Genome-wide mutagenesis of Zea mays L. using RescueMu transposons

Maize mutant screens: from classical methods to new CRISPR-based approaches

Mutations play a pivotal role in shaping the trajectory and outcomes of a species evolution and domestication. Maize (Zea mays) has been a major staple crop and model for genetic research for more than 100 yr. With the arrival of site-directed mutagenesis and genome editing (GE) driven by the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), maize mutational research is once again in the spotlight. If we combine the powerful physiological and genetic characteristics of maize with the already available and ever increasing toolbox of CRISPR-Cas, prospects for its future trait engineering are very promising. This review aimed to give an overview of the progression and learnings of maize screening studies analyzing forward genetics, natural variation and reverse genetics to focus on recent GE approaches. We will highlight how each strategy and resource has contributed to our understanding of maize natural and induced trait variability and how this information could be used to design the next generation of mutational screenings.

Mutator transposon insertions within maize genes often provide a novel outward reading promoter

The highly active family of Mutator (Mu) DNA transposons has been widely used for forward and reverse genetics in maize. There are examples of Mu-suppressible alleles that result in conditional phenotypic effects based on the activity of Mu. Phenotypes from these Mu-suppressible mutations are observed in Mu-active genetic backgrounds, but absent when Mu activity is lost. For some Mu-suppressible alleles, phenotypic suppression likely results from an outward-reading promoter within Mu that is only active when the autonomous Mu element is silenced or lost. We isolated 35 Mu alleles from the UniformMu population that represent insertions in 24 different genes. Most of these mutant alleles are due to insertions within gene coding sequences, but several 5' UTR and intron insertions were included. RNA-seq and de novo transcript assembly were utilized to document the transcripts produced from 33 of these Mu insertion alleles. For 20 of the 33 alleles, there was evidence of transcripts initiating within the Mu sequence reading through the gene. This outward-reading promoter activity was detected in multiple types of Mu elements and does not depend on the orientation of Mu. Expression analyses of Mu-initiated transcripts revealed the Mu promoter often provides gene expression levels and patterns that are similar to the wild-type gene. These results suggest the Mu promoter may represent a minimal promoter that can respond to gene cis-regulatory elements. Findings from this study have implications for maize researchers using the UniformMu population, and more broadly highlight a strategy for transposons to co-exist with their host.

Independent evolution of transposase and TIRs facilitated by recombination between Mutator transposons from divergent clades in maize

Nearly all eukaryotes carry DNA transposons of the Robertson's Mutator (Mu) superfamily, a widespread source of genome instability and genetic variation. Despite their pervasive impact on host genomes, much remains unknown about the evolution of these transposons. Transposase recognition of terminal inverted repeats (TIRs) is thought to drive and constrain coevolution of MuDR transposase genes and TIRs. To address the extent of this relationship and its impact, we compared separate phylogenies of TIRs and MuDR gene sequences from Mu elements in the maize genome. Five major clades were identified. As expected, most Mu elements were bound by highly similar TIRs from the same clade (homomorphic type). However, a subset of elements contained dissimilar TIRs derived from divergent clades. These "heteromorphs" typically occurred in multiple copies indicating active transposition in the genome. In addition, analysis of internal sequences showed that exchanges between elements having divergent TIRs produced new mudra and mudrb gene combinations. In several instances, TIR homomorphs had been regenerated within a heteromorph clade with retention of distinctive internal MuDR sequence combinations. Results reveal that recombination between divergent clades facilitates independent evolution of transposase (mudra), transposase-binding targets (TIRs), and capacity for insertion (mudrb) of active Mu elements. This mechanism would be enhanced by the preference of Mu insertions for recombination-rich regions near the 5' ends of genes. We suggest that cycles of recombination give rise to alternating homo- and heteromorph forms that enhance the diversity on which selection for Mu fitness can operate.

Boundary domain genes were recruited to suppress bract growth and promote branching in maize

Grass inflorescence development is diverse and complex and involves sophisticated but poorly understood interactions of genes regulating branch determinacy and leaf growth. Here, we use a combination of transcript profiling and genetic and phylogenetic analyses to investigate tasselsheath1 (tsh1) and tsh4, two maize genes that simultaneously suppress inflorescence leaf growth and promote branching. We identify a regulatory network of inflorescence leaf suppression that involves the phase change gene tsh4 upstream of tsh1 and the ligule identity gene liguleless2 (lg2). We also find that a series of duplications in the tsh1 gene lineage facilitated its shift from boundary domain in nongrasses to suppressed inflorescence leaves of grasses. Collectively, these results suggest that the boundary domain genes tsh1 and lg2 were recruited to inflorescence leaves where they suppress growth and regulate a nonautonomous signaling center that promotes inflorescence branching, an important component of yield in cereal grasses.

A Novel Mutator-Like Transposable Elements With Unusual Structure and Recent Transpositions in Barley (Hordeum vulgare)

Mutator-like transposable elements (MULEs) represent a unique superfamily of DNA transposons as they can capture host genes and cause higher frequency of mutations in some eukaryotes. Despite their essential roles in plant evolution and functional genomics, MULEs are not fully understood yet in many important crops including barley (Hordeum vulgare). In this study, we analyzed the barley genome and identified a new mutator transposon Hvu_Abermu. This transposon is present at extremely high copy number in barley and shows unusual structure as it contains three open reading frames (ORFs) including one ORF (ORF1) encoding mutator transposase protein and one ORF (ORFR) showing opposite transcriptional orientation. We identified homologous sequences of Hvu_Abermu in both monocots and dicots and grouped them into a large mutator family named Abermu. Abermu transposons from different species share significant sequence identity, but they exhibit distinct sequence structures. Unlike the transposase proteins which are highly conserved between Abermu transposons from different organisms, the ORFR-encoded proteins are quite different from distant species. Phylogenetic analysis indicated that Abermu transposons shared closer evolutionary relationships with the maize MuDR transposon than other reported MULEs. We also found phylogenetic incongruence for the Abermu transposons identified in rice and its wild species implying the possibility of horizontal transfer of transposon. Further comparison indicated that over 200 barley genes contain Abermu-related sequences. We analyzed the barley pan genomes and detected polymorphic Hvu_Abermu transposons between the sequenced 23 wild and cultivated barley genomes. Our efforts identified a novel mutator transposon and revealed its recent transposition activity, which may help to develop genetic tools for barley and other crops.

BonnMu: A Sequence-Indexed Resource of Transposon-Induced Maize Mutations for Functional Genomics Studies

Sequence-indexed insertional libraries in maize (Zea mays) are fundamental resources for functional genetics studies. Here, we constructed a Mutator (Mu) insertional library in the B73 inbred background designated BonnMu A total of 1,152 Mu-tagged F2-families were sequenced using the Mu-seq approach. We detected 225,936 genomic Mu insertion sites and 41,086 high quality germinal Mu insertions covering 16,392 of the annotated maize genes (37% of the B73v4 genome). On average, each F2-family of the BonnMu libraries captured 37 germinal Mu insertions in genes of the Filtered Gene Set (FGS). All BonnMu insertions and phenotypic seedling photographs of Mu-tagged F2-families can be accessed via MaizeGDB.org Downstream examination of 137,410 somatic and germinal insertion sites revealed that 50% of the tagged genes have a single hotspot, targeted by Mu By comparing our BonnMu (B73) data to the UniformMu (W22) library, we identified conserved insertion hotspots between different genetic backgrounds. Finally, the vast majority of BonnMu and UniformMu transposons was inserted near the transcription start site of genes. Remarkably, 75% of all BonnMu insertions were in closer proximity to the transcription start site (distance: 542 bp) than to the start codon (distance: 704 bp), which corresponds to open chromatin, especially in the 5' region of genes. Our European sequence-indexed library of Mu insertions provides an important resource for functional genetics studies of maize.

A Sequence-Indexed Mutator Insertional Library for Maize Functional Genomics Study

Sequence-indexed insertional libraries are important resources for functional gene study in model plants. However, the maize (Zea mays) UniformMu library covers only 36% of the annotated maize genes. Here, we generated a new sequence-indexed maize Mutator insertional library named ChinaMu through high-throughput sequencing of enriched Mu-tagged sequences. A total of 2,581 Mu F2 lines were analyzed, and 311,924 nonredundant Mu insertion sites were obtained. Based on experimental validation, ChinaMu contains about 97,000 germinal Mu insertions, about twice as many as UniformMu. About two-thirds (66,565) of the insertions are high-quality germinal insertions (positive rate > 90%), 89.6% of which are located in genic regions. Furthermore, 45.7% (20,244) of the 44,300 annotated maize genes are effectively tagged and about two-thirds (13,425) of these genes harbor multiple insertions. We tested the utility of ChinaMu using pentatricopeptide repeat (PPR) genes. For published PPR genes with defective kernel phenotypes, 17 out of 20 were tagged, 11 of which had the previously reported mutant phenotype. For 16 unstudied PPR genes with both Mu insertions and defective kernel phenotypes, 6 contained insertions that cosegregated with the mutant phenotype. Our sequence-indexed Mu insertional library provides an important resource for functional genomics study in maize.

CsMYB60 is a key regulator of flavonols and proanthocyanidans that determine the colour of fruit spines in cucumber

Spine colour is an important fruit quality trait that influences the commercial value of cucumber (Cucumis sativus). However, little is known about the metabolites and the regulatory mechanisms of their biosynthesis in black spine varieties. In this study, we determined that the pigments of black spines are flavonoids, including flavonols and proanthocyanidins (PAs). We identified CsMYB60 as the best candidate for the previously identified B (Black spine) locus. Expression levels of CsMYB60 and the key genes involved in flavonoid biosynthesis were higher in black-spine inbred lines than that in white-spine lines at different developmental stages. The insertion of a Mutator-like element (CsMULE) in the second intron of CsMYB60 decreased its expression in a white-spine line. Transient overexpression assays indicated that CsMYB60 is a key regulatory gene and Cs4CL is a key structural gene in the pigmentation of black spines. In addition, the DNA methylation level in the CsMYB60 promoter was much lower in the black-spine line compared with white-spine line. The CsMULE insert may decrease the expression level of CsMYB60, causing hindered synthesis of flavonols and PAs in cucumber fruit spines.

Control of sexuality by the sk1-encoded UDP-glycosyltransferase of maize

Sex determination in maize involves the production of staminate and pistillate florets from an initially bisexual floral meristem. Pistil elimination in staminate florets requires jasmonic acid signaling, and functional pistils are protected by the action of the silkless 1 (sk1) gene. The sk1 gene was identified and found to encode a previously uncharacterized family 1 uridine diphosphate glycosyltransferase that localized to the plant peroxisomes. Constitutive expression of an sk1 transgene protected all pistils in the plant, causing complete feminization, a gain-of-function phenotype that operates by blocking the accumulation of jasmonates. The segregation of an sk1 transgene was used to effectively control the production of pistillate and staminate inflorescences in maize plants.

Mutator and MULE Transposons

The Mutator system of transposable elements (TEs) is a highly mutagenic family of transposons in maize. Because they transpose at high rates and target genic regions, these transposons can rapidly generate large numbers of new mutants, which has made the Mutator system a favored tool for both forward and reverse mutagenesis in maize. Low copy number versions of this system have also proved to be excellent models for understanding the regulation and behavior of Class II transposons in plants. Notably, the availability of a naturally occurring locus that can heritably silence autonomous Mutator elements has provided insights into the means by which otherwise active transposons are recognized and silenced. This chapter will provide a review of the biology, regulation, evolution and uses of this remarkable transposon system, with an emphasis on recent developments in our understanding of the ways in which this TE system is recognized and epigenetically silenced as well as recent evidence that Mu-like elements (MULEs) have had a significant impact on the evolution of plant genomes.

Get the jump - Do 3'UTRs protect transposable elements from silencing?

Eukaryotic genomes contain large numbers of transposable elements and repetitive sequences that are subjected to silencing through epigenetic mechanisms. These involve primarily DNA methylation, chromatin modifications and small RNA. It is known that these transposable elements can affect the expression of neighboring genes; however, little is known about how transposable element silencing depends on the general chromosomal environment at the insertion site. Taking advantage of the vast genomic resources available in Arabidopsis thaliana, a recent report begins to unravel these interactions by identifying insertion sites of one specific MULE element, AtMu1c across the A. thaliana lineage. Among over 30 insertion sites analyzed, a correlation between the loss of epigenetic silencing and the insertion into the 3'end of protein coding genes was found. Here, we discuss details, implications and potential mechanisms of these findings that may be applicable to a much wider set of transposable elements and across diverse species.

Transposition of a rice Mutator-like element in the yeast Saccharomyces cerevisiae

Mutator-like transposable elements (MULEs) are widespread in plants and are well known for their high transposition activity as well as their ability to duplicate and amplify host gene fragments. Despite their abundance and importance, few active MULEs have been identified. In this study, we demonstrated that a rice (Oryza sativa) MULE, Os3378, is capable of excising and reinserting in yeast (Saccharomyces cerevisiae), suggesting that yeast harbors all the host factors for the transposition of MULEs. The transposition activity induced by the wild-type transposase is low but can be altered by modification of the transposase sequence, including deletion, fusion, and substitution. Particularly, fusion of a fluorescent protein to the transposase enhanced the transposition activity, representing another approach to manipulate transposases. Moreover, we identified a critical region in the transposase where the net charge of the amino acids seems to be important for activity. Finally, transposition efficiency is also influenced by the element and its flanking sequences (i.e., small elements are more competent than their large counterparts). Perfect target site duplication is favorable, but not required, for precise excision. In addition to the potential application in functional genomics, this study provides the foundation for further studies of the transposition mechanism of MULEs.

Identification of the maize gravitropism gene lazy plant1 by a transposon-tagging genome resequencing strategy

Since their initial discovery, transposons have been widely used as mutagens for forward and reverse genetic screens in a range of organisms. The problems of high copy number and sequence divergence among related transposons have often limited the efficiency at which tagged genes can be identified. A method was developed to identity the locations of Mutator (Mu) transposons in the Zea mays genome using a simple enrichment method combined with genome resequencing to identify transposon junction fragments. The sequencing library was prepared from genomic DNA by digesting with a restriction enzyme that cuts within a perfectly conserved motif of the Mu terminal inverted repeats (TIR). Paired-end reads containing Mu TIR sequences were computationally identified and chromosomal sequences flanking the transposon were mapped to the maize reference genome. This method has been used to identify Mu insertions in a number of alleles and to isolate the previously unidentified lazy plant1 (la1) gene. The la1 gene is required for the negatively gravitropic response of shoots and mutant plants lack the ability to sense gravity. Using bioinformatic and fluorescence microscopy approaches, we show that the la1 gene encodes a cell membrane and nuclear localized protein. Our Mu-Taq method is readily adaptable to identify the genomic locations of any insertion of a known sequence in any organism using any sequencing platform.

Mu-seq: sequence-based mapping and identification of transposon induced mutations

Mutations tagged by transposon insertions can be readily mapped and identified in organisms with sequenced genomes. Collections of such mutants allow a systematic analysis of gene function, and can be sequence-indexed to build invaluable resources. Here we present Mu-seq (Mutant-seq), a high-throughput NextGen sequencing method for harnessing high-copy transposons. We illustrate the efficacy of Mu-seq by applying it to the Robertson's Mutator system in a large population of maize plants. A single Mu-seq library, for example, constructed from 576 different families (2304 plants), enabled 4, 723 novel, germinal, transposon insertions to be detected, identified, and mapped with single base-pair resolution. In addition to the specificity, efficiency, and reproducibility of Mu-seq, a key feature of this method is its adjustable scale that can accommodate simultaneous profiling of transposons in thousands of individuals. We also describe a Mu-seq bioinformatics framework tailored to high-throughput, genome-wide, and population-wide analysis of transposon insertions.

TED, an autonomous and rare maize transposon of the mutator superfamily with a high gametophytic excision frequency

Mutator (Mu) elements, one of the most diverse superfamilies of DNA transposons, are found in all eukaryotic kingdoms, but are particularly numerous in plants. Most of the present knowledge on the transposition behavior of this superfamily comes from studies of the maize (Zea mays) Mu elements, whose transposition is mediated by the autonomous Mutator-Don Robertson (MuDR) element. Here, we describe the maize element TED (for Transposon Ellen Dempsey), an autonomous cousin that differs significantly from MuDR. Element excision and reinsertion appear to require both proteins encoded by MuDR, but only the single protein encoded by TED. Germinal excisions, rare with MuDR, are common with TED, but arise in one of the mitotic divisions of the gametophyte, rather than at meiosis. Instead, transposition-deficient elements arise at meiosis, suggesting that the double-strand breaks produced by element excision are repaired differently in mitosis and meiosis. Unlike MuDR, TED is a very low-copy transposon whose number and activity do not undergo dramatic changes upon inbreeding or outcrossing. Like MuDR, TED transposes mostly to unlinked sites and can form circular transposition products. Sequences closer to TED than to MuDR were detected only in the grasses, suggesting a rather recent evolutionary split from a common ancestor.

Cytological characterization and allelism testing of anther developmental mutants identified in a screen of maize male sterile lines

Proper regulation of anther differentiation is crucial for producing functional pollen, and defects in or absence of any anther cell type result in male sterility. To deepen understanding of processes required to establish premeiotic cell fate and differentiation of somatic support cell layers a cytological screen of maize male-sterile mutants has been conducted which yielded 42 new mutants including 22 mutants with premeiotic cytological defects (increasing this class fivefold), 7 mutants with postmeiotic defects, and 13 mutants with irregular meiosis. Allelism tests with known and new mutants confirmed new alleles of four premeiotic developmental mutants, including two novel alleles of msca1 and single new alleles of ms32, ms8, and ocl4, and two alleles of the postmeiotic ms45. An allelic pair of newly described mutants was found. Premeiotic mutants are now classified into four categories: anther identity defects, abnormal anther structure, locular wall defects and premature degradation of cell layers, and/or microsporocyte collapse. The range of mutant phenotypic classes is discussed in comparison with developmental genetic investigation of anther development in rice and Arabidopsis to highlight similarities and differences between grasses and eudicots and within the grasses.

Regulation of the Mutator system of transposons in maize

The Mutator system has proved to be an invaluable tool for elucidating gene function via insertional mutagenesis. Its high copy number, high transposition frequency, relative lack of insertion specificity, and ease of use has made it the preferred method for gene tagging in maize. Recent advances in high throughput sequencing of insertion sites, combined with the availability of large numbers of pre-mutagenized and sequence-indexed stocks, ensure that this resource will only be more useful in the years ahead. Muk is a locus that can silence Mu-active lines, making it possible to ameliorate the phenotypic effects of high numbers of active Mu transposons and reduce the copy number of these elements during introgressions.

Using MuDR/Mu transposons in directed tagging strategies

An introduction to MuDR/Mu transposons as mutagens is provided along with protocols for using these elements to tag maize genes. Selection for retention of Mutator activity is described as well as details for establishing and screening tagging populations efficiently.

A highly conserved, small LTR retrotransposon that preferentially targets genes in grass genomes

LTR retrotransposons are often the most abundant components of plant genomes and can impact gene and genome evolution. Most reported LTR retrotransposons are large elements (>4 kb) and are most often found in heterochromatic (gene poor) regions. We report the smallest LTR retrotransposon found to date, only 292 bp. The element is found in rice, maize, sorghum and other grass genomes, which indicates that it was present in the ancestor of grass species, at least 50-80 MYA. Estimated insertion times, comparisons between sequenced rice lines, and mRNA data indicate that this element may still be active in some genomes. Unlike other LTR retrotransposons, the small LTR retrotransposons (SMARTs) are distributed throughout the genomes and are often located within or near genes with insertion patterns similar to MITEs (miniature inverted repeat transposable elements). Our data suggests that insertions of SMARTs into or near genes can, in a few instances, alter both gene structures and gene expression. Further evidence for a role in regulating gene expression, SMART-specific small RNAs (sRNAs) were identified that may be involved in gene regulation. Thus, SMARTs may have played an important role in genome evolution and genic innovation and may provide a valuable tool for gene tagging systems in grass.

Mu killer-Mediated and Spontaneous Silencing of Zea mays Mutator Family Transposable Elements Define Distinctive Paths of Epigenetic Inactivation

Mu killer contains a partial inverted duplication of the mudrA transposase gene and two copies of the terminal inverted repeat A (TIRA) region of the master MuDR element of maize. Mu killer can effectively silence single copy MuDR/Mu lines, and it is proposed that a ∼4 kb hairpin RNA is generated by read through transcription from a flanking gene and that this transcript serves as a substrate for siRNA production. Mu killer was sequenced, except for a recalcitrant portion in the center of the locus, and shown to be co-linear with mudrA as originally proposed. The ability of the dominant Mu killer locus to silence a standard high copy number MuDR/Mu transposon line was evaluated. After two generations of exposure, about three quarters of individuals were silenced indicating reasonable effectiveness as measured by the absence of mudrA transposase transcripts. Mu killer individuals that effectively silenced MuDR expressed two short antisense transcripts. In contrast, Mu killer individuals that failed to silence MuDR expressed multiple sense transcripts, derived from read through transcription initiating in a flanking gene, but no antisense transcripts were detected.

Whole organ, venation and epidermal cell morphological variations are correlated in the leaves of Arabidopsis mutants

Despite the large number of genes known to affect leaf shape or size, we still have a relatively poor understanding of how leaf morphology is established. For example, little is known about how cell division and cell expansion are controlled and coordinated within a growing leaf to eventually develop into a laminar organ of a definite size. To obtain a global perspective of the cellular basis of variations in leaf morphology at the organ, tissue and cell levels, we studied a collection of 111 non-allelic mutants with abnormally shaped and/or sized leaves, which broadly represent the mutational variations in Arabidopsis thaliana leaf morphology not associated with lethality. We used image-processing techniques on these mutants to quantify morphological parameters running the gamut from the palisade mesophyll and epidermal cells to the venation, whole leaf and rosette levels. We found positive correlations between epidermal cell size and leaf area, which is consistent with long-standing Avery's hypothesis that the epidermis drives leaf growth. In addition, venation parameters were positively correlated with leaf area, suggesting that leaf growth and vein patterning share some genetic controls. Positional cloning of the genes affected by the studied mutations will eventually establish functional links between genotypes, molecular functions, cellular parameters and leaf phenotypes.

Phenotypic and genomic analyses of a fast neutron mutant population resource in soybean

Mutagenized populations have become indispensable resources for introducing variation and studying gene function in plant genomics research. In this study, fast neutron (FN) radiation was used to induce deletion mutations in the soybean (Glycine max) genome. Approximately 120,000 soybean seeds were exposed to FN radiation doses of up to 32 Gray units to develop over 23,000 independent M2 lines. Here, we demonstrate the utility of this population for phenotypic screening and associated genomic characterization of striking and agronomically important traits. Plant variation was cataloged for seed composition, maturity, morphology, pigmentation, and nodulation traits. Mutants that showed significant increases or decreases in seed protein and oil content across multiple generations and environments were identified. The application of comparative genomic hybridization (CGH) to lesion-induced mutants for deletion mapping was validated on a midoleate x-ray mutant, M23, with a known FAD2-1A (for fatty acid desaturase) gene deletion. Using CGH, a subset of mutants was characterized, revealing deletion regions and candidate genes associated with phenotypes of interest. Exome resequencing and sequencing of PCR products confirmed FN-induced deletions detected by CGH. Beyond characterization of soybean FN mutants, this study demonstrates the utility of CGH, exome sequence capture, and next-generation sequencing approaches for analyses of mutant plant genomes. We present this FN mutant soybean population as a valuable public resource for future genetic screens and functional genomics research.

An Ac transposon system based on maize chromosome 4S for isolating long-distance-transposed Ac tags in the maize genome

Transposon tagging is an important tool for gene isolation and functional studies. In maize, several transposon-tagging systems have been developed, mostly using Activator/Dissociation (Ac/Ds) and Mutator systems. Here, we establish another Ac-based transposon system with the donor Ac tightly linked with sugary1 (su1) on maize chromosome 4S. Newly transposed Ac (tr-Acs) were detected based on a negative dosage effect, and long-distance-transposed Ac events were identified and isolated from the donor Ac by a simple backcross scheme. In this study, we identified 208 independent long-distance-transposed Ac lines. Thirty-one flanking sequences of these tr-Acs were isolated and localized in the maize genome. As found in previous studies, the tr-Acs preferentially inserted into genic sequences. The distribution of tr-Acs is not random. In our study, the tr-Acs preferentially transposed into chromosomes 1, 2, 9 and 10. We discuss the preferential distribution of tr-Acs from Ac systems. Our system is complementary to two other Ac-based regional-mutagenesis systems in maize, and the combined use of these systems will achieve an even and high-density distribution of Ac elements throughout the maize genome for functional-genomics studies.

[Construction of a methylation filtration library in Hevea brasiliensis.]

In order to enrich gene encoding region of Hevea brasiliensis, a methylation filtration library was constructed using Escherichia coli McrBC restriction-modification system. The titers of the non-amplified library and the amplified library were 2.6×106 pfu/ml and 9.0×109, respectively. The rate of positive clones was 86.4%. The lengths of inserted DNA sequence ranged from 1 kb to 2.5 kb and the average size of inserts was 1.2 kb. One hundred clones were selected randomly for sequencing, resulting in splicing out of 81 non-redundant sequences, including 6 contigs and 75 singlets. The redundancy was 17.35%. Blast analysis showed that 39.5% of non-redundant sequences were homologous with the Nr database, 14.81% with the EST database, and 32.1% were unknown sequences. Some sequences were related genes for flowering, insect and disease resistance. Therefore, the rubber tree methylation library is helpful for discovery and cloning of functional genes.

Distributed simple sequence repeat markers for efficient mapping from maize public mutagenesis populations

The genome sequence of the B73 maize inbred enables map-based cloning of genetic variants underlying phenotypes. In parallel to sequencing efforts, multiple public mutagenesis resources are being developed predominantly in the W22 and B73 inbreds. Efficient platforms to map mutants in these genetic backgrounds would aid molecular genetic analysis of the public resources. We screened 505 simple sequence repeat markers for polymorphisms between the B73, Mo17, and W22 inbreds. Using common thermocycling conditions, 47.1% of the markers showed co-dominant polymorphisms in at least one pair of inbreds. Based on these results, we identified 85 distributed markers for mapping in all three inbred pairs. For each inbred pair, the distributed set has 64-71 polymorphic markers with a mean distance of 27-29 cM between markers. The distributed markers give nearly complete coverage of the genetic map for each inbred pair. We demonstrate the utility of the marker set for efficient placement of mutants on the maize genetic map with an example mapping experiment of a seed mutant from the UniformMu mutagenesis resource. We conclude that these distributed molecular markers enable rapid mapping of phenotypic variants from public mutagenesis populations.

Genome-wide distribution of transposed Dissociation elements in maize

The maize (Zea mays) transposable element Dissociation (Ds) was mobilized for large-scale genome mutagenesis and to study its endogenous biology. Starting from a single donor locus on chromosome 10, over 1500 elements were distributed throughout the genome and positioned on the maize physical map. Genetic strategies to enrich for both local and unlinked insertions were used to distribute Ds insertions. Global, regional, and local insertion site trends were examined. We show that Ds transposed to both linked and unlinked sites and displayed a nonuniform distribution on the genetic map around the donor r1-sc:m3 locus. Comparison of Ds and Mutator insertions reveals distinct target preferences, which provide functional complementarity of the two elements for gene tagging in maize. In particular, Ds displays a stronger preference for insertions within exons and introns, whereas Mutator insertions are more enriched in promoters and 5'-untranslated regions. Ds has no strong target site consensus sequence, but we identified properties of the DNA molecule inherent to its local structure that may influence Ds target site selection. We discuss the utility of Ds for forward and reverse genetics in maize and provide evidence that genes within a 2- to 3-centimorgan region flanking Ds insertions will serve as optimal targets for regional mutagenesis.

Genetic resources for maize cell wall biology

Grass species represent a major source of food, feed, and fiber crops and potential feedstocks for biofuel production. Most of the biomass is contributed by cell walls that are distinct in composition from all other flowering plants. Identifying cell wall-related genes and their functions underpins a fundamental understanding of growth and development in these species. Toward this goal, we are building a knowledge base of the maize (Zea mays) genes involved in cell wall biology, their expression profiles, and the phenotypic consequences of mutation. Over 750 maize genes were annotated and assembled into gene families predicted to function in cell wall biogenesis. Comparative genomics of maize, rice (Oryza sativa), and Arabidopsis (Arabidopsis thaliana) sequences reveal differences in gene family structure between grass species and a reference eudicot species. Analysis of transcript profile data for cell wall genes in developing maize ovaries revealed that expression within families differed by up to 100-fold. When transcriptional analyses of developing ovaries before pollination from Arabidopsis, rice, and maize were contrasted, distinct sets of cell wall genes were expressed in grasses. These differences in gene family structure and expression between Arabidopsis and the grasses underscore the requirement for a grass-specific genetic model for functional analyses. A UniformMu population proved to be an important resource in both forward- and reverse-genetics approaches to identify hundreds of mutants in cell wall genes. A forward screen of field-grown lines by near-infrared spectroscopic screen of mature leaves yielded several dozen lines with heritable spectroscopic phenotypes. Pyrolysis-molecular beam mass spectrometry confirmed that several nir mutants had altered carbohydrate-lignin compositions.

Mu transposon insertion sites and meiotic recombination events co-localize with epigenetic marks for open chromatin across the maize genome

The Mu transposon system of maize is highly active, with each of the ∼50–100 copies transposing on average once each generation. The approximately one dozen distinct Mu transposons contain highly similar ∼215 bp terminal inverted repeats (TIRs) and generate 9-bp target site duplications (TSDs) upon insertion. Using a novel genome walking strategy that uses these conserved TIRs as primer binding sites, Mu insertion sites were amplified from Mu stocks and sequenced via 454 technology. 94% of ∼965,000 reads carried Mu TIRs, demonstrating the specificity of this strategy. Among these TIRs, 21 novel Mu TIRs were discovered, revealing additional complexity of the Mu transposon system. The distribution of >40,000 non-redundant Mu insertion sites was strikingly non-uniform, such that rates increased in proportion to distance from the centromere. An identified putative Mu transposase binding consensus site does not explain this non-uniformity. An integrated genetic map containing more than 10,000 genetic markers was constructed and aligned to the sequence of the maize reference genome. Recombination rates (cM/Mb) are also strikingly non-uniform, with rates increasing in proportion to distance from the centromere. Mu insertion site frequencies are strongly correlated with recombination rates. Gene density does not fully explain the chromosomal distribution of Mu insertion and recombination sites, because pronounced preferences for the distal portion of chromosome are still observed even after accounting for gene density. The similarity of the distributions of Mu insertions and meiotic recombination sites suggests that common features, such as chromatin structure, are involved in site selection for both Mu insertion and meiotic recombination. The finding that Mu insertions and meiotic recombination sites both concentrate in genomic regions marked with epigenetic marks of open chromatin provides support for the hypothesis that open chromatin enhances rates of both Mu insertion and meiotic recombination.

Genomic insertion sites of Mu transposons were amplified and sequenced via next generation technology, revealing more than 40,000 non-redundant Mu insertion sites that are non-uniformly distributed across the maize genome and within genes. Along chromosomes, frequencies of Mu transposon insertions are strongly correlated with recombination rates. Although both Mu and recombination occur preferentially in genes, gene density does not fully explain these patterns. Instead, the finding that Mu insertions and meiotic recombination sites both concentrate in genomic regions marked with epigenetic marks of open chromatin provides support for the hypothesis that open chromatin enhances rates of both Mu insertion and meiotic recombination.

DLA-based strategies for cloning insertion mutants: cloning the gl4 locus of maize using Mu transposon tagged alleles

Digestion-ligation-amplification (DLA), a novel adaptor-mediated PCR-based method that uses a single-stranded oligo as the adaptor, was developed to overcome difficulties of amplifying unknown sequences flanking known DNA sequences in large genomes. DLA specifically overcomes the problems associated with existing methods for amplifying genomic sequences flanking Mu transposons, including high levels of nonspecific amplification. Two DLA-based strategies, MuClone and DLA-454, were developed to isolate Mu-tagged alleles. MuClone allows for the amplification of subsets of the numerous Mu transposons in the genome, using unique three-nucleotide tags at the 3' ends of primers, simplifying the identification of flanking sequences that cosegregate with mutant phenotypes caused by Mu insertions. DLA-454, which combines DLA with 454 pyrosequencing, permits the efficient cloning of genes for which multiple independent insertion alleles are available without the need to develop segregating populations. The utility of each approach was validated by independently cloning the gl4 (glossy4) gene. Mutants of gl4 lack the normal accumulation of epicuticular waxes. The gl4 gene is a homolog of the Arabidopsis CUT1 gene, which encodes a condensing enzyme involved in the synthesis of very-long-chain fatty acids, which are precursors of epicuticular waxes.

[Advances of gene enrichment in plant genome]

The genome size varies greatly in higher plants. Repetitive sequences account for most of the large plant genomes while low-copy or single copy genic sequences, referred to as gene space, take up only a small portion of the genomes. Considering the large amount of repetitive sequences, it is a great challenge to obtain genic sequences using high-throughout methods in non-model plants bearing large genomes. Currently, several approaches have been developed for gene enrichment on a genome-wide scale, such as cDNA library, methylation filtration library, high Cot library and transposon tagging. Here, we reviewed the technical principles, advantages and disadvantages of these methods, as well as the recent development of methylation filtration technology. An in-depth discussion was performed for selection of one method or combination of methods according to the research objectives and plant materials, especially for plants with large genomes.

Mutator transposon activity reprograms the transcriptomes and proteomes of developing maize anthers

Despite the high conservation of anther gene expression patterns across maize lines, Mu transposition programmed by transcriptionally active MuDR results in a 25% change in the transcriptome, monitored over 90 h of immature anther development, without altering the morphology, anatomy or pace of development. Most transcriptome changes are stage specific: cases of suppression of normal transcripts and ectopic activation are equally represented. Protein abundance changes were validated for numerous metabolic enzymes, and highlight the increased carbon and reactive oxygen management in Mutator anthers. Active Mutator lines appear to experience chronic stress, on a par with abiotic treatments that stimulate early flowering. Despite the diversity of acclimation responses, anther development progresses normally, in contrast to male-sterile mutants that disrupt anther cell fate or function completely, and cause fewer transcriptome changes. The early flowering phenotype ultimately confers an advantage in Mu element transmission.

Plasmid rescue: recovery of flanking genomic sequences from transgenic transposon insertion sites

The transgenic RescueMu lines were designed for and successfully used in our maize gene discovery project. The pBluescript-containing RescueMu transposon can be readily recovered by a procedure called plasmid rescue. Plasmid rescue is a technique for recovering bacterial plasmids from transgenic eukaryotic genomic DNA. Total maize DNA was first digested with restriction enzyme(s), ligated, and then transformed into E. coli cells. Colonies were recovered under selection against the antibiotic marker(s) in the transgene vector. Ampicillin or carbenicillin was used for RescueMu transgene recovery. The flanking genomic sequences at RescueMu insertion sites were simultaneously captured and then sequenced using RescueMu-readout primers. Genomic DNA from an individual plant or from pooled samples of up to approximately 50 plants could be used in a single rescue. Because the majority of transgenic constructs currently used in flowering plants were made in the form of plasmids, this protocol could therefore be adapted by and useful to researchers involved in other transgenic work and be versatile for characterizing transgene loci.

Phenome analysis in plant species using loss-of-function and gain-of-function mutants

Analysis of genetic mutations is one of the most effective ways to investigate gene function. We now have methods that allow for mass production of mutant lines and cells in a variety of model species. Recently, large numbers of mutant lines have been generated by both 'loss-of-function' and 'gain-of-function' techniques. In parallel, phenotypic information covering various mutant resources has been acquired and released in web-based databases. As a result, significant progress in comprehensive phenotype analysis is being made through the use of these tools. Arabidopsis and rice are two major model plant species in which genome sequencing projects have been completed. Arabidopsis is the most widely used experimental plant, with a large number of mutant resources and several examples of systematic phenotype analysis. Rice is a major crop species and is used as a model plant, with an increasing number of mutant resources. Other plant species are also being employed in functional genetics research. In this review, the present status of mutant resources for large-scale studies of gene function in plant research and the current perspective on using loss-of-function and gain-of-function mutants in phenome research will be discussed.

Mutator transposon in maize and MULEs in the plant genome

Mutator (Mu) is by far the most mutagenic plant transposon. The high frequency of transposition and the tendency to insert into low copy sequences for such transposon have made it the primary means by which genes are mutagenized in maize (Zea mays L.). Mus like elements (MULEs) are widespread among angiosperms and multiple-diverged functional variants can be present in a single genome. MULEs often capture genetic sequences. These Pack-MuLEs can mobilize thousands of gene fragments, which may have had a significant impact on host genome evolution. There is also evidence that MULEs can move between reproductively isolated species. Here we present an overview of the discovery, features and utility of Mu transposon. Classification of Mu elements and future directions of related research are also discussed. Understanding Mu will help us elucidate the dynamic genome.

Nonradioactive genomic DNA blots for detection of low abundant sequences in transgenic maize

Sensitive and reproducible genotyping tools are fundamental in interpreting and substantiating genetic data. In cases where alternative assays like PCR are not applicable, a sensitive genomic Southern protocol is needed. Our maize gene discovery work using the RescueMu transgenic lines was such a task. The direct proof of each new germinal insertion event can be assessed only on a genomic DNA hybridization analysis, and therefore we developed the following protocol to screen efficiently through hundreds up to thousands of samples in a relatively short time. The DNA extraction protocol was scaled to accommodate samples processed in a microcentrifuge with consistent yield of approximately 50 microg of high molecular weight DNA. A trained person can easily process several hundred samples in a few days. Once the DNA is extracted, final results can be obtained routinely within a week on approximately 100 or more samples, depending on the capacity of the electrophoresis and hybridization apparatus available. Under our optimized conditions, the method described below generates blots with high sensitivity and low background even after repeated stripping and reprobing. Single to low-copy transgenes as well as maize genomic sequences can be detected consistently. The nonradioactive DNA probes employed are not only safer, compared to the conventional radioactive probes, but also greatly shorten the exposure time. Confident estimation of copy number - as good as quantitative PCR - and visualization of transgene complexity are just a few more advantages of this protocol.

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.

Testing the palindromic target site model for DNA transposon insertion using the Drosophila melanogaster P-element

Understanding the molecular mechanisms that influence transposable element target site preferences is a fundamental challenge in functional and evolutionary genomics. Large-scale transposon insertion projects provide excellent material to study target site preferences in the absence of confounding effects of post-insertion evolutionary change. Growing evidence from a wide variety of prokaryotes and eukaryotes indicates that DNA transposons recognize staggered-cut palindromic target site motifs (TSMs). Here, we use over 10 000 accurately mapped P-element insertions in the Drosophila melanogaster genome to test predictions of the staggered-cut palindromic target site model for DNA transposon insertion. We provide evidence that the P-element targets a 14-bp palindromic motif that can be identified at the primary sequence level, which predicts the local spacing, hotspots and strand orientation of P-element insertions. Intriguingly, we find that the although P-element destroys the complete 14-bp target site upon insertion, the terminal three nucleotides of the P-element inverted repeats complement and restore the original TSM, suggesting a mechanistic link between transposon target sites and their terminal inverted repeats. Finally, we discuss how the staggered-cut palindromic target site model can be used to assess the accuracy of genome mappings for annotated P-element insertions.

Maize and sorghum: genetic resources for bioenergy grasses

The highly photosynthetic-efficient C4 grasses, such as switchgrass (Panicum virgatum), Miscanthus (Miscanthusxgiganteus), sorghum (Sorghum bicolor) and maize (Zea mays), are expected to provide abundant and sustainable resources of lignocellulosic biomass for the production of biofuels. A deeper understanding of the synthesis, deposition and hydrolysis of the distinctive cell walls of grasses is crucial to gain genetic control of traits that contribute to biomass yield and quality. With a century of genetic investigations and breeding success, recently completed genome sequences, well-characterized cell wall compositions, and a close evolutionary relationship with future bioenergy perennial grasses, we propose that maize and sorghum are key model systems for gene discovery relating to biomass yield and quality in the bioenergy grasses.

The art and design of genetic screens: maize

Maize (Zea mays) is an excellent model for basic research. Genetic screens have informed our understanding of developmental processes, meiosis, epigenetics and biochemical pathways--not only in maize but also in other cereal crops. We discuss the forward and reverse genetic screens that are possible in this organism, and emphasize the available tools. Screens exploit the well-studied behaviour of transposon systems, and the distinctive chromosomes allow an integration of cytogenetics into mutagenesis screens and analyses. The imminent completion of the maize genome sequence provides the essential resource to move seamlessly from gene to phenotype and back.

A drastic reduction in DOF1 transcript levels does not affect C4-specific gene expression in maize

The transcription factor DOF1 has been suggested to regulate photosynthetic gene expression in maize. By screening a RescueMu transposon-tagged mutant library, we identified a maize mutant with a transposon integration in the Dof1 gene 16 bp upstream of the transcription initiation site (TIS). Sequencing of the Dof1 promoter region revealed an unusual promoter structure missing any typical elements. Homozygous (ho) mutant lines were generated by selfing and subsequent PCR and DNA gel blot analyses. The transposon integration reduced Dof1 transcript levels to less than 20% compared to the wild-type and overlapping RT-PCR systems revealed that these transcripts were not initiated from the native transcription start site. Dof1 transcripts transiently accumulate in wild-type plants after illumination of darkened seedlings, but this accumulation cannot be observed in mutant lines. However, the time-course of transcript accumulation from the C(4)-specific phosphoenolpyruvate carboxylase (PEPC) gene, a possible target of DOF1, is not altered. Moreover, no impact on the steady-state levels of five additional transcripts involved in C(4)-metabolism can be observed. The contents of amino acids, glucose, and malate as well as the carbon to nitrogen ratio in the leaves remained unchanged when comparing wild-type and mutant plants. Our data question the importance of DOF1 in the control of photosynthetic gene expression in maize.

Sequence-indexed mutations in maize using the UniformMu transposon-tagging population

BACKGROUND

Gene knockouts are a critical resource for functional genomics. In Arabidopsis, comprehensive knockout collections were generated by amplifying and sequencing genomic DNA flanking insertion mutants. These Flanking Sequence Tags (FSTs) map each mutant to a specific locus within the genome. In maize, FSTs have been generated using DNA transposons. Transposable elements can generate unstable insertions that are difficult to analyze for simple knockout phenotypes. Transposons can also generate somatic insertions that fail to segregate in subsequent generations.

RESULTS

Transposon insertion sites from 106 UniformMu FSTs were tested for inheritance by locus-specific PCR. We confirmed 89% of the FSTs to be germinal transposon insertions. We found no evidence for somatic insertions within the 11% of insertion sites that were not confirmed. Instead, this subset of insertion sites had errors in locus-specific primer design due to incomplete or low-quality genomic sequences. The locus-specific PCR assays identified a knockout of a 6-phosphogluconate dehydrogenase gene that co-segregates with a seed mutant phenotype. The mutant phenotype linked to this knockout generates novel hypotheses about the role for the plastid-localized oxidative pentose phosphate pathway during grain-fill.

CONCLUSION

We show that FSTs from the UniformMu population identify stable, germinal insertion sites in maize. Moreover, we show that these sequence-indexed mutations can be readily used for reverse genetic analysis. We conclude from these data that the current collection of 1,882 non-redundant insertion sites from UniformMu provide a genome-wide resource for reverse genetics.

Cytological visualization of DNA transposons and their transposition pattern in somatic cells of maize

Global genomic analysis of transposable element distributions of both natural (En/Spm, Ac-Ds, and MuDR/Mu) and modified (RescueMu) types was performed by fluorescence in situ hybridization (FISH) on somatic chromosomes coupled with karyotyping of each chromosome. In lines without an active transposable element, the locations of silent En/Spm, Ac-Ds, and MuDR/Mu elements were visualized, revealing variation in copy number and position among lines but no apparent locational bias. The ability to detect single elements was validated by using previously mapped active Ac elements. Somatic transpositions were documented in plants containing an engineered Mutator element, RescueMu, via use of the karyotyping system. By analyzing the RescueMu lines, we found that transposition of RescueMu in root-tip cells follows the cut-and-paste type of transposition. This work demonstrates the utility of FISH and karyotyping in the study of transposon activity and its consequences.

Differential methylation of genes and repeats in land plants

The hypomethylated fraction of plant genomes is usually enriched in genes and can be selectively cloned using methylation filtration (MF). Therefore, MF has been used as a gene enrichment technology in sorghum and maize, where gene enrichment was proportional to genome size. Here we apply MF to a broad variety of plant species spanning a wide range of genome sizes. Differential methylation of genic and non-genic sequences was observed in all species tested, from non-vascular to vascular plants, but in some cases, such as wheat and pine, a lower than expected level of enrichment was observed. Remarkably, hexaploid wheat and pine show a dramatically large number of gene-like sequences relative to other plants. In hexaploid wheat, this apparent excess of genes may reflect an abundance of methylated pseudogenes, which may thus be more prevalent in recent polyploids.

Comparative plant genomics resources at PlantGDB

PlantGDB (http://www.plantgdb.org/) is a database of plant molecular sequences. Expressed sequence tag (EST) sequences are assembled into contigs that represent tentative unique genes. EST contigs are functionally annotated with information derived from known protein sequences that are highly similar to the putative translation products. Tentative Gene Ontology terms are assigned to match those of the similar sequences identified. Genome survey sequences are assembled similarly. The resulting genome survey sequence contigs are matched to ESTs and conserved protein homologs to identify putative full-length open reading frame-containing genes, which are subsequently provisionally classified according to established gene family designations. For Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), the exon-intron boundaries for gene structures are annotated by spliced alignment of ESTs and full-length cDNAs to their respective complete genome sequences. Unique genome browsers have been developed to present all available EST and cDNA evidence for current transcript models (for Arabidopsis, see the AtGDB site at http://www.plantgdb.org/AtGDB/; for rice, see the OsGDB site at http://www.plantgdb.org/OsGDB/). In addition, a number of bioinformatic tools have been integrated at PlantGDB that enable researchers to carry out sequence analyses on-site using both their own data and data residing within the database.

Reduced representation sequencing: a success in maize and a promise for other plant genomes

Plant, and particularly cereal genomes, are challenging to sequence due to their large size and high repetitive DNA content. Gene-enrichment strategies are alternative or complementary approaches to complete genome sequencing that yield, rapidly and inexpensively, useful sequence data from large and complex genomes. The maize genome is large (2.7 Gbp) and contains large amounts of conserved repetitive elements. Furthermore, the high allelic diversity found between maize inbred lines may necessitate sequencing several inbred lines in order to recover the maize "gene pool". Two gene-enrichment approaches, methylation filtration (MF) and high C(o)t (HC) sequencing have been tested in maize and their ability to sample the gene space has been examined. Combined with other genomic sequencing strategies, gene-enriched genomic sequencing is a practical way to examine the maize gene pool, to order and orient the genic sequences on the genome, and to enable investigation of gene content of other complex plant genomes.