Identification of an active <Emphasis Type="Italic">Mutator</Emphasis>-like element (MULE) in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis>)

Transposable elements (TEs) represent an important fraction of plant genomes and play a significant role in gene and genome evolution. Among all TE superfamilies discovered in plants, Mutator from maize (Zea mays) is the most active and mutagenic element. Mutator-like elements (MULEs) were identified in a wide range of plants. However, only few active MULEs have been reported, and the transposition mechanism of the elements is still poorly understood. In this study, an active MULE named Os3378 was discovered in rice (Oryza sativa) by a combination of computational and experimental approaches. The four newly identified Os3378 elements share more than 98% sequence identity between each other, and all of them encode transposases without any deletion derivatives, indicating their capability of autonomous transposition. Os3378 is present in the rice species with AA genome type but is absent in other non-AA genome species. A new insertion of Os3378 was identified in a rice somaclonal mutant Z418, and the element remained active in the descendants of the mutant for more than ten generations. Both germinal and somatic excision events of Os3378 were observed, and no footprint was detected after excision. Furthermore, the occurrence of somatic excision of Os3378 appeared to be associated with plant developmental stages and tissue types. Taken together, Os3378 is a unique active element in rice, which provides a valuable resource for further studying of transposition mechanism and evolution of MULEs.     

Transposition behavior of nonautonomous a <Emphasis Type="Italic">hAT</Emphasis> superfamily transposon <Emphasis Type="Italic">nDart</Emphasis> in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Transposable elements (TEs) have a significant impact on the evolution of gene function and genome structures. An endogenous nonautonomous transposable element nDart was discovered in an albino mutant that had an insertion in the Mg-protoporphyrin IX methyltransferase gene in rice. In this study, we elucidated the transposition behavior of nDart, the frequency of nDart transposition and characterized the footprint of nDart. Novel independent nDart insertions in backcrossed progenies were detected by DNA blotting analysis. In addition, germinal excision of nDart occurred at very low frequency compared with that of somatic excision, 0–13.3%, in the nDart1-4(3-2) and nDart1-A loci by a locus-specific PCR strategy. A total of 253 clones from somatic excision at five nDart loci in 10 varieties were determined. nDart rarely caused deletions beyond target site duplication (TSD). The footprint of nDart contained few transversions of nucleotides flanking to both sides of the TSD. The predominant footprint of nDart was an 8-bp addition. Precise excision of nDart was detected at a rate of only 2.2%, which occurred at two loci among the five loci examined. Furthermore, the results in this study revealed that a highly conserved mechanism of transposition is involved between maize Ac/Ds and rice Dart/nDart, which are two-component transposon systems of the hAT superfamily transposons in plant species.     

<Emphasis Type="Italic">hATpin</Emphasis>, a Family of MITE-like <Emphasis Type="Italic">hAT</Emphasis> Mobile Elements Conserved in Diverse Plant Species that Forms Highly Stable Secondary Structures

We identified a 178 bp mobile DNA element in lettuce with characteristic CGAGC/GCTCG repeats in the subterminal regions. This element has terminal inverted repeats and 8-bp target site duplications typical of the hAT superfamily of class II mobile elements, but its small size and potential to form a single-stranded stable hairpin-like secondary structure suggest that it is related to MITE elements. In silico searches for related elements identified 252 plant sequences with 8-bp target site duplications and sequence similarity in their terminal and subterminal regions. Some of these sequences were predicted to encode transposases and may be autonomous elements; these constituted a separate clade within the phylogram of hAT transposases. We demonstrate that the CGAGC/GCTCG pentamer maximizes the hairpin stability compared to any other pentamer with the same C + G content, and the secondary structures of these elements are more stable than for most MITEs. We named these elements collectively as hATpin elements because of the hAT similarity and their hairpin structures. The nearly complete rice genome sequence and the highly advanced genome annotation allowed us to localize most rice elements and to deduce insertion preferences. hATpin elements are distributed on all chromosomes, but with significant bias for chromosomes 1 and 10 and in regions of moderate gene density. This family of class II mobile elements is found primarily in monocot species, but is also present in dicot species. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Direct determination of the effects of genotype and extreme temperature on the transposition of <Emphasis Type="Italic">roo</Emphasis> in long-term mutation accumulation lines of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Transposable elements (TEs) are mobile repetitive DNA sequences that constitute a structurally dynamic component of genomes. In order to understand the dynamics of TEs it is necessary to have information about the control of transposition and its dependence of environmental factors. After a great deal of previous work on transposition conducted on long-term mutation accumulation (MA) lines of Drosophila melanogaster started in 1987, only roo out of 16 families was found active in this genotype. Here we test the effect of the modification of the genetic background by introducing a Cy chromosome, and the effect of extreme temperature (28°C) on the transposition rate of roo. Thermal stress did not affect the transposition rate, whereas the presence of a Cy chromosome in heterozygosis lowered it. There was an excess of insertions in the X chromosome, with respect to autosomes, and in the proximal and distal regions of chromosome arms that can be interpreted as target site preference. One of the control lines became highly unstable with mean insertion and excision rates of 3.0 × 10⁻³ and 8.5 × 10⁻⁴, respectively. Instability arose spontaneously during generations of mutation accumulation, and can be attributed to “de novo” mutation. Transposition in the unstable line could be directly studied on the progeny of individual males and females, from where we deduced that transposition occurs mainly, if not exclusively, in males, with a rate of 1.125 insertions per gamete. In situ hybridization with an LTR probe showed that most excisions (12 out of 14) were precise. Our data show the prominent role of genotype in transposition control and can explain rapid turnovers in the genome without increasing the number of copies.     

Transpositional reactivation of the <Emphasis Type="Italic">Dart</Emphasis> transposon family in rice lines derived from introgressive hybridization with <Emphasis Type="Italic">Zizania latifolia</Emphasis>

Background  

It is widely recognized that interspecific hybridization may induce "genome shock", and lead to genetic and epigenetic instabilities in the resultant hybrids and/or backcrossed introgressants. A prominent component involved in the genome shock is reactivation of cryptic transposable elements (TEs) in the hybrid genome, which is often associated with alteration in the elements' epigenetic modifications like cytosine DNA methylation. We have previously reported that introgressants derived from hybridization between Oryza sativa (rice) and Zizania latifolia manifested substantial methylation re-patterning and rampant mobilization of two TEs, a copia retrotransposon Tos17 and a MITE mPing. It was not known however whether other types of TEs had also been transpositionally reactivated in these introgressants, their relevance to alteration in cytosine methylation, and their impact on expression of adjacent cellular genes.     

New insights into <Emphasis Type="Italic">Oryza</Emphasis> genome evolution: high gene colinearity and differential retrotransposon amplification

An ∼247-kb genomic region from FF genome of wild rice Oryza brachyantha, possessing the smallest Oryza genome, was compared to the orthologous ∼450-kb region from AA genome, O. sativa L. ssp. japonica. 37 of 38 genes in the orthologous regions are shared between japonica and O. brachyantha. Analyses of nucleotide substitution in coding regions suggest the two genomes diverged ∼10 million years ago. Comparisons of transposable elements (TEs) reveal that the density of DNA TEs in O. brachyantha is comparable to O. sativa; however, the density of RNA TEs is dramatically lower. The genomic fraction of RNA TEs in japonica is two times greater than in O. brachyantha. Differences, particularly in RNA TEs, in this region and in BAC end sequences from five wild and two cultivated Oryza species explain major genome size differences between sativa and brachyantha. Gene expression analyses of three ObDREB1 genes in the sequenced region indicate orthologous genes retain similar expression patterns following cold stress. Our results demonstrate that size and number of RNA TEs play a major role in genomic differentiation and evolution in Oryza. Additionally, distantly related O. brachyantha shares colinearity with O. sativa, offering opportunities to use comparative genomics to explore the genetic diversity of wild species to improve cultivated rice. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Data deposition: Sequence data from this article were deposited with GenBank Library under accession number DQ810282. Shibo Zhang and Yong Qiang Gu contributed equally to the work     

Analysis of the human <Emphasis Type="Italic">Alu</Emphasis> Ye lineage

Background  

Alu elements are short (~300 bp) interspersed elements that amplify in primate genomes through a process termed retroposition. The expansion of these elements has had a significant impact on the structure and function of primate genomes. Approximately 10 % of the mass of the human genome is comprised of Alu elements, making them the most abundant short interspersed element (SINE) in our genome. The majority of Alu amplification occurred early in primate evolution, and the current rate of Alu retroposition is at least 100 fold slower than the peak of amplification that occurred 30–50 million years ago. Alu elements are therefore a rich source of inter- and intra-species primate genomic variation.     

Molecular evolution of the clustered <Emphasis Type="Italic">MIC</Emphasis>-<Emphasis Type="Italic">3</Emphasis> multigene family of <Emphasis Type="Italic">Gossypium</Emphasis> species

The Gossypium MIC-3 (Meloidogyne Induced Cotton-3) gene family is of great interest for molecular evolutionary studies because of its uniqueness to Gossypium species, multi-gene content, clustered localization, and root-knot nematode resistance-associated features. Molecular evolution of the MIC-3 gene family was studied in 15 tetraploid and diploid Gossypium genotypes that collectively represent seven phylogenetically distinct genomes. Synonymous (d_S) and non-synonymous (d_N) nucleotide substitution rates suggest that the second of the two exons of the MIC-3 genes has been under strong positive selection pressure, while the first exon has been under strong purifying selection to preserve function. Based on nucleotide substitution rates, we conclude that MIC-3 genes are evolving by a birth-and-death process and that a ‘gene amplification’ mechanism has helped to retain all duplicate copies, which best fits with the “bait and switch” model of R-gene evolution. The data indicate MIC-3 gene duplication events occurred at various rates, once per 1 million years (MY) in the allotetraploids, once per ~2 MY in the A/F genome clade, and once per ~8 MY in the D-genome clade. Variations in the MIC-3 gene family seem to reflect evolutionary selection for increased functional stability, while also expanding the capacity to develop novel “switch” pockets for responding to diverse pests and pathogens. Such evolutionary roles are congruent with the hypothesis that members of this unique resistance gene family provide fitness advantages in Gossypium.     

Identification and genomic distribution of<Emphasis Type="Italic"> gypsy</Emphasis> like retrotransposons in<Emphasis Type="Italic"> Citrus</Emphasis> and<Emphasis Type="Italic"> Poncirus</Emphasis>

Transposable elements might be importantly involved in citrus genetic instability and genome evolution. The presence of gypsy like retrotransposons, their heterogeneity and genomic distribution in Citrus and Poncirus, have been investigated. Eight clones containing part of the POL coding region of gypsy like retrotransposons have been isolated from a commercial variety of Citrus clementina, one of the few sexual species in Citrus. Four of the eight clones might correspond to active elements given that they present all the conserved motifs described in the literature as essential for activity, no in-frame stop codon and no frame-shift mutation. High homology has been found between some of these citrus elements and retroelements within a resistance-gene cluster from potato, another from Poncirus trifoliata and two putative resistance polyproteins from rice. Nested copies of gypsy like elements are scattered along the Citrus and Poncirus genomes. The results on genomic distribution show that these elements were introduced before the divergence of both genera and evolved separately thereafter. IRAPs based on gypsy and copia types of retrotransposons seem to distribute differently, therefore gypsy based IRAPs prove a new, complementary set of molecular markers in Citrus to study and map genetic variability, especially for disease resistance. Similarly to copia-derived IRAPs, the number of copies and heterozygosity values found for gypsy derived IRAPs are lower in Poncirus than in Citrus aurantium, which is less apomictic and the most usual rootstock for clementines until 1970.Communicated by C. Möllers     

Both sense and antisense strands of the LTR of the <Emphasis Type="Italic">Schistosoma mansoni Pao</Emphasis>-like retrotransposon <Emphasis Type="Italic">Sinbad</Emphasis> drive luciferase expression

Long terminal repeat (LTR) retrotransposons, mobile genetic elements comprising substantial proportions of many eukaryotic genomes, are so named for the presence of LTRs, direct repeats about 250–600 bp in length flanking the open reading frames that encode the retrotransposon enzymes and structural proteins. LTRs include promotor functions as well as other roles in retrotransposition. LTR retrotransposons, including the Gypsy-like Boudicca and the Pao/BEL-like Sinbad elements, comprise a substantial proportion of the genome of the human blood fluke, Schistosoma mansoni. In order to deduce the capability of specific copies of Boudicca and Sinbad LTRs to function as promotors, these LTRs were investigated analytically and experimentally. Sequence analysis revealed the presence of TATA boxes, canonical polyadenylation signals, and direct inverted repeats within the LTRs of both the Boudicca and Sinbad retrotransposons. Inserted in the reporter plasmid pGL3, the LTR of Sinbad drove firefly luciferase activity in HeLa cells in its forward and inverted orientation. In contrast, the LTR of Boudicca did not drive luciferase activity in HeLa cells. The ability of the Sinbad LTR to transcribe in both its forward and inverted orientation represents one of few documented examples of bidirectional promotor function.     

Characterization of <Emphasis Type="Italic">EamaT1</Emphasis>, a member of <Emphasis Type="Italic">maT</Emphasis> family of transposable elements from the earthworm <Emphasis Type="Italic">Eisenia andrei</Emphasis> (Annelida,Oligochaeta)

The maT family is a unique clade within the Tc1-mariner superfamily, and their distribution is to date known as being limited to invertebrates. A novel transposon named EamaT1 is described from the genome of the earthworm Eisenia andrei. The full sized EamaT1 was obtained by degenerate and inverse PCR-based amplification. Sequence analysis of multiple copies of the EamaT1, which consisted of 0.9 and 1.4 kb elements, showed that the consensual EamaT1 with inverted terminal repeats (ITRs) of 69 bp was 1,422 bp long and flanked by a duplicated TA dinucleotide. The EamaT1 is present in approximately 120–250 copies per diploid genome but undergoes an inactivation process as a result of accumulating multiple mutations and is nonfunctional. The open reading frame (ORF) of the EamaT1 consensus encoding 356 amino acid sequences of transposase contained a DD37D signature and a conserved paired-like DNA binding motif for the transposition mechanism. The result of ITRs comparison confirmed their consensus terminal sequences (5′-CAGGGTG-3′) and AT-rich region on the internal bases for ITRs-transposase interaction.     

The <Emphasis Type="Italic">Caenorhabditis</Emphasis><Emphasis Type="Italic">briggsae</Emphasis> genome contains active <Emphasis Type="Italic">CbmaT1</Emphasis> and <Emphasis Type="Italic">Tcb1</Emphasis> transposons

The maT clade of transposons is a group of transposable elements intermediate in sequence and predicted protein structure to mariner and Tc transposons, with a distribution thus far limited to a few invertebrate species. We present evidence, based on searches of publicly available databases, that the nematode Caenorhabditis briggsae has several maT-like transposons, which we have designated as CbmaT elements, dispersed throughout its genome. We also describe two additional transposon sequences that probably share their evolutionary history with the CbmaT transposons. One resembles a fold back variant of a CbmaT element, with long (380-bp) inverted terminal repeats (ITRs) that show a high degree (71%) of identity to CbmaT1. The other, which shares only the 26-bp ITR sequences with one of the CbmaT variants, is present in eight nearly identical copies, but does not have a transposase gene and may therefore be cross mobilised by a CbmaT transposase. Using PCR-based mobility assays, we show that CbmaT1 transposons are capable of excising from the C. briggsae genome. CbmaT1 excised approximately 500 times less frequently than Tcb1 in the reference strain AF16, but both CbmaT1 and Tcb1 excised at extremely high frequencies in the HK105 strain. The HK105 strain also exhibited a high frequency of spontaneous induction of unc-22 mutants, suggesting that it may be a mutator strain of C. briggsae.     

Comparative polytene chromosome maps of <Emphasis Type="Italic">D. montana</Emphasis> and <Emphasis Type="Italic">D. virilis</Emphasis>

Chromosomal inversion polymorphism was characterized in Finnish Drosophila montana populations. A total of 14 polymorphic inversions were observed in Finnish D. montana of which nine had not been described before. The number of polymorphic inversions in each chromosome was not significantly different from that expected, assuming equal chance of occurrence in the euchromatic genome. There was, however, no correlation between the number of polymorphic inversions and that of fixed inversions in each chromosome. Therefore, a simple neutral model does not explain the evolutionary dynamics of inversions. Furthermore, in contrast to results obtained by others, no significant correlation was found between the two transposable elements (TEs) Penelope and Ulysses and inversion breakpoints in D. montana. This result suggests that these TEs were not involved in the creation of the polymorphic inversions seen in D. montana. A comparative analysis of D. montana and Drosophila virilis polytene chromosomes 4 and 5 was performed with D. virilis bacteriophage P1 clones, thus completing the comparative studies of the two species.     

Polymorphism of canonical and noncanonical <Emphasis Type="Italic">gypsy</Emphasis> sequences in different species of <Emphasis Type="Italic">Drosophila </Emphasis><Emphasis Type="Italic">melanogaster</Emphasis> subgroup: possible evolutionary relations

Mobile genetic elements constitute a substantial part of eukaryotic genome and play an important role in its organization and functioning. Co-evolution of retrotransposons and their hosts resulted in the establishment of control systems employing mechanisms of RNA interference that seem to be impossible to evade. However, “active” copies of endogenous retrovirus gypsy escape cellular control in some cases, while its evolutionary elder “inactive” variants do not. To clarify the evolutionary relationship between “active” and “inactive” gypsy we combined two approaches: the analysis of gypsy sequences, isolated from G32 Drosophila melanogaster strain and from different Drosophila species of the melanogaster subgroup, as well as the study of databases, available on the Internet. No signs of “intermediate” (between “active” and “inactive”) gypsy form were found in GenBank, and four full-size G32 gypsy copies demonstrated a convergence that presumably involves gene conversion. No “active” gypsy were revealed among PCR generated gypsy ORF3 sequences from the various Drosophila species indicating that “active” gypsy appeared in some population of D. melanogaster and then started to spread out. Analysis of sequences flanking gypsy variants in G32 revealed their predominantly heterochromatic location. Discrepancy between the structure of actual gypsy sites in G32 and corresponding sequences in database might indicate significant inter-strain heterochromatin diversity. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The <Emphasis Type="Italic">Salmonella enterica</Emphasis> Pan-genome

Salmonella enterica is divided into four subspecies containing a large number of different serovars, several of which are important zoonotic pathogens and some show a high degree of host specificity or host preference. We compare 45 sequenced S. enterica genomes that are publicly available (22 complete and 23 draft genome sequences). Of these, 35 were found to be of sufficiently good quality to allow a detailed analysis, along with two Escherichia coli strains (K-12 substr. DH10B and the avian pathogenic E. coli (APEC O1) strain). All genomes were subjected to standardized gene finding, and the core and pan-genome of Salmonella were estimated to be around 2,800 and 10,000 gene families, respectively. The constructed pan-genomic dendrograms suggest that gene content is often, but not uniformly correlated to serotype. Any given Salmonella strain has a large stable core, whilst there is an abundance of accessory genes, including the Salmonella pathogenicity islands (SPIs), transposable elements, phages, and plasmid DNA. We visualize conservation in the genomes in relation to chromosomal location and DNA structural features and find that variation in gene content is localized in a selection of variable genomic regions or islands. These include the SPIs but also encompass phage insertion sites and transposable elements. The islands were typically well conserved in several, but not all, isolates—a difference which may have implications in, e.g., host specificity.     

Mining <Emphasis Type="Italic">Xanthomonas</Emphasis> and <Emphasis Type="Italic">Streptomyces</Emphasis> genomes for new pectinase-encoding sequences and their heterologous expression in <Emphasis Type="Italic">Escherichia coli</Emphasis>

Microbial genome sequencing has left a legacy of annotated yet uncharacterized genes or open reading frames, activities that may have useful applications in health and/or the environment. We are interested in the discovery and characterization of potentially new pectinolytic activities for the enzymatic retting of natural bast fibers such as hemp and flax. A highlight in this study is the discovery of a cold-active pectate lyase among five pectate-lyase-encoding sequences and two polygalacturonase-encoding sequences that we have cloned from the genomes of Xanthomonas campestris pv. campestris and Streptomyces coelicolor A3(2). Heterologous expression of these sequences as active pectate lyases and polygalacturonases required their subcloning in Escherichia coli Rosetta™ cells. The most active recombinant pectate lyase (XcPL NP_638163), a cold-active pectate lyase (XcPL NP_636037), and a polygalacturonase (XcPG NP_638805) were purified to near homogeneity and their kinetic parameters were determined. A significant amount of pectin degradation products was shown to be released by the two pectate lyases but not the polygalacturonase when hemp fiber pectin was used as substrate. Results of this study showed that genome data mining, besides an economical approach to new gene acquisition, may uncover new findings such as the discovery of a cold-active pectate-lyase-encoding sequence from X. campestris, a mesophilic microorganism.     

Strong Expression and Conserved Regulation of <Emphasis Type="Italic">ACT2</Emphasis> in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> and <Emphasis Type="Italic">Physcomitrella patens</Emphasis>

Arabidopsis ACT2 represents an ancient class of vegetative plant actins and is strongly and constitutively expressed in almost all Arabidopsis sporophyte vegetative tissues. Using the beta glucuronidase report system, the studies showed that ACT2 5′ regulatory region was significantly more active than CaMV 35S promoter in Arabidopsis seedlings and gametophyte vegetative tissues of Physcomitrella patens. Its activity was also observed in rice and maize seedlings. Thus, the ACT2 5′ regulatory region could potentially serve as a strong regulator to express a transgene in divergent plant species. ACT2 5′ regulatory region contained 15 conserved sequence elements, an ancient intron in its 5′ un-translated region (5′ UTR), and a purine-rich stretch followed by a pyrimidine-rich stretch (PuPy). Mutagenesis and deletion analysis illustrated that some of the conserved sequence elements and the region containing PuPy sequences played regulatory roles in Arabidopsis. Interestingly, mutation of the conserved elements did not lead a dramatic change in the activity of ACT2 5′ regulatory region. The ancient intron in ACT2 5′ UTR was required for its strong expression in both Arabidopsis and P. patens, but did not fully function as a canonical intron. Thus, it was likely that some of the conserved sequence elements and gene structures had been preserved in ACT2 5′ regulatory region over the course of land plant evolution partly due to their functional importance. The studies provided additional evidences that identification of evolutionarily conserved features in non-coding region might be used as an efficient strategy to predict gene regulatory elements.