STRUCTURAL EVOLUTION OF AFRICAN CICHLID GENOMES

An unanswered question in biology is how the evolution of genome structure supports or accompanies diversification and speciation on different time scales. African cichlid fishes are a well-documented system ideal for studying rapid evolution, due to their phenotypic diversity and high number of speciation events over the last several million years. I generated two de novo genome assemblies of the riverine cichlid Oreochromis niloticus (tilapia) and the Lake Malawi cichlid Metriaclima zebra using high-coverage long-read sequencing data and anchored the assemblies to chromosomes using several genetic and physical maps, to produce two high-quality anchored references. By comparing these chromosome-scale assemblies to integrated recombination, transcriptome, and resequencing data of multiple genera and species, I identified and characterized many large novel genome rearrangement events. These rearrangements included multiple novel sex-determination inversions, several metacentric-acrocentric karyotype differences via centromere assembly and placement, and wide regions of suppressed recombination in genera- and species-level crosses of Lake Malawi cichlids. Karyotype evolution in cichlids was further analyzed with long-read sequencing, specifically revealing the complex structure and content of a highly repetitive supernumerary chromosome present in some but not all individuals of a population across a wide range of eukaryotes, including many cichlid species. These supernumerary "B" chromosomes are shown to be limited to female Lake Malawi cichlids and have a unique evolutionary history with B chromosomes present in Lake Victorian cichlids male and females. This work reveals how structural genomic changes impact a rapidly evolving clade, while providing high-quality resources for the community, a context for previous genetic studies, and a robust platform for future genome research in cichlids

An improved genome reference for the African cichlid, Metriaclima zebra

Funding for Open Access provided by the UMD Libraries Open Access Publishing Fund.Background: Problems associated with using draft genome assemblies are well documented and have become more pronounced with the use of short read data for de novo genome assembly. We set out to improve the draft genome assembly of the African cichlid fish, Metriaclima zebra, using a set of Pacific Biosciences SMRT sequencing reads corresponding to 16.5x coverage of the genome. Here we characterize the improvements that these long reads allowed us to make to the state-of-the-art draft genome previously assembled from short read data. Results: Our new assembly closed 68 % of the existing gaps and added 90.6Mbp of new non-gap sequence to the existing draft assembly of M. zebra. Comparison of the new assembly to the sequence of several bacterial artificial chromosome clones confirmed the accuracy of the new assembly. The closure of sequence gaps revealed thousands of new exons, allowing significant improvement in gene models. We corrected one known misassembly, and identified and fixed other likely misassemblies. 63.5 Mbp (70 %) of the new sequence was classified as repetitive and the new sequence allowed for the assembly of many more transposable elements. Conclusions: Our improvements to the M. zebra draft genome suggest that a reasonable investment in long reads could greatly improve many comparable vertebrate draft genome assemblies.https://doi.org/10.1186/s12864-015-1930-

Structure and decay of a proto-Y region in Tilapia, Oreochromis niloticus

Funding for Open Access provided by the UMD Libraries Open Access Publishing Fund.Sex-determination genes drive the evolution of adjacent chromosomal regions. Sexually antagonistic selection favors the accumulation of inversions that reduce recombination in regions adjacent to the sex-determination gene. Once established, the clonal inheritance of sex-linked inversions leads to the accumulation of deleterious alleles, repetitive elements and a gradual decay of sex-linked genes. This in turn creates selective pressures for the evolution of mechanisms that compensate for the unequal dosage of gene expression. Here we use whole genome sequencing to characterize the structure of a young sex chromosome and quantify sex-specific gene expression in the developing gonad. We found an 8.8 Mb block of strong differentiation between males and females that corresponds to the location of a previously mapped sex-determiner on linkage group 1 of Oreochromis niloticus. Putatively disruptive mutations are found in many of the genes within this region. We also found a significant female-bias in the expression of genes within the block of differentiation compared to those outside the block of differentiation. Eight candidate sex-determination genes were identified within this region. This study demonstrates a block of differentiation on linkage group 1, suggestive of an 8.8 Mb inversion encompassing the sex-determining locus. The enrichment of female-biased gene expression inside the proposed inversion suggests incomplete dosage compensation. This study helps establish a model for studying the early-to-intermediate stages of sex chromosome evolution.https://doi.org/10.1186/1471-2164-15-97

Comparative analysis of a sex chromosome from the blackchin tilapia, Sarotherodon melanotheron

Background Inversions and other structural polymorphisms often reduce the rate of recombination between sex chromosomes, making it impossible to fine map sex-determination loci using traditional genetic mapping techniques. Here we compare distantly related species of tilapia that each segregate an XY system of sex-determination on linkage group 1. We use whole genome sequencing to identify shared sex-patterned polymorphisms, which are candidates for the ancestral sex-determination mutation. Results We found that Sarotherodon melanotheron segregates an XY system on LG1 in the same region identified in Oreochromis niloticus. Both species have higher densities of sex-patterned SNPs, as well as elevated number of ancestral copy number variants in this region when compared to the rest of the genome, but the pattern of differentiation along LG1 differs between species. The number of sex-patterned SNPs shared by the two species is small, but larger than expected by chance, suggesting that a novel Y-chromosome arose just before the divergence of the two species. We identified a shared sex-patterned SNP that alters a Gata4 binding site near Wilms tumor protein that might be responsible for sex-determination. Conclusions Shared sex-patterned SNPs, insertions and deletions suggest an ancestral sex-determination system that is common to both S. melanotheron and O. niloticus. Functional analyses are needed to evaluate shared SNPs near candidate genes that might play a role in sex-determination of these species. Interspecific variation in the sex chromosomes of tilapia species provides an excellent model system for understanding the evolution of vertebrate sex chromosomes. (Résumé d'auteur

Mapping of pigmentation QTL on an anchored genome assembly of the cichlid fish, Metriaclima zebra

Pigmentation patterns are one of the most recognizable phenotypes across the animal kingdom. They play an important role in camouflage, communication, mate recognition and mate choice. Most progress on understanding the genetics of pigmentation has been achieved via mutational analysis, with relatively little work done to understand variation in natural populations. Pigment patterns vary dramatically among species of cichlid fish from Lake Malawi, and are thought to be important in speciation. In this study, we crossed two species, Metriaclima zebra and M. mbenjii, that differ in several aspects of their body and fin color. We genotyped 798 SNPs in 160 F2 male individuals to construct a linkage map that was used to identify quantitative trait loci (QTL) associated with the pigmentation traits of interest. We also used the linkage map to anchor portions of the M. zebra genome assembly. We constructed a linkage map consisting of 834 markers in 22 linkage groups that spanned over 1,933 cM. QTL analysis detected one QTL each for dorsal fin xanthophores, caudal fin xanthophores, and pelvic fin melanophores. Dorsal fin and caudal fin xanthophores share a QTL on LG12, while pelvic fin melanophores have a QTL on LG11. We used the mapped markers to anchor 66.5% of the M. zebra genome assembly. Within each QTL interval we identified several candidate genes that might play a role in pigment cell development. This is one of a few studies to identify QTL for natural variation in fish pigmentation. The QTL intervals we identified did not contain any pigmentation genes previously identified by mutagenesis studies in other species. We expect that further work on these intervals will identify new genes involved in pigment cell development in natural populations

An EST resource for tilapia based on 17 normalized libraries and assembly of 116,899 sequence tags

<p>Abstract</p> <p>Background</p> <p>Large collections of expressed sequence tags (ESTs) are a fundamental resource for analysis of gene expression and annotation of genome sequences. We generated 116,899 ESTs from 17 normalized and two non-normalized cDNA libraries representing 16 tissues from tilapia, a cichlid fish widely used in aquaculture and biological research.</p> <p>Results</p> <p>The ESTs were assembled into 20,190 contigs and 36,028 singletons for a total of 56,218 unique sequences and a total assembled length of 35,168,415 bp. Over the whole project, a unique sequence was discovered for every 2.079 sequence reads. 17,722 (31.5%) of these unique sequences had significant BLAST hits (e-value < 10^-10) to the UniProt database.</p> <p>Conclusion</p> <p>Normalization of the cDNA pools with double-stranded nuclease allowed us to efficiently sequence a large collection of ESTs. These sequences are an important resource for studies of gene expression, comparative mapping and annotation of the forthcoming tilapia genome sequence.</p

Modelling the effects of past and future climate on the risk of bluetongue emergence in Europe

Vector-borne diseases are among those most sensitive to climate because the ecology of vectors and the development rate of pathogens within them are highly dependent on environmental conditions. Bluetongue (BT), a recently emerged arboviral disease of ruminants in Europe, is often cited as an illustration of climate's impact on disease emergence, although no study has yet tested this association. Here, we develop a framework to quantitatively evaluate the effects of climate on BT's emergence in Europe by integrating high-resolution climate observations and model simulations within a mechanistic model of BT transmission risk. We demonstrate that a climate-driven model explains, in both space and time, many aspects of BT's recent emergence and spread, including the 2006 BT outbreak in northwest Europe which occurred in the year of highest projected risk since at least 1960. Furthermore, the model provides mechanistic insight into BT's emergence, suggesting that the drivers of emergence across Europe differ between the South and the North. Driven by simulated future climate from an ensemble of 11 regional climate models, the model projects increase in the future risk of BT emergence across most of Europe with uncertainty in rate but not in trend. The framework described here is adaptable and applicable to other diseases, where the link between climate and disease transmission risk can be quantified, permitting the evaluation of scale and uncertainty in climate change's impact on the future of such diseases

Transcriptome characterization via 454 pyrosequencing of the annelid Pristina leidyi, an emerging model for studying the evolution of regeneration

Background: The naid annelids contain a number of species that vary in their ability to regenerate lost body parts, making them excellent candidates for evolution of regeneration studies. However, scant sequence data exists to facilitate such studies. We constructed a cDNA library from the naid Pristina leidyi, a species that is highly regenerative and also reproduces asexually by fission, using material from a range of regeneration and fission stages for our library. We then sequenced the transcriptome of P. leidyi using 454 technology. Results: 454 sequencing produced 1,550,174 reads with an average read length of 376 nucleotides. Assembly of 454 sequence reads resulted in 64,522 isogroups and 46,679 singletons for a total of 111,201 unigenes in this transcriptome. We estimate that over 95% of the transcripts in our library are present in our transcriptome. 17.7% of isogroups had significant BLAST hits to the UniProt database and these include putative homologs of a number of genes relevant to regeneration research. Although many sequences are incomplete, the mean sequence length of transcripts (isotigs) is 707 nucleotides. Thus, many sequences are large enough to be immediately useful for downstream applications such as gene expression analyses. Using in situ hybridization, we show that two Wnt/β-catenin pathway genes (homologs of frizzled and β-catenin) present in our transcriptome are expressed in the regeneration blastema of P. leidyi, demonstrating the usefulness of this resource for regeneration research. Conclusions: 454 sequencing is a rapid and efficient approach for identifying large numbers of genes in an organism that lacks a sequenced genome. This transcriptome dataset will be a valuable resource for molecular analyses of regeneration in P. leidyi and will serve as a starting point for comparisons to non-regenerating naids. It also contributes significantly to the still limited genomic resources available for annelids and lophotrochozoans more generally.https://doi.org/10.1186/1471-2164-13-28