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

Cis-regulatory elements (CREs) are critical in regulating gene expression, and yet understanding of CRE evolution remains challenging. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa, integrating data from 103,911 nuclei representing 126 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-derived 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 silencer-like CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of plant cell-type-specific CRE evolution.

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.

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 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.

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.

RelocaTE2: a high resolution transposable element insertion site mapping tool for population resequencing.

BACKGROUND: Transposable element (TE) polymorphisms are important components of population genetic variation. The functional impacts of TEs in gene regulation and generating genetic diversity have been observed in multiple species, but the frequency and magnitude of TE variation is under appreciated. Inexpensive and deep sequencing technology has made it affordable to apply population genetic methods to whole genomes with methods that identify single nucleotide and insertion/deletion polymorphisms. However, identifying TE polymorphisms, particularly transposition events or non-reference insertion sites can be challenging due to the repetitive nature of these sequences, which hamper both the sensitivity and specificity of analysis tools. METHODS: We have developed the tool RelocaTE2 for identification of TE insertion sites at high sensitivity and specificity. RelocaTE2 searches for known TE sequences in whole genome sequencing reads from second generation sequencing platforms such as Illumina. These sequence reads are used as seeds to pinpoint chromosome locations where TEs have transposed. RelocaTE2 detects target site duplication (TSD) of TE insertions allowing it to report TE polymorphism loci with single base pair precision. RESULTS AND DISCUSSION: The performance of RelocaTE2 is evaluated using both simulated and real sequence data. RelocaTE2 demonstrate high level of sensitivity and specificity, particularly when the sequence coverage is not shallow. In comparison to other tools tested, RelocaTE2 achieves the best balance between sensitivity and specificity. In particular, RelocaTE2 performs best in prediction of TSDs for TE insertions. Even in highly repetitive regions, such as those tested on rice chromosome 4, RelocaTE2 is able to report up to 95% of simulated TE insertions with less than 0.1% false positive rate using 10-fold genome coverage resequencing data. RelocaTE2 provides a robust solution to identify TE insertion sites and can be incorporated into analysis workflows in support of describing the complete genotype from light coverage genome sequencing.

Functional characterization of the active Mutator-like transposable element, Muta1 from the mosquito Aedes aegypti.

BACKGROUND: Mutator-like transposable elements (MULEs) are widespread with members in fungi, plants, and animals. Most of the research on the MULE superfamily has focused on plant MULEs where they were discovered and where some are extremely active and have significant impact on genome structure. The maize MuDR element has been widely used as a tool for both forward and reverse genetic studies because of its high transposition rate and preference for targeting genic regions. However, despite being widespread, only a few active MULEs have been identified, and only one, the rice Os3378, has demonstrated activity in a non-host organism. RESULTS: Here we report the identification of potentially active MULEs in the mosquito Aedes aegypti. We demonstrate that one of these, Muta1, is capable of excision and reinsertion in a yeast transposition assay. Element reinsertion generated either 8 bp or 9 bp target site duplications (TSDs) with no apparent sequence preference. Mutagenesis analysis of donor site TSDs in the yeast assay indicates that their presence is important for precise excision and enhanced transposition. Site directed mutagenesis of the putative DDE catalytic motif and other conserved residues in the transposase protein abolished transposition activity. CONCLUSIONS: Collectively, our data indicates that the Muta1 transposase of Ae. aegypti can efficiently catalyze both excision and reinsertion reactions in yeast. Mutagenesis analysis reveals that several conserved amino acids, including the DDE triad, play important roles in transposase function. In addition, donor site TSD also impacts the transposition of Muta1.

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.

The use of RelocaTE and unassembled short reads to produce high-resolution snapshots of transposable element generated diversity in rice.

Transposable elements (TEs) are dynamic components of genomes that often vary in copy number among members of the same species. With the advent of next-generation sequencing TE insertion-site polymorphism can be examined at an unprecedented level of detail when combined with easy-to-use bioinformatics software. Here we report a new tool, RelocaTE, that rapidly identifies specific TE insertions that are either polymorphic or shared between a reference and unassembled next-generation sequencing reads. Furthermore, a novel companion tool, CharacTErizer, exploits the depth of coverage to classify genotypes of nonreference insertions as homozygous, heterozygous or, when analyzing an active TE family, as rare somatic insertion or excision events. It does this by comparing the numbers of RelocaTE aligned reads to reads that map to the same genomic position without the TE. Although RelocaTE and CharacTErizer can be used for any TE, they were developed to analyze the very active mPing element which is undergoing massive amplification in specific strains of Oryza sativa (rice). Three individuals of one of these strains, A123, were resequenced and analyzed for mPing insertion site polymorphisms. The majority of mPing insertions found (~97%) are not present in the reference, and two siblings from a self-crossed of this strain were found to share only ~90% of their insertions. Private insertions are primarily heterozygous but include both homozygous and predicted somatic insertions. The reliability of the predicted genotypes was validated by polymerase chain reaction.

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.

Transposing from the laboratory to the classroom to generate authentic research experiences for undergraduates.

Large lecture classes and standardized laboratory exercises are characteristic of introductory biology courses. Previous research has found that these courses do not adequately convey the process of scientific research and the excitement of discovery. Here we propose a model that provides beginning biology students with an inquiry-based, active learning laboratory experience. The Dynamic Genome course replicates a modern research laboratory focused on eukaryotic transposable elements where beginning undergraduates learn key genetics concepts, experimental design, and molecular biological skills. Here we report on two key features of the course, a didactic module and the capstone original research project. The module is a modified version of a published experiment where students experience how virtual transposable elements from rice (Oryza sativa) are assayed for function in transgenic Arabidopsis thaliana. As part of the module, students analyze the phenotypes and genotypes of transgenic plants to determine the requirements for transposition. After mastering the skills and concepts, students participate in an authentic research project where they use computational analysis and PCR to detect transposable element insertion site polymorphism in a panel of diverse maize strains. As a consequence of their engagement in this course, students report large gains in their ability to understand the nature of research and demonstrate that they can apply that knowledge to independent research projects.

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.

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.