Plant transposable elements: where genetics meets genomics

Gene duplication, population genomics, and species-level differentiation within a tropical mountain shrub

Gene duplication leads to paralogy, which complicates the de novo assembly of genotyping-by-sequencing (GBS) data. The issue of paralogous genes is exacerbated in plants, because they are particularly prone to gene duplication events. Paralogs are normally filtered from GBS data before undertaking population genomics or phylogenetic analyses. However, gene duplication plays an important role in the functional diversification of genes and it can also lead to the formation of postzygotic barriers. Using populations and closely related species of a tropical mountain shrub, we examine 1) the genomic differentiation produced by putative orthologs, and 2) the distribution of recent gene duplication among lineages and geography. We find high differentiation among populations from isolated mountain peaks and species-level differentiation within what is morphologically described as a single species. The inferred distribution of paralogs among populations is congruent with taxonomy and shows that GBS could be used to examine recent gene duplication as a source of genomic differentiation of nonmodel species.

Illumina TruSeq synthetic long-reads empower de novo assembly and resolve complex, highly-repetitive transposable elements

High-throughput DNA sequencing technologies have revolutionized genomic analysis, including the de novo assembly of whole genomes. Nevertheless, assembly of complex genomes remains challenging, in part due to the presence of dispersed repeats which introduce ambiguity during genome reconstruction. Transposable elements (TEs) can be particularly problematic, especially for TE families exhibiting high sequence identity, high copy number, or complex genomic arrangements. While TEs strongly affect genome function and evolution, most current de novo assembly approaches cannot resolve long, identical, and abundant families of TEs. Here, we applied a novel Illumina technology called TruSeq synthetic long-reads, which are generated through highly-parallel library preparation and local assembly of short read data and which achieve lengths of 1.5–18.5 Kbp with an extremely low error rate (0.03% per base). To test the utility of this technology, we sequenced and assembled the genome of the model organism Drosophila melanogaster (reference genome strain y; cn, bw, sp) achieving an N50 contig size of 69.7 Kbp and covering 96.9% of the euchromatic chromosome arms of the current reference genome. TruSeq synthetic long-read technology enables placement of individual TE copies in their proper genomic locations as well as accurate reconstruction of TE sequences. We entirely recovered and accurately placed 4,229 (77.8%) of the 5,434 annotated transposable elements with perfect identity to the current reference genome. As TEs are ubiquitous features of genomes of many species, TruSeq synthetic long-reads, and likely other methods that generate long-reads, offer a powerful approach to improve de novo assemblies of whole genomes.

Transposon-based tagging: IRAP, REMAP, and iPBS

Retrotransposons are a major component of virtually all eukaryotic genomes, which makes them useful as molecular markers. Various molecular marker systems have been developed that exploit the ubiquitous nature of these genetic elements and their property of stable integration into dispersed chromosomal loci that are polymorphic within species. To detect polymorphisms for retrotransposon insertions, marker systems generally rely on PCR amplification between the retrotransposon termini and some component of flanking genomic DNA. The main methods of IRAP, REMAP, RBIP, and SSAP all detect the polymorphic sites at which the retrotransposon DNA is integrated into the genome. Marker systems exploiting these methods can be easily developed and are inexpensively deployed in the absence of extensive genome sequence data. Here, we describe protocols for the IRAP, REMAP, and iPBS techniques, including methods for PCR amplification with a single primer or with two primers, and agarose gel electrophoresis of the product using optimal electrophoresis buffers; we also describe iPBS techniques for the rapid isolation of retrotransposon termini and full-length elements.

Retrotransposon-based molecular markers for analysis of genetic diversity within the Genus Linum

SSAP method was used to study the genetic diversity of 22 Linum species from sections Linum, Adenolinum, Dasylinum, Stellerolinum, and 46 flax cultivars. All the studied flax varieties were distinguished using SSAP for retrotransposons FL9 and FL11. Thus, the validity of SSAP method was demonstrated for flax marking, identification of accessions in genebank collections, and control during propagation of flax varieties. Polymorphism of Fl1a, Fl1b, and Cassandra insertions were very low in flax varieties, but these retrotransposons were successfully used for the investigation of Linum species. Species clusterization based on SSAP markers was in concordance with their taxonomic division into sections Dasylinum, Stellerolinum, Adenolinum, and Linum. All species of sect. Adenolinum clustered apart from species of sect. Linum. The data confirmed the accuracy of the separation in these sections. Members of section Linum are not as closely related as members of other sections, so taxonomic revision of this section is desirable. L. usitatissimum accessions genetically distant from modern flax cultivars were revealed in our work. These accessions are of utmost interest for flax breeding and introduction of new useful traits into flax cultivars. The chromosome localization of Cassandra retrotransposon in Linum species was determined.

Analysis of correlation structures in the Synechocystis PCC6803 genome

Transfer of nucleotide strings in the Synechocystis sp. PCC6803 genome is investigated to exhibit periodic and non-periodic correlation structures by using the recurrence plot method and the phase space reconstruction technique. The periodic correlation structures are generated by periodic transfer of several substrings in long periodic or non-periodic nucleotide strings embedded in the coding regions of genes. The non-periodic correlation structures are generated by non-periodic transfer of several substrings covering or overlapping with the coding regions of genes. In the periodic and non-periodic transfer, some gaps divide the long nucleotide strings into the substrings and prevent their global transfer. Most of the gaps are either the replacement of one base or the insertion/reduction of one base. In the reconstructed phase space, the points generated from two or three steps for the continuous iterative transfer via the second maximal distance can be fitted by two lines. It partly reveals an intrinsic dynamics in the transfer of nucleotide strings. Due to the comparison of the relative positions and lengths, the substrings concerned with the non-periodic correlation structures are almost identical to the mobile elements annotated in the genome. The mobile elements are thus endowed with the basic results on the correlation structures.

Polyploidy and novelty: Gottlieb's legacy

Nearly four decades ago, Roose & Gottlieb (Roose & Gottlieb 1976 Evolution 30, 818-830. (doi:10.2307/2407821)) showed that the recently derived allotetraploids Tragopogon mirus and T. miscellus combined the allozyme profiles of their diploid parents (T. dubius and T. porrifolius, and T. dubius and T. pratensis, respectively). This classic paper addressed the link between genotype and biochemical phenotype and documented enzyme additivity in allopolyploids. Perhaps more important than their model of additivity, however, was their demonstration of novelty at the biochemical level. Enzyme multiplicity-the production of novel enzyme forms in the allopolyploids-can provide an extensive array of polymorphism for a polyploid individual and may explain, for example, the expanded ranges of polyploids relative to their diploid progenitors. In this paper, we extend the concept of evolutionary novelty in allopolyploids to a range of genetic and ecological features. We observe that the dynamic nature of polyploid genomes-with alterations in gene content, gene number, gene arrangement, gene expression and transposon activity-may generate sufficient novelty that every individual in a polyploid population or species may be unique. Whereas certain combinations of these features will undoubtedly be maladaptive, some unique combinations of newly generated variation may provide tremendous evolutionary potential and adaptive capabilities.

Retrotransposon-based molecular markers for assessment of genomic diversity

Retrotransposons (RTs) are major components of most eukaryotic genomes. They are ubiquitous, dispersed throughout the genome, and their abundance correlates with genome size. Their copy-and-paste lifestyle in the genome consists of three molecular steps involving transcription of an RNA copy from the genomic RT, followed by reverse transcription to generate cDNA, and finally, reintegration into a new location in the genome. This process leads to new genomic insertions without excision of the original element. The target sites of insertions are relatively random and independent for different taxa; however, some elements cluster together in 'repeat seas' or have a tendency to cluster around the centromeres and telomeres. The structure and copy number of retrotransposon families are strongly influenced by the evolutionary history of the host genome. Molecular markers play an essential role in all aspects of genetics and genomics, and RTs represent a powerful tool compared with other molecular and morphological markers. All features of integration activity, persistence, dispersion, conserved structure and sequence motifs, and high copy number suggest that RTs are appropriate genomic features for building molecular marker systems. To detect polymorphisms for RTs, marker systems generally rely on the amplification of sequences between the ends of the RT, such as (long-terminal repeat)-retrotransposons and the flanking genomic DNA. Here, we review the utility of some commonly used PCR retrotransposon-based molecular markers, including inter-primer binding sequence (IPBS), sequence-specific amplified polymorphism (SSAP), retrotransposon-based insertion polymorphism (RBIP), inter retrotransposon amplified polymorphism (IRAP), and retrotransposon-microsatellite amplified polymorphism (REMAP).

Mutation of a major CG methylase in rice causes genome-wide hypomethylation, dysregulated genome expression, and seedling lethality

Cytosine methylation at CG sites ((m)CG) plays critical roles in development, epigenetic inheritance, and genome stability in mammals and plants. In the dicot model plant Arabidopsis thaliana, methyltransferase 1 (MET1), a principal CG methylase, functions to maintain (m)CG during DNA replication, with its null mutation resulting in global hypomethylation and pleiotropic developmental defects. Null mutation of a critical CG methylase has not been characterized at a whole-genome level in other higher eukaryotes, leaving the generality of the Arabidopsis findings largely speculative. Rice is a model plant of monocots, to which many of our important crops belong. Here we have characterized a null mutant of OsMet1-2, the major CG methylase in rice. We found that seeds homozygous for OsMet1-2 gene mutation (OsMET1-2(-/-)), which directly segregated from normal heterozygote plants (OsMET1-2(+/-)), were seriously maldeveloped, and all germinated seedlings underwent swift necrotic death. Compared with wild type, genome-wide loss of (m)CG occurred in the mutant methylome, which was accompanied by a plethora of quantitative molecular phenotypes including dysregulated expression of diverse protein-coding genes, activation and repression of transposable elements, and altered small RNA profiles. Our results have revealed conservation but also distinct functional differences in CG methylases between rice and Arabidopsis.

Epigenetic control of plant senescence and linked processes

Senescence processes are part of the plant developmental programme. They involve reprogramming of gene expression and are under the control of a complex regulatory network closely linked to other developmental and stress-responsive pathways. Recent evidence indicates that leaf senescence is regulated via epigenetic mechanisms. In the present review, the epigenetic control of plant senescence is discussed in the broader context of environment-sensitive plant development. The review outlines the concept of epigenetic control of interconnected regulatory pathways steering stress responses and plant development. Besides giving an overview of techniques used in the field, it summarizes recent findings on global alterations in chromatin structure, histone and DNA modifications, and ATP-dependent chromatin remodelling during plant senescence and linked processes.

Transposon variation by order during allopolyploidisation between Brassica oleracea and Brassica rapa

Although many studies have shown that transposable element (TE) activation is induced by hybridisation and polyploidisation in plants, much less is known on how different types of TE respond to hybridisation, and the impact of TE-associated sequences on gene function. We investigated the frequency and regularity of putative transposon activation for different types of TE, and determined the impact of TE-associated sequence variation on the genome during allopolyploidisation. We designed different types of TE primers and adopted the Inter-Retrotransposon Amplified Polymorphism (IRAP) method to detect variation in TE-associated sequences during the process of allopolyploidisation between Brassica rapa (AA) and Brassica oleracea (CC), and in successive generations of self-pollinated progeny. In addition, fragments with TE insertions were used to perform Blast2GO analysis to characterise the putative functions of the fragments with TE insertions. Ninety-two primers amplifying 548 loci were used to detect variation in sequences associated with four different orders of TE sequences. TEs could be classed in ascending frequency into LTR-REs, TIRs, LINEs, SINEs and unknown TEs. The frequency of novel variation (putative activation) detected for the four orders of TEs was highest from the F1 to F2 generations, and lowest from the F2 to F3 generations. Functional annotation of sequences with TE insertions showed that genes with TE insertions were mainly involved in metabolic processes and binding, and preferentially functioned in organelles. TE variation in our study severely disturbed the genetic compositions of the different generations, resulting in inconsistencies in genetic clustering. Different types of TE showed different patterns of variation during the process of allopolyploidisation.

Efficient screening of long terminal repeat retrotransposons that show high insertion polymorphism via high-throughput sequencing of the primer binding site

Retrotransposons have been used frequently for the development of molecular markers by using their insertion polymorphisms among cultivars, because multiple copies of these elements are dispersed throughout the genome and inserted copies are inherited genetically. Although a large number of long terminal repeat (LTR) retrotransposon families exist in the higher eukaryotic genomes, the identification of families that show high insertion polymorphism has been challenging. Here, we performed an efficient screening of these retrotransposon families using an Illumina HiSeq2000 sequencing platform with comprehensive LTR library construction based on the primer binding site (PBS), which is located adjacent to the 5' LTR and has a motif that is universal and conserved among LTR retrotransposon families. The paired-end sequencing library of the fragments containing a large number of LTR sequences and their insertion sites was sequenced for seven strawberry (Fragaria × ananassa Duchesne) cultivars and one diploid wild species (Fragaria vesca L.). Among them, we screened 24 families with a "unique" insertion site that appeared only in one cultivar and not in any others, assuming that this type of insertion should have occurred quite recently. Finally, we confirmed experimentally the selected LTR families showed high insertion polymorphisms among closely related cultivars.

Spy: a new group of eukaryotic DNA transposons without target site duplications

Class 2 or DNA transposons populate the genomes of most eukaryotes and like other mobile genetic elements have a profound impact on genome evolution. Most DNA transposons belong to the cut-and-paste types, which are relatively simple elements characterized by terminal-inverted repeats (TIRs) flanking a single gene encoding a transposase. All eukaryotic cut-and-paste transposons so far described are also characterized by target site duplications (TSDs) of host DNA generated upon chromosomal insertion. Here, we report a new group of evolutionarily related DNA transposons called Spy, which also include TIRs and DDE motif-containing transposase but surprisingly do not create TSDs upon insertion. Instead, Spy transposons appear to transpose precisely between 5'-AAA and TTT-3' host nucleotides, without duplication or modification of the AAATTT target sites. Spy transposons were identified in the genomes of diverse invertebrate species based on transposase homology searches and structure-based approaches. Phylogenetic analyses indicate that Spy transposases are distantly related to IS5, ISL2EU, and PIF/Harbinger transposases. However, Spy transposons are distinct from these and other DNA transposon superfamilies by their lack of TSD and their target site preference. Our findings expand the known diversity of DNA transposons and reveal a new group of eukaryotic DDE transposases with unusual catalytic properties.

BrassicaTED - a public database for utilization of miniature transposable elements in Brassica species

Background

MITE, TRIM and SINEs are miniature form transposable elements (mTEs) that are ubiquitous and dispersed throughout entire plant genomes. Tens of thousands of members cause insertion polymorphism at both the inter- and intra- species level. Therefore, mTEs are valuable targets and resources for development of markers that can be utilized for breeding, genetic diversity and genome evolution studies. Taking advantage of the completely sequenced genomes of Brassica rapa and B. oleracea, characterization of mTEs and building a curated database are prerequisite to extending their utilization for genomics and applied fields in Brassica crops.

Findings

We have developed BrassicaTED as a unique web portal containing detailed characterization information for mTEs of Brassica species. At present, BrassicaTED has datasets for 41 mTE families, including 5894 and 6026 members from 20 MITE families, 1393 and 1639 members from 5 TRIM families, 1270 and 2364 members from 16 SINE families in B. rapa and B. oleracea, respectively. BrassicaTED offers different sections to browse structural and positional characteristics for every mTE family. In addition, we have added data on 289 MITE insertion polymorphisms from a survey of seven Brassica relatives. Genes with internal mTE insertions are shown with detailed gene annotation and microarray-based comparative gene expression data in comparison with their paralogs in the triplicated B. rapa genome. This database also includes a novel tool, K BLAST (Karyotype BLAST), for clear visualization of the locations for each member in the B. rapa and B. oleracea pseudo-genome sequences.

Conclusions

BrassicaTED is a newly developed database of information regarding the characteristics and potential utility of mTEs including MITE, TRIM and SINEs in B. rapa and B. oleracea. The database will promote the development of desirable mTE-based markers, which can be utilized for genomics and breeding in Brassica species. BrassicaTED will be a valuable repository for scientists and breeders, promoting efficient research on Brassica species. BrassicaTED can be accessed at http://im-crop.snu.ac.kr/BrassicaTED/index.php.

A rapid and enhanced DNA detection method for crop cultivar discrimination

In many crops species, the development of a rapid and precise cultivar discrimination system has been required for plant breeding and patent protection of plant cultivars and agricultural products. Here, we successfully evaluated strawberry cultivars via a novel method, namely, the single tag hybridization (STH) chromatographic printed array strip (PAS) using the PCR products of eight genomic regions. In a previous study, we showed that genotyping of eight genomic regions derived from FaRE1 retrotransposon insertion site enabled to discriminate 32 strawberry cultivars precisely, however, this method required agarose/acrylamide gel electrophoresis, thus has the difficulty for practical application. In contrast, novel DNA detection method in this study has some great advantages over standard DNA detection methods, including agarose/acrylamide gel electrophoresis, because it produces signals for DNA detection with dramatically higher sensitivity in a shorter time without any preparation or staining of a gel. Moreover, this method enables the visualization of multiplex signals simultaneously in a single reaction using several independent amplification products. We expect that this novel method will become a rapid and convenient cultivar screening assay for practical purposes, and will be widely applied to various situations, including laboratory research, and on-site inspection of plant cultivars and agricultural products.

Efficient DNA fingerprinting based on the targeted sequencing of active retrotransposon insertion sites using a bench-top high-throughput sequencing platform

In many crop species, DNA fingerprinting is required for the precise identification of cultivars to protect the rights of breeders. Many families of retrotransposons have multiple copies throughout the eukaryotic genome and their integrated copies are inherited genetically. Thus, their insertion polymorphisms among cultivars are useful for DNA fingerprinting. In this study, we conducted a DNA fingerprinting based on the insertion polymorphisms of active retrotransposon families (Rtsp-1 and LIb) in sweet potato. Using 38 cultivars, we identified 2,024 insertion sites in the two families with an Illumina MiSeq sequencing platform. Of these insertion sites, 91.4% appeared to be polymorphic among the cultivars and 376 cultivar-specific insertion sites were identified, which were converted directly into cultivar-specific sequence-characterized amplified region (SCAR) markers. A phylogenetic tree was constructed using these insertion sites, which corresponded well with known pedigree information, thereby indicating their suitability for genetic diversity studies. Thus, the genome-wide comparative analysis of active retrotransposon insertion sites using the bench-top MiSeq sequencing platform is highly effective for DNA fingerprinting without any requirement for whole genome sequence information. This approach may facilitate the development of practical polymerase chain reaction-based cultivar diagnostic system and could also be applied to the determination of genetic relationships.

mPing: The bursting transposon

Though transposable elements (TEs) have been considered as an efficient source of evolution, it has never been possible to test this hypothesis because most of TE insertions had occurred millions of years ago, or because currently active TEs have very few copies in a host genome. However, mPing, the first active DNA transposon in rice, was revealed to hold a key to answer this question. mPing has attained high copy numbers and still retained very high activity in a traditional rice strain, which enabled direct observation of behavior and impact of a bursting TE. A comprehensive analysis of mPing insertion sites has revealed it avoids exons but prefers promoter regions and thus moderately affects transcription of neighboring genes. Some of the mPing insertions have introduced possibly useful expression profile to adjacent genes that indicated TE's potential in de novo formation of gene regulatory network.

Recurrent horizontal transfers of Chapaev transposons in diverse invertebrate and vertebrate animals

Horizontal transfer (HT) of a transposable element (TE) into a new genome is regarded as an important force to drive genome variation and biological innovation. In addition, HT also plays an important role in the persistence of TEs in eukaryotic genomes. Here, we provide the first documented example for the repeated HT of three families of Chapaev transposons in a wide range of animal species, including mammals, reptiles, jawed fishes, lampreys, insects, and in an insect bracovirus. Multiple alignments of the Chapaev transposons identified in these species revealed extremely high levels of nucleotide sequence identity (79-99%), which are inconsistent with vertical evolution given the deep divergence time separating these host species. Rather, the discontinuous distribution amongst species and lack of purifying selection acting on these transposons strongly suggest that they were independently and horizontally transferred into these species lineages. The detection of Chapaev transposons in an insect bracovirus indicated that these viruses might act as a possible vector for the horizontal spread of Chapaev transposons. One of the Chapaev families was also shared by lampreys and some of their common hosts (such as sturgeon and paddlefish), which suggested that parasite-host interaction might facilitate HTs.

cDNA Microarray Analysis of Rice Anther Genes under Chilling Stress at the Microsporogenesis Stage Revealed Two Genes with DNA TransposonCastawayin the 5′-Flanking Region

Rice is most chilling sensitive at the onset of microspore release. Chilling treatment at this stage causes male sterility. The gene expression profile during the microspore development process under chilling stress was revealed using a microarray that included 8,987 rice cDNAs. As many as 160 cDNAs were up- or down-regulated by chilling during the microspore release stage. RT-PCR analysis of 5 genes confirmed the microarray results. We identified 3 novel genes whose expression levels were remarkably changed by chilling in rice anther. A new cis element that includes a DNA transposon Castaway sequence was found in the 5' upstream region of two genes which were conspicuously down-regulated by chilling temperatures in rice anther.

Analysis of transposable elements in the genome of Asparagus officinalis from high coverage sequence data

Asparagus officinalis is an economically and nutritionally important vegetable crop that is widely cultivated and is used as a model dioecious species to study plant sex determination and sex chromosome evolution. To improve our understanding of its genome composition, especially with respect to transposable elements (TEs), which make up the majority of the genome, we performed Illumina HiSeq2000 sequencing of both male and female asparagus genomes followed by bioinformatics analysis. We generated 17 Gb of sequence (12×coverage) and assembled them into 163,406 scaffolds with a total cumulated length of 400 Mbp, which represent about 30% of asparagus genome. Overall, TEs masked about 53% of the A. officinalis assembly. Majority of the identified TEs belonged to LTR retrotransposons, which constitute about 28% of genomic DNA, with Ty1/copia elements being more diverse and accumulated to higher copy numbers than Ty3/gypsy. Compared with LTR retrotransposons, non-LTR retrotransposons and DNA transposons were relatively rare. In addition, comparison of the abundance of the TE groups between male and female genomes showed that the overall TE composition was highly similar, with only slight differences in the abundance of several TE groups, which is consistent with the relatively recent origin of asparagus sex chromosomes. This study greatly improves our knowledge of the repetitive sequence construction of asparagus, which facilitates the identification of TEs responsible for the early evolution of plant sex chromosomes and is helpful for further studies on this dioecious plant.

Environmental stress activation of plant long-terminal repeat retrotransposons

Genomic retrotransposons (RTs) are major components of most plant genomes. They spread throughout the genomes by a process termed retrotransposition, which consists of reverse transcription and reinsertion of the copied element into a new genomic location (a copy-and-paste system). Abiotic and biotic stresses activate long-terminal repeat (LTR) RTs in photosynthetic eukaryotes from algae to angiosperms. LTR RTs could represent a threat to the integrity of host genomes because of their activity and mutagenic potential by epigenetic regulation. Host genomes have developed mechanisms to control the activity of the retroelements and their mutagenic potential. Some LTR RTs escape these defense mechanisms, and maintain their ability to be activated and transpose as a result of biotic or abiotic stress stimuli. These stimuli include pathogen infection, mechanical damage, in vitro tissue culturing, heat, drought and salt stress, generation of doubled haploids, X-ray irradiation and many others. Reactivation of LTR RTs differs between different plant genomes. The expression levels of reactivated RTs are influenced by the transcriptional and post-transcriptional gene silencing mechanisms (e.g. DNA methylation, heterochromatin formation and RNA interference). Moreover, the insertion of RTs (e.g. Triticum aestivum L. Wis2-1A) into or next to coding regions of the host genome can generate changes in the expression of adjacent host genes of the host. In this paper, we review the ways that plant genomic LTR RTs are activated by environmental stimuli to affect restructuring and diversification of the host genome.

Genomic localization of AtRE1 and AtRE2, copia-type retrotransposons, in natural variants of Arabidopsis thaliana

Retrotransposons are ubiquitous components of plant genomes. They affect genome organization, and can also affect the expression patterns of neighboring genes. Retrotransposons are therefore important elements for changing genomic information. To understand the evolution of the Arabidopsis genome, we examined the distribution of certain retrotransposons, AtRE1s and AtRE2s, in the genomes of 12 natural variants (accessions) of Arabidopsis thaliana. AtRE1 and AtRE2 are copia-type retrotransposons that are potentially active. Their copy numbers are low, and they are absent from the genomes of some accessions. We detected four loci with AtRE1s inserted in six accessions, and one locus with an insertion of a solo-LTR-like sequence derived from AtRE1 in two accessions. Seven loci with AtRE2s inserted were detected on eight accessions. These loci were distributed in euchromatic regions of chromosomes 1, 2, 3, and 4. The AtRE1 and AtRE2 sequences at some loci identified in this study have not been recorded in the database of the 1001 Genome project. The sequences of AtRE1s and those of AtRE2s in different accessions and at different loci were highly conserved. There was a complete or almost complete conservation of sequences of both long terminal repeats in each AtRE1 and in each AtRE2. These results suggest that AtRE1 and AtRE2 appeared quite recently in the Arabidopsis genome. Furthermore, sequence comparisons of AtRE1 and AtRE2 loci among accessions revealed the possibility that large deletions containing entire sequences of AtRE1 and AtRE2 have occurred in some accessions.

Genome-wide comparative analysis of 20 miniature inverted-repeat transposable element families in Brassica rapa and B. oleracea

Miniature inverted-repeat transposable elements (MITEs) are ubiquitous, non-autonomous class II transposable elements. Here, we conducted genome-wide comparative analysis of 20 MITE families in B. rapa, B. oleracea, and Arabidopsis thaliana. A total of 5894 and 6026 MITE members belonging to the 20 families were found in the whole genome pseudo-chromosome sequences of B. rapa and B. oleracea, respectively. Meanwhile, only four of the 20 families, comprising 573 members, were identified in the Arabidopsis genome, indicating that most of the families were activated in the Brassica genus after divergence from Arabidopsis. Copy numbers varied from 4 to 1459 for each MITE family, and there was up to 6-fold variation between B. rapa and B. oleracea. In particular, analysis of intact members showed that whereas eleven families were present in similar copy numbers in B. rapa and B. oleracea, nine families showed copy number variation ranging from 2- to 16-fold. Four of those families (BraSto-3, BraTo-3, 4, 5) were more abundant in B. rapa, and the other five (BraSto-1, BraSto-4, BraTo-1, 7 and BraHAT-1) were more abundant in B. oleracea. Overall, 54% and 51% of the MITEs resided in or within 2 kb of a gene in the B. rapa and B. oleracea genomes, respectively. Notably, 92 MITEs were found within the CDS of annotated genes, suggesting that MITEs might play roles in diversification of genes in the recently triplicated Brassica genome. MITE insertion polymorphism (MIP) analysis of 289 MITE members showed that 52% and 23% were polymorphic at the inter- and intra-species levels, respectively, indicating that there has been recent MITE activity in the Brassica genome. These recently activated MITE families with abundant MIP will provide useful resources for molecular breeding and identification of novel functional genes arising from MITE insertion.

Impact of the competition between mating types on the cultivation of Tuber melanosporum: Romeo and Juliet and the matter of space and time

Major breakthroughs in our understanding of the life cycles of the symbiotic ascomycetes belonging to the genus Tuber have occurred over the last several years. A number of Tuber species produce edible fruiting bodies, known as truffles, that are marketed worldwide. A better understanding of the basic biological characteristics of Tuber spp. is likely to have tremendous practical relevance for their cultivation. Tuber melanosporum produces the most valuable black truffles and its genome has been recently sequenced. This species is now serving as a model for studying the biology of truffles. Here, we review recent progress in the understanding of sexual reproduction modalities in T. melanosporum. The practical relevance of these findings is outlined. In particular, the discoveries that T. melanosporum is heterothallic and that strains of different mating types compete to persist on the roots of host plants suggest that the spatial and temporal distributional patterns of strains of different mating types are key determinants of truffle fructification. The spatial segregation of the two mating types in areas where T. melanosporum occurs likely limits truffle production. Thus, host plant inoculation techniques and agronomic practices that might be pursued to manage T. melanosporum orchards with a balanced presence of the two mating partners are described.

The precarious prokaryotic chromosome

Evolutionary selection for optimal genome preservation, replication, and expression should yield similar chromosome organizations in any type of cells. And yet, the chromosome organization is surprisingly different between eukaryotes and prokaryotes. The nuclear versus cytoplasmic accommodation of genetic material accounts for the distinct eukaryotic and prokaryotic modes of genome evolution, but it falls short of explaining the differences in the chromosome organization. I propose that the two distinct ways to organize chromosomes are driven by the differences between the global-consecutive chromosome cycle of eukaryotes and the local-concurrent chromosome cycle of prokaryotes. Specifically, progressive chromosome segregation in prokaryotes demands a single duplicon per chromosome, while other "precarious" features of the prokaryotic chromosomes can be viewed as compensations for this severe restriction.

Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species

Hot pepper (Capsicum annuum), one of the oldest domesticated crops in the Americas, is the most widely grown spice crop in the world. We report whole-genome sequencing and assembly of the hot pepper (Mexican landrace of Capsicum annuum cv. CM334) at 186.6× coverage. We also report resequencing of two cultivated peppers and de novo sequencing of the wild species Capsicum chinense. The genome size of the hot pepper was approximately fourfold larger than that of its close relative tomato, and the genome showed an accumulation of Gypsy and Caulimoviridae family elements. Integrative genomic and transcriptomic analyses suggested that change in gene expression and neofunctionalization of capsaicin synthase have shaped capsaicinoid biosynthesis. We found differential molecular patterns of ripening regulators and ethylene synthesis in hot pepper and tomato. The reference genome will serve as a platform for improving the nutritional and medicinal values of Capsicum species.

Major repeat components covering one-third of the ginseng (Panax ginseng C.A. Meyer) genome and evidence for allotetraploidy

Ginseng (Panax ginseng) is a famous medicinal herb, but the composition and structure of its genome are largely unknown. Here we characterized the major repeat components and inspected their distribution in the ginseng genome. By analyzing three repeat-rich bacterial artificial chromosome (BAC) sequences from ginseng, we identified complex insertion patterns of 34 long terminal repeat retrotransposons (LTR-RTs) and 11 LTR-RT derivatives accounting for more than 80% of the BAC sequences. The LTR-RTs were classified into three Ty3/gypsy (PgDel, PgTat and PgAthila) and two Ty1/Copia (PgTork and PgOryco) families. Mapping of 30-Gbp Illumina whole-genome shotgun reads to the BAC sequences revealed that these five LTR-RT families occupy at least 34% of the ginseng genome. The Ty3/Gypsy families were predominant, comprising 74 and 33% of the BAC sequences and the genome, respectively. In particular, the PgDel family accounted for 29% of the genome and presumably played major roles in enlargement of the size of the ginseng genome. Fluorescence in situ hybridization (FISH) revealed that the PgDel1 elements are distributed throughout the chromosomes along dispersed heterochromatic regions except for ribosomal DNA blocks. The intensity of the PgDel2 FISH signals was biased toward 24 out of 48 chromosomes. Unique gene probes showed two pairs of signals with different locations, one pair in subtelomeric regions on PgDel2-rich chromosomes and the other in interstitial regions on PgDel2-poor chromosomes, demonstrating allotetraploidy in ginseng. Our findings promote understanding of the evolution of the ginseng genome and of that of related species in the Araliaceae.

Dicer-like 3 produces transposable element-associated 24-nt siRNAs that control agricultural traits in rice

Transposable elements (TEs) and repetitive sequences make up over 35% of the rice (Oryza sativa) genome. The host regulates the activity of different TEs by different epigenetic mechanisms, including DNA methylation, histone H3K9 methylation, and histone H3K4 demethylation. TEs can also affect the expression of host genes. For example, miniature inverted repeat TEs (MITEs), dispersed high copy-number DNA TEs, can influence the expression of nearby genes. In plants, 24-nt small interfering RNAs (siRNAs) are mainly derived from repeats and TEs. However, the extent to which TEs, particularly MITEs associated with 24-nt siRNAs, affect gene expression remains elusive. Here, we show that the rice Dicer-like 3 homolog OsDCL3a is primarily responsible for 24-nt siRNA processing. Impairing OsDCL3a expression by RNA interference caused phenotypes affecting important agricultural traits; these phenotypes include dwarfism, larger flag leaf angle, and fewer secondary branches. We used small RNA deep sequencing to identify 535,054 24-nt siRNA clusters. Of these clusters, ∼82% were OsDCL3a-dependent and showed significant enrichment of MITEs. Reduction of OsDCL3a function reduced the 24-nt siRNAs predominantly from MITEs and elevated expression of nearby genes. OsDCL3a directly targets genes involved in gibberellin and brassinosteroid homeostasis; OsDCL3a deficiency may affect these genes, thus causing the phenotypes of dwarfism and enlarged flag leaf angle. Our work identifies OsDCL3a-dependent 24-nt siRNAs derived from MITEs as broadly functioning regulators for fine-tuning gene expression, which may reflect a conserved epigenetic mechanism in higher plants with genomes rich in dispersed repeats or TEs.

MAKER-P: a tool kit for the rapid creation, management, and quality control of plant genome annotations

We have optimized and extended the widely used annotation engine MAKER in order to better support plant genome annotation efforts. New features include better parallelization for large repeat-rich plant genomes, noncoding RNA annotation capabilities, and support for pseudogene identification. We have benchmarked the resulting software tool kit, MAKER-P, using the Arabidopsis (Arabidopsis thaliana) and maize (Zea mays) genomes. Here, we demonstrate the ability of the MAKER-P tool kit to automatically update, extend, and revise the Arabidopsis annotations in light of newly available data and to annotate pseudogenes and noncoding RNAs absent from The Arabidopsis Informatics Resource 10 build. Our results demonstrate that MAKER-P can be used to manage and improve the annotations of even Arabidopsis, perhaps the best-annotated plant genome. We have also installed and benchmarked MAKER-P on the Texas Advanced Computing Center. We show that this public resource can de novo annotate the entire Arabidopsis and maize genomes in less than 3 h and produce annotations of comparable quality to those of the current The Arabidopsis Information Resource 10 and maize V2 annotation builds.

Analysis of T-DNA/Host-Plant DNA Junction Sequences in Single-Copy Transgenic Barley Lines

Sequencing across the junction between an integrated transfer DNA (T-DNA) and a host plant genome provides two important pieces of information. The junctions themselves provide information regarding the proportion of T-DNA which has integrated into the host plant genome, whilst the transgene flanking sequences can be used to study the local genetic environment of the integrated transgene. In addition, this information is important in the safety assessment of GM crops and essential for GM traceability. In this study, a detailed analysis was carried out on the right-border T-DNA junction sequences of single-copy independent transgenic barley lines. T-DNA truncations at the right-border were found to be relatively common and affected 33.3% of the lines. In addition, 14.3% of lines had rearranged construct sequence after the right border break-point. An in depth analysis of the host-plant flanking sequences revealed that a significant proportion of the T-DNAs integrated into or close to known repetitive elements. However, this integration into repetitive DNA did not have a negative effect on transgene expression.

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.

Annotation and sequence diversity of transposable elements in common bean (Phaseolus vulgaris)

Common bean (Phaseolus vulgaris) is an important legume crop grown and consumed worldwide. With the availability of the common bean genome sequence, the next challenge is to annotate the genome and characterize functional DNA elements. Transposable elements (TEs) are the most abundant component of plant genomes and can dramatically affect genome evolution and genetic variation. Thus, it is pivotal to identify TEs in the common bean genome. In this study, we performed a genome-wide transposon annotation in common bean using a combination of homology and sequence structure-based methods. We developed a 2.12-Mb transposon database which includes 791 representative transposon sequences and is available upon request or from www.phytozome.org. Of note, nearly all transposons in the database are previously unrecognized TEs. More than 5,000 transposon-related expressed sequence tags (ESTs) were detected which indicates that some transposons may be transcriptionally active. Two Ty1-copia retrotransposon families were found to encode the envelope-like protein which has rarely been identified in plant genomes. Also, we identified an extra open reading frame (ORF) termed ORF2 from 15 Ty3-gypsy families that was located between the ORF encoding the retrotransposase and the 3'LTR. The ORF2 was in opposite transcriptional orientation to retrotransposase. Sequence homology searches and phylogenetic analysis suggested that the ORF2 may have an ancient origin, but its function is not clear. These transposon data provide a useful resource for understanding the genome organization and evolution and may be used to identify active TEs for developing transposon-tagging system in common bean and other related genomes.

Comparative genomic paleontology across plant kingdom reveals the dynamics of TE-driven genome evolution

Long terminal repeat-retrotransposons (LTR-RTs) are the most abundant class of transposable elements (TEs) in plants. They strongly impact the structure, function, and evolution of their host genome, and, in particular, their role in genome size variation has been clearly established. However, the dynamics of the process through which LTR-RTs have differentially shaped plant genomes is still poorly understood because of a lack of comparative studies. Using a new robust and automated family classification procedure, we exhaustively characterized the LTR-RTs in eight plant genomes for which a high-quality sequence is available (i.e., Arabidopsis thaliana, A. lyrata, grapevine, soybean, rice, Brachypodium dystachion, sorghum, and maize). This allowed us to perform a comparative genome-wide study of the retrotranspositional landscape in these eight plant lineages from both monocots and dicots. We show that retrotransposition has recurrently occurred in all plant genomes investigated, regardless their size, and through bursts, rather than a continuous process. Moreover, in each genome, only one or few LTR-RT families have been active in the recent past, and the difference in genome size among the species studied could thus mostly be accounted for by the extent of the latest transpositional burst(s). Following these bursts, LTR-RTs are efficiently eliminated from their host genomes through recombination and deletion, but we show that the removal rate is not lineage specific. These new findings lead us to propose a new model of TE-driven genome evolution in plants.

Unintended compositional changes in genetically modified (GM) crops: 20 years of research

The compositional equivalency between genetically modified (GM) crops and nontransgenic comparators has been a fundamental component of human health safety assessment for 20 years. During this time, a large amount of information has been amassed on the compositional changes that accompany both the transgenesis process and traditional breeding methods; additionally, the genetic mechanisms behind these changes have been elucidated. After two decades, scientists are encouraged to objectively assess this body of literature and determine if sufficient scientific uncertainty still exists to continue the general requirement for these studies to support the safety assessment of transgenic crops. It is concluded that suspect unintended compositional effects that could be caused by genetic modification have not materialized on the basis of this substantial literature. Hence, compositional equivalence studies uniquely required for GM crops may no longer be justified on the basis of scientific uncertainty.

Genetics and consequences of crop domestication

Phenotypic variation has been manipulated by humans during crop domestication, which occurred primarily between 3000 and 10000 years ago in the various centers of origin around the world. The process of domestication has profound consequences on crops, where the domesticate has moderately reduced genetic diversity relative to the wild ancestor across the genome, and severely reduced diversity for genes targeted by domestication. The question that remains is whether reduction in genetic diversity has affected crop production today. A case study in maize ( Zea mays ) demonstrates the application of understanding relationships between genetic diversity and phenotypic diversity in the wild ancestor and the domesticate. As an outcrossing species, maize has tremendous genetic variation. The complementary combination of genome-wide association mapping (GWAS) approaches, large HapMap data sets, and germplasm resources is leading to important discoveries of the relationship between genetic diversity and phenotypic variation and the impact of domestication on trait variation.

Isolation and characterization of novel Ty1-copia-like retrotransposons from lily

Species of the genus Lilium are well known for their large genomes. Although expansion of noncoding repeated DNA is believed to account for this genome size, retroelement del Ty3-gypsy is the only one described so far in the genus Lilium. We isolated Ty1-copia elements from Lilium longiflorum and named them LIREs (lily retrotransposons). The long terminal repeats, primer binding site, and polypurine tract sequences are highly similar among the LIRE elements, indicating that they are in the same lineage. Although the protein-coding regions were highly decayed, the sequence motifs of the integrase, reverse transcriptase, and RNase H domains were identifiable as belonging to the order of Ty1-copia elements. Phylogenetic analysis and primer binding site sequences revealed that these elements belonged to the Ale lineage among the six lineages of plant Ty1-copia elements. Base substitutions in the long terminal repeats estimated that the integration times of the LIRE Ty1-copia elements were between 0.7 and 5.5 mya. In situ hybridization showed that the LIRE elements were present in all the chromosomes of L. longiflorum and L. lancifolium, but absent in centromeres, telomeres, and 45S rRNA sites in both species. The LIRE elements were present very abundantly in species of the genus Lilium, but absent in other genera of the family Liliaceae, implying that the LIRE elements might have contributed to the expansion of the genome in the genus Lilium.

Mobile DNA distributions refine the phylogeny of "matsutake" mushrooms, Tricholoma sect. Caligata

"Matsutake" mushrooms are formed by several species of Tricholoma sect. Caligata distributed across the northern hemisphere. A phylogenetic analysis of matsutake based on virtually neutral mutations in DNA sequences resolved robust relationships among Tricholoma anatolicum, Tricholoma bakamatsutake, Tricholoma magnivelare, Tricholoma matsutake, and Tricholoma sp. from Mexico (=Tricholoma sp. Mex). However, relationships among these matsutake and other species, such as Tricholoma caligatum and Tricholoma fulvocastaneum, were ambiguous. We, therefore, analyzed genomic copy numbers of σ marY1 , marY1, and marY2N retrotransposons by comparing them with the single-copy mobile DNA megB1 using real-time polymerase chain reaction (PCR) to clarify matsutake phylogeny. We also examined types of megB1-associated domains, composed of a number of poly (A) and poly (T) reminiscent of RNA-derived DNA elements among these species. Both datasets resolved two distinct groups, one composed of T. bakamatsutake, T. fulvocastaneum, and T. caligatum that could have diverged earlier and the other comprising T. magnivelare, Tricholoma sp. Mex, T. anatolicum, and T. matsutake that could have evolved later. In the first group, T. caligatum was the closest to the second group, followed by T. fulvocastaneum and T. bakamatsutake. Within the second group, T. magnivelare was clearly differentiated from the other species. The data suggest that matsutake underwent substantial evolution between the first group, mostly composed of Fagaceae symbionts, and the second group, comprised only of Pinaceae symbionts, but diverged little within each groups. Mobile DNA markers could be useful in resolving difficult phylogenies due to, for example, closely spaced speciation events.

Lessons from natural and artificial polyploids in higher plants

Polyploidy in higher plants is a major source of genetic novelty upon which selection may act to drive evolution, as evidenced by the widespread success of polyploid species in the wild. However, research into the effects of polyploidy can be confounded by the entanglement of several processes: genome duplication, hybridisation (allopolyploidy is frequent in plants) and subsequent evolution. The discovery of the chemical agent colchicine, which can be used to produce artificial polyploids on demand, has enabled scientists to unravel these threads and understand the complex genomic changes involved in each. We present here an overview of lessons learnt from studies of natural and artificial polyploids, and from comparisons between the 2, covering basic cellular and metabolic consequences through to alterations in epigenetic gene regulation, together with 2 in-depth case studies in Senecio and Glycine. See also the sister article focusing on animals by Arai and Fujimoto in this themed issue.

Identification of genetic elements associated with EPSPs gene amplification

Weed populations can have high genetic plasticity and rapid responses to environmental selection pressures. For example, 100-fold amplification of the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene evolved in the weed species Amaranthus palmeri to confer resistance to glyphosate, the world's most important herbicide. However, the gene amplification mechanism is unknown. We sequenced the EPSPS gene and genomic regions flanking EPSPS loci in A. palmeri, and searched for mobile genetic elements or repetitive sequences. The EPSPS gene was 10,229 bp, containing 8 exons and 7 introns. The gene amplification likely proceeded through a DNA-mediated mechanism, as introns exist in the amplified gene copies and the entire amplified sequence is at least 30 kb in length. Our data support the presence of two EPSPS loci in susceptible (S) A. palmeri, and that only one of these was amplified in glyphosate-resistant (R) A. palmeri. The EPSPS gene amplification event likely occurred recently, as no sequence polymorphisms were found within introns of amplified EPSPS copies from R individuals. Sequences with homology to miniature inverted-repeat transposable elements (MITEs) were identified next to EPSPS gene copies only in R individuals. Additionally, a putative Activator (Ac) transposase and a repetitive sequence region were associated with amplified EPSPS genes. The mechanism controlling this DNA-mediated amplification remains unknown. Further investigation is necessary to determine if the gene amplification may have proceeded via DNA transposon-mediated replication, and/or unequal recombination between different genomic regions resulting in replication of the EPSPS gene.

Genome-wide prediction of nucleosome occupancy in maize reveals plant chromatin structural features at genes and other elements at multiple scales

The nucleosome is a fundamental structural and functional chromatin unit that affects nearly all DNA-templated events in eukaryotic genomes. It is also a biochemical substrate for higher order, cis-acting gene expression codes and the monomeric structural unit for chromatin packaging at multiple scales. To predict the nucleosome landscape of a model plant genome, we used a support vector machine computational algorithm trained on human chromatin to predict the nucleosome occupancy likelihood (NOL) across the maize (Zea mays) genome. Experimentally validated NOL plots provide a novel genomic annotation that highlights gene structures, repetitive elements, and chromosome-scale domains likely to reflect regional gene density. We established a new genome browser (http://www.genomaize.org) for viewing support vector machine-based NOL scores. This annotation provides sequence-based comprehensive coverage across the entire genome, including repetitive genomic regions typically excluded from experimental genomics data. We find that transposable elements often displayed family-specific NOL profiles that included distinct regions, especially near their termini, predicted to have strong affinities for nucleosomes. We examined transcription start site consensus NOL plots for maize gene sets and discovered that most maize genes display a typical +1 nucleosome positioning signal just downstream of the start site but not upstream. This overall lack of a -1 nucleosome positioning signal was also predicted by our method for Arabidopsis (Arabidopsis thaliana) genes and verified by additional analysis of previously published Arabidopsis MNase-Seq data, revealing a general feature of plant promoters. Our study advances plant chromatin research by defining the potential contribution of the DNA sequence to observed nucleosome positioning and provides an invariant baseline annotation against which other genomic data can be compared.

Loss-of-function of a ubiquitin-related modifier promotes the mobilization of the active MITE mPing

Miniature inverted-repeat transposable elements (MITEs) are widespread in both prokaryotic and eukaryotic genomes, where their copy numbers can attain several thousands. Little is known, however, about the genetic factor(s) affecting their transpositions. Here, we show that disruption of a gene encoding ubiquitin-like protein markedly enhances the transposition activity of a MITE mPing in intact rice plants without any exogenous stresses. We found that the transposition activity of mPing is far higher in the lines harboring a non-functional allele at the Rurm1 (Rice ubiquitin-related modifier-1) locus than in the wild-type line. Although the alteration of cytosine methylation pattern triggers the activation of transposable elements under exogenous stress conditions, the methylation degrees in the whole genome, the mPing-body region, and the mPing-flanking regions of the non-functional Rurm1 line were unchanged. This study provides experimental evidence for one of the models of genome shock theory that genetic accidents within cells enhance the transposition activities of transposable elements.

Rice epigenomics and epigenetics: challenges and opportunities

During recent years rice genome-wide epigenomic information such as DNA methylation and histone modifications, which are important for genome activity has been accumulated. The function of a number of rice epigenetic regulators has been studied, many of which are found to be involved in a diverse range of developmental and stress-responsive pathways. Analysis of epigenetic variations among different rice varieties indicates that epigenetic modification may lead to inheritable phenotypic variation. Characterizing phenotypic consequences of rice epigenomic variations and the underlining chromatin mechanism and identifying epialleles related to important agronomic traits may provide novel strategies to enhance agronomically favorable traits and grain productivity in rice.

Metagenomic analysis of healthy and white plague-affected Mussismilia braziliensis corals

Coral health is under threat throughout the world due to regional and global stressors. White plague disease (WP) is one of the most important threats affecting the major reef builder of the Abrolhos Bank in Brazil, the endemic coral Mussismilia braziliensis. We performed a metagenomic analysis of healthy and WP-affected M. braziliensis in order to determine the types of microbes associated with this coral species. We also optimized a protocol for DNA extraction from coral tissues. Our taxonomic analysis revealed Proteobacteria, Bacteroidetes, Firmicutes, Cyanobacteria, and Actinomycetes as the main groups in all healthy and WP-affected corals. Vibrionales, members of the Cytophaga-Flavobacterium-Bacteroides complex, Rickettsiales, and Neisseriales were more abundant in the WP-affected corals. Diseased corals also had more eukaryotic metagenomic sequences identified as Alveolata and Apicomplexa. Our results suggest that WP disease in M. braziliensis is caused by a polymicrobial consortium.

Characterization of a new high copy Stowaway family MITE, BRAMI-1 in Brassica genome

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) are expected to play important roles in evolution of genes and genome in plants, especially in the highly duplicated plant genomes. Various MITE families and their roles in plants have been characterized. However, there have been fewer studies of MITE families and their potential roles in evolution of the recently triplicated Brassica genome.

RESULTS

We identified a new MITE family, BRAMI-1, belonging to the Stowaway super-family in the Brassica genome. In silico mapping revealed that 697 members are dispersed throughout the euchromatic regions of the B. rapa pseudo-chromosomes. Among them, 548 members (78.6%) are located in gene-rich regions, less than 3 kb from genes. In addition, we identified 516 and 15 members in the 470 Mb and 15 Mb genomic shotgun sequences currently available for B. oleracea and B. napus, respectively. The resulting estimated copy numbers for the entire genomes were 1440, 1464 and 2490 in B. rapa, B. oleracea and B. napus, respectively. Concurrently, only 70 members of the related Arabidopsis ATTIRTA-1 MITE family were identified in the Arabidopsis genome. Phylogenetic analysis revealed that BRAMI-1 elements proliferated in the Brassica genus after divergence from the Arabidopsis lineage. MITE insertion polymorphism (MIP) was inspected for 50 BRAMI-1 members, revealing high levels of insertion polymorphism between and within species of Brassica that clarify BRAMI-1 activation periods up to the present. Comparative analysis of the 71 genes harbouring the BRAMI-1 elements with their non-insertion paralogs (NIPs) showed that the BRAMI-1 insertions mainly reside in non-coding sequences and that the expression levels of genes with the elements differ from those of their NIPs.

CONCLUSION

A Stowaway family MITE, named as BRAMI-1, was gradually amplified and remained present in over than 1400 copies in each of three Brassica species. Overall, 78% of the members were identified in gene-rich regions, and it is assumed that they may contribute to the evolution of duplicated genes in the highly duplicated Brassica genome. The resulting MIPs can serve as a good source of DNA markers for Brassica crops because the insertions are highly dispersed in the gene-rich euchromatin region and are polymorphic between or within species.

CHH islands: de novo DNA methylation in near-gene chromatin regulation in maize

Small RNA-mediated regulation of chromatin structure is an important means of suppressing unwanted genetic activity in diverse plants, fungi, and animals. In plants specifically, 24-nt siRNAs direct de novo methylation to repetitive DNA, both foreign and endogenous, in a process known as RNA-directed DNA methylation (RdDM). Many components of the de novo methylation machinery have been identified recently, including multiple RNA polymerases, but specific genetic features that trigger methylation remain poorly understood. By applying whole-genome bisulfite sequencing to maize, we found that transposons close to cellular genes (particularly within 1 kb of either a gene start or end) are strongly associated with de novo methylation, as evidenced both by 24-nt siRNAs and by methylation specifically in the CHH sequence context. In addition, we found that the major classes of transposons exhibited a gradient of CHH methylation determined by proximity to genes. Our results further indicate that intergenic chromatin in maize exists in two major forms that are distinguished based on proximity to genes-one form marked by dense CG and CHG methylation and lack of transcription, and one marked by CHH methylation and activity of multiple forms of RNA polymerase. The existence of the latter, which we call CHH islands, may have implications for how cellular gene expression could be coordinated with immediately adjacent transposon repression in a large genome with a complex organization of genes interspersed in a landscape of transposons.

Comparative genomics of parasitic silkworm microsporidia reveal an association between genome expansion and host adaptation

Background

Microsporidian Nosema bombycis has received much attention because the pébrine disease of domesticated silkworms results in great economic losses in the silkworm industry. So far, no effective treatment could be found for pébrine. Compared to other known Nosema parasites, N. bombycis can unusually parasitize a broad range of hosts. To gain some insights into the underlying genetic mechanism of pathological ability and host range expansion in this parasite, a comparative genomic approach is conducted. The genome of two Nosema parasites, N. bombycis and N. antheraeae (an obligatory parasite to undomesticated silkworms Antheraea pernyi), were sequenced and compared with their distantly related species, N. ceranae (an obligatory parasite to honey bees).

Results

Our comparative genomics analysis show that the N. bombycis genome has greatly expanded due to the following three molecular mechanisms: 1) the proliferation of host-derived transposable elements, 2) the acquisition of many horizontally transferred genes from bacteria, and 3) the production of abundnant gene duplications. To our knowledge, duplicated genes derived not only from small-scale events (e.g., tandem duplications) but also from large-scale events (e.g., segmental duplications) have never been seen so abundant in any reported microsporidia genomes. Our relative dating analysis further indicated that these duplication events have arisen recently over very short evolutionary time. Furthermore, several duplicated genes involving in the cytotoxic metabolic pathway were found to undergo positive selection, suggestive of the role of duplicated genes on the adaptive evolution of pathogenic ability.

Conclusions

Genome expansion is rarely considered as the evolutionary outcome acting on those highly reduced and compact parasitic microsporidian genomes. This study, for the first time, demonstrates that the parasitic genomes can expand, instead of shrink, through several common molecular mechanisms such as gene duplication, horizontal gene transfer, and transposable element expansion. We also showed that the duplicated genes can serve as raw materials for evolutionary innovations possibly contributing to the increase of pathologenic ability. Based on our research, we propose that duplicated genes of N. bombycis should be treated as primary targets for treatment designs against pébrine.

Miniature inverted-repeat transposable elements: discovery, distribution, and activity

Eukaryotic organisms have dynamic genomes, with transposable elements (TEs) as a major contributing factor. Although the large autonomous TEs can significantly shape genomic structures during evolution, genomes often harbor more miniature nonautonomous TEs that can infest genomic niches where large TEs are rare. In spite of their cut-and-paste transposition mechanisms that do not inherently favor copy number increase, miniature inverted-repeat transposable elements (MITEs) are abundant in eukaryotic genomes and exist in high copy numbers. Based on the large number of MITE families revealed in previous studies, accurate annotation of MITEs, particularly in newly sequenced genomes, will identify more genomes highly rich in these elements. Novel families identified from these analyses, together with the currently known families, will further deepen our understanding of the origins, transposase sources, and dramatic amplification of these elements.