The unusual dRemp retrotransposon is abundant, highly mutagenic, and mobilized only in the second pollen mitosis of some maize lines.

The frequent mutations recovered recently from the pollen of select maize lines resulted from the meiotic mobilization of specific low-copy number long-terminal repeat (LTR) retrotransposons, which differ among lines. Mutations that arise at male meiosis produce kernels with concordant mutant phenotypes in both endosperm and embryo because the two sperms that participate in double fertilization are genetically identical. Those are in a majority. However, a small minority of kernels with a mutant endosperm carry a nonconcordant normal embryo, pointing to a postmeiotic or microgametophytic origin. In this study, we have identified the basis for those nonconcordant mutations. We find that all are produced by transposition of a defective LTR retrotransposon that we have termed dRemp (defective retroelement mobile in pollen). This element has several unique properties. Unlike the mutagenic LTR retrotransposons identified previously, dRemp is present in hundreds of copies in all sequenced lines. It seems to transpose only at the second pollen mitosis because all dRemp insertion mutants are nonconcordant yet recoverable in either the endosperm or the embryo. Although it does not move in most lines, dRemp is highly mobile in the Corn Belt inbred M14, identified earlier by breeders as being highly unstable. Lastly, it can be recovered in an array of structures, ranging from solo LTRs to tandem dRemp repeats containing several internal LTRs, suggestive of extensive recombination during retrotransposition. These results shed further light on the spontaneous mutation process and on the possible basis for inbred instability in maize.

Necrotic upper tips1 mimics heat and drought stress and encodes a protoxylem-specific transcription factor in maize.

Maintaining sufficient water transport during flowering is essential for proper organ growth, fertilization, and yield. Water deficits that coincide with flowering result in leaf wilting, necrosis, tassel browning, and sterility, a stress condition known as "tassel blasting." We identified a mutant, necrotic upper tips1 (nut1), that mimics tassel blasting and drought stress and reveals the genetic mechanisms underlying these processes. The nut1 phenotype is evident only after the floral transition, and the mutants have difficulty moving water as shown by dye uptake and movement assays. These defects are correlated with reduced protoxylem vessel thickness that indirectly affects metaxylem cell wall integrity and function in the mutant. nut1 is caused by an Ac transposon insertion into the coding region of a unique NAC transcription factor within the VND clade of Arabidopsis NUT1 localizes to the developing protoxylem of root, stem, and leaf sheath, but not metaxylem, and its expression is induced by flowering. NUT1 downstream target genes function in cell wall biosynthesis, apoptosis, and maintenance of xylem cell wall thickness and strength. These results show that maintaining protoxylem vessel integrity during periods of high water movement requires the expression of specialized, dynamically regulated transcription factors within the vasculature.

Spontaneous mutations in maize pollen are frequent in some lines and arise mainly from retrotranspositions and deletions.

While studying spontaneous mutations at the maize bronze (bz) locus, we made the unexpected discovery that specific low-copy number retrotransposons are mobile in the pollen of some maize lines, but not of others. We conducted large-scale genetic experiments to isolate new bz mutations from several Bz stocks and recovered spontaneous stable mutations only in the pollen parent in reciprocal crosses. Most of the new stable bz mutations resulted from either insertions of low-copy number long terminal repeat (LTR) retrotransposons or deletions, the same two classes of mutations that predominated in a collection of spontaneous wx mutations [Wessler S (1997) The Mutants of Maize, pp 385-386]. Similar mutations were recovered at the closely linked sh locus. These events occurred with a frequency of 2-4 × 10-5 in two lines derived from W22 and in 4Co63, but not at all in B73 or Mo17, two inbreds widely represented in Corn Belt hybrids. Surprisingly, the mutagenic LTR retrotransposons differed in the active lines, suggesting differences in the autonomous element make-up of the lines studied. Some active retrotransposons, like Hopscotch, Magellan, and Bs2, a Bs1 variant, were described previously; others, like Foto and Focou in 4Co63, were not. By high-throughput sequencing of retrotransposon junctions, we established that retrotranposition of Hopscotch, Magellan, and Bs2 occurs genome-wide in the pollen of active lines, but not in the female germline or in somatic tissues. We discuss here the implications of these results, which shed light on the source, frequency, and nature of spontaneous mutations in maize.

Maize defective kernel mutant generated by insertion of a Ds element in a gene encoding a highly conserved TTI2 cochaperone.

We have used the newly engineered transposable element Dsg to tag a gene that gives rise to a defective kernel (dek) phenotype. Dsg requires the autonomous element Ac for transposition. Upon excision, it leaves a short DNA footprint that can create in-frame and frameshift insertions in coding sequences. Therefore, we could create alleles of the tagged gene that confirmed causation of the dek phenotype by the Dsg insertion. The mutation, designated dek38-Dsg, is embryonic lethal, has a defective basal endosperm transfer (BETL) layer, and results in a smaller seed with highly underdeveloped endosperm. The maize dek38 gene encodes a TTI2 (Tel2-interacting protein 2) molecular cochaperone. In yeast and mammals, TTI2 associates with two other cochaperones, TEL2 (Telomere maintenance 2) and TTI1 (Tel2-interacting protein 1), to form the triple T complex that regulates DNA damage response. Therefore, we cloned the maize Tel2 and Tti1 homologs and showed that TEL2 can interact with both TTI1 and TTI2 in yeast two-hybrid assays. The three proteins regulate the cellular levels of phosphatidylinositol 3-kinase-related kinases (PIKKs) and localize to the cytoplasm and the nucleus, consistent with known subcellular locations of PIKKs. dek38-Dsg displays reduced pollen transmission, indicating TTI2's importance in male reproductive cell development.

Rolling-circle amplification of centromeric Helitrons in plant genomes.

The unusual eukaryotic Helitron transposons can readily capture host sequences and are, thus, evolutionarily important. They are presumed to amplify by rolling-circle replication (RCR) because some elements encode predicted proteins homologous to RCR prokaryotic transposases. In support of this replication mechanism, it was recently shown that transposition of a bat Helitron generates covalently closed circular intermediates. Another strong prediction is that RCR should generate tandem Helitron concatemers, yet almost all Helitrons identified to date occur as solo elements in the genome. To investigate alternative modes of Helitron organization in present-day genomes, we have applied the novel computational tool HelitronScanner to 27 plant genomes and have uncovered numerous tandem arrays of partially decayed, truncated Helitrons in all of them. Strikingly, most of these Helitron tandem arrays are interspersed with other repeats in centromeres. Many of these arrays have multiple Helitron 5' ends, but a single 3' end. The number of repeats in any one array can range from a handful to several hundreds. We propose here an RCR model that conforms to the present Helitron landscape of plant genomes. Our study provides strong evidence that plant Helitrons amplify by RCR and that the tandemly arrayed replication products accumulate mostly in centromeres.

Polarized gene conversion at the bz locus of maize.

Nucleotide diversity is greater in maize than in most organisms studied to date, so allelic pairs in a hybrid tend to be highly polymorphic. Most recombination events between such pairs of maize polymorphic alleles are crossovers. However, intragenic recombination events not associated with flanking marker exchange, corresponding to noncrossover gene conversions, predominate between alleles derived from the same progenitor. In these dimorphic heterozygotes, the two alleles differ only at the two mutant sites between which recombination is being measured. To investigate whether gene conversion at the bz locus is polarized, two large diallel crossing matrices involving mutant sites spread across the bz gene were performed and more than 2,500 intragenic recombinants were scored. In both diallels, around 90% of recombinants could be accounted for by gene conversion. Furthermore, conversion exhibited a striking polarity, with sites located within 150 bp of the start and stop codons converting more frequently than sites located in the middle of the gene. The implications of these findings are discussed with reference to recent data from genome-wide studies in other plants.

HelitronScanner uncovers a large overlooked cache of Helitron transposons in many plant genomes.

Transposons make up the bulk of eukaryotic genomes, but are difficult to annotate because they evolve rapidly. Most of the unannotated portion of sequenced genomes is probably made up of various divergent transposons that have yet to be categorized. Helitrons are unusual rolling circle eukaryotic transposons that often capture gene sequences, making them of considerable evolutionary importance. Unlike other DNA transposons, Helitrons do not end in inverted repeats or create target site duplications, so they are particularly challenging to identify. Here we present HelitronScanner, a two-layered local combinational variable (LCV) tool for generalized Helitron identification that represents a major improvement over previous identification programs based on DNA sequence or structure. HelitronScanner identified 64,654 Helitrons from a wide range of plant genomes in a highly automated way. We tested HelitronScanner's predictive ability in maize, a species with highly heterogeneous Helitron elements. LCV scores for the 5' and 3' termini of the predicted Helitrons provide a primary confidence level and element copy number provides a secondary one. Newly identified Helitrons were validated by PCR assays or by in silico comparative analysis of insertion site polymorphism among multiple accessions. Many new Helitrons were identified in model species, such as maize, rice, and Arabidopsis, and in a variety of organisms where Helitrons had not been reported previously to our knowledge, leading to a major upward reassessment of their abundance in plant genomes. HelitronScanner promises to be a valuable tool in future comparative and evolutionary studies of this major transposon superfamily.

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.

Variation in allelic expression associated with a recombination hotspot in Zea mays.

Gene expression is a complex process, requiring precise spatial and temporal regulation of transcription factor activity; however, modifications of individual cis- and trans-acting modules can be molded by natural selection to create a sizeable number of novel phenotypes. Results from decades of research indicate that developmental and phenotypic divergence among eukaryotic organisms is driven primarily by variation in levels of gene expression that are dictated by mutations, either in structural or regulatory regions, of genes. The relative contributions and interplay of cis- and trans-acting regulatory factors to this evolutionary process, however, remain poorly understood. Analysis of eight genes in the Bz1-Sh1 interval of Zea mays (maize) indicates significant allele-specific expression biases in at least one tissue for all genes, ranging from 1.3-fold to 36-fold. All detected effects were cis-regulatory in nature, although genetic background may also influence the level of expression bias and tissue specificity for some allelic combinations. Most allelic pairs exhibited the same direction and approximate intensity of bias across all four tissues; however, a subset of allelic pairs show alternating dominance across different tissue types or variation in the degree of bias in different tissues. In addition, the genes showing the most striking levels of allelic bias co-localize with a previously described recombination hotspot in this region, suggesting a naturally occurring genetic mechanism for creating regulatory variability for a subset of plant genes that may ultimately lead to evolutionary diversification.

The spectrum and frequency of self-inflicted and host gene mutations produced by the transposon Ac in maize.

The autonomous transposon Activator (Ac) is a powerful mutagen. Ac-induced mutations range from small footprints of host sequences to large rearrangements of transposon or host sequences. These mutations arise by different repair mechanisms of the double-strand break produced by Ac excision: footprints by nonhomologous end joining and rearrangements by various mechanisms, including DNA replication repair. Footprints greatly outnumber other mutations, masking them because they usually share a nonfunctional phenotype. To determine the spectrum and frequencies of host and self-mutations generated by Ac, we used an allele harboring Ac in the 5' untranslated region bronze (bz). In this system, simple excisions produce purple revertants, whereas deletions of host or transposon sequences produce stable bronze (bz-s) mutants. Internal and terminal deletions of Ac predominated among the 72 bz-s derivatives. Most internal deletions (52 of 54) behaved as nonautonomous Dissociation (Ds) elements. All nine terminal deletions or fractured Ac (fAc) elements had rearrangements of adjacent host sequences. Most Ds and fAc deletion junctions displayed microhomologies and contained filler DNA from nearby sequences, suggesting an origin by DNA repair synthesis followed by microhomology-mediated end joining. All mutations occurred more frequently in pollen, where one in 200 grains carried new Ds or fAc elements.

Dynamic evolution of bz orthologous regions in the Andropogoneae and other grasses.

Genome structure exhibits remarkable plasticity within Zea mays. To examine how haplotype structure has evolved within the Andropogoneae tribe, we have analyzed the bz gene-rich region of maize (Zea mays), the Zea teosintes mays ssp. mexicana, luxurians and diploperennis, Tripsacum dactyloides, Coix lacryma-jobi and Sorghum propinquum. We sequenced and annotated BAC clones from these species and re-annotated the orthologous Sorghum bicolor region. Gene colinearity in the region is well conserved within the genus Zea. However, the orthologous regions of Coix and Sorghum exhibited several micro-rearrangements relative to Zea, including addition, truncation and deletion of genes. The stc1 gene, involved in the production of a terpenoid insect defense signal, is evolving particularly fast, and its progressive disappearance from some species is occurring by microhomology-mediated recombination. LTR retrotransposons are the main contributors to the dynamic evolution of the bz region. Common transposon insertion sites occur among haplotypes from different Zea mays sub-species, but not outside the species. As in Zea, different patterns of interspersion between genes and retrotransposons are observed in Sorghum. We estimate that the mean divergence times between maize and Tripsacum, Coix and Sorghum are 8.5, 12.1 and 12.4 million years ago, respectively, and that between Coix and Sorghum is 9.3 million years ago. A comparison of the bz orthologous regions of Zea, Sorghum and Coix with those of Brachypodium, Setaria and Oryza allows us to infer how the region has evolved by addition and deletion of genes in the approximately 50 million years since these genera diverged from a common progenitor.

InsertionMapper: a pipeline tool for the identification of targeted sequences from multidimensional high throughput sequencing data.

BACKGROUND: The advent of next-generation high-throughput technologies has revolutionized whole genome sequencing, yet some experiments require sequencing only of targeted regions of the genome from a very large number of samples. These regions can be amplified by PCR and sequenced by next-generation methods using a multidimensional pooling strategy. However, there is at present no available generalized tool for the computational analysis of target-enriched NGS data from multidimensional pools. RESULTS: Here we present InsertionMapper, a pipeline tool for the identification of targeted sequences from multidimensional high throughput sequencing data. InsertionMapper consists of four independently working modules: Data Preprocessing, Database Modeling, Dimension Deconvolution and Element Mapping. We illustrate InsertionMapper with an example from our project 'New reverse genetics resources for maize', which aims to sequence-index a collection of 15,000 independent insertion sites of the transposon Ds in maize. Identified sequences are validated by PCR assays. This pipeline tool is applicable to similar scenarios requiring analysis of the tremendous output of short reads produced in NGS sequencing experiments of targeted genome sequences. CONCLUSIONS: InsertionMapper is proven efficacious for the identification of target-enriched sequences from multidimensional high throughput sequencing data. With adjustable parameters and experiment configurations, this tool can save great computational effort to biologists interested in identifying their sequences of interest within the huge output of modern DNA sequencers. InsertionMapper is freely accessible at https://sourceforge.net/p/insertionmapper and http://bo.csam.montclair.edu/du/insertionmapper.

Haplotype structure strongly affects recombination in a maize genetic interval polymorphic for Helitron and retrotransposon insertions.

We have asked here how the remarkable variation in maize haplotype structure affects recombination. We compared recombination across a genetic interval of 9S in 2 highly dissimilar heterozygotes that shared 1 parent. The genetic interval in the common haplotype is approximately 100 kb long and contains 6 genes interspersed with gene-fragment-bearing Helitrons and retrotransposons that, together, comprise 70% of its length. In one heterozygote, most intergenic insertions are homozygous, although polymorphic, enabling us to determine whether any recombination junctions fall within them. In the other, most intergenic insertions are hemizygous and, thus, incapable of homologous recombination. Our analysis of the frequency and distribution of recombination in the interval revealed that: (i) Most junctions were circumscribed to the gene space, where they showed a highly nonuniform distribution. In both heterozygotes, more than half of the junctions fell in the stc1 gene, making it a clear recombination hotspot in the region. However, the genetic size of stc1 was 2-fold lower when flanked by a hemizygous 25-kb retrotransposon cluster. (ii) No junctions fell in the hypro1 gene in either heterozygote, making it a genic recombination coldspot. (iii) No recombination occurred within the gene fragments borne on Helitrons nor within retrotransposons, so neither insertion class contributes to the interval's genetic length. (iv) Unexpectedly, several junctions fell in an intergenic region not shared by all 3 haplotypes. (v) In general, the ability of a sequence to recombine correlated inversely with its methylation status. Our results show that haplotypic structural variability strongly affects the frequency and distribution of recombination events in maize.

The polychromatic Helitron landscape of the maize genome.

Maize Helitron transposons are intriguing because of their notable ability to capture gene fragments and move them around the genome. To document more extensively their variability and their contribution to the remarkable genome structure variation of present-day maize, we have analyzed their composition, copy number, timing of insertion, and chromosomal distribution. First, we searched 2.4 Gb of sequences generated by the Maize Genome Sequencing Project with our HelitronFinder program. We identified 2,791 putative nonautonomous Helitrons and manually curated a subset of 272. The predicted Helitrons measure 11.9 kb on average and carry from zero to nine gene fragments, captured from 376 different genes. Although the diversity of Helitron gene fragments in maize is greater than in other species, more than one-third of annotated Helitrons carry fragments derived from just one of two genes. Most members in these two subfamilies inserted in the genome less than one million years ago. Second, we conducted a BLASTN search of the maize sequence database with queries from two previously described agenic Helitrons not detected by HelitronFinder. Two large subfamilies of Helitrons or Helitron-related transposons were identified. One subfamily, termed Cornucopious, consists of thousands of copies of an approximately 1.0-kb agenic Helitron that may be the most abundant transposon in maize. The second subfamily consists of >150 copies of a transposon-like sequence, termed Heltir, that has terminal inverted repeats resembling Helitron 3' termini. Nonautonomous Helitrons make up at least 2% of the maize genome and most of those tested show +/- polymorphisms among modern inbred lines.

The complete Ac/Ds transposon family of maize.

BACKGROUND: The nonautonomous maize Ds transposons can only move in the presence of the autonomous element Ac. They comprise a heterogeneous group that share 11-bp terminal inverted repeats (TIRs) and some subterminal repeats, but vary greatly in size and composition. Three classes of Ds elements can cause mutations: Ds-del, internal deletions of the 4.6-kb Ac element; Ds1, ~400-bp in size and sharing little homology with Ac, and Ds2, variably-sized elements containing about 0.5 kb from the Ac termini and unrelated internal sequences. Here, we analyze the entire complement of Ds-related sequences in the genome of the inbred B73 and ask whether additional classes of Ds-like (Ds-l) elements, not uncovered genetically, are mobilized by Ac. We also compare the makeup of Ds-related sequences in two maize inbreds of different origin. RESULTS: We found 903 elements with 11-bp Ac/Ds TIRs flanked by 8-bp target site duplications. Three resemble Ac, but carry small rearrangements. The others are much shorter, once extraneous insertions are removed. There are 331 Ds1 and 39 Ds2 elements, many of which are likely mobilized by Ac, and two novel classes of Ds-l elements. Ds-l3 elements lack subterminal homology with Ac, but carry transposase gene fragments, and represent decaying Ac elements. There are 44 such elements in B73. Ds-l4 elements share little similarity with Ac outside of the 11-bp TIR, have a modal length of ~1 kb, and carry filler DNA which, in a few cases, could be matched to gene fragments. Most Ds-related elements in B73 (486/903) fall in this class. None of the Ds-l elements tested responded to Ac. Only half of Ds insertion sites examined are shared between the inbreds B73 and W22. CONCLUSIONS: The majority of Ds-related sequences in maize correspond to Ds-l elements that do not transpose in the presence of Ac. Unlike actively transposing elements, many Ds-l elements are inserted in repetitive DNA, where they probably become methylated and begin to decay. The filler DNA present in most elements is occasionally captured from genes, a rare feature in transposons of the hAT superfamily to which Ds belongs. Maize inbreds of different origin are highly polymorphic in their DNA transposon makeup.

Give-and-take: interactions between DNA transposons and their host plant genomes.

Recent genome sequencing efforts have revealed how extensively transposable elements (TEs) have contributed to the shaping of present day plant genomes. DNA transposons associate preferentially with the euchromatic or genic component of plant genomes and have had the opportunity to interact intimately with the genes of the plant host. These interactions have resulted in TEs acquiring host sequences, forming chimeric genes through exon shuffling, replacing regulatory sequences, mobilizing genes around the genome, and contributing genes to the host. The close interaction of transposons with genes has also led to the evolution of intricate cellular mechanisms for silencing transposon activity. Transposons have thus become important subjects of study in understanding epigenetic regulation and, in cases where transposons have amplified to high numbers, how to escape that regulation.

Excision of Helitron transposons in maize.

Helitrons are novel transposons discovered by bioinformatic analysis of eukaryotic genome sequences. They are believed to move by rolling circle (RC) replication because their predicted transposases are homologous to those of bacterial RC transposons. We report here evidence of somatic Helitron excision in maize, an unexpected finding suggesting that Helitrons can exhibit an excisive mode of transposition.

Organization and variability of the maize genome.

With a size approximating that of the human genome, the maize genome is about to become the largest plant genome yet sequenced. Contributing to that size are a whole-genome duplication event and a retrotransposition explosion that produced a large amount of repetitive DNA. This DNA is greatly under-represented in cDNA collections, so analysis of the maize transcriptome has been an expedient way of assessing the gene content of maize. Over 2 million maize cDNA sequences are now available, making maize the third most widely studied organism, behind mouse and man. To date, the sequencing of large-sized DNA clones has been largely driven by the genetic interests of different investigators. The recent construction of a physical map that is anchored to the genetic map will aid immensely in the maize genome-sequencing effort. However, studies showing that the repetitive DNA component is highly polymorphic among maize inbred lines point to the need to sample vertically a few specific regions of the genome to evaluate the extent and importance of this variability.

Computational prediction and molecular confirmation of Helitron transposons in the maize genome.

BACKGROUND: Helitrons represent a new class of transposable elements recently uncovered in plants and animals. One remarkable feature of Helitrons is their ability to capture gene sequences, which makes them of considerable potential evolutionary importance. However, because Helitrons lack the typical structural features of other DNA transposable elements, identifying them is a challenge. Currently, most researchers identify Helitrons manually by comparing sequences. With the maize whole genome sequencing project underway, an automated computational Helitron searching tool is needed. The characterization of Helitron activities in maize needs to be addressed in order to better understand the impact of Helitrons on the organization of the genome. RESULTS: We developed and implemented a heuristic searching algorithm in PERL for identifying Helitrons. Our HelitronFinder program will (i) take FASTA-formatted DNA sequences as input and identify the hairpin looping patterns, and (ii) exploit the consensus 5' and 3' end sequences of known Helitrons to identify putative ends. We randomly selected five predicted Helitrons from the program's high quality output for molecular verification. Four out of the five predicted Helitrons were confirmed by PCR assays and DNA sequencing in different maize inbred lines. The HelitronFinder program identified two head-to-head dissimilar Helitrons in a maize BAC sequence. CONCLUSION: We have identified 140 new Helitron candidates in maize with our computational tool HelitronFinder by searching maize DNA sequences currently available in GenBank. Four out of five candidates were confirmed to be real by empirical methods, thus validating the predictions of HelitronFinder. Additional points to emerge from our study are that Helitrons do not always insert at an AT dinucleotide in the host sequences, that they can insert immediately adjacent to an existing Helitron, and that their movement may cause changes in the flanking region, such as deletions.

The maize W22 genome provides a foundation for functional genomics and transposon biology.

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.