Search

Multicopy sequences evolve adaptations for increasing their copy number within nuclei. The activities of multicopy sequences under constraints imposed by cellular and organismal selection result in a rich intranuclear ecology in germline cells. Mitochondrial and ribosomal DNA are managed as domestic herds subject to selective breeding by the genes of the single-copy genome. Transposable elements lead a peripatetic existence in which they must continually move to new sites to keep ahead of inactivating mutations at old sites and undergo exponential outbreaks when the production of new copies exceeds the rate of inactivation of old copies. Centromeres become populated by repeats that do little harm. Organisms with late sequestration of germ cells tend to evolve more "junk" in their genomes than organisms with early sequestration of germ cells. [ABSTRACT FROM AUTHOR]

Ribosomal DNA (rDNA) has a vital role in ribosome biogenesis as it contains the genes that encode ribosomal RNA (rRNA) separated by intergenic spacers (IGSs). The rRNA genes occur in hundreds to tens of thousands of copies per haploid genome in eukaryotes and are generally highly conserved with low variation within species. Due to the repetitive nature and large size of rDNA arrays, detecting intraindividual variation can be difficult. In this study, we use whole-genome sequences of 169 Daphnia pulex individuals from 10 natural populations to measure the copy number and sequence variation in rDNA. This revealed that variation in rDNA copy number between individuals spans an order of magnitude. We further observed a substantial level of sequence variation within individual genomes. As expected, single-nucleotide polymorphisms occurred in regions of lower functional constraint such as the IGS and expansion segments of the rRNA genes. The presence of strong linkage disequilibrium among variants facilitated identification of haplotypes within each population. Although there was evidence of recombination among haplotypes from different populations, it is insufficient to eliminate linkage disequilibrium within populations. Estimating copy number and haplotype diversity within individuals revealed that the level of intraindividual sequence variation is not strongly correlated with copy number. The observed patterns of variation highlight a complex evolutionary history of rDNA in D. pulex. Future research should explore the functional implications of rDNA copy number and sequence variation on organismal phenotypes. [ABSTRACT FROM AUTHOR]

Although both are salient features of genomes, at first glance ribosomal DNAs and transposable elements are genetic elements with not much in common: whereas ribosomal DNAs are mainly viewed as housekeeping genes that uphold all prime genome functions, transposable elements are generally portrayed as selfish and disruptive. These opposing characteristics are also mirrored in other attributes: organization in tandem (ribosomal DNAs) versus organization in a dispersed manner (transposable elements); evolution in a concerted manner (ribosomal DNAs) versus evolution by diversification (transposable elements); and activity that prolongs genomic stability (ribosomal DNAs) versus activity that shortens it (transposable elements). Re-visiting relevant instances in which ribosomal DNA–transposable element interactions have been reported, we note that both repeat types share at least four structural and functional hallmarks: (1) they are repetitive DNAs that shape genomes in evolutionary timescales, (2) they exchange structural motifs and can enter co-evolution processes, (3) they are tightly controlled genomic stress sensors playing key roles in senescence/aging, and (4) they share common epigenetic marks such as DNA methylation and histone modification. Here, we give an overview of the structural, functional, and evolutionary characteristics of both ribosomal DNAs and transposable elements, discuss their roles and interactions, and highlight trends and future directions as we move forward in understanding ribosomal DNA–transposable element associations. [ABSTRACT FROM AUTHOR]

The 5S rRNA genes are among the most conserved nucleotide sequences across all species. Similar to the 5S preservation we observe the occurrence of 5S-related nonautonomous retrotransposons, so-called Cassandras. Cassandras harbor highly conserved 5S rDNA-related sequences within their long terminal repeats, advantageously providing them with the 5S internal promoter. However, the dynamics of Cassandra retrotransposon evolution in the context of 5S rRNA gene sequence information and structural arrangement are still unclear, especially: (1) do we observe repeated or gradual domestication of the highly conserved 5S promoter by Cassandras and (2) do changes in 5S organization such as in the linked 35S-5S rDNA arrangements impact Cassandra evolution? Here, we show evidence for gradual co-evolution of Cassandra sequences with their corresponding 5S rDNAs. To follow the impact of 5S rDNA variability on Cassandra TEs, we investigate the Asteraceae family where highly variable 5S rDNAs, including 5S promoter shifts and both linked and separated 35S-5S rDNA arrangements have been reported. Cassandras within the Asteraceae mirror 5S rDNA promoter mutations of their host genome, likely as an adaptation to the host's specific 5S transcription factors and hence compensating for evolutionary changes in the 5S rDNA sequence. Changes in the 5S rDNA sequence and in Cassandras seem uncorrelated with linked/separated rDNA arrangements. We place all these observations into the context of angiosperm 5S rDNA-Cassandra evolution, discuss Cassandra's origin hypotheses (single or multiple) and Cassandra's possible impact on rDNA and plant genome organization, giving new insights into the interplay of ribosomal genes and transposable elements. [ABSTRACT FROM AUTHOR]

Introduction: The Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing nonself peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods. Methods: Here, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHCIIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages. Results: The analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and β2 chains of the MHC-II heterodimer. We detected no evidence for an overall (genewide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptidebinding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes. Discussion: Our results suggest that specific pairing of MHC-II α and β chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research. [ABSTRACT FROM AUTHOR]

IntroductionThe Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing non-self peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods.MethodsHere, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHC-IIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages.ResultsThe analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and β2 chains of the MHC-II heterodimer. We detected no evidence for an overall (gene-wide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptide-binding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes.DiscussionOur results suggest that specific pairing of MHC-II α and β chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research.

Ubiquitin is an extraordinarily highly conserved 76 amino acid protein encoded by three different types of gene, where the primary translation products are fusions either of ubiquitin with one of two ribosomal proteins (RPs) or of multiple ubiquitin monomers from head to tail. Here, we investigate the evolution of ubiquitin genes in mammalian malaria parasites (Plasmodium species). The ubiquitin encoded by the RPS27a fusion gene is highly divergent, as previously found in a variety of protists. However, we also find that two other forms of divergent ubiquitin sequence, each previously thought to be extremely rare, have arisen recently during the divergence of Plasmodium subgenera. On two occasions, in two distinct lineages, the ubiquitin encoded by the RPL40 fusion gene has rapidly diverged. In addition, in one of these lineages, the polyubiquitin genes have undergone a single codon insertion, previously considered a unique feature of Rhizaria. There has been disagreement whether the multiple ubiquitin coding repeats within a genome exhibit concerted evolution or undergo a birth-and-death process; the Plasmodium ubiquitin genes show clear signs of concerted evolution, including the spread of this codon insertion to multiple repeats within the polyubiquitin gene. [ABSTRACT FROM AUTHOR]

In most animals, mitochondrial DNA is strictly maternally inherited and non-recombining. One exception to this pattern is called doubly uniparental inheritance (DUI), a phenomenon involving the independent transmission of female and male mitochondrial genomes. DUI is known only from the molluskan class Bivalvia. The phylogenetic distribution of male-transmitted mitochondrial DNA (M mtDNA) in bivalves is consistent with several evolutionary scenarios, including multiple independent gains, losses, and varying degrees of recombination with female-transmitted mitochondrial DNA (F mtDNA). In this study, we use phylogenetic methods to test M mtDNA origination hypotheses and infer the prevalence of mitochondrial recombination in bivalves with DUI. Phylogenetic modeling using site concordance factors supported a single origin of M mtDNA in bivalves coupled with recombination acting over long evolutionary timescales. Ongoing mitochondrial recombination is present in Mytilida and Venerida, which results in a pattern of concerted evolution of F mtDNA and M mtDNA. Mitochondrial recombination could be favored to offset the deleterious effects of asexual inheritance and maintain mitonuclear compatibility across tissues. Cardiida and Unionida have gone without recent recombination, possibly due to an extension of the COX2 gene in male mitochondrial DNA. The loss of recombination could be connected to the role of M mtDNA in sex determination or sexual development. Our results support that recombination events may occur throughout the mitochondrial genomes of DUI species. Future investigations may reveal more complex patterns of inheritance of recombinants, which could explain the retention of signal for a single origination of M mtDNA in protein-coding genes. [ABSTRACT FROM AUTHOR]

Guanine and cytosine (GC) content is a fundamental component of genetic diversity and essential for phylogenetic analyses. However, the GC content of the ribosomal internal transcribed spacer 2 (ITS2) remains unknown, despite the fact that ITS2 is a widely used phylogenetic marker. Here, the ITS2 was high-throughput sequenced from 29 Corydalis species, and their GC contents were comparatively investigated in the context of ITS2's characteristic secondary structure and concerted evolution. Our results showed that the GC contents of ITS2 were 131% higher than those of their adjacent 5.8S regions, suggesting that ITS2 underwent GC-biased evolution. These GCs were distributed in a heterogeneous manner in the ITS2 secondary structure, with the paired regions being 130% larger than the unpaired regions, indicating that GC is chosen for thermodynamic stability. In addition, species with homogeneous ITS2 sequences were always GC-rich, supporting GC-biased gene conversion (gBGC), which occurred with ITS2's concerted evolution. The RNA substitution model inferred also showed a GC preference among base pair transformations, which again supports gBGC. Overall, structurally based GC investigation reveals that ITS2 evolves under structural stability and gBGC selection, significantly increasing its GC content. [ABSTRACT FROM AUTHOR]

Heat-shock (HS) assays to understand the connection between standing inversion variation and evolutionary response to climate change in Drosophila subobscura found that "warm-climate" inversion O3+4 exhibits non-HS levels of Hsp70 protein like those of "cold-climate" OST after HS induction. This was unexpected, as overexpression of Hsp70 can incur multiple fitness costs. To understand the genetic basis of this finding, we have determined the genomic sequence organization of the Hsp70 family in four different inversions, including OST , O3+4 , O3+4+8 and O3+4+16 , using as outgroups the remainder of the subobscura species subgroup, namely Drosophila madeirensis and Drosophila guanche. We found (i) in all the assayed lines, the Hsp70 family resides in cytological locus 94A and consists of only two genes, each with four HS elements (HSEs) and three GAGA sites on its promoter. Yet, in OST , the family is comparatively more compact; (ii) the two Hsp70 copies evolve in concert through gene conversion, except in D. guanche; (iii) within D. subobscura, the rate of concerted evolution is strongly structured by inversion, being higher in OST than in O3+4 ; and (iv) in D. guanche, the two copies accumulated multiple differences, including a newly evolved "gap-type" HSE2. The absence of concerted evolution in this species may be related to a long-gone-unnoticed observation that it lacks Hsp70 HS response, perhaps because it has evolved within a narrow thermal range in an oceanic island. Our results point to a previously unrealized link between inversions and concerted evolution, with potentially major implications for understanding genome evolution.

Heat-shock (HS) assays to understand the connection between standing inversion variation and evolutionary response to climate change in Drosophila subobscura found that "warm-climate" inversion O3+4 exhibits non-HS levels of Hsp70 protein like those of "cold-climate" OST after HS induction. This was unexpected, as overexpression of Hsp70 can incur multiple fitness costs. To understand the genetic basis of this finding, we have determined the genomic sequence organization of the Hsp70 family in four different inversions, including OST , O3+4 , O3+4+8 and O3+4+16 , using as outgroups the remainder of the subobscura species subgroup, namely Drosophila madeirensis and Drosophila guanche. We found (i) in all the assayed lines, the Hsp70 family resides in cytological locus 94A and consists of only two genes, each with four HS elements (HSEs) and three GAGA sites on its promoter. Yet, in OST , the family is comparatively more compact; (ii) the two Hsp70 copies evolve in concert through gene conversion, except in D. guanche; (iii) within D. subobscura, the rate of concerted evolution is strongly structured by inversion, being higher in OST than in O3+4 ; and (iv) in D. guanche, the two copies accumulated multiple differences, including a newly evolved "gap-type" HSE2. The absence of concerted evolution in this species may be related to a long-gone-unnoticed observation that it lacks Hsp70 HS response, perhaps because it has evolved within a narrow thermal range in an oceanic island. Our results point to a previously unrealized link between inversions and concerted evolution, with potentially major implications for understanding genome evolution.

Constitutive heterochromatin, consisting of repetitive sequences, diverges very rapidly; therefore, its nucleotide sequences and chromosomal distributions are often largely different, even between closely related species. The chromosome C-banding patterns of two Gerbillinae species, Meriones unguiculatus and Gerbillus perpallidus, vary greatly, even though they belong to the same subfamily. To understand the evolution of C-positive heterochromatin in these species, we isolated highly repetitive sequences, determined their nucleotide sequences, and characterized them using chromosomal and filter hybridization. We obtained a centromeric repeat (MUN-HaeIII) and a chromosome 13-specific repeat (MUN-EcoRI) from M. unguiculatus. We also isolated a centromeric/pericentromeric repeat (GPE-MBD) and an interspersed-type repeat that was predominantly amplified in the X and Y chromosomes (GPE-EcoRI) from G. perpallidus. GPE-MBD was found to contain a 17-bp motif that is essential for binding to the centromere-associated protein CENP-B. This indicates that it may play a role in the formation of a specified structure and/or function of centromeres. The nucleotide sequences of the three sequence families, except GPE-EcoRI, were conserved only in Gerbillinae. GPE-EcoRI was derived from the long interspersed nuclear elements 1 retrotransposon and showed sequence homology throughout Muridae and Cricetidae species, indicating that the repeat sequence occurred at least in the common ancestor of Muridae and Cricetidae. Due to a lack of assembly data of highly repetitive sequences constituting heterochromatin in whole-genome sequences of vertebrate species published to date, the knowledge obtained in this study provides useful information for a deep understanding of the evolution of repetitive sequences in not only rodents but also in mammals. [ABSTRACT FROM AUTHOR]

Background and Aims Polyploidy is an important process that often generates genomic diversity within lineages, but it can also cause changes that result in loss of genomic material. Island lineages, while often polyploid, typically show chromosomal stasis but have not been investigated in detail regarding smaller-scale gene loss. Our aim was to investigate post-polyploidization genome dynamics in a chromosomally stable lineage of Malvaceae endemic to New Zealand. Methods We determined chromosome numbers and used fluorescence in situ hybridization to localize 18S and 5S rDNA. Gene sequencing of 18S rDNA, the internal transcribed spacers (ITS) with intervening 5.8S rDNA, and a low-copy nuclear gene, GBSSI-1 , was undertaken to determine if gene loss occurred in the New Zealand lineage following polyploidy. Key Results The chromosome number for all species investigated was 2 n  = 42, with the first published report for the monotypic Australian genus Asterotrichion. The five species investigated all had two 5S rDNA signals localized interstitially on the long arm of one of the largest chromosome pairs. All species, except Plagianthus regius , had two 18S rDNA signals localized proximally on the short arm of one of the smallest chromosome pairs. Plagianthus regius had two additional 18S rDNA signals on a separate chromosome, giving a total of four. Sequencing of nuclear ribosomal 18S rDNA and the ITS cistron indicated loss of historical ribosomal repeats. Phylogenetic analysis of a low-copy nuclear gene, GBSSI-1 , indicated that some lineages maintained three copies of the locus, while others have lost one or two copies. Conclusions Although island endemic lineages show chromosomal stasis, with no additional changes in chromosome number, they may undergo smaller-scale processes of gene loss and concerted evolution ultimately leading to further genome restructuring and downsizing. [ABSTRACT FROM AUTHOR]

Orb-web weaving spiders use a variety of silk types for particular tasks, and each silk type is composed of at least two spider silk proteins (spidroins). In the early stage of divergence, however, the molecular evolutionary processes act on spidroin variants are still unclear because of a lack of knowledge for full-length paralogous and orthologous gene sequences among closely related species. Here, we present two complete gene sequences encoding the tubuliform spidroin TuSp1 variants (TuSp1-v2 and TuSp1-v3) from orb-weaving spider Neoscona theisi. Both N. theisi TuSp1-v2 and TuSp1-v3 genes contain a single enormous exon (14,139 bp for TuSp1-v2 and 13,152 bp for TuSp1-v3) and dozens of tandemly arrayed repeats (25 repeats for TuSp1-v2 and 23 repeats for TuSp1-v3) with extreme intragenic homogenization. The pattern of expression for these two spidroins revealed that the level of TuSp1-v3 mRNA is ~3-fold higher than that of TuSp1-v2 in tubuliform gland. Phylogenetic analyses of spidroins not only show the occurrence of a gene duplication event for TuSp1-v2 and TuSp1-v3 in the common ancestor of the Neoscona and Araneus lineage but reinforce the role of concerted evolution for the extreme homogenization of TuSp1 repeats. [ABSTRACT FROM AUTHOR]

The greater long-tailed hamster (Tscherskia triton, Cricetinae) has a unique karyotype (2n = 28), containing 11 pairs of acrocentric chromosomes with large C-band-positive centromeric heterochromatin blocks. To understand the origin and evolutionary process of heterochromatin in this species, we isolated novel families of chromosome site-specific highly repetitive DNA sequences from TaqI-digested genomic DNA and then characterized them by chromosome in situ and filter hybridization. The TaqI-families of repetitive sequences were classified into 2 types by their genome organization and chromosomal distribution: the 110-bp repeated sequence organized in large tandem arrays (as satellite DNA), localized to centromeric C-positive heterochromatin of acrocentric autosomes (chromosomes 1–11) and submetacentric X chromosome, and the 405-bp repeated sequence that was composed of 30–32-bp internal repeats, distributed in the pericentromeric region on the short arms of X and Y chromosomes. The repetitive sequences did not cross-hybridize with genomic DNA of any genera of Cricetinae (Mesocricetus, Cricetulus, and Phodopus). These results suggest that the 110-bp and 405-bp repeats rapidly diverged in the lineage of T. triton, evolving in a concerted manner among autosomes and X chromosome and within X and Y chromosomes, respectively. The 110-bp centromeric repeat contained a 17-bp motif in which 9 bases are essential for binding with the centromere-associated protein CENP-B, suggesting the possibility that the 110-bp major satellite DNA carrying the 17-bp motif may have a role in the formation of specified structure and/or function of centromeres in T. triton. [ABSTRACT FROM AUTHOR]

Intragenomic variation is the molecular variation within the genome among repetitive DNA. As a multigene family, nuclear ribosomal DNA (rDNA) has been widely used in fungal taxonomy for their ease in amplification and suitable variability to attain various levels of taxonomic resolution. At the intraspecific level, rDNA is believed to be under concerted evolution and the internal transcribed spacers (ITS) region is actually accepted as a universal barcoding marker for fungi. However, documentation of intragenomic variation of rDNA indicated that it can be problematic in species delimitation and identification. Fungal taxonomic studies have not generally taken into account the intragenomic variation of rDNA in a systematic manner. In this review, our objective is to address the definition, the origin and the mechanisms for maintenance of intragenomic variation, as well as its implication in the domain of fungal molecular taxonomy, particularly for species delimitation, identification and DNA barcoding. With advanced sequencing technologies (second and third generations), we also addressed how these technologies can be used to study the intragenomic variation of rDNA and also how the intragenomic variation will impact on DNA barcoding via high-throughput sequencing. • nrDNA is the most used genomic region for fungal identification and barcoding. • But intragenomic variation can be problematics for identification and barcoding. • In silico survey of intragenomic variation of ITS was undertaken. • High-throughput sequencing (HTS) is effective for studying the intragenomic variation. • HTS-based identification of fungi must take into account the intragenomic variation. [ABSTRACT FROM AUTHOR]

The paper nautilus or greater argonaut, Argonauta argo , is a species of octopods which is characterized by its pelagic lifestyle and by the presence of a protective spiral-shaped shell-like eggcase in females. To reveal the genomic background of how the species adapted to the pelagic lifestyle and acquired its shell-like eggcase, we sequenced the draft genome of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters, and some gene clusters that could probably be related to the pelagic lifestyle, such as reflectin , tyrosinase , and opsin. The gene models also revealed several homologous genes related to calcified shell formation in Conchiferan mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the shell-less octopus, as well as Nautilus , which has a true outer shell. Therefore, the draft genome sequence of Arg. argo presented here has helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes to form an important, converging extended phenotypic structure such as the shell and the shell-like eggcase. Additionally, it allows us to explore the evolution of from benthic to pelagic lifestyles in cephalopods and octopods. [ABSTRACT FROM AUTHOR]

Homologous centromeres compete for segregation to the secondary oocyte nucleus at female meiosis I. Centromeric repeats also compete with each other to populate centromeres in mitotic cells of the germline and have become adapted to use the recombinational machinery present at centromeres to promote their own propagation. Repeats are not needed at centromeres, rather centromeres appear to be hospitable habitats for the colonization and proliferation of repeats. This is probably an indirect consequence of two distinctive features of centromeric DNA. Centromeres are subject to breakage by the mechanical forces exerted by microtubules and meiotic crossing‐over is suppressed. Centromeric proteins acting in trans are under selection to mitigate the costs of centromeric repeats acting in cis. Collateral costs of mitotic competition at centromeres may help to explain the high rates of aneuploidy observed in early human embryos. [ABSTRACT FROM AUTHOR]

Apomictic genera have a complex evolutionary history, including reticulate hybridization, polyploidization and variation in reproduction modes. Restrictions of functional meiosis in polyploid apomictic taxa considerably hamper gene conversion, leaving footprints of past hybridization events in their genomes, which are masked by the homogenization of tandemly organized arrays of nuclear ribosomal DNA in sexual species. Dandelions (Taraxacum) have a highly complex reticulate evolutionary history. Detailed knowledge remains hazy; earlier investigations of rDNA sequences have uncovered the complex pattern derived from evolution but without direct evidence for any particular processes. We investigated the position and number of 45S and 5S rDNA loci in 38 Taraxacum taxa (covering different reproduction modes, geographical regions and putative phylogenetic groups) using fluorescent in situ hybridization (FISH) and measured genome size and GC content. The ITS1-5.8S-ITS2 regions of four sexual and five apomictic taxa were sequenced to investigate inter- and intra-individual variation. Most species considerably differ in the chromosome positions of loci and karyotype patterns, but mostly share the same number of studied loci (45S:5S) in a 1:2 ratio per haploid genome (x = 8), with six exceptions (up to a 4:2 ratio). Genome size (2C) varies sixfold and with GC content partly distinguishes major evolutionary groups. Sexual taxa show limited variation in sequenced regions (with two to eight variants), but apomictic taxa vary significantly (with 20–36 variants). Extensive reticulate evolution in Taraxacum and subsequent phenomena such as genome repatterning and non-effective concerted evolution are probably the cause of the dynamic nature of Taraxacum karyotypes and the large variation in genome size and rDNA sequences. [ABSTRACT FROM AUTHOR]

The complete mitochondrial genome was sequenced in a mangrove rivulus Kryptolebias marmoratus from western Florida using next-generation sequencing. The 17 329 bp-long genome was identical in length and 99.8% similar to a previously published genome of this species from a specimen of unknown geographic origin. Gene arrangement in K. marmoratus is similar to other cyprinodontiform fishes, except for the presence of a second copy of the control region inserted upstream of the nad1 gene.

Abstract Background Vertebrate brain structure is characterised not only by relative consistency in scaling between components, but also by many examples of divergence from these general trends.. Alternative hypotheses explain these patterns by emphasising either ‘external’ processes, such as coordinated or divergent selection, or ‘internal’ processes, like developmental coupling among brain regions. Although these hypotheses are not mutually exclusive, there is little agreement over their relative importance across time or how that importance may vary across evolutionary contexts. Results We introduce an agent-based model to simulate brain evolution in a ‘bare-bones’ system and examine dependencies between variables shaping brain evolution. We show that ‘concerted’ patterns of brain evolution do not, in themselves, provide evidence for developmental coupling, despite these terms often being treated as synonymous in the literature. Instead, concerted evolution can reflect either functional or developmental integration. Our model further allows us to clarify conditions under which such developmental coupling, or uncoupling, is potentially adaptive, revealing support for the maintenance of both mechanisms in neural evolution. Critically, we illustrate how the probability of deviation from concerted evolution depends on the cost/benefit ratio of neural tissue, which increases when overall brain size is itself under constraint. Conclusions We conclude that both developmentally coupled and uncoupled brain architectures can provide adaptive mechanisms, depending on the distribution of selection across brain structures, life history and costs of neural tissue. However, when constraints also act on overall brain size, heterogeneity in selection across brain structures will favour region specific, or mosaic, evolution. Regardless, the respective advantages of developmentally coupled and uncoupled brain architectures mean that both may persist in fluctuating environments. This implies that developmental coupling is unlikely to be a persistent constraint, but could evolve as an adaptive outcome to selection to maintain functional integration.

Guanine and cytosine (GC) content is a fundamental component of genetic diversity and essential for phylogenetic analyses. However, the GC content of the ribosomal internal transcribed spacer 2 (ITS2) remains unknown, despite the fact that ITS2 is a widely used phylogenetic marker. Here, the ITS2 was high-throughput sequenced from 29 Corydalis species, and their GC contents were comparatively investigated in the context of ITS2’s characteristic secondary structure and concerted evolution. Our results showed that the GC contents of ITS2 were 131% higher than those of their adjacent 5.8S regions, suggesting that ITS2 underwent GC-biased evolution. These GCs were distributed in a heterogeneous manner in the ITS2 secondary structure, with the paired regions being 130% larger than the unpaired regions, indicating that GC is chosen for thermodynamic stability. In addition, species with homogeneous ITS2 sequences were always GC-rich, supporting GC-biased gene conversion (gBGC), which occurred with ITS2’s concerted evolution. The RNA substitution model inferred also showed a GC preference among base pair transformations, which again supports gBGC. Overall, structurally based GC investigation reveals that ITS2 evolves under structural stability and gBGC selection, significantly increasing its GC content.

Summary: We present here the complete mtDNA genome (mitogenome) of Heterorhabditis bacteriophora , an important biological control agent of soil-dwelling insect pests in agriculture and horticulture. This is the first description of a mitogenome for a member of the family Heterorhabditidae. The genome contains the typical chromadorean complement of 12 protein-coding genes, 22 tRNA genes and two rRNA genes. All genes are transcribed in the same direction and have a nucleotide composition high in A and T. For the entire genome, the nucleotide contents are 47.02% (T), 28.81% (A), 16.10% (G), 8.08% (C) and 75.83% (AT). Heterorhabditis bacteriophora has a unique, idiosyncratic gene arrangement. It differs from that of Caenorhabditis elegans in having a block of seven genes: trnQ-trnF-cytb-trnL1-cox3-trnT-nad4 translocated to a position between nad3 and nad5 , as well as having a change in the position of the four tRNA block gene cluster, trnC-trnM-trnD-trnG , where trnC and trnM have switched places and trnD and trnG have translocated between nad4 and nad5 genes. The H. bacteriophora mitogenome is 18 128 bp long, and thus is ca 4 kb larger than the mitogenomes of most chromadoreans. This relatively large genome is due to the presence of five non-coding regions (NCR): NCR1 (114 bp), NCR2 (159 bp), NCR3 (498 bp), NCR4 (1917 bp) and NCR5 (2154 bp), which make up 26.7% of the genome. The NCR5 had the highest A + T content of 83.47% indicating that this region is the likely AT-rich control region. The complete 498 bp NCR3 sequence is duplicated in NCR4 and in NCR5 (the putative AT-rich control region). Such an organisation has not been reported previously in nematode mtDNA. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Neotropical fishes have highly diversified karyotypic and genomic characteristics and present many diverse sex chromosome systems, with various degrees of sex chromosome differentiation. Knowledge on their sex-specific composition and evolution, however, is still limited. Satellite DNAs (satDNAs) are tandemly repeated sequences with pervasive genomic distribution and distinctive evolutionary pathways, and investigating satDNA content might shed light into how genome architecture is organized in fishes and in their sex chromosomes. The present study investigated the satellitome of Megaleporinus elongatus , a freshwater fish with a proposed Z1Z1Z2Z2/Z1W1Z2W2 multiple sex chromosome system that encompasses a highly heterochromatic and differentiated W1 chromosome. The species satellitome comprises of 140 different satDNA families, including previously isolated sequences and new families found in this study. This diversity is remarkable considering the relatively low proportion that satDNAs generally account for the M. elongatus genome (around only 5%). Differences between the sexes in regards of satDNA content were also evidenced, as these sequences are 14% more abundant in the female genome. The occurrence of sex-biased signatures of satDNA evolution in the species is tightly linked to satellite enrichment associated with W1 in females. Although both sexes share practically all satDNAs, the overall massive amplification of only a few of them accompanied the W1 differentiation. We also investigated the expansion and diversification of the two most abundant satDNAs of M. elongatus, MelSat01-36 and MelSat02-26, both highly amplified sequences in W1 and, in MelSat02-26's case, also harbored by Z2 and W2 chromosomes. We compared their occurrences in M. elongatus and the sister species M. macrocephalus (with a standard ZW sex chromosome system) and concluded that both satDNAs have led to the formation of highly amplified arrays in both species; however, they formed species-specific organization on female-restricted sex chromosomes. Our results show how satDNA composition is highly diversified in M. elongatus , in which their accumulation is significantly contributing to W1 differentiation and not satDNA diversity per se. Also, the evolutionary behavior of these repeats may be associated with genome plasticity and satDNA variability between the sexes and between closely related species, influencing how seemingly homeologous heteromorphic sex chromosomes undergo independent satDNA evolution. [ABSTRACT FROM AUTHOR]

The rearrangement of 37 genes with one control region, firstly identified in Gallus gallus mitogenome, is believed to be ancestral for all Aves. However, mitogenomic sequences obtained in recent years revealed that many avian mitogenomes contain duplicated regions that were omitted in previous genomic versions. Their evolution and mechanism of duplication are still poorly understood. The order of Accipitriformes is especially interesting in this context because its representatives contain a duplicated control region in various stages of degeneration. Therefore, we applied an appropriate PCR strategy to look for duplications within the mitogenomes of the early diverged species Sagittarius serpentarius and Cathartiformes, which is a sister order to Accipitriformes. The analyses revealed the same duplicated gene order in all examined taxa and the common ancestor of these groups. The duplicated regions were subjected to gradual degeneration and homogenization during concerted evolution. The latter process occurred recently in the species of Cathartiformes as well as in the early diverged lineages of Accipitriformes, that is, Sagittarius serpentarius and Pandion haliaetus. However, in other lineages, that is, Pernis ptilorhynchus , as well as representatives of Aegypiinae, Aquilinae, and five related subfamilies of Accipitriformes (Accipitrinae, Circinae, Buteoninae, Haliaeetinae, and Milvinae), the duplications were evolving independently for at least 14–47 Myr. Different portions of control regions in Cathartiformes showed conflicting phylogenetic signals indicating that some sections of these regions were homogenized at a frequency higher than the rate of speciation, whereas others have still evolved separately. [ABSTRACT FROM AUTHOR]

Neotropical fishes have highly diversified karyotypic and genomic characteristics and present many diverse sex chromosome systems, with various degrees of sex chromosome differentiation. Knowledge on their sex-specific composition and evolution, however, is still limited. Satellite DNAs (satDNAs) are tandemly repeated sequences with pervasive genomic distribution and distinctive evolutionary pathways, and investigating satDNA content might shed light into how genome architecture is organized in fishes and in their sex chromosomes. The present study investigated the satellitome of Megaleporinus elongatus, a freshwater fish with a proposed Z1Z1Z2Z2/Z1W1Z2W2 multiple sex chromosome system that encompasses a highly heterochromatic and differentiated W1 chromosome. The species satellitome comprises of 140 different satDNA families, including previously isolated sequences and new families found in this study. This diversity is remarkable considering the relatively low proportion that satDNAs generally account for the M. elongatus genome (around only 5%). Differences between the sexes in regards of satDNA content were also evidenced, as these sequences are 14% more abundant in the female genome. The occurrence of sex-biased signatures of satDNA evolution in the species is tightly linked to satellite enrichment associated with W1 in females. Although both sexes share practically all satDNAs, the overall massive amplification of only a few of them accompanied the W1 differentiation. We also investigated the expansion and diversification of the two most abundant satDNAs of M. elongatus, MelSat01-36 and MelSat02-26, both highly amplified sequences in W1 and, in MelSat02-26’s case, also harbored by Z2 and W2 chromosomes. We compared their occurrences in M. elongatus and the sister species M. macrocephalus (with a standard ZW sex chromosome system) and concluded that both satDNAs have led to the formation of highly amplified arrays in both species; however, they formed species-specific organization on female-restricted sex chromosomes. Our results show how satDNA composition is highly diversified in M. elongatus, in which their accumulation is significantly contributing to W1 differentiation and not satDNA diversity per se. Also, the evolutionary behavior of these repeats may be associated with genome plasticity and satDNA variability between the sexes and between closely related species, influencing how seemingly homeologous heteromorphic sex chromosomes undergo independent satDNA evolution.

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat‐enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long‐read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth‐and‐death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra‐ and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus. [ABSTRACT FROM AUTHOR]

Copyright of Chinese Journal of Applied Entomology is the property of Chinese Journal of Applied Entomology, Editorial Department and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

BackgroundSpider silks are spectacular examples of phenotypic diversity arising from adaptive molecular evolution. An individual spider can produce an array of specialized silks, with the majority of constituent silk proteins encoded by members of the spidroin gene family. Spidroins are dominated by tandem repeats flanked by short, non-repetitive N- and C-terminal coding regions. The remarkable mechanical properties of spider silks have been largely attributed to the repeat sequences. However, the molecular evolutionary processes acting on spidroin terminal and repetitive regions remain unclear due to a paucity of complete gene sequences and sampling of genetic variation among individuals. To better understand spider silk evolution, we characterize a complete aciniform spidroin gene from an Argiope orb-weaving spider and survey aciniform gene fragments from congeneric individuals.ResultsWe present the complete aciniform spidroin (AcSp1) gene from the silver garden spider Argiope argentata (Aar_AcSp1), and document multiple AcSp1 loci in individual genomes of A. argentata and the congeneric A. trifasciata and A. aurantia. We find that Aar_AcSp1 repeats have >98% pairwise nucleotide identity. By comparing AcSp1 repeat amino acid sequences between Argiope species and with other genera, we identify regions of conservation over vast amounts of evolutionary time. Through a PCR survey of individual A. argentata, A. trifasciata, and A. aurantia genomes, we ascertain that AcSp1 repeats show limited variation between species whereas terminal regions are more divergent. We also find that average dN/dS across codons in the N-terminal, repetitive, and C-terminal encoding regions indicate purifying selection that is strongest in the N-terminal region.ConclusionsUsing the complete A. argentata AcSp1 gene and spidroin genetic variation between individuals, this study clarifies some of the molecular evolutionary processes underlying the spectacular mechanical attributes of aciniform silk. It is likely that intragenic concerted evolution and functional constraints on A. argentata AcSp1 repeats result in extreme repeat homogeneity. The maintenance of multiple AcSp1 encoding loci in Argiope genomes supports the hypothesis that Argiope spiders require rapid and efficient protein production to support their prolific use of aciniform silk for prey-wrapping and web-decorating. In addition, multiple gene copies may represent the early stages of spidroin diversification.

BackgroundSpider silks are spectacular examples of phenotypic diversity arising from adaptive molecular evolution. An individual spider can produce an array of specialized silks, with the majority of constituent silk proteins encoded by members of the spidroin gene family. Spidroins are dominated by tandem repeats flanked by short, non-repetitive N- and C-terminal coding regions. The remarkable mechanical properties of spider silks have been largely attributed to the repeat sequences. However, the molecular evolutionary processes acting on spidroin terminal and repetitive regions remain unclear due to a paucity of complete gene sequences and sampling of genetic variation among individuals. To better understand spider silk evolution, we characterize a complete aciniform spidroin gene from an Argiope orb-weaving spider and survey aciniform gene fragments from congeneric individuals.ResultsWe present the complete aciniform spidroin (AcSp1) gene from the silver garden spider Argiope argentata (Aar_AcSp1), and document multiple AcSp1 loci in individual genomes of A. argentata and the congeneric A. trifasciata and A. aurantia. We find that Aar_AcSp1 repeats have >98% pairwise nucleotide identity. By comparing AcSp1 repeat amino acid sequences between Argiope species and with other genera, we identify regions of conservation over vast amounts of evolutionary time. Through a PCR survey of individual A. argentata, A. trifasciata, and A. aurantia genomes, we ascertain that AcSp1 repeats show limited variation between species whereas terminal regions are more divergent. We also find that average dN/dS across codons in the N-terminal, repetitive, and C-terminal encoding regions indicate purifying selection that is strongest in the N-terminal region.ConclusionsUsing the complete A. argentata AcSp1 gene and spidroin genetic variation between individuals, this study clarifies some of the molecular evolutionary processes underlying the spectacular mechanical attributes of aciniform silk. It is likely that intragenic concerted evolution and functional constraints on A. argentata AcSp1 repeats result in extreme repeat homogeneity. The maintenance of multiple AcSp1 encoding loci in Argiope genomes supports the hypothesis that Argiope spiders require rapid and efficient protein production to support their prolific use of aciniform silk for prey-wrapping and web-decorating. In addition, multiple gene copies may represent the early stages of spidroin diversification.

Abstract Sea anemones have a wide array of toxic compounds (peptide toxins found in their venom) which have potential uses as therapeutics. To date, the majority of studies characterizing toxins in sea anemones have been restricted to species from the superfamily, Actinioidea. No highly complete draft genomes are currently available for this superfamily, however, highlighting our limited understanding of the genes encoding toxins in this important group. Here we have sequenced, assembled, and annotated a draft genome for Actinia tenebrosa. The genome is estimated to be approximately 255 megabases, with 31,556 protein‐coding genes. Quality metrics revealed that this draft genome matches the quality and completeness of other model cnidarian genomes, including Nematostella, Hydra, and Acropora. Phylogenomic analyses revealed strong conservation of the Cnidaria and Hexacorallia core‐gene set. However, we found that lineage‐specific gene families have undergone significant expansion events compared with shared gene families. Enrichment analysis performed for both gene ontologies, and protein domains revealed that genes encoding toxins contribute to a significant proportion of the lineage‐specific genes and gene families. The results make clear that the draft genome of A. tenebrosa will provide insight into the evolution of toxins and lineage‐specific genes, and provide an important resource for the discovery of novel biological compounds.

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

Molecular evolution of ribosomal DNA can be highly dynamic. Hundreds to thousands of copies in the genome are subject to concerted evolution, which homogenizes sequence variants to different degrees. If well homogenized, sequences are suitable for phylogeny reconstruction; if not, sequence polymorphism has to be handled appropriately. Here we investigate non-coding rDNA sequences (ITS/ETS, 5S-NTS) along with the chromosomal organization of their respective loci (45S and 5S rDNA) in diploids of the Hieraciinae. The subtribe consists of genera Hieracium, Pilosella, Andryala, and Hispidella and has a complex evolutionary history characterized by ancient intergeneric hybridization, allele sharing among species, and incomplete lineage sorting. Direct or cloned Sanger sequences and phased alleles derived from Illumina genome sequencing were subjected to phylogenetic analyses. Patterns of homogenization and tree topologies based on the three regions were compared. In contrast to most other plant groups, 5S-NTS sequences were generally better homogenized than ITS and ETS sequences. A novel case of ancient intergeneric hybridization between Hispidella and Hieracium was inferred, and some further incongruences between the trees were found, suggesting independent evolution of these regions. In some species, homogenization of ITS/ETS and 5S-NTS sequences proceeded in different directions although the 5S rDNA locus always occurred on the same chromosome with one 45S rDNA locus. The ancestral rDNA organization in the Hieraciinae comprised 4 loci of 45S rDNA in terminal positions and 2 loci of 5S rDNA in interstitial positions per diploid genome. In Hieracium, some deviations from this general pattern were found (3, 6, or 7 loci of 45S rDNA; three loci of 5S rDNA). Some of these deviations concerned intraspecific variation, and most of them occurred at the tips of the tree or independently in different lineages. This indicates that the organization of rDNA loci is more dynamic than the evolution of sequences contained in them and that locus number is therefore largely unsuitable to inform about species relationships in Hieracium. No consistent differences in the degree of sequence homogenization and the number of 45S rDNA loci were found, suggesting interlocus concerted evolution.

Background: Vertebrate brain structure is characterised not only by relative consistency in scaling between components, but also by many examples of divergence from these general trends.. Alternative hypotheses explain these patterns by emphasising either 'external' processes, such as coordinated or divergent selection, or 'internal' processes, like developmental coupling among brain regions. Although these hypotheses are not mutually exclusive, there is little agreement over their relative importance across time or how that importance may vary across evolutionary contexts. Results: We introduce an agent-based model to simulate brain evolution in a 'bare-bones' system and examine dependencies between variables shaping brain evolution. We show that 'concerted' patterns of brain evolution do not, in themselves, provide evidence for developmental coupling, despite these terms often being treated as synonymous in the literature. Instead, concerted evolution can reflect either functional or developmental integration. Our model further allows us to clarify conditions under which such developmental coupling, or uncoupling, is potentially adaptive, revealing support for the maintenance of both mechanisms in neural evolution. Critically, we illustrate how the probability of deviation from concerted evolution depends on the cost/benefit ratio of neural tissue, which increases when overall brain size is itself under constraint. Conclusions: We conclude that both developmentally coupled and uncoupled brain architectures can provide adaptive mechanisms, depending on the distribution of selection across brain structures, life history and costs of neural tissue. However, when constraints also act on overall brain size, heterogeneity in selection across brain structures will favour region specific, or mosaic, evolution. Regardless, the respective advantages of developmentally coupled and uncoupled brain architectures mean that both may persist in fluctuating environments. This implies that developmental coupling is unlikely to be a persistent constraint, but could evolve as an adaptive outcome to selection to maintain functional integration. [ABSTRACT FROM AUTHOR]

Molecular evolution of ribosomal DNA can be highly dynamic. Hundreds to thousands of copies in the genome are subject to concerted evolution, which homogenizes sequence variants to different degrees. If well homogenized, sequences are suitable for phylogeny reconstruction; if not, sequence polymorphism has to be handled appropriately. Here we investigate non-coding rDNA sequences (ITS/ETS, 5S-NTS) along with the chromosomal organization of their respective loci (45S and 5S rDNA) in diploids of the Hieraciinae. The subtribe consists of genera Hieracium , Pilosella , Andryala , and Hispidella and has a complex evolutionary history characterized by ancient intergeneric hybridization, allele sharing among species, and incomplete lineage sorting. Direct or cloned Sanger sequences and phased alleles derived from Illumina genome sequencing were subjected to phylogenetic analyses. Patterns of homogenization and tree topologies based on the three regions were compared. In contrast to most other plant groups, 5S-NTS sequences were generally better homogenized than ITS and ETS sequences. A novel case of ancient intergeneric hybridization between Hispidella and Hieracium was inferred, and some further incongruences between the trees were found, suggesting independent evolution of these regions. In some species, homogenization of ITS/ETS and 5S-NTS sequences proceeded in different directions although the 5S rDNA locus always occurred on the same chromosome with one 45S rDNA locus. The ancestral rDNA organization in the Hieraciinae comprised 4 loci of 45S rDNA in terminal positions and 2 loci of 5S rDNA in interstitial positions per diploid genome. In Hieracium , some deviations from this general pattern were found (3, 6, or 7 loci of 45S rDNA; three loci of 5S rDNA). Some of these deviations concerned intraspecific variation, and most of them occurred at the tips of the tree or independently in different lineages. This indicates that the organization of rDNA loci is more dynamic than the evolution of sequences contained in them and that locus number is therefore largely unsuitable to inform about species relationships in Hieracium. No consistent differences in the degree of sequence homogenization and the number of 45S rDNA loci were found, suggesting interlocus concerted evolution. [ABSTRACT FROM AUTHOR]

Spider silk fibers have impressive mechanical properties and are primarily composed of highly repetitive structural proteins (termed spidroins) encoded by a single gene family. Most characterized spidroin genes are incompletely known because of their extreme size (typically >9 kb) and repetitiveness, limiting understanding of the evolutionary processes that gave rise to their unusual gene architectures. The only complete spidroin genes characterized thus far form the dragline in the Western black widow, Latrodectus hesperus. Here, we describe the first complete gene sequence encoding the aciniform spidroin AcSp1, the primary component of spider prey-wrapping fibers. L. hesperus AcSp1 contains a single enormous (∼19 kb) exon. The AcSp1 repeat sequence is exceptionally conserved between two widow species (∼94% identity) and between widows and distantly related orb-weavers (∼30% identity), consistent with a history of strong purifying selection on its amino acid sequence. Furthermore, the 16 repeats (each 371-375 amino acids long) found in black widow AcSp1 are, on average, >99% identical at the nucleotide level. A combination of stabilizing selection on amino acid sequence, selection on silent sites, and intragenic recombination likely explains the extreme homogenization of AcSp1 repeats. In addition, phylogenetic analyses of spidroin paralogs support a gene duplication event occurring concomitantly with specialization of the aciniform glands and the tubuliform glands, which synthesize egg-case silk. With repeats that are dramatically different in length and amino acid composition from dragline spidroins, our L. hesperus AcSp1 expands the knowledge base for developing silk-based biomimetic technologies.

Due to their universal presence and high sequence conservation, ribosomal RNA (rRNA) sequences are used widely in phylogenetics for inferring evolutionary relationships between microbes and in metagenomics for analyzing the composition of microbial communities. Most microbial genomes encode multiple copies of rRNA genes to supply cells with sufficient capacity for protein synthesis. These copies typically undergo concerted evolution that keeps their sequences identical, or nearly so, due to gene conversion, a type of intragenomic recombination that changes one copy of a homologous sequence to exactly match another. Widely varying rates of rRNA gene conversion have previously been estimated by comparative genomics methods and using genetic reporter assays. To more directly measure rates of rRNA intragenomic recombination, we sequenced the seven Escherichia coli rRNA operons in 15 lineages that were evolved for ~13,750 generations with frequent single-cell bottlenecks that reduce the effects of selection. We identified 38 gene conversion events and estimated an overall rate of intragenomic recombination within the 16S and 23S genes between rRNA copies of 3.6×10-4 per genome per generation or 8.6×10-6 per rRNA operon per homologous donor operon per generation. This rate varied only slightly from random expectations at different sites within the rRNA genes and between rRNA operons located at different positions in the genome. Our accurate estimate of the rate of rRNA gene conversions fills a gap in our quantitative understanding of how ribosomal sequences and other multicopy elements diversify and homogenize during microbial genome evolution [ABSTRACT FROM AUTHOR]

Background and Aims The ribosomal DNA (rDNA) gene family, encoding ribosomal RNA (rRNA), has long been regarded as an archetypal example illustrating the model of concerted evolution. However, controversy is arising, as rDNA in many eukaryotic species has been proved to be polymorphic. Here, a metagenomic strategy was applied to detect the intragenomic polymorphism as well as the evolutionary patterns of 26S rDNA across the genus Camellia. Methods Degenerate primer pairs were designed to amplify the 26S rDNA fragments from different Camellia species. The amplicons were then paired-end sequenced on the Illumina MiSeq platform. Key Results An extremely high level of rDNA polymorphism existed universally in Camellia. However, functional rDNA was still the major component of the family, and was relatively conserved among different Camellia species. Sequence variations mainly came from rRNA pseudogenes and favoured regions that are rich in GC. Specifically, some rRNA pseudogenes have existed in the genome for a long time, and have even experienced several expansion events, which has greatly enriched the abundance of rDNA polymorphism. Conclusions Camellia represents a group in which rDNA is subjected to a mixture of concerted and birth-and-death evolution. Some rRNA pseudogenes may still have potential functions. Conversely, when released from selection constraint, they can evolve in the direction of decreasing GC content and structural stability through a methylation-induced process, and finally be eliminated from the genome. [ABSTRACT FROM AUTHOR]

The Pannonian endemic species Vincetoxicum pannonicum was described from specimens collected in Hungary and occurs at only few locations. It is considered “vulnerable” according to the International Red List. The chromosome set was reported to be tetraploid, and the species was hypothesized to be an allotetraploid hybrid of the Balkan species V. fuscatum and the Adriatic species hirundinaria subsp. adriaticum. We investigated the origin of V. pannonicum using molecular phylogenetic methods by separately analyzing the multicopy nuclear ribosomal internal transcribed spacer (nrITS) and the plastid-encoded trnH-psbA DNA regions and by evaluating discrepancies between the produced gene trees. Paralogs in the nrITS region clustered in two main groups, one of which was closest to V. fuscatum, and the other included V. hirundinaria subsp. adriaticum. According to trnH-psbA sequences, V. pannonicum and V. hirundinaria subsp. adriaticum formed a single group. Our results show that V. pannonicum diversified because of hybrid speciation, in which V. fuscatum was the pollen donor. We discovered a similar placement of V. maeoticum, which suggests a further hybridization event between V. fuscatum and a species of the V. hirundinaria group. Our genome-size estimate indicates almost sixfold larger genome size in V. pannonicum compared to the maternal diploid parent, suggesting hexaploidy; however, V. pannonicum is tetraploid. This may suggest cytological diploidization in the allopolyploid V. pannonicum. We observed substantial genetic distance between V. hirundinaria subsp. adriaticum and all other subspecies of V. hirundinaria, and we therefore propose that V. adriaticum should be regarded as a separate species.

Background Species in the genus Morus (Moraceae) are deciduous woody plants of great economic importance. The classification and phylogenetic relationships of Morus, especially the abundant mulberry resources in China, is still undetermined. Internal transcribed spacer (ITS) regions are among the most widely used molecular markers in phylogenetic analyses of angiosperms. However, according to the previous phylogenetic analyses of ITS sequences, most of the mulberry accessions collected in China were grouped into the largest clade lacking for phylogenetic resolution. Compared with functional ITS sequences, ITS pseudogenes show higher sequence diversity, so they can provide useful phylogenetic information. Methods We sequenced the ITS regions and the chloroplast DNA regions TrnL-TrnF and TrnT-TrnL from 33 mulberry accessions, and performed phylogenetic analyses to explore the evolution of mulberry. Results We found ITS pseudogenes in 11 mulberry accessions. In the phylogenetic tree constructed from ITS sequences, clade B was separated into short-type sequence clades (clades 1 and 2), and a long-type sequence clade (clade 3). Pseudogene sequences were separately clustered into two pseudogroups, designated as pseudogroup 1 and pseudogroup 2. The phylogenetic tree generated from cpDNA sequences also separated clade B into two clades. Conclusions Two species were separated in clade B. The existence of three connection patterns and incongruent distribution patterns between the phylogenetic trees generated from cpDNA and ITS sequences suggested that the ITS pseudogene sequences connect with genetic information from the female progenitor. Hybridization has played important roles in the evolution of mulberry, resulting in low resolution of the phylogenetic analysis based on ITS sequences. An evolutionary pattern illustrating the evolution history of mulberry is proposed. These findings have significance for the conservation of local mulberry resources. Polyploidy, hybridization, and concerted evolution have all played the roles in the evolution of ITS sequences in mulberry. This study will expand our understanding of mulberry evolution.

Local gene duplication is a prominent mechanism of gene copy number expansion. Elucidating the mechanisms by which local duplicates arise is necessary in understanding the evolution of genomes and their host organisms. Chromosome one of Arabidopsis thaliana contains an 81-gene array subdivided into 27 triplet units (t-units), with each t-unit containing three pre-transfer RNA genes. We utilized phylogenetic tree reconstructions and comparative genomics to order the events leading to the array’s formation, and propose a model using unequal crossing-over as the primary mechanism of array formation. The model is supported by additional phylogenetic information from intergenic spacer sequences separating each t-unit, comparative analysis to an orthologous array of 12 t-units in the sister taxa Arabidopsis lyrata, and additional modeling using a stochastic simulation of orthologous array divergence. Lastly, comparative phylogenetic analysis demonstrates that the two orthologous t-unit arrays undergo concerted evolution within each taxa and are likely fluctuating in copy number under neutral evolutionary drift. These findings hold larger implications for future research concerning gene and genome evolution.

Satellite DNAs (satDNAs) are a ubiquitous feature of eukaryotic genomes and are usually the major components of constitutive heterochromatin. The 1.688 satDNA, also known as the 359 bp satellite, is one of the most abundant repetitive sequences in Drosophila melanogaster and has been linked to several different biological functions. We investigated the presence and evolution of the 1.688 satDNA in 16 Drosophila genomes. We find that the 1.688 satDNA family is much more ancient than previously appreciated, being shared among part of the melanogaster group that diverged from a common ancestor ~27 Mya. We found that the 1.688 satDNA family has two major subfamilies spread throughout Drosophila phylogeny (~360 bp and ~190 bp). Phylogenetic analysis of ~10,000 repeats extracted from 14 of the species revealed that the 1.688 satDNA family is present within heterochromatin and euchromatin. A high number of euchromatic repeats are gene proximal, suggesting the potential for local gene regulation. Notably, heterochromatic copies display concerted evolution and a species-specific pattern, whereas euchromatic repeats display a more typical evolutionary pattern, suggesting that chromatin domains may influence the evolution of these sequences. Overall, our data indicate the 1.688 satDNA as the most perduring satDNA family described in Drosophila phylogeny to date. Our study provides a strong foundation for future work on the functional roles of 1.688 satDNA across many Drosophila species. [ABSTRACT FROM AUTHOR]

Spiders produce a diversity of silk fibers from multiple morphologically distinct silk glands for specific tasks, and each silk type primarily composed of one or more particular silk proteins encoded by silk gene family members believed to generated by duplication and divergence of ancient silk genes. Egg case silks spun from tubuliform glands are used to construct the tough outer structure of egg cases, are important for their reproduction. Here we present two novel complete TuSp1 sequences from orb weaving spider Araneus ventricosus. Alignment of the two spidroin iterated repeats showed both extreme intragenic homogenization. The pairwise Ka/Ks analysis revealed the terminal and repetitive regions for three TuSp1 loci including the reported TuSp1 gene are all under purifying selection. Phylogenetic analysis showed the two new TuSp1 variants could derive from recent duplication events. [ABSTRACT FROM AUTHOR]