Excision of Helitron Transposons in Maize

Helitrons are novel transposons discovered by bioinformatic analysis of eukaryotic genome sequences. They are believed to move by rolling circle (RC) replication because their predicted transposases are homologous to those of bacterial RC transposons. We report here evidence of somatic Helitron excision in maize, an unexpected finding suggesting that Helitrons can exhibit an excisive mode of transposition.     

Evolution of Genome Size and Complexity in Pinus

Background

Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood.

Methodology/Principal Findings

Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA.

Conclusions/Significance

Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes.     

Complete Genome Analysis of Three Acinetobacter baumannii Clinical Isolates in China for Insight into the Diversification of Drug Resistance Elements

Background

The emergence and rapid spreading of multidrug-resistant Acinetobacter baumannii strains has become a major health threat worldwide. To better understand the genetic recombination related with the acquisition of drug-resistant elements during bacterial infection, we performed complete genome analysis on three newly isolated multidrug-resistant A. baumannii strains from Beijing using next-generation sequencing technology.

Methodologies/Principal Findings

Whole genome comparison revealed that all 3 strains share some common drug resistant elements including carbapenem-resistant bla _OXA-23 and tetracycline (tet) resistance islands, but the genome structures are diversified among strains. Various genomic islands intersperse on the genome with transposons and insertions, reflecting the recombination flexibility during the acquisition of the resistant elements. The blood-isolated BJAB07104 and ascites-isolated BJAB0868 exhibit high similarity on their genome structure with most of the global clone II strains, suggesting these two strains belong to the dominant outbreak strains prevalent worldwide. A large resistance island (RI) of about 121-kb, carrying a cluster of resistance-related genes, was inserted into the ATPase gene on BJAB07104 and BJAB0868 genomes. A 78-kb insertion element carrying tra-locus and bla _OXA-23 island, can be either inserted into one of the tniB gene in the 121-kb RI on the chromosome, or transformed to conjugative plasmid in the two BJAB strains. The third strains of this study, BJAB0715, which was isolated from spinal fluid, exhibit much more divergence compared with above two strains. It harbors multiple drug-resistance elements including a truncated AbaR-22-like RI on its genome. One of the unique features of this strain is that it carries both bla _OXA-23 and bla _OXA-58 genes on its genome. Besides, an Acinetobacter lwoffii adeABC efflux element was found inserted into the ATPase position in BJAB0715.

Conclusions

Our comparative analysis on currently completed Acinetobacter baumannii genomes revealed extensive and dynamic genome organizations, which may facilitate the bacteria to acquire drug-resistance elements into their genomes.     

Interspecies Insertion Polymorphism Analysis Reveals Recent Activity of Transposable Elements in Extant Coelacanths

Coelacanths are lobe-finned fish represented by two extant species, Latimeria chalumnae in South Africa and Comoros and L. menadoensis in Indonesia. Due to their intermediate phylogenetic position between ray-finned fish and tetrapods in the vertebrate lineage, they are of great interest from an evolutionary point of view. In addition, extant specimens look similar to 300 million-year-old fossils; because of their apparent slowly evolving morphology, coelacanths have been often described as « living fossils ». As an underlying cause of such a morphological stasis, several authors have proposed a slow evolution of the coelacanth genome. Accordingly, sequencing of the L. chalumnae genome has revealed a globally low substitution rate for protein-coding regions compared to other vertebrates. However, genome and gene evolution can also be influenced by transposable elements, which form a major and dynamic part of vertebrate genomes through their ability to move, duplicate and recombine. In this work, we have searched for evidence of transposition activity in coelacanth genomes through the comparative analysis of orthologous genomic regions from both Latimeria species. Comparison of 5.7 Mb (0.2%) of the L. chalumnae genome with orthologous Bacterial Artificial Chromosome clones from L. menadoensis allowed the identification of 27 species-specific transposable element insertions, with a strong relative contribution of CR1 non-LTR retrotransposons. Species-specific homologous recombination between the long terminal repeats of a new coelacanth endogenous retrovirus was also detected. Our analysis suggests that transposon activity is responsible for at least 0.6% of genome divergence between both Latimeria species. Taken together, this study demonstrates that coelacanth genomes are not evolutionary inert: they contain recently active transposable elements, which have significantly contributed to post-speciation genome divergence in Latimeria.     

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

Background

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.

Results

We examine the relationship between intron length and levels of evolutionary constraint by analyzing inter-specific divergence at 225 intron fragments in Drosophila melanogaster and Drosophila simulans, sampled from a broad distribution of intron lengths. We document a strongly negative correlation between intron length and divergence. Interestingly, we also find that divergence in introns is negatively correlated with GC content. This relationship does not account for the correlation between intron length and divergence, however, and may simply reflect local variation in mutational rates or biases.

Conclusion

Short introns make up only a small fraction of total intronic DNA in the genome. Our finding that long introns evolve more slowly than average implies that, while the majority of introns in the Drosophila genome may experience little or no selective constraint, most intronic DNA in the genome is likely to be evolving under considerable constraint. Our results suggest that functional elements may be ubiquitous within longer introns and that these introns may have a more general role in regulating gene expression than previously appreciated. Our finding that GC content and divergence are negatively correlated in introns has important implications for the interpretation of the correlation between divergence and levels of codon bias observed in Drosophila.     

Organization and evolution of transposable elements along the bread wheat chromosome 3B

Background

The 17 Gb bread wheat genome has massively expanded through the proliferation of transposable elements (TEs) and two recent rounds of polyploidization. The assembly of a 774 Mb reference sequence of wheat chromosome 3B provided us with the opportunity to explore the impact of TEs on the complex wheat genome structure and evolution at a resolution and scale not reached so far.

Results

We develop an automated workflow, CLARI-TE, for TE modeling in complex genomes. We delineate precisely 56,488 intact and 196,391 fragmented TEs along the 3B pseudomolecule, accounting for 85% of the sequence, and reconstruct 30,199 nested insertions. TEs have been mostly silent for the last one million years, and the 3B chromosome has been shaped by a succession of bursts that occurred between 1 to 3 million years ago. Accelerated TE elimination in the high-recombination distal regions is a driving force towards chromosome partitioning. CACTAs overrepresented in the high-recombination distal regions are significantly associated with recently duplicated genes. In addition, we identify 140 CACTA-mediated gene capture events with 17 genes potentially created by exon shuffling and show that 19 captured genes are transcribed and under selection pressure, suggesting the important role of CACTAs in the recent wheat adaptation.

Conclusion

Accurate TE modeling uncovers the dynamics of TEs in a highly complex and polyploid genome. It provides novel insights into chromosome partitioning and highlights the role of CACTA transposons in the high level of gene duplication in wheat.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0546-4) contains supplementary material, which is available to authorized users.     

Ancient DNA and the tropics: a rodent's tale

Most genetic studies of Holocene fauna have been performed with ancient samples from dry and cold regions, in which preservation of fossils is facilitated and molecular damage is reduced. Ancient DNA work from tropical regions has been precluded owing to factors that limit DNA preservation (e.g. temperature, hydrolytic damage). We analysed ancient DNA from rodent jawbones identified as Ototylomys phyllotis, found in Holocene and Late Pleistocene stratigraphic layers from Loltún, a humid tropical cave located in the Yucatan peninsula. We extracted DNA and amplified six short overlapping fragments of the cytochrome b gene, totalling 666 bp, which represents an unprecedented success considering tropical ancient DNA samples. We performed genetic, phylogenetic and divergence time analyses, combining sequences from ancient and modern O. phyllotis, in order to assess the ancestry of the Loltún samples. Results show that all ancient samples fall into a unique clade that diverged prior to the divergence of the modern O. phyllotis, supporting it as a distinct Pleistocene form of the Ototylomys genus. Hence, this rodent''s tale suggests that the sister group to modern O. phyllotis arose during the Miocene–Pliocene, diversified during the Pleistocene and went extinct in the Holocene.     

Pif1 Helicases and the Evidence for a Prokaryotic Origin of Helitrons

Helitrons are the only group of rolling-circle transposons that encode a transposase with a helicase domain (Hel), which belongs to the Pif1 family. Because Pif1 helicases are important components of eukaryotic genomes, it has been suggested that Hel domains probably originated after a host eukaryotic Pif1 gene was captured by a Helitron ancestor. However, the few analyses exploring the evolution of Helitron transposases (RepHel) have focused on its Rep domain, which is also present in other mobile genetic elements. Here, we used phylogenetic and nonmetric multidimensional scaling analyses to investigate the relationship between Hel domains and Pif1-like helicases from a variety of organisms. Our results reveal that Hel domains are only distantly related to genomic helicases from eukaryotes and prokaryotes, and thus are unlikely to have originated from a captured Pif1 gene. Based on this evidence, and on recent studies indicating that Rep domains are more closely related to rolling-circle plasmids and phages, we suggest that Helitrons are descendants of a RepHel-encoding prokaryotic plasmid element that invaded eukaryotic genomes before the radiation of its major groups. We discuss how a Pif1-like helicase domain might have favored the transposition of Helitrons in eukaryotes beyond simply unwinding DNA intermediates. Finally, we demonstrate that some examples in the literature describing genomic helicases from eukaryotes actually consist of Hel domains from Helitrons, a finding that underscores how transposons can hamper the analysis of eukaryotic genes. This investigation also revealed that two groups of land plants appear to have lost genomic Pif1 helicases independently.     

Palindromic repeated sequences (PRSs) in the mitochondrial genome of rice: evidence for their insertion after divergence of the genus Oryza from the other Gramineae

We have identified a family of small repeated sequences (from 60 to 66 bp in length) in the mitochondrial genome of rice (Oryza sativa cv. Nipponbare). There are at least ten copies of these sequences and they are distributed throughout the mitochondrial genome. Each is potentially capable of forming a stem-and-loop structure and we have designated them PRSs (palindromic repeated sequences). Their features are reminiscent of the small dispersed repeats in the mitochondrial DNA (mtDNA) of some lower eukaryotes, such as Saccharomyces cerevisiae, Neurospora crassa and Chlamydomonas reinhardtii. Some of the PRSs of rice mtDNA are located in the intron of the gene for ribosomal protein S3 (rps3) and in the flanking sequence of the gene for chloroplast-like tRNA^Asn (trnN). An analysis of PCR-amplified fragments of these regions from the DNA of some Gramineae suggests that the PRSs were inserted into these regions of the Oryza mtDNA after the divergence of Oryza from the other Gramineae.     

Quantifying the relative contributions of the X chromosome,autosomes, and mitochondrial genome to local adaptation*

During local adaptation with gene flow, some regions of the genome are inherently more responsive to selection than others. Recent theory predicts that X‐linked genes should disproportionately contribute to local adaptation relative to other genomic regions, yet this prediction remains to be tested. We carried out a multigeneration crossing scheme, using two cline‐end populations of Drosophila melanogaster, to estimate the relative contributions of the X chromosome, autosomes, and mitochondrial genome to divergence in four traits involved in local adaptation (wing size, resistance to heat, desiccation, and starvation stresses). We found that the mitochondrial genome and autosomes contributed significantly to clinal divergence in three of the four traits. In contrast, the X made no significant contribution to divergence in these traits. Given the small size of the mitochondrial genome, our results indicate that it plays a surprisingly large role in clinal adaptation. In contrast, the X, which represents roughly 20% of the Drosophila genome, contributes negligibly—a pattern that conflicts with theoretical predictions. These patterns reinforce recent work implying a central role of mitochondria in climatic adaptation, and suggest that different genomic regions may play fundamentally different roles in processes of divergence with gene flow.     

Mono-allelic retrotransposon insertion addresses epigenetic transcriptional repression in human genome

Background

Retrotransposons have been extensively studied in plants and animals and have been shown to have an impact on human genome dynamics and evolution. Their ability to move within genomes gives retrotransposons to affect genome instability.

Methods

we examined the polymorphic inserted AluYa5, evolutionary young Alu, in the progesterone receptor gene to determine the effects of Alu insertion on molecular environment. We used mono-allelic inserted cell lines which carry both Alu-present and Alu-absent alleles. To determine the epigenetic change and gene expression, we performed restriction enzyme digestion, Pyrosequencing, and Chromatin Immunoprecipitation.

Results

We observed that the polymorphic insertion of evolutionally young Alu causes increasing levels of DNA methylation in the surrounding genomic area and generates inactive histone tail modifications. Consequently the Alu insertion deleteriously inactivates the neighboring gene expression.

Conclusion

The mono-allelic Alu insertion cell line clearly showed that polymorphic inserted repetitive elements cause the inactivation of neighboring gene expression, bringing aberrant epigenetic changes.     

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (d_XY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.     

Contrasting evolutionary patterns and target specificities among three Tourist-like MITE families in the maize genome

Miniature inverted-repeat transposable elements (MITEs) are short, non autonomous DNA elements that are widespread and abundant in plant genomes. The high sequence and size conservation observed in many MITE families suggest that they have spread recently throughout their respective host genomes. Here we present a maize genome wide analysis of three Tourist-like MITE families, mPIF, and two previously uncharacterized families, ZmV1 and Zead8. We undertook a bioinformatic analysis of MITE insertion sites, developed methyl-sensitive transposon display (M-STD) assays to estimate the associated level of CpG methylation at MITE flanking regions, and conducted a population genetics approach to investigate MITE patterns of expansion. Our results reveal that the three MITE families insert into genomic regions that present specific molecular features: they are preferentially AT rich, present low level of cytosine methylation as compared to the LTR retrotransposon Grande, and target site duplications are flanked by large and conserved palindromic sequences. Moreover, the analysis of MITE distances from predicted genes shows that 73% of 263 copies are inserted at less than 5 kb from the nearest predicted gene, and copies from Zead8 family are significantly more abundant upstream of genes. By employing a population genetic approach we identified contrasting patterns of expansion among the three MITE families. All elements seem to have inserted roughly 1 million years ago but ZmV1 and Zead8 families present evidences for activity of several master copies within the last 0.4 Mya.     

Divergence of Gene Body DNA Methylation and Evolution of Plant Duplicate Genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.     

Characterization of Imcrop, a Mutator-like MITE family in the rice genome

Sequence comparisons of ammonium transporter 1?C2 genes (OsAMT1-2) in different rice accessions revealed a MITE insertion in the upstream region of the gene. The 391-bp MITE, classified as a Mutator superfamily member and named Imcrop, included terminal inverted repeat (TIR) and 9-bp target site duplication (TSD) sequences. We identified 151 Imcrop elements dispersed on 12 chromosomes of the japonica reference genome. Of these, 12.6% were found in genic regions and 33.1% were located within 1.5 kb of annotated rice genes. We constructed comparative insertion maps with 111 and 102 intact Imcrop elements in the japonica and indica reference genomes, respectively. The Imcrop elements showed relatively even distribution across all chromosomes although their frequency was higher on chromosomes 1, 3, and 4 in both genomes. Seventy seven Imcrop elements were detected in both subspecies, whereas 34 and 25 insertions were found only in the japonica or indica genome, respectively. We compared insertion polymorphisms of 19 Imcrop elements found inside genes in 48 Korean rice cultivars, consisting of 42 japonica and six Tongil-types (indica-japonica cross). Thirteen insertions were common to all cultivars indicating these elements were present before indica-japonica divergence. The six other elements showed insertion polymorphisms among accessions, showing their recent insertion history or no critical positive effect of their insertion on the rice genome.     

Genetic Insertions and Diversification of the PolB-Type DNA Polymerase (gp43) of T4-Related Phages

In Escherichia coli phage T4 and many of its phylogenetic relatives, gene 43 consists of a single cistron that encodes a PolB family (PolB-type) DNA polymerase. We describe the divergence of this phage gene and its protein product (gp43) (gene product 43) among 26 phylogenetic relatives of T4 and discuss our observations in the context of diversity among the widely distributed PolB enzymes in nature. In two T4 relatives that grow in Aeromonas salmonicida phages 44RR and 25, gene 43 is fragmented by different combinations of three distinct types of DNA insertion elements: (a) a short intercistronic untranslated sequence (IC-UTS) that splits the polymerase gene into two cistrons, 43A and 43B, corresponding to N-terminal (gp43A) and C-terminal (gp43B) protein products; (b) a freestanding homing endonuclease gene (HEG) inserted between the IC-UTS and the 43B cistron; and (c) a group I intron in the 43B cistron. Phage 25 has all three elements, whereas phage 44RR has only the IC-UTS. We present evidence that (a) the split gene of phage 44RR encodes a split DNA polymerase consisting of a complex between gp43A and gp43B subunits; (b) the putative HEG encodes a double-stranded DNA endonuclease that specifically cleaves intron-free homologues of the intron-bearing 43B site; and (c) the group I intron is a self-splicing RNA. Our results suggest that some freestanding HEGs can mediate the homing of introns that do not encode their own homing enzymes. The results also suggest that different insertion elements can converge on a polB gene and evolve into a single integrated system for lateral transfer of polB genetic material. We discuss the possible pathways for the importation of such insertion elements into the genomes of T4-related phages.