Genomic diversity generated by a transposable element burst in a rice recombinant inbred population.

Genomes of all characterized higher eukaryotes harbor examples of transposable element (TE) bursts-the rapid amplification of TE copies throughout a genome. Despite their prevalence, understanding how bursts diversify genomes requires the characterization of actively transposing TEs before insertion sites and structural rearrangements have been obscured by selection acting over evolutionary time. In this study, rice recombinant inbred lines (RILs), generated by crossing a bursting accession and the reference Nipponbare accession, were exploited to characterize the spread of the very active Ping/mPing family through a small population and the resulting impact on genome diversity. Comparative sequence analysis of 272 individuals led to the identification of over 14,000 new insertions of the mPing miniature inverted-repeat transposable element (MITE), with no evidence for silencing of the transposase-encoding Ping element. In addition to new insertions, Ping-encoded transposase was found to preferentially catalyze the excision of mPing loci tightly linked to a second mPing insertion. Similarly, structural variations, including deletion of rice exons or regulatory regions, were enriched for those with break points at one or both ends of linked mPing elements. Taken together, these results indicate that structural variations are generated during a TE burst as transposase catalyzes both the high copy numbers needed to distribute linked elements throughout the genome and the DNA cuts at the TE ends known to dramatically increase the frequency of recombination.

Tracking the genome-wide outcomes of a transposable element burst over decades of amplification.

To understand the success strategies of transposable elements (TEs) that attain high copy numbers, we analyzed two pairs of rice (Oryza sativa) strains, EG4/HEG4 and A119/A123, undergoing decades of rapid amplification (bursts) of the class 2 autonomous Ping element and the nonautonomous miniature inverted repeat transposable element (MITE) mPing Comparative analyses of whole-genome sequences of the two strain pairs validated that each pair has been maintained for decades as inbreds since divergence from their respective last common ancestor. Strains EG4 and HEG4 differ by fewer than 160 SNPs and a total of 264 new mPing insertions. Similarly, strains A119 and A123 exhibited about half as many SNPs (277) as new mPing insertions (518). Examination of all other potentially active TEs in these genomes revealed only a single new insertion out of ∼40,000 loci surveyed. The virtual absence of any new TE insertions in these strains outside the mPing bursts demonstrates that the Ping/mPing family gradually attains high copy numbers by maintaining activity and evading host detection for dozens of generations. Evasion is possible because host recognition of mPing sequences appears to have no impact on initiation or maintenance of the burst. Ping is actively transcribed, and both Ping and mPing can transpose despite methylation of terminal sequences. This finding suggests that an important feature of MITE success is that host recognition does not lead to the silencing of the source of transposase.

Transposition of Mutator-like transposable elements (MULEs) resembles hAT and Transib elements and V(D)J recombination.

Mutator-like transposable elements (MULEs) are widespread across fungal, plant and animal species. Despite their abundance and importance as genetic tools in plants, the transposition mechanism of the MULE superfamily was previously unknown. Discovery of the Muta1 element from Aedes aegypti and its successful transposition in yeast facilitated the characterization of key steps in Muta1 transposition. Here we show that purified transposase binds specifically to the Muta1 ends and catalyzes excision through double strand breaks (DSB) and the joining of newly excised transposon ends with target DNA. In the process, the DSB forms hairpin intermediates on the flanking DNA side. Analysis of transposase proteins containing site-directed mutations revealed the importance of the conserved DDE motif and a W residue. The transposition pathway resembles that of the V(D)J recombination reaction and the mechanism of hAT and Transib transposases including the importance of the conserved W residue in both MULEs and hATs. In addition, yeast transposition and in vitro assays demonstrated that the terminal motif and subterminal repeats of the Muta1 terminal inverted repeat also influence Muta1 transposition. Collectively, our data provides new insights to understand the evolutionary relationships between MULE, hAT and Transib elements and the V(D)J recombinase.

Fine-scale variation in meiotic recombination in Mimulus inferred from population shotgun sequencing.

Meiotic recombination rates can vary widely across genomes, with hotspots of intense activity interspersed among cold regions. In yeast, hotspots tend to occur in promoter regions of genes, whereas in humans and mice, hotspots are largely defined by binding sites of the positive-regulatory domain zinc finger protein 9. To investigate the detailed recombination pattern in a flowering plant, we use shotgun resequencing of a wild population of the monkeyflower Mimulus guttatus to precisely locate over 400,000 boundaries of historic crossovers or gene conversion tracts. Their distribution defines some 13,000 hotspots of varying strengths, interspersed with cold regions of undetectably low recombination. Average recombination rates peak near starts of genes and fall off sharply, exhibiting polarity. Within genes, recombination tracts are more likely to terminate in exons than in introns. The general pattern is similar to that observed in yeast, as well as in positive-regulatory domain zinc finger protein 9-knockout mice, suggesting that recombination initiation described here in Mimulus may reflect ancient and conserved eukaryotic mechanisms.

Unexpected consequences of a sudden and massive transposon amplification on rice gene expression.

High-copy-number transposable elements comprise the majority of eukaryotic genomes where they are major contributors to gene and genome evolution. However, it remains unclear how a host genome can survive a rapid burst of hundreds or thousands of insertions because such bursts are exceedingly rare in nature and therefore difficult to observe in real time. In a previous study we reported that in a few rice strains the DNA transposon mPing was increasing its copy number by approximately 40 per plant per generation. Here we exploit the completely sequenced rice genome to determine 1,664 insertion sites using high-throughput sequencing of 24 individual rice plants and assess the impact of insertion on the expression of 710 genes by comparative microarray analysis. We find that the vast majority of transposable element insertions either upregulate or have no detectable effect on gene transcription. This modest impact reflects a surprising avoidance of exon insertions by mPing and a preference for insertion into 5' flanking sequences of genes. Furthermore, we document the generation of new regulatory networks by a subset of mPing insertions that render adjacent genes stress inducible. As such, this study provides evidence for models first proposed previously for the involvement of transposable elements and other repetitive sequences in genome restructuring and gene regulation.

The catalytic domain of all eukaryotic cut-and-paste transposase superfamilies.

Cut-and-paste DNA transposable elements are major components of eukaryotic genomes and are grouped into superfamilies (e.g., hAT, P) based on sequence similarity of the element-encoded transposase. The transposases from several superfamilies possess a protein domain containing an acidic amino acid triad (DDE or DDD) that catalyzes the "cut and paste" transposition reaction. However, it was unclear whether this domain was shared by the transposases from all superfamilies. Through multiple-alignment of transposase sequences from a diverse collection of previously identified and recently annotated elements from a wide range of organisms, we identified the putative DDE/D triad for all superfamilies. Furthermore, we identified additional highly conserved amino acid residues or motifs within the DDE/D domain that together form a "signature string" that is specific to each superfamily. These conserved residues or motifs were exploited as phylogenetic characters to infer evolutionary relationships among all superfamilies. The phylogenetic analysis revealed three major groups that were not previously discerned and led us to revise the classification of several currently recognized superfamilies. Taking the data together, this study suggests that all eukaryotic cut-and-paste transposable element superfamilies have a common evolutionary origin and establishes a phylogenetic framework for all future cut-and-paste transposase comparisons.

Evolution of plant cell-type-specific cis -regulatory elements.

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species ( Zea mays, Sorghum bicolor, Panicum miliaceum , and Urochloa fusca ). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1 epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Comparison of class 2 transposable elements at superfamily resolution reveals conserved and distinct features in cereal grass genomes.

BACKGROUND: Class 2 transposable elements (TEs) are the predominant elements in and around plant genes where they generate significant allelic diversity. Using the complete sequences of four grasses, we have performed a novel comparative analysis of class 2 TEs. To ensure consistent comparative analyses, we re-annotated class 2 TEs in Brachypodium distachyon, Oryza sativa (rice), Sorghum bicolor and Zea mays and assigned them to one of the five cut-and-paste superfamilies found in plant genomes (Tc1/mariner, PIF/Harbinger, hAT, Mutator, CACTA). We have focused on noncoding elements because of their abundance, and compared superfamily copy number, size and genomic distribution as well as correlation with the level of nearby gene expression. RESULTS: Our comparison revealed both unique and conserved features. First, the average length or size distribution of elements in each superfamily is largely conserved, with the shortest always being Tc1/mariner elements, followed by PIF/Harbinger, hAT, Mutator and CACTA. This order also holds for the ratio of the copy numbers of noncoding to coding elements. Second, with the exception of CACTAs, noncoding TEs are enriched within and flanking genes, where they display conserved distribution patterns, having the highest peak in the promoter region. Finally, our analysis of microarray data revealed that genes associated with Tc1/mariner and PIF/Harbinger noncoding elements have significantly higher expression levels than genes without class 2 TEs. In contrast, genes with CACTA elements have significantly lower expression than genes without class 2 TEs. CONCLUSIONS: We have achieved the most comprehensive annotation of class 2 TEs to date in these four grass genomes. Comparative analysis of this robust dataset led to the identification of several previously unknown features of each superfamily related to copy number, element size, genomic distribution and correlation with the expression levels of nearby genes. These results highlight the importance of distinguishing TE superfamilies when assessing their impact on gene and genome evolution.

The rice miniature inverted repeat transposable element mPing is an effective insertional mutagen in soybean.

Insertional mutagenesis of legume genomes such as soybean (Glycine max) should aid in identifying genes responsible for key traits such as nitrogen fixation and seed quality. The relatively low throughput of soybean transformation necessitates the use of a transposon-tagging strategy where a single transformation event will produce many mutations over a number of generations. However, existing transposon-tagging tools being used in legumes are of limited utility because of restricted transposition (Ac/Ds: soybean) or the requirement for tissue culture activation (Tnt1: Medicago truncatula). A recently discovered transposable element from rice (Oryza sativa), mPing, and the genes required for its mobilization, were transferred to soybean to determine if it will be an improvement over the other available transposon-tagging tools. Stable transformation events in soybean were tested for mPing transposition. Analysis of mPing excision at early and late embryo developmental stages revealed increased excision during late development in most transgenic lines, suggesting that transposition is developmentally regulated. Transgenic lines that produced heritable mPing insertions were identified, with the plants from the highest activity line producing at least one new insertion per generation. Analysis of the mPing insertion sites in the soybean genome revealed that features displayed in rice were retained including transposition to unlinked sites and a preference for insertion within 2.5 kb of a gene. Taken together these findings indicate that mPing has the characteristics necessary for an effective transposon-tagging resource.

MITE-Hunter: a program for discovering miniature inverted-repeat transposable elements from genomic sequences.

Miniature inverted-repeat transposable elements (MITEs) are a special type of Class 2 non-autonomous transposable element (TE) that are abundant in the non-coding regions of the genes of many plant and animal species. The accurate identification of MITEs has been a challenge for existing programs because they lack coding sequences and, as such, evolve very rapidly. Because of their importance to gene and genome evolution, we developed MITE-Hunter, a program pipeline that can identify MITEs as well as other small Class 2 non-autonomous TEs from genomic DNA data sets. The output of MITE-Hunter is composed of consensus TE sequences grouped into families that can be used as a library file for homology-based TE detection programs such as RepeatMasker. MITE-Hunter was evaluated by searching the rice genomic database and comparing the output with known rice TEs. It discovered most of the previously reported rice MITEs (97.6%), and found sixteen new elements. MITE-Hunter was also compared with two other MITE discovery programs, FINDMITE and MUST. Unlike MITE-Hunter, neither of these programs can search large genomic data sets including whole genome sequences. More importantly, MITE-Hunter is significantly more accurate than either FINDMITE or MUST as the vast majority of their outputs are false-positives.

Detailed analysis of a contiguous 22-Mb region of the maize genome.

Most of our understanding of plant genome structure and evolution has come from the careful annotation of small (e.g., 100 kb) sequenced genomic regions or from automated annotation of complete genome sequences. Here, we sequenced and carefully annotated a contiguous 22 Mb region of maize chromosome 4 using an improved pseudomolecule for annotation. The sequence segment was comprehensively ordered, oriented, and confirmed using the maize optical map. Nearly 84% of the sequence is composed of transposable elements (TEs) that are mostly nested within each other, of which most families are low-copy. We identified 544 gene models using multiple levels of evidence, as well as five miRNA genes. Gene fragments, many captured by TEs, are prevalent within this region. Elimination of gene redundancy from a tetraploid maize ancestor that originated a few million years ago is responsible in this region for most disruptions of synteny with sorghum and rice. Consistent with other sub-genomic analyses in maize, small RNA mapping showed that many small RNAs match TEs and that most TEs match small RNAs. These results, performed on approximately 1% of the maize genome, demonstrate the feasibility of refining the B73 RefGen_v1 genome assembly by incorporating optical map, high-resolution genetic map, and comparative genomic data sets. Such improvements, along with those of gene and repeat annotation, will serve to promote future functional genomic and phylogenomic research in maize and other grasses.

TARGeT: a web-based pipeline for retrieving and characterizing gene and transposable element families from genomic sequences.

Gene families compose a large proportion of eukaryotic genomes. The rapidly expanding genomic sequence database provides a good opportunity to study gene family evolution and function. However, most gene family identification programs are restricted to searching protein databases where data are often lagging behind the genomic sequence data. Here, we report a user-friendly web-based pipeline, named TARGeT (Tree Analysis of Related Genes and Transposons), which uses either a DNA or amino acid 'seed' query to: (i) automatically identify and retrieve gene family homologs from a genomic database, (ii) characterize gene structure and (iii) perform phylogenetic analysis. Due to its high speed, TARGeT is also able to characterize very large gene families, including transposable elements (TEs). We evaluated TARGeT using well-annotated datasets, including the ascorbate peroxidase gene family of rice, maize and sorghum and several TE families in rice. In all cases, TARGeT rapidly recapitulated the known homologs and predicted new ones. We also demonstrated that TARGeT outperforms similar pipelines and has functionality that is not offered elsewhere.

A Course-Based Undergraduate Research Experience in CRISPR-Cas9 Experimental Design to Support Reverse Genetic Studies in Arabidopsis thaliana.

Gene-editing tools such as CRISPR-Cas9 have created unprecedented opportunities for genetic studies in plants and animals. We designed a course-based undergraduate research experience (CURE) to train introductory biology students in the concepts and implementation of gene-editing technology as well as develop their soft skills in data management and scientific communication. We present two versions of the course that can be implemented with twice-weekly meetings over a 5-week period. In the remote-learning version, students performed homology searches, designed guide RNAs (gRNAs) and primers, and learned the principles of molecular cloning. This version is appropriate when access to laboratory equipment or in-person instruction is limited, such as during closures that have occurred in response to the COVID-19 pandemic. In person, students designed gRNAs, cloned CRISPR-Cas9 constructs, and performed genetic transformation of Arabidopsis thaliana. Students learned how to design effective gRNA pairs targeting their assigned gene with an 86% success rate. Final exams tested students' ability to apply knowledge of an unfamiliar genome database to characterize gene structure and to properly design gRNAs. Average final exam scores of ∼73% and ∼84% for in-person and remote-learning CUREs, respectively, indicated that students met learning outcomes. The highly parallel nature of the CURE makes it possible to target dozens to hundreds of genes, depending on the number of sections. Applying this approach in a sensitized mutant background enables focused reverse genetic screens for genetic suppressors or enhancers. The course can be adapted readily to other organisms or projects that employ gene editing.

Pack-MULE transposable elements mediate gene evolution in plants.

Mutator-like transposable elements (MULEs) are found in many eukaryotic genomes and are especially prevalent in higher plants. In maize, rice and Arabidopsis a few MULEs were shown to carry fragments of cellular genes. These chimaeric elements are called Pack-MULEs in this study. The abundance of MULEs in rice and the availability of most of the genome sequence permitted a systematic analysis of the prevalence and nature of Pack-MULEs in an entire genome. Here we report that there are over 3,000 Pack-MULEs in rice containing fragments derived from more than 1,000 cellular genes. Pack-MULEs frequently contain fragments from multiple chromosomal loci that are fused to form new open reading frames, some of which are expressed as chimaeric transcripts. About 5% of the Pack-MULEs are represented in collections of complementary DNA. Functional analysis of amino acid sequences and proteomic data indicate that some captured gene fragments might be functional. Comparison of the cellular genes and Pack-MULE counterparts indicates that fragments of genomic DNA have been captured, rearranged and amplified over millions of years. Given the abundance of Pack-MULEs in rice and the widespread occurrence of MULEs in all characterized plant genomes, gene fragment acquisition by Pack-MULEs might represent an important new mechanism for the evolution of genes in higher plants.

Transposition of the rice miniature inverted repeat transposable element mPing in Arabidopsis thaliana.

An active miniature inverted repeat transposable element (MITE), mPing, was discovered by computer-assisted analysis of rice genome sequence. The mPing element is mobile in rice cell culture and in a few rice strains where it has been amplified to >1,000 copies during recent domestication. However, determination of the transposase source and characterization of the mechanism of transposition have been hampered by the high copy number of mPing and the presence of several candidate autonomous elements in the rice genome. Here, we report that mPing is active in Arabidopsis thaliana, where its transposition is catalyzed by three sources of transposase from rice: the autonomous Ping and Pong elements and by a cDNA derived from a Ping transcript. In addition to transposase, the product of a second element-encoded ORF of unknown function is also required for mPing transposition. Excision of mPing in A. thaliana is usually precise, and transposed copies usually insert into unlinked sites in the genome that are preferentially in or near genes. As such, this will be a valuable assay system for the dissection of MITE transposition and a potentially powerful tagging system for gene discovery in eukaryotes.

Tracking the origin of two genetic components associated with transposable element bursts in domesticated rice.

Transposable elements (TEs) shape genome evolution through periodic bursts of amplification. In this study prior knowledge of the mPing/Ping/Pong TE family is exploited to track their copy numbers and distribution in genome sequences from 3,000 accessions of domesticated Oryza sativa (rice) and the wild progenitor Oryza rufipogon. We find that mPing bursts are restricted to recent domestication and is likely due to the accumulation of two TE components, Ping16A and Ping16A_Stow, that appear to be critical for mPing hyperactivity. Ping16A is a variant of the autonomous element with reduced activity as shown in a yeast transposition assay. Transposition of Ping16A into a Stowaway element generated Ping16A_Stow, the only Ping locus shared by all bursting accessions, and shown here to correlate with high mPing copies. Finally, we show that sustained activity of the mPing/Ping family in domesticated rice produced the components necessary for mPing bursts, not the loss of epigenetic regulation.

Efficient capture of unique sequences from eukaryotic genomes.

Cot-based cloning and sequencing (CBCS), a synthesis of Cot analysis, DNA cloning and high-throughput sequencing, promises to accelerate the study of eukaryotic genomes. In particular, CBCS will (1) permit efficient gene discovery in species with substantial quantities of repetitive DNA, (2) allow the sequence complexity (i.e. all the unique sequence information) of large genomes to be elucidated at a fraction of the cost of shotgun sequencing, and (3) enhance genome sequencing efforts by facilitating capture of low-copy sequences not secured by EST sequencing. CBCS should accelerate comparative genomics research, especially in large genomes such as those of many crops.

QTL analysis of floral traits in Louisiana iris hybrids.

The formation of hybrid zones between nascent species is a widespread phenomenon. The evolutionary consequences of hybridization are influenced by numerous factors, including the action of natural selection on quantitative trait variation. Here we examine how the genetic basis of floral traits of two species of Louisiana Irises affects the extent of quantitative trait variation in their hybrids. Quantitative trait locus (QTL) mapping was used to assess the size (magnitude) of phenotypic effects of individual QTL, the degree to which QTL for different floral traits are colocalized, and the occurrence of mixed QTL effects. These aspects of quantitative genetic variation would be expected to influence (1) the number of genetic steps (in terms of QTL substitutions) separating the parental species phenotypes; (2) trait correlations; and (3) the potential for transgressive segregation in hybrid populations. Results indicate that some Louisiana Iris floral trait QTL have large effects and QTL for different traits tend to colocalize. Transgressive variation was observed for six of nine traits, despite the fact that mixed QTL effects influence few traits. Overall, our QTL results imply that the genetic basis of floral morphology and color traits might facilitate the maintenance of phenotypic divergence between Iris fulva and Iris brevicaulis, although a great deal of phenotypic variation was observed among hybrids.

The transposable element landscape of the model legume Lotus japonicus.

The largest component of plant and animal genomes characterized to date is transposable elements (TEs). The availability of a significant amount of Lotus japonicus genome sequence has permitted for the first time a comprehensive study of the TE landscape in a legume species. Here we report the results of a combined computer-assisted and experimental analysis of the TEs in the 32.4 Mb of finished TAC clones. While computer-assisted analysis facilitated a determination of TE abundance and diversity, the availability of complete TAC sequences permitted identification of full-length TEs, which facilitated the design of tools for genomewide experimental analysis. In addition to containing all TE types found in previously characterized plant genomes, the TE component of L. japonicus contained several surprises. First, it is the second species (after Oryza sativa) found to be rich in Pack-MULEs, with >1000 elements that have captured and amplified gene fragments. In addition, we have identified what appears to be a legume-specific MULE family that was previously identified only in fungal species. Finally, the L. japonicus genome contains many hundreds, perhaps thousands of Sireviruses: Ty1/copia-like elements with an extra ORF. Significantly, several of the L. japonicus Sireviruses have recently amplified and may still be actively transposing.