Comparative analysis reveals signatures of differentiation amid genomic polymorphism in Lake Malawi cichlids
© Loh et al.; licensee BioMed Central Ltd. 2008
Received: 25 April 2008
Accepted: 10 July 2008
Published: 10 July 2008
Cichlid fish from East Africa are remarkable for phenotypic and behavioral diversity on a backdrop of genomic similarity. In 2006, the Joint Genome Institute completed low coverage survey sequencing of the genomes of five phenotypically and ecologically diverse Lake Malawi species. We report a computational and comparative analysis of these data that provides insight into the mechanisms that make closely related species different from one another.
We produced assemblies for the five species ranging in aggregate length from 68 to 79 megabase pairs, identified putative orthologs for more than 12,000 human genes, and predicted more than 32,000 cross-species single nucleotide polymorphisms (SNPs). Nucleotide diversity was lower than that found among laboratory strains of the zebrafish. We collected around 36,000 genotypes to validate a subset of SNPs within and among populations and across multiple individuals of about 75 Lake Malawi species. Notably, there were no fixed differences observed between focal species nor between major lineages. Roughly 3% to 5% of loci surveyed are statistical outliers for genetic differentiation (FST) within species, between species, and between major lineages. Outliers for FST are candidate genes that may have experienced a history of natural selection in the Malawi lineage.
We present a novel genome sequencing strategy, which is useful when evolutionary diversity is the question of interest. Lake Malawi cichlids are phenotypically and behaviorally diverse, but they appear genetically like a subdivided population. The unique structure of Lake Malawl cichlid genomes should facilitate conceptually new experiments, employing SNPs to identity genotype-phenotype association, using the entire species flock as a mapping panel.
Cichlid fishes from the East African Rift lakes Victoria, Tanganyika, and Malawi represent a preeminent example of replicated and rapid evolutionary radiation . This group of fishes is a significant model of the evolutionary process and the coding of genotype to phenotype, largely because tremendous diversity has evolved in a short period of time among lineages with similar genomes [2–4]. Recently evolved cichlid species segregate ancestral polymorphism [5, 6] and may exchange genes [7, 8]. Numerous genomic resources have been developed for East African cichlids (many of which are summarized by the Cichlid Genome Consortium ). These include the following: genetic linkage maps for tilapia [10–12] and Lake Malawi species [10, 13]; fingerprinted bacterial artificial chromosome libraries ; expressed sequence tag sequences for Lake Tanganyika and Lake Victoria cichlids ; and first-generation microarrays [16, 17]. Many studies have used these resources to study cichlid population genetics, molecular ecology, and phylogeny (for review [18, 19]). Recent reports have capitalized on the diversity among East African cichlids to study the evolution and genetic basis of many traits, including behavior , olfaction , pigmentation [22–24], vision [25, 26], sex determination [24, 27], the brain , and craniofacial development [10, 13, 29].
In 2006, under the auspices of the Community Sequencing Program, the Joint Genome Institute (JGI) completed low coverage survey sequencing of the genomes of five Lake Malawi species. Species were chosen to maximize the morphological, behavioral, and genetic diversity among the Malawi species flock. This represents a novel genome project. Low coverage sequencing is now a routine strategy to uncover functional or 'constrained' genomic elements . The rationale is as follows; one compares genome sequences of distantly related organisms (for example, shark, diverse mammals) with that of a reference (for instance, human, mouse), and outliers of similarity will be observed against the background expectation of divergence [31–34]. Our interests in diversity suggest a conceptually similar but logically reversed research objective. When the background expectation is similarity, how does one use low coverage genome sequencing to detect that which makes organisms distinct?
Here, we report computational and comparative analyses of survey sequence data to address the question of diversity. We had four major goals: to produce a low coverage assembly for each of the five Lake Malawi species; to identify orthologs of vertebrate genes in these data; to predict single nucleotide polymorphisms (SNPs) segregating between species; and to use SNPs to evaluate the degree of genomic polymorphism and divergence at different evolutionary scales. Consequently, we produced assemblies for the five species ranging in aggregate length from 68 to 79 megabases (Mb), identified putative orthologs for more than 12,000 human genes, and predicted more than 32,000 cross-species segregating sites (with about 2,700 located in genic regions). We genotyped a set of these SNPs within and between Lake Malawi cichlid lineages and demonstrate signatures of differentiation on the background of similarity and polymorphism. Our work should facilitate further understanding of evolutionary processes in the species flocks of East African cichlids. Moreover, the approach we outline should be broadly applicable in other lineages where phenotypic and behavioral diversity has evolved in a short window of evolutionary time.
First-pass genomic assembly statistics for five Lake Malawi cichlid species
Total number of contigs in assembly
Total length (bases)
Genome coveragea (%)
Mean trace length (bases)
Shortest contig length (bases)
Longest contig length (bases)
Mean contig length (bases)
Q25 contig length (bases)
Q50 (median) contig length (bases)
Q75 contig length (bases)
Total genic length (bases)
We noted that these first-pass assemblies were 'over-assembled' by roughly a factor of 2 when compared with theoretical expectations . Theory suggests that random shotgun sequencing of single copy DNA, at 15% coverage of a 1.1 gigabase genome, will result in an assembly length of about 153 Mb. We reasoned that our assemblies might be shorter than expected because multicopy elements were grouped as if they were single copy sequence. Given the theoretical expectation (again for 15% coverage of a 1.1 gigabase genome) that individual bases should only be sequenced a maximum of four to five times, we examined whether contigs were built from five or more trace sequences contributing overlapping bases. We observed that about 10 Mb of each first-pass assembly were derived from such contigs, and excluded these data from subsequent analyses (for example SNP prediction [see below]). Notably, individual sequences contributing to these 'high trace number' contigs were not identified by RepeatMasker but did sometimes have Basic Local Alignment Search Tool (BLAST) matches to putative repetitive elements (for example, pol polyprotein, reverse transcriptase). Because of the keen interest in repetitive DNA families in cichlids  and other organisms , we have retained alignments of these 'high trace number' contigs and have marked them as such (see Additional data files 3 and 4).
Gene content and coverage
To establish the extent of gene content and coverage present in each assembly, we carried out BLASTX similarity searches (10-10 E value cutoff) for each of the five assemblies against a reference human proteome (RefSeq proteins). The average proportion of putative genic sequence amounted to 3.9% of the available genomes. The MZ assembly contained the highest gene coverage, possessing genic loci that were significantly similar to approximately 5,240 unique human proteins. The remaining four species yielded approximately similar numbers ranging from 5,020 to 5,170 genes. It must be noted, however, that most of these genes are highly fragmented and incomplete, because of low coverage of the assembly. In all, a total of 36% (12,211 genes out of 34,180; see Additional data file 2) of the reference human proteome could be identified in one or more of the cichlid species.
Clustering and alignment
Roughly 1% of the alignments (294 alignments) showed percentages of variable sites above 2% (about tenfold higher than the average). It is impossible to know, given the low coverage of the sequenced genomes, whether these represent orthologous but divergent regions of cichlid genomes or the alignment of paralogous sequence. We therefore retained these alignments, and included a calculation of polymorphism for each alignment (see Additional data file 3), for the consideration of researchers using these data. For example, alignment 108,866 contains sequence with similarity to asteroid homolog 1, with 8% of sites variable and a majority of replacement polymorphism. Given the lack of functional information about this novel signaling protein (first described in Drosophila ), this alignment provides useful information even if (and perhaps because) it includes paralogous loci. Another 12% of the alignments (2,119 total) contained individual species contigs that had consensus base positions derived from five or more trace sequences (see above).
For all subsequent analyses, we excluded 2,413 alignments that exhibited a high percentage of variable sites and/or higher than expected coverage. More than 11.6 million bases of multiple species alignments remain, of which roughly 1.06 Mb were inferred as genic. This included 10,902,011 (986,506 genic) bases of two-species alignments, 721,049 (75,371 genic) bases of three-species alignments, 27,951 (2,898 genic) bases of four-species alignments, and 877 (193 genic) bases of alignments containing all five species.
Further analysis of these 11.6 million bases of multiple alignments identified a total of 32,417 (0.28%) cross-species SNPs. In order to classify the quality of an identified variable site, a polymorphism quality score (PQS) was defined, corresponding to the first digit of the lowest Phrap quality score among the nucleotides of the different species present at the polymorphic site (for example, a polymorphic site between four species with base quality scores of 34, 45, 46, and 50 would be assigned a PQS of 3). In total, 4,468 (13.8%) variable sites had a PQS of 5 or higher, 7,952 (24.5%) had a PQS of 4, 8,236 (25.4%) a PQS of 3, and the remaining 11,761 (36.3%) had a PQS of 2. PQS for each variable site are provided on the alignments described in Additional data file 3 (also available online ). Nucleotide diversity (Watterson's θw) averaged over two-, three-, and four-species alignments was 0.00257. Roughly 8% of all polymorphic sites (2,709) were located within the putative genic regions identified earlier. Alignments with fish and human proteins provided us with the phase information required to further classify these into 1,066 synonymous and 1,643 nonsynonymous SNPs. Summaries of all alignments containing genic and nongenic polymorphisms are provided in Additional data files 3 and 4.
In order to investigate the pair-wise differences between any two of the five species, all sequence alignment segments with two or more species were broken up into all possible pair-wise alignments; this resulted in 1.06 to 1.55 Mb of alignment per pair. We then calculated the Jukes-Cantor distance between species pairs. The three shortest distances were between LF and MZ (0.229%), followed by MA/MZ (0.232%) and LF/MA (0.241%), and the greatest was between LF and RE (0.288%). These genetic distances include both within-species polymorphism and the fixed differences between species. Currently, there is no exhaustive estimate of within-species polymorphism for Malawi cichlids. Unpublished data from our own group (Streelman JT) indicates that for LF and MZ, within-species diversity (π) may be as high as 0.2%. Thus, the percentage of fixed genetic differences is likely to be extremely small in this assemblage (see following sections).
Finally, we calculated the ratio of replacement to synonymous substitutions (Ka/Ks) for concatenated genic alignments among all pairs of species. We used concatenated sequences because each segment represented only a small fraction of a gene, with only few nonsynonymous and synonymous sites. Ka/Ks ranged from 0.380 in MC/LF to 0.562 in LF/MA. These numbers are greater than the ratios found between Fugu and Tetraodon (0.127 to 0.144 ). Such high Ka/Ks values may indicate that positive selection, driven by adaptive radiation, is prevalent in cichlid fishes. However, given the expectation of few fixed differences between groups, this topic should be revisited with more data on the levels of segregating and fixed nucleotide substitutions among lineages.
Validation and generality of SNPs
We genotyped 96 SNPs in 384 Lake Malawi cichlid samples using Beckman Coulter SNPstream™ technology (Beckman Coulter, Inc., Fullerton, CA). The SNPs were partitioned into three categories to help us evaluate the comparative success rate of automated SNP prediction. First, we included 13 positive controls: genes previously sequenced by others [3, 25] and by us (Streelman JT, unpublished data), with expected variation in Malawi cichlids. Positive controls included genes involved in morphogenesis (otx1, otx2, and pax9), pigmentation (mitf, ednrb, and aim1), and visual sensitivity (opsins rh1, sws1, lws, sws2a, and sws2b). Next, we genotyped 59 SNPs identified using the automated procedure described in this report. We selected these SNPs to represent a range of PQS (from 2 to 5) and a variety of sequence types (genic, nongenic with a BLAST match < e-100 to Tetraodon, and nongenic with no BLAST match). Finally, we wished to compare our automated SNP selection to a manual approach. Therefore, we included an additional 24 SNPs identified by manual inspection of BLAST matches between single JGI traces and Tetraodon chromosome 11; we have previously shown Tetraodon 11 to share orthologs with cichlid chromosome 5 . Note that these SNPs were most often not discovered by our automated procedure because they originated in single traces that did not meet percentage quality cutoffs and/or they did not align into comparative contigs because of overlap cutoffs.
Our validation strategy sought to document the general use and segregation of these markers among Lake Malawi cichlids. Given recent divergence times among species (some as recent as 1,000 years ), we expected that SNPs might segregate throughout the assemblage. Therefore, Malawi samples comprised about ten individuals from each of ten populations of MZ and LF, as well as one to five individuals of 77 additional species (25 of which were rock-dwelling mbuna). Taxa were included to represent the morphological, functional, and behavioral diversity of the Malawi lineage, which may contain more than 800 species .
SNP genotyping success categorized by detection method
SNP detection method
Number of genotyped loci
Number of polymorphic loci
Number of fixed loci
Number of failed loci
Successful SNP detection (%)
SNP genotyping success categorized by polymorphic quality score
Polymorphic quality score
Number of genotyped loci
Number of polymorphic loci
Number of fixed/failed loci
Successful SNP detection (%)
Genetic polymorphism and divergence at multiple scales
Strikingly, among all 68 loci showing polymorphism, no SNP locus was alternately fixed between LF and MZ, or between rock-dwelling mbuna and non-mbuna. We thus sought to investigate the degree of polymorphism versus divergence at multiple evolutionary scales.
Genetic clustering and ancestry
For comparison, additional analyses were performed setting the predefined number of genetic clusters to from two to five. When set to two genetic clusters, species were accurately classified as mbuna or non-mbuna. At settings of four or five, the program was unable to yield stable classification results between replicate runs. Thus, these latter three sets of analyses (data not shown) did not provide any further insights into the genetic lineages of Malawi cichlids.
African cichlid fishes are important models of evolutionary diversification in form and function . They are singularly remarkable for the extent of phenotypic and behavioral diversity on a backdrop of genomic similarity. Lake Malawi is home to the most species-rich assemblage of African cichlids; as many as 800 to 1,000 species are thought to have evolved from a common ancestor during the past 500,000 to 1 million years ago . These recently formed species segregate ancestral polymorphism and exchange genes by hybridization [5, 7, 49]. Such circumstances present both opportunities and challenges for understanding evolutionary history and biological diversity. Opportunistically, researchers have used molecular markers across studies to interrogate the genetic basis of phenotypic differentiation [13, 22, 24, 29]. This approach views Malawi cichlid species as natural mutants screened for function by natural selection, with essentially identical ancestral genomes honed by contrasting historical processes. By contrast, the task of reconstructing a phylogeny of species has been hindered by the very same phenomena of genomic similarity and mosaicism [2, 3]; even the promising approach of Amplified Fragment Length Polymorphism (AFLP) does not provide strong resolution of the relationships among genera [23, 48, 50, 51]. The data we present here should provide new resources and perspectives for cichlid evolutionary genomics.
Cichlid species exhibit genomic polymorphism
Lake Malawi cichlid species sequenced by the JGI embody the phylogenetic, morphological, and behavioral diversity found within the assemblage. Rhamphochromis esox (RE) is a large (about 0.5 m) pelagic predator that represents one of the basal lineages of the species flock [3, 4, 47]. Mchenga conophorus (MC) is a sand-dwelling species that breeds on leks, where males construct 'bowers' to attract females. Melanochromis auratus (MA), Maylandia zebra (MZ), and Labeotropheus fuelleborni (LF) are rock-dwelling (mbuna) species that differ in color pattern, trophic ecology, body shape, and craniofacial morphology (pictures of these and others are available online ).
Our data confirm the conclusions from previous genetic analyses on a smaller scale; Lake Malawi species are genetically similar. Nucleotide diversity observed among the five cichlid species (Watterson's θw = 0.26%) is less than that found among laboratory strains of the zebrafish Danio rerio (Watterson's θw = 0.48% ). Although overall nucleotide diversity is less than that observed in Danio, the ratio of replacement to silent change is nearly fivefold higher in the Lake Malawi genomes. Such a result might suggest that East African cichlid evolution is characterized by adaptive molecular evolution, as has been indicated in a few instances [25, 54], or a relaxation of purifying selection attributable to small effective population size. However, we should view this estimate of Ka/Ks with caution because of one of the remarkable features of these data (see below). Variable sites identified from cross-species alignments are not substitutions fixed between species. The Ka/Ks approach to identifying selection may be largely inappropriate for such young species where ancestral alleles segregate as polymorphisms.
The pattern of variation observed across the approximately 75 species genotyped in this study demonstrates that biallelic polymorphisms segregate widely throughout the Malawi species flock. SNPs segregate within and between MZ and LF populations, as well as within and among mbuna species and other lineages. No SNP locus surveyed is alternately fixed in LF versus MZ, nor between mbuna and non-mbuna. Remarkably, the degree of genetic differentiation (FST) within species is roughly equivalent to that between species and to that between major lineages. Lake Malawi cichlid species are mosaics of ancestrally polymorphic genomes. Add to this a propensity of recently diverged species to exchange genes , and Malawi cichlids present a case of complex and dynamic evolutionary diversification, where recombination and the sorting of ancestral polymorphism may be more important than new mutation as sources of genetic variation. Despite allele sharing, SNP frequencies contain a clear signal of ancestry for the entire flock. Rock-dwelling mbuna comprise a genetic cluster, as do pelagic and sand-dwelling species, in addition to Rhamphochromis. Notably, Astatotilapia calliptera, one of a few nonendemic haplochromines in Lake Malawi, appears to retain a reservoir of ancestral polymorphisms from which mbuna and non-mbuna genomes have emerged.
Genomic polymorphism and the divergence of Malawi cichlids
Our hierarchical sampling design allows us to consider whether there are loci exhibiting extreme genetic differentiation against the background of shared polymorphism within species, between species, and between major lineages. Strikingly, regardless of the evolutionary scale, statistical outliers comprise approximately 3% to 5% of loci surveyed. Opsin loci lws, rh1, and sws1 are differentiated among populations of LF and MZ, adding to reports that opsin polymorphisms are associated with population-specific color patterns or visual environments .
SNPs in csrp1, β-catenin, and ptc2 exhibit greater than expected differentiation between LF and MZ. Csrp1 (cysteine-rich protein) is a vertebrate LIM-domain family member acting in the noncanonical WNT pathway, expressed in gut, intestine, and cardiac mesoderm . β-catenin acts to transduce signals in the canonical WNT pathway  and is expressed in developing cichlid fins, dentitions, brains, and lateral lines (Fraser GJ, Streelman JT, unpublished data). Patched is a receptor for sonic hedgehog ; both areexpressed in developing cichlid dentitions, jaws, and brains (Fraser GJ, Sylvester JB, Streelman JT, unpublished data). A SNP in irx1 nearly perfectly differentiates rock-dwelling mbuna from the remainder of the Malawi species flock. Irx1 acts to position the boundary between the telencephalon and the posterior forebrain . Finally, a SNP located between contactin 3 and ncam L1 exhibits differentiation between mbuna and non-mbuna lineages; these genes are linked in other genomes and functionally interact to pattern dendritic branching in the neocortex . Taken together, differentiated loci are interesting in the context of cichlid diversification because they affect the phenotypes that vary among lineages: color and vision [25, 26], guts , dentitions [13, 62], jaws [10, 29], and brains .
Discovery for evolutionary biology
There are obvious challenges when attempting to extract information from low coverage genomic sequence, and also obvious payoffs [31–34]. Most previous studies have used this information for species-specific discovery (for example, dog breeds) or broad evolutionary comparisons with respect to a reference genome (for example, dog-human, shark-human, or cat-mammal). Our goals in the present analysis stem from the unique characteristics of Lake Malawi cichlids; these are biologic species that behave genetically like a single subdivided population. Therefore, our biggest challenge was to devise a strategy that retains information from these low coverage survey sequences (75% genomic coverage spread over five closely related species), but minimizes error and bias in assembly and cross-species alignment for SNP identification. For example, we excluded many contigs because they appeared to be over-assembled, and we excluded multi-species alignments if they exceeded a polymorphism threshold. The over-assembly problem limits the coverage of these genomes in relation to expectation; this phenomenon, observed in the cat genome and in simulation, has complex and varying causes and has yet to be fully resolved . It is likely to be mitigated to some degree by comparison with a higher coverage reference sequence. The power of the data we present comes from the broad utility of the genic sequences and SNPs we have identified for many questions in genomic evolutionary biology.
Our analyses identified about 12,000 Lake Malawi cichlid sequences with similarity to human and fish proteins. This is a significant advance in our understanding of cichlid genomic content. To put this in context, approximately 13,500 unique expressed sequence tags, from three different East African cichlids, represent the sum total of such publicly released sequences . Our contribution roughly doubles the available data.
The approximately 32,000 (2,700 genic) SNPs we identified should provide a wealth of molecular markers for studies of population genetics and molecular ecology, linkage and quantitative trait locus mapping, association mapping, and phylogeny. We convert about 70% of predicted SNPs to polymorphic markers; this percentage is comparable to that of other studies from white spruce (74% to 85%, depending on quality cutoffs ), zebrafish (65% ), and cow (43% ). We have shown these biallelic markers to be of general use, many segregating across the major cichlid lineages of Lake Malawi. We used the SNPs to assign Malawi species to ancestral genetic clusters, and this approach should hold promise for similar questions of genetic structure that span the population versus species continuum. It is important to note that early runs of this analysis, with fewer SNP loci, resulted in stable results with more individuals showing mosaic genomes. This suggests that careful consideration should be given to the number of polymorphic loci necessary to yield confidence in evolutionary interpretation. As more SNP loci (with known genome coordinates) are assayed, it will be possible to compute and compare ancestry proportions across scales (for example, genome versus chromosome versus gene cluster).
Notably, we have used the background level of genomic similarity and polymorphism to identify loci that may have experienced a history of selection within species, between species and between major lineages. Because SNP markers are co-dominant, easy to genotype, reliable and reproducible from laboratory to laboratory, and readily mapped in silico (NHGRI will sequence a related cichlid, the tilapia, to 7-fold draft assembly coverage in 2008), they are likely to complement microsatellites and AFLP for most applications in cichlid evolutionary genomics. Given the unique mosaic structure of Lake Malawl cichlid genomes, it is exciting to envision experiments employing SNPs to identity genotype-phenotype associations, using the entire species flock as a mapping panel. Finally, as sequencing costs continue to drop, the approach we outline here should prove applicable to those studying evolutionary and phenotypic diversity among closely related species .
Materials and methods
Individuals of Mchenga conophorus (MC), Labeotropheus fuelleborni (LF), Melanochromis auratus (MA), Maylandia zebra (MZ), and Rhamphochromis esox (RE) were sampled from the wild during an expedition to Malawi in 2005. Specimens prepared for survey sequencing by the JGI were collected from Mazinzi Reef (MZ), Domwe Island (LF and MA), and Otter Point (MC and RE), all of which are locales in the southeastern portion of the lake. High-quality DNA was extracted and prepared in the laboratory of TDK.
Trace sequences generated by the JGI for MC, LF, MA, MZ, and RE, together with their sequence quality scores, were downloaded (6 May 2007) from the National Center for Biotechnology Information (NCBI) Trace Archive. The dataset for each species consisted of an average of about 152,000 individual trace reads with total read lengths ranging from 137 to 185 million bases. Detailed sequence statistics for each species are provided in Additional data file 1.
Sequence preprocessing and assembly
The trace and quality sequences were first pre-processed for assembly by masking out all possible vector sequences available from the NCBI UniVec vector sequence database (downloaded 6 May 2007). The vector masking was performed using the cross_match.pl perl script provided by the Phred-Phrap package . In order to reduce the computational complexity and time required for the final assembly, repeat sequences were masked before assembly using RepeatMasker version 3.1.8 (Smit AFA, Hubley R and Green P, unpublished data) in conjunction with the latest repeatmasker libraries from RepBase Update . Bases with sequencing quality score of less than 20 were also masked. The actual assembly of each species' trace sequences into contiguous sequences (contigs) was then performed using the Phrap version 0.990329 assembly program from the Phred-Phrap package. Contigs with more than 80% low quality bases (defined as <20 assembly quality score) were removed from the assembly. This whole genome shotgun project has been deposited at DDBJ/EMBL/GenBank under the project accessions ABPJ00000000 (MC), ABPK00000000 (LF), ABPL00000000 (MA), ABPM00000000 (MZ), and ABPN00000000 (RE). The versions described in this paper are the first versions: ABPJ01000000, ABPK01000000, ABPL01000000, ABPM01000000, and ABPN01000000.
Similarity search and alignment
Orthologous genomic contig pairs were first identified using reciprocal BLASTN similarity searches with a strict E-value cutoff of 10-100, performed across the sequence contigs of all possible species pairs. To reduce spurious ortholog assignments, putative ortholog contig pairs were only retained if their regions of high sequence similarity formed good end-to-end overlaps (defined as within 100 bases of the 5' end or 30 bases from the 3' end of a sequence) or overlap more than 80% of the shorter contig. Although some of the filtered regions could represent biologically relevant loci where recombination or translocations might have occurred, we decided to remove them from this analysis. Contig pair assignments were then passed to an algorithm that created clusters of contigs whereby each contig within the cluster must be related to all other contigs in the cluster through one or more putatively orthologous relations.
Each cluster of contigs was then individually aligned using Phrap, resulting in a continuous alignment tiling path where each alignment position may consist of a base from any one or up to all five cichlid species (Figure 1). Segregating sites were then identified from alignment positions with high quality bases (>20 score) from two or more species. A PQS was defined, corresponding to the first digit of the lowest Phrap quality score among the nucleotides of the different species present at the polymorphic site (for example, a polymorphic site between four species with base quality scores of 34, 45, 46, and 50 would be assigned a PQS of 3). To compare the extent of nucleotide diversity among the five cichlid species, we calculated Watterson's theta (θw ). This measure takes into account the number of variable positions and the sample size analyzed. Our data violate the assumption of an infinite, interbreeding population, but we chose this metric to in order to make direct comparisons to similar measures from study of other genomes (for example, zebrafish).
Protein-coding sequence identification
Cichlid protein coding sequences were inferred based on similarity searches to known protein databases of fishes and humans. BLASTX searches with E-value cutoff of 10-10 were performed for the each cichlid genomic assembly as well as the overall consensus sequence of the cluster alignments, against a protein database made up of all GenBank Actinopterygii (ray-finned fishes) sequences (downloaded 2 June 2007; 163,471 entries) and all human RefSeq proteins (downloaded 25 June 2007; 34,180 sequences). The alignment with the highest scoring hit for each genomic locus was then used as a reference to determine the coding strand and phase of the protein-coding cichlid locus.
Evolutionary sequence divergence among JGI species
All cluster alignment segments with contributing bases from two or more species were split into pairwise alignments (each two, three, four, or five species alignment position can be split into one, three, six, or ten pair-wise alignments respectively). Pair-wise alignments within each of the ten possible species pair combinations (MC-LF, MC-MA, MC-MZ, MC-RE, LF-MA, LF-MZ, LF-RE, MA-MZ, MA-RE, and MZ-RE) were then concatenated and the number of substitutions counted. Jukes-Cantor correction for multiple substitutions was applied to these direct distance measurements . Pair-wise alignments consisting of only genic sequences were obtained from multi-species cluster alignment segments in a manner similar to that described above. The DNAStatistics package of Bioperl  was then used to calculate the Ka/Ks values of pair-wise alignments.
Genotyping and validation of SNPs
We genotyped 96 SNPs in 364 diverse Lake Malawi cichlid samples. These SNPs included 13 positive controls, 59 loci from the automated procedure described in this report, and an additional 24 loci chosen manually by BLAST of individual traces to the Tetraodon genome (see main text for further description). The GenomeLab SNPstream Genotyping System Software Suite v2.3 (Beckman Coulter, Inc.) was used for experimental setup, data uploading, image analysis, genotype calling and QC review, at Emory University's Center for Medical Genomics. In brief, marker panel data (multiplexed SNP panel designed by SNPstream's Primer Design Engine website ) were first uploaded to the SNPstream database using the PlateExplorer application software. Also uploaded was the Process Group Data containing all test sample information generated through a Laboratory Information Management System (Nautilus 2002; Thermo Fisher Scientific, Waltham, MA, USA). An on-board CCD camera of the SNPstream Imager took two snapshot images of each well of the 384-well tag array, one under a blue excitation laser and the other under a green excitation laser. Image application software was used to analyze the captured images to detect spots, overlay an alignment grid, and determine spot intensity. The fluorescent pixel intensity data for each SNP under the two channels, representing the relative abundance of the two alleles, were uploaded to the database. The GetGenos application software was used to calculate and generate a Log(B+G) versus B/(B+G) plot, where B and G were the pixel intensities under the blue and green channels, respectively, for each sample and each SNP. Next, automated genotype calling was accomplished using the QCReview application software based on a number of criteria (for instance, signal baseline, clustering pattern of the three genotypes, and Hardy-Weinberg score). A genotype summary was generated using the Report application software.
Genetic differentiation within and among lineages
Locus-specific FST  was calculated using FSTAT version 126.96.36.199  for three evolutionary scales: within LF and MZ; between LF and MZ; and between mbuna and non-mbuna. We determined that a SNP locus was a statistical outlier using the empirical distribution of FST values. FST outliers exceed the sum of the upper quartile value and 1.5 times the interquartile range.
We used a Bayesian method (STRUCTURE v.2.2 ) to determine how well our SNP genotypes assigned individuals to evolutionary lineages. We chose to define the number of K genetic clusters in accord with previous research showing about three major evolutionary groups of Lake Malawi cichlids [3–5, 47]. Note that we do not intend this to mean that three is the best supported estimate of K in these data; our rationale is rather to demonstrate how individual genomes are composites (or not) of the major evolutionary lineages found in the lake. Thus, we used the admixture model to estimate q, the proportion of each genome derived from each of K genetic clusters. For comparison, we also ran analyses with K set to two, four, or five (not shown). Each run of the program included 50,000 cycles of burn-in and run length of 50,000 steps. Multiple runs were conducted to ensure reliability and consistency of results.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is a table of trace sequence statistics of five Lake Malawi cichlid species. Additional data file 2 is a list of human gene homologs found in the five cichlid species. Additional data file 3 is a list of alignments and polymorphic sites. Additional data file 4 is a list of alignments with BLAST hits to fish and humans. Additional data file 5 is a table of major allele frequencies for biallelic SNPs surveyed across Lake Malawi cichlid populations and species.
Basic Local Alignment Search Tool
Joint Genome Institute
ratio of replacement to synonymous substitutions
National Center for Biotechnology Information
polymorphism quality score
single nucleotide polymorphism.
We thank members of the Streelman laboratory, Karen Carleton, and two anonymous reviewers for comments on previous drafts of the manuscript. The research is supported by grants from the NSF (IOS 0546423), NIH (R21 DE017182), and Alfred P Sloan Foundation (BR-4499) to JTS. Drs Karen Carleton and Federica DiPalma extracted high-quality DNA from the five species of Malawi cichlid. Library construction and sequencing was performed by the JGI under the auspices of the US Department of Energy's Office of Science, Biological and Environmental Research Program and by the University of California, Lawrence Livermore National Laboratory under contract number W-7405-Eng-48, Lawrence Berkeley National Laboratory under contract number DE-AC03-76SF00098 and Los Alamos National Laboratory under contract number W-7405-ENG-36.
- Kocher TD: Adaptive evolution and explosive speciation: the cichlid fish model. Nat Rev Genet. 2004, 5: 288-298. 10.1038/nrg1316.PubMedView ArticleGoogle Scholar
- Won Y-J, Sivasundar A, Wang Y, Hey J: On the origin of Lake Malawi cichlid species. Proc Natl Acad Sci USA. 2005, 102: 6581-6586. 10.1073/pnas.0502127102.PubMedPubMed CentralView ArticleGoogle Scholar
- Won Y-J, Wang Y, Sivasundar A, Raincrow J, Hey J: Nuclear gene variation and molecular dating of the cichlid species flock of Lake Malawi. Mol Biol Evol. 2006, 23: 828-837. 10.1093/molbev/msj101.PubMedView ArticleGoogle Scholar
- Hulsey CD, Mims MC, Streelman JT: Do constructional constraints influence cichlid craniofacial diversification?. Proc Biol Sci. 2007, 274: 1867-1875. 10.1098/rspb.2007.0444.PubMedPubMed CentralView ArticleGoogle Scholar
- Moran P, Kornfield I: Retention of ancestral polymorphism in the Mbuna species flock of Lake Malawi. Mol Biol Evol. 1993, 10: 1015-1029.Google Scholar
- Nagl S, Tichy H, Mayer WE, Takahata N, Klein J: Persistence of neutral polymorphisms in Lake Victoria cichlid fish. Proc Natl Acad Sci USA. 1998, 24: 14238-14243. 10.1073/pnas.95.24.14238.View ArticleGoogle Scholar
- Smith PF, Konings A, Kornfield I: Hybrid origin of a cichlid population in Lake Malawi: implications for genetic variation and species diversity. Mol Ecol. 2003, 12: 2497-2504. 10.1046/j.1365-294X.2003.01905.x.PubMedView ArticleGoogle Scholar
- Seehausen O: Hybridization and adaptive radiation. Trends Ecol Evol. 2004, 19: 198-207. 10.1016/j.tree.2004.01.003.PubMedView ArticleGoogle Scholar
- Cichlid Genome Consortium. [http://www.cichlidgenome.org]
- Albertson RC, Streelman JT, Kocher TD: Directional selection has shaped the oral jaws of Lake Malawi cichlid fishes. Proc Natl Acad Sci USA. 2003, 100: 5252-5257. 10.1073/pnas.0930235100.PubMedPubMed CentralView ArticleGoogle Scholar
- Kocher TD, Lee W-J, Sobolewska H, Penman D, McAndrew B: A genetic linkage map of the cichlid fish, the tilapia (Oreochromis niloticus). Genetics. 1998, 148: 1225-1232.PubMedPubMed CentralGoogle Scholar
- Carleton KL, Streelman JT, Lee B-Y, Garnhart N, Kidd MR, Kocher TD: Rapid isolation of CA microsatellites from the cichlid genome. Anim Genet. 2002, 33: 140-144. 10.1046/j.1365-2052.2002.00817.x.PubMedView ArticleGoogle Scholar
- Streelman JT, Albertson RC: Evolution of novelty in the cichlid dentition. J Exp Zoolog B Mol Dev Evol. 2006, 306: 216-226. 10.1002/jez.b.21101.View ArticleGoogle Scholar
- Katigiri T, Kidd CE, Tomasino E, Davis JT, Wishon C, Stern JE, Carleton KL, Howe AE, Kocher TD: A BAC-based physical map of the Nile tilapia genome 89. BMC Genomics. 2005, 9: 89-10.1186/1471-2164-6-89.View ArticleGoogle Scholar
- The Gene Index Project. [http://compbio.dfci.harvard.edu/tgi]
- Kijimoto T, Watanabe M, Fujimura K, Nakazawa M, Murakami Y, Kuratani S, Kohara Y, Gojobori T, Okada N: cimp, a novel astacin family metalloproteinase gene from East African cichlids, is differentially expressed between species during growth. Mol Biol Evol. 2005, 22: 1649-1660. 10.1093/molbev/msi159.PubMedView ArticleGoogle Scholar
- Renn SC, Aubin-Horth N, Hofmann HA: Biologically meaningful expression profiling across species using heterologous hybridization to a cDNA microarray. BMC Genomics. 2004, 6: 42-10.1186/1471-2164-5-42.View ArticleGoogle Scholar
- Kornfield I, Smith PF: African cichlid fishes: model systems for evolutionary biology. Ann Rev Ecol Evol Syst. 2000, 31: 163-196.View ArticleGoogle Scholar
- Genner MJ, Turner GF: The mbuna cichlids of Lake Malawi: a model for rapid speciation and adaptive radiation. Fish Fisheries. 2005, 6: 1-34. 10.1111/j.1467-2679.2005.00173.x.View ArticleGoogle Scholar
- Aubin-Horth N, Desjardins JK, Martei YM, Balshine S, Hofmann HA: Masculinized dominant females in a cooperatively breeding species. Mol Ecol. 2007, 16: 1349-1358. 10.1111/j.1365-294X.2007.03249.x.PubMedView ArticleGoogle Scholar
- Blais J, Rico C, van Oosterhout C, Cable J, Turner GF, Bernatchez L: MHC adaptive divergence between closely related and sympatric African cichlids. PloS ONE. 2007, 2: e734-10.1371/journal.pone.0000734.PubMedPubMed CentralView ArticleGoogle Scholar
- Streelman JT, Albertson RC, Kocher TD: Genome mapping of the orange blotch colour pattern in cichlid fishes. Mol Ecol. 2003, 12: 2465-2471. 10.1046/j.1365-294X.2003.01920.x.PubMedView ArticleGoogle Scholar
- Allender CJ, Seehausen O, Knight ME, Turner GF, Maclean N: Divergent selection during speciation of Lake Malawi cichlid fishes inferred from parallel radiations in nuptual coloration. Proc Natl Acad Sci USA. 2003, 100: 14074-14079. 10.1073/pnas.2332665100.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee B-Y, Lee W-J, Streelman JT, Carleton KL, Howe AE, Hulata G, Slettan A, Stern JE, Terai Y, Kocher TD: A second-generation genetic linkage map of tilapia (Oreochromis spp). Genetics. 2005, 170: 237-244. 10.1534/genetics.104.035022.PubMedPubMed CentralView ArticleGoogle Scholar
- Spady TC, Seehausen O, Loew ER, Jordan RC, Kocher TD, Carleton KL: Adaptive molecular evolution in the opsin genes of rapidly speciating cichlid fishes. Mol Biol Evol. 2005, 22: 1412-1422. 10.1093/molbev/msi137.PubMedView ArticleGoogle Scholar
- Parry JW, Carleton KL, Spady T, Carboo A, Hunt DM, Bowmaker JK: Mix and match color vision: tuning spectral sensitivity by differential gene expression in Lake Malawi cichlids. Curr Biol. 2005, 15: 1734-1739. 10.1016/j.cub.2005.08.010.PubMedView ArticleGoogle Scholar
- Lee B-Y, Hulata G, Kocher TD: Two unlinked loci controlling the sex of blue tilapia (Oreochromis aureus). Heredity. 2004, 92: 543-549. 10.1038/sj.hdy.6800453.PubMedView ArticleGoogle Scholar
- Huber R, van Staaden MJ, Kaufman LS, Liem KF: Microhabitat use, trophic patterns, and the evolution of brain structure in African cichlids. Brain Behav Evol. 1997, 50: 167-182. 10.1159/000113330.PubMedView ArticleGoogle Scholar
- Albertson RC, Streelman JT, Kocher TD, Yelick PC: Integration and evolution of the cichlid mandible: the molecular basis of alternative feeding strategies. Proc Natl Acad Sci USA. 2005, 102: 16287-16292. 10.1073/pnas.0506649102.PubMedPubMed CentralView ArticleGoogle Scholar
- Margulies EH, Birney E: Approaches to comparative sequence analysis: towards a functional view of vertebrate genomes. Nat Rev Genet. 2008, 9: 303-313. 10.1038/nrg2185.PubMedView ArticleGoogle Scholar
- Kirkness EF, Bafna V, Halpern AL, Levy S, Remington K, Rusch DB, Delcher AL, Pop M, Wang W, Fraser CM, Venter JC: The dog genome: survey sequencing and comparative analysis. Science. 2003, 310: 1898-1903. 10.1126/science.1086432.View ArticleGoogle Scholar
- Venkatesh B, Kirkness EF, Loh YH, Halpern AL, Lee AP, Johnson J, Dandona N, Viswanathan LD, Tay A, Venter JC, Strausberg RL, Brenner S: Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome. PloS Biology. 2007, 5: e101-10.1371/journal.pbio.0050101.PubMedPubMed CentralView ArticleGoogle Scholar
- Pontius JU, Mullikin JC, Smith DR, Agencourt Sequencing Team, Lindblad-Toh K, Gnerre S, Clamp M, Chang J, Stephens R, Neelam B, Volfovsky N, Schäffer AA, Agarwala R, Narfström K, Murphy WJ, Giger U, Roca AL, Antunes A, Menotti-Raymond M, Yuhki N, Pecon-Slattery J, Johnson WE, Bourque G, Tesler G, NISC Comparative Sequencing Program, O'Brien SJ: Initial sequence and comparative analysis of the cat. Genome Res. 2007, 17: 1675-1689. 10.1101/gr.6380007.PubMedPubMed CentralView ArticleGoogle Scholar
- Margulies EH, Vinson JP, NISC Comparative Sequencing Program, Miller W, Jaffe DB, Lindblad-Toh K, Chang JL, Green ED, Lander ES, Mullikin JC, Clamp M: An initial strategy for the systematic identification of functional elements in the human genome by low-redundancy comparative sequencing. Proc Natl Acad Sci USA. 2005, 102: 4795-4800. 10.1073/pnas.0409882102.PubMedPubMed CentralView ArticleGoogle Scholar
- Gregory TR, Nicol JA, Tamm H, Kullman K, Leitch IJ, Murray BG, Kapraun DF, Greilhuber J, Bennett MD: Eukaryotic genome size database. Nucleic Acids Res. 2007, 35: D332-D338. 10.1093/nar/gkl828.PubMedPubMed CentralView ArticleGoogle Scholar
- Lander ES, Waterman MS: Genomic mapping by fingerprinting random clones: a mathematical analysis. Genomics. 1988, 2: 231-239. 10.1016/0888-7543(88)90007-9.PubMedView ArticleGoogle Scholar
- Takahashi K, Okada N: Mosaic structure and retropositional dynamics during evolution of subfamilies of short interspersed elements in African cichlids. Mol Biol Evol. 2002, 19: 1303-1312.PubMedView ArticleGoogle Scholar
- Jordan IK, Rogozin IB, Glazko GV, Koonin EV: Origin of a substantial fraction of human regulatory sequences from transposable elements. Trends Genet. 2003, 19: 68-72. 10.1016/S0168-9525(02)00006-9.PubMedView ArticleGoogle Scholar
- Kotarski MA, Leonard DA, Bennett SA, Bishop CP, Wahn SD, Sedore SA, Shrader M: The Drosophila gene asteroid encodes a novel protein and displays dosage-sensitive interactions with Star and Egfr. Genome. 1998, 41: 295-302. 10.1139/gen-41-2-295.PubMedView ArticleGoogle Scholar
- Georgia Tech Streelman Lab: Online Cichlid Resources. [http://cichlids.biology.gatech.edu]
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, Nicaud S, Jaffe D, Fisher S, Lutfalla G, Dossat C, Segurens B, Dasilva C, Salanoubat M, Levy M, Boudet N, Castellano S, Anthouard V, Jubin C, Castelli V, Katinka M, Vacherie B, Biémont C, Skalli Z, Cattolico L, Poulain J, et al: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431: 946-957. 10.1038/nature03025.PubMedView ArticleGoogle Scholar
- Turner GF, Seehausen O, Knight ME, Allender CF, Robinson RL: How many species of cichlid fishes are there in African lakes?. Mol Ecol. 2001, 10: 793-806. 10.1046/j.1365-294x.2001.01200.x.PubMedView ArticleGoogle Scholar
- Danley PD, Markert JA, Arnegard ME, Kocher TD: Divergence with gene flow in the rock-dwelling cichlids of Lake Malawi. Evolution. 2000, 54: 1725-1737.PubMedView ArticleGoogle Scholar
- Streelman JT, Peichel CL, Parichy DM: Developmental genetics of adaptation in fishes: the case for novelty. Ann Rev Ecol Evol Syst. 2007, 38: 655-681. 10.1146/annurev.ecolsys.38.091206.095537.View ArticleGoogle Scholar
- Arnegard ME, Markert JA, Danley PD, Stauffer JR, Ambali AJ, Kocher TD: Population structure and colour variation of the cichlid fish Labeotropheus fuelleborni along a recently formed archipelago of rocky habitat patches in southern Lake Malawi. Proc Biol Sci. 1999, 266: 119-130. 10.1098/rspb.1999.0611.PubMed CentralView ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
- Kocher TD, Conroy JA, McKaye KR, Stauffer JR, Lockwood SF: Evolution of NADH dehydrogenase subunit 2 in East African cichlid fish. Mol Phylogenet Evol. 1995, 4: 420-432. 10.1006/mpev.1995.1039.PubMedView ArticleGoogle Scholar
- Seehausen O, Koetsier E, Schneider MV, Chapman LJ, Chapman CA, Knight ME, Turner GF, van Alphen JJM, Bills R: Nuclear markers reveal unexpected genetic variation and a Congolese-Nilotic origin of the Lake Victoria cichlid species flock. Proc Biol Sci. 2003, 270: 129-137. 10.1098/rspb.2002.2153.PubMedPubMed CentralView ArticleGoogle Scholar
- Streelman JT, Gmyrek SL, Kidd MR, Kidd CE, Robinson RL, Hert E, Ambali AJ, Kocher TD: Hybridization and contemporary evolution in an introduced cichlid fish from Lake Malawi National Park. Mol Ecol. 2004, 13: 2471-2479. 10.1111/j.1365-294X.2004.02240.x.PubMedView ArticleGoogle Scholar
- Albertson RC, Markert JA, Danley PD, Kocher TD: Phylogeny of a rapidly evolving clade: the cichlid fishes of Lake Malawi, East Africa. Proc Natl Acad Sci USA. 1999, 96: 5107-5110. 10.1073/pnas.96.9.5107.PubMedPubMed CentralView ArticleGoogle Scholar
- Kidd MR, Kidd CE, Kocher TD: Axes of differentiation in the bower-building cichlids of Lake Malawi. Mol Ecol. 2006, 15: 459-478. 10.1111/j.1365-294X.2005.02787.x.PubMedView ArticleGoogle Scholar
- The Cichlid Fishes of Lake Malawi, Africa. [http://malawicichlids.com]
- Guryev V, Koudijs MJ, Berezikov E, Johnson SL, Plasterk RHA, van Eeden FJM, Cuppen E: Genetic variation in the zebrafish. Genome Res. 2006, 16: 491-497. 10.1101/gr.4791006.PubMedPubMed CentralView ArticleGoogle Scholar
- Terai Y, Morikawa N, Okada N: The evolution of the pro-domain of bone morphogenetic protein 4 (Bmp4) in an explosively speciated lineage of East African cichlid fishes. Mol Biol Evol. 2002, 19: 1628-1632.PubMedView ArticleGoogle Scholar
- Carleton KL, Parry JW, Bowmaker JK, Hunt DM, Seehausen O: Color vision and speciation in Lake Victoria cichlids of the genus Pundamilia. Mol Ecol. 2005, 14: 4341-4353.PubMedView ArticleGoogle Scholar
- Miyasaka KY, Kida YS, Sato T, Minami M, Ogura T: Csrp1 regulates dynamic cell movements of the mesendoderm and cardiac mesoderm through interactions with Dishevelled and Diversin. Proc Natl Acad Sci USA. 2007, 104: 11274-11279. 10.1073/pnas.0702000104.PubMedPubMed CentralView ArticleGoogle Scholar
- Chenn A, Walsh CA: Regulation of cerebral cortical size by control of cell cycle exit in neural precursors. Science. 2002, 297: 365-369. 10.1126/science.1074192.PubMedView ArticleGoogle Scholar
- Koudijs MJ, den Broeder MJ, Groot E, van Eeden GF: Genetic analysis of the two zebrafish patched homologues identifies novel roles for the hedgehog signaling pathway. BMC Dev Biol. 2008, 8: 15-10.1186/1471-213X-8-15.PubMedPubMed CentralView ArticleGoogle Scholar
- Scholpp S, Foucher I, Staudt N, Peukert D, Lumsden A, Houart C: Otx1l, Otx2 and Irx1b establish and position the ZLI in the diencephalon. Development. 2007, 134: 3167-3176. 10.1242/dev.001461.PubMedView ArticleGoogle Scholar
- Ye H, Tan YL, Ponniah S, Takeda Y, Wang SQ, Schachner M, Watanabe K, Pallen CJ, Xiao ZC: Neural recognition molecules CHL1 and NB-3 regulate apical dentrite orientation in the neocortex via PTP alpha. EMBO J. 2008, 27: 188-200. 10.1038/sj.emboj.7601939.PubMedPubMed CentralView ArticleGoogle Scholar
- Reinthal PN: The feeding habits of a group of herbivorous rock-dwelling fishes from Lake Malawi, Africa. Env Biol Fishes. 1990, 27: 215-233. 10.1007/BF00001674.View ArticleGoogle Scholar
- Fraser GJ, Bloomquist RF, Streelman JT: A periodic pattern generator for dental diversity. BMC Biology. 2008, 6: 32-10.1186/1741-7007-6-32.PubMedPubMed CentralView ArticleGoogle Scholar
- Green P: 2 × genomes - does depth matter?. Genome Res. 2007, 17: 1547-1549. 10.1101/gr.7050807.PubMedView ArticleGoogle Scholar
- Pavy N, Parsons LS, Paule C, MacKay J, Bousquet J: Automated SNP detection from a large collection of white spruce expressed sequences: contributing factors and approaches for the categorization of SNPs. BMC Genomics. 2006, 7: 174-10.1186/1471-2164-7-174.PubMedPubMed CentralView ArticleGoogle Scholar
- Moon S, Shin HD, Cheong HS, Cho HY, Namgoong S, Kim EM, Han CS, Kim H: BcSNPdb: Bovine coding region single nucleotide polymorphisms located proximal to quantitative trait loci. J Biochem Mol Biol. 2007, 40: 95-99.PubMedView ArticleGoogle Scholar
- Ewing B, Hiller L, Wendl M, Green P: Basecalling of automated sequencer traces using Phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.PubMedView ArticleGoogle Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467. 10.1159/000084979.PubMedView ArticleGoogle Scholar
- Watterson GA: On the number of segregating sites in genetic models without recombination. Theor Pop Biol. 1975, 7: 256-276. 10.1016/0040-5809(75)90020-9.View ArticleGoogle Scholar
- Jukes TH, Cantor CR: Evolution of protein molecules. Mammalian Protein Metabolism. Edited by: Munro HN. 1969, New York, NY: Academic Press, 21-132.View ArticleGoogle Scholar
- Bioperl. [http://www.bioperl.org]
- Beckman Coulter Autoprimer.com. [http://www.autoprimer.com]
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.Google Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.