Sanzhen Liu

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.

Heritable Epigenetic Variation among Maize Inbreds

Epigenetic variation describes heritable differences that are not attributable to changes in DNA sequence. There is the potential for pure epigenetic variation that occurs in the absence of any genetic change or for more complex situations that involve both genetic and epigenetic differences. Methylation of cytosine residues provides one mechanism for the inheritance of epigenetic information. A genome-wide profiling of DNA methylation in two different genotypes of Zea mays (ssp. mays), an organism with a complex genome of interspersed genes and repetitive elements, allowed the identification and characterization of examples of natural epigenetic variation. The distribution of DNA methylation was profiled using immunoprecipitation of methylated DNA followed by hybridization to a high-density tiling microarray. The comparison of the DNA methylation levels in the two genotypes, B73 and Mo17, allowed for the identification of approximately 700 differentially methylated regions (DMRs). Several of these DMRs occur in genomic regions that are apparently identical by descent in B73 and Mo17 suggesting that they may be examples of pure epigenetic variation. The methylation levels of the DMRs were further studied in a panel of near-isogenic lines to evaluate the stable inheritance of the methylation levels and to assess the contribution of cis- and trans- acting information to natural epigenetic variation. The majority of DMRs that occur in genomic regions without genetic variation are controlled by cis-acting differences and exhibit relatively stable inheritance. This study provides evidence for naturally occurring epigenetic variation in maize, including examples of pure epigenetic variation that is not conditioned by genetic differences. The epigenetic differences are variable within maize populations and exhibit relatively stable trans-generational inheritance. The detected examples of epigenetic variation, including some without tightly linked genetic variation, may contribute to complex trait variation.

Gene Mapping via Bulked Segregant RNA-Seq (BSR-Seq)

Bulked segregant analysis (BSA) is an efficient method to rapidly and efficiently map genes responsible for mutant phenotypes. BSA requires access to quantitative genetic markers that are polymorphic in the mapping population. We have developed a modification of BSA (BSR-Seq) that makes use of RNA-Seq reads to efficiently map genes even in populations for which no polymorphic markers have been previously identified. Because of the digital nature of next-generation sequencing (NGS) data, it is possible to conduct de novo SNP discovery and quantitatively genotype BSA samples by analyzing the same RNA-Seq data using an empirical Bayesian approach. In addition, analysis of the RNA-Seq data provides information on the effects of the mutant on global patterns of gene expression at no extra cost. In combination these results greatly simplify gene cloning experiments. To demonstrate the utility of this strategy BSR-Seq was used to clone the glossy3 (gl3) gene of maize. Mutants of the glossy loci exhibit altered accumulation of epicuticular waxes on juvenile leaves. By subjecting the reference allele of gl3 to BSR-Seq, we were able to map the gl3 locus to an ∼2 Mb interval. The single gene located in the ∼2 Mb mapping interval whose expression was down-regulated in the mutant pool was subsequently demonstrated to be the gl3 gene via the analysis of multiple independent transposon induced mutant alleles. The gl3 gene encodes a putative myb transcription factor, which directly or indirectly affects the expression of a number of genes involved in the biosynthesis of very-long-chain fatty acids.

High-Throughput Genetic Mapping of Mutants via Quantitative Single Nucleotide Polymorphism Typing

Advances in next-generation sequencing technology have facilitated the discovery of single nucleotide polymorphisms (SNPs). Sequenom-based SNP-typing assays were developed for 1359 maize SNPs identified via comparative next-generation transcriptomic sequencing. Approximately 75% of these SNPs were successfully converted into genetic markers that can be scored reliably and used to generate a SNP-based genetic map by genotyping recombinant inbred lines from the intermated B73 × Mo17 population. The quantitative nature of Sequenom-based SNP assays led to the development of a time- and cost-efficient strategy to genetically map mutants via quantitative bulked segregant analysis. This strategy was used to rapidly map the loci associated with several dozen recessive mutants. Because a mutant can be mapped using as few as eight multiplexed sets of SNP assays on a bulk of as few as 20 mutant F2 individuals, this strategy is expected to be widely adopted for mapping in many species.

B73-Mo17 Near-Isogenic Lines Demonstrate Dispersed Structural Variation in Maize

Recombinant inbred lines developed from the maize (Zea mays ssp. mays) inbreds B73 and Mo17 have been widely used to discover quantitative trait loci controlling a wide variety of phenotypic traits and as a resource to produce high-resolution genetic maps. These two parents were used to produce a set of near-isogenic lines (NILs) with small regions of introgression into both backgrounds. A novel array-based genotyping platform was used to score genotypes of over 7,000 loci in 100 NILs with B73 as the recurrent parent and 50 NILs with Mo17 as the recurrent parent. This population contains introgressions that cover the majority of the maize genome. The set of NILs displayed an excess of residual heterozygosity relative to the amount expected based on their pedigrees, and this excess residual heterozygosity is enriched in the low-recombination regions near the centromeres. The genotyping platform provided the ability to survey copy number variants that exist in more copies in Mo17 than in B73. The majority of these Mo17-specific duplications are located in unlinked positions throughout the genome. The utility of this population for the discovery and validation of quantitative trait loci was assessed through analysis of plant height variation.

Mendelian and Non-Mendelian Regulation of Gene Expression in Maize

Transcriptome variation plays an important role in affecting the phenotype of an organism. However, an understanding of the underlying mechanisms regulating transcriptome variation in segregating populations is still largely unknown. We sought to assess and map variation in transcript abundance in maize shoot apices in the intermated B73×Mo17 recombinant inbred line population. RNA–based sequencing (RNA–seq) allowed for the detection and quantification of the transcript abundance derived from 28,603 genes. For a majority of these genes, the population mean, coefficient of variation, and segregation patterns could be predicted by the parental expression levels. Expression quantitative trait loci (eQTL) mapping identified 30,774 eQTL including 96 trans-eQTL “hotspots,” each of which regulates the expression of a large number of genes. Interestingly, genes regulated by a trans-eQTL hotspot tend to be enriched for a specific function or act in the same genetic pathway. Also, genomic structural variation appeared to contribute to cis-regulation of gene expression. Besides genes showing Mendelian inheritance in the RIL population, we also found genes whose expression level and variation in the progeny could not be predicted based on parental difference, indicating that non-Mendelian factors also contribute to expression variation. Specifically, we found 145 genes that show patterns of expression reminiscent of paramutation such that all the progeny had expression levels similar to one of the two parents. Furthermore, we identified another 210 genes that exhibited unexpected patterns of transcript presence/absence. Many of these genes are likely to be gene fragments resulting from transposition, and the presence/absence of their transcripts could influence expression levels of their ancestral syntenic genes. Overall, our results contribute to the identification of novel expression patterns and broaden the understanding of transcriptional variation in plants.

The maize brown midrib2 (bm2) gene encodes a methylenetetrahydrofolate reductase that contributes to lignin accumulation

The midribs of maize brown midrib (bm) mutants exhibit a reddish-brown color associated with reductions in lignin concentration and alterations in lignin composition. Here, we report the mapping, cloning, and functional and biochemical analyses of the bm2 gene. The bm2 gene was mapped to a small region of chromosome 1 that contains a putative methylenetetrahydrofolate reductase (MTHFR) gene, which is down-regulated in bm2 mutant plants. Analyses of multiple Mu-induced bm2-Mu mutant alleles confirmed that this constitutively expressed gene is bm2. Yeast complementation experiments and a previously published biochemical characterization show that the bm2 gene encodes a functional MTHFR. Quantitative RT-PCR analyses demonstrated that the bm2 mutants accumulate substantially reduced levels of bm2 transcript. Alteration of MTHFR function is expected to influence accumulation of the methyl donor S-adenosyl-l-methionine (SAM). Because SAM is consumed by two methyltransferases in the lignin pathway (Ye et al., 1994), the finding that bm2 encodes a functional MTHFR is consistent with its lignin phenotype. Consistent with this functional assignment of bm2, the expression patterns of genes in a variety of SAM-dependent or -related pathways, including lignin biosynthesis, are altered in the bm2 mutant. Biochemical assays confirmed that bm2 mutants accumulate reduced levels of lignin with altered composition compared to wild-type. Hence, this study demonstrates a role for MTHFR in lignin biosynthesis.

Histone lysine methyltransferase SDG8 is involved in brassinosteroid-regulated gene expression in Arabidopsis thaliana.

The plant steroid hormones, brassinosteroids (BRs), play important roles in plant growth, development, and responses to environmental stresses. BRs signal through receptors localized to the plasma membrane and other signaling components to regulate the BES1/BZR1 family of transcription factors, which modulates the expression of thousands of genes. How BES1/BZR1 and their interacting proteins function to regulate the large number of genes are not completely understood. Here we report that histone lysine methyltransferase SDG8, implicated in histone 3 lysine 36 di- and trimethylation (H3K36me2 and me3), is involved in BR-regulated gene expression. BES1 interacts with SDG8, directly or indirectly through IWS1, a transcription elongation factor involved in BR-regulated gene expression. The knockout mutant sdg8 displays a reduced growth phenotype with compromised BR responses. Global gene expression studies demonstrated that, while BR regulates about 5000 genes in wild-type plants, the hormone regulates fewer than 700 genes in sdg8 mutant. In addition, more than half of BR-regulated genes are differentially affected in sdg8 mutant. A Chromatin Immunoprecipitation (ChIP) experiment showed that H3K36me3 is reduced in BR-regulated genes in the sdg8 mutant. Based on these results, we propose that SDG8 plays an essential role in mediating BR-regulated gene expression. Our results thus reveal a major mechanism by which histone modifications dictate hormonal regulation of gene expression.

The Maize glossy13 Gene, Cloned via BSR-Seq and Seq-Walking Encodes a Putative ABC Transporter Required for the Normal Accumulation of Epicuticular Waxes

Aerial plant surfaces are covered by epicuticular waxes that among other purposes serve to control water loss. Maize glossy mutants originally identified by their “glossy” phenotypes exhibit alterations in the accumulation of epicuticular waxes. By combining data from a BSR-Seq experiment and the newly developed Seq-Walking technology, GRMZM2G118243 was identified as a strong candidate for being the glossy13 gene. The finding that multiple EMS-induced alleles contain premature stop codons in GRMZM2G118243, and the one knockout allele of gl13, validates the hypothesis that gene GRMZM2G118243 is gl13. Consistent with this, GRMZM2G118243 is an ortholog of AtABCG32 (Arabidopsis thaliana), HvABCG31 (barley) and OsABCG31 (rice), which encode ABCG subfamily transporters involved in the trans-membrane transport of various secondary metabolites. We therefore hypothesize that gl13 is involved in the transport of epicuticular waxes onto the surfaces of seedling leaves.

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.