Chromosome evolution in African cichlid fish: contributions from the physical mapping of repeated DNAs

Evidence for selfing in a vertebrate from whole-genome sequencing

A growing number of recent genomic studies report asexual parthenogenetic reproduction in a wide range of taxa, including vertebrate species from the reptile, bird, and fish lineages. Yet, self-fertilization (selfing) has been recorded only in a single vertebrate, the mangrove killifish Kryptolebias marmoratus In cichlid fishes, sex determination is notably diverse and can be influenced by the environment, and sequential hermaphroditism has been reported for some species. Here, we present evidence for a case of facultative selfing in the cichlid fish Benitochromis nigrodorsalis, which is otherwise known as biparentally reproducing ovophilic mouthbrooder from Western Africa. Our laboratory observations revealed that a wild-caught individual produced repeatedly viable offspring in absence of a mating partner. By analyzing genome-wide single-nucleotide polymorphism (SNP) data, we compare that individual and two of its offspring to shed light on its reproductive mode. First, our results confirm uniparental reproduction. Second, overall heterozygosity is reduced in the offspring compared with outbred individuals. Retained maternal heterozygosity in the offspring is ∼51%, which is close to the theoretically expected value of a heterozygosity reduction of 50% by selfing. Heterozygosity patterns along individual chromosomes do not point to alternative parthenogenetic reproductive mechanisms like automixis by terminal or central fusion. Facultative selfing may represent an adaptive strategy ensuring reproduction when mating partners are absent and, hence, contribute to the cichlids' enormous evolutionary success.

Fast satellite DNA evolution in Nothobranchius annual killifishes

Satellite DNA (satDNA) is a rapidly evolving class of tandem repeats, with some monomers being involved in centromere organization and function. To identify repeats associated with (peri)centromeric regions, we investigated satDNA across Southern and Coastal clades of African annual killifishes of the genus Nothobranchius. Molecular cytogenetic and bioinformatic analyses revealed that two previously identified satellites, designated here as NkadSat01-77 and NfurSat01-348, are associated with (peri)centromeres only in one lineage of the Southern clade. NfurSat01-348 was, however, additionally detected outside centromeres in three members of the Coastal clade. We also identified a novel satDNA, NrubSat01-48, associated with (peri)centromeres in N. foerschi, N. guentheri, and N. rubripinnis. Our findings revealed fast turnover of satDNA associated with (peri)centromeres and different trends in their evolution in two clades of the genus Nothobranchius.

Chromosome diversity in Buthidae and Chactidae scorpions from Brazilian fauna: Diploid number and distribution of repetitive DNA sequences

In this work, we analyzed cytogenetically eight Chactidae and Buthidae, including the localization of repetitive DNA sequences. The chactids possess monocentric chromosomes and the highest diploid numbers (2n=50 in Brotheas amazonicus, 2n=36 in Chactopsis amazonica, 2n=30 in Neochactas sp.) when compared with buthids (2n=10 in Tityus bahiensis, 2n=14 in Tityus apiacas and Tityus metuendus, 2n=18 in Tityus aba, 2n=26 in Ischnotelson peruassu). The localization of rDNA genes and (TTAGG)n sequences exhibited a conserved pattern of two terminal/subterminal ribosomal cistrons and terminal telomere signals. However, the comparison between the data of C-banding, DAPI after FISH and Cot-DNA fraction indicated a variable quantity and distribution of these regions, as follow: (i) positive heterochromatin and Cot-DNA signals (B. amazonicus and I. peruassu), (ii) small blocks of heterochromatin with large Cot-DNA signals (T. metuendus), (iii) positive heterochromatic regions and absence of Cot-DNA signals (T. aba and T. apiacas), and (iv) negative heterochromatin and Cot-DNA signals (T. bahiensis). Therefore, our results revealed that there still is not a clear relation between quantity of heterochromatin and presence of monocentric or holocentric chromosomes and occurrence of chromosomal rearrangements, indicating that repetitive regions in scorpions must be analyzed using different cytogenetic approaches.

Environmental Cues and Mechanisms Underpinning Sex Change in Fish

Fishes are the only vertebrates that undergo sex change during their lifetime, but even within this group, a unique reproductive strategy is displayed by only 1.5% of the teleosts. This lability in alternating sexual fate is the result of the simultaneous suppression and activation of opposing male and female networks. Here, we provide a brief review summarizing recent advances in our understanding of the environmental cues that trigger sex change and their perception, integration, and translation into molecular cascades that convert the sex of an individual. We particularly focus on molecular events underpinning the complex behavioral and morphological transformation involved in sex change, dissecting the main molecular players and regulatory networks that shape the transformation of one sex into the opposite. We show that histological changes and molecular pathways governing gonadal reorganization are better described than the neuroendocrine basis of sex change and that, despite important advances, information is lacking for the majority of hermaphrodite species. We highlight significant gaps in our knowledge of how sex change takes place and suggest future research directions.

Origin of a Giant Sex Chromosome

Chromosome size and morphology vary within and among species, but little is known about the proximate or ultimate causes of these differences. Cichlid fish species in the tribe Oreochromini share an unusual giant chromosome that is ∼3 times longer than the other chromosomes. This giant chromosome functions as a sex chromosome in some of these species. We test two hypotheses of how this giant sex chromosome may have evolved. The first hypothesis proposes that it evolved by accumulating repetitive elements as recombination was reduced around a dominant sex determination locus, as suggested by canonical models of sex chromosome evolution. An alternative hypothesis is that the giant sex chromosome originated via the fusion of an autosome with a highly repetitive B chromosome, one of which carried a sex determination locus. We test these hypotheses using comparative analysis of chromosome-scale cichlid and teleost genomes. We find that the giant sex chromosome consists of three distinct regions based on patterns of recombination, gene and transposable element content, and synteny to the ancestral autosome. The WZ sex determination locus encompasses the last ∼105 Mb of the 134-Mb giant chromosome. The last 47 Mb of the giant chromosome shares no obvious homology to any ancestral chromosome. Comparisons across 69 teleost genomes reveal that the giant sex chromosome contains unparalleled amounts of endogenous retroviral elements, immunoglobulin genes, and long noncoding RNAs. The results favor the B chromosome fusion hypothesis for the origin of the giant chromosome.

Genetic Variation and Hybridization in Evolutionary Radiations of Cichlid Fishes

Evolutionary radiations are responsible for much of the variation in biodiversity across taxa. Cichlid fishes are well known for spectacular evolutionary radiations, as they have repeatedly evolved into large and phenotypically diverse arrays of species. Cichlid genomes carry signatures of past events and, at the same time, are the substrate for ongoing evolution. We survey genome-wide data and the available literature covering 438 cichlid populations (412 species) across multiple radiations to synthesize information about patterns and sharing of genetic variation. Nucleotide diversity within species is low in cichlids, with 92% of surveyed populations having less diversity than the median value found in other vertebrates. Divergence within radiations is also low, and a large proportion of variation is shared among species due to incomplete lineage sorting and widespread hybridization. Population genetics therefore provides a suitable conceptual framework for evolutionary genomic studies of cichlid radiations. We focus in detail on the roles of hybridization in shaping the patterns of genetic variation and in promoting cichlid diversification.

High-quality chromosome-level genomes of two tilapia species reveal their evolution of repeat sequences and sex chromosomes

Tilapias are one of the most farmed fishes that are coined as "aquatic chicken" by the food industry. Nile tilapia and blue tilapia exhibit very recent transition of sex chromosome systems since their divergence approximately five million years ago, making them a great model for elucidating the molecular and evolutionary mechanisms of sex chromosome turnovers. Studies of their sex-determining pathways are also critical for developing genetic sex control in aquaculture. We report here the newly produced genomes of Nile tilapia and blue tilapia that integrate long-read sequencing and chromatin conformation data. The two nearly complete genomes have anchored over 97% of the sequences into linkage groups (LGs), and assembled majorities of complex repetitive regions including telomeres, centromeres and rDNA clusters. In particular, we inferred two episodes of repeat expansion at LG3 respectively in the ancestor of cichlids and that of tilapias. The consequential large heterochromatic region concentrated at one end of LG3 comprises tandem arrays of mRNA and small RNA genes, among which we have identified a candidate female determining gene Paics in blue tilapia. Paics shows female-specific patterns of single-nucleotide variants, copy numbers and expression patterns in gonads during early gonadogenesis. Our work provides a very important genomic resource for functional studies of cichlids, and suggested that unequal distribution of repeat content that impacts the local recombination rate might make some chromosomes more likely to become sex chromosomes.

Genomics to accelerate genetic improvement in tilapia

Selective breeding of tilapia populations started in the early 1990s and over the past three decades tilapia has become one of the most important farmed freshwater species, being produced in more than 125 countries around the globe. Although genome assemblies have been available since 2011, most of the tilapia industry still depends on classical selection techniques using mass spawning or pedigree information to select for growth traits with reported genetic gains of up to 20% per generation. The involvement of international breeding companies and research institutions has resulted in the rapid development and application of genomic resources in the last few years. GWAS and genomic selection are expected to contribute to uncovering the genetic variants involved in economically relevant traits and increasing the genetic gain in selective breeding programs, respectively. Developments over the next few years will probably focus on achieving a deep understanding of genetic architecture of complex traits, as well as accelerating genetic progress in the selection for growth-, quality- and robustness-related traits. Novel phenotyping technologies (i.e. phenomics), lower-cost whole-genome sequencing approaches, functional genomics and gene editing tools will be crucial in future developments for the improvement of tilapia aquaculture.

Adaptive Radiation from a Chromosomal Perspective: Evidence of Chromosome Set Stability in Cichlid Fishes (Cichlidae: Teleostei) from the Barombi Mbo Lake, Cameroon

Cichlid fishes are the subject of scientific interest because of their rapid adaptive radiation, resulting in extensive ecological and taxonomic diversity. In this study, we examined 11 morphologically distinct cichlid species endemic to Barombi Mbo, the largest crater lake in western Cameroon, namely Konia eisentrauti, Konia dikume, Myaka myaka, Pungu maclareni, Sarotherodon steinbachi, Sarotherodon lohbergeri, Sarotherodon linnellii, Sarotherodon caroli, Stomatepia mariae, Stomatepia pindu, and Stomatepia mongo. These species supposedly evolved via sympatric ecological speciation from a common ancestor, which colonized the lake no earlier than one million years ago. Here we present the first comparative cytogenetic analysis of cichlid species from Barombi Mbo Lake using both conventional (Giemsa staining, C-banding, and CMA₃/DAPI staining) and molecular (fluorescence in situ hybridization with telomeric, 5S, and 28S rDNA probes) methods. We observed stability on both macro and micro-chromosomal levels. The diploid chromosome number was 2n = 44, and the karyotype was invariably composed of three pairs of meta/submetacentric and 19 pairs of subtelo/acrocentric chromosomes in all analysed species, with the same numbers of rDNA clusters and distribution of heterochromatin. The results suggest the evolutionary stability of chromosomal set; therefore, the large-scale chromosomal rearrangements seem to be unlikely associated with the sympatric speciation in Barombi Mbo.

Extensive Chromosomal Reorganization in Apistogramma Fishes (Cichlidae, Cichlinae) Fits the Complex Evolutionary Diversification of the Genus

Neotropical cichlid fishes are one of the most diversified and evolutionarily successful species assemblages. Extremely similar forms and intraspecific polychromatism present challenges for the taxonomy of some of these groups. Several species complexes have a largely unknown origin and unresolved evolutionary processes. Dwarf cichlids of the genus Apistogramma, comprising more than a hundred species, exhibit intricate taxonomic and biogeographic patterns, with both allopatric and sympatric distributions. However, karyotype evolution and the role of chromosomal changes in Apistogramma are still unknown. In the present study, nine South American Apistogramma species were analyzed using conventional cytogenetic methods and the mapping of repetitive DNA sequences [18S rDNA, 5S rDNA, and (TTAGGG)n] by fluorescence in situ hybridization (FISH). Our results showed that Apistogramma has unique cytogenetic characteristics in relation to closely related groups, such as a reduced 2n and a large number of bi-armed chromosomes. Interspecific patterns revealed a scenario of remarkable karyotypic changes, including a reduction of 2n, the occurrence of B-chromosomes and evolutionary dynamic of rDNA tandem repeats. In addition to the well-known pre-zygotic reproductive isolation, the karyotype reorganization in the genus suggests that chromosomal changes could act as postzygotic barriers in areas where Apistogramma congeners overlap.

Chromosome-scale assemblies reveal the structural evolution of African cichlid genomes

BACKGROUND

African cichlid fishes are well known for their rapid radiations and are a model system for studying evolutionary processes. Here we compare multiple, high-quality, chromosome-scale genome assemblies to elucidate the genetic mechanisms underlying cichlid diversification and study how genome structure evolves in rapidly radiating lineages.

RESULTS

We re-anchored our recent assembly of the Nile tilapia (Oreochromis niloticus) genome using a new high-density genetic map. We also developed a new de novo genome assembly of the Lake Malawi cichlid, Metriaclima zebra, using high-coverage Pacific Biosciences sequencing, and anchored contigs to linkage groups (LGs) using 4 different genetic maps. These new anchored assemblies allow the first chromosome-scale comparisons of African cichlid genomes. Large intra-chromosomal structural differences (∼2-28 megabase pairs) among species are common, while inter-chromosomal differences are rare (<10 megabase pairs total). Placement of the centromeres within the chromosome-scale assemblies identifies large structural differences that explain many of the karyotype differences among species. Structural differences are also associated with unique patterns of recombination on sex chromosomes. Structural differences on LG9, LG11, and LG20 are associated with reduced recombination, indicative of inversions between the rock- and sand-dwelling clades of Lake Malawi cichlids. M. zebra has a larger number of recent transposable element insertions compared with O. niloticus, suggesting that several transposable element families have a higher rate of insertion in the haplochromine cichlid lineage.

CONCLUSION

This study identifies novel structural variation among East African cichlid genomes and provides a new set of genomic resources to support research on the mechanisms driving cichlid adaptation and speciation.

Insertion Hot Spots of DIRS1 Retrotransposon and Chromosomal Diversifications among the Antarctic Teleosts Nototheniidae

By their faculty to transpose, transposable elements are known to play a key role in eukaryote genomes, impacting both their structuration and remodeling. Their integration in targeted sites may lead to recombination mechanisms involved in chromosomal rearrangements. The Antarctic fish family Nototheniidae went through several waves of species radiations. It is a suitable model to study transposable element (TE)-mediated mechanisms associated to genome and chromosomal diversifications. After the characterization of Gypsy (GyNoto), Copia (CoNoto), and DIRS1 (YNoto) retrotransposons in the genomes of Nototheniidae (diversity, distribution, conservation), we focused on their chromosome location with an emphasis on the three identified nototheniid radiations (the Trematomus, the plunderfishes, and the icefishes). The strong intrafamily TE conservation and wide distribution across species of the whole family suggest an ancestral acquisition with potential secondary losses in some lineages. GyNoto and CoNoto (including Hydra and GalEa clades) mostly produced interspersed signals along chromosomal arms. On the contrary, insertion hot spots accumulating in localized regions (mainly next to centromeric and pericentromeric regions) highlighted the potential role of YNoto in chromosomal diversifications as facilitator of the fusions which occurred in many nototheniid lineages, but not of the fissions.

Development and Validation of 58K SNP-Array and High-Density Linkage Map in Nile Tilapia (O. niloticus)

Despite being the second most important aquaculture species in the world accounting for 7.4% of global production in 2015, tilapia aquaculture has lacked genomic tools like SNP-arrays and high-density linkage maps to improve selection accuracy and accelerate genetic progress. In this paper, we describe the development of a genotyping array containing more than 58,000 SNPs for Nile tilapia (Oreochromis niloticus). SNPs were identified from whole genome resequencing of 32 individuals from the commercial population of the Genomar strain, and were selected for the SNP-array based on polymorphic information content and physical distribution across the genome using the Orenil1.1 genome assembly as reference sequence. SNP-performance was evaluated by genotyping 4991 individuals, including 689 offspring belonging to 41 full-sib families, which revealed high-quality genotype data for 43,588 SNPs. A preliminary genetic linkage map was constructed using Lepmap2 which in turn was integrated with information from the O_niloticus_UMD1 genome assembly to produce an integrated physical and genetic linkage map comprising 40,186 SNPs distributed across 22 linkage groups (LGs). Around one-third of the LGs showed a different recombination rate between sexes, with the female being greater than the male map by a factor of 1.2 (1632.9 to 1359.6 cM, respectively), with most LGs displaying a sigmoid recombination profile. Finally, the sex-determining locus was mapped to position 40.53 cM on LG23, in the vicinity of the anti-Müllerian hormone (amh) gene. These new resources has the potential to greatly influence and improve the genetic gain when applying genomic selection and surpass the difficulties of efficient selection for invasively measured traits in Nile tilapia.

Unusual Diversity of Sex Chromosomes in African Cichlid Fishes

African cichlids display a remarkable assortment of jaw morphologies, pigmentation patterns, and mating behaviors. In addition to this previously documented diversity, recent studies have documented a rich diversity of sex chromosomes within these fishes. Here we review the known sex-determination network within vertebrates, and the extraordinary number of sex chromosomes systems segregating in African cichlids. We also propose a model for understanding the unusual number of sex chromosome systems within this clade.

Mobilization of retrotransposons as a cause of chromosomal diversification and rapid speciation: the case for the Antarctic teleost genus Trematomus

Background

The importance of transposable elements (TEs) in the genomic remodeling and chromosomal rearrangements that accompany lineage diversification in vertebrates remains the subject of debate. The major impediment to understanding the roles of TEs in genome evolution is the lack of comparative and integrative analyses on complete taxonomic groups. To help overcome this problem, we have focused on the Antarctic teleost genus Trematomus (Notothenioidei: Nototheniidae), as they experienced rapid speciation accompanied by dramatic chromosomal diversity. Here we apply a multi-strategy approach to determine the role of large-scale TE mobilization in chromosomal diversification within Trematomus species.

Results

Despite the extensive chromosomal rearrangements observed in Trematomus species, our measurements revealed strong interspecific genome size conservation. After identifying the DIRS1, Gypsy and Copia retrotransposon superfamilies in genomes of 13 nototheniid species, we evaluated their diversity, abundance (copy numbers) and chromosomal distribution. Four families of DIRS1, nine of Gypsy, and two of Copia were highly conserved in these genomes; DIRS1 being the most represented within Trematomus genomes. Fluorescence in situ hybridization mapping showed preferential accumulation of DIRS1 in centromeric and pericentromeric regions, both in Trematomus and other nototheniid species, but not in outgroups: species of the Sub-Antarctic notothenioid families Bovichtidae and Eleginopsidae, and the non-notothenioid family Percidae.

Conclusions

In contrast to the outgroups, High-Antarctic notothenioid species, including the genus Trematomus, were subjected to strong environmental stresses involving repeated bouts of warming above the freezing point of seawater and cooling to sub-zero temperatures on the Antarctic continental shelf during the past 40 millions of years (My). As a consequence of these repetitive environmental changes, including thermal shocks; a breakdown of epigenetic regulation that normally represses TE activity may have led to sequential waves of TE activation within their genomes. The predominance of DIRS1 in Trematomus species, their transposition mechanism, and their strategic location in “hot spots” of insertion on chromosomes are likely to have facilitated nonhomologous recombination, thereby increasing genomic rearrangements. The resulting centric and tandem fusions and fissions would favor the rapid lineage diversification, characteristic of the nototheniid adaptive radiation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4714-x) contains supplementary material, which is available to authorized users.

High-Throughput Sequencing and Linkage Mapping of a Clownfish Genome Provide Insights on the Distribution of Molecular Players Involved in Sex Change

Clownfishes are an excellent model system for investigating the genetic mechanism governing hermaphroditism and socially-controlled sex change in their natural environment because they are broadly distributed and strongly site-attached. Genomic tools, such as genetic linkage maps, allow fine-mapping of loci involved in molecular pathways underlying these reproductive processes. In this study, a high-density genetic map of Amphiprion bicinctus was constructed with 3146 RAD markers in a full-sib family organized in 24 robust linkage groups which correspond to the haploid chromosome number of the species. The length of the map was 4294.71 cM, with an average marker interval of 1.38 cM. The clownfish linkage map showed various levels of conserved synteny and collinearity with the genomes of Asian and European seabass, Nile tilapia and stickleback. The map provided a platform to investigate the genomic position of genes with differential expression during sex change in A. bicinctus. This study aims to bridge the gap of genome-scale information for this iconic group of species to facilitate the study of the main gene regulatory networks governing social sex change and gonadal restructuring in protandrous hermaphrodites.

A high quality assembly of the Nile Tilapia (Oreochromis niloticus) genome reveals the structure of two sex determination regions

Background

Tilapias are the second most farmed fishes in the world and a sustainable source of food. Like many other fish, tilapias are sexually dimorphic and sex is a commercially important trait in these fish. In this study, we developed a significantly improved assembly of the tilapia genome using the latest genome sequencing methods and show how it improves the characterization of two sex determination regions in two tilapia species.

Results

A homozygous clonal XX female Nile tilapia (Oreochromis niloticus) was sequenced to 44X coverage using Pacific Biosciences (PacBio) SMRT sequencing. Dozens of candidate de novo assemblies were generated and an optimal assembly (contig NG50 of 3.3Mbp) was selected using principal component analysis of likelihood scores calculated from several paired-end sequencing libraries. Comparison of the new assembly to the previous O. niloticus genome assembly reveals that recently duplicated portions of the genome are now well represented. The overall number of genes in the new assembly increased by 27.3%, including a 67% increase in pseudogenes. The new tilapia genome assembly correctly represents two recent vasa gene duplication events that have been verified with BAC sequencing. At total of 146Mbp of additional transposable element sequence are now assembled, a large proportion of which are recent insertions. Large centromeric satellite repeats are assembled and annotated in cichlid fish for the first time. Finally, the new assembly identifies the long-range structure of both a ~9Mbp XY sex determination region on LG1 in O. niloticus, and a ~50Mbp WZ sex determination region on LG3 in the related species O. aureus.

Conclusions

This study highlights the use of long read sequencing to correctly assemble recent duplications and to characterize repeat-filled regions of the genome. The study serves as an example of the need for high quality genome assemblies and provides a framework for identifying sex determining genes in tilapia and related fish species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3723-5) contains supplementary material, which is available to authorized users.

Chromosomal Mapping of Repetitive DNA Sequences in the Genus Bryconamericus (Characidae) and DNA Barcoding to Differentiate Populations

The mapping of repetitive DNA sites by fluorescence in situ hybridization has been widely used for karyotype studies in different species of fish, especially when dealing with related species or even genera presenting high chromosome variability. This study analyzed three populations of Bryconamericus, with diploid number preserved, but with different karyotype formulae. Bryconamericus ecai, from the Forquetinha river/RS, presented three new cytotypes, increasing the number of karyotype forms to seven in this population. Other two populations of Bryconamericus sp. from the Vermelho stream/PR and Cambuta river/PR exhibited interpopulation variation. The chromosome mapping of rDNA sites revealed unique markings among the three populations, showing inter- and intrapopulation variability located in the terminal region. The molecular analysis using DNA barcoding complementing the cytogenetic analysis also showed differentiation among the three populations. The U2 small nuclear DNA repetitive sequence exhibited conserved features, being located in the interstitial region of a single chromosome pair. This is the first report on its occurrence in the genus Bryconamericus. Data obtained revealed a karyotype variability already assigned to the genus, along with polymorphism of ribosomal sites, demonstrating that this group of fish can be undergoing a divergent evolutionary process, constituting a substantive model for studies of chromosomal evolution.

Identification of the sex-determining region in flathead grey mullet (Mugil cephalus)

Elucidation of the sex-determination mechanism in flathead grey mullet (Mugil cephalus) is required to exploit its economic potential by production of genetically determined monosex populations and application of hormonal treatment to parents rather than to the marketed progeny. Our objective was to construct a first-generation linkage map of the M. cephalus in order to identify the sex-determining region and sex-determination system. Deep-sequencing data of a single male was assembled and aligned to the genome of Nile tilapia (Oreochromis niloticus). A total 245 M. cephalus microsatellite markers were designed, spanning the syntenic tilapia genome assembly at intervals of 10 Mb. In the mapping family of full-sib progeny, 156 segregating markers were used to construct a first-generation linkage map of 24 linkage groups (LGs), corresponding to the number of chromosomes. The linkage map spanned approximately 1200 cM with an average inter-marker distance of 10.6 cM. Markers segregating on LG9 in two independent mapping families showed nearly complete concordance with gender (R²  = 0.95). The sex determining locus was fine mapped within an interval of 8.6 cM on LG9. The sex of offspring was determined only by the alleles transmitted from the father, thus indicating an XY sex-determination system.

Integrated cytogenetics and genomics analysis of transposable elements in the Nile tilapia, Oreochromis niloticus

Integration of cytogenetics and genomics has become essential to a better view of architecture and function of genomes. Although the advances on genomic sequencing have contributed to study genes and genomes, the repetitive DNA fraction of the genome is still enigmatic and poorly understood. Among repeated DNAs, transposable elements (TEs) are major components of eukaryotic chromatin and their investigation has been hindered even after the availability of whole sequenced genomes. The cytogenetic mapping of TEs in chromosomes has proved to be of high value to integrate information from the micro level of nucleotide sequence to a cytological view of chromosomes. Different TEs have been cytogenetically mapped in cichlids; however, neither details about their genomic arrangement nor appropriated copy number are well defined by these approaches. The current study integrates TEs distribution in Nile tilapia Oreochromis niloticus genome based on cytogenetic and genomics/bioinformatics approach. The results showed that some elements are not randomly distributed and that some are genomic dependent on each other. Moreover, we found extensive overlap between genomics and cytogenetics data and that tandem duplication may be the major mechanism responsible for the genomic dynamics of TEs here analyzed. This paper provides insights in the genomic organization of TEs under an integrated view based on cytogenetics and genomics.

A microsatellite-based genetic linkage map and putative sex-determining genomic regions in Lake Victoria cichlids

Cichlid fishes in East Africa have undergone extensive adaptive radiation, which has led to spectacular diversity in their morphology and ecology. To date, genetic linkage maps have been constructed for several tilapias (riverine), Astatotilapia burtoni (Lake Tanganyika), and hybrid lines of Lake Malawi cichlids to facilitate genome-wide comparative analyses. In the present study, we constructed a genetic linkage map of the hybrid line of Lake Victoria cichlids, so that maps of cichlids from all the major areas of East Africa will be available. The genetic linkage map shown here is derived from the F2 progeny of an interspecific cross between Haplochromis chilotes and Haplochromis sauvagei and is based on 184 microsatellite and two single-nucleotide polymorphism (SNP) markers. Most of the microsatellite markers used in the present study were originally designed for other genetic linkage maps, allowing us to directly compare each linkage group (LG) among different cichlid groups. We found 25 LGs, the total length of which was 1133.2cM with an average marker spacing of about 6.09cM. Our subsequent linkage mapping analysis identified two putative sex-determining loci in cichlids. Interestingly, one of these two loci is located on cichlid LG5, on which the female heterogametic ZW locus and several quantitative trait loci (QTLs) related to adaptive evolution have been reported in Lake Malawi cichlids. We also found that V1R1 and V1R2, candidate genes for the fish pheromone receptor, are located very close to the recently detected sex-determining locus on cichlid LG5. The genetic linkage map study presented here may provide a valuable foundation for studying the chromosomal evolution of East African cichlids and the possible role of sex chromosomes in generating their genomic diversity.

Comparative cytogenetics of Neotropical cichlid fishes (Nannacara, Ivanacara and Cleithracara) indicates evolutionary reduction of diploid chromosome numbers

A comparative cytogenetic analysis was carried out in five species of a monophyletic clade of neotropical Cichlasomatine cichlids, namely Cleithracara maronii Steindachner, 1881, Ivanacara adoketa (Kullander & Prada-Pedreros, 1993), Nannacara anomala Regan, 1905, N. aureocephalus Allgayer, 1983 and N. taenia Regan, 1912. Karyotypes and other chromosomal characteristics were revealed by CDD banding and mapped onto the phylogenetic hypothesis based on molecular analyses of four genes, namely cyt b, 16S rRNA, S7 and RAG1. The diploid numbers of chromosomes ranged from 44 to 50, karyotypes were composed predominantly of monoarmed chromosomes and one to three pairs of CMA3 signal were observed. The results showed evolutionary reduction in this monophyletic clade and the cytogenetic mechanisms (fissions/fusions) were hypothesized and discussed.

Differential expression of a retrotransposable element, Rex6, in Colossoma macropomum fish from different Amazonian environments

Transposable elements (TEs) are DNA sequences that have the ability to move and replicate within the genomes. TEs can be classified according to their intermediates of transposition, RNA (retrotransposons) or DNA. In some aquatic organisms, it has been observed that environmental factors such as pH, temperature and pollution may stimulate differential transcription and mobilization of retrotransposons. In light of this information, the present study sought to evaluate the expression of Rex6 TE transcripts in Colossoma macropomum, which is a very commercially exploited fish in Brazil. In order to establish a comparative analysis using real-time PCR, the samples were collected from Amazonian rivers with different physical and chemical characteristics (distinguished by clear water and black water). Quantitative RT-PCR analyses revealed a differential pattern of expression between tissues collected from different types of water (clear and black waters). When it came to the hepatic and muscle tissues sampled, the levels of Rex6 transcripts were significantly different between the two Amazonian water types. These results suggest that environmental conditions operate differently in the regulation of Rex6 transcription in C. macropomum, results which have implications in the reshaping of the genome against environmental variations.

Construction of chromosome markers from the Lake Victoria cichlid Paralabidochromis chilotes and their application to comparative mapping

Cichlid fishes in the African Great Lakes are known as a spectacular example of adaptive radiation in vertebrates. Four linkage maps have been constructed to identify the genes responsible for adaptation and speciation, and the genetic linkages of those genes are assumed to play an important role during adaptive evolution. However, it is difficult to analyze such linkages because the linkage groups of one species do not match well with those of the other species. Chromosome markers are a powerful tool for the direct identification of linkage homology between different species. We used information about the linkage map of the Lake Malawi cichlid (Labeotropheus fuelleborni/Metriaclima zebra) to isolate bacterial artificial chromosome (BAC) clones from the BAC library of Paralabidochromis chilotes, Lake Victoria. We identified 18 of 22 P. chilotes chromosomes by single- and multi-color BAC fluorescence in situ hybridization using 19 BAC clones. Comparative mapping with the chromosome markers of P. chilotes in Astatotilapia burtoni (2n = 40) from Lake Tanganyika revealed the chromosome rearrangements that have occurred in this lineage. These chromosome markers will be useful for delineating the process of genome and chromosome evolution in African species.

Evolutionary dynamics of retrotransposable elements Rex1, Rex3 and Rex6 in neotropical cichlid genomes

Background

Transposable elements (TEs) have the potential to produce broad changes in the genomes of their hosts, acting as a type of evolutionary toolbox and generating a collection of new regulatory and coding sequences. Several TE classes have been studied in Neotropical cichlids; however, the information gained from these studies is restricted to the physical chromosome mapping, whereas the genetic diversity of the TEs remains unknown. Therefore, the genomic organization of the non-LTR retrotransposons Rex1, Rex3, and Rex6 in five Amazonian cichlid species was evaluated using physical chromosome mapping and DNA sequencing to provide information about the role of TEs in the evolution of cichlid genomes.

Results

Physical mapping revealed abundant TE clusters dispersed throughout the chromosomes. Furthermore, several species showed conspicuous clusters accumulation in the centromeric and terminal portions of the chromosomes. These TE chromosomal sites are associated with both heterochromatic and euchromatic regions. A higher number of Rex1 clusters were observed among the derived species. The Rex1 and Rex3 nucleotide sequences were more conserved in the basal species than in the derived species; however, this pattern was not observed in Rex6. In addition, it was possible to observe conserved blocks corresponding to the reverse transcriptase fragment of the Rex1 and Rex3 clones and to the endonuclease of Rex6.

Conclusion

Our data showed no congruence between the Bayesian trees generated for Rex1, Rex3 and Rex6 of cichlid species and phylogenetic hypothesis described for the group. Rex1 and Rex3 nucleotide sequences were more conserved in the basal species whereas Rex6 exhibited high substitution rates in both basal and derived species. The distribution of Rex elements in cichlid genomes suggests that such elements are under the action of evolutionary mechanisms that lead to their accumulation in particular chromosome regions, mostly in heterochromatins.

Chromosomal mapping of microsatellite repeats in the rock bream fish Oplegnathus fasciatus, with emphasis of their distribution in the neo-Y chromosome

Despite the theoretical and experimental progress, our understanding on sex chromosome differentiation is still diagrammatic. The accumulation of repetitive DNA sequences is believed to occur in early stages of such differentiation. As fish species present a wide range of sex chromosome systems they are excellent models to examine the differentiation of these chromosomes. In the present study, the chromosomal distribution of 9 mono-, di- and tri-nucleotide microsatellites were analyzed using fluorescence in situ hybrization (FISH) in rock bream fish (Oplegnathus fasciatus), which is characterized by an X1X2Y sex chromosome system. Generally, the males and females exhibited the same autosomal pattern of distribution for a specific microsatellite probe. The male specific Y chromosome displays a specific amount of distinct microsatellites repeats along both arms. However, the accumulation of these repetitive sequences was not accompanied by a huge heterochromatinization process. The present data provide new insights into the chromosomal constitution of the multiple sex chromosomes and allow further investigations on the true role of the microsatellite repeats in the differentiation process of this sex system.

Evolutionary dynamics of rRNA gene clusters in cichlid fish

Background

Among multigene families, ribosomal RNA (rRNA) genes are the most frequently studied and have been explored as cytogenetic markers to study the evolutionary history of karyotypes among animals and plants. In this report, we applied cytogenetic and genomic methods to investigate the organization of rRNA genes among cichlid fishes. Cichlids are a group of fishes that are of increasing scientific interest due to their rapid and convergent adaptive radiation, which has led to extensive ecological diversity.

Results

The present paper reports the cytogenetic mapping of the 5S rRNA genes from 18 South American, 22 African and one Asian species and the 18S rRNA genes from 3 African species. The data obtained were comparatively analyzed with previously published information related to the mapping of rRNA genes in cichlids. The number of 5S rRNA clusters per diploid genome ranged from 2 to 15, with the most common pattern being the presence of 2 chromosomes bearing a 5S rDNA cluster. Regarding 18S rDNA mapping, the number of sites ranged from 2 to 6, with the most common pattern being the presence of 2 sites per diploid genome. Furthermore, searching the Oreochromis niloticus genome database led to the identification of a total of 59 copies of 5S rRNA and 38 copies of 18S rRNA genes that were distributed in several genomic scaffolds. The rRNA genes were frequently flanked by transposable elements (TEs) and spread throughout the genome, complementing the FISH analysis that detect only clustered copies of rRNA genes.

Conclusions

The organization of rRNA gene clusters seems to reflect their intense and particular evolutionary pathway and not the evolutionary history of the associated taxa. The possible role of TEs as one source of rRNA gene movement, that could generates the spreading of ribosomal clusters/copies, is discussed. The present paper reinforces the notion that the integration of cytogenetic data and genomic analysis provides a more complete picture for understanding the organization of repeated sequences in the genome.

Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish

BACKGROUND

The availability of a large number of recently sequenced vertebrate genomes opens new avenues to integrate cytogenetics and genomics in comparative and evolutionary studies. Cytogenetic mapping can offer alternative means to identify conserved synteny shared by distinct genomes and also to define genome regions that are still not fine characterized even after wide-ranging nucleotide sequence efforts. An efficient way to perform comparative cytogenetic mapping is based on BAC clones mapping by fluorescence in situ hybridization. In this report, to address the knowledge gap on the genome evolution in cichlid fishes, BAC clones of an Oreochromis niloticus library covering the linkage groups (LG) 1, 3, 5, and 7 were mapped onto the chromosomes of 9 African cichlid species. The cytogenetic mapping data were also integrated with BAC-end sequences information of O. niloticus and comparatively analyzed against the genome of other fish species and vertebrates.

RESULTS

The location of BACs from LG1, 3, 5, and 7 revealed a strong chromosomal conservation among the analyzed cichlid species genomes, which evidenced a synteny of the markers of each LG. Comparative in silico analysis also identified large genomic blocks that were conserved in distantly related fish groups and also in other vertebrates.

CONCLUSIONS

Although it has been suggested that fishes contain plastic genomes with high rates of chromosomal rearrangements and probably low rates of synteny conservation, our results evidence that large syntenic chromosome segments have been maintained conserved during evolution, at least for the considered markers. Additionally, our current cytogenetic mapping efforts integrated with genomic approaches conduct to a new perspective to address important questions involving chromosome evolution in fishes.

Genomic content and new insights on the origin of the B chromosome of the cichlid fish Astatotilapia latifasciata

B chromosomes are additional chromosomes widely studied in a diversity of eukaryotic groups, including fungi, plants and animals, but their origin, evolution and possible functions are not clearly understood. To further understand the genomic content and the evolutionary history of B chromosomes, classical and molecular cytogenetic analyses were conducted in the cichlid fish Astatotilapia latifasciata, which harbor 1–2 B chromosomes. Through cytogenetic mapping of several probes, including transposable elements, rRNA genes, a repeated DNA genomic fraction (C0t - 1 DNA), whole genome probes (comparative genomic hybridization), and BAC clones from Oreochromis niloticus, we found similarities between the B chromosome and the 1st chromosome pair and chromosomes harboring rRNA genes. Based on the cytogenetic mapping data, we suggest the B chromosome may have evolved from a small chromosomal fragment followed by the invasion of the proto-B chromosome by several repeated DNA families.

A high-resolution map of the Nile tilapia genome: a resource for studying cichlids and other percomorphs

Background

The Nile tilapia (Oreochromis niloticus) is the second most farmed fish species worldwide. It is also an important model for studies of fish physiology, particularly because of its broad tolerance to an array of environments. It is a good model to study evolutionary mechanisms in vertebrates, because of its close relationship to haplochromine cichlids, which have undergone rapid speciation in East Africa. The existing genomic resources for Nile tilapia include a genetic map, BAC end sequences and ESTs, but comparative genome analysis and maps of quantitative trait loci (QTL) are still limited.

Results

We have constructed a high-resolution radiation hybrid (RH) panel for the Nile tilapia and genotyped 1358 markers consisting of 850 genes, 82 markers corresponding to BAC end sequences, 154 microsatellites and 272 single nucleotide polymorphisms (SNPs). From these, 1296 markers could be associated in 81 RH groups, while 62 were not linked. The total size of the RH map is 34,084 cR₃₅₀₀ and 937,310 kb. It covers 88% of the entire genome with an estimated inter-marker distance of 742 Kb. Mapping of microsatellites enabled integration to the genetic map. We have merged LG8 and LG24 into a single linkage group, and confirmed that LG16-LG21 are also merged. The orientation and association of RH groups to each chromosome and LG was confirmed by chromosomal in situ hybridizations (FISH) of 55 BACs. Fifty RH groups were localized on the 22 chromosomes while 31 remained small orphan groups. Synteny relationships were determined between Nile tilapia, stickleback, medaka and pufferfish.

Conclusion

The RH map and associated FISH map provide a valuable gene-ordered resource for gene mapping and QTL studies. All genetic linkage groups with their corresponding RH groups now have a corresponding chromosome which can be identified in the karyotype. Placement of conserved segments indicated that multiple inter-chromosomal rearrangements have occurred between Nile tilapia and the other model fishes. These maps represent a valuable resource for organizing the forthcoming genome sequence of Nile tilapia, and provide a foundation for evolutionary studies of East African cichlid fishes.

Characterisation of the chromosome fusions in Oreochromis karongae

Oreochromis karongae, one of the "chambo" tilapia species from Lake Malawi, has a karyotype of 2n = 38, making it one of the few species investigated to differ from the typical tilapia karyotype (2n = 44). The O. karongae karyotype consists of one large subtelocentric pair of chromosomes, four medium-sized pairs (three subtelocentric and one submetacentric) and 14 small pairs. The five largest pairs could be distinguished from each other on the basis of size, morphology and a series of fluorescence in situ hybridisation (FISH) probes. The largest pair is easily distinguished on the basis of size and a chromosome 1 (linkage group 3) bacterial artificial chromosome (BAC) FISH probe from Oreochromis niloticus. BAC clones from O. niloticus chromosome 2 (linkage group 7) hybridised to one of the medium-sized subtelocentric chromosome pairs (no. 5) of O. karongae, distinguishing the ancestral medium-sized pair from the three other medium-sized chromosome pairs (nos. 2, 3 and 4) that appear to have resulted from fusions. SATA repetitive DNA hybridised to the centromeres of all 19 chromosome pairs and also revealed the locations of the relic centromeres in the three fused pairs. Telomeric (TTAGGG)(n) repeats were identified in the telomeres of all chromosomes, and an interstitial telomeric site (ITS) was identified in three chromosomal pairs (no. 2, 3 and 4). Additionally, two ITS sites were identified in the largest chromosome pair (pair 1), confirming the origin of this chromosome from three ancestral chromosomes. SATA and ITS sites allowed the orientation of the fusions in pairs 2, 3 and 4, which all appear to have been in different orientations (q-q, p-q and p-p, respectively). One of these fusions (O. karongae chromosome pair no. 2) involves a small chromosome (equivalent to linkage group 1), which in O. niloticus carries the main sex-determining gene. 4',6-Diamidino-2-phenyloindole staining of the synaptonemal complex in male O. karongae revealed the presumptive positions of the kinetochores, which correspond well to the centromeric positions observed in the mitotic karyotype.