Songbird germline-restricted chromosome as a potential arena of genetic conflicts

Genetic conflicts can arise between components of the genome with different inheritance strategies. The germline-restricted chromosome (GRC) of songbirds shows unusual mitotic and meiotic transmission compared with the rest of the genome. It is excluded from somatic cells and maintained only in the germline. It is usually present in one copy in the male germline and eliminated during spermatogenesis, while in the female germline, it usually occurs in two copies and behaves as a regular chromosome. Here, we review what is known about the GRC's evolutionary history, genetic content, and expression and discuss how it may be involved in different types of genetic conflicts. Finally, we interrogate the potential role of the GRC in songbird germline development, highlighting several unsolved mysteries.

Rapid gene content turnover on the germline-restricted chromosome in songbirds

The germline-restricted chromosome (GRC) of songbirds represents a taxonomically widespread example of programmed DNA elimination. Despite its apparent indispensability, we still know very little about the GRC's genetic composition, function, and evolutionary significance. Here we assemble the GRC in two closely related species, the common and thrush nightingale. In total we identify 192 genes across the two GRCs, with many of them present in multiple copies. Interestingly, the GRC appears to be under little selective pressure, with the genetic content differing dramatically between the two species and many GRC genes appearing to be pseudogenized fragments. Only one gene, cpeb1, has a complete coding region in all examined individuals of the two species and shows no copy number variation. The acquisition of this gene by the GRC corresponds with the earliest estimates of the GRC origin, making it a good candidate for the functional indispensability of the GRC in songbirds.

Micro germline-restricted chromosome in blue tits: evidence for meiotic functions

The germline-restricted chromosome (GRC) is likely present in all songbird species, but differs widely in size and gene content. This extra chromosome has been described as either a microchromosome with only limited basic gene content or a macrochromosome with enriched gene functions related to female gonad and embryo development. Here, we assembled, annotated and characterized the first micro-GRC in the blue tit (Cyanistes caeruleus) using high-fidelity long-read sequencing data. Although some genes on the blue tit GRC show signals of pseudogenization, others potentially have important functions, either currently or in the past. We highlight the GRC gene paralog BMP15, which is among the highest expressed GRC genes both in blue tits and in zebra finches (Taeniopygia guttata), and is known to play a role in oocyte and follicular maturation in other vertebrates. The GRC genes of the blue tit are further enriched for functions related to the synaptonemal complex. We found a similar functional enrichment when analyzing published data on GRC genes from two nightingale species (Luscinia spp.). We hypothesize that these genes play a role in maintaining standard maternal inheritance or in recombining maternal and paternal GRCs during potential episodes of biparental inheritance.

Mendelian nightmares: the germline-restricted chromosome of songbirds

Germline-restricted chromosomes (GRCs) are accessory chromosomes that occur only in germ cells. They are eliminated from somatic cells through programmed DNA elimination during embryo development. GRCs have been observed in several unrelated animal taxa and show peculiar modes of non-Mendelian inheritance and within-individual elimination. Recent cytogenetic and phylogenomic evidence suggests that a GRC is present across the species-rich songbirds, but absent in non-passerine birds, implying that over half of all 10,500 bird species have extensive germline/soma genome differences. Here, we review recent insights gained from genomic, transcriptomic, and cytogenetic approaches with regard to the genetic content, phylogenetic distribution, and inheritance of the songbird GRC. While many questions remain unsolved in terms of GRC inheritance, elimination, and function, we discuss plausible scenarios and future directions for understanding this widespread form of programmed DNA elimination.

Satellitome comparison of two oedipodine grasshoppers highlights the contingent nature of satellite DNA evolution

BACKGROUND

The full catalog of satellite DNA (satDNA) within a same genome constitutes the satellitome. The Library Hypothesis predicts that satDNA in relative species reflects that in their common ancestor, but the evolutionary mechanisms and pathways of satDNA evolution have never been analyzed for full satellitomes. We compare here the satellitomes of two Oedipodine grasshoppers (Locusta migratoria and Oedaleus decorus) which shared their most recent common ancestor about 22.8 Ma ago.

RESULTS

We found that about one third of their satDNA families (near 60 in every species) showed sequence homology and were grouped into 12 orthologous superfamilies. The turnover rate of consensus sequences was extremely variable among the 20 orthologous family pairs analyzed in both species. The satDNAs shared by both species showed poor association with sequence signatures and motives frequently argued as functional, except for short inverted repeats allowing short dyad symmetries and non-B DNA conformations. Orthologous satDNAs frequently showed different FISH patterns at both intra- and interspecific levels. We defined indices of homogenization and degeneration and quantified the level of incomplete library sorting between species.

CONCLUSIONS

Our analyses revealed that satDNA degenerates through point mutation and homogenizes through partial turnovers caused by massive tandem duplications (the so-called satDNA amplification). Remarkably, satDNA amplification increases homogenization, at intragenomic level, and diversification between species, thus constituting the basis for concerted evolution. We suggest a model of satDNA evolution by means of recursive cycles of amplification and degeneration, leading to mostly contingent evolutionary pathways where concerted evolution emerges promptly after lineages split.

Occasional paternal inheritance of the germline-restricted chromosome in songbirds

Songbirds have one special accessory chromosome, the so-called germline-restricted chromosome (GRC), which is only present in germline cells and absent from all somatic tissues. Earlier work on the zebra finch (Taeniopygia guttata castanotis) showed that the GRC is inherited only through the female line-like the mitochondria-and is eliminated from the sperm during spermatogenesis. Here, we show that the GRC has the potential to be paternally inherited. Confocal microscopy using GRC-specific fluorescent in situ hybridization probes indicated that a considerable fraction of sperm heads (1 to 19%) in zebra finch ejaculates still contained the GRC. In line with these cytogenetic data, sequencing of ejaculates revealed that individual males from two families differed strongly and consistently in the number of GRCs in their ejaculates. Examining a captive-bred male hybrid of the two zebra finch subspecies (T. g. guttata and T. g. castanotis) revealed that the mitochondria originated from a castanotis mother, whereas the GRC came from a guttata father. Moreover, analyzing GRC haplotypes across nine castanotis matrilines, estimated to have diverged for up to 250,000 y, showed surprisingly little variability among GRCs. This suggests that a single GRC haplotype has spread relatively recently across all examined matrilines. A few diagnostic GRC mutations that arose since this inferred spreading suggest that the GRC has continued to jump across matriline boundaries. Our findings raise the possibility that certain GRC haplotypes could selfishly spread through the population via occasional paternal transmission, thereby outcompeting other GRC haplotypes that were limited to strict maternal inheritance, even if this was partly detrimental to organismal fitness.

Out of patterns, the euchromatic B chromosome of the grasshopper Abracris flavolineata is not enriched in high-copy repeats

In addition to the normal set of standard (A) chromosomes, some eukaryote species harbor supernumerary (B) chromosomes. In most cases, B chromosomes show differential condensation with respect to A chromosomes and display dark C-bands of heterochromatin, and some of them are highly enriched in repetitive DNA. Here we perform a comprehensive NGS (next-generation sequencing) analysis of the repeatome in the grasshopper Abracris flavolineata aimed at uncovering the molecular composition and origin of its B chromosome. Our results have revealed that this B chromosome shows a DNA repeat content highly similar to the DNA repeat content observed for euchromatic (non-C-banded) regions of A chromosomes. Moreover, this B chromosome shows little enrichment for high-copy repeats, with only a few elements showing overabundance in B-carrying individuals compared to the 0B individuals. Consequently, the few satellite DNAs (satDNAs) mapping on the B chromosome were mostly restricted to its centromeric and telomeric regions, and they displayed much smaller bands than those observed on the A chromosomes. Our data support the intraspecific origin of the B chromosome from the longest autosome by misdivision, isochromosome formation, and additional restructuring, with accumulation of specific repeats in one or both B chromosome arms, yielding a submetacentric B. Finally, the absence of B-specific satDNAs, which are frequent in other species, along with its euchromatic nature, suggest that this B chromosome arose recently and might still be starting a heterochromatinization process. On this basis, it could be a good model to investigate the initial steps of B chromosome evolution.

Transposable element landscapes illuminate past evolutionary events in the endangered fern Vandenboschia speciosa

Vandenboschia speciosa is an endangered tetraploid fern species with a large genome (10.5 Gb). Its geographical distribution is characterized by disjoined tertiary flora refuges, with relict populations that survived past climate crises. Here, we analyzed the transposable elements (TEs) and found that they comprise approximately 76% of the V. speciosa genome, thus being the most abundant type of DNA sequence in this gigantic genome. The V. speciosa genome is composed of 51% and 5.6% of Class I and Class II elements, respectively. LTR retrotransposons were the most abundant TEs in this species (at least 42% of the genome), followed by non-LTR retrotransposons, which constituted at least 8.7% of the genome of this species. We introduce an additional analysis to identify the nature of non-annotated elements (19% of the genome). A BLAST search of the non-annotated contigs against the V. speciosa TE database allowed for the identification of almost half of them, which were most likely diverged sequence variants of the annotated TEs. In general, the TE composition in V. speciosa resembles the TE composition in seed plants. In addition, repeat landscapes revealed three episodes of amplification for all TEs, most likely due to demographic changes associated with past climate crises.

Long-term persistence of supernumerary B chromosomes in multiple species of Astyanax fish

BACKGROUND

Eukaryote genomes frequently harbor supernumerary B chromosomes in addition to the "standard" A chromosome set. B chromosomes are thought to arise as byproducts of genome rearrangements and have mostly been considered intraspecific oddities. However, their evolutionary transcendence beyond species level has remained untested.

RESULTS

Here we reveal that the large metacentric B chromosomes reported in several fish species of the genus Astyanax arose in a common ancestor at least 4 million years ago. We generated transcriptomes of A. scabripinnis and A. paranae 0B and 1B individuals and used these assemblies as a reference for mapping all gDNA and RNA libraries to quantify coverage differences between B-lacking and B-carrying genomes. We show that the B chromosomes of A. scabripinnis and A. paranae share 19 protein-coding genes, of which 14 and 11 were also present in the B chromosomes of A. bockmanni and A. fasciatus, respectively. Our search for B-specific single-nucleotide polymorphisms (SNPs) identified the presence of B-derived transcripts in B-carrying ovaries, 80% of which belonged to nobox, a gene involved in oogenesis regulation. Importantly, the B chromosome nobox paralog is expressed > 30× more than the A chromosome paralog. This indicates that the normal regulation of this gene is altered in B-carrying females, which could potentially facilitate B inheritance at higher rates than Mendelian law prediction.

CONCLUSIONS

Taken together, our results demonstrate the long-term survival of B chromosomes despite their lack of regular pairing and segregation during meiosis and that they can endure episodes of population divergence leading to species formation.

A Long-Term Conserved Satellite DNA That Remains Unexpanded in Several Genomes of Characiformes Fish Is Actively Transcribed

Eukaryotic genomes contain large amounts of repetitive DNA sequences, such as tandemly repeated satellite DNAs (satDNAs). These sequences are highly dynamic and tend to be genus- or species-specific due to their particular evolutionary pathways, although there are few unusual cases of conserved satDNAs over long periods of time. Here, we used multiple approaches to reveal that an satDNA named CharSat01-52 originated in the last common ancestor of Characoidei fish, a superfamily within the Characiformes order, ∼140-78 Ma, whereas its nucleotide composition has remained considerably conserved in several taxa. We show that 14 distantly related species within Characoidei share the presence of this satDNA, which is highly amplified and clustered in subtelomeric regions in a single species (Characidium gomesi), while remained organized as small clusters in all the other species. Defying predictions of the molecular drive of satellite evolution, CharSat01-52 shows similar values of intra- and interspecific divergence. Although we did not provide evidence for a specific functional role of CharSat01-52, its transcriptional activity was demonstrated in different species. In addition, we identified short tandem arrays of CharSat01-52 embedded within single-molecule real-time long reads of Astyanax paranae (536 bp-3.1 kb) and A. mexicanus (501 bp-3.9 kb). Such arrays consisted of head-to-tail repeats and could be found interspersed with other sequences, inverted sequences, or neighbored by other satellites. Our results provide a detailed characterization of an old and conserved satDNA, challenging general predictions of satDNA evolution.

High dynamism for neo-sex chromosomes: satellite DNAs reveal complex evolution in a grasshopper

A common characteristic of sex chromosomes is the accumulation of repetitive DNA, which accounts for their diversification and degeneration. In grasshoppers, the X0 sex-determining system in males is considered ancestral. However, in some species, derived variants like neo-XY in males evolved several times independently by Robertsonian translocation. This is the case of Ronderosia bergii, in which further large pericentromeric inversion in the neo-Y also took place, making this species particularly interesting for investigating sex chromosome evolution. Here, we characterized the satellite DNAs (satDNAs) and transposable elements (TEs) of the species to investigate the quantitative differences in repeat composition between male and female genomes putatively associated with sex chromosomes. We found a total of 53 satDNA families and 56 families of TEs. The satDNAs were 13.5% more abundant in males than in females, while TEs were just 1.02% more abundant in females. These results imply differential amplification of satDNAs on neo-Y chromosome and a minor role of TEs in sex chromosome differentiation. We showed highly differentiated neo-XY sex chromosomes owing to major amplification of satDNAs in neo-Y. Furthermore, chromosomal mapping of satDNAs suggests high turnover of neo-sex chromosomes in R. bergii at the intrapopulation level, caused by multiple paracentric inversions, amplifications, and transpositions. Finally, the species is an example of the action of repetitive DNAs in the generation of variability for sex chromosomes after the suppression of recombination, and helps understand sex chromosome evolution at the intrapopulation level.

Eight Million Years of Satellite DNA Evolution in Grasshoppers of the Genus Schistocerca Illuminate the Ins and Outs of the Library Hypothesis

Satellite DNA (satDNA) is an abundant class of tandemly repeated noncoding sequences, showing high rate of change in sequence, abundance, and physical location. However, the mechanisms promoting these changes are still controversial. The library model was put forward to explain the conservation of some satDNAs for long periods, predicting that related species share a common collection of satDNAs, which mostly experience quantitative changes. Here, we tested the library model by analyzing three satDNAs in ten species of Schistocerca grasshoppers. This group represents a valuable material because it diversified during the last 7.9 Myr across the American continent from the African desert locust (Schistocerca gregaria), and this thus illuminates the direction of evolutionary changes. By combining bioinformatic and cytogenetic, we tested whether these three satDNA families found in S. gregaria are also present in nine American species, and whether differential gains and/or losses have occurred in the lineages. We found that the three satDNAs are present in all species but display remarkable interspecies differences in their abundance and sequences while being highly consistent with genus phylogeny. The number of chromosomal loci where satDNA is present was also consistent with phylogeny for two satDNA families but not for the other. Our results suggest eminently chance events for satDNA evolution. Several evolutionary trends clearly imply either massive amplifications or contractions, thus closely fitting the library model prediction that changes are mostly quantitative. Finally, we found that satDNA amplifications or contractions may influence the evolution of monomer consensus sequences and by chance playing a major role in driftlike dynamics.

Programmed DNA elimination of germline development genes in songbirds

In some eukaryotes, germline and somatic genomes differ dramatically in their composition. Here we characterise a major germline–soma dissimilarity caused by a germline-restricted chromosome (GRC) in songbirds. We show that the zebra finch GRC contains >115 genes paralogous to single-copy genes on 18 autosomes and the Z chromosome, and is enriched in genes involved in female gonad development. Many genes are likely functional, evidenced by expression in testes and ovaries at the RNA and protein level. Using comparative genomics, we show that genes have been added to the GRC over millions of years of evolution, with embryonic development genes bicc1 and trim71 dating to the ancestor of songbirds and dozens of other genes added very recently. The somatic elimination of this evolutionarily dynamic chromosome in songbirds implies a unique mechanism to minimise genetic conflict between germline and soma, relevant to antagonistic pleiotropy, an evolutionary process underlying ageing and sexual traits.

Songbirds have extensive germline–soma genome differences due to developmental elimination of a germline-specific chromosome (GRC). Here, the authors show that the GRC contains dozens of expressed developmental genes, some of which have been on the GRC since the ancestor of all songbirds.

Differential Expression of Genes Related to Sexual Determination Can Modify the Reproductive Cycle of Astyanax scabripinnis (Characiformes: Characidae) in B Chromosome Carrier Individuals

The species complex Astyanax scabripinnis is one of the most studied with respect to origin, distribution, and frequency of B chromosomes, and is considered a model organism for evolutionary studies. Research using population inferences about the occurrence and frequency of the B chromosome shows seasonal variation between sexes, which is associated with the presence of this supernumerary element. We hypothesized that the B chromosome could influence the sex ratio of these animals. Based on this assumption, the present work aimed to investigate if differences exist among levels of gene expression with qRT-PCR of the amh (associated with testicular differentiation) and foxl2a (associated with ovarian differentiation) genes between B-carrier and non-B-carrier individuals. The results showed that for the amh gene, the difference in expression between animals with B chromosomes was not accentuated compared to that in animals without this chromosome. Expression of foxl2a in B-carrier females, however, was reduced by 73.56% compared to females that lacked the B chromosome. Males had no difference in expression of the amh and foxl2a genes between carriers and non-carriers of the B chromosome. Results indicate that the presence of B chromosomes is correlated with the differential expression of sex-associated genes. An analysis of these results integrated with data from other studies on the reproductive cycle in the same species reveals that this difference in expression may be expanding the reproductive cycle of the species.

Satellite DNA content of B chromosomes in the characid fish Characidium gomesi supports their origin from sex chromosomes

The origin of supernumerary (B) chromosomes is clearly conditioned by their ancestry from the standard (A) chromosomes. Sequence similarity between A and B chromosomes is thus crucial to determine B chromosome origin. For this purpose, we compare here the DNA sequences from A and B chromosomes in the characid fish Characidium gomesi using two main approaches. First, we found 59 satellite DNA (satDNA) families constituting the satellitome of this species and performed FISH analysis for 18 of them. This showed the presence of six satDNAs on the B chromosome: one shared with sex chromosomes and autosomes, two shared with sex chromosomes, one shared with autosomes and two being B-specific. This indicated that B chromosomes most likely arose from the sex chromosomes. Our second approach consisted of the analysis of five repetitive DNA families: 18S and 5S ribosomal DNA (rDNA), the H3 histone gene, U2 snDNA and the most abundant satDNA (CgoSat01-184) on DNA obtained from microdissected B chromosomes and from B-lacking genomes. PCR and sequence analysis of these repetitive sequences was successful for three of them (5S rDNA, H3 histone gene and CgoSat01-184), and sequence comparison revealed that DNA sequences obtained from the B chromosomes displayed higher identity with C. gomesi genomic DNA than with those obtained from other Characidium species. Taken together, our results support the intraspecific origin of B chromosomes in C. gomesi and point to sex chromosomes as B chromosome ancestors, which raises interesting prospects for future joint research on the genetic content of sex and B chromosomes in this species.

DNA metabarcoding of Neotropical ichthyoplankton: Enabling high accuracy with lower cost

Knowledge of ichthyoplankton dynamics is extremely important for conservation management as it can provide information about preferential spawning sites, reproductive period, migratory routes and recruitment success, which can be used to guide management and conservation efforts. However, identification of the eggs and larvae of Neotropical freshwater fish is a difficult task. DNA barcodes have emerged as an alternative and highly accurate approach for species identification, but DNA barcoding can be time-consuming and costly. To solve this problem, we aimed to develop a simple protocol based on DNA metabarcoding, to investigate whether it is possible to detect and quantify all species present in a pool of organisms. To do this, 230 larvae were cut in half, one half was sequenced by the Sanger technique and the other half was used to compose six arrays with a pool of larvae that were sequenced using a next-generation technique (NGS). The results of the Sanger sequencing allowed the identification of almost all larvae at species level, and the results from NGS showed high accuracy in species detection, ranging from 83% to 100%, with an average of 95% in all samples. No false positives were detected. The frequency of organisms in the two methods was positively correlated (Pearson), with low variation among species. In conclusion, this protocol represents a considerable advance in ichthyoplankton studies, allowing a rapid, cost-effective, quali-quantitative approach that improves the accuracy of identification.

Satellitome landscape analysis of Megaleporinus macrocephalus (Teleostei, Anostomidae) reveals intense accumulation of satellite sequences on the heteromorphic sex chromosome

The accumulation of repetitive DNA sequences on the sex-limited W or Y chromosomes is a well-known process that is likely triggered by the suppression of recombination between the sex chromosomes, which leads to major differences in their sizes and genetic content. Here, we report an analysis conducted on the satellitome of Megaleporinus macrocephalus that focuses specifically on the satDNAs that have been shown to have higher abundances in females and are putatively located on the W chromosome in this species. We characterized 164 satellite families in M. macrocephalus, which is, by far, the most satellite-rich species discovered to date. Subsequently, we mapped 30 satellites, 22 of which were located on the W chromosome, and 14 were shown to exist only on the W chromosome. Finally, we report two simple, quick and reliable methods that can be used for sex identification in M. macrocephalus individuals using fin clips or scales, which could be applicable to future studies conducted in the field of aquaculture.

Characterization of the satellitome in lower vascular plants: the case of the endangered fern Vandenboschia speciosa

BACKGROUND AND AIMS

Vandenboschia speciosa is a highly vulnerable fern species, with a large genome (10.5 Gb). Haploid gametophytes and diploid sporophytes are perennial, can reproduce vegetatively, and certain populations are composed only of independent gametophytes. These features make this fern a good model: (1) for high-throughput analysis of satellite DNA (satDNA) to investigate possible evolutionary trends in satDNA sequence features; (2) to determine the relative contribution of satDNA and other repetitive DNAs to its large genome; and (3) to analyse whether the reproduction mode or phase alternation between long-lasting haploid and diploid stages influences satDNA abundance or divergence.

METHODS

We analysed the repetitive fraction of the genome of this species in three different populations (one comprised only of independent gametophytes) using Illumina sequencing and bioinformatic analysis with RepeatExplorer and satMiner.

KEY RESULTS

The satellitome of V. speciosa is composed of 11 satDNA families, most of them showing a short repeat length and being A + T rich. Some satDNAs had complex repeats composed of sub-repeats, showing high similarity to shorter satDNAs. Three families had particular structural features and highly conserved motifs. SatDNA only amounts to approx. 0.4 % of its genome. Likewise, microsatellites do not represent more than 2 %, but transposable elements (TEs) represent approx. 50 % of the sporophytic genomes. We found high resemblance in satDNA abundance and divergence between both gametophyte and sporophyte samples from the same population and between populations.

CONCLUSIONS

(1) Longer (and older) satellites in V. speciosa have a higher A + T content and evolve from shorter ones and, in some cases, microsatellites were a source of new satDNAs; (2) the satellitome does not explain the huge genome size in this species while TEs are the major repetitive component of the V. speciosa genome and mostly contribute to its large genome; and (3) reproduction mode or phase alternation between gametophytes and sporophytes does not entail accumulation or divergence of satellites.

Complete mitochondrial genome of the Western Capercaillie Tetrao urogallus (Phasianidae, Tetraoninae)

The Western Capercaillie (Tetrao urogallus) is a galliform bird of boreal climax forests from Scandinavia to eastern Siberia, with a fragmented population in southwestern Europe. We extracted the DNA of T. urogallus aquitanicus and obtained the complete mitochondrial genome (mitogenome) sequence by combining Illumina and Sanger sequencing sequence data. The mitochondrial genome of T. urogallus is 16,683 bp long and is very similar to that of Lyrurus tetrix (16,677 bp). The T. urogallus mitogenome contains the normal 13 protein-coding genes (PCGs), 22 transfer RNAs, 2 ribosomal RNAs, and the control region. The number, order, and orientation of the mitochondrial genes are the same as in L. tetrix and in other species of the same and other bird families. The three domains of the control region contained conserved sequences (ETAS; CSBs), boxes (F, E, D, C, B, BS box), the putative origin of replication of the H-strand (OH) and bidirectional promoters of translation (LSP/HSP).

Full plastome sequence of the fern Vandenboschia speciosa (Hymenophyllales): structural singularities and evolutionary insights

We provide here the first full chloroplast genome sequence, i.e., the plastome, for a species belonging to the fern order Hymenophyllales. The phylogenetic position of this order within leptosporangiate ferns, together with the general scarcity of information about fern plastomes, places this research as a valuable study on the analysis of the diversity of plastomes throughout fern evolution. Gene content of V. speciosa plastome was similar to that in most ferns, although there were some characteristic gene losses and lineage-specific differences. In addition, an important number of genes required U to C RNA editing for proper protein translation and two genes showed start codons alternative to the canonical AUG (AUA). Concerning gene order, V. speciosa shared the specific 30-kb inversion of euphyllophytes plastomes and the 3.3-kb inversion of fern plastomes, keeping the ancestral gene order shared by eusporangiate and early leptosporangiate ferns. Conversely, V. speciosa has expanded IR regions comprising the rps7, rps12, ndhB and trnL genes in addition to rRNA and other tRNA genes, a condition shared with several eusporangiate ferns, lycophytes and hornworts, as well as most seed plants.

High-throughput analysis of satellite DNA in the grasshopper Pyrgomorpha conica reveals abundance of homologous and heterologous higher-order repeats

Satellite DNA (satDNA) constitutes an important fraction of repetitive DNA in eukaryotic genomes, but it is barely known in most species. The high-throughput analysis of satDNA in the grasshopper Pyrgomorpha conica revealed 87 satDNA variants grouped into 76 different families, representing 9.4% of the genome. Fluorescent in situ hybridization (FISH) analysis of the 38 most abundant satDNA families revealed four different patterns of chromosome distribution. Homology search between the 76 satDNA families showed the existence of 15 superfamilies, each including two or more families, with the most abundant superfamily representing more than 80% of all satDNA found in this species. This also revealed the presence of two types of higher-order repeats (HORs), one showing internal homologous subrepeats, as conventional HORs, and an additional type showing non-homologous internal subrepeats, the latter arising by the combination of a given satDNA family with a non-annotated sequence, or with telomeric DNA. Interestingly, the heterologous subrepeats included in these HORs showed higher divergence within the HOR than outside it, suggesting that heterologous HORs show poor homogenization, in high contrast with conventional (homologous) HORs. Finally, heterologous HORs can show high differences in divergence between their constituent subrepeats, suggesting the possibility of regional homogenization.

First complete female mitochondrial genome in four bivalve species genus Donax and their phylogenetic relationships within the Veneroida order

BACKGROUND

Four species of the genus Donax (D. semistriatus, D. trunculus, D. variegatus and D. vittatus) are common on Iberian Peninsula coasts. Nevertheless, despite their economic importance and overexploitation, scarce genetic resources are available. In this work, we newly determined the complete mitochondrial genomes of these four representatives of the family Donacidae, with the aim of contributing to unveil phylogenetic relationships within the Veneroida order, and of developing genetic markers being useful in wedge clam identification and authentication, and aquaculture stock management.

PRINCIPAL FINDINGS

The complete female mitochondrial genomes of the four species vary in size from 17,044 to 17,365 bp, and encode 13 protein-coding genes (including the atp8 gene), 2 rRNAs and 22 tRNAs, all located on the same strand. A long non-coding region was identified in each of the four Donax species between cob and cox2 genes, presumably corresponding to the Control Region. The Bayesian and Maximum Likelihood phylogenetic analysis of the Veneroida order indicate that all four species of Donax form a single clade as a sister group of other bivalves within the Tellinoidea superfamily. However, although Tellinoidea is actually monophyletic, none of its families are monophyletic.

CONCLUSIONS

Sequencing of complete mitochondrial genomes provides highly valuable information to establish the phylogenetic relationships within the Veneroida order. Furthermore, we provide here significant genetic resources for further research and conservation of this commercially important fishing resource.

Quantitative sequence characterization for repetitive DNA content in the supernumerary chromosome of the migratory locust

Repetitive DNA is a major component in most eukaryotic genomes but is ignored in most genome sequencing projects. Here, we report the quantitative composition in repetitive DNA for a supernumerary (B) chromosome, in the migratory locust (Locusta migratoria), by Illumina sequencing of genomic DNA from B-carrying and B-lacking individuals and DNA obtained from a microdissected B chromosome, as well as the physical mapping of some elements. B chromosome DNA of 94.9% was repetitive, in high contrast with the 64.1% of standard (A) chromosomes. B chromosomes are enriched in satellite DNA (satDNA) (65.2% of B-DNA), with a single satellite (LmiSat02-176) comprising 55% of the B. Six satDNAs were visualized by FISH on the B chromosome, and the only A chromosome carrying all these satellites was autosome 9, pointing to this chromosome, along with autosome 8 (which shares histone genes with the B) as putative ancestors of the B chromosome. We found several transposable elements (TEs) showing nucleotidic variation specific to B-carrying individuals, which was also present in B-carrying transcriptomes. Remarkably, an interstitial region of the B chromosome included a 17 kb chimera composed of 29 different TEs, suggesting reiterative TE insertion in this B chromosome region.

A Glimpse into the Satellite DNA Library in Characidae Fish (Teleostei, Characiformes)

Satellite DNA (satDNA) is an abundant fraction of repetitive DNA in eukaryotic genomes and plays an important role in genome organization and evolution. In general, satDNA sequences follow a concerted evolutionary pattern through the intragenomic homogenization of different repeat units. In addition, the satDNA library hypothesis predicts that related species share a series of satDNA variants descended from a common ancestor species, with differential amplification of different satDNA variants. The finding of a same satDNA family in species belonging to different genera within Characidae fish provided the opportunity to test both concerted evolution and library hypotheses. For this purpose, we analyzed here sequence variation and abundance of this satDNA family in ten species, by a combination of next generation sequencing (NGS), PCR and Sanger sequencing, and fluorescence in situ hybridization (FISH). We found extensive between-species variation for the number and size of pericentromeric FISH signals. At genomic level, the analysis of 1000s of DNA sequences obtained by Illumina sequencing and PCR amplification allowed defining 150 haplotypes which were linked in a common minimum spanning tree, where different patterns of concerted evolution were apparent. This also provided a glimpse into the satDNA library of this group of species. In consistency with the library hypothesis, different variants for this satDNA showed high differences in abundance between species, from highly abundant to simply relictual variants.

Protein-coding genes in B chromosomes of the grasshopper Eyprepocnemis plorans

For many years, parasitic B chromosomes have been considered genetically inert elements. Here we show the presence of ten protein-coding genes in the B chromosome of the grasshopper Eyprepocnemis plorans. Four of these genes (CIP2A, GTPB6, KIF20A, and MTG1) were complete in the B chromosome whereas the six remaining (CKAP2, CAP-G, HYI, MYCB2, SLIT and TOP2A) were truncated. Five of these genes (CIP2A, CKAP2, CAP-G, KIF20A, and MYCB2) were significantly up-regulated in B-carrying individuals, as expected if they were actively transcribed from the B chromosome. This conclusion is supported by three truncated genes (CKAP2, CAP-G and MYCB2) which showed up-regulation only in the regions being present in the B chromosome. Our results indicate that B chromosomes are not so silenced as was hitherto believed. Interestingly, the five active genes in the B chromosome code for functions related with cell division, which is the main arena where B chromosome destiny is played. This suggests that B chromosome evolutionary success can lie on its gene content.

The complete mitochondrial genome sequence of Astyanax paranae(Teleostei: characiformes)

Complete mitochondrial genome of the characiform fish Astyanax paranae was characterized in the present study. The whole mitogenome was 16,707 bp long and consisted of 13 protein-coding genes, 22 tRNAs, 2 rRNAs genes, a control region and origin of light-strand replication. The gene order is similar to those of the congeneric blind cavefish A. mexicanus. The nucleotide content of A. paranae mitogenome was 29.5% for A, 27.6% for T, 15.8% for G, 27.1% for C. Nucleotide identity between A. paranae and A. mexicanus across all the 37 genic regions ranged from 74.9% (ND2) to 90.3% (COX3).

High-throughput analysis of the satellitome illuminates satellite DNA evolution

Satellite DNA (satDNA) is a major component yet the great unknown of eukaryote genomes and clearly underrepresented in genome sequencing projects. Here we show the high-throughput analysis of satellite DNA content in the migratory locust by means of the bioinformatic analysis of Illumina reads with the RepeatExplorer and RepeatMasker programs. This unveiled 62 satDNA families and we propose the term "satellitome" for the whole collection of different satDNA families in a genome. The finding that satDNAs were present in many contigs of the migratory locust draft genome indicates that they show many genomic locations invisible by fluorescent in situ hybridization (FISH). The cytological pattern of five satellites showing common descent (belonging to the SF3 superfamily) suggests that non-clustered satDNAs can become into clustered through local amplification at any of the many genomic loci resulting from previous dissemination of short satDNA arrays. The fact that all kinds of satDNA (micro- mini- and satellites) can show the non-clustered and clustered states suggests that all these elements are mostly similar, except for repeat length. Finally, the presence of VNTRs in bacteria, showing similar properties to non-clustered satDNAs in eukaryotes, suggests that this kind of tandem repeats show common properties in all living beings.

Origin of B chromosomes in the genus Astyanax (Characiformes, Characidae) and the limits of chromosome painting

Eukaryote genomes are frequently burdened with the presence of supernumerary (B) chromosomes. Their origin is frequently investigated by chromosome painting, under the hypothesis that sharing the repetitive DNA sequences contained in the painting probes is a sign of common descent. However, the intragenomic mobility of many anonymous DNA sequences contained in these probes (e.g., transposable elements) adds high uncertainty to this conclusion. Here we test the validity of chromosome painting to investigate B chromosome origin by comparing its results for seven B chromosome types in two fish species genus Astyanax, with those obtained (1) by means of the physical mapping of 18S ribosomal DNA (rDNA), H1 histone genes, the As51 satellite DNA and the (AC)15 microsatellite, and (2) by comparing the nucleotide sequence of one of these families (ITS regions from ribosomal DNA) between genomic DNA from B-lacking individuals in both species and the microdissected DNA from two metacentric B chromosomes found in these same species. Intra- and inter-specific painting suggested that all B chromosomes that were assayed shared homologous DNA sequences among them, as well as with a variable number of A chromosomes in each species. This finding would be consistent with a common origin for all seven B chromosomes analyzed. By contrast, the physical mapping of repetitive DNA sequences failed to give support to this hypothesis, as no more than two B-types shared a given repetitive DNA. Finally, sequence analysis of the ITS regions suggested that at least some of the B chromosomes could have had a common origin.

Uncovering the Ancestry of B Chromosomes in Moenkhausia sanctaefilomenae (Teleostei, Characidae)

B chromosomes constitute a heterogeneous mixture of genomic parasites that are sometimes derived intraspecifically from the standard genome of the host species, but result from interspecific hybridization in other cases. The mode of origin determines the DNA content, with the B chromosomes showing high similarity with the A genome in the first case, but presenting higher similarity with a different species in the second. The characid fish Moenkhausia sanctaefilomenae harbours highly invasive B chromosomes, which are present in all populations analyzed to date in the Parana and Tietê rivers. To investigate the origin of these B chromosomes, we analyzed two natural populations: one carrying B chromosomes and the other lacking them, using a combination of molecular cytogenetic techniques, nucleotide sequence analysis and high-throughput sequencing (Illumina HiSeq2000). Our results showed that i) B chromosomes have not yet reached the Paranapanema River basin; ii) B chromosomes are mitotically unstable; iii) there are two types of B chromosomes, the most frequent of which is lightly C-banded (similar to euchromatin in A chromosomes) (B₁), while the other is darkly C-banded (heterochromatin-like) (B₂); iv) the two B types contain the same tandem repeat DNA sequences (18S ribosomal DNA, H3 histone genes, MS3 and MS7 satellite DNA), with a higher content of 18S rDNA in the heterochromatic variant; v) all of these repetitive DNAs are present together only in the paracentromeric region of autosome pair no. 6, suggesting that the B chromosomes are derived from this A chromosome; vi) the two B chromosome variants show MS3 sequences that are highly divergent from each other and from the 0B genome, although the B₂-derived sequences exhibit higher similarity with the 0B genome (this suggests an independent origin of the two B variants, with the less frequent, B₂ type presumably being younger); and vii) the dN/dS ratio for the H3.2 histone gene is almost 4–6 times higher for B chromosomes than for A chromosome sequences, suggesting that purifying selection is relaxed for the DNA sequences located on the B chromosomes, presumably because they are mostly inactive.

DNA barcoding of Iberian Peninsula and North Africa Tawny Owls Strix aluco suggests the Strait of Gibraltar as an important barrier for phylogeography

Eight subspecies have been proposed within the Tawny Owl (Strix aluco) species. However, recent molecular data have challenged this view, encouraging further work in this species complex. Here we reevaluated the taxonomic status between the North-Western African Tawny Owl, S. a. mauritanica, and its closest Iberian Tawny Owl population (from the S. a. sylvatica to S. a. aluco clade) separated by the Strait of Gibraltar. The Tawny Owl is a non-migratory and territorial species, and juvenile dispersal is restricted to a few kilometers around the natal site. This limited dispersal and the barrier imposed by the Strait of Gibraltar predicted a strong differentiation between the two populations. We tested this using DNA barcoding, Bayesian phylogenetic and species delimitation analysis. We found that an 81.1% of variation is due to the intergroups variation. In addition, the inter-intraspecific distances distribution revealed a barcoding gap among the two subspecies. Also, posterior probabilities and the P_AB value allowed to reject the hypothesis that observed degree of distinctiveness is due to random coalescence processes. These findings clearly support the Strait of Gibraltar as an isolating barrier for this species. The subspecific status is confirmed and species status is even suggested for S. a. mauritanica.

Non-random expression of ribosomal DNA units in a grasshopper showing high intragenomic variation for the ITS2 region

We analyse intragenomic variation of the ITS2 internal transcribed spacer of ribosomal DNA (rDNA) in the grasshopper Eyprepocnemis plorans, by means of tagged PCR 454 amplicon sequencing performed on both genomic DNA (gDNA) and RNA-derived complementary DNA (cDNA), using part of the ITS2 flanking coding regions (5.8S and 28S rDNA) as an internal control for sequencing errors. Six different ITS2 haplotypes (i.e. variants for at least one nucleotide in the complete ITS2 sequence) were found in a single population, one of them (Hap4) being specific to a supernumerary (B) chromosome. The analysis of both gDNA and cDNA from the same individuals provided an estimate of the expression efficiency of the different haplotypes. We found random expression (i.e. about similar recovery in gDNA and cDNA) for three haplotypes (Hap1, Hap2 and Hap5), but significant underexpression for three others (Hap3, Hap4 and Hap6). Hap4 was the most extremely underexpressed and, remarkably, it showed the lowest sequence conservation for the flanking 5.8-28S coding regions in the gDNA reads but the highest conservation (100%) in the cDNA ones, suggesting the preferential expression of mutation-free rDNA units carrying this ITS2 haplotype. These results indicate that the ITS2 region of rDNA is far from complete homogenization in this species, and that the different rDNA units are not expressed at random, with some of them being severely downregulated.

Intragenomic distribution of RTE retroelements suggests intrachromosomal movement

Much is known about the abundance of transposable elements (TEs) in eukaryotic genomes, but much is still unknown on their behaviour within cells. We employ here a combination of cytological, molecular and genomic approaches providing information on the intragenomic distribution and behaviour of non-long terminal repeat (LTR) retrotransposon-like elements (RTE). We microdissected every chromosome in a single first meiotic metaphase cell of the grasshopper Eyprepocnemis plorans and polymerase chain reaction (PCR) amplified a fragment of the RTE reverse transcriptase gene with specific primers. PCR products were cloned and 139 clones were sequenced. Analysis of molecular variance (AMOVA) showed significant intragenomic structure for these elements, with 4.6 % of molecular variance being found between chromosomes. A maximum likelihood tree built with the RTE sequences revealed the frequent presence of two or more elements showing very high similarity and being located on the same chromosome, thus suggesting intrachromosome movement. The 454 pyrosequencing of genomic DNA gave strong support to the microdissection results and provided evidence for the existence of 5' truncated elements. Our results thus indicate a tendency of RTE elements to reinsert into the same chromosome from where they were transcribed, which could be achieved if retrotranscription and insertion takes place immediately after transcription.

Genomics of ecological adaptation in cactophilic Drosophila

Cactophilic Drosophila species provide a valuable model to study gene–environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii–D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii–D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.

A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs

The desert locust (Schistocerca gregaria) has been used as material for numerous cytogenetic studies. Its genome size is estimated to be 8.55 Gb of DNA comprised in 11 autosomes and the X chromosome. Its X0/XX sex chromosome determinism therefore results in females having 24 chromosomes whereas males have 23. Surprisingly, little is known about the DNA content of this locust's huge chromosomes. Here, we use the Feulgen Image Analysis Densitometry and C-banding techniques to respectively estimate the DNA quantity and heterochromatin content of each chromosome. We also identify three satellite DNAs using both restriction endonucleases and next-generation sequencing. We then use fluorescent in situ hybridization to determine the chromosomal location of these satellite DNAs as well as that of six tandem repeat DNA gene families. The combination of the results obtained in this work allows distinguishing between the different chromosomes not only by size, but also by the kind of repetitive DNAs that they contain. The recent publication of the draft genome of the migratory locust (Locusta migratoria), the largest animal genome hitherto sequenced, invites for sequencing even larger genomes. S. gregaria is a pest that causes high economic losses. It is thus among the primary candidates for genome sequencing. But this species genome is about 50 % larger than that of L. migratoria, and although next-generation sequencing currently allows sequencing large genomes, sequencing it would mean a greater challenge. The chromosome sizes and markers provided here should not only help planning the sequencing project and guide the assembly but would also facilitate assigning assembled linkage groups to actual chromosomes.

Next generation sequencing and FISH reveal uneven and nonrandom microsatellite distribution in two grasshopper genomes

Simple sequence repeats (SSRs), also known as microsatellites, are one of the prominent DNA sequences shaping the repeated fraction of eukaryotic genomes. In spite of their profuse use as molecular markers for a variety of genetic and evolutionary studies, their genomic location, distribution, and function are not yet well understood. Here we report the first thorough joint analysis of microsatellite motifs at both genomic and chromosomal levels in animal species, by a combination of 454 sequencing and fluorescent in situ hybridization (FISH) techniques performed on two grasshopper species. The in silico analysis of the 454 reads suggested that microsatellite expansion is not driving size increase of these genomes, as SSR abundance was higher in the species showing the smallest genome. However, the two species showed the same uneven and nonrandom location of SSRs, with clear predominance of dinucleotide motifs and association with several types of repetitive elements, mostly histone gene spacers, ribosomal DNA intergenic spacers (IGS), and transposable elements (TEs). The FISH analysis showed a dispersed chromosome distribution of microsatellite motifs in euchromatic regions, in coincidence with chromosome location patterns previously observed for many mobile elements in these species. However, some SSR motifs were clustered, especially those located in the histone gene cluster.

Delimiting the origin of a B chromosome by FISH mapping, chromosome painting and DNA sequence analysis in Astyanax paranae (Teleostei, Characiformes)

Supernumerary (B) chromosomes have been shown to contain a wide variety of repetitive sequences. For this reason, fluorescent in situ hybridisation (FISH) is a useful tool for ascertaining the origin of these genomic elements, especially when combined with painting from microdissected B chromosomes. In order to investigate the origin of B chromosomes in the fish species Astyanax paranae, these two approaches were used along with PCR amplification of specific DNA sequences obtained from the B chromosomes and its comparison with those residing in the A chromosomes. Remarkably, chromosome painting with the one-arm metacentric B chromosome probe showed hybridization signals on entire B chromosome, while FISH mapping revealed the presence of H1 histone and 18S rDNA genes symmetrically placed in both arms of the B chromosome. These results support the hypothesis that the B chromosome of A. paranae is an isochromosome. Additionally, the chromosome pairs Nos. 2 or 23 are considered the possible B chromosome ancestors since both contain syntenic H1 and 18S rRNA sequences. The analysis of DNA sequence fragments of the histone and rRNA genes obtained from the microdissected B chromosomes showed high similarity with those obtained from 0B individuals, which supports the intraspecific origin of B chromosomes in A. paranae. Finally, the population hereby analysed showed a female-biased B chromosome presence suggesting that B chromosomes in this species could influence sex determinism.

A mobile insulator system to detect and disrupt cis-regulatory landscapes in vertebrates

In multicellular organisms, cis-regulation controls gene expression in space and time. Despite the essential implication of cis-regulation in the development and evolution of organisms and in human diseases, our knowledge about regulatory sequences largely derives from analyzing their activity individually and outside their genomic context. Indeed, the contribution of these sequences to the expression of their target genes in their genomic context is still largely unknown. Here we present a novel genetic screen designed to visualize and interrupt gene regulatory landscapes in vertebrates. In this screen, based on the random insertion of an engineered Tol2 transposon carrying a strong insulator separating two fluorescent reporter genes, we isolated hundreds of zebrafish lines containing insertions that disrupt the cis-regulation of tissue-specific expressed genes. We therefore provide a new easy-to-handle tool that will help to disrupt and chart the regulatory activity spread through the vast noncoding regions of the vertebrate genome.

DNA amount of X and B chromosomes in the grasshoppers Eyprepocnemis plorans and Locusta migratoria

We analyzed the DNA amount in X and B chromosomes of 2 XX/X0 grasshopper species (Eyprepocnemis plorans and Locusta migratoria), by means of Feulgen image analysis densitometry (FIAD), using previous estimates in L. migratoria as standard (5.89 pg). We first analyzed spermatids of 0B males and found a bimodal distribution of integrated optical densities (IODs), suggesting that one peak corresponded to +X and the other to -X spermatids. The difference between the 2 peaks corresponded to the X chromosome DNA amount, which was 1.28 pg in E. plorans and 0.80 pg in L. migratoria. In addition, the +X peak in E. plorans gave an estimate of the C-value in this species (10.39 pg). We next analyzed diplotene cells from 1B males in E. plorans and +B males in L. migratoria (a species where Bs are mitotically unstable and no integer B number can be defined for an individual) and measured B chromosome IOD relative to X chromosome IOD, within the same cell, taking advantage of the similar degree of condensation for both positively heteropycnotic chromosomes at this meiotic stage. From this proportion, we estimated the DNA amount for 3 different B chromosome variants found in individuals from 3 E. plorans Spanish populations (0.54 pg for B1 from Saladares, 0.51 pg for B2 from Salobreña and 0.64 for B24 from Torrox). Likewise, we estimated the DNA amount of the B chromosome in L. migratoria to be 0.15 pg. To automate measurements, we wrote a GPL3 licensed Python program (pyFIA). We discuss the utility of the present approach for estimating X and B chromosome DNA amount in a variety of situations, and the meaning of the DNA amount estimates for X and B chromosomes in these 2 species.