Ribosomal DNA in the grasshopper Podisma pedestris: escape from concerted evolution

Structure and Evolution of Ribosomal Genes of Insect Chromosomes

Currently, clusters of 45S and 5S ribosomal DNA (rDNA) have been studied in about 1000 and 100 species of the class Insecta, respectively. Although the number of insect species with known 45S rDNA clusters (also referred to as nucleolus-organizing regions, or NORs) constitutes less than 0.1 percent of the described members of this enormous group, certain conclusions can already be drawn. Since haploid karyotypes with single 45S and 5S rDNA clusters predominate in both basal and derived insect groups, this character state is apparently ancestral for the class Insecta in general. Nevertheless, the number, chromosomal location, and other characteristics of both 45S and 5S rDNA sites substantially vary across different species, and sometimes even within the same species. There are several main factors and molecular mechanisms that either maintain these parameters or alter them on the short-term and/or long-term scale. Chromosome structure (i.e., monocentric vs. holokinetic chromosomes), excessive numbers of rRNA gene copies per cluster, interactions with transposable elements, pseudogenization, and meiotic recombination are perhaps the most important among them.

Intragenomic rDNA variation - the product of concerted evolution, mutation, or something in between?

The classical model of concerted evolution states that hundreds to thousands of ribosomal DNA (rDNA) units undergo homogenization, making the multiple copies of the individual units more uniform across the genome than would be expected given mutation frequencies and gene redundancy. While the universality of this over 50-year-old model has been confirmed in a range of organisms, advanced high throughput sequencing techniques have also revealed that rDNA homogenization in many organisms is partial and, in rare cases, even apparently failing. The potential underpinning processes leading to unexpected intragenomic variation have been discussed in a number of studies, but a comprehensive understanding remains to be determined. In this work, we summarize information on variation or polymorphisms in rDNAs across a wide range of taxa amongst animals, fungi, plants, and protists. We discuss the definition and description of concerted evolution and describe whether incomplete concerted evolution of rDNAs predominantly affects coding or non-coding regions of rDNA units and if it leads to the formation of pseudogenes or not. We also discuss the factors contributing to rDNA variation, such as interspecific hybridization, meiotic cycles, rDNA expression status, genome size, and the activity of effector genes involved in genetic recombination, epigenetic modifications, and DNA editing. Finally, we argue that a combination of approaches is needed to target genetic and epigenetic phenomena influencing incomplete concerted evolution, to give a comprehensive understanding of the evolution and functional consequences of intragenomic variation in rDNA.

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.

Guess who? Taxonomic problems in the genus Eiseniella revisited by integrated approach

Eiseniella neapolitana is a semi-aquatic, diploid earthworm that for many years was related to the cosmopolitan species Eiseniella tetraedra and even considered a subspecies of it. Norealidys andaluciana was described in Spain and is usually synonymized with E. neapolitana. We collected 69 specimens from Italy, Spain, and Cyprus and studied five molecular markers (COI, 16S, 28S, 12S, and ND1) and their morphology to solve this taxonomic problem. Phylogenetic analyses reveal the possible existence of two separate genera confounded under the name Eiseniella, but the study of more molecular markers and species of the genus would be necessary to confirm this. Therefore, the synonymy between Eiseniella and Norealidys is maintained. Various genetic analyses, including species delimitation, confirm the separation between E. neapolitana and E. andaluciana (= N. andaluciana) and excluded that E. neapolitana is a subspecies of E. tetraedra. The resemblance in external appearance despite clear genetic differences of the three species could be explained by convergent adaptation to the aquatic habitat. Despite the expected low haplotype diversity based on the 28S gene, we found a surprisingly high variability in the E. andaluciana (= N. andaluciana) population in Spain. However, its stable predicted secondary structure and its high content of G + C reject the presence of a pseudogene.

Varying strength of selection contributes to the intragenomic diversity of rRNA genes

Ribosome biogenesis in eukaryotes is supported by hundreds of ribosomal RNA (rRNA) gene copies that are encoded in the ribosomal DNA (rDNA). The multiple copies of rRNA genes are thought to have low sequence diversity within one species. Here, we present species-wide rDNA sequence analysis in Saccharomyces cerevisiae that challenges this view. We show that rDNA copies in this yeast are heterogeneous, both among and within isolates, and that many variants avoided fixation or elimination over evolutionary time. The sequence diversity landscape across the rDNA shows clear functional stratification, suggesting different copy-number thresholds for selection that contribute to rDNA diversity. Notably, nucleotide variants in the most conserved rDNA regions are sufficiently deleterious to exhibit signatures of purifying selection even when present in only a small fraction of rRNA gene copies. Our results portray a complex evolutionary landscape that shapes rDNA sequence diversity within a single species and reveal unexpectedly strong purifying selection of multi-copy genes.

Polymorphism of 16s rRNA Gene: Any Effect on the Biomolecular Quantitation of the Honey Bee (Apis mellifera L., 1758) Pathogen Nosema ceranae?

The microsporidian Nosema ceranae is a severe threat to the western honey bee Apis mellifera, as it is responsible for nosemosis type C, which leads the colonies to dwindle and collapse. Infection quantification is essential to clinical and research aims. Assessment is made often with molecular assays based on rRNA genes, which are present in the N. ceranae genome as multiple and polymorphic copies. This study aims to compare two different methods of Real-Time PCR (qPCR), respectively relying on the 16S rRNA and Hsp70 genes, the first of which is described as a multiple and polymorphic gene. Young worker bees, hatched in the laboratory and artificially inoculated with N. ceranae spores, were incubated at 33 °C and subject to different treatment regimens. Samples were taken post-infection and analyzed with both qPCR methods. Compared to Hsp70, the 16S rRNA method systematically detected higher abundance. Straightforward conversion between the two methods is made impossible by erratic 16s rRNA/Hsp70 ratios. The 16s rRNA polymorphism showed an increase around the inoculated dose, where a higher prevalence of ungerminated spores was expected due to the treatment effects. The possible genetic background of that irregular distribution is discussed in detail. The polymorphic nature of 16S rRNA showed to be a limit in the infection quantification. More reliably, the N. ceranae abundance can be assessed in honey bee samples with methods based on the single-copy gene Hsp70.

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. .

Intragenomic sequence variations in the second internal transcribed spacer (ITS2) ribosomal DNA of the malaria vector Anopheles stephensi

Second Internal Transcribed Spacer (ITS2) ribosomal DNA (rDNA) sequence is a widely used molecular marker for species-identification or -delimitation due to observed concerted evolution which is believed to homogenize rDNA copies in an interbreeding population. However, intra-specific differences in ITS2 of Anopheles stephensi have been reported. This study reports the presence of intragenomic sequence variation in the ITS2-rDNA of An. stephensi and hypothesizes that observed intra-specific differences in this species may have resulted due to ambiguous DNA sequence-chromatogram resulting from intragenomic heterogeneity. Anopheles stephensi collected from different parts of India were sequenced for complete ITS2 and the variable region of 28S-rDNA (d1-d3 domains). Intragenomic variations were found in ITS2 region of all An. stephensi sequenced, but no such variation was observed in d1 to d3 domains of 28S-rDNA. Cloning and sequencing of ITS2 through the d3 domain of the 28S region of rDNA from representative samples from northern, central, and southern India confirmed the presence of intragenomic variation in ITS2 due to transitions at three loci and two bp indel in a di-nucleotide microsatellite locus. Multiple haplotypes were observed in ITS2 raised from such variations. Due to the absence of detectable intragenomic sequence variation in the d1 to d3 domain of 28S rDNA of An. stephensi, this region can serve as an ideal reference sequence for taxonomic and phylogenetic studies. The presence of intragenomic variation in rDNA should be carefully examined before using this as a molecular marker for species delimitation or phylogenetic analyses.

It Is Just a Matter of Time: Balancing Homologous Recombination and Non-homologous End Joining at the rDNA Locus During Meiosis

Ribosomal RNA genes (rDNAs) are located in large domains of hundreds of rDNA units organized in a head-to-tail manner. The proper and stable inheritance of rDNA clusters is of paramount importance for survival. Yet, these highly repetitive elements pose a potential risk to the genome since they can undergo non-allelic exchanges. Here, we review the current knowledge of the organization of the rDNA clusters in Arabidopsis thaliana and their stability during meiosis. Recent findings suggest that during meiosis, all rDNA loci are embedded within the nucleolus favoring non-homologous end joining (NHEJ) as a repair mechanism, while DNA repair via homologous recombination (HR) appears to be a rare event. We propose a model where (1) frequent meiotic NHEJ events generate abundant single nucleotide polymorphisms and insertions/deletions within the rDNA, resulting in a heterogeneous population of rDNA units and (2) rare HR events dynamically change rDNA unit numbers, only to be observed in large populations over many generations. Based on the latest efforts to delineate the entire rDNA sequence in A. thaliana, we discuss evidence supporting this model. The results compiled so far draw a surprising picture of rDNA sequence heterogeneity between individual units. Furthermore, rDNA cluster sizes have been recognized as relatively stable when observing less than 10 generations, yet emerged as major determinant of genome size variation between different A. thaliana ecotypes. The sequencing efforts also revealed that transcripts from the diverse rDNA units yield heterogenous ribosome populations with potential functional implications. These findings strongly motivate further research to understand the mechanisms that maintain the metastable state of rDNA loci.

Intragenomic internal transcribed spacer 2 variation in a genus of parasitoid wasps (Hymenoptera: Braconidae): implications for accurate species delimitation and phylogenetic analysis

A recent DNA barcoding study of Australian microgastrines (Hymenoptera: Braconidae) sought to use next-generation sequencing of the cytochrome c oxidase subunit 1 (COI) barcoding gene region, the wingless (WG) gene and the internal transcribed spacer 2 (ITS2) to delimit molecular species in a highly diverse group of parasitic wasps. Large intragenomic distances between ITS2 variants, often larger than the average interspecific variation, caused difficulties in using ITS2 for species delimitation in both threshold and tree-based approaches, and the gene was not included in the reported results of the previous DNA barcoding study. We here report on the intragenomic, and the intra- and interspecies, variation in ITS2in the microgastrine genus Diolcogasterto further investigate the value of ITS2as a marker for species delimitation and phylogenetics of the Microgastrinae. Distinctive intragenomic variant patterns were found in different species of Diolcogaster, with some species possessing a single major variant, and others possessing many divergent variants. Characterizing intragenomic variation of ITS2is critical as it is a widely used marker in hymenopteran phylogenetics and species delimitation, and large intragenomic distances such as those found in this study may obscure phylogenetic signal.

The Birth and Death of Toxins with Distinct Functions: A Case Study in the Sea Anemone Nematostella

The cnidarian Nematostella vectensis has become an established lab model, providing unique opportunities for venom evolution research. The Nematostella venom system is multimodal: involving both nematocytes and ectodermal gland cells, which produce a toxin mixture whose composition changes throughout the life cycle. Additionally, their modes of interaction with predators and prey vary between eggs, larvae, and adults, which is likely shaped by the dynamics of the venom system.

Nv1 is a major component of adult venom, with activity against arthropods (through specific inhibition of sodium channel inactivation) and fish. Nv1 is encoded by a cluster of at least 12 nearly identical genes that were proposed to be undergoing concerted evolution. Surprisingly, we found that Nematostella venom includes several Nv1 paralogs escaping a pattern of general concerted evolution, despite belonging to the Nv1-like family. Here, we show two of these new toxins, Nv4 and Nv5, are lethal for zebrafish larvae but harmless to arthropods, unlike Nv1. Furthermore, unlike Nv1, the newly identified toxins are expressed in early life stages. Using transgenesis and immunostaining, we demonstrate that Nv4 and Nv5 are localized to ectodermal gland cells in larvae.

The evolution of Nv4 and Nv5 can be described either as neofunctionalization or as subfunctionalization. Additionally, the Nv1-like family includes several pseudogenes being an example of nonfunctionalization and venom evolution through birth-and-death mechanism. Our findings reveal the evolutionary history for a toxin radiation and point toward the ecological function of the novel toxins constituting a complex cnidarian venom.

Remarkable variation of ribosomal DNA organization and copy number in gnetophytes, a distinct lineage of gymnosperms

INTRODUCTION

Gnetophytes, comprising the genera Ephedra, Gnetum and Welwitschia, are an understudied, enigmatic lineage of gymnosperms with a controversial phylogenetic relationship to other seed plants. Here we examined the organization of ribosomal DNA (rDNA) across representative species.

METHODS

We applied high-throughput sequencing approaches to isolate and reconstruct rDNA units and to determine their intragenomic homogeneity. In addition, fluorescent in situ hybridization and Southern blot hybridization techniques were used to reveal the chromosome and genomic organization of rDNA.

KEY RESULTS

The 5S and 35S rRNA genes were separate (S-type) in Gnetum montanum, Gnetum gnemon and Welwitschia mirabilis and linked (L-type) in Ephedra altissima. There was considerable variability in 5S rDNA abundance, ranging from as few as ~4000 (W. mirabilis) to >100 000 (G. montanum) copies. A similar large variation was also observed in 5S rDNA locus numbers (two to 16 sites per diploid cell). 5S rRNA pseudogenes were interspersed between functional genes forming a single unit in E. altissima and G. montanum. Their copy number was comparable or even higher than that of functional 5S rRNA genes. In E. altissima internal transcribed spacers of 35S rDNA were long and intrinsically repetitive while in G. montanum and W. mirabilis they were short without the subrepeats.

CONCLUSIONS

Gnetophytes are distinct from other gymnosperms and angiosperms as they display surprisingly large variability in rDNA organization and rDNA copy and locus numbers between genera, with no relationship between copy numbers and genome sizes apparent. Concerted evolution of 5S rDNA units seems to have led to the amplification of 5S pseudogenes in both G. montanum and E. altissima. Evolutionary patterns of rDNA show both gymnosperm and angiosperm features underlining the diversity of the group.

Genotyping of individual Ceratonova shasta (Cnidaria: Myxosporea) myxospores reveals intra-spore ITS-1 variation and invalidates the distinction of genotypes II and III

Genotypes of the myxosporean parasite Ceratonova shasta are defined by the number of ATC repeats in the parasite's ribosomal DNA internal transcribed spacer region 1. These genotypes correlate with specific salmonid fish hosts. We observed coho salmon (Oncorhynchus kisutch) and rainbow trout (Oncorhynchus mykiss) with mixtures of genotypes II and III, and assumed that this was a consequence of fish having an aggregate infection from multiple individual parasites. We hypothesized that although multiple ITS copies are present within a parasite spore, the DNA sequences of these copies are identical, and thus individual C. shasta spores are a single genotype. We tested this by extracting and sequencing DNA from individual myxospores. We trialed three approaches for in-tube DNA extraction; digestion with proteinase K was superior to simply rehydrating spores, or incubation in the buffer. Sequences from 14 myxospores were each a mixture of genotypes II and III. Therefore, intra-genomic ribosomal DNA variants exist within individual parasite spores, and II and III should no longer be regarded as discrete C. shasta genotypes. This single-spore genotyping approach will be a useful tool for testing validity of other C. shasta genotypes, and for correctly matching genotype with phenotype for mixed infections of other myxozoan species.

Evolutionary trends in animal ribosomal DNA loci: introduction to a new online database

Ribosomal DNA (rDNA) loci encoding 5S and 45S (18S-5.8S-28S) rRNAs are important components of eukaryotic chromosomes. Here, we set up the animal rDNA database containing cytogenetic information about these loci in 1343 animal species (264 families) collected from 542 publications. The data are based on in situ hybridisation studies (both radioactive and fluorescent) carried out in major groups of vertebrates (fish, reptiles, amphibians, birds, and mammals) and invertebrates (mostly insects and mollusks). The database is accessible online at www.animalrdnadatabase.com . The median number of 45S and 5S sites was close to two per diploid chromosome set for both rDNAs despite large variation (1-74 for 5S and 1-54 for 45S sites). No significant correlation between the number of 5S and 45S rDNA loci was observed, suggesting that their distribution and amplification across the chromosomes follow independent evolutionary trajectories. Each group, irrespective of taxonomic classification, contained rDNA sites at any chromosome location. However, the distal and pericentromeric positions were the most prevalent (> 75% karyotypes) for 45S loci, while the position of 5S loci was more variable. We also examined potential relationships between molecular attributes of rDNA (homogenisation and expression) and cytogenetic parameters such as rDNA positions, chromosome number, and morphology.

Variations in 5S rDNAs in diploid and tetraploid offspring of red crucian carp × common carp

BACKGROUND

The allotetraploid hybrid fish (4nAT) that was created in a previous study through an intergeneric cross between red crucian carp (Carassius auratus red var., ♀) and common carp (Cyprinus carpio L., ♂) provided an excellent platform to investigate the effect of hybridization and polyploidization on the evolution of 5S rDNA. The 5S rDNAs of paternal common carp were made up of a coding sequence (CDS) and a non-transcribed spacer (NTS) unit, and while the 5S rDNAs of maternal red crucian carp contained a CDS and a NTS unit, they also contained a variable number of interposed regions (IPRs). The CDSs of the 5S rDNAs in both parental fishes were conserved, while their NTS units seemed to have been subjected to rapid evolution.

RESULTS

The diploid hybrid 2nF₁ inherited all the types of 5S rDNAs in both progenitors and there were no signs of homeologous recombination in the 5S rDNAs of 2nF₁ by sequencing of PCR products. We obtained two segments of 5S rDNA with a total length of 16,457 bp from allotetraploid offspring 4nAT through bacterial artificial chromosome (BAC) sequencing. Using this sequence together with the 5S rDNA sequences amplified from the genomic DNA of 4nAT, we deduced that the 5S rDNAs of 4nAT might be inherited from the maternal progenitor red crucian carp. Additionally, the IPRs in the 5S rDNAs of 4nAT contained A-repeats and TA-repeats, which was not the case for the IPRs in the 5S rDNAs of 2nF₁. We also detected two signals of a 200-bp fragment of 5S rDNA in the chromosomes of parental progenitors and hybrid progenies by fluorescence in situ hybridization (FISH).

CONCLUSIONS

We deduced that during the evolution of 5S rDNAs in different ploidy hybrid fishes, interlocus gene conversion events and tandem repeat insertion events might occurred in the process of polyploidization. This study provided new insights into the relationship among the evolution of 5S rDNAs, hybridization and polyploidization, which were significant in clarifying the genome evolution of polyploid fish.

Contrasting Patterns of rDNA Homogenization within the Zygosaccharomyces rouxii Species Complex

Arrays of repetitive ribosomal DNA (rDNA) sequences are generally expected to evolve as a coherent family, where repeats within such a family are more similar to each other than to orthologs in related species. The continuous homogenization of repeats within individual genomes is a recombination process termed concerted evolution. Here, we investigated the extent and the direction of concerted evolution in 43 yeast strains of the Zygosaccharomyces rouxii species complex (Z. rouxii, Z. sapae, Z. mellis), by analyzing two portions of the 35S rDNA cistron, namely the D1/D2 domains at the 5’ end of the 26S rRNA gene and the segment including the internal transcribed spacers (ITS) 1 and 2 (ITS regions). We demonstrate that intra-genomic rDNA sequence variation is unusually frequent in this clade and that rDNA arrays in single genomes consist of an intermixing of Z. rouxii, Z. sapae and Z. mellis-like sequences, putatively evolved by reticulate evolutionary events that involved repeated hybridization between lineages. The levels and distribution of sequence polymorphisms vary across rDNA repeats in different individuals, reflecting four patterns of rDNA evolution: I) rDNA repeats that are homogeneous within a genome but are chimeras derived from two parental lineages via recombination: Z. rouxii in the ITS region and Z. sapae in the D1/D2 region; II) intra-genomic rDNA repeats that retain polymorphisms only in ITS regions; III) rDNA repeats that vary only in their D1/D2 domains; IV) heterogeneous rDNA arrays that have both polymorphic ITS and D1/D2 regions. We argue that an ongoing process of homogenization following allodiplodization or incomplete lineage sorting gave rise to divergent evolutionary trajectories in different strains, depending upon temporal, structural and functional constraints. We discuss the consequences of these findings for Zygosaccharomyces species delineation and, more in general, for yeast barcoding.

Variation in the number of nucleoli and incomplete homogenization of 18S ribosomal DNA sequences in leaf cells of the cultivated Oriental ginseng (Panax ginseng Meyer)

BACKGROUND

Wild ginseng, Panax ginseng Meyer, is an endangered species of medicinal plants. In the present study, we analyzed variations within the ribosomal DNA (rDNA) cluster to gain insight into the genetic diversity of the Oriental ginseng, P. ginseng, at artificial plant cultivation.

METHODS

The roots of wild P. ginseng plants were sampled from a nonprotected natural population of the Russian Far East. The slides were prepared from leaf tissues using the squash technique for cytogenetic analysis. The 18S rDNA sequences were cloned and sequenced. The distribution of nucleotide diversity, recombination events, and interspecific phylogenies for the total 18S rDNA sequence data set was also examined.

RESULTS

In mesophyll cells, mononucleolar nuclei were estimated to be dominant (75.7%), while the remaining nuclei contained two to four nucleoli. Among the analyzed 18S rDNA clones, 20% were identical to the 18S rDNA sequence of P. ginseng from Japan, and other clones differed in one to six substitutions. The nucleotide polymorphism was more expressed at the positions 440-640 bp, and distributed in variable regions, expansion segments, and conservative elements of core structure. The phylogenetic analysis confirmed conspecificity of ginseng plants cultivated in different regions, with two fixed mutations between P. ginseng and other species.

CONCLUSION

This study identified the evidences of the intragenomic nucleotide polymorphism in the 18S rDNA sequences of P. ginseng. These data suggest that, in cultivated plants, the observed genome instability may influence the synthesis of biologically active compounds, which are widely used in traditional medicine.

Contrasting evolutionary patterns of 28S and ITS rRNA genes reveal high intragenomic variation in Cephalenchus (Nematoda): Implications for species delimitation

Concerted evolution is often assumed to be the evolutionary force driving multi-family genes, including those from ribosomal DNA (rDNA) repeat, to complete homogenization within a species, although cases of non-concerted evolution have been also documented. In this study, sequence variation of 28S and ITS ribosomal RNA (rRNA) genes in the genus Cephalenchus is assessed at three different levels, intragenomic, intraspecific, and interspecific. The findings suggest that not all Cephalenchus species undergo concerted evolution. High levels of intraspecific polymorphism, mostly due to intragenomic variation, are found in Cephalenchus sp1 (BRA-01). Secondary structure analyses of both rRNA genes and across different species show a similar substitution pattern, including mostly compensatory (CBC) and semi-compensatory (SBC) base changes, thus suggesting the functionality of these rRNA copies despite the variation found in some species. This view is also supported by low sequence variation in the 5.8S gene in relation to the flanking ITS-1 and ITS-2 as well as by the existence of conserved motifs in the former gene. It is suggested that potential cross-fertilization in some Cephalenchus species, based on inspection of female reproductive system, might contribute to both intragenomic and intraspecific polymorphism of their rRNA genes. These results reinforce the potential implications of intragenomic and intraspecific genetic diversity on species delimitation, especially in biodiversity studies based solely on metagenetic approaches. Knowledge of sequence variation will be crucial for accurate species diversity estimation using molecular methods.

Nucleolar dominance and maternal control of 45S rDNA expression

Using a system of interspecies hybrids, trihybrids, and recombinants with varying proportions of genomes from three distinct Xenopus species, we provide evidence for de novo epigenetic silencing of paternal 45 S ribosomal ribonucleic acid (rRNA) genes and their species-dependent expression dominance that escapes transcriptional inactivation after homologous recombination. The same pattern of imprinting is maintained in the offspring from mothers being genetic males (ZZ) sex-reversed to females, indicating that maternal control of ribosomal deoxyribonucleic acid (rDNA) expression is not sex-chromosome linked. Nucleolar dominance (nucleolus underdevelopment) in Xenopus hybrids appears to be associated with a major non-Mendelian reduction in the number of 45 S rDNA gene copies rather than a specific pattern of their expression. The loss of rRNA gene copies in F1 hybrids was non-random with respect to the parental species, with the transcriptionally dominant variant preferentially removed from hybrid zygotes. This dramatic disruption in the structure and function of 45 S rDNA impacts transcriptome patterns of small nucleolar RNAs and messenger RNAs, with genes from the ribosome and oxidative stress pathways being among the most affected. Unorthodoxies of rDNA inheritance and expression may be interpreted as hallmarks of genetic conflicts between parental genomes, as well as defensive epigenetic mechanisms employed to restore genome integrity.

Astonishing 35S rDNA diversity in the gymnosperm species Cycas revoluta Thunb

In all eukaryotes, the highly repeated 35S ribosomal DNA (rDNA) sequences encoding 18S-5.8S-26S ribosomal RNA (rRNA) typically show high levels of intragenomic uniformity due to homogenisation processes, leading to concerted evolution of 35S rDNA repeats. Here, we compared 35S rDNA divergence in several seed plants using next generation sequencing and a range of molecular and cytogenetic approaches. Most species showed similar 35S rDNA homogeneity indicating concerted evolution. However, Cycas revoluta exhibits an extraordinary diversity of rDNA repeats (nucleotide sequence divergence of different copies averaging 12 %), influencing both the coding and non-coding rDNA regions nearly equally. In contrast, its rRNA transcriptome was highly homogeneous suggesting that only a minority of genes (<20 %) encode functional rRNA. The most common SNPs were C > T substitutions located in symmetrical CG and CHG contexts which were also highly methylated. Both functional genes and pseudogenes appear to cluster on chromosomes. The extraordinary high levels of 35S rDNA diversity in C. revoluta, and probably other species of cycads, indicate that the frequency of repeat homogenisation has been much lower in this lineage, compared with all other land plant lineages studied. This has led to the accumulation of methylation-driven mutations and pseudogenisation. Potentially, the reduced homology between paralogs prevented their elimination by homologous recombination, resulting in long-term retention of rDNA pseudogenes in the genome.

A Multiplex Real-Time PCR Assay to Diagnose and Separate Helicoverpa armigera and H. zea (Lepidoptera: Noctuidae) in the New World

The Old World bollworm, Helicoverpa armigera (Hübner), and the corn earworm, H. zea (Boddie), are two of the most important agricultural pests in the world. Diagnosing these two species is difficult-adults can only be separated with a complex dissection, and larvae cannot be identified to species using morphology, necessitating the use of geographic origin for identification in most instances. With the discovery of H. armigera in the New World, identification of immature Helicoverpa based on origin is no longer possible because H. zea also occurs in all of the geographic regions where H. armigera has been discovered. DNA barcoding and restriction fragment length polymorphism (RFLP) analyses have been reported in publications to distinguish these species, but these methods both require post-PCR processing (i.e., DNA sequencing or restriction digestion) to complete. We report the first real-time PCR assay to distinguish these pests based on two hydrolysis probes that bind to a segment of the internal transcribed spacer region 2 (ITS2) amplified using a single primer pair. One probe targets H. armigera, the second probe targets H. zea, and a third probe that targets a conserved segment of 18S rDNA is used as a control of DNA quality. The assay can be completed in 50 minutes when using isolated DNA and is successfully tested on larvae intercepted at ports of entry and adults captured during domestic surveys. We demonstrate that the assay can be run in triplex with no negative effects on sensitivity, can be run using alternative real-time PCR reagents and instruments, and does not cross react with other New World Heliothinae.

Comparison of five methods for delimitating species in Ophion Fabricius, a diverse genus of parasitoid wasps (Hymenoptera, Ichneumonidae)

DNA taxonomy has been proposed as a method to quickly assess diversity and species limits in highly diverse, understudied taxa. Here we use five methods for species delimitation and two genetic markers (COI and ITS2) to assess species diversity within the parasitoid genus, Ophion. We searched for compensatory base changes (CBC's) in ITS2, and determined that they are too rare to be of practical use in delimiting species in this genus. The other four methods used both COI and ITS2, and included distance-based (threshold analysis and ABGD) and tree-based (GMYC and PTP) models. We compared the results of these analyses to each other under various parameters and tested their performance with respect to 11 Nearctic species/morphospecies and 15 described Palearctic species. We also computed barcode accumulation curves of COI sequences to assess the completeness of sampling. The species count was highly variable depending on the method and parameters used, ranging from 47 to 168 species, with more conservative estimates of 89-121 species. Despite this range, many of the Nearctic test species were fairly robust with respect to method. We concluded that while there was often good congruence between methods, GMYC and PTP were less reliant on arbitrary parameters than the other two methods and more easily applied to genetic markers other than COI. However, PTP was less successful at delimiting test species than was GMYC. All methods, as well as the barcode accumulation curves, indicate that several Palearctic species remain undescribed and that we have scarcely begun to appreciate the Nearctic diversity within this genus.

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

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

Extent and divergence of heteroplasmy of the DNA barcoding region in Anapodisma miramae (Orthoptera: Acrididae)

A partial sequence of the mitochondrial cytochrome oxidase subunit I (COI) gene is widely used as a molecular marker for species identification in animals, also termed a DNA barcode. However, the presence of more than one sequence type in a single individual, also known as heteroplasmy, is one of the shortcomings of barcode identification. In this study, we examined the extent and divergence of COI heteroplasmy, including nuclear-encoded mitochondrial pseudogenes (NUMTs), at the genomic-DNA level from 13 insect species including orthopteran Anapodisma miramae, and a long fragment of mitochondrial DNA and cDNA from A. miramae as templates. When multiple numbers of clones originated from genomic DNA were sequenced, heteroplasmy was prevalent in all species and NUMTs were observed in five species. Long fragment DNA (∼13.5 kb) also is a source of heteroplasmic amplification, but the divergent haplotypes and NUMTs obtained from genomic DNA were not detected in A. miramae. On the other hand, cDNA was relatively heteroplasmy-free. Consistently, one dominant haplotype was always obtained from the genomic DNA-origin clones in all species and also from the long fragment- and cDNA-origin clones in the two tested individuals of A. miramae. Furthermore, the dominant haplotype was identical in sequence, regardless of the DNA source in A. miramae. Thus, one possible solution to avoid the barcoding problem in relationship to heteroplasmy could be the acquisition of multiple numbers of barcoding sequences to determine a dominant haplotype that can be assigned as barcoding sequence for a given species.

A nuclear ribosomal DNA pseudogene in triatomines opens a new research field of fundamental and applied implications in Chagas disease

A pseudogene, designated as "ps(5.8S+ITS-2)", paralogous to the 5.8S gene and internal transcribed spacer (ITS)-2 of the nuclear ribosomal DNA (rDNA), has been recently found in many triatomine species distributed throughout North America, Central America and northern South America. Among characteristics used as criteria for pseudogene verification, secondary structures and free energy are highlighted, showing a lower fit between minimum free energy, partition function and centroid structures, although in given cases the fit only appeared to be slightly lower. The unique characteristics of "ps(5.8S+ITS-2)" as a processed or retrotransposed pseudogenic unit of the ghost type are reviewed, with emphasis on its potential functionality compared to the functionality of genes and spacers of the normal rDNA operon. Besides the technical problem of the risk for erroneous sequence results, the usefulness of "ps(5.8S+ITS-2)" for specimen classification, phylogenetic analyses and systematic/taxonomic studies should be highlighted, based on consistence and retention index values, which in pseudogenic sequence trees were higher than in functional sequence trees. Additionally, intraindividual, interpopulational and interspecific differences in pseudogene amount and the fact that it is a pseudogene in the nuclear rDNA suggests a potential relationships with fitness, behaviour and adaptability of triatomine vectors and consequently its potential utility in Chagas disease epidemiology and control.

Evolutionary relationships in the medicinally important genus Fritillaria L. (Liliaceae)

Fritillaria (Liliaceae) is a genus of approximately 140 species of bulbous perennial plants that includes taxa of both horticultural and medicinal importance. As well as being commercially valuable, Fritillaria species have attracted attention because of their exceptionally large genome sizes, with all values recorded to date in excess of 30Gb. Despite such interest in the genus, phylogenetic relationships between the majority of species have remained untested. Here we present the first phylogenetic reconstruction of relationships to encompass most of the currently recognised species diversity in the genus. Three regions of the plastid genome were sequenced in 117 individuals of Fritillaria, representing 92 species (c. 66% of the genus) and in representatives of nine other genera of Liliaceae. Eleven low-copy nuclear gene regions were also screened in selected species for their potential utility. Phylogenetic analysis of a combined plastid dataset using maximum parsimony and Bayesian inference provided support for the monophyly of the majority of currently recognised subgenera. However, subgenus Fritillaria, which is by far the largest of the subgenera and includes the most important species used in traditional Chinese medicine, is found to be polyphyletic. Moreover, several taxa that were represented by multiple individuals show evidence of species non-monophyly. The Japanese endemic subgenus Japonica, which contains the species with the largest recorded genome size for any diploid plant, is resolved as sister to the predominantly Middle Eastern and Central Asian subgenus Rhinopetalum. Whilst relationships between most of the major Fritillaria lineages can now be resolved, our results also highlight the need for data from additional independently evolving loci; an endeavour that may be particularly challenging in light of the huge nuclear genomes found in these plants.

New tools and insights to assist with the molecular identification of Simulium guianense s.l., main Onchocerca volvulus vector within the highland areas of the Amazonia onchocerciasis focus

Following the success of the Onchocerciasis Elimination Programme for the Americas (OEPA), there is now just one Latin American onchocerciasis focus where onchocerciasis transmission is described as 'on-going:' the Amazonia Onchocerciasis focus. In the hyperendemic highland areas of the Amazonia focus, Simulium guianense s.l. Wise are the most important vectors of the disease. Populations of S. guianense s.l. are, however, known to vary in their cytogenetics and in a range of behaviours, including in their biting habits. In the hypoendemic lowland areas of the Amazonia focus, for example, S. guianense s.l. are generally regarded as zoophilic and consequently unimportant to disease transmission. Robust tools, to discriminate among various populations of S. guianense s.l. have, however, not yet been developed. In the work reported here, we have assessed the utility of a ribosomal DNA sequence fragment spanning the nuclear ribosomal ITS-1, ITS-2 and 5.8S sequence regions and a ∼850 nucleotide portion of the mitochondrial cytochrome oxidase gene (CO1) for species-level identification and for resolving the within species substructuring. We report here how we have generated 78 CO1 sequences from a rich set of both zoophilic and anthropophilic populations of S. guianense s.l. that were collected from eight sites that are broadly distributed across Brazil. Consistent with previous findings, our analysis supports the genetic isolation of Simulium litobranchium from S. guianense s.l. In contrast with previous findings, however, our results did not provide support for the divergence of the two species prior to the radiation of S. guianense s.l. In our analysis of the S. guianense s.l. ribosomal DNA sequence trace files we generated, we provide clear evidence of multiple within-specimen single nucleotide polymorphisms and indels suggesting that S. guianense s.l. ribosomal DNA is not a good target for conventional DNA barcoding. This is the first report of S. guianense s.l. within individual ribosomal DNA variation and thus the first evidence that the species is not subject to the normal effects of concerted evolution. Collectively, these data illustrate the need for diverse sampling in the development of robust molecular tools for vector identification and suggest that ribosomal DNA might be able to assist with resolving S. guianense s.l. species substructuring that C01 barcoding has hitherto failed to.

Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans

Wide arrays of repetitive DNA sequences form an important part of eukaryotic genomes. These repeats appear to evolve as coherent families, where repeats within a family are more similar to each other than to other orthologous representatives in related species. The continuous homogenization of repeats, through selective and non-selective processes, is termed concerted evolution. Ascertaining the level of variation between repeats is crucial to determining which evolutionary model best explains the homogenization observed for these sequences. Here, for the grasshopper Eyprepocnemis plorans, we present the analysis of intragenomic diversity for two repetitive DNA sequences (a satellite DNA (satDNA) and the 45S rDNA) resulting from the independent microdissection of several chromosomes. Our results show different homogenization patterns for these two kinds of paralogous DNA sequences, with a high between-chromosome structure for rDNA but no structure at all for the satDNA. This difference is puzzling, considering the adjacent localization of the two repetitive DNAs on paracentromeric regions in most chromosomes. The disparate homogenization patterns detected for these two repetitive DNA sequences suggest that several processes participate in the concerted evolution in E. plorans, and that these mechanisms might not work as genome-wide processes but rather as sequence-specific ones.

Amblyomma cajennense (Fabricius, 1787) (Acari: Ixodidae), the Cayenne tick: phylogeography and evidence for allopatric speciation

Background

Amblyomma cajennense F. is one of the best known and studied ticks in the New World because of its very wide distribution, its economical importance as pest of domestic ungulates, and its association with a variety of animal and human pathogens. Recent observations, however, have challenged the taxonomic status of this tick and indicated that intraspecific cryptic speciation might be occurring. In the present study, we investigate the evolutionary and demographic history of this tick and examine its genetic structure based on the analyses of three mitochondrial (12SrDNA, d-loop, and COII) and one nuclear (ITS2) genes. Because A. cajennense is characterized by a typical trans-Amazonian distribution, lineage divergence dating is also performed to establish whether genetic diversity can be linked to dated vicariant events which shaped the topology of the Neotropics.

Results

Total evidence analyses of the concatenated mtDNA and nuclear + mtDNA datasets resulted in well-resolved and fully congruent reconstructions of the relationships within A. cajennense. The phylogenetic analyses consistently found A. cajennense to be monophyletic and to be separated into six genetic units defined by mutually exclusive haplotype compositions and habitat associations. Also, genetic divergence values showed that these lineages are as distinct from each other as recognized separate species of the same genus. The six clades are deeply split and node dating indicates that they started diverging in the middle-late Miocene.

Conclusions

Behavioral differences and the results of laboratory cross-breeding experiments had already indicated that A. cajennense might be a complex of distinct taxonomic units. The combined and congruent mitochondrial and nuclear genetic evidence from this study reveals that A. cajennense is an assembly of six distinct species which have evolved separately from each other since at least 13.2 million years ago (Mya) in the earliest and 3.3 Mya in the latest lineages. The temporal and spatial diversification modes of the six lineages overlap the phylogeographical history of other organisms with similar extant trans-Amazonian distributions and are consistent with the present prevailing hypothesis that Neotropical diversity often finds its origins in the Miocene, after the Andean uplift changed the topology and consequently the climate and ecology of the Neotropics.

Intra-genomic variation in the ribosomal repeats of nematodes

Ribosomal loci represent a major tool for investigating environmental diversity and community structure via high-throughput marker gene studies of eukaryotes (e.g. 18S rRNA). Since the estimation of species' abundance is a major goal of environmental studies (by counting numbers of sequences), understanding the patterns of rRNA copy number across species will be critical for informing such high-throughput approaches. Such knowledge is critical, given that ribosomal RNA genes exist within multi-copy repeated arrays in a genome. Here we measured the repeat copy number for six nematode species by mapping the sequences from whole genome shotgun libraries against reference sequences for their rRNA repeat. This revealed a 6-fold variation in repeat copy number amongst taxa investigated, with levels of intragenomic variation ranging from 56 to 323 copies of the rRNA array. By applying the same approach to four C. elegans mutation accumulation lines propagated by repeated bottlenecking for an average of ~400 generations, we find on average a 2-fold increase in repeat copy number (rate of increase in rRNA estimated at 0.0285-0.3414 copies per generation), suggesting that rRNA repeat copy number is subject to selection. Within each Caenorhabditis species, the majority of intragenomic variation found across the rRNA repeat was observed within gene regions (18S, 28S, 5.8S), suggesting that such intragenomic variation is not a product of selection for rRNA coding function. We find that the dramatic variation in repeat copy number among these six nematode genomes would limit the use of rRNA in estimates of organismal abundance. In addition, the unique pattern of variation within a single genome was uncorrelated with patterns of divergence between species, reflecting a strong signature of natural selection for rRNA function. A better understanding of the factors that control or affect copy number in these arrays, as well as their rates and patterns of evolution, will be critical for informing estimates of global biodiversity.

Next generation sequencing analysis reveals a relationship between rDNA unit diversity and locus number in Nicotiana diploids

BACKGROUND

Tandemly arranged nuclear ribosomal DNA (rDNA), encoding 18S, 5.8S and 26S ribosomal RNA (rRNA), exhibit concerted evolution, a pattern thought to result from the homogenisation of rDNA arrays. However rDNA homogeneity at the single nucleotide polymorphism (SNP) level has not been detailed in organisms with more than a few hundred copies of the rDNA unit. Here we study rDNA complexity in species with arrays consisting of thousands of units.

METHODS

We examined homogeneity of genic (18S) and non-coding internally transcribed spacer (ITS1) regions of rDNA using Roche 454 and/or Illumina platforms in four angiosperm species, Nicotiana sylvestris, N. tomentosiformis, N. otophora and N. kawakamii. We compared the data with Southern blot hybridisation revealing the structure of intergenic spacer (IGS) sequences and with the number and distribution of rDNA loci.

RESULTS AND CONCLUSIONS

In all four species the intragenomic homogeneity of the 18S gene was high; a single ribotype makes up over 90% of the genes. However greater variation was observed in the ITS1 region, particularly in species with two or more rDNA loci, where >55% of rDNA units were a single ribotype, with the second most abundant variant accounted for >18% of units. IGS heterogeneity was high in all species. The increased number of ribotypes in ITS1 compared with 18S sequences may reflect rounds of incomplete homogenisation with strong selection for functional genic regions and relaxed selection on ITS1 variants. The relationship between the number of ITS1 ribotypes and the number of rDNA loci leads us to propose that rDNA evolution and complexity is influenced by locus number and/or amplification of orphaned rDNA units at new chromosomal locations.

Discovery of multiple IGS haplotypes within genotypes of Puccinia striiformis

In a search for specific molecular markers for population analysis of Puccinia striiformis f. sp. tritici, the ribosomal DNA (rDNA) intergenic spacer (IGS) 1 region (rDNA-IGS1, between the 28S and the 5S rDNA genes) was amplified, cloned, and sequenced. It was found to exhibit multiple bands and length polymorphism. Surprisingly, single isolates were found to possess between three to five different IGS1 haplotypes. Bands were cloned and sequenced, and two highly variable regions (α and β) were found between conserved regions, with repeat units interspersed in both types of regions. There were 14 different repeat units, and these were sometimes grouped further into four combinations of repeat units, with a few individual nucleotides (A or C) inserted between the repeats. Among three geographically dispersed isolates, the variable region α was divided into eight types, and the variable region β was divided into two types based on repeat units. Most of the 14 repeat units were shared by the variable and the conserved regions. Among the three isolates, there were a total of 12 IGS1 haplotypes, but some of these were shared between isolates such that there were only eight unique haplotypes. The occurrence of multiple haplotypes within single isolates may be useful for analyzing the population structure, tracking the origin of annual epidemics and providing insights into evolutionary biology of this pathogen.

Evaluation of the internal transcribed spacer 2 (ITS2) as a molecular marker for phylogenetic inference using sequence and secondary structure information in blow flies (Diptera: Calliphoridae)

The internal transcribed spacer 2 (ITS2) is a small non-coding region located inside the nuclear ribosomal DNA cluster. ITS2 sequence variability is thought to be appropriate to differentiate species and for phylogenetic reconstructions analyses, which can be further improved if structural information is considered. We evaluated the potential of ITS2 as a molecular marker for phylogenetic inference in Calliphoridae (Diptera: Brachycera) using a broad range of inference methods and different substitution models, accounting or not for structural information. Sequence analyses revealed a hierarchically organized pattern of sequence variation and a small level of nucleotide substitution saturation. Intragenomic variation due to small sequence repeats was found mainly in the most variable domain (IV), but it has no significant impact on the phylogenetic signal at the species level. Inferred secondary structures revealed that GC pairs are more frequently found flanking bulges and loops regions in more conserved domains, thus ensuring structure stability. In the phylogenetic analyses, the use of substitution models accounting for structural information significantly improves phylogenetic inference in both neighbour-joining and Bayesian analyses, although the former provides limited resolution for dealing with highly divergent sequences. For Bayesian analyses, a significant improvement in likelihood was observed when considering structure information, although with small changes in topology and overall support, probably reflecting better evolutionary rates estimates. Based on these findings, ITS2 is a suitable molecular marker for phylogenetic analyses in Calliphoridae, at both species and generic level.

Low evolutionary diversification in a widespread and abundant uncultured protist (MAST-4)

Recent culture-independent studies of marine planktonic protists have unveiled a large diversity at all phylogenetic scales and the existence of novel groups. MAST-4 represents one of these novel uncultured lineages, and it is composed of small (~2 μm) bacterivorous eukaryotes that are widely distributed in marine systems. MAST-4 accounts for a significant fraction of the marine heterotrophic flagellates at the global level, playing key roles in the marine ecological network. In this study, we investigated the diversity of MAST-4, aiming to assess its limits and structure. Using ribosomal DNA (rDNA) sequences obtained in this study (both pyrosequencing reads and clones with large rDNA operon coverage), complemented with GenBank sequences, we show that MAST-4 is composed of only five main clades, which are well supported by small subunit and large subunit phylogenies. The differences in the conserved regions of the internal transcribed spacers 1 and 2 (ITS1 and ITS2) secondary structures strongly suggest that these five clades are different biological species. Based on intraclade divergence, ITS secondary structures and comparisons of ITS1 and ITS2 trees, we did not find evidence of more than one species within clade A, whereas as many as three species might be present within other clades. Overall, the genetic divergence of MAST-4 was surprisingly low for an organism with a global population size estimated to be around 10(24), indicating a very low evolutionary diversification within the group.

The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays

Background

Ribosomal 5S genes are well known for the critical role they play in ribosome folding and functionality. These genes are thought to evolve in a concerted fashion, with high rates of homogenization of gene copies. However, the majority of previous analyses regarding the evolutionary process of rDNA repeats were conducted in invertebrates and plants. Studies have also been conducted on vertebrates, but these analyses were usually restricted to the 18S, 5.8S and 28S rRNA genes. The recent identification of divergent 5S rRNA gene paralogs in the genomes of elasmobranches and teleost fishes indicate that the eukaryotic 5S rRNA gene family has a more complex genomic organization than previously thought. The availability of new sequence data from lower vertebrates such as teleosts and elasmobranches enables an enhanced evolutionary characterization of 5S rDNA among vertebrates.

Results

We identified two variant classes of 5S rDNA sequences in the genomes of Potamotrygonidae stingrays, similar to the genomes of other vertebrates. One class of 5S rRNA genes was shared only by elasmobranches. A broad comparative survey among 100 vertebrate species suggests that the 5S rRNA gene variants in fishes originated from rounds of genome duplication. These variants were then maintained or eliminated by birth-and-death mechanisms, under intense purifying selection. Clustered multiple copies of 5S rDNA variants could have arisen due to unequal crossing over mechanisms. Simultaneously, the distinct genome clusters were independently homogenized, resulting in the maintenance of clusters of highly similar repeats through concerted evolution.

Conclusions

We believe that 5S rDNA molecular evolution in fish genomes is driven by a mixed mechanism that integrates birth-and-death and concerted evolution.

Incomplete homogenization of 18 S ribosomal DNA coding regions in Arabidopsis thaliana

Background

As a result of concerted evolution, coding regions of ribosomal DNA sequences are highly conserved within species and variation is generally thought to be limited to a few nucleotides. However, rDNA sequence variation has not been systematically examined in plant genomes, including that of the model plant Arabidopsis thaliana whose genome was the first to be sequenced.

Findings

Both genomic and transcribed 18 S sequences were sampled and revealed that most deviation from the consensus sequence was limited to single nucleotide substitutions except for a variant with a 270 bp deletion from position 456 to 725 in Arabidopsis numbering. The deletion maps to the functionally important and highly conserved 530 loop or helix18 in the structure of E. coli 16 S. The expression of the deletion variant is tightly controlled during developmental growth stages. Transcripts were not detectable in young seedlings but could be amplified from RNA extracts of mature leaves, stems, flowers and roots of Arabidopsis thaliana ecotype Columbia. We also show polymorphism for the deletion variant among four Arabidopsis ecotypes examined.

Conclusion

Despite a strong purifying selection that might be expected against functionally impaired rDNAs, the newly identified variant is maintained in the Arabidopsis genome. The expression of the variant and the polymorphism displayed by Arabidopsis ecotypes suggest a transition state in concerted evolution.

Polymorphism and recombination for rDNA in the putatively asexual microsporidian Nosema ceranae, a pathogen of honeybees

Nosema ceranae is currently one of the major pathogens of honeybees, related to the worldwide colony losses phenomenon. The genotyping of strains based on ribosomal DNA (rDNA) can be misleading if the repeated units are not identical. The analysis of cloned rDNA fragments containing the intergenic spacer (IGS) and part of the rDNA small-subunit (SSU) gene, from N. ceranae isolates from different European and Central Asia populations, revealed a high diversity of sequences. The variability involved single-nucleotide polymorphisms and insertion/deletions, resulting in 79 different haplotypes. Two sequences from the same isolate could be as different as any pair of sequences from different samples; in contrast, identical haplotypes were also found in very different geographical origins. Consequently, haplotypes cannot be organized in a consistent phylogenetic tree, clearly indicating that rDNA is not a reliable marker for the differentiation of N. ceranae strains. The results indicate that recombination between different sequences may produce new variants, which is quite surprising in microsporidia, usually considered to have an asexual mode of reproduction. The diversity of sequences and their geographical distribution indicate that haplotypes of different lineages may occasionally be present in a same cell and undergo homologue recombination, therefore suggesting a sexual haplo-diploid cycle.

Frequent silencing of rDNA loci on the univalent-forming genomes contrasts with their stable expression on the bivalent-forming genomes in polyploid dogroses (Rosa sect. Caninae)

The polyploid species in Rosa section Caninae (2n=21, 28 or 35) are characterized by an unusual reproductive system known as odd (or asymmetric) meiosis. Only two chromosome sets form bivalents in meiosis, whereas the remaining chromosomes are transmitted as univalents through the female germline. Evolution of ribosomal rRNA genes (rDNA) does not seem to be significantly affected by interlocus homogenization in dogroses. As a consequence, most species contain several rDNA families falling into two main clades (beta and gamma) thought to be differentially distributed between bivalent and univalent chromosomes, respectively. Here, we have investigated expression of rRNA gene families in five pentaploid species (R. canina, R. rubiginosa, R. dumalis, R. sherardii and R. caesia, 2n=35) and in one tetraploid (R. mollis, 2n=28). Using extensive sequencing of ITS clones and cleaved amplified polymorphism sequence (CAPS) analysis, we found that the beta-family was constitutively expressed in all species. However, there was large variation in the expression patterns of families constituting the gamma-clade. In addition, a single family can be active in one species, whereas silenced in another. The data show that the families on bivalent-forming chromosomes dominate rDNA expression in all dogrose species. We hypothesize that genes on bivalent genomes are stably expressed, whereas those on univalent genomes undergo variable levels of epigenetic silencing. Nonetheless, mosaic expression of univalent genomes could contribute to phenotypic variation between the species.

ITS2 secondary structure improves phylogeny estimation in a radiation of blue butterflies of the subgenus Agrodiaetus (Lepidoptera: Lycaenidae: Polyommatus )

BACKGROUND

Current molecular phylogenetic studies of Lepidoptera and most other arthropods are predominantly based on mitochondrial genes and a limited number of nuclear genes. The nuclear genes, however, generally do not provide sufficient information for young radiations. ITS2 , which has proven to be an excellent nuclear marker for similarly aged radiations in other organisms like fungi and plants, is only rarely used for phylogeny estimation in arthropods, although universal primers exist. This is partly due to difficulties in the alignment of ITS2 sequences in more distant taxa. The present study uses ITS2 secondary structure information to elucidate the phylogeny of a species-rich young radiation of arthropods, the butterfly subgenus Agrodiaetus. One aim is to evaluate the efficiency of ITS2 to resolve the phylogeny of the subgenus in comparison with COI , the most important mitochondrial marker in arthropods. Furthermore, we assess the use of compensatory base changes in ITS2 for the delimitation of species and discuss the prospects of ITS2 as a nuclear marker for barcoding studies.

RESULTS

In the butterfly family Lycaenidae, ITS2 secondary structure enabled us to successfully align sequences of different subtribes in Polyommatini and produce a Profile Neighbour Joining tree of this tribe, the resolution of which is comparable to phylogenetic trees obtained with COI+COII . The subgenus Agrodiaetus comprises 6 major clades which are in agreement with COI analyses. A dispersal-vicariance analysis (DIVA) traced the origin of most Agrodiaetus clades to separate biogeographical areas in the region encompassing Eastern Anatolia, Transcaucasia and Iran.

CONCLUSIONS

With the inclusion of secondary structure information, ITS2 appears to be a suitable nuclear marker to infer the phylogeny of young radiations, as well as more distantly related genera within a diverse arthropod family. Its phylogenetic signal is comparable to the mitochondrial marker COI . Compensatory base changes are very rare within Polyommatini and cannot be used for species delimitation. The implementation of secondary structure information into character-based phylogenetic methods is suggested to further improve the versatility of this marker in phylogenetic studies.

Invasions, DNA barcodes, and rapid biodiversity assessment using ants of Mauritius

BACKGROUND

Using an understudied taxon (Hymenoptera, Formicidae) found on a tropical island (Mauritius) where native flora and fauna have been threatened by 400 years of habitat modification and introduced species, we tested whether estimated incidences of diversity and complementarity were similar when measured by standard morphological alpha-taxonomy or phylogenetic diversity (PD) based on a standardized mitochondrial barcode and corroborating nuclear marker.

RESULTS

We found that costs related to site loss (considered loss of evolutionary history measured as loss of barcode PD) were not significantly different from predictions made either a) using standard morphology-based taxonomy, or b) measured using a nuclear marker. Integrating morphology and barcode results permitted us to identify a case of initially morphologically-cryptic variation as a new and endemic candidate species. However, barcode estimates of the relative importance of each site or network of sites were dramatically affected when the species in question was known to be indigenous or introduced.

CONCLUSION

This study goes beyond a mere demonstration of the rapid gains possible for diversity assessment using a standardized DNA barcode. Contextualization of these gains with ecological and natural history information is necessary to calibrate this wealth of standardized information. Without such an integrative approach, critical opportunities to advance knowledge will be missed.