Development of chromosomal markers based on next-generation sequencing: the B chromosome of the cichlid fish Astatotilapia latifasciata as a model

Comprehending the dynamism of B chromosomes in their journey towards becoming unselfish

Investigated for more than a century now, B chromosomes (Bs) research has come a long way from Bs being considered parasitic or neutral to becoming unselfish and bringing benefits to their hosts. B chromosomes exist as accessory chromosomes along with the standard A chromosomes (As) across eukaryotic taxa. Represented singly or in multiple copies, B chromosomes are largely heterochromatic but also contain euchromatic and organellar segments. Although B chromosomes are derived entities, they follow their species-specific evolutionary pattern. B chromosomes fail to pair with the standard chromosomes during meiosis and vary in their number, size, composition and structure across taxa and ensure their successful transmission through non-mendelian mechanisms like mitotic, pre-meiotic, meiotic or post-meiotic drives, unique non-disjunction, self-pairing or even imparting benefits to the host when they lack drive. B chromosomes have been associated with cellular processes like sex determination, pathogenicity, resistance to pathogens, phenotypic effects, and differential gene expression. With the advancements in B-omics research, novel insights have been gleaned on their functions, some of which have been associated with the regulation of gene expression of A chromosomes through increased expression of miRNAs or differential expression of transposable elements located on them. The next-generation sequencing and emerging technologies will further likely unravel the cellular, molecular and functional behaviour of these enigmatic entities. Amidst the extensive fluidity shown by B chromosomes in their structural and functional attributes, we perceive that the existence and survival of B chromosomes in the populations most likely seem to be a trade-off between the drive efficiency and adaptive significance versus their adverse effects on reproduction.

First characterization of PIWI-interacting RNA clusters in a cichlid fish with a B chromosome

BACKGROUND

B chromosomes are extra elements found in several eukaryote species. Usually, they do not express a phenotype in the host. However, advances in bioinformatics over the last decades have allowed us to describe several genes and molecular functions related to B chromosomes. These advances enable investigations of the relationship between the B chromosome and the host to understand how this element has been preserved in genomes. However, considering that transposable elements (TEs) are highly abundant in this supernumerary chromosome, there is a lack of knowledge concerning the dynamics of TE control in B-carrying cells. Thus, the present study characterized PIWI-interacting RNA (piRNA) clusters and pathways responsible for silencing the mobilization of TEs in gonads of the cichlid fish Astatotilapia latifasciata carrying the B chromosome.

RESULTS

Through small RNA-seq and genome assembly, we predicted and annotated piRNA clusters in the A. latifasciata genome for the first time. We observed that these clusters had biased expression related to sex and the presence of the B chromosome. Furthermore, three piRNA clusters, named curupira, were identified in the B chromosome. Two of them were expressed exclusively in gonads of samples with the B chromosome. The composition of these curupira sequences was derived from LTR, LINE, and DNA elements, representing old and recent transposition events in the A. latifasciata genome and the B chromosome. The presence of the B chromosome also affected the expression of piRNA pathway genes. The mitochondrial cardiolipin hydrolase-like (pld6) gene is present in the B chromosome, as previously reported, and an increase in its expression was detected in gonads with the B chromosome.

CONCLUSIONS

Due to the high abundance of TEs in the B chromosome, it was possible to investigate the origin of piRNA from these jumping genes. We hypothesize that the B chromosome has evolved its own genomic guardians to prevent uncontrolled TE mobilization. Furthermore, we also detected an expression bias in the presence of the B chromosome over A. latifasciata piRNA clusters and pathway genes.

Meiotic behavior, transmission and active genes of B chromosomes in the cichlid Astatotilapia latifasciata: new clues about nature, evolution and maintenance of accessory elements

Supernumerary B chromosomes (Bs) are dispensable genetic elements widespread in eukaryotes and are poorly understood mainly in relation to mechanisms of maintenance and transmission. The cichlid Astatotilapia latifasciata can harbor Bs in a range of 0 (named B -) and 1-2 (named B +). The B in A. latifasciata is rich in several classes of repetitive DNA sequences, contains protein coding genes, and affects hosts in diverse ways, including sex-biased effects. To advance in the knowledge about the mechanisms of maintenance and transmission of B chromosomes in A. latifasciata, here, we studied the meiotic behavior in males and transmission rates of A. latifasciata B chromosome. We also analyzed structurally and functionally the predicted B chromosome copies of the cell cycle genes separin-like, tubb1-like and kif11-like. We identified in the meiotic structure relative to the B chromosome the presence of proteins associated with Synaptonemal Complex organization (SMC3, SYCP1 and SYCP3) and found that the B performs self-pairing. These data suggest that isochromosome formation was a step during B chromosome evolution and this element is in a stage of diversification of the two arms keeping the self-pairing behavior to protect the A chromosome complement of negative effects of recombination. Moreover, we observed no occurrence of B-drive and confirmed the presence of cell cycle genes copies in the B chromosome and their transcription in encephalon, muscle and gonads, which can indicates beneficial effects to hosts and contribute to B maintenance.

Differential expression of miRNAs in the presence of B chromosome in the cichlid fish Astatotilapia latifasciata

BACKGROUND

B chromosomes (Bs) are extra elements observed in diverse eukaryotes, including animals, plants and fungi. Although Bs were first identified a century ago and have been studied in hundreds of species, their biology is still enigmatic. Recent advances in omics and big data technologies are revolutionizing the B biology field. These advances allow analyses of DNA, RNA, proteins and the construction of interactive networks for understanding the B composition and behavior in the cell. Several genes have been detected on the B chromosomes, although the interaction of B sequences and the normal genome remains poorly understood.

RESULTS

We identified 727 miRNA precursors in the A. latifasciata genome, 66% which were novel predicted sequences that had not been identified before. We were able to report the A. latifasciata-specific miRNAs and common miRNAs identified in other fish species. For the samples carrying the B chromosome (B⁺), we identified 104 differentially expressed (DE) miRNAs that are down or upregulated compared to samples without B chromosome (B^-) (p < 0.05). These miRNAs share common targets in the brain, muscle and gonads. These targets were used to construct a protein-protein-miRNA network showing the high interaction between the targets of differentially expressed miRNAs in the B⁺ chromosome samples. Among the DE-miRNA targets there are protein-coding genes reported for the B chromosome that are present in the protein-protein-miRNA network. Additionally, Gene Ontology (GO) terms related to nuclear matrix organization and response to stimulus are exclusive to DE miRNA targets of B⁺ samples.

CONCLUSIONS

This study is the first to report the connection of B chromosomes and miRNAs in a vertebrate species. We observed that the B chromosome impacts the miRNAs expression in several tissues and these miRNAs target several mRNAs involved with important biological processes.

Structure and Sequence of the Sex Determining Locus in Two Wild Populations of Nile Tilapia

In domesticated strains of the Nile tilapia, phenotypic sex has been linked to genetic variants on linkage groups 1, 20 and 23. This diversity of sex-loci might reflect a naturally polymorphic sex determination system in Nile tilapia, or it might be an artefact arising from the process of domestication. Here, we searched for sex-determiners in wild populations from Kpandu, Lake Volta (Ghana-West Africa), and from Lake Koka (Ethiopia-East Africa) that have not been subjected to any genetic manipulation. We analysed lab-reared families using double-digest Restriction Associated DNA sequencing (ddRAD) and analysed wild-caught males and females with pooled whole-genome sequencing (WGS). Strong sex-linked signals were found on LG23 in both populations, and sex-linked signals with LG3 were observed in Kpandu samples. WGS uncovered blocks of high sequence coverage, suggesting the presence of B chromosomes. We confirmed the existence of a tandem amh duplication in LG23 in both populations and determined its breakpoints between the oaz1 and dot1l genes. We found two common deletions of ~5 kb in males and confirmed the presence of both amhY and amh∆Y genes. Males from Lake Koka lack both the previously reported 234 bp deletion and the 5 bp frameshift-insertion that creates a premature stop codon in amh∆Y.

De novo genome assembly of the cichlid fish Astatotilapia latifasciata reveals a higher level of genomic polymorphism and genes related to B chromosomes

Supernumerary B chromosomes (Bs) are accessory elements to the regular chromosome set (As) and have been observed in a huge diversity of eukaryotic species. Although extensively investigated, the biological significance of Bs remains enigmatic. Here, we present de novo genome assemblies for the cichlid fish Astatotilapia latifasciata, a well-known model to study Bs. High coverage data with Illumina sequencing was obtained for males and females with 0B (B-), 1B, and 2B (B+) chromosomes to provide information regarding the diversity among these genomes. The draft assemblies comprised 771 Mb for the B- genome and 781 Mb for the B+ genome. Comparative analysis of the B+ and B- assemblies reveals syntenic discontinuity, duplicated blocks and several insertions, deletions, and inversions indicative of rearrangements in the B+ genome. Hundreds of transposable elements and 1546 protein coding sequences were annotated in the duplicated B+ regions. Our work contributes a list of thousands of genes harbored on the B chromosome, with functions in several biological processes, including the cell cycle.

Epigenetic DNA Modifications Are Correlated With B Chromosomes and Sex in the Cichlid Astatotilapia latifasciata

Supernumerary B chromosomes are dispensable elements found in several groups of eukaryotes, and their impacts in host organisms are not clear. The cichlid fish Astatotilapia latifasciata presents one or two large metacentric B chromosomes. These elements affect the transcription of several classes of RNAs. Here, we evaluated the epigenetic DNA modification status of B chromosomes using immunocytogenetics and assessed the impact of B chromosome presence on the global contents of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) and the molecular mechanisms underlying these variations. We found that the B chromosome of A. latifasciata has an active pattern of DNA epimarks, and its presence promotes the loss of 5mC in gonads of females with B chromosome (FB+) and promotes the loss of 5hmC in the muscle of males with the B element (MB+). Based on the transcriptional quantification of DNA modification genes (dnmt, tet, and tdg) and their candidate regulators (idh genes, microRNAs, and long non-coding RNAs) and on RNA-protein interaction prediction, we suggest the occurrence of passive demethylation in gonads of FB+ and 5hmC loss by Tet inhibition or by 5hmC oxidation in MB+ muscle. We suggest that these results can also explain the previously reported variations in the transcription levels of several classes of RNA depending on B chromosome presence. The DNA modifications detected here are also influenced by sex. Although the correlation between B chromosomes and sex has been previously reported, it remains unexplained. The B chromosome of A. latifasciata seems to be active and impacts cell physiology in a very complex way, including at the epigenetic level.

Landscape of Transposable Elements Focusing on the B Chromosome of the Cichlid Fish Astatotilapia latifasciata

B chromosomes (Bs) are supernumerary elements found in many taxonomic groups. Most B chromosomes are rich in heterochromatin and composed of abundant repetitive sequences, especially transposable elements (TEs). B origin is generally linked to the A-chromosome complement (A). The first report of a B chromosome in African cichlids was in Astatotilapia latifasciata, which can harbor 0, 1, or 2 Bs Classical cytogenetic studies found high a TE content on this B chromosome. In this study, we aimed to understand TE composition and expression in the A. latifasciata genome and its relation to the B chromosome. We used bioinformatics analysis to explore the genomic organization of TEs and their composition on the B chromosome. The bioinformatics findings were validated by fluorescent in situ hybridization (FISH) and real-time PCR (qPCR). A. latifasciata has a TE content similar to that of other cichlid fishes and several expanded elements on its B chromosome. With RNA sequencing data (RNA-seq), we showed that all major TE classes are transcribed in the brain, muscle, and male and female gonads. An evaluation of TE transcription levels between B- and B+ individuals showed that few elements are differentially expressed between these groups and that the expanded B elements are not highly transcribed. Putative silencing mechanisms may act on the B chromosome of A. latifasciata to prevent the adverse consequences of repeat transcription and mobilization in the genome.

The hnRNP Q-like gene is retroinserted into the B chromosomes of the cichlid fish Astatotilapia latifasciata

B chromosomes are dispensable elements observed in many eukaryotic species, including the African cichlid Astatotilapia latifasciata, which might have one or two B chromosomes. Although there have been many studies focused on the biology of these chromosomes, questions about the evolution, maintenance, and potential effects of these chromosomes remain. Here, we identified a variant form of the hnRNP Q-like gene inserted into the B chromosome of A. latifasciata that is characterized by a high copy number and intron-less structure. The absence of introns and presence of transposable elements with a reverse transcriptase domain flanking hnRNP Q-like sequences suggest that this gene was retroinserted into the B chromosome. RNA-Seq analysis did not show that the B variant retroinserted copies are transcriptionally active. However, RT-qPCR results showed variations in the canonical hnRNP Q-like copy expression levels among exons, tissues, sex, and B presence/absence. Although the patterns of transcription are not well understood, the exons of the B retrocopies were overexpressed, and a bias for female B+ expression was also observed. These results suggest that retroinsertion is an additional and important mechanism contributing to B chromosome formation. Furthermore, these findings indicate a bias towards female differential expression of B chromosome sequences, suggesting that B chromosomes and sex determination are somehow associated in cichlids.