Retrotranspositions in orthologous regions of closely related grass species

Calcium-Binding Protein and Polymorphism in Musa spp. Somaclones Resistant to Fusarium oxysporum

The fresh fruits of 'Grande Naine' (Cavendish AAA-Musa spp.) dominate the world market, especially in countries with a population in a situation of social vulnerability. However, Fusarium wilt, caused by the fungus Fusarium oxysporum f.sp. cubense race 4 Subtropical (Foc ST4), emerges as a serious threat to banana production, requiring the development of resistant cultivars based on biotechnological strategies, such as the induction of mutation in tissue culture. This study aimed to identify and characterize genetic variation in somaclones resistant to Fusarium oxysporum f.sp. cubense subtropical race 4 (Foc ST4), derived from 'Grand Naine' bananas, by molecular markers based on retrotransposons IRAP (Inter-retrotransposon Amplified Polymorphism) and REMAP (Retrotransposon-Microsatellite Amplified Polymorphism). Nine combinations of IRAP and six combinations of REMAP primers were used. The low number of polymorphic bands did not allow for genetic diversity studies; however, ten polymorphic bands between the somaclones and control were sequenced. Of these, three presented good base calling and were aligned, namely, 1AF, 2AF, and 3AF bands. Only the 1AF band presented function related to stress response with homology to a calcium-binding protein. These proteins act early in plant infection as secondary messengers activated by pathogen-associated molecular patterns (PAMPs), initiating the cascade of plant defense signals. The fact that this band is present in all somaclones reinforces previous assessments of their resistance to Foc ST4. The use of markers IRAP and REMAP produced polymorphic bands that can, through future primer design and field validations, accelerate the identification of resistant banana genotypes for use in banana genetic breeding programs.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Genome-wide characterization of LTR retrotransposons in the non-model deep-sea annelid Lamellibrachia luymesi

Background

Long Terminal Repeat retrotransposons (LTR retrotransposons) are mobile genetic elements composed of a few genes between terminal repeats and, in some cases, can comprise over half of a genome’s content. Available data on LTR retrotransposons have facilitated comparative studies and provided insight on genome evolution. However, data are biased to model systems and marine organisms, including annelids, have been underrepresented in transposable elements studies. Here, we focus on genome of Lamellibrachia luymesi, a vestimentiferan tubeworm from deep-sea hydrocarbon seeps, to gain knowledge of LTR retrotransposons in a deep-sea annelid.

Results

We characterized LTR retrotransposons present in the genome of L. luymesi using bioinformatic approaches and found that intact LTR retrotransposons makes up about 0.1% of L. luymesi genome. Previous characterization of the genome has shown that this tubeworm hosts several known LTR-retrotransposons. Here we describe and classify LTR retrotransposons in L. luymesi as within the Gypsy, Copia and Bel-pao superfamilies. Although, many elements fell within already recognized families (e.g., Mag, CSRN1), others formed clades distinct from previously recognized families within these superfamilies. However, approximately 19% (41) of recovered elements could not be classified. Gypsy elements were the most abundant while only 2 Copia and 2 Bel-pao elements were present. In addition, analysis of insertion times indicated that several LTR-retrotransposons were recently transposed into the genome of L. luymesi, these elements had identical LTR’s raising possibility of recent or ongoing retrotransposon activity.

Conclusions

Our analysis contributes to knowledge on diversity of LTR-retrotransposons in marine settings and also serves as an important step to assist our understanding of the potential role of retroelements in marine organisms. We find that many LTR retrotransposons, which have been inserted in the last few million years, are similar to those found in terrestrial model species. However, several new groups of LTR retrotransposons were discovered suggesting that the representation of LTR retrotransposons may be different in marine settings. Further study would improve understanding of the diversity of retrotransposons across animal groups and environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07749-1.

Transposons play an important role in the evolution and diversification of centromeres among closely related species

Centromeres are important chromosomal regions necessary for eukaryotic cell segregation and replication. Due to high amounts of tandem repeats and transposons, centromeres have been difficult to sequence in most multicellular organisms, thus their sequence structure and evolution are poorly understood. In this study, we analyzed transposons in the centromere 8 (Cen8) from the African cultivated rice (O. glaberrima) and two subspecies of the Asian cultivated rice (O. sativa), indica and japonica. We detected much higher transposon contents (>69%) in centromere regions than in the whole genomes of O. sativa ssp. japonica and O. glaberrima (~35%). We compared the three Cen8s and identified numerous recent insertions of transposons that were frequently organized into multiple-layer nested blocks, similar to nested transposons in maize. Except for the Hopi retrotransposon, all LTR retrotransposons were shared but exhibit different abundances amongst the three Cen8s. Even though a majority of the transposons were located in intergenic regions, some gene-related transposons were found and may be involved in gene diversification. Chromatin immunoprecipitated (ChIP) data analysis revealed that 165 families from both Class I and Class II transposons were found in CENH3-associated chromatin sequences. These results indicate essential roles for transposons in centromeres and that the rapid divergence of the Cen8 sequences between the two cultivated rice species was primarily caused by recent transposon insertions.

Transposable elements and genome size variations in plants

Although the number of protein-coding genes is not highly variable between plant taxa, the DNA content in their genomes is highly variable, by as much as 2,056-fold from a 1C amount of 0.0648 pg to 132.5 pg. The mean 1C-value in plants is 2.4 pg, and genome size expansion/contraction is lineage-specific in plant taxonomy. Transposable element fractions in plant genomes are also variable, as low as ~3% in small genomes and as high as ~85% in large genomes, indicating that genome size is a linear function of transposable element content. Of the 2 classes of transposable elements, the dynamics of class 1 long terminal repeat (LTR) retrotransposons is a major contributor to the 1C value differences among plants. The activity of LTR retrotransposons is under the control of epigenetic suppressing mechanisms. Also, genome-purging mechanisms have been adopted to counter-balance the genome size amplification. With a wealth of information on whole-genome sequences in plant genomes, it was revealed that several genome-purging mechanisms have been employed, depending on plant taxa. Two genera, Lilium and Fritillaria, are known to have large genomes in angiosperms. There were twice times of concerted genome size evolutions in the family Liliaceae during the divergence of the current genera in Liliaceae. In addition to the LTR retrotransposons, non-LTR retrotransposons and satellite DNAs contributed to the huge genomes in the two genera by possible failure of genome counter-balancing mechanisms.

Phylogenetic analysis of the expansion of the MATH-BTB gene family in the grasses

MATH-BTB proteins are known to act as substrate-specific adaptors of cullin3 (CUL3)-based ubiquitin E3 ligases to target protein for ubiquitination. In a previous study we reported the presence of 31 MATH-BTB genes in the maize genome and determined the regulatory role of the MATH-BTB protein MAB1 during meiosis to mitosis transition. In contrast to maize, there are only 6 homologous genes in the model plant Arabidopsis, while this family has largely expanded in grasses. Here, we report a phylogenetic analysis of the MATH-BTB gene family in 9 land plant species including various mosses, eudicots, and grasses. We extend a previous classification of the plant MATH-BTB family and additionally arrange the expanded group into 5 grass-specific clades. Synteny studies indicate that expansion occurred to a large extent due to local gene duplications. Expression studies of 3 closely related MATH-BTB genes in maize (MAB1-3) indicate highly specific expression pattern. In summary, this work provides a solid base for further studies comparing genetic and functional information of the MATH-BTB family especially in the grasses.

Co-evolution of plant LTR-retrotransposons and their host genomes

Transposable elements (TEs), particularly, long terminal repeat retrotransposons (LTR-RTs), are the most abundant DNA components in all plant species that have been investigated, and are largely responsible for plant genome size variation. Although plant genomes have experienced periodic proliferation and/or recent burst of LTR-retrotransposons, the majority of LTR-RTs are inactivated by DNA methylation and small RNA-mediated silencing mechanisms, and/or were deleted/truncated by unequal homologous recombination and illegitimate recombination, as suppression mechanisms that counteract genome expansion caused by LTR-RT amplification. LTR-RT DNA is generally enriched in pericentromeric regions of the host genomes, which appears to be the outcomes of preferential insertions of LTR-RTs in these regions and low effectiveness of selection that purges LTR-RT DNA from these regions relative to chromosomal arms. Potential functions of various TEs in their host genomes remain blurry; nevertheless, LTR-RTs have been recognized to play important roles in maintaining chromatin structures and centromere functions and regulation of gene expressions in their host genomes.

Epiallele biogenesis in maize

We have correlated cytosine methylation of two epialleles, P1-rr and P1-pr, with variation in gene expression and therefore phenotype. The p1 gene in maize encodes a transcription factor that controls phlobaphene pigment accumulation in floral tissues. While cytosine methylation was assayed in various regions spanning 17 kb, the only difference in DNA methylation pattern between the expressed P1-rr allele and the silenced P1-pr allele was detected in a region that consists of a complex arrangement of transposons and adjacent repeats. This region, which comprises the distal enhancer element of P1-rr, is hypermethylated in P1-pr compared to P1-rr. Based on other precedents, we hypothesize that DNA methylation spreads from the transposable elements into the flanking P1-rr enhancer, thereby transcriptionally silencing the gene. Interestingly, P1-pr is reactivated in mutants of the dominant epigenetic modifier Ufo1. DNA methylation in the distal enhancer sequence is significantly reduced, which inversely correlates with increased transcript levels and pigmentation in P1-pr Ufo1 plants. If in general DNA methylation spreads from transposons into adjacent sequences containing regulatory elements for neighboring genes, the corresponding genes could be silenced by chance. Given the large amount of transposable elements in the maize genome, epialleles may be far more frequent than previously estimated.

Mutator System Derivatives Isolated from Sugarcane Genome Sequence

Mutator-like transposase is the most represented transposon transcript in the sugarcane transcriptome. Phylogenetic reconstructions derived from sequenced transcripts provided evidence that at least four distinct classes exist (I-IV) and that diversification among these classes occurred early in Angiosperms, prior to the divergence of Monocots/Eudicots. The four previously described classes served as probes to select and further sequence six BAC clones from a genomic library of cultivar R570. A total of 579,352 sugarcane base pairs were produced from these "Mutator system" BAC containing regions for further characterization. The analyzed genomic regions confirmed that the predicted structure and organization of the Mutator system in sugarcane is composed of two true transposon lineages, each containing a specific terminal inverted repeat and two transposase lineages considered to be domesticated. Each Mutator transposase class displayed a particular molecular structure supporting lineage specific evolution. MUSTANG, previously described domesticated genes, are located in syntenic regions across Sacharineae and, as expected for a host functional gene, posses the same gene structure as in other Poaceae. Two sequenced BACs correspond to hom(eo)logous locus with specific retrotransposon insertions that discriminate sugarcane haplotypes. The comparative studies presented, add information to the Mutator systems previously identified in the maize and rice genomes by describing lineage specific molecular structure and genomic distribution pattern in the sugarcane genome. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s12042-012-9104-y) contains supplementary material, which is available to authorized users.

Genetic and epigenetic dynamics of a retrotransposon after allopolyploidization of wheat

Allopolyploidy, or the combination of two or more distinct genomes in one nucleus, is usually accompanied by radical genomic changes involving transposable elements (TEs). The dynamics of TEs after an allopolyploidization event are poorly understood. In this study, we analyzed the methylation state and genetic rearrangements of a high copied, newly amplified terminal-repeat retrotransposon in miniature (TRIM) family in wheat termed Veju. We found that Veju insertion sites underwent massive methylation changes in the first four generations of a newly formed wheat allohexaploid. Hypomethylation or hypermethylation occurred in ∼43% of the tested insertion sites; while hypomethylation was significantly predominant in the first three generations of the newly formed allohexaploid, hypermethylation became predominant in the subsequent generation. In addition, we determined that the methylation state of Veju long terminal repeats (LTRs) might be correlated with the deletion and/or insertion of the TE. While most of the methylation changes and deletions of Veju occurred in the first generation of the newly formed allohexaploid, most Veju insertions were seen in the second generation. Finally, using quantitative PCR, we quantitatively assessed the genome composition of Veju in the newly formed allohexaploid and found that up to 50% of Veju LTRs were deleted in the first generation. Retrotransposition bursts in subsequent generations, however, led to increases in Veju elements. In light of these findings, the underlying mechanisms of TRIM rearrangements are discussed.

Genome Size Dynamics and Evolution in Monocots

Monocot genomic diversity includes striking variation at many levels. This paper compares various genomic characters (e.g., range of chromosome numbers and ploidy levels, occurrence of endopolyploidy, GC content, chromosome packaging and organization, genome size) between monocots and the remaining angiosperms to discern just how distinctive monocot genomes are. One of the most notable features of monocots is their wide range and diversity of genome sizes, including the species with the largest genome so far reported in plants. This genomic character is analysed in greater detail, within a phylogenetic context. By surveying available genome size and chromosome data it is apparent that different monocot orders follow distinctive modes of genome size and chromosome evolution. Further insights into genome size-evolution and dynamics were obtained using statistical modelling approaches to reconstruct the ancestral genome size at key nodes across the monocot phylogenetic tree. Such approaches reveal that while the ancestral genome size of all monocots was small ( pg), there have been several major increases and decreases during monocot evolution. In addition, notable increases in the rates of genome size-evolution were found in Asparagales and Poales compared with other monocot lineages.

Derepression of the plant Chromovirus LORE1 induces germline transposition in regenerated plants

Transposable elements represent a large proportion of the eukaryotic genomes. Long Terminal Repeat (LTR) retrotransposons are very abundant and constitute the predominant family of transposable elements in plants. Recent studies have identified chromoviruses to be a widely distributed lineage of Gypsy elements. These elements contain chromodomains in their integrases, which suggests a preference for insertion into heterochromatin. In turn, this preference might have contributed to the patterning of heterochromatin observed in host genomes. Despite their potential importance for our understanding of plant genome dynamics and evolution, the regulatory mechanisms governing the behavior of chromoviruses and their activities remain largely uncharacterized. Here, we report a detailed analysis of the spatio-temporal activity of a plant chromovirus in the endogenous host. We examined LORE1a, a member of the endogenous chromovirus LORE1 family from the model legume Lotus japonicus. We found that this chromovirus is stochastically de-repressed in plant populations regenerated from de-differentiated cells and that LORE1a transposes in the male germline. Bisulfite sequencing of the 5′ LTR and its surrounding region suggests that tissue culture induces a loss of epigenetic silencing of LORE1a. Since LTR promoter activity is pollen specific, as shown by the analysis of transgenic plants containing an LTR::GUS fusion, we conclude that male germline-specific LORE1a transposition in pollen grains is controlled transcriptionally by its own cis-elements. New insertion sites of LORE1a copies were frequently found in genic regions and show no strong insertional preferences. These distinctive novel features of LORE1 indicate that this chromovirus has considerable potential for generating genetic and epigenetic diversity in the host plant population. Our results also define conditions for the use of LORE1a as a genetic tool.

In contrast to animals, where germline differentiation initiates early in embryogenesis, germline differentiation in plants starts in the adult phase during reproductive development. Transpositions of transposable elements in both somatic and gametic cells can be transmitted to the next generation. As a result, plant genomes may contain transposable elements exhibiting a variety of tissue-specific activities. Thus far, the spatio-temporal activity of LTR retrotransposons, the most abundant class of transposable elements in plants, has not been well characterized. Here, we report a detailed analysis of the spatio-temporal transposition pattern of a plant LTR retrotransposon in the endogenous system. Using the model legume Lotus japonicus, we found that LORE1a, a member of the chromovirus LORE1 family that belongs to the Gypsy superfamily, was epigenetically de-repressed via tissue culture. Activation was stochastic and derepression was maintained in regenerated plants. This feature made it possible to trace the original spatio-temporal activity of the retrotransposon in the intact plants. We determined that the plant chromovirus retrotransposes mainly in the male germline, without obvious insertional preferences for chromosomal regions. This finding suggests that the tissue specificity of transposable elements should be taken into account when considering their impact on the host genome dynamics and evolution.

Small RNAs originated from pseudogenes: cis- or trans-acting?

Pseudogenes are significant components of eukaryotic genomes, and some have acquired novel regulatory roles. To date, no study has characterized rice pseudogenes systematically or addressed their impact on the structure and function of the rice genome. In this genome-wide study, we have identified 11,956 non-transposon-related rice pseudogenes, most of which are from gene duplications. About 12% of the rice protein-coding genes, half of which are in singleton families, have a pseudogene paralog. Interestingly, we found that 145 of these pseudogenes potentially gave rise to antisense small RNAs after examining approximately 1.5 million small RNAs from developing rice grains. The majority (>50%) of these antisense RNAs are 24-nucleotides long, a feature often seen in plant repeat-associated small interfering RNAs (siRNAs) produced by RNA-dependent RNA polymerase (RDR2) and Dicer-like protein 3 (DCL3), suggesting that some pseudogene-derived siRNAs may be implicated in repressing pseudogene transcription (i.e., cis-acting). Multiple lines of evidence, however, indicate that small RNAs from rice pseudogenes might also function as natural antisense siRNAs either by interacting with the complementary sense RNAs from functional parental genes (38 cases) or by forming double-strand RNAs with transcripts of adjacent paralogous pseudogenes (2 cases) (i.e., trans-acting). Further examinations of five additional small RNA libraries revealed that pseudogene-derived antisense siRNAs could be produced in specific rice developmental stages or physiological growth conditions, suggesting their potentially important roles in normal rice development. In summary, our results show that pseudogenes derived from protein-coding genes are prevalent in the rice genome, and a subset of them are strong candidates for producing small RNAs with novel regulatory roles. Our findings suggest that pseudogenes of exapted functions may be a phenomenon ubiquitous in eukaryotic organisms.

Evolution of teleost fish retroviruses: characterization of new retroviruses with cellular genes

The interactions between retroviruses and their hosts can be of a beneficial or detrimental nature. Some endogenous retroviruses are involved in development, while others cause disease. The Genome Parsing Suite (GPS) is a software tool to track and trace all Retroid agents in any sequenced genome (M. A. McClure et al., Genomics 85:512-523, 2005). Using the GPS, the retroviral content was assessed in four model teleost fish. Eleven new species of fish retroviruses are identified and characterized. The reverse transcriptase protein sequences were used to reconstruct a fish retrovirus phylogeny, thereby, significantly expanding the epsilon-retrovirus family. Most of these novel retroviruses encode additional genes, some of which are homologous to cellular genes that would confer viral advantage. Although the fish divergence is much more ancient, retroviruses began infecting fish genomes approximately 4 million years ago.

Change of gene structure and function by non-homologous end-joining, homologous recombination, and transposition of DNA

An important objective in genome research is to relate genome structure to gene function. Sequence comparisons among orthologous and paralogous genes and their allelic variants can reveal sequences of functional significance. Here, we describe a 379-kb region on chromosome 1 of maize that enables us to reconstruct chromosome breakage, transposition, non-homologous end-joining, and homologous recombination events. Such a high-density composition of various mechanisms in a small chromosomal interval exemplifies the evolution of gene regulation and allelic diversity in general. It also illustrates the evolutionary pace of changes in plants, where many of the above mechanisms are of somatic origin. In contrast to animals, somatic alterations can easily be transmitted through meiosis because the germline in plants is contiguous to somatic tissue, permitting the recovery of such chromosomal rearrangements. The analyzed region contains the P1-wr allele, a variant of the genetically well-defined p1 gene, which encodes a Myb-like transcriptional activator in maize. The P1-wr allele consists of eleven nearly perfect P1-wr 12-kb repeats that are arranged in a tandem head-to-tail array. Although a technical challenge to sequence such a structure by shotgun sequencing, we overcame this problem by subcloning each repeat and ordering them based on nucleotide variations. These polymorphisms were also critical for recombination and expression analysis in presence and absence of the trans-acting epigenetic factor Ufo1. Interestingly, chimeras of the p1 and p2 genes, p2/p1 and p1/p2, are framing the P1-wr cluster. Reconstruction of sequence amplification steps at the p locus showed the evolution from a single Myb-homolog to the multi-gene P1-wr cluster. It also demonstrates how non-homologous end-joining can create novel gene fusions. Comparisons to orthologous regions in sorghum and rice also indicate a greater instability of the maize genome, probably due to diploidization following allotetraploidization.

Methylation-sensitive linking libraries enhance gene-enriched sequencing of complex genomes and map DNA methylation domains

Background

Many plant genomes are resistant to whole-genome assembly due to an abundance of repetitive sequence, leading to the development of gene-rich sequencing techniques. Two such techniques are hypomethylated partial restriction (HMPR) and methylation spanning linker libraries (MSLL). These libraries differ from other gene-rich datasets in having larger insert sizes, and the MSLL clones are designed to provide reads localized to "epigenetic boundaries" where methylation begins or ends.

Results

A large-scale study in maize generated 40,299 HMPR sequences and 80,723 MSLL sequences, including MSLL clones exceeding 100 kb. The paired end reads of MSLL and HMPR clones were shown to be effective in linking existing gene-rich sequences into scaffolds. In addition, it was shown that the MSLL clones can be used for anchoring these scaffolds to a BAC-based physical map. The MSLL end reads effectively identified epigenetic boundaries, as indicated by their preferential alignment to regions upstream and downstream from annotated genes. The ability to precisely map long stretches of fully methylated DNA sequence is a unique outcome of MSLL analysis, and was also shown to provide evidence for errors in gene identification. MSLL clones were observed to be significantly more repeat-rich in their interiors than in their end reads, confirming the correlation between methylation and retroelement content. Both MSLL and HMPR reads were found to be substantially gene-enriched, with the SalI MSLL libraries being the most highly enriched (31% align to an EST contig), while the HMPR clones exhibited exceptional depletion of repetitive DNA (to ~11%). These two techniques were compared with other gene-enrichment methods, and shown to be complementary.

Conclusion

MSLL technology provides an unparalleled approach for mapping the epigenetic status of repetitive blocks and for identifying sequences mis-identified as genes. Although the types and natures of epigenetic boundaries are barely understood at this time, MSLL technology flags both approximate boundaries and methylated genes that deserve additional investigation. MSLL and HMPR sequences provide a valuable resource for maize genome annotation, and are a uniquely valuable complement to any plant genome sequencing project. In order to make these results fully accessible to the community, a web display was developed that shows the alignment of MSLL, HMPR, and other gene-rich sequences to the BACs; this display is continually updated with the latest ESTs and BAC sequences.

Dynamics and differential proliferation of transposable elements during the evolution of the B and A genomes of wheat

Transposable elements (TEs) constitute >80% of the wheat genome but their dynamics and contribution to size variation and evolution of wheat genomes (Triticum and Aegilops species) remain unexplored. In this study, 10 genomic regions have been sequenced from wheat chromosome 3B and used to constitute, along with all publicly available genomic sequences of wheat, 1.98 Mb of sequence (from 13 BAC clones) of the wheat B genome and 3.63 Mb of sequence (from 19 BAC clones) of the wheat A genome. Analysis of TE sequence proportions (as percentages), ratios of complete to truncated copies, and estimation of insertion dates of class I retrotransposons showed that specific types of TEs have undergone waves of differential proliferation in the B and A genomes of wheat. While both genomes show similar rates and relatively ancient proliferation periods for the Athila retrotransposons, the Copia retrotransposons proliferated more recently in the A genome whereas Gypsy retrotransposon proliferation is more recent in the B genome. It was possible to estimate for the first time the proliferation periods of the abundant CACTA class II DNA transposons, relative to that of the three main retrotransposon superfamilies. Proliferation of these TEs started prior to and overlapped with that of the Athila retrotransposons in both genomes. However, they also proliferated during the same periods as Gypsy and Copia retrotransposons in the A genome, but not in the B genome. As estimated from their insertion dates and confirmed by PCR-based tracing analysis, the majority of differential proliferation of TEs in B and A genomes of wheat (87 and 83%, respectively), leading to rapid sequence divergence, occurred prior to the allotetraploidization event that brought them together in Triticum turgidum and Triticum aestivum, <0.5 million years ago. More importantly, the allotetraploidization event appears to have neither enhanced nor repressed retrotranspositions. We discuss the apparent proliferation of TEs as resulting from their insertion, removal, and/or combinations of both evolutionary forces.

Genome-wide analysis of transposon insertion polymorphisms reveals intraspecific variation in cultivated rice

Insertions and precise eliminations of transposable elements generated numerous transposon insertion polymorphisms (TIPs) in rice (Oryza sativa). We observed that TIPs represent more than 50% of large insertions and deletions (>100 bp) in the rice genome. Using a comparative genomic approach, we identified 2,041 TIPs between the genomes of two cultivars, japonica Nipponbare and indica 93-11. We also identified 691 TIPs between Nipponbare and indica Guangluai 4 in the 23-Mb collinear regions of chromosome 4. Among them, retrotransposon-based insertion polymorphisms were used to reveal the evolutionary relationships of these three cultivars. Our conservative estimates suggest that the TIPs generated approximately 14% of the genomic DNA sequence differences between subspecies indica and japonica. It was also found that more than 10% of TIPs were located in expressed gene regions, representing an important source of genetic variation. Transcript evidence implies that these TIPs induced a series of genetic differences between two subspecies, including interrupting host genes, creating different expression forms, drastically changing intron length, and affecting expression levels of adjacent genes. These analyses provide genome-wide insights into evolutionary history and genetic variation of rice.

LTR retrotransposon landscape in Medicago truncatula: more rapid removal than in rice

Background

Long terminal repeat retrotransposons (LTR elements) are ubiquitous Eukaryotic TEs that transpose through RNA intermediates. Accounting for significant proportion of many plant genomes, LTR elements have been well established as one of the major forces underlying the evolution of plant genome size, structure and function. The accessibility of more than 40% of genomic sequences of the model legume Medicago truncatula (Mt) has made the comprehensive study of its LTR elements possible.

Results

We use a newly developed tool LTR_FINDER to identify LTR retrotransposons in the Mt genome and detect 526 full-length elements as well as a great number of copies related to them. These elements constitute about 9.6% of currently available genomic sequences. They are classified into 85 families of which 64 are reported for the first time. The majority of the LTR retrotransposons belong to either Copia or Gypsy superfamily and the others are categorized as TRIMs or LARDs by their length. We find that the copy-number of Copia-like families is 3 times more than that of Gypsy-like ones but the latter contribute more to the genome. The analysis of PBS and protein-coding domain structure of the LTR families reveals that they tend to use only 4–5 types of tRNAs and many families have quite conservative ORFs besides known TE domains. For several important families, we describe in detail their abundance, conservation, insertion time and structure. We investigate the amplification-deletion pattern of the elements and find that the detectable full-length elements are relatively young and most of them were inserted within the last 0.52 MY. We also estimate that more than ten million bp of the Mt genomic sequences have been removed by the deletion of LTR elements and the removal of the full-length structures in Mt has been more rapid than in rice.

Conclusion

This report is the first comprehensive description and analysis of LTR retrotransposons in the Mt genome. Many important novel LTR families were discovered and their evolution is elucidated. Our results may outline the LTR retrotransposon landscape of the model legume.

Physical and genetic structure of the maize genome reflects its complex evolutionary history

Maize (Zea mays L.) is one of the most important cereal crops and a model for the study of genetics, evolution, and domestication. To better understand maize genome organization and to build a framework for genome sequencing, we constructed a sequence-ready fingerprinted contig-based physical map that covers 93.5% of the genome, of which 86.1% is aligned to the genetic map. The fingerprinted contig map contains 25,908 genic markers that enabled us to align nearly 73% of the anchored maize genome to the rice genome. The distribution pattern of expressed sequence tags correlates to that of recombination. In collinear regions, 1 kb in rice corresponds to an average of 3.2 kb in maize, yet maize has a 6-fold genome size expansion. This can be explained by the fact that most rice regions correspond to two regions in maize as a result of its recent polyploid origin. Inversions account for the majority of chromosome structural variations during subsequent maize diploidization. We also find clear evidence of ancient genome duplication predating the divergence of the progenitors of maize and rice. Reconstructing the paleoethnobotany of the maize genome indicates that the progenitors of modern maize contained ten chromosomes.

As a cash crop and a model biological system, maize is of great public interest. To facilitate maize molecular breeding and its basic biology research, we built a high-resolution physical map with two different fingerprinting methods on the same set of bacterial artificial chromosome clones. The physical map was integrated to a high-density genetic map and further serves as a framework for the maize genome-sequencing project. Comparative genomics showed that the euchromatic regions between rice and maize are very conserved. Physically we delimited these conserved regions and thus detected many genome rearrangements. We defined extensively the duplication blocks within the maize genome. These blocks allowed us to reconstruct the chromosomes of the maize progenitor. We detected that maize genome has experienced two rounds of genome duplications, an ancient one before maize–rice divergence and a recent one after tetraploidization.

Transcriptome Analysis of ESTs from a Chaetognath Reveals a Deep-Branching Clade of Retrovirus-Like Retrotransposons

Chaetognaths constitute a small marine phylum exhibiting several characteristic which are highly unusual in animal genomes, including two classes of both rRNA and protein ribosomal genes. As in this phylum presence of retrovirus-like elements has never been documented, analysis of a published expressed sequence tag (EST) collection of the chaetognath Spadella cephaloptera has been made. Twelve sequences representing transcript sections of reverse transcriptase domain of active retrotransposons were isolated from~11,000 ESTs. Five of them are originated from Gypsy retrovirus-like elements, whereas the other are transcripts from a Bel-Pao LTR-retrotransposon, a Penelope-like element and LINE retrotransposons. Moreover, a part of a putative integrase has also been found. Phylogenetic analyses suggest a deep-branching clade of the retrovirus-like elements, which is in agreement with the probably Cambrian origin of the phylum. Moreover, retrotransposons have not been found in telomeric-like transcripts which are probably constituted by both vertebrate and arthropod canonical repeats.

The evolution of transposon repeat-induced point mutation in the genome of Colletotrichum cereale: reconciling sex, recombination and homoplasy in an ''asexual" pathogen

Mobile transposable elements are among the primary drivers of the evolution of eukaryotic genomes. For fungi, repeat-induced point mutation (RIP) silencing minimizes deleterious effects of transposons by mutating multicopy DNA during meiosis. In this study we identify five transposon species from the mitosporic fungus Colletotrichum cereale and report the signature pattern of RIP acting in a lineage-specific manner on 21 of 35 unique transposon copies, providing the first evidence for sexual recombination for this species. Sequence analysis of genomic populations of the retrotransposon Ccret2 showed repeated rounds of RIP mutation acting on different copies of the element. In the RIPped Ccret2 population, there were multiple inferences of incongruence primarily attributed to RIP-induced homoplasy. This study supports the view that the sequence variability of transposon populations in filamentous fungi reflects the activities of evolutionary processes that fall outside of typical phylogenetic or population genetic reconstructions.

New insights into Oryza genome evolution: high gene colinearity and differential retrotransposon amplification

An approximately 247-kb genomic region from FF genome of wild rice Oryza brachyantha, possessing the smallest Oryza genome, was compared to the orthologous approximately 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 approximately 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.

Genome-wide comparative analysis of copia retrotransposons in Triticeae, rice, and Arabidopsis reveals conserved ancient evolutionary lineages and distinct dynamics of individual copia families

Although copia retrotransposons are major components of all plant genomes, the evolutionary relationships between individual copia families and between elements from different plant species are only poorly studied. We used 20 copia families from the large-genome plants barley and wheat to identify 46 families of homologous copia elements from rice and 22 from Arabidopsis, two plant species with much smaller genomes. In total, 599 copia elements were analyzed. Phylogenetic analysis showed that copia elements from the four species can be classified into six ancient lineages that existed before the divergence of monocots and dicots. The six lineages show a surprising degree of conservation in sequence organization and other characteristics across species. Additionally, the phylogenetic data suggest at least one case of horizontal gene transfer between the Arabidopsis and rice lineages. Insertion time estimates for 522 high-copy elements showed that retrotransposons from rice were active at different times in waves of activity lasting 0.5-2 million years, depending on the family, whereas elements from wheat and barley had longer periods of activity. We estimated that half of the rice copia elements are truncated or otherwise rearranged after approximately 790,000 yr, which is almost twice the half-life of Arabidopsis elements. In contrast, wheat and barley copia elements appear to have a massively longer half-life, beyond our ability to estimate from the available data. These findings suggest that genome size can be explained by the specific rate of DNA removal from the genome and the length of active periods of retrotransposon families.

Analysis of EST sequences suggests recent origin of allotetraploid colonial and creeping bentgrasses

Advances in plant genomics have permitted the analysis of several members of the grass family, including the major domesticated species, and provided new insights into the evolution of the major crops on earth. Two members, colonial bentgrass (Agrostis capillaris L.) and creeping bentgrass (A. stolonifera L.) have only recently been domesticated and provide an interesting case of polyploidy and comparison to crops that have undergone human selection for thousands of years. As an initial step of characterizing these genomes, we have sampled roughly 10% of their gene content, thereby also serving as a starting point for the construction of their physical and genetic maps. Sampling mRNA from plants subjected to environmental stress showed a remarkable increase in transcription of transposable elements. Both colonial and creeping bentgrass are allotetraploids and are considered to have one genome in common, designated the A2 genome. Analysis of conserved genes present among the ESTs suggests the colonial and creeping bentgrass A2 genomes diverged from a common ancestor approximately 2.2 million years ago (MYA), thereby providing an enhanced evolutionary zoom in respect to the origin of maize, which formed 4.8 MYA, and tetraploid wheat, which formed only 0.5 MYA and is the progenitor of domesticated hexaploid wheat.

Comparison of orthologous loci from small grass genomes Brachypodium and rice: implications for wheat genomics and grass genome annotation

Brachypodium sylvaticum and Brachypodium distachyon were recently proposed as new model plants because of their small genomes and their phylogenetic position between rice and Triticeae crops. We sequenced a 371-kb region in B. sylvaticum, the largest genomic sequence available so far from this species, providing quantitative data on gene conservation, collinearity and phylogeny. We compared it with orthologous regions from rice and wheat. Brachypodium and wheat show perfect macro-collinearity of genetic markers, whereas rice contains an approximately 220-kb inversion. Rice contains almost twice as many genes as Brachypodium in the region studied, whereas wheat has about 40% more. Through comparative annotation, we identified alternative transcripts and improved the annotation for several rice genes, indicating that approximately 15% of rice genes might require re-annotation. Surprisingly, our data suggest that 10-15% of functional sequences in small grass genomes may not encode any proteins. From available genomic and expressed sequence tag sequences, we estimated Brachypodium to have diverged from wheat about 35-40 Mya, significantly more recently than the divergence of rice and wheat. However, our data also indicate that orthologous regions from Brachypodium and wheat differ considerably in gene content, thus the Brachypodium genome sequence probably cannot replace genomic studies in the large Triticeae genomes.

Maize haplotype with a helitron-amplified cytidine deaminase gene copy

Background

Genetic maps are based on recombination of orthologous gene sequences between different strains of the same species. Therefore, it was unexpected to find extensive non-collinearity of genes between different inbred strains of maize. Interestingly, disruption of gene collinearity can be caused among others by a rolling circle-type copy and paste mechanism facilitated by Helitrons. However, understanding the role of this type of gene amplification has been hampered by the lack of finding intact gene sequences within Helitrons.

Results

By aligning two haplotypes of the z1C1 locus of maize we found a Helitron that contains two genes, one encoding a putative cytidine deaminase and one a hypothetical protein with part of a 40S ribosomal protein. The cytidine deaminase gene, called ZmCDA3, has been copied from the ZmCDA1 gene on maize chromosome 7 about 4.5 million years ago (mya) after maize was formed by whole-genome duplication from two progenitors. Inbred lines contain gene copies of both progenitors, the ZmCDA1 and ZmCDA2 genes. Both genes diverged when the progenitors of maize split and are derived from the same progenitor as the rice OsCDA1 gene. The ZmCDA1 and ZmCDA2 genes are both transcribed in leaf and seed tissue, but transcripts of the paralogous ZmCDA3 gene have not been found yet. Based on their protein structure the maize CDA genes encode a nucleoside deaminase that is found in bacterial systems and is distinct from the mammalian RNA and/or DNA modifying enzymes.

Conclusion

The conservation of a paralogous gene sequence encoding a cytidine deaminase gene over 4.5 million years suggests that Helitrons could add functional gene sequences to new chromosomal positions and thereby create new haplotypes. However, the function of such paralogous gene copies cannot be essential because they are not present in all maize strains. However, it is interesting to note that maize hybrids can outperform their inbred parents. Therefore, certain haplotypes may function only in combination with other haplotypes or under specialized environmental conditions.