Chicken rRNA Gene Cluster Structure

Description of the complete rDNA repeat unit structure of Coturnixjaponica Temminck et Schlegel, 1849 (Aves)

Ribosomal RNA (18S, 5.8S, 28S) gene clusters in genomes form regions that consist of multiple tandem repeats. They are located on a single or several pairs of chromosomes and play an important role in the formation of the nucleolus responsible for the assembly of ribosome subunits. The rRNA gene cluster sequences are widely used for taxonomic studies, however at present, complete information on the avian rDNA repeat unit structure including intergenic spacer sequence is available only for the chicken (Gallusgallusdomesticus Linnaeus, 1758). The GC enrichment and high-order repeats peculiarities within the intergenic spacer described for the chicken rDNA cluster may be responsible for these failures. The karyotype of the Japanese quail (Coturnixjaponica Temminck et Schlegel, 1849) deserves close attention because, unlike most birds, it has three pairs of nucleolar organizer bearing chromosomes, two of which are microchromosomes enriched in repeating elements and heterochromatin that carry translocated terminal nucleolar organizers. Here we assembled and annotated the complete Japanese quail ribosomal gene cluster sequence of 21166 base pairs (GenBank under the registration tag BankIt2509210 CoturnixOK523374). This is the second deciphered avian rDNA cluster after the chicken. Despite the revealed high similarity with the chicken corresponding sequence, it has a number of specific features, which include a slightly lower degree of GC content and the presence of bendable elements in the content of both the transcribed spacer I and the non-transcribed intergenic spacer.

Alternative probe hybridization buffers for target RNA depletion and viral sequence recovery in NGS for poultry samples

Non-targeted next generation sequencing (NGS) is widely applied to identify the diversity of pathogens in field samples. However, abundance of host RNA (especially rRNA) and other environmental nucleic acids can reduce the abundance of pathogen specific reads of interest, reduce depth of coverage and increase surveillance costs. We presently deplete chicken- and selected bacterial-specific rRNAs in poultry field RNA samples with complementary DNA probes in a commercially available probe hybridization buffer followed by digestion of the RNA:DNA hybrids with RNase H. Because the current buffer is an expensive special order reagent of proprietary composition, we tested in-house and other commercially available buffers and identified a viable alternative that yields equivalent host rRNA depletion and viral-specific reads in poultry samples as the current special order reagent but at a reduced cost.

Single Copies of the 5S rRNA Inserted into 45S rDNA Intergenic Spacers in the Genomes of Nototheniidae (Perciformes, Actinopterygii)

In the vast majority of Animalia genomes, the 5S rRNA gene repeats are located on chromosomes outside of the 45S rDNA arrays of the nucleolar organiser (NOR). We analysed the genomic databases available and found that a 5S rDNA sequence is inserted into the intergenic spacer (IGS) between the 45S rDNA repeats in ten species of the family Nototheniidae (Perciformes, Actinopterigii). We call this sequence the NOR-5S rRNA gene. Along with Testudines and Crocodilia, this is the second case of a close association between four rRNA genes within one repetitive unit in deuterostomes. In both cases, NOR-5S is oriented opposite the 45S rDNA. None of the three nucleotide substitutions compared to the canonical 5S rRNA gene influenced the 5S rRNA secondary structure. In transcriptomes of the Patagonian toothfish, we only found NOR-5S rRNA reads in ovaries and early embryos, but not in testis or somatic tissues of adults. Thus, we consider the NOR-5S gene to be a maternal-type 5S rRNA template. The colocalization of the 5S and 45S ribosomal genes appears to be essential for the equimolar production of all four rRNAs in the species that show rDNA amplification during oogenesis. Most likely, the integration of 5S and NOR rRNA genes occurred prior to Nototheniidae lineage diversification.

Molecular identification and genetic diversity of equine ocular setariasis in Thailand based on the COI, 12S rDNA, and ITS1 regions

Equine ocular setariasis is mainly caused by Setaria digitata, and the identification of this filarial nematode is based on morphology. However, morphological characterization alone is insufficient for the detection and differentiation of S. digitata from its congeners. In Thailand, the molecular detection of S. digitata is lacking and its genetic diversity is still unknown. This study aimed to phylogenetically characterize equine S. digitata from Thailand based on sequences derived from the mitochondrial cytochrome c oxidase subunit 1 (COI), the mitochondrial small subunit ribosomal DNA (12S rDNA), the nuclear internal transcribed spacer 1 (ITS1) and Wolbachia surface protein (wsp). Five samples of S. digitata were characterized, submitted to the NCBI database, and used for phylogenetic analysis as well as the assessment of similarity, entropy, and haplotype diversity. Phylogenetic analyses revealed that the S. digitata Thai strain was similar to S. digitata from China and Sri Lanka, with 99 to 100% similarity. The entropy and haplotype diversity indicated that the S. digitata Thai isolate was conserved and closely related to S. digitata worldwide. This is the first report on the molecular detection of equine ocular setariasis caused by S. digitata in Thailand.

Chromosome-Length Assembly of the Baikal Seal (Pusa sibirica) Genome Reveals a Historically Large Population Prior to Isolation in Lake Baikal

Pusa sibirica, the Baikal seal, is the only extant, exclusively freshwater, pinniped species. The pending issue is, how and when they reached their current habitat—the rift lake Baikal, more than three thousand kilometers away from the Arctic Ocean. To explore the demographic history and genetic diversity of this species, we generated a de novo chromosome-length assembly, and compared it with three closely related marine pinniped species. Multiple whole genome alignment of the four species compared with their karyotypes showed high conservation of chromosomal features, except for three large inversions on chromosome VI. We found the mean heterozygosity of the studied Baikal seal individuals was relatively low (0.61 SNPs/kbp), but comparable to other analyzed pinniped samples. Demographic reconstruction of seals revealed differing trajectories, yet remarkable variations in Ne occurred during approximately the same time periods. The Baikal seal showed a significantly more severe decline relative to other species. This could be due to the difference in environmental conditions encountered by the earlier populations of Baikal seals, as ice sheets changed during glacial–interglacial cycles. We connect this period to the time of migration to Lake Baikal, which occurred ~3–0.3 Mya, after which the population stabilized, indicating balanced habitat conditions.

The Gallus gallus RJF reference genome reveals an MHCY haplotype organized in gene blocks that contain 107 loci including 45 specialized, polymorphic MHC class I loci, 41 C-type lectin-like loci, and other loci amid hundreds of transposable elements

MHCY is a second major histocompatibility complex-like gene region in chickens originally identified by the presence of major histocompatibility complex class I-like and class II-like gene sequences. Up to now, the MHCY gene region has been poorly represented in genomic sequence data. A high density of repetitive sequence and multiple members of several gene families prevented the accurate assembly of short-read sequence data for MHCY. Identified here by single-molecule real-time sequencing sequencing of BAC clones for the Gallus gallus Red Jungle Fowl reference genome are 107 MHCY region genes (45 major histocompatibility complex class I-like, 41 c-type-lectin-like, 8 major histocompatibility complex class IIβ, 8 LENG9-like, 4 zinc finger protein loci, and a single only zinc finger-like locus) located amid hundreds of retroelements within 4 contigs representing the region. Sequences obtained for nearby ribosomal RNA genes have allowed MHCY to be precisely mapped with respect to the nucleolar organizer region. Gene sequences provide insights into the unusual structure of the MHCY class I molecules. The MHCY class I loci are polymorphic and group into 22 types based on predicted amino acid sequences. Some MHCY class I loci are full-length major histocompatibility complex class I genes. Others with altered gene structure are considered gene candidates. The amino acid side chains at many of the polymorphic positions in MHCY class I are directed away rather than into the antigen-binding groove as is typical of peptide-binding major histocompatibility complex class I molecules. Identical and nearly identical blocks of genomic sequence contribute to the observed multiplicity of identical MHCY genes and the large size (>639 kb) of the Red Jungle Fowl MHCY haplotype. Multiple points of hybridization observed in fluorescence in situ hybridization suggest that the Red Jungle Fowl MHCY haplotype is made up of linked, but physically separated genomic segments. The unusual gene content, the evidence of highly similar duplicated segments, and additional evidence of variation in haplotype size distinguish polymorphic MHCY from classical polymorphic major histocompatibility complex regions.

45S rDNA Repeats of Turtles and Crocodiles Harbor a Functional 5S rRNA Gene Specifically Expressed in Oocytes

In most eukaryotic genomes, tandemly repeated copies of 5S rRNA genes are clustered outside the nucleolus organizer region (NOR), which normally encodes three other major rRNAs: 18S, 5.8S, and 28S. Our analysis of turtle rDNA sequences has revealed a 5S rDNA insertion into the NOR intergenic spacer in antisense orientation. The insertion (hereafter called NOR-5S rRNA gene) has a length of 119 bp and coexists with the canonical 5S rDNA clusters outside the NOR. Despite the ∼20% nucleotide difference between the two 5S gene sequences, their internal control regions for RNA polymerase III are similar. Using the turtle Trachemys scripta as a model species, we showed the NOR-5S rDNA specific expression in oocytes. This expression is concurrent with the NOR rDNA amplification during oocyte growth. We show that in vitellogenic oocytes, the NOR-5S rRNA prevails over the canonical 5S rRNA in the ribosomes, suggesting a role of modified ribosomes in oocyte-specific translation. The orders Testudines and Crocodilia seem to be the only taxa of vertebrates with such a peculiar rDNA organization. We speculate that the amplification of the 5S rRNA genes as a part of the NOR DNA during oogenesis provides a dosage balance between transcription of all the four ribosomal RNAs while producing a maternal pool of extra ribosomes. We further hypothesize that the NOR-5S rDNA insertion appeared in the Archelosauria clade during the Permian period and was lost later in the ancestors of Aves.

Analyses of the Updated "Animal rDNA Loci Database" with an Emphasis on Its New Features

We report on a major update to the animal rDNA loci database, which now contains cytogenetic information for 45S and 5S rDNA loci in more than 2600 and 1000 species, respectively.The data analyses show the following: (i) A high variability in 5S and 45S loci numbers, with both showing 50-fold or higher variability. However, karyotypes with an extremely high number of loci were rare, and medians generally converged to two 5S sites and two 45S rDNA sites per diploid genome. No relationship was observed between the number of 5S and 45S loci. (ii) The position of 45S rDNA on sex chromosomes was relatively frequent in some groups, particularly in arthropods (14% of karyotypes). Furthermore, 45S rDNA was almost exclusively located in microchromosomes when these were present (in birds and reptiles). (iii) The proportion of active NORs (positively stained with silver staining methods) progressively decreased with an increasing number of 45S rDNA loci, and karyotypes with more than 12 loci showed, on average, less than 40% of active loci. In conclusion, the updated version of the database provides some new insights into the organization of rRNA genes in chromosomes. We expect that its updated content will be useful for taxonomists, comparative cytogeneticists, and evolutionary biologists. .

Ageing European lobsters (Homarus gammarus) using DNA methylation of evolutionarily conserved ribosomal DNA

Crustaceans are notoriously difficult to age because of their indeterminate growth and the moulting of their exoskeleton throughout life. The poor knowledge of population age structure in crustaceans therefore hampers accurate assessment of population dynamics and consequently sustainable fisheries management. Quantification of DNA methylation of the evolutionarily conserved ribosomal DNA (rDNA) may allow for age prediction across diverse species. However, the rDNA epigenetic clock remains to be tested in crustaceans, despite its potential to inform both ecological and evolutionary understanding, as well as conservation and management practices. Here, patterns of rDNA methylation with age were measured across 5154 bp of rDNA corresponding to 355 quality-filtered loci in the economically important European lobster (Homarus gammarus). Across 0- to 51-month-old lobsters (n = 155), there was a significant linear relationship between age and percentage rDNA methylation in claw tissue at 60% of quality-filtered loci (n = 214). An Elastic Net regression model using 46 loci allowed for the accurate and precise age estimation of individuals (R 2 = 0.98; standard deviation = 1.6 months). Applying this ageing model to antennal DNA from wild lobsters of unknown age (n = 38) resulted in predicted ages that are concordant with estimates of minimum size at age in the wild (mean estimated age = 40.1 months; range 32.8-55.7 months). Overall, the rDNA epigenetic clock shows potential as a novel, nonlethal ageing technique for European lobsters. However, further validation is required across a wider range of known-age individuals and tissue types before the model can be used in fisheries management.

Repeat Sequence Mapping Shows Different W Chromosome Evolutionary Pathways in Two Caprimulgiformes Families

Although birds belonging to order Caprimulgiformes show extensive karyotype variation, data concerning their genomic organization is still scarce, as most studies have presented only results obtained from conventional staining analyses. Nevertheless, some interesting findings have been observed, such as the W chromosome of the Common Potoo, Nyctibius griseus (2n = 86), which has the same morphology and size of the Z chromosome, a rare feature in Neognathae birds. Hence, we aimed to investigate the process by which the W chromosome of this species was enlarged. For that, we analyzed comparatively the chromosome organization of the Common Potoo and the Scissor-tailed Nightjar, Hydropsalis torquata (2n = 74), which presents the regular differentiated sex chromosomes, by applying C-banding, G-banding and mapping of repetitive DNAs (microsatellite repeats and 18S rDNA). Our results showed an accumulation of constitutive heterochromatin in the W chromosome of both species. However, 9 out of 11 microsatellite sequences hybridized in the large W chromosome in the Common Potoo, while none of them hybridized in the W chromosome of the Scissor-tailed Nightjar. Therefore, we can conclude that the accumulation of microsatellite sequences, and consequent increase in constitutive heterochromatin, was responsible for the enlargement of the W chromosome in the Common Potoo. Based on these results, we conclude that even though these two species belong to the same order, their W chromosomes have gone through different evolutionary histories, with an extra step of accumulation of repetitive sequences in the Common Potoo.

The distribution of 45S rDNA sites in bird chromosomes suggests multiple evolutionary histories

The distribution of 45S rDNA cluster in avian karyotypes varies in different aspects, such as position, number of bearer chromosomes, and bearers being macro- or microchromosomes. The present study investigated the patterns of variation in the 45S rDNA-bearer chromosomes of birds in order to understand the evolutionary dynamics of the cluster configuration and its contribution to the evolution of bird karyotypes. A total of 73 bird species were analyzed, including both published data and species for which rDNA-FISH was conducted for the first time. In most birds, the 45S rDNA clusters were located in a single pair of microchromosomes. Hence, the location of 45S rDNA in macrochromosomes, observed only in Neognathae species, seems to be a derived state, probably the result of chromosomal fusion between microchromosomes and distinct macrochromosomes. Additionally, the 45S rDNA was observed in multiple microchromosomes in different branches of the bird phylogeny, suggesting recurrence of dispersion processeses, such as duplications and translocations. Overall, this study indicated that the redistribution of the 45S rDNA sites in bird chromosomes followed different evolutionary trajectories with respect to each lineage of the class Aves.

Comparative Cytogenetics of Four Sea Turtle Species (Cheloniidae): G-Banding Pattern and in situ Localization of Repetitive DNA Units

Sea turtles are considered flagship species for marine biodiversity conservation and are considered to be at varying risk of extinction globally. Cases of hybridism have been reported in sea turtles, but chromosomal analyses are limited to classical karyotype descriptions and a few molecular cytogenetic studies. In order to compare karyotypes and understand evolutive mechanisms related to chromosome dif­ferentiation in this group, <i>Chelonia mydas</i>, <i>Caretta caretta</i>, <i>Eretmochelys imbricata</i>, and <i>Lepidochelys olivacea</i> were cytogenetically characterized in the present study. When the obtained cytogenetic data were compared with the putative ancestral Cryptodira karyotype, the studied species showed the same diploid number (2n) of 56 chromosomes, with some variations in chromosomal morphology (karyotypic formula) and minor changes in longitudinal band locations. In situ localization using a 18S ribosomal DNA probe indicated a homeologous microchromosome pair bearing a 45S ribosomal DNA locus and size heteromorphism in all 4 species. Interstitial telomeric sites were identified in a microchromosome pair in <i>C. mydas</i> and <i>C. caretta</i>. The data showed that interspecific variations occurred in chromosomal sets among the Cheloniidae species, in addition to other Cryptodira karyotypes. These variations generated lineage-specific karyotypic diversification in sea turtles, which will have considerable implications for hybrid recognition and for the study, the biology, ecology, and evolutionary history of regional and global populations. Furthermore, we demonstrated that some chromosome rearrangements occurred in sea turtle species, which is in conflict with the hypothesis of conserved karyotypes in this group.

Structure of the intergenic spacers in chicken ribosomal DNA

BACKGROUND

Ribosomal DNA (rDNA) repeats are situated in the nucleolus organizer regions (NOR) of chromosomes and transcribed into rRNA for ribosome biogenesis. Thus, they are an essential component of eukaryotic genomes. rDNA repeat units consist of rRNA gene clusters that are transcribed into single pre-rRNA molecules, each separated by intergenic spacers (IGS) that contain regulatory elements for rRNA gene cluster transcription. Because of their high repeat content, rDNA sequences are usually absent from genome assemblies. In this work, we used the long-read sequencing technology to describe the chicken IGS and fill the knowledge gap on rDNA sequences of one of the key domesticated animals.

METHODS

We used the long-read PacBio RSII technique to sequence the BAC clone WAG137G04 (Wageningen BAC library) known to contain chicken NOR elements and the HGAP workflow software suit to assemble the PacBio RSII reads. Whole-genome sequence contigs homologous to the chicken rDNA repetitive unit were identified based on the Gallus_gallus-5.0 assembly with BLAST. We used the Geneious 9.0.5 and Mega software, maximum likelihood method and Chickspress project for sequence evolution analysis, phylogenetic tree construction and analysis of the raw transcriptome data.

RESULTS

Three complete IGS sequences in the White Leghorn chicken genome and one IGS sequence in the red junglefowl contig AADN04001305.1 (Gallus_gallus-5.0) were detected. They had various lengths and contained three groups of tandem repeats (some of them being very GC rich) that form highly organized arrays. Initiation and termination sites of rDNA transcription were located within small and large unique regions (SUR and LUR), respectively. No functionally significant sites were detected within the tandem repeat sequences.

CONCLUSIONS

Due to the highly organized GC-rich repeats, the structure of the chicken IGS differs from that of IGS in human, apes, Xenopus or fish rDNA. However, the chicken IGS shares some molecular organization features with that of the turtles, which are other representatives of the Sauropsida clade that includes birds and reptiles. Our current results on the structure of chicken IGS together with the previously reported ribosomal gene cluster sequence provide sufficient data to consider that the complete chicken rDNA sequence is assembled with confidence in terms of molecular DNA organization.

Sequence properties of certain GC rich avian genes, their origins and absence from genome assemblies: case studies

Background

More and more eukaryotic genomes are sequenced and assembled, most of them presented as a complete model in which missing chromosomal regions are filled by Ns and where a few chromosomes may be lacking. Avian genomes often contain sequences with high GC content, which has been hypothesized to be at the origin of many missing sequences in these genomes. We investigated features of these missing sequences to discover why some may not have been integrated into genomic libraries and/or sequenced.

Results

The sequences of five red jungle fowl cDNA models with high GC content were used as queries to search publicly available datasets of Illumina and Pacbio sequencing reads. These were used to reconstruct the leptin, TNFα, MRPL52, PCP2 and PET100 genes, all of which are absent from the red jungle fowl genome model. These gene sequences displayed elevated GC contents, had intron sizes that were sometimes larger than non-avian orthologues, and had non-coding regions that contained numerous tandem and inverted repeat sequences with motifs able to assemble into stable G-quadruplexes and intrastrand dyadic structures. Our results suggest that Illumina technology was unable to sequence the non-coding regions of these genes. On the other hand, PacBio technology was able to sequence these regions, but with dramatically lower efficiency than would typically be expected.

Conclusions

High GC content was not the principal reason why numerous GC-rich regions of avian genomes are missing from genome assembly models. Instead, it is the presence of tandem repeats containing motifs capable of assembling into very stable secondary structures that is likely responsible.

Comparative analysis of the ribosomal DNA repeat unit (rDNA) of Perna viridis (Linnaeus, 1758) and Perna canaliculus (Gmelin, 1791)

Perna viridis and P. canaliculus are economically and ecologically important species of shellfish. In this study, the complete ribosomal DNA (rDNA) unit sequences of these species were determined for the first time. The gene order, 18S rRNA-internal transcribed spacer (ITS) 1-5.8S rRNA-ITS2-28S rRNA-intergenic spacer (IGS), was similar to that observed in other eukaryotes. The lengths of the P. viridis and P. canaliculus rDNA sequences ranged from 8,432 to 8,616 bp and from 7,597 to 7,610 bp, respectively, this variability was mainly attributable to the IGS region. The putative transcription termination site and initiation site were confirmed. Perna viridis and P. canaliculus rDNA contained two (length: 93 and 40 bp) and one (length: 131 bp) repeat motifs, respectively. Individual intra-species differences mainly involved the copy number of repeat units. In P. viridis, three cytosine-guanine (CpG) sites with sizes of 440, 1,075 and 537 bp were found to cover nearly the entire IGS sequence, whereas in P. canaliculus, two CpG islands with sizes of 361 and 484 bp were identified. The phylogenetic trees constructed with maximum likelihood and neighbour-joining methods and based on ITS sequences were identical and included three major clusters. Species of the same genus were easily clustered together.

Evolution of ribosomal internal transcribed spacers in Deuterostomia

Sequences of ribosomal internal transcribed spacers (ITSs) are of great importance to molecular phylogenetics and DNA barcoding, but remain unstudied in some large taxa of Deuterostomia. We have analyzed complete ITS1 and ITS2 sequences in 62 species from 16 Deuterostomia classes, with ITS sequences in 24 species from 11 classes initially obtained using unannotated contigs and raw read sequences. A general tendency for both ITS length and GC-content increase from interior to superior Deuterostomia taxa, a uniform GC-content in both ITSs within the same species, thymine content decrease in sense DNA sequences of both ITSs are shown. A possible role of GC-based gene conversion in Deuterostomia ITS evolutionary changes is hypothesized. The first example of non-LTR retrotransposon insertion into ITS sequence in Deuterostomia is described in turtle Geochelone nigra. The roles of mobile genetic element insertions in the evolution of ITS sequences in some Sauropsida taxa are discussed as well.

Comparative Cytogenetics between Two Important Songbird, Models: The Zebra Finch and the Canary

Songbird species (order Passeriformes, suborder Oscines) are important models in various experimental fields spanning behavioural genomics to neurobiology. Although the genomes of some songbird species were sequenced recently, the chromosomal organization of these species is mostly unknown. Here we focused on the two most studied songbird species in neuroscience, the zebra finch (Taeniopygia guttata) and the canary (Serinus canaria). In order to clarify these issues and also to integrate chromosome data with their assembled genomes, we used classical and molecular cytogenetics in both zebra finch and canary to define their chromosomal homology, localization of heterochromatic blocks and distribution of rDNA clusters. We confirmed the same diploid number (2n = 80) in both species, as previously reported. FISH experiments confirmed the occurrence of multiple paracentric and pericentric inversions previously found in other species of Passeriformes, providing a cytogenetic signature for this order, and corroborating data from in silico analyses. Additionally, compared to other Passeriformes, we detected differences in the zebra finch karyotype concerning the morphology of some chromosomes, in the distribution of 5S rDNA clusters, and an inversion in chromosome 1.

Ribosomal RNA gene functioning in avian oogenesis

Despite long-term exploration into ribosomal RNA gene functioning during the oogenesis of various organisms, many intriguing problems remain unsolved. In this review, we describe nucleolus organizer region (NOR) activity in avian oocytes. Whereas oocytes from an adult avian ovary never reveal the formation of the nucleolus in the germinal vesicle (GV), an ovary from juvenile birds possesses both nucleolus-containing and non-nucleolus-containing oocytes. The evolutionary diversity of oocyte NOR functioning and the potential non-rRNA-related functions of the nucleolus in oocytes are also discussed.