The recent evolution of human L1 retrotransposons

L1 elements are the most successful retrotransposons in mammals and are responsible for at least 30% of human DNA. Far from being indolent genomic parasites, L1 elements have evolved and amplified rapidly during human evolution. Indeed during just the last 25 million years (MY) five distinct L1 families have emerged and generated tens of thousands of copies. The most recently evolved human specific L1 family is currently active and L1 copies have been accumulating in the human genome at about the same rate per generation as the currently active L1 families in Old World rats and mice. At times during the last 25 MY L1 activity constituted a significant enough genetic load to be subject to negative selection. During these same times, and in apparent response to the host, L1 underwent adaptive evolution. Understanding the molecular basis for these evolutionary changes should help illuminate one of the least understood but most important aspects of L1 biology, namely the extent and nature of the interaction between L1 and its host.     

Phylogenomic analysis of chromoviruses

Genome sequences of model organisms provide a unique opportunity to obtain insight into the complete diversity of any transposable element (TE) group. A limited number of chromoviruses, the chromodomain containing genus of Metaviridae, is known from plant, fungal and vertebrate genomes. By searching diverse eukaryotic genome databases, we have found a surprisingly large number of new, structurally intact and highly conserved chromoviral elements, greatly exceeding the number of previously known chromoviruses. In this study, we examined the diversity, origin and evolution of chromoviruses in Eukaryota. Chromoviral diversity in plants, fungi and vertebrates, as shown by phylogenetic analyses, was found to be much greater than previously expected. A novel centromere-specific chromoviral lineage was found to be widespread and highly conserved in all seed plants. The age of chromoviruses has been significantly extended by finding their representatives in the most basal plant lineages (green and red algae), in Heterokonta (oomycetes) and in Cercozoa (plasmodiophorids). The evolutionary origin of chromoviruses has been found to be no earlier than in Cercozoa, since none can be found in the basal eukaryotic lineages, despite the extensive genome data. The evolutionary dynamics of chromoviruses can be explained by a strict vertical transmission in plants and fungi, while in Metazoa it is more complex. The currently available genome data clearly show that chromoviruses are the most widespread and one of the oldest Metaviridae clade.     

Phylogenomic analysis of type 1 NADH:Quinone oxidoreductase

We performed phylogenomic analysis of the catalytic core of NADH:quinone oxidoreductases of type 1 (NDH-1). Analysis of phylogenetic trees, as constructed for the core subunits of NDH-1, revealed fundamental differences in their topologies. In the case of four putatively homologous ion-carrying membrane subunits, the trees for the NuoH and NuoN subunits contained separate archaeal clades, whereas subunits NuoL and NuoM were characterized by multiple archaeal clades spread among bacterial branches. Large, separate clades, which united sequences belonging to different archaeal subdomains, were also found for cytoplasmic subunits NuoD and NuoB, homologous to the large and small subunits of nickel-iron hydrogenases. A smaller such clade was also shown for subunit NuoC. Based on these data, we suggest that the ancestral NDH-1 complex could be present already at the stage of the Last Universal Cellular Ancestor (LUCA). Ancestral forms of membrane subunits NuoN and NuoH and cytoplasmic subunits NuoD, NuoB, and, perhaps NuoC, may have formed a membrane complex that operated as an ion-translocating membrane hydrogenase. After the complex attained the ability to reduce membrane quinones, gene duplications could yield the subunits NuoL and NuoM, which enabled translocation of additional ions.     

Phylogenomic analysis of the uracil-DNA glycosylase superfamily

The spontaneous deamination of cytosine produces uracil mispaired with guanine in DNA, which will produce a mutation, unless repaired. In all domains of life, uracil-DNA glycosylases (UDGs) are responsible for the elimination of uracil from DNA. Thus, UDGs contribute to the integrity of the genetic information and their loss results in mutator phenotypes. We are interested in understanding the role of UDG genes in the evolutionary variation of the rate and the spectrum of spontaneous mutations. To this end, we determined the presence or absence of the five main UDG families in more than 1,000 completely sequenced genomes and analyzed their patterns of gene loss and gain in eubacterial lineages. We observe nonindependent patterns of gene loss and gain between UDG families in Eubacteria, suggesting extensive functional overlap in an evolutionary timescale. Given that UDGs prevent transitions at G:C sites, we expected the loss of UDG genes to bias the mutational spectrum toward a lower equilibrium G + C content. To test this hypothesis, we used phylogenetically independent contrasts to compare the G + C content at intergenic and 4-fold redundant sites between lineages where UDG genes have been lost and their sister clades. None of the main UDG families present in Eubacteria was associated with a higher G + C content at intergenic or 4-fold redundant sites. We discuss the reasons of this negative result and report several features of the evolution of the UDG superfamily with implications for their functional study. uracil-DNA glycosylase, mutation rate evolution, mutational bias, GC content, DNA repair, mutator gene.     

Identifying related L1 retrotransposons by analyzing 3' transduced sequences

Background  

A large fraction of the human genome is attributable to L1 retrotransposon sequences. Not only do L1s themselves make up a significant portion of the genome, but L1-encoded proteins are thought to be responsible for the transposition of other repetitive elements and processed pseudogenes. In addition, L1s can mobilize non-L1, 3'-flanking DNA in a process called 3' transduction. Using computational methods, we collected DNA sequences from the human genome for which we have high confidence of their mobilization through L1-mediated 3' transduction.     

Phylogenomic analysis of the cystatin superfamily in eukaryotes and prokaryotes

Background  

The cystatin superfamily comprises cysteine protease inhibitors that play key regulatory roles in protein degradation processes. Although they have been the subject of many studies, little is known about their genesis, evolution and functional diversification. Our aim has been to obtain a comprehensive insight into their origin, distribution, diversity, evolution and classification in Eukaryota, Bacteria and Archaea.     

Phylogenomic analysis clarifies the evolutionary origin of Coffea arabica

Interspecific hybridization events have played a major role in plant speciation, yet the evolutionary origin of hybrid species often remains enigmatic. Here, we inferred the evolutionary origin of the allotetraploid species Coffea arabica, which is widely cultivated for Arabica coffee production. We estimated genetic distances between C. arabica and all species that are known to be closely related to C. arabica using genotyping-by-sequencing (GBS) data. In addition, we reconstructed a time-calibrated multilabeled phylogenetic tree of 24 species to estimate the age of the C. arabica hybridization event. Ancestral states of self-compatibility were also inferred to shed new light on the evolution of self-compatibility in Coffea. Coffea canephora and C. eugenioides were confirmed as the putative progenitor species of C. arabica. These species most likely hybridized between 1.08 million and 543 000 years ago, coinciding with periods of environmental upheaval, which may have induced range shifts of the progenitor species that facilitated the emergence of C. arabica.     

Phylogenomic analysis of echinoderm class relationships supports Asterozoa

While some aspects of the phylogeny of the five living echinoderm classes are clear, the position of the ophiuroids (brittlestars) relative to asteroids (starfish), echinoids (sea urchins) and holothurians (sea cucumbers) is controversial. Ophiuroids have a pluteus-type larva in common with echinoids giving some support to an ophiuroid/echinoid/holothurian clade named Cryptosyringida. Most molecular phylogenetic studies, however, support an ophiuroid/asteroid clade (Asterozoa) implying either convergent evolution of the pluteus or reversals to an auricularia-type larva in asteroids and holothurians. A recent study of 10 genes from four of the five echinoderm classes used ‘phylogenetic signal dissection’ to separate alignment positions into subsets of (i) suboptimal, heterogeneously evolving sites (invariant plus rapidly changing) and (ii) the remaining optimal, homogeneously evolving sites. Along with most previous molecular phylogenetic studies, their set of heterogeneous sites, expected to be more prone to systematic error, support Asterozoa. The homogeneous sites, in contrast, support an ophiuroid/echinoid grouping, consistent with the cryptosyringid clade, leading them to posit homology of the ophiopluteus and echinopluteus. Our new dataset comprises 219 genes from all echinoderm classes; analyses using probabilistic Bayesian phylogenetic methods strongly support Asterozoa. The most reliable, slowly evolving quartile of genes also gives highest support for Asterozoa; this support diminishes in second and third quartiles and the fastest changing quartile places the ophiuroids close to the root. Using phylogenetic signal dissection, we find heterogenous sites support an unlikely grouping of Ophiuroidea + Holothuria while homogeneous sites again strongly support Asterozoa. Our large and taxonomically complete dataset finds no support for the cryptosyringid hypothesis; in showing strong support for the Asterozoa, our preferred topology leaves the question of homology of pluteus larvae open.     

RNase L restricts the mobility of engineered retrotransposons in cultured human cells

Retrotransposons are mobile genetic elements, and their mobility can lead to genomic instability. Retrotransposon insertions are associated with a diverse range of sporadic diseases, including cancer. Thus, it is not a surprise that multiple host defense mechanisms suppress retrotransposition. The 2′,5′-oligoadenylate (2-5A) synthetase (OAS)-RNase L system is a mechanism for restricting viral infections during the interferon antiviral response. Here, we investigated a potential role for the OAS-RNase L system in the restriction of retrotransposons. Expression of wild type (WT) and a constitutively active form of RNase L (NΔ385), but not a catalytically inactive RNase L mutant (R667A), impaired the mobility of engineered human LINE-1 (L1) and mouse intracisternal A-type particle retrotransposons in cultured human cells. Furthermore, WT RNase L, but not an inactive RNase L mutant (R667A), reduced L1 RNA levels and subsequent expression of the L1-encoded proteins (ORF1p and ORF2p). Consistently, confocal immunofluorescent microscopy demonstrated that WT RNase L, but not RNase L R667A, prevented formation of L1 cytoplasmic foci. Finally, siRNA-mediated depletion of endogenous RNase L in a human ovarian cancer cell line (Hey1b) increased the levels of L1 retrotransposition by ∼2-fold. Together, these data suggest that RNase L might function as a suppressor of structurally distinct retrotransposons.     

Phylogenomic analysis of the receptor-like proteins of rice and Arabidopsis

The tomato (Lycopersicon esculentum) Cf-9 resistance gene encodes the first characterized member of the plant receptor-like protein (RLP) family. Other RLPs such as CLAVATA2 and TOO MANY MOUTHS are known to regulate development. The domain structure of RLPs consists of extracellular leucine-rich repeats, a transmembrane helix, and a short cytoplasmic region. Here, we identify 90 RLPs in rice (Oryza sativa) and compare them with functionally characterized RLPs from different plant species and with 56 Arabidopsis (Arabidopsis thaliana) RLPs, including the downy mildew resistance protein RPP27. Many RLPs cluster into four distinct superclades, three of which include RLPs known to be involved in plant defense. Sequence comparisons reveal diagnostic amino acid residues that may specify different molecular functions in different RLP subtypes. This analysis of rice RLPs thus identified at least 73 candidate resistance genes and four genes potentially involved in development. Due to the synteny between rice and other Gramineae, this analysis should provide valuable tools for experimental studies in rice and other cereals.     

Phylogenomic analysis reveals dynamic evolutionary history of the Drosophila heterochromatin protein 1 (HP1) gene family

Heterochromatin is the gene-poor, satellite-rich eukaryotic genome compartment that supports many essential cellular processes. The functional diversity of proteins that bind and often epigenetically define heterochromatic DNA sequence reflects the diverse functions supported by this enigmatic genome compartment. Moreover, heterogeneous signatures of selection at chromosomal proteins often mirror the heterogeneity of evolutionary forces that act on heterochromatic DNA. To identify new such surrogates for dissecting heterochromatin function and evolution, we conducted a comprehensive phylogenomic analysis of the Heterochromatin Protein 1 gene family across 40 million years of Drosophila evolution. Our study expands this gene family from 5 genes to at least 26 genes, including several uncharacterized genes in Drosophila melanogaster. The 21 newly defined HP1s introduce unprecedented structural diversity, lineage-restriction, and germline-biased expression patterns into the HP1 family. We find little evidence of positive selection at these HP1 genes in both population genetic and molecular evolution analyses. Instead, we find that dynamic evolution occurs via prolific gene gains and losses. Despite this dynamic gene turnover, the number of HP1 genes is relatively constant across species. We propose that karyotype evolution drives at least some HP1 gene turnover. For example, the loss of the male germline-restricted HP1E in the obscura group coincides with one episode of dramatic karyotypic evolution, including the gain of a neo-Y in this lineage. This expanded compendium of ovary- and testis-restricted HP1 genes revealed by our study, together with correlated gain/loss dynamics and chromosome fission/fusion events, will guide functional analyses of novel roles supported by germline chromatin.     

Molecular evolution of two lineages of L1 (LINE-1) retrotransposons in the california mouse,Peromyscus californicus.

The large number of L1 [long interspersed elements (LINE)-1] sequences found in the genome is due to the insertion of copies of the retrotransposon over evolutionary time. The majority of copies appear to be replicates of a few active, or "master" templates. A continual replacement of master templates over time gives rise to lineages distinguishable by their own unique set of shared-sequence variants. A previous analysis of L1 sequences in deer mice, Peromyscus maniculatus and P. leucopus, revealed two active L1 lineages, marked by different rates of evolution, whose most recent common ancestor predates the expansion of the Peromyscus species. Here we exploit lineage-specific, shared-sequence variants to reveal a paucity of Lineage 2 sequences in at least one species, P. californicus. The dearth of Lineage 2 copies in P. californicus suggests that Lineage 2 may have been unproductive until after the most recent common ancestor of P. californicus and P. maniculatus. We also show that Lineage 1 appears to have a higher rate of evolution in P. maniculatus relative to either P. californicus or P. leucopus. As a phylogenetic tool, L1 lineage-specific variants support a close affinity between P. californicus and P. eremicus relative to the other species examined.