<title>Background</title> <p>Cichlid fishes from the East African Rift lakes Victoria, Tanganyika, and Malawi represent a preeminent example of replicated and rapid evolutionary radiation [##REF##15131652##1##]. 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 [##REF##15851665##2##, ####REF##16461358##3##, ##REF##17519189##4####17519189##4##]. Recently evolved cichlid species segregate ancestral polymorphism [##UREF##0##5##,##REF##9826684##6##] and may exchange genes [##REF##12919487##7##,##REF##16701254##8##]. Numerous genomic resources have been developed for East African cichlids (many of which are summarized by the Cichlid Genome Consortium [##UREF##1##9##]). These include the following: genetic linkage maps for tilapia [##REF##12704237##10##, ####REF##9539437##11##, ##REF##12047227##12####12047227##12##] and Lake Malawi species [##REF##12704237##10##,##REF##16555305##13##]; fingerprinted bacterial artificial chromosome libraries [##UREF##2##14##]; expressed sequence tag sequences for Lake Tanganyika and Lake Victoria cichlids [##UREF##3##15##]; and first-generation microarrays [##REF##15858202##16##,##REF##15238158##17##]. Many studies have used these resources to study cichlid population genetics, molecular ecology, and phylogeny (for review [##UREF##4##18##,##UREF##5##19##]). Recent reports have capitalized on the diversity among East African cichlids to study the evolution and genetic basis of many traits, including behavior [##REF##17391260##20##], olfaction [##REF##17710134##21##], pigmentation [##REF##12919484##22##, ####REF##14614144##23##, ##REF##15716505##24####15716505##24##], vision [##REF##15772376##25##,##REF##16213819##26##], sex determination [##REF##15716505##24##,##REF##15100706##27##], the brain [##REF##9288416##28##], and craniofacial development [##REF##12704237##10##,##REF##16555305##13##,##REF##16251275##29##].</p> <p>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 [##REF##18347593##30##]. 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 [##REF##14512627##31##, ####REF##17407382##32##, ##REF##17975172##33##, ##REF##15778292##34####15778292##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?</p> <p>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.</p>
<title>Materials and methods</title> <title>Samples</title> <p>Individuals of <italic>Mchenga conophorus </italic>(MC), <italic>Labeotropheus fuelleborni </italic>(LF), <italic>Melanochromis auratus </italic>(MA), <italic>Maylandia zebra </italic>(MZ), and <italic>Rhamphochromis esox </italic>(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.</p> <title>Trace sequences</title> <p>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.</p> <title>Sequence preprocessing and assembly</title> <p>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 [##REF##9521921##66##]. 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 [##REF##16093699##67##]. 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 &lt;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.</p> <title>Similarity search and alignment</title> <p>Orthologous genomic contig pairs were first identified using reciprocal BLASTN similarity searches with a strict E-value cutoff of 10<sup>-100</sup>, 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.</p> <p>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 ##FIG##0##1##). Segregating sites were then identified from alignment positions with high quality bases (&gt;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 (θ<sub>w </sub>[##REF##1145509##68##]). 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).</p> <title>Protein-coding sequence identification</title> <p>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<sup>-10 </sup>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 <italic>Actinopterygii </italic>(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.</p> <title>Evolutionary sequence divergence among JGI species</title> <p>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 [##UREF##11##69##]. 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 [##UREF##12##70##] was then used to calculate the K<sub>a</sub>/K<sub>s </sub>values of pair-wise alignments.</p> <title>Genotyping and validation of SNPs</title> <p>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 <italic>Tetraodon </italic>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 [##UREF##13##71##]) 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.</p> <title>Genetic differentiation within and among lineages</title> <p>Locus-specific F<sub>ST </sub>[##UREF##14##72##] was calculated using FSTAT version 2.9.3.2 [##UREF##15##73##] 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 F<sub>ST </sub>values. F<sub>ST </sub>outliers exceed the sum of the upper quartile value and 1.5 times the interquartile range.</p> <title>Genomic assignment</title> <p>We used a Bayesian method (STRUCTURE v.2.2 [##REF##10835412##46##]) 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 [##REF##16461358##3##, ####REF##17519189##4##, ##UREF##0##5####0##5##,##REF##8747298##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.</p>
<title>Results</title> <title>Sequence assembly</title> <p>Trace sequences of five Lake Malawi cichlid species, namely <italic>Mchenga conophorus </italic>(MC; formerly genus <italic>Copadichromis</italic>), <italic>Labeotropheus fuelleborni </italic>(LF), <italic>Melanochromis auratus </italic>(MA), <italic>Maylandia zebra </italic>(MZ; formerly genus <italic>Metriaclima</italic>) and <italic>Rhamphochromis esox </italic>(RE), were downloaded from the GenBank Trace Archive and assembled into contiguous (contig) sequences. The average cichlid genome is 1.1 × 10<sup>9 </sup>bases [##REF##17090588##35##], so the traces represent a sequence coverage of 12-17% for each of the five species (see Additional data file 1). Through several quality filtering and assembly steps (see Materials and methods [below]), the resultant genomic assemblies of the five cichlid species yielded an average of 60,862 contigs with a mean length of 1,193 bases per contig. The total first-pass assembly sequence length for each species ranged from 68,238,634 bases (MA) to 79,168,277 bases (MZ), or about 7% of an average cichlid genome. Assembly statistics are shown in Table ##TAB##0##1##.</p> <p>We noted that these first-pass assemblies were 'over-assembled' by roughly a factor of 2 when compared with theoretical expectations [##REF##3294162##36##]. 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 [##REF##12140242##37##] and other organisms [##REF##12547512##38##], we have retained alignments of these 'high trace number' contigs and have marked them as such (see Additional data files 3 and 4).</p> <title>Gene content and coverage</title> <p>To establish the extent of gene content and coverage present in each assembly, we carried out BLASTX similarity searches (10<sup>-10 </sup>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.</p> <title>Clustering and alignment</title> <p>We obtained 25,458 clusters of putatively orthologous sequences, which were individually assembled into multi-species alignments for subsequent comparative analyses. Genic regions, as identified by similarity searches to known human and fish genes, were marked onto each alignment. Figure ##FIG##0##1## illustrates a typical example of one such alignment.</p> <p>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 <italic>Drosophila </italic>[##REF##9644838##39##]), 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).</p> <p>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.</p> <title>Segregating sites</title> <p>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 [##UREF##6##40##]). Nucleotide diversity (Watterson's θ<sub>w</sub>) 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.</p> <p>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).</p> <p>Finally, we calculated the ratio of replacement to synonymous substitutions (K<sub>a</sub>/K<sub>s</sub>) 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. K<sub>a</sub>/K<sub>s </sub>ranged from 0.380 in MC/LF to 0.562 in LF/MA. These numbers are greater than the ratios found between <italic>Fugu </italic>and <italic>Tetraodon </italic>(0.127 to 0.144 [##REF##15496914##41##]). Such high K<sub>a</sub>/K<sub>s </sub>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.</p> <title>Validation and generality of SNPs</title> <p>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 [##REF##16461358##3##,##REF##15772376##25##] and by us (Streelman JT, unpublished data), with expected variation in Malawi cichlids. Positive controls included genes involved in morphogenesis (<italic>otx1</italic>, <italic>otx2</italic>, and <italic>pax9</italic>), pigmentation (<italic>mitf</italic>, <italic>ednrb</italic>, and <italic>aim1</italic>), and visual sensitivity (opsins <italic>rh1</italic>, <italic>sws1</italic>, <italic>lws</italic>, <italic>sws2a</italic>, and <italic>sws2b</italic>). 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 &lt; e<sup>-100 </sup>to <italic>Tetraodon</italic>, 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 <italic>Tetraodon </italic>chromosome 11; we have previously shown <italic>Tetraodon </italic>11 to share orthologs with cichlid chromosome 5 [##REF##16555305##13##]. 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.</p> <p>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 [##REF##15851665##2##]), 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 [##REF##11298988##42##].</p> <p>Ten out of 13 (about 77%) positive controls gave reliable genotypes and were variable across the dataset. For the 59 SNPs predicted by our automated procedure, 11 were fixed (no variation) in all samples, indicating an error in sequencing (or genotyping), an error in prediction, or the presence of a low frequency allele in the sequenced samples. Six predicted SNPs did not produce data reliable enough for genotype calls. The remaining 42 loci from automated predictions (about 71%) were polymorphic across the dataset. For 24 SNPs predicted using manual similarity searches, four were fixed and four failed reliability for genotype calls, with the remaining 16 loci (about 67%) showing polymorphism (Table ##TAB##1##2##). Twelve out of 20 (60%) predicted SNPs with PQS of 3 or less were successful, whereas 30 out of 39 (76%) predictions with PQS of at least 4 yielded polymorphisms (Table ##TAB##2##3##). There is evidence of ascertainment bias in our genotypic data (see Additional data file 5). For example, three SNP loci (Aln100674, Aln114498, and Aln102321) exhibit alleles unique to <italic>Rhamphochromis</italic>. Similarly, SNPs predicted from comparisons of RE and mbuna (LF, MA, and MZ) are sometimes fixed in mbuna. Polymorphisms predicted from comparisons of mbuna taxa are more likely to vary within LF and MZ populations and across mbuna species.</p> <title>Genetic polymorphism and divergence at multiple scales</title> <p>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.</p> <p>The data (Additional data file 5) support the previously reported population structures in MZ [##REF##11108599##43##,##UREF##7##44##] and LF [##UREF##8##45##], as well as the genetic distinction between these species (MC Mims, unpublished data). For example, mean genetic differentiation (F<sub>ST</sub>) in MZ is 0.148 and in LF is 0.271. Mean F<sub>ST </sub>between LF and MZ was 0.215, and between mbuna (25 species) and non-mbuna (52 species) it was 0.224, demonstrating that most genetic variation segregates within and not between lineages, regardless of evolutionary scale. Nevertheless, these distributions of F<sub>ST </sub>yielded statistical outliers, which exhibit greater than average genetic differentiation (Figure ##FIG##1##2##). Four loci were found to be statistical outliers for F<sub>ST </sub>among MZ and LF populations. In MZ the opsin loci <italic>lws </italic>(F<sub>ST </sub>= 0.514), <italic>sws1 </italic>(0.572) and <italic>rh1 </italic>(0.733), and in LF the opsin locus <italic>rh1 </italic>(0.853) exhibit differentiation between populations. Between LF and MZ, three loci were identified as outliers: a nonsynonymous polymorphism in <italic>csrp1 </italic>(F<sub>ST </sub>= 0.893), a synonymous polymorphism in <italic>β-catenin </italic>(Aln101106_1089; F<sub>ST </sub>= 0.904), and an intronic polymorphism in <italic>ptc2 </italic>(Aln100281_1741; F<sub>ST </sub>= 0.863). Two statistical outliers were identified for F<sub>ST </sub>between rock-dwelling mbuna and non-mbuna groups: a nonsynonymous polymorphism in <italic>irx1 </italic>(Aln102504_1609; F<sub>ST </sub>= 0.984), and a nongenic polymorphism (Aln103534_280; F<sub>ST </sub>= 0.919) in sequence with similarity to pufferfish and stickleback genomes between <italic>contactin 3 </italic>and <italic>ncam L1</italic>.</p> <title>Genetic clustering and ancestry</title> <p>To further visualize the segregation of SNPs across the Malawi cichlid flock, we utilized a Bayesian approach that assigns individuals to a predefined number of genetic clusters [##REF##10835412##46##]. Specifically, we were interested in how species would be assigned to major Malawi cichlid lineages identified in previous studies [##REF##16461358##3##,##REF##17519189##4##,##REF##8747298##47##]. There are three such groups supported by the majority of molecular data: the rock-dwelling mbuna; pelagic and sand-dwelling species; and a group comprised of <italic>Rhamphochromis</italic>, <italic>Diplotaxodon</italic>, and other deep-water taxa. Analysis of 68 SNP loci accurately classifies species to respective lineages (Figure ##FIG##2##3##). For instance, all species considered mbuna (blue) cluster with other mbuna, to the exclusion of other groups; species thought to represent the earliest divergence within the species flock (<italic>Rhamphochromis</italic>) clustered together as a separate group (green); all remaining non-mbuna species formed the third group (red). Notably, deepwater genera <italic>Diplotaxodon </italic>and <italic>Pallidochromis </italic>contain individuals with mosaic genomes (red and green) and <italic>Astatotilapia calliptera</italic>, a nonendemic species and possible Malawi ancestor [##REF##12590750##48##] combines mbuna and non-mbuna genomes.</p> <p>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.</p>
<title>Discussion</title> <p>African cichlid fishes are important models of evolutionary diversification in form and function [##UREF##7##44##]. 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 [##REF##11298988##42##]. These recently formed species segregate ancestral polymorphism and exchange genes by hybridization [##UREF##0##5##,##REF##12919487##7##,##REF##15245419##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 [##REF##16555305##13##,##REF##12919484##22##,##REF##15716505##24##,##REF##16251275##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 [##REF##15851665##2##,##REF##16461358##3##]; even the promising approach of Amplified Fragment Length Polymorphism (AFLP) does not provide strong resolution of the relationships among genera [##REF##14614144##23##,##REF##12590750##48##,##REF##10220426##50##,##REF##16448413##51##]. The data we present here should provide new resources and perspectives for cichlid evolutionary genomics.</p> <title>Cichlid species exhibit genomic polymorphism</title> <p>Lake Malawi cichlid species sequenced by the JGI embody the phylogenetic, morphological, and behavioral diversity found within the assemblage. <italic>Rhamphochromis esox </italic>(RE) is a large (about 0.5 m) pelagic predator that represents one of the basal lineages of the species flock [##REF##16461358##3##,##REF##17519189##4##,##REF##8747298##47##]. <italic>Mchenga conophorus </italic>(MC) is a sand-dwelling species that breeds on leks, where males construct 'bowers' to attract females. <italic>Melanochromis auratus </italic>(MA), <italic>Maylandia zebra </italic>(MZ), and <italic>Labeotropheus fuelleborni </italic>(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 [##UREF##9##52##]).</p> <p>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 θ<sub>w </sub>= 0.26%) is less than that found among laboratory strains of the zebrafish <italic>Danio rerio </italic>(Watterson's θ<sub>w </sub>= 0.48% [##REF##16533913##53##]). Although overall nucleotide diversity is less than that observed in <italic>Danio</italic>, 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 [##REF##15772376##25##,##REF##12200490##54##], or a relaxation of purifying selection attributable to small effective population size. However, we should view this estimate of K<sub>a</sub>/K<sub>s </sub>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 K<sub>a</sub>/K<sub>s </sub>approach to identifying selection may be largely inappropriate for such young species where ancestral alleles segregate as polymorphisms.</p> <p>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 (F<sub>ST</sub>) 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 [##REF##15851665##2##], 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 <italic>Rhamphochromis</italic>. Notably, <italic>Astatotilapia calliptera</italic>, 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.</p> <title>Genomic polymorphism and the divergence of Malawi cichlids</title> <p>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 <italic>lws</italic>, <italic>rh1</italic>, and <italic>sws1 </italic>are differentiated among populations of LF and MZ, adding to reports that opsin polymorphisms are associated with population-specific color patterns or visual environments [##REF##16313597##55##].</p> <p>SNPs in <italic>csrp1</italic>, <italic>β-catenin</italic>, and <italic>ptc2 </italic>exhibit greater than expected differentiation between LF and MZ. <italic>Csrp1 </italic>(cysteine-rich protein) is a vertebrate LIM-domain family member acting in the noncanonical WNT pathway, expressed in gut, intestine, and cardiac mesoderm [##REF##17592114##56##]. <italic>β-catenin </italic>acts to transduce signals in the canonical WNT pathway [##REF##12130776##57##] 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 [##REF##18284698##58##]; both areexpressed in developing cichlid dentitions, jaws, and brains (Fraser GJ, Sylvester JB, Streelman JT, unpublished data). A SNP in <italic>irx1 </italic>nearly perfectly differentiates rock-dwelling mbuna from the remainder of the Malawi species flock. <italic>Irx1 </italic>acts to position the boundary between the telencephalon and the posterior forebrain [##REF##17670791##59##]. Finally, a SNP located between <italic>contactin 3 </italic>and <italic>ncam L1 </italic>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 [##REF##18046458##60##]. Taken together, differentiated loci are interesting in the context of cichlid diversification because they affect the phenotypes that vary among lineages: color and vision [##REF##15772376##25##,##REF##16213819##26##], guts [##UREF##10##61##], dentitions [##REF##16555305##13##,##REF##18625062##62##], jaws [##REF##12704237##10##,##REF##16251275##29##], and brains [##REF##9288416##28##].</p> <title>Discovery for evolutionary biology</title> <p>There are obvious challenges when attempting to extract information from low coverage genomic sequence, and also obvious payoffs [##REF##14512627##31##, ####REF##17407382##32##, ##REF##17975172##33##, ##REF##15778292##34####15778292##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 [##REF##17975171##63##]. 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.</p> <p>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 [##UREF##3##15##]. Our contribution roughly doubles the available data.</p> <p>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 [##REF##16824208##64##]), zebrafish (65% [##REF##16533913##53##]), and cow (43% [##REF##17244488##65##]). 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).</p> <p>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 [##UREF##7##44##].</p>
<p>This is an open access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <p>Low coverage survey sequencing shows that although Lake Malawi cichlids are phenotypically and behaviorally diverse, they appear genetically like a subdivided population.</p> <title>Background</title> <p>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.</p> <title>Results</title> <p>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 (F<sub>ST</sub>) within species, between species, and between major lineages. Outliers for F<sub>ST </sub>are candidate genes that may have experienced a history of natural selection in the Malawi lineage.</p> <title>Conclusion</title> <p>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.</p>
<title>Abbreviations</title> <p>BLAST, Basic Local Alignment Search Tool; F<sub>ST</sub>, genetic differentiation; JGI, Joint Genome Institute; K<sub>a</sub>/K<sub>s</sub>, ratio of replacement to synonymous substitutions; LF, <italic>Labeotropheus fuelleborni</italic>; MA, <italic>Melanochromis auratus</italic>; Mb, megabases; MC, <italic>Mchenga conophorus</italic>; MZ, <italic>Maylandia zebra</italic>; NCBI, National Center for Biotechnology Information; PQS, polymorphism quality score; RE, <italic>Rhamphochromis esox</italic>; SNP, single nucleotide polymorphism.</p> <title>Authors' contributions</title> <p>YHL, JTS, SVY, and TDK conceived the idea and designed the study. YHL, LSK, and MCM performed the research. YHL and JTS analyzed the data and drafted the manuscript. All authors read and approved the final manuscript.</p> <title>Additional data files</title> <p>The following additional data are available with the online version of this paper. Additional data file ##SUPPL##0##1## is a table of trace sequence statistics of five Lake Malawi cichlid species. Additional data file ##SUPPL##1##2## is a list of human gene homologs found in the five cichlid species. Additional data file ##SUPPL##2##3## is a list of alignments and polymorphic sites. Additional data file ##SUPPL##3##4## is a list of alignments with BLAST hits to fish and humans. Additional data file ##SUPPL##4##5## is a table of major allele frequencies for biallelic SNPs surveyed across Lake Malawi cichlid populations and species.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>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.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Alignment of a typical cluster of orthologous sequences. <bold>(a) </bold>Overall alignment of assembly contigs from three different cichlid species with alignment positions indicated. <bold>(b) </bold>Expanded detail of nucleotide alignment. Filled pink block shows the expanded alignment corresponding to dotted red box in panel a. Filled blue block shows the alignment of corresponding species' traces that made up the assembly sequences. Lower case nucleotides have base quality scores under 20. Dashes '-' represent sequence unavailability. Asterisks '*' represent gaps inserted into the sequences. Dots '·' represent identity in alignment. Cap '^' represents segregating site. Alignment positions shown after consensus sequence. Polymorphism quality score shown below A-G single nucleotide polymorphism site.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Box-and-whisker plots of F<sub>ST </sub>values. F<sub>ST </sub>values were calculated for the following: within MZ, within LF, LF versus MZ, and Mbuna versus non-Mbuna. Upper and lower box bounds represent 75th and 25th percentiles, respectively. The solid lines within boxes represent the median value. Whiskers mark the furthest points from the median that are not classified as outliers. Unfilled circles represent outliers that are more than 1.5 times the interquartile range higher than the upper box bound. F<sub>ST</sub>, genetic differentiation; LF, <italic>Labeotropheus fuelleborni</italic>; MA, <italic>Melanochromis auratus</italic>; Mb, megabases; MC, <italic>Mchenga conophorus</italic>; MZ, <italic>Maylandia zebra</italic>.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Bayesian assignment of Lake Malawi cichlids to different evolutionary lineages. We show the contribution to each individual genome (q, which ranges from 0% to 100%) from each of K = 3 predefined genetic clusters (blue, red, and green), for data derived from single nucleotide polymorphisms (SNPs) in Tables 2 and 3. Note that this method predefines the number but not the identity of genetic clusters. Species names are written once; multiple individuals from species are grouped together (for example, four individuals of <italic>Pseudotropheus crabro</italic>). Species considered mbuna (blue) cluster with other mbuna, to the exclusion of other groups; species thought to represent the earliest divergence within the species flock (<italic>Rhamphochromis</italic>) clustered together as a separate group (green); and all remaining non-mbuna species formed the third group (red).</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>First-pass genomic assembly statistics for five Lake Malawi cichlid species</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left">MC</td><td align="left">LF</td><td align="left">MA</td><td align="left">MZ</td><td align="left">RE</td></tr></thead><tbody><tr><td align="left">Total number of contigs in assembly</td><td align="left">61,923</td><td align="left">58,245</td><td align="left">63,297</td><td align="left">65,094</td><td align="left">55,751</td></tr><tr><td align="left">Total length (bases)</td><td align="left">73,425,564</td><td align="left">70,858,381</td><td align="left">68,238,634</td><td align="left">79,168,277</td><td align="left">71,295,074</td></tr><tr><td align="left">Genome coverage<sup>a </sup>(%)</td><td align="left">6.68</td><td align="left">6.44</td><td align="left">6.20</td><td align="left">7.20</td><td align="left">6.48</td></tr><tr><td align="left">Mean trace length (bases)</td><td align="left">1,055</td><td align="left">1,092</td><td align="left">991</td><td align="left">1,145</td><td align="left">1,153</td></tr><tr><td align="left">Shortest contig length (bases)</td><td align="left">50</td><td align="left">50</td><td align="left">50</td><td align="left">50</td><td align="left">50</td></tr><tr><td align="left">Longest contig length (bases)</td><td align="left">19,632</td><td align="left">17,437</td><td align="left">21,601</td><td align="left">15,371</td><td align="left">21,351</td></tr><tr><td align="left">Mean contig length (bases)</td><td align="left">1,186</td><td align="left">1,217</td><td align="left">1,078</td><td align="left">1,216</td><td align="left">1,279</td></tr><tr><td align="left">Q25 contig length (bases)</td><td align="left">759</td><td align="left">846</td><td align="left">783</td><td align="left">805</td><td align="left">934</td></tr><tr><td align="left">Q50 (median) contig length (bases)</td><td align="left">966</td><td align="left">1,063</td><td align="left">949</td><td align="left">1,163</td><td align="left">1,113</td></tr><tr><td align="left">Q75 contig length (bases)</td><td align="left">1,403</td><td align="left">1,355</td><td align="left">1,102</td><td align="left">1,417</td><td align="left">1,407</td></tr><tr><td colspan="6"><hr/></td></tr><tr><td align="left">Total genic length (bases)</td><td align="left">2,863,110 (3.9%)</td><td align="left">2,841,933 (4.0%)</td><td align="left">2,761,941 (4.0%)</td><td align="left">2,851,968 (3.6%)</td><td align="left">2,797,548 (3.9%)</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>SNP genotyping success categorized by detection method</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">SNP detection method</td><td align="left">Control genes</td><td align="left">Automated</td><td align="left">Manual BLAST</td></tr></thead><tbody><tr><td align="left">Number of genotyped loci</td><td align="left">13</td><td align="left">59</td><td align="left">24</td></tr><tr><td align="left">Number of polymorphic loci</td><td align="left">10</td><td align="left">42</td><td align="left">16</td></tr><tr><td align="left">Number of fixed loci</td><td align="left">3</td><td align="left">11</td><td align="left">4</td></tr><tr><td align="left">Number of failed loci</td><td align="left">0</td><td align="left">6</td><td align="left">4</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left">Successful SNP detection (%)</td><td align="left">76.9</td><td align="left">71.2</td><td align="left">66.7</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>SNP genotyping success categorized by polymorphic quality score</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Polymorphic quality score</td><td align="left">2</td><td align="left">3</td><td align="left">4</td><td align="left">5</td></tr></thead><tbody><tr><td align="left">Number of genotyped loci</td><td align="left">5</td><td align="left">15</td><td align="left">28</td><td align="left">11</td></tr><tr><td align="left">Number of polymorphic loci</td><td align="left">2</td><td align="left">10</td><td align="left">24</td><td align="left">6</td></tr><tr><td align="left">Number of fixed/failed loci</td><td align="left">3</td><td align="left">5</td><td align="left">4</td><td align="left">5</td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="left">Successful SNP detection (%)</td><td align="left">40</td><td align="left">66.7</td><td align="left">85.7</td><td align="left">54.5</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional data file 1</title><p>Presented is a table of trace sequence statistics of five Lake Malawi cichlid species.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional data file 2</title><p>Presented is a list of human gene homologs found in the five cichlid species.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional data file 3</title><p>Presented is a list of alignments and polymorphic sites.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional data file 4</title><p>Presented is a list of alignments with BLAST hits to fish and humans.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional data file 5</title><p>Presented is a table of major allele frequencies for biallelic SNPs surveyed across Lake Malawi cichlid populations and species.</p></caption></supplementary-material>
<table-wrap-foot><p><sup>a</sup>Using an average cichlid genome size of 1.1 × 10<sup>9 </sup>bases. LF, <italic>Labeotropheus fuelleborni</italic>; MA, <italic>Melanochromis auratus</italic>; MC, <italic>Mchenga conophorus</italic>; MZ, <italic>Maylandia zebra</italic>; RE, <italic>Rhamphochromis esox</italic>; Q25, 25<sup>th </sup>percentile; Q50, median or 50<sup>th </sup>percentile; Q75, 75<sup>th </sup>percentile.</p></table-wrap-foot> <table-wrap-foot><p>BLAST, Basic Local Alignment Search Tool; SNP, single nucleotide polymorphism.</p></table-wrap-foot> <table-wrap-foot><p>SNP, single nucleotide polymorphism.</p></table-wrap-foot>
<graphic xlink:href="gb-2008-9-7-r113-1"/> <graphic xlink:href="gb-2008-9-7-r113-2"/> <graphic xlink:href="gb-2008-9-7-r113-3"/>
<media xlink:href="gb-2008-9-7-r113-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-7-r113-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-7-r113-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-7-r113-S4.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-7-r113-S5.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Moran", "Kornfield"], "given-names": ["P", "I"], "article-title": ["Retention of ancestral polymorphism in the Mbuna species flock of Lake Malawi."], "source": ["Mol Biol Evol"], "year": ["1993"], "volume": ["10"], "fpage": ["1015"], "lpage": ["1029"]}, {"article-title": ["Cichlid Genome Consortium"]}, {"surname": ["Katigiri", "Kidd", "Tomasino", "Davis", "Wishon", "Stern", "Carleton", "Howe", "Kocher"], "given-names": ["T", "CE", "E", "JT", "C", "JE", "KL", "AE", "TD"], "article-title": ["A BAC-based physical map of the Nile tilapia genome 89."], "source": ["BMC Genomics"], "year": ["2005"], "volume": ["9"], "fpage": ["89"], "pub-id": ["10.1186/1471-2164-6-89"]}, {"article-title": ["The Gene Index Project"]}, {"surname": ["Kornfield", "Smith"], "given-names": ["I", "PF"], "article-title": ["African cichlid fishes: model systems for evolutionary biology."], "source": ["Ann Rev Ecol Evol Syst"], "year": ["2000"], "volume": ["31"], "fpage": ["163"], "lpage": ["196"]}, {"surname": ["Genner", "Turner"], "given-names": ["MJ", "GF"], "article-title": ["The mbuna cichlids of Lake Malawi: a model for rapid speciation and adaptive radiation."], "source": ["Fish Fisheries"], "year": ["2005"], "volume": ["6"], "fpage": ["1"], "lpage": ["34"], "pub-id": ["10.1111/j.1467-2679.2005.00173.x"]}, {"article-title": ["Georgia Tech Streelman Lab: Online Cichlid Resources"]}, {"surname": ["Streelman", "Peichel", "Parichy"], "given-names": ["JT", "CL", "DM"], "article-title": ["Developmental genetics of adaptation in fishes: the case for novelty."], "source": ["Ann Rev Ecol Evol Syst"], "year": ["2007"], "volume": ["38"], "fpage": ["655"], "lpage": ["681"], "pub-id": ["10.1146/annurev.ecolsys.38.091206.095537"]}, {"surname": ["Arnegard", "Markert", "Danley", "Stauffer", "Ambali", "Kocher"], "given-names": ["ME", "JA", "PD", "JR", "AJ", "TD"], "article-title": ["Population structure and colour variation of the cichlid fish "], "italic": ["Labeotropheus fuelleborni "], "source": ["Proc Biol Sci"], "year": ["1999"], "volume": ["266"], "fpage": ["119"], "lpage": ["130"], "pub-id": ["10.1098/rspb.1999.0611"]}, {"article-title": ["The Cichlid Fishes of Lake Malawi, Africa"]}, {"surname": ["Reinthal"], "given-names": ["PN"], "article-title": ["The feeding habits of a group of herbivorous rock-dwelling fishes from Lake Malawi, Africa."], "source": ["Env Biol Fishes"], "year": ["1990"], "volume": ["27"], "fpage": ["215"], "lpage": ["233"], "pub-id": ["10.1007/BF00001674"]}, {"surname": ["Jukes", "Cantor", "Munro HN"], "given-names": ["TH", "CR"], "article-title": ["Evolution of protein molecules."], "source": ["Mammalian Protein Metabolism"], "year": ["1969"], "publisher-name": ["New York, NY: Academic Press"], "fpage": ["21"], "lpage": ["132"]}, {"article-title": ["Bioperl"]}, {"article-title": ["Beckman Coulter Autoprimer.com"]}, {"surname": ["Weir", "Cockerham"], "given-names": ["BS", "CC"], "article-title": ["Estimating F-statistics for the analysis of population structure."], "source": ["Evolution"], "year": ["1984"], "volume": ["38"], "fpage": ["1358"], "lpage": ["1370"], "pub-id": ["10.2307/2408641"]}, {"surname": ["Goudet"], "given-names": ["J"], "article-title": ["FSTAT (Version 1.2): A computer program to calculate F-statistics."], "source": ["J Hered"], "year": ["1995"], "volume": ["86"], "fpage": ["485"], "lpage": ["486"]}]
Genome Biol. 2008 Jul 10; 9(7):R113