Molecular evidence for genetic exchanges among ribosomal genes on nonhomologous chromosomes in man and apes

Functional conservation of specialized ribosomes bearing genome-encoded variant rRNAs in Vibrio species

Heterogeneity of ribosomal RNA (rRNA) sequences has recently emerged as a mechanism that can lead to subpopulations of specialized ribosomes. Our previous study showed that ribosomes containing highly divergent rRNAs expressed from the rrnI operon (I-ribosomes) can preferentially translate a subset of mRNAs such as hspA and tpiA in the Vibrio vulnificus CMCP6 strain. Here, we explored the functional conservation of I-ribosomes across Vibrio species. Exogenous expression of the rrnI operon in another V. vulnificus strain, MO6-24/O, and in another Vibrio species, V. fischeri (strain MJ11), decreased heat shock susceptibility by upregulating HspA expression. In addition, we provide direct evidence for the preferential synthesis of HspA by I-ribosomes in the V. vulnificus MO6-24/O strain. Furthermore, exogenous expression of rrnI in V. vulnificus MO6-24/O cells led to higher mortality of infected mice when compared to the wild-type (WT) strain and a strain expressing exogenous rrnG, a redundant rRNA gene in the V. vulnificus CMCP6 strain. Our findings suggest that specialized ribosomes bearing heterogeneous rRNAs play a conserved role in translational regulation among Vibrio species. This study shows the functional importance of rRNA heterogeneity in gene expression control by preferential translation of specific mRNAs, providing another layer of specialized ribosome system.

The 18S rRNA genes of Haemoproteus (Haemosporida, Apicomplexa) parasites from European songbirds with remarks on improved parasite diagnostics

BACKGROUND

The nuclear ribosomal RNA genes of Plasmodium parasites are assumed to evolve according to a birth-and-death model with new variants originating by duplication and others becoming deleted. For some Plasmodium species, it has been shown that distinct variants of the 18S rRNA genes are expressed differentially in vertebrate hosts and mosquito vectors. The central aim was to evaluate whether avian haemosporidian parasites of the genus Haemoproteus also have substantially distinct 18S variants, focusing on lineages belonging to the Haemoproteus majoris and Haemoproteus belopolskyi species groups.

METHODS

The almost complete 18S rRNA genes of 19 Haemoproteus lineages of the subgenus Parahaemoproteus, which are common in passeriform birds from the Palaearctic, were sequenced. The PCR products of 20 blood and tissue samples containing 19 parasite lineages were subjected to molecular cloning, and ten clones in mean were sequenced each. The sequence features were analysed and phylogenetic trees were calculated, including sequence data published previously from eight additional Parahaemoproteus lineages. The geographic and host distribution of all 27 lineages was visualised as CytB haplotype networks and pie charts. Based on the 18S sequence data, species-specific oligonucleotide probes were designed to target the parasites in host tissue by in situ hybridization assays.

RESULTS

Most Haemoproteus lineages had two or more variants of the 18S gene like many Plasmodium species, but the maximum distances between variants were generally lower. Moreover, unlike in most mammalian and avian Plasmodium species, the 18S sequences of all but one parasite lineage clustered into reciprocally monophyletic clades. Considerably distinct 18S clusters were only found in Haemoproteus tartakovskyi hSISKIN1 and Haemoproteus sp. hROFI1. The presence of chimeric 18S variants in some Haemoproteus lineages indicates that their ribosomal units rather evolve in a semi-concerted fashion than according to a strict model of birth-and-death evolution.

CONCLUSIONS

Parasites of the subgenus Parahaemoproteus contain distinct 18S variants, but the intraspecific variability is lower than in most mammalian and avian Plasmodium species. The new 18S data provides a basis for more thorough investigations on the development of Haemoproteus parasites in host tissue using in situ hybridization techniques targeting specific parasite lineages.

On the Roles of the Nuclear Non-Coding RNA-Dependent Membrane-Less Organelles in the Cellular Stress Response

At the beginning of the 21st century, it became obvious that radical changes had taken place in the concept of living matter and, in particular, in the concept of the organization of intracellular space. The accumulated data testify to the essential importance of phase transitions of biopolymers (first of all, intrinsically disordered proteins and RNA) in the spatiotemporal organization of the intracellular space. Of particular interest is the stress-induced reorganization of the intracellular space. Examples of organelles formed in response to stress are nuclear A-bodies and nuclear stress bodies. The formation of these organelles is based on liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) and non-coding RNA. Despite their overlapping composition and similar mechanism of formation, these organelles have different functional activities and physical properties. In this review, we will focus our attention on these membrane-less organelles (MLOs) and describe their functions, structure, and mechanism of formation.

Genetic characteristics and growth patterns of the hybrid grouper derived from the hybridization of Epinephelus fuscoguttatus (female) × Epinephelus polyphekadion (male)

Hybridization is one of the primary methods used to cultivate farmed grouper species. The hybrid grouper derived from crossing Epinephelus fuscoguttatus (♀) and E. polyphekadion (♂) exhibits growth superiority over its parents. The genetic characteristics and growth patterns of the hybrid grouper have not yet been defined. This study confirms the ploidy level of the hybrid grouper (2n = 48) using chromosome count analysis and flow cytometry. The 5S rDNA family was used to evaluate genetic diversity. Only one 5S class (~400 bp) was detected in the hybrid grouper, which could be used to distinguish between two different types based on nucleotide sequences, likely representing homologous unit classes from the female and male parental species. Growth patterns of 5-8-month-old hybrid groupers were also monitored. In this phase, a positive allometric growth pattern in body mass with total length was found. Body height and body mass were significantly correlated based on correlation and path coefficient, suggesting that body height could serve as an excellent index to increase body mass. These results aid our understanding of the genetic evolution of the hybrid grouper and inform the development of improved rearing techniques.

Transcript-specific selective translation by specialized ribosomes bearing genome-encoded heterogeneous rRNAs in V. vulnificus CMCP6

Ribosomes composed of genome-encoded heterogeneous rRNAs are implicated in the rapid adaptation of bacterial cells to environmental changes. A previous study showed that ribosomes bearing the most heterogeneous rRNAs expressed from the rrnI operon (I-ribosomes) are implicated in the preferential translation of a subset of mRNAs, including hspA and tpiA, in Vibrio vulnificus CMCP6. In this study, we show that HspA nascent peptides were predominantly bound to I-ribosomes. Specifically, I-ribosomes were enriched more than two-fold in ribosomes that were pulled down by immunoprecipitation of HspA peptides compared with the proportion of I-ribosomes in crude ribosomes and ribosomes pulled down by immunoprecipitation of RNA polymerase subunit β peptides in the wild-type (WT) and rrnI-completed strains. Other methods that utilized the incorporation of an affinity tag in 23S rRNA or chimeric rRNA tethering 16S and 23S rRNAs, which generated specialized functional ribosomes in Escherichia coli, did not result in functional I-ribosomes in V. vulnificus CMCP6. This study provides direct evidence of the preferential translation of hspA mRNA by I-ribosomes.

5S Ribosomal DNA of Genus Solanum: Molecular Organization, Evolution, and Taxonomy

The Solanum genus, being one of the largest among high plants, is distributed worldwide and comprises about 1,200 species. The genus includes numerous agronomically important species such as Solanum tuberosum (potato), Solanum lycopersicum (tomato), and Solanum melongena (eggplant) as well as medical and ornamental plants. The huge Solanum genus is a convenient model for research in the field of molecular evolution and structural and functional genomics. Clear knowledge of evolutionary relationships in the Solanum genus is required to increase the effectiveness of breeding programs, but the phylogeny of the genus is still not fully understood. The rapidly evolving intergenic spacer region (IGS) of 5S rDNA has been successfully used for inferring interspecific relationships in several groups of angiosperms. Here, combining cloning and sequencing with bioinformatic analysis of genomic data available in the SRA database, we evaluate the molecular organization and diversity of IGS for 184 accessions, representing 137 species of the Solanum genus. It was found that the main mechanisms of IGS molecular evolution was step-wise accumulation of single base substitution or short indels, and that long indels and multiple base substitutions, which arose repeatedly during evolution, were mostly not conserved and eliminated. The reason for this negative selection seems to be association between indels/multiple base substitutions and pseudogenization of 5S rDNA. Comparison of IGS sequences allowed us to reconstruct the phylogeny of the Solanum genus. The obtained dendrograms are mainly congruent with published data: same major and minor clades were found. However, relationships between these clades and position of some species (S. cochoae, S. clivorum, S. macrocarpon, and S. spirale) were different from those of previous results and require further clarification. Our results show that 5S IGS represents a convenient molecular marker for phylogenetic studies on the Solanum genus. In particular, the simultaneous presence of several structural variants of rDNA in the genome enables the detection of reticular evolution, especially in the largest and economically most important sect. Petota. The origin of several polyploid species should be reconsidered.

Specialised ribosomes as versatile regulators of gene expression

The ribosome has long been thought to be a homogeneous cellular machine that constitutively and globally synthesises proteins from mRNA. However, recent studies have revealed that ribosomes are highly heterogeneous, dynamic macromolecular complexes with specialised roles in translational regulation in many organisms across the kingdoms. In this review, we summarise the current understanding of ribosome heterogeneity and the specialised functions of heterogeneous ribosomes. We also discuss specialised translation systems that utilise orthogonal ribosomes.

DNA Environment of Centromeres and Non-Homologous Chromosomes Interactions in Mouse

Although the pericentromeric regions of chromosomes that are enriched in tandemly repeated satellite DNA represent a significant part of eukaryotic genomes, they remain understudied, which is mainly due to interdisciplinary knowledge gaps. Recent studies suggest their important role in genome regulation, karyotype stability, and evolution. Thus, the idea of satellite DNA as a junk part of the genome has been refuted. The integration of data regarding molecular composition, chromosome behaviour, and the details of the in situ organization of pericentromeric regions is of great interest. The objective of this work was a cytogenetic analysis of the interactions between pericentromeric regions from non-homologous chromosomes in mouse spermatocytes using immuno-FISH. We analysed two events: the associations between centromeric regions of the X chromosome and autosomes and the associations between the centromeric regions of the autosomal bivalents that form chromocenters. We concluded that the X chromosome forms temporary synaptic associations with different autosomes in early meiotic prophase I, which can normally be found until the pachytene-diplotene, without signs of pachytene arrest. These associations are formed between the satellite-DNA-rich centromeric regions of the X chromosome and different autosomes but do not involve the satellite-DNA-poor centromeric region of the Y chromosome. We suggest the hypothetical model of X chromosome competitive replacement from such associations during synaptic correction. We showed that the centromeric region of the X chromosome in association remains free of γH2Ax-dependent chromatin inactivation, while the Y chromosome is completely inactivated. This finding highlights the predominant role of associations between satellite DNA-rich regions of different chromosomes, including the X chromosome. We suppose that X-autosomal transient associations are a manifestation of an additional synaptic disorder checkpoint. These associations are normally corrected before the late diplotene stage. We revealed that the intense spreading conditions that were applied to the spermatocyte I nuclei did not lead to the destruction of stretched chromatin fibers of elongated chromocenters enriched in satellite DNA. The tight associations that we revealed between the pericentromeric regions of different autosomal bivalents and the X chromosome may represent the basis for a mechanism for maintaining the repeats stability in the autosomes and in the X chromosome. The consequences of our findings are discussed.

Concerted evolution of ribosomal DNA: Somatic peace amid germinal strife: Intranuclear and cellular selection maintain the quality of rRNA

Most eukaryotes possess many copies of rDNA. Organismal selection alone cannot maintain rRNA function because the effects of mutations in one rDNA are diluted by the presence of many other rDNAs. rRNA quality is maintained by processes that increase homogeneity of rRNA within, and heterogeneity among, germ cells thereby increasing the effectiveness of cellular selection on ribosomal function. A successful rDNA repeat will possess adaptations for spreading within tandem arrays by intranuclear selection. These adaptations reside in the non-coding regions of rDNA. Single-copy genes are predicted to manage processes of intranuclear and cellular selection in the germline to maintain the quality of rRNA expressed in somatic cells of future generations.

Regulatory roles of nucleolus organizer region-derived long non-coding RNAs

The nucleolus is the largest sub-nuclear domain, serving primarily as the place for ribosome biogenesis. A delicately regulated function of the nucleolus is vital to the cell not only for maintaining proper protein synthesis but is also tightly associated with responses to different types of cellular stresses. Recently, several long non-coding RNAs (lncRNAs) were found to be part of the regulatory network that modulate nucleolar functions. Several of these lncRNAs are encoded in the ribosomal DNA (rDNA) repeats or are transcribed from the genomic regions that are located near the nucleolus organizer regions (NORs). In this review, we first discuss the current understanding of the sequence of the NORs and variations between different NORs. We then focus on the NOR-derived lncRNAs in mammalian cells and their functions in rRNA transcription and the organization of nucleolar structure under different cellular conditions. The identification of these lncRNAs reveals great potential of the NORs in harboring novel genes involved in the regulation of nucleolar functions.

The role of internal transcribed spacer 2 secondary structures in classifying mycoparasitic Ampelomyces

Many fungi require specific growth conditions before they can be identified. Direct environmental DNA sequencing is advantageous, although for some taxa, specific primers need to be used for successful amplification of molecular markers. The internal transcribed spacer region is the preferred DNA barcode for fungi. However, inter- and intra-specific distances in ITS sequences highly vary among some fungal groups; consequently, it is not a solely reliable tool for species delineation. Ampelomyces, mycoparasites of the fungal phytopathogen order Erysiphales, can have ITS genetic differences up to 15%; this may lead to misidentification with other closely related unknown fungi. Indeed, Ampelomyces were initially misidentified as other pycnidial mycoparasites, but subsequent research showed that they differ in pycnidia morphology and culture characteristics. We investigated whether the ITS2 nucleotide content and secondary structure was different between Ampelomyces ITS2 sequences and those unrelated to this genus. To this end, we retrieved all ITS sequences referred to as Ampelomyces from the GenBank database. This analysis revealed that fungal ITS environmental DNA sequences are still being deposited in the database under the name Ampelomyces, but they do not belong to this genus. We also detected variations in the conserved hybridization model of the ITS2 proximal 5.8S and 28S stem from two Ampelomyces strains. Moreover, we suggested for the first time that pseudogenes form in the ITS region of this mycoparasite. A phylogenetic analysis based on ITS2 sequences-structures grouped the environmental sequences of putative Ampelomyces into a different clade from the Ampelomyces-containing clades. Indeed, when conducting ITS2 analysis, resolution of genetic distances between Ampelomyces and those putative Ampelomyces improved. Each clade represented a distinct consensus ITS2 S2, which suggested that different pre-ribosomal RNA (pre-rRNA) processes occur across different lineages. This study recommends the use of ITS2 S2s as an important tool to analyse environmental sequencing and unveiling the underlying evolutionary processes.

Horizontally Acquired nrDNAs Persist in Low Amounts in Host Hordeum Genomes and Evolve Independently of Native nrDNA

Nuclear ribosomal DNA (nrDNA) has displayed extraordinary dynamics during the evolution of plant species. However, the patterns and evolutionary significance of nrDNA array expansion or contraction are still relatively unknown. Moreover, only little is known of the fate of minority nrDNA copies acquired between species via horizontal transfer. The barley genus Hordeum (Poaceae) represents a good model for such a study, as species of section Stenostachys acquired nrDNA via horizontal transfer from at least five different panicoid genera, causing long-term co-existence of native (Hordeum-like) and non-native (panicoid) nrDNAs. Using quantitative PCR, we investigated copy number variation (CNV) of nrDNA in the diploid representatives of the genus Hordeum. We estimated the copy number of the foreign, as well as of the native ITS types (ribotypes), and followed the pattern of their CNV in relation to the genus' phylogeny, species' genomes size and the number of nrDNA loci. For the native ribotype, we encountered an almost 19-fold variation in the mean copy number among the taxa analysed, ranging from 1689 copies (per 2C content) in H. patagonicum subsp. mustersii to 31342 copies in H. murinum subsp. glaucum. The copy numbers did not correlate with any of the genus' phylogeny, the species' genome size or the number of nrDNA loci. The CNV was high within the recognised groups (up to 13.2 × in the American I-genome species) as well as between accessions of the same species (up to 4×). Foreign ribotypes represent only a small fraction of the total number of nrDNA copies. Their copy numbers ranged from single units to tens and rarely hundreds of copies. They amounted, on average, to between 0.1% (Setaria ribotype) and 1.9% (Euclasta ribotype) of total nrDNA. None of the foreign ribotypes showed significant differences with respect to phylogenetic groups recognised within the sect. Stenostachys. Overall, no correlation was found between copy numbers of native and foreign nrDNAs suggesting the sequestration and independent evolution of native and non-native nrDNA arrays. Therefore, foreign nrDNA in Hordeum likely poses a dead-end by-product of horizontal gene transfer events.

Evidence of two deeply divergent co-existing mitochondrial genomes in the Tuatara reveals an extremely complex genomic organization

Animal mitochondrial genomic polymorphism occurs as low-level mitochondrial heteroplasmy and deeply divergent co-existing molecules. The latter is rare, known only in bivalvian mollusks. Here we show two deeply divergent co-existing mt-genomes in a vertebrate through genomic sequencing of the Tuatara (Sphenodon punctatus), the sole-representative of an ancient reptilian Order. The two molecules, revealed using a combination of short-read and long-read sequencing technologies, differ by 10.4% nucleotide divergence. A single long-read covers an entire mt-molecule for both strands. Phylogenetic analyses suggest a 7-8 million-year divergence between genomes. Contrary to earlier reports, all 37 genes typical of animal mitochondria, with drastic gene rearrangements, are confirmed for both mt-genomes. Also unique to vertebrates, concerted evolution drives three near-identical putative Control Region non-coding blocks. Evidence of positive selection at sites linked to metabolically important transmembrane regions of encoded proteins suggests these two mt-genomes may confer an adaptive advantage for an unusually cold-tolerant reptile.

Parallel Concerted Evolution of Ribosomal Protein Genes in Fungi and Its Adaptive Significance

Ribosomal protein (RP) genes encode structural components of ribosomes, the cellular machinery for protein synthesis. A single functional copy has been maintained in most of 78–80 RP families in animals due to evolutionary constraints imposed by gene dosage balance. Some fungal species have maintained duplicate copies in most RP families. The mechanisms by which the RP genes were duplicated and maintained and their functional significance are poorly understood. To address these questions, we identified all RP genes from 295 fungi and inferred the timing and nature of gene duplication events for all RP families. We found that massive duplications of RP genes have independently occurred by different mechanisms in three distantly related lineages: budding yeasts, fission yeasts, and Mucoromycota. The RP gene duplicates in budding yeasts and Mucoromycota were mainly created by whole genome duplication events. However, duplicate RP genes in fission yeasts were likely generated by retroposition, which is unexpected considering their dosage sensitivity. The sequences of most RP paralogs have been homogenized by repeated gene conversion in each species, demonstrating parallel concerted evolution, which might have facilitated the retention of their duplicates. Transcriptomic data suggest that the duplication and retention of RP genes increased their transcript abundance. Physiological data indicate that increased ribosome biogenesis allowed these organisms to rapidly consume sugars through fermentation while maintaining high growth rates, providing selective advantages to these species in sugar-rich environments.

The Helvella corium species aggregate in Nordic countries - phylogeny and species delimitation

Mycologists have always been curious about the elaborate morphotypes and shapes of species of the genus Helvella. The small, black, cupulate Helvella specimens have mostly been assigned to Helvella corium, a broadly defined morpho-species. Recent phylogenetic analyses, however, have revealed an aggregate of species hidden under this name. We performed a multispecies coalescent analysis to re-assess species limits and evolutionary relationships of the Helvella corium species aggregate in the Nordic countries. To achieve this, we used morphology and phylogenetic evidence from five loci – heat shock protein 90 (hsp), translation elongation factor 1-alpha (tef), RNA polymerase II (rpb2), and the 5.8S and large subunit (LSU) of the nuclear ribosomal DNA. All specimens under the name Helvella corium in the larger university fungaria of Norway, Sweden and Denmark were examined and barcoded, using partial hsp and/or rpb2 as the preferential secondary barcodes in Helvella. Additional fresh specimens were collected in three years (2015–2018) to obtain in vivo morphological data to aid in species discrimination. The H. corium species aggregate consists of seven phylogenetically distinct species, nested in three divergent lineages, i.e. H. corium, H. alpina and H. pseudoalpina sp. nov. in the /alpina-corium lineage, H. alpestris, H. macrosperma and H. nannfeldtii in the /alpestris-nannfeldtii lineage, and H. alpicola as a weakly supported sister to the /alpestris-nannfeldtii lineage. Among the seven species, the ribosomal loci expressed substantial variation in evolutionary rates, suggesting care in the use of these regions alone in delimitation of Helvella species. Altogether, 469 out of 496 available fungarium specimens were successfully barcoded.

Chromosome Dynamics Regulating Genomic Dispersion and Alteration of Nucleolus Organizer Regions (NORs)

The nucleolus organizer regions (NORs) demonstrate differences in genomic dispersion and transcriptional activity among all organisms. I postulate that such differences stem from distinct genomic structures and their interactions from chromosome observations using fluorescence in situ hybridization and silver nitrate staining methods. Examples in primates and Australian bulldog ants indicate that chromosomal features indeed play a significant role in determining the properties of NORs. In primates, rDNA arrays that are located on the short arm of acrocentrics frequently form reciprocal associations ("affinity"), but they lack such associations ("non-affinity") with other repeat arrays-a binary molecular effect. These "rules" of affinity vs. non-affinity are extrapolated from the chromosomal configurations of meiotic prophase. In bulldog ants, genomic dispersions of rDNA loci expand much more widely following an increase in the number of acrocentric chromosomes formed by centric fission. Affinity appears to be a significantly greater force: associations likely form among rDNA and heterochromatin arrays of acrocentrics-thus, more acrocentrics bring about more rDNA loci. The specific interactions among NOR-related genome structures remain unclear and require further investigation. Here, I propose that there are limited and non-limited genomic dispersion systems that result from genomic affinity rules, inducing specific chromosomal configurations that are related to NORs.

Widespread prevalence but contrasting patterns of intragenomic rRNA polymorphisms in nematodes: Implications for phylogeny, species delimitation and life history inference

Ribosomal RNA genes have long been a favoured locus in phylogenetic and metabarcoding studies. Within a genome, rRNA loci are organized as tandem repeated arrays and the copies are homogenized through the process of concerted evolution. However, some level of rRNA variation (intragenomic polymorphism) is known to persist and be maintained in the genomes of many species. In nematode worms, the extent of rRNA polymorphism (RP) across species and the evolutionary and life history factors that contribute to the maintenance of intragenomic RP is largely unknown. Here, we present an extensive analysis across 30 terrestrial nematode species representing a range of free-living and parasitic taxa isolated worldwide. Our results indicate that RP is common and widespread, ribosome function appears to be maintained despite mutational changes, and intragenomic variants are stable in the genome and neutrally evolving. However, levels of variation were varied widely across rRNA locus and species, with some taxa observed to lack RP entirely. Higher levels of RP were significantly correlated with shorter generation time and high reproductive rates, and population-level factors may play a role in the geographic and phylogenetic structuring of rRNA variants observed in genera such as Rotylenchulus and Pratylenchus. Although RP did not dramatically impact the clustering and recovery of taxa in mock metabarcoding analyses, the present study has significant implications for global biodiversity estimates of nematode species derived from environmental rRNA amplicon studies, as well as our understanding of the evolutionary and ecological factors shaping genetic diversity across the nematode Tree of Life.

The nuclear 18S ribosomal DNAs of avian haemosporidian parasites

Background

Plasmodium species feature only four to eight nuclear ribosomal units on different chromosomes, which are assumed to evolve independently according to a birth-and-death model, in which new variants originate by duplication and others are deleted throughout time. Moreover, distinct ribosomal units were shown to be expressed during different developmental stages in the vertebrate and mosquito hosts. Here, the 18S rDNA sequences of 32 species of avian haemosporidian parasites are reported and compared to those of simian and rodent Plasmodium species.

Methods

Almost the entire 18S rDNAs of avian haemosporidians belonging to the genera Plasmodium (7), Haemoproteus (9), and Leucocytozoon (16) were obtained by PCR, molecular cloning, and sequencing ten clones each. Phylogenetic trees were calculated and sequence patterns were analysed and compared to those of simian and rodent malaria species. A section of the mitochondrial CytB was also sequenced.

Results

Sequence patterns in most avian Plasmodium species were similar to those in the mammalian parasites with most species featuring two distinct 18S rDNA sequence clusters. Distinct 18S variants were also found in Haemoproteus tartakovskyi and the three Leucocytozoon species, whereas the other species featured sets of similar haplotypes. The 18S rDNA GC-contents of the Leucocytozoon toddi complex and the subgenus Parahaemoproteus were extremely high with 49.3% and 44.9%, respectively. The 18S sequences of several species from all three genera showed chimeric features, thus indicating recombination.

Conclusion

Gene duplication events leading to two diverged main sequence clusters happened independently in at least six out of seven avian Plasmodium species, thus supporting evolution according to a birth-and-death model like proposed for the ribosomal units of simian and rodent Plasmodium species. Patterns were similar in the 18S rDNAs of the Leucocytozoon toddi complex and Haemoproteus tartakovskyi. However, the 18S rDNAs of the other species seem to evolve in concerted fashion like in most eukaryotes, but the presence of chimeric variants indicates that the ribosomal units rather evolve in a semi-concerted manner. The new data may provide a basis for studies testing whether differential expression of distinct 18S rDNA also occurs in avian Plasmodium species and related haemosporidian parasites.

Structural variations of subterminal satellite blocks and their source mechanisms as inferred from the meiotic configurations of chimpanzee chromosome termini

African great apes have large constitutive heterochromatin (C-band) blocks in subtelomeric regions of the majority of their chromosomes, but humans lack these. Additionally, the chimpanzee meiotic cell division process demonstrates unique partial terminal associations in the first meiotic prophase (pachytene). These are likely formed as a result of interaction among subtelomeric C-band blocks. We thus conducted an extensive study to define the features in the subtelomeric heterochromatic regions of chimpanzee chromosomes undergoing mitotic metaphase and meiotic cell division. Molecular cytogenetic analyses with probes of both subterminal satellite DNA (a main component of C-band) and rDNA demonstrated principles of interaction among DNA arrays. The results suggest that homologous and ectopic recombination through persistent subtelomeric associations (post-bouquet association observed in 32% of spermatocytes in the pachytene stage) appears to create variability in heterochromatin patterns and simultaneously restrain subtelomeric genome polymorphisms. That is, the meeting of non-homologous chromosome termini sets the stage for ectopic pairing which, in turn, is the mechanism for generating variability and genomic dispersion of subtelomeric C-band blocks through a system of concerted evolution. Comparison between the present study and previous reports indicated that the chromosomal distribution rate of sutelomeric regions seems to have antagonistic correlation with arm numbers holding subterminal satellite blocks in humans, chimpanzees, and gorillas. That is, the increase of subterminal satellite blocks probably reduces genomic diversity in the subtelomeric regions. The acquisition vs. loss of the subtelomeric C-band blocks is postulated as the underlying engine of this chromosomal differentiation yielded by meiotic chromosomal interaction.

Analysis of genetic diversity and phylogeny of Philosamia ricini (Lepidoptera: Saturniidae) by using RAPD and internal transcribed spacer DNA1

The Indian Eri silkworm, Philosamia ricini Hutt, a commercial silk producing insect, is believed to have originated in the Brahmaputra valley of Assam. In this study, the genetic diversity and phylogenetic relationships of six morphs of Eri silkworm viz. white plain, white zebra, white spotted, blue plain, blue zebra and blue spotted collected from different geographical locations of North-East India were investigated by using random amplified polymorphic DNA (RAPD) and the first internal transcribed spacer region (ITS1). This study revealed a low genetic diversity among the morphs of Eri silkworm. Twenty-eight random primers generated 199 bands. Out of these, 112 were polymorphic (56.28%) with an average of 7.1 bands per primer. The genetic similarity matrix ranged from 0.56 to 0.99. The morphs collected from same geographical area shared the same cluster in the dendrogram. The genetic diversity in case of ITS1 sequences (2.19%) was found to be less as compared to RAPD. The ITS1 sequences of the morphs collected from same geographical area were found to be identical. The information generated in this study will help in conservation and effective breeding program to improve its productivity.

Implications of sequence variation on the evolution of rRNA

The evolution of the multi-copy family of ribosomal RNA (rRNA) genes is unique in regard to its genetics and genome evolution. Paradoxically, rRNA genes are highly homogenized within and between individuals, yet they are globally distinct between species. Here, we discuss the implications for models of rRNA gene evolution in light of our recent discoveries that ribosomes bearing rRNA sequence variants can affect gene expression and physiology and that intra-individual rRNA alleles exhibit both context- and tissue-specific expression.

Comparative analysis of C-heterochromatin, ribosomal and telomeric DNA markers in chromosomes of Pamphagidae grasshoppers from Morocco

The karyotypes and the localization of C-bands, clusters of ribosomal DNA and telomeric repeats of 10 species of the family Pamphagidae from Morocco are described for the first time. The species studied belong to the subfamilies Pamphaginae and Thrinchinae. All species have karyotypes consisting of 19 and 20 acrocentric chromosomes and X0/XX sex chromosome system in males and females, respectively (2n♂=19, NF=19; 2n♀=20, NF=20). Despite the karyotype conservatism, we revealed differences in the location and size of C-heterochromatin blocks and ribosomal DNA clusters. A comparative analysis of these differences shows that karyotype divergences in this group is connected not to structural chromosome rearrangements, but to the evolution of repetitive DNA.

Variation in human chromosome 21 ribosomal RNA genes characterized by TAR cloning and long-read sequencing

Despite the key role of the human ribosome in protein biosynthesis, little is known about the extent of sequence variation in ribosomal DNA (rDNA) or its pre-rRNA and rRNA products. We recovered ribosomal DNA segments from a single human chromosome 21 using transformation-associated recombination (TAR) cloning in yeast. Accurate long-read sequencing of 13 isolates covering ∼0.82 Mb of the chromosome 21 rDNA complement revealed substantial variation among tandem repeat rDNA copies, several palindromic structures and potential errors in the previous reference sequence. These clones revealed 101 variant positions in the 45S transcription unit and 235 in the intergenic spacer sequence. Approximately 60% of the 45S variants were confirmed in independent whole-genome or RNA-seq data, with 47 of these further observed in mature 18S/28S rRNA sequences. TAR cloning and long-read sequencing enabled the accurate reconstruction of multiple rDNA units and a new, high-quality 44 838 bp rDNA reference sequence, which we have annotated with variants detected from chromosome 21 of a single individual. The large number of variants observed reveal heterogeneity in human rDNA, opening up the possibility of corresponding variations in ribosome dynamics.

Coumarin Content, Morphological Variation, and Molecular Phylogenetics of Melilotus

Melilotus albus and Melilotus officinalis are widely used in forage production and herbal medicine due to the biological activity of their coumarins, which have many biological and pharmacological activities, including anti-HIV and anti-tumor effects. To comprehensively evaluate M. albus and M. officinalis coumarin content (Cou), morphological variation, and molecular phylogeny, we examined the Cou, five morphological traits and the molecular characterization based on the trnL-F spacer and internal transcribed spacer (ITS) regions of 93 accessions. Significant (p < 0.05) variation was observed in the Cou and all five morphological traits in both species. Analysis of population differentiation (Pst) of the phenotypic traits showed that powdery mildew resistance (PMR) had the greatest Pst, meaning that this trait demonstrated the largest genetic differentiation among the accessions. The Pst values of dry matter yield (DMY) and Cou were relatively high. Biplot analysis identified accessions with higher DMY and higher and lower Cou. Analysis of molecular sequence variation identified seven haplotypes of the trnL-F spacer and 13 haplotypes of the ITS region. Based on haplotype and sequence analyses, the genetic variation of M. officinalis was higher than that of M. albus. Additionally, ITS sequence analysis showed that the variation among accessions was larger than that among species across three geographical areas: Asia, Europe, and North America. Similarly, variation among accessions for both the trnL-F and ITS sequences were larger than the differences between the geographical areas. Our results indicate that there has been considerable gene flow between the two Melilotus species. Our characterization of Cou and the morphological and genetic variations of these two Melilotus species may provide useful insights into germplasm improvement to enhance DMY and Cou.

Taxonomy of Cyrtochilum-alliance (Orchidaceae) in the light of molecular and morphological data

BACKGROUND

The generic separateness and specific composition of the orchid genus Cyrtochilum was discussed for almost two centuries. Over the years several smaller taxa were segregated from this taxon, but their separateness was recently questioned based on molecular studies outcomes. The aim of our study was to revise concepts of morphological-based generic delimitation in Cyrtochilum-alliance and to compare it with the results of genetic analysis. We used phylogenetic framework in combination with phenetical analysis to provide proposal of the generic delimitation within Cyrtochilum-alliance. Two molecular markers, ITS and matK were used to construct phylogenetic tree. A total of over 5000 herbarium specimens were included in the morphological examination and the phenetical analysis included 29 generative and vegetative characters.

RESULTS

Comparative morphology of the previously recognized genera: Buesiella, Dasyglossum, Neodryas, Rusbyella, Siederella and Trigonochilum is presented. A new species within the the latter genus is described. Fourteen new combinations are proposed. The key to the identification of the genera of the Cyrtochilum-alliance and morphological characteristics of each genus are provided.

CONCLUSIONS

A total of six separated genera are recognized within Cyrtochilum-alliance. The reasons of the incompatibility between morphological differences observed within studied taxa and phylogenetic tree are argued and the taxonomic implications of such inconsistency, resulting in fragmentation or lumping of taxonomic units, are discussed.

Hybridization between two cryptic filamentous brown seaweeds along the shore: analysing pre- and postzygotic barriers in populations of individuals with varying ploidy levels

We aimed to study the importance of hybridization between two cryptic species of the genus Ectocarpus, a group of filamentous algae with haploid-diploid life cycles that include the principal genetic model organism for the brown algae. In haploid-diploid species, the genetic structure of the two phases of the life cycle can be analysed separately in natural populations. Such life cycles provide a unique opportunity to estimate the frequency of hybrid genotypes in diploid sporophytes and meiotic recombinant genotypes in haploid gametophytes allowing the effects of reproductive barriers preventing fertilization or preventing meiosis to be untangle. The level of hybridization between E. siliculosus and E. crouaniorum was quantified along the European coast. Clonal cultures (568 diploid, 336 haploid) isolated from field samples were genotyped using cytoplasmic and nuclear markers to estimate the frequency of hybrid genotypes in diploids and recombinant haploids. We identified admixed individuals using microsatellite loci, classical assignment methods and a newly developed Bayesian method (XPloidAssignment), which allows the analysis of populations that exhibit variations in ploidy level. Over all populations, the level of hybridization was estimated at 8.7%. Hybrids were exclusively observed in sympatric populations. More than 98% of hybrids were diploids (40% of which showed signs of aneuploidy) with a high frequency of rare alleles. The near absence of haploid recombinant hybrids demonstrates that the reproductive barriers are mostly postzygotic and suggests that abnormal chromosome segregation during meiosis following hybridization of species with different genome sizes could be a major cause of interspecific incompatibility in this system.

Epistatic and allelic interactions control expression of ribosomal RNA gene clusters in Arabidopsis thaliana

BACKGROUND

Ribosomal RNA (rRNA) accounts for the majority of the RNA in eukaryotic cells, and is encoded by hundreds to thousands of nearly identical gene copies, only a subset of which are active at any given time. In Arabidopsis thaliana, 45S rRNA genes are found in two large ribosomal DNA (rDNA) clusters and little is known about the contribution of each to the overall transcription pattern in the species.

RESULTS

By taking advantage of genome sequencing data from the 1001 Genomes Consortium, we characterize rRNA gene sequence variation within and among accessions. Notably, variation is not restricted to the pre-rRNA sequences removed during processing, but it is also present within the highly conserved ribosomal subunits. Through linkage mapping we assign these variants to a particular rDNA cluster unambiguously and use them as reporters of rDNA cluster-specific expression. We demonstrate that rDNA cluster-usage varies greatly among accessions and that rDNA cluster-specific expression and silencing is controlled via genetic interactions between entire rDNA cluster haplotypes (alleles).

CONCLUSIONS

We show that rRNA gene cluster expression is controlled via complex epistatic and allelic interactions between rDNA haplotypes that apparently regulate the entire rRNA gene cluster. Furthermore, the sequence polymorphism we discovered implies that the pool of rRNA in a cell may be heterogeneous, which could have functional consequences.

Secondary structure of nrDNA Internal Transcribed Spacers as a useful tool to align highly divergent species in phylogenetic studies

Recently, it has been suggested that internal transcribed spacer (ITS) sequences are under selective constraints to preserve their secondary structure. Here, we investigate the patterns of the ITS nucleotide and secondary structure conservation across the Passiflora L. genus to evaluate the potential use of secondary structure data as a helpful tool for the alignment in taxonomically complex genera. Considering the frequent use of ITS, this study also presents a perspective on future analyses in other plant groups. The ITS1 and ITS2 sequences presented significant differences for mean values of the lowest energy state (LES) and for number of hairpins in different Passiflora subgenera. Statistical analyses for the subgenera separately support significant differences between the LES values and the total number of secondary structures for ITS. In order to evaluate whether the LES values of ITS secondary structures were related to selective constraints, we compared these results among 120 ITS sequences from Passiflora species and 120 randomly generated sequences. These analyses indicated that Passiflora ITS sequences present characteristics of a region under selective constraint to maintain the secondary structure showing to be a promising tool to improve the alignments and identify sites with non-neutral substitutions or those correlated evolutionary steps.

Evolutional dynamics of 45S and 5S ribosomal DNA in ancient allohexaploid Atropa belladonna

BACKGROUND

Polyploid hybrids represent a rich natural resource to study molecular evolution of plant genes and genomes. Here, we applied a combination of karyological and molecular methods to investigate chromosomal structure, molecular organization and evolution of ribosomal DNA (rDNA) in nightshade, Atropa belladonna (fam. Solanaceae), one of the oldest known allohexaploids among flowering plants. Because of their abundance and specific molecular organization (evolutionarily conserved coding regions linked to variable intergenic spacers, IGS), 45S and 5S rDNA are widely used in plant taxonomic and evolutionary studies.

RESULTS

Molecular cloning and nucleotide sequencing of A. belladonna 45S rDNA repeats revealed a general structure characteristic of other Solanaceae species, and a very high sequence similarity of two length variants, with the only difference in number of short IGS subrepeats. These results combined with the detection of three pairs of 45S rDNA loci on separate chromosomes, presumably inherited from both tetraploid and diploid ancestor species, example intensive sequence homogenization that led to substitution/elimination of rDNA repeats of one parent. Chromosome silver-staining revealed that only four out of six 45S rDNA sites are frequently transcriptionally active, demonstrating nucleolar dominance. For 5S rDNA, three size variants of repeats were detected, with the major class represented by repeats containing all functional IGS elements required for transcription, the intermediate size repeats containing partially deleted IGS sequences, and the short 5S repeats containing severe defects both in the IGS and coding sequences. While shorter variants demonstrate increased rate of based substitution, probably in their transition into pseudogenes, the functional 5S rDNA variants are nearly identical at the sequence level, pointing to their origin from a single parental species. Localization of the 5S rDNA genes on two chromosome pairs further supports uniparental inheritance from the tetraploid progenitor.

CONCLUSIONS

The obtained molecular, cytogenetic and phylogenetic data demonstrate complex evolutionary dynamics of rDNA loci in allohexaploid species of Atropa belladonna. The high level of sequence unification revealed in 45S and 5S rDNA loci of this ancient hybrid species have been seemingly achieved by different molecular mechanisms.

Interstitial NORs, Fragile Sites, and Chromosome Evolution: A Not So Simple Relationship - The Example of Melolontha melolontha and Genus Protaetia (Coleoptera: Scarabaeidae)

In the present study, the origin of recurrent rearrangements involving chromosome 6 in 3.2% of cells of Melolontha melolontha (Coleoptera, Scarabaeidae) was investigated. Various chromosome staining techniques, including C-banding, Giemsa and silver staining, as well as fluorescence in situ hybridization with a human 28S rDNA probe, were applied to M. melolontha chromosome spreads. In addition, related species of the genera Melolontha and Protaetia were studied. On chromosome 6 of M. melolontha, there is a fragile site-like structure which corresponds to an interstitial nucleolus organizer region (NOR). Despite this instability, the NOR remains unique and interstitial in this species, as well as in the other species studied. It is proposed that the intercalary position of the NOR both facilitates the detection of its fragile site-like instability and correlates with its relative stability during evolution. We explain this apparent paradox by strong counter-selection for imbalances of the chromosome fragment distal to the interstitial NORs, which would recurrently occur in the progeny of translocation carriers. Thus, the frequent telomeric position of the NORs in most animal and plant taxa would have no functional rationale but would be the consequence of selection against the meiotic transmission of chromosome imbalances.

The Antarctic Circumpolar Current as a diversification trigger for deep-sea octocorals

BACKGROUND

Antarctica is surrounded by the Antarctic Circumpolar Current (ACC), the largest and strongest current in the world. Despite its potential importance for shaping biogeographical patterns, the distribution and connectivity of deep-sea populations across the ACC remain poorly understood. In this study we conducted the first assessment of phylogeographical patterns in deep-sea octocorals in the South Pacific and Southern Ocean, specifically a group of closely related bottlebrush octocorals (Primnoidae: Tokoprymno and Thourella), as a test case to study the effect of the ACC on the population structure of brooding species. We assessed the degree to which the ACC constitutes a barrier to gene flow between northern and southern populations and whether the onset of diversification of these corals coincides with the origin of the ACC (Oligocene-Miocene boundary).

RESULTS

Based on DNA sequences of two nuclear genes from 80 individuals and a combination of phylogeographic model-testing approaches we found a phylogenetic break corresponding to the spatial occurrence of the ACC. We also found significant genetic structure among our four regional populations. However, we uncovered shared haplotypes among certain population pairs, suggesting long-distance, asymmetrical migration. Our divergence time analyses indicated that the separation of amphi-ACC populations took place during the Middle Miocene around 12.6 million years ago, i.e., after the formation of the ACC.

CONCLUSION

We suggest that the ACC constitutes a semi-permeable barrier to these deep-sea octocorals capable of separating and structuring populations, while allowing short periods of gene flow. The fluctuations in latitudinal positioning of the ACC during the Miocene likely contributed to the diversification of these octocorals. Additionally, we provide evidence that the populations from each of our four sampling regions could actually constitute different species.

Molecular Differentiation of the African Yellow Fever Vector Aedes bromeliae (Diptera: Culicidae) from Its Sympatric Non-vector Sister Species, Aedes lilii

INTRODUCTION

Yellow fever continues to be a problem in sub-Saharan Africa with repeated epidemics occurring. The mosquito Aedes bromeliae is a major vector of yellow fever, but it cannot be readily differentiated from its non-vector zoophilic sister species Ae. lilii using morphological characters. Genetic differences have been reported between anthropophilic Ae. bromeliae and zoophilic Ae. lilii and between forest and domestic populations. However, due to the application of different molecular markers and non-overlapping populations employed in previous studies, interpretation of species delimitation is unclear.

METHODOLOGY/PRINCIPLE FINDINGS

DNA sequences were generated from specimens of Ae. simpsoni s.l. from the Republic of Benin, Tanzania and Uganda for two nuclear genes apolipophorin 2 (apoLp2) and cytochrome p450 (CYPJ92), the ribosomal internal transcribed spacer region (ITS) and the mitochondrial cytochrome c oxidase (COI) barcoding region. Nuclear genes apoLp2 and CYPJ92 were unable to differentiate between species Ae. bromeliae and Ae. lilii due to ancestral lineage sorting, while ITS sequence data provided clear topological separation on a phylogeny. The standard COI barcoding region was shown to be subject to species introgression and unable to clearly distinguish the two taxa. Here we present a reliable direct PCR-based method for differentiation of the vector species Ae. bromeliae from its isomorphic, sympatric and non-biomedically important sister taxon, Ae. lilii, based on the ITS region. Using molecular species verification, we describe novel immature habitats for Ae. lilii and report both sympatric and allopatric populations. Whereas only Ae. lilii is found in the Republic of Benin and only Ae. bromeliae in Tanzania, both species are sympatric in Uganda.

CONCLUSIONS/SIGNIFICANCE

Our accurate identification method will allow informed distribution and detailed ecological studies that will facilitate assessment of arboviral disease risk and development of future targeted vector control.

Patterns of Gene Conversion in Duplicated Yeast Histones Suggest Strong Selection on a Coadapted Macromolecular Complex

We find evidence for interlocus gene conversion in five duplicated histone genes from six yeast species. The sequences of these duplicated genes, surviving from the ancient genome duplication, show phylogenetic patterns inconsistent with the well-resolved orthology relationships inferred from a likelihood model of gene loss after the genome duplication. Instead, these paralogous genes are more closely related to each other than any is to its nearest ortholog. In addition to simulations supporting gene conversion, we also present evidence for elevated rates of radical amino acid substitutions along the branches implicated in the conversion events. As these patterns are similar to those seen in ribosomal proteins that have undergone gene conversion, we speculate that in cases where duplicated genes code for proteins that are a part of tightly interacting complexes, selection may favor the fixation of gene conversion events in order to maintain high protein identities between duplicated copies.

Rapid identification of Helicoverpa armigera and Helicoverpa zea (Lepidoptera: Noctuidae) using ribosomal RNA internal transcribed spacer 1

Rapid identification of invasive species is crucial for deploying management strategies to prevent establishment. Recent Helicoverpa armigera (Hübner) invasions and subsequent establishment in South America has increased the risk of this species invading North America. Morphological similarities make differentiation of H. armigera from the native Helicoverpa zea (Boddie) difficult. Characteristics of adult male genitalia and nucleotide sequence differences in mitochondrial DNA are two of the currently available methods to differentiate these two species. However, current methods are likely too slow to be employed as rapid detection methods. In this study, conserved differences in the internal transcribed spacer 1 (ITS1) of the ribosomal RNA genes were used to develop species-specific oligonucleotide primers that amplified ITS1 fragments of 147 and 334 bp from H. armigera and H. zea, respectively. An amplicon (83 bp) from a conserved region of 18S ribosomal RNA subunit served as a positive control. Melting temperature differences in ITS1 amplicons yielded species-specific dissociation curves that could be used in high resolution melt analysis to differentiate the two Helicoverpa species. In addition, a rapid and inexpensive procedure for obtaining amplifiable genomic DNA from a small amount of tissue was identified. Under optimal conditions, the process was able to detect DNA from one H. armigera leg in a pool of 25 legs. The high resolution melt analysis combined with rapid DNA extraction could be used as an inexpensive method to genetically differentiate large numbers of H. armigera and H. zea using readily available reagents.

Rare Events of Intragenus and Intraspecies Horizontal Transfer of the 16S rRNA Gene

Horizontal gene transfer (HGT) of operational genes has been widely reported in prokaryotic organisms. However, informational genes such as those involved in transcription and translation processes are very difficult to be horizontally transferred, as described by Woese’s complexity hypothesis. Here, we analyzed all of the completed prokaryotic genome sequences (2,143 genomes) in the NCBI (National Center for Biotechnology Information) database, scanned for genomes with high intragenomic heterogeneity of 16S rRNA gene copies, and explored potential HGT events of ribosomal RNA genes based on the phylogeny, genomic organization, and secondary structures of the ribosomal RNA genes. Our results revealed 28 genomes with relatively high intragenomic heterogeneity of multiple 16S rRNA gene copies (lowest pairwise identity <98.0%), and further analysis revealed HGT events and potential donors of the heterogeneous copies (such as HGT from Chlamydia suis to Chlamydia trachomatis) and mutation events of some heterogeneous copies (such as Streptococcus suis JS14). Interestingly, HGT of the 16S rRNA gene only occurred at intragenus or intraspecies levels, which is quite different from the HGT of operational genes. Our results improve our understanding regarding the exchange of informational genes.

Interspecific Phylogenic Relationships within Genus Melilotus Based on Nuclear and Chloroplast DNA

Melilotus comprises 19 species, while the phylogenetic relationships between species remain unclear. In the present work, three chloroplast genes, rbcL, matK, trnL-F, and one nuclear region, ITS (internal transcribed spacer) belonging to 48 populations of 18 species of Melilotus were sequenced and phylogenetic trees were constructed to study their interspecific relationships. Based on the phylogenetic tree generated in this study using rbcL analysis, the Melilotus genus is clearly monophyletic in the legume family. Both Bayesian and maximum-parsimony approaches were used to analyze the data. The nrDNA ITS provided more informative characteristics (9.8%) than cpDNA (3.0%). Melilotus contains two closely related groups, clade I and clade II. M. spicatus, M. indicus and M. segetalis have a close relationship. M. infestus, M. siculus and M. sulcatus are closely related. The comparing between molecular phylogeny and flower color classification in Melilotus showed that the flower color is not much informative for phylogenetics of this genus.

Interlocus gene conversion explains at least 2.7% of single nucleotide variants in human segmental duplications

Background

Interlocus gene conversion (IGC) is a recombination-based mechanism that results in the unidirectional transfer of short stretches of sequence between paralogous loci. Although IGC is a well-established mechanism of human disease, the extent to which this mutagenic process has shaped overall patterns of segregating variation in multi-copy regions of the human genome remains unknown. One expected manifestation of IGC in population genomic data is the presence of one-to-one paralogous SNPs that segregate identical alleles.

Results

Here, I use SNP genotype calls from the low-coverage phase 3 release of the 1000 Genomes Project to identify 15,790 parallel, shared SNPs in duplicated regions of the human genome. My approach for identifying these sites accounts for the potential redundancy of short read mapping in multi-copy genomic regions, thereby effectively eliminating false positive SNP calls arising from paralogous sequence variation. I demonstrate that independent mutation events to identical nucleotides at paralogous sites are not a significant source of shared polymorphisms in the human genome, consistent with the interpretation that these sites are the outcome of historical IGC events. These putative signals of IGC are enriched in genomic contexts previously associated with non-allelic homologous recombination, including clear signals in gene families that form tandem intra-chromosomal clusters.

Conclusions

Taken together, my analyses implicate IGC, not point mutation, as the mechanism generating at least 2.7 % of single nucleotide variants in duplicated regions of the human genome.

Electronic supplementary material

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

Concerted copy number variation balances ribosomal DNA dosage in human and mouse genomes

Tandemly repeated ribosomal DNA (rDNA) arrays are among the most evolutionary dynamic loci of eukaryotic genomes. The loci code for essential cellular components, yet exhibit extensive copy number (CN) variation within and between species. CN might be partly determined by the requirement of dosage balance between the 5S and 45S rDNA arrays. The arrays are nonhomologous, physically unlinked in mammals, and encode functionally interdependent RNA components of the ribosome. Here we show that the 5S and 45S rDNA arrays exhibit concerted CN variation (cCNV). Despite 5S and 45S rDNA elements residing on different chromosomes and lacking sequence similarity, cCNV between these loci is strong, evolutionarily conserved in humans and mice, and manifested across individual genotypes in natural populations and pedigrees. Finally, we observe that bisphenol A induces rapid and parallel modulation of 5S and 45S rDNA CN. Our observations reveal a novel mode of genome variation, indicate that natural selection contributed to the evolution and conservation of cCNV, and support the hypothesis that 5S CN is partly determined by the requirement of dosage balance with the 45S rDNA array. We suggest that human disease variation might be traced to disrupted rDNA dosage balance in the genome.

Mitochondrial and nuclear phylogenetic analysis with Sanger and next-generation sequencing shows that, in Área de Conservación Guanacaste, northwestern Costa Rica, the skipper butterfly named Urbanus belli (family Hesperiidae) comprises three morphologically cryptic species

Background

Skipper butterflies (Hesperiidae) are a relatively well-studied family of Lepidoptera. However, a combination of DNA barcodes, morphology, and natural history data has revealed several cryptic species complexes within them. Here, we investigate three DNA barcode lineages of what has been identified as Urbanus belli (Hesperiidae, Eudaminae) in Área de Conservación Guanacaste (ACG), northwestern Costa Rica.

Results

Although no morphological traits appear to distinguish among the three, congruent nuclear and mitochondrial lineage patterns show that “Urbanus belli” in ACG is a complex of three sympatric species. A single strain of Wolbachia present in two of the three cryptic species indicates that Urbanus segnestami Burns (formerly Urbanus belliDHJ01), Urbanus bernikerni Burns (formerly Urbanus belliDHJ02), and Urbanus ehakernae Burns (formerly Urbanus belliDHJ03) may be biologically separated by Wolbachia, as well as by their genetics. Use of parallel sequencing through 454-pyrosequencing improved the utility of ITS2 as a phylogenetic marker and permitted examination of the intra- and interlineage relationships of ITS2 variants within the species complex. Interlineage, intralineage and intragenomic compensatory base pair changes were discovered in the secondary structure of ITS2.

Conclusion

These findings corroborate the existence of three cryptic species. Our confirmation of a novel cryptic species complex, initially suggested by DNA barcode lineages, argues for using a multi-marker approach coupled with next-generation sequencing for exploration of other suspected species complexes.

Molecular organization of the 25S-18S rDNA IGS of Fagus sylvatica and Quercus suber: a comparative analysis

The 35S ribosomal DNA (rDNA) units, repeated in tandem at one or more chromosomal loci, are separated by an intergenic spacer (IGS) containing functional elements involved in the regulation of transcription of downstream rRNA genes. In the present work, we have compared the IGS molecular organizations in two divergent species of Fagaceae, Fagus sylvatica and Quercus suber, aiming to comprehend the evolution of the IGS sequences within the family. Self- and cross-hybridization FISH was done on representative species of the Fagaceae. The IGS length variability and the methylation level of 18 and 25S rRNA genes were assessed in representatives of three genera of this family: Fagus, Quercus and Castanea. The intergenic spacers in Beech and Cork Oak showed similar overall organizations comprising putative functional elements needed for rRNA gene activity and containing a non-transcribed spacer (NTS), a promoter region, and a 5′-external transcribed spacer. In the NTS: the sub-repeats structure in Beech is more organized than in Cork Oak, sharing some short motifs which results in the lowest sequence similarity of the entire IGS; the AT-rich region differed in both spacers by a GC-rich block inserted in Cork Oak. The 5′-ETS is the region with the higher similarity, having nonetheless different lengths. FISH with the NTS-5′-ETS revealed fainter signals in cross-hybridization in agreement with the divergence between genera. The diversity of IGS lengths revealed variants from ∼2 kb in Fagus, and Quercus up to 5.3 kb in Castanea, and a lack of correlation between the number of variants and the number of rDNA loci in several species. Methylation of 25S Bam HI site was confirmed in all species and detected for the first time in the 18S of Q. suber and Q. faginea. These results provide important clues for the evolutionary trends of the rDNA 25S-18S IGS in the Fagaceae family.

Robertsonian chromosomes and the nuclear architecture of mouse meiotic prophase spermatocytes

Background

The nuclear architecture of meiotic prophase spermatocytes is based on higher-order patterns of spatial associations among chromosomal domains from different bivalents. The meiotic nuclear architecture depends on the chromosome characteristics and consequently is prone to modification by chromosomal rearrangements. In this work, we consider Mus domesticus spermatocytes with diploid chromosome number 2n = 40, all telocentric, and investigate a possible modification of the ancestral nuclear architecture due to the emergence of derived Rb chromosomes, which may be present in the homozygous or heterozygous condition.

Results

In the 2n = 40 spermatocyte nuclei random associations mediated by pericentromeric heterochromatin among the 19 telocentric bivalents ocurr at the nuclear periphery. The observed frequency of associations among them, made distinguishable by specific probes and FISH, seems to be the same for pairs that may or may not form Rb chromosomes. In the homozygote Rb 2n = 24 spermatocytes, associations also mediated by pericentromeric heterochromatin occur mainly between the three telocentric or the eight metacentric bivalents themselves. In heterozygote Rb 2n = 32 spermatocytes all heterochromatin is localized at the nuclear periphery, yet associations are mainly observed among the three telocentric bivalents and between the asynaptic axes of the trivalents.

Conclusions

The Rb chromosomes pose sharp restrictions for interactions in the 2n = 24 and 2n = 32 spermatocytes, as compared to the ample possibilities for interactions between bivalents in the 2n = 40 spermatocytes. Undoubtedly the emergence of Rb chromosomes changes the ancestral nuclear architecture of 2n = 40 spermatocytes since they establish new types of interactions among chromosomal domains, particularly through centromeric and heterochromatic regions at the nuclear periphery among telocentric and at the nuclear center among Rb metacentric ones.