Cheng-Ting Yeh

Maize inbreds exhibit high levels of copy number variation (CNV) and presence/absence variation (PAV) in genome content

Following the domestication of maize over the past ∼10,000 years, breeders have exploited the extensive genetic diversity of this species to mold its phenotype to meet human needs. The extent of structural variation, including copy number variation (CNV) and presence/absence variation (PAV), which are thought to contribute to the extraordinary phenotypic diversity and plasticity of this important crop, have not been elucidated. Whole-genome, array-based, comparative genomic hybridization (CGH) revealed a level of structural diversity between the inbred lines B73 and Mo17 that is unprecedented among higher eukaryotes. A detailed analysis of altered segments of DNA conservatively estimates that there are several hundred CNV sequences among the two genotypes, as well as several thousand PAV sequences that are present in B73 but not Mo17. Haplotype-specific PAVs contain hundreds of single-copy, expressed genes that may contribute to heterosis and to the extraordinary phenotypic diversity of this important crop.

Parallel domestication of the Shattering1 genes in cereals

A key step during crop domestication is the loss of seed shattering. Here, we show that seed shattering in sorghum is controlled by a single gene, Shattering1 (Sh1), which encodes a YABBY transcription factor. Domesticated sorghums harbor three different mutations at the Sh1 locus. Variants at regulatory sites in the promoter and intronic regions lead to a low level of expression, a 2.2-kb deletion causes a truncated transcript that lacks exons 2 and 3, and a GT-to-GG splice-site variant in the intron 4 results in removal of the exon 4. The distributions of these non-shattering haplotypes among sorghum landraces suggest three independent origins. The function of the rice ortholog (OsSh1) was subsequently validated with a shattering-resistant mutant, and two maize orthologs (ZmSh1-1 and ZmSh1-5.1+ZmSh1-5.2) were verified with a large mapping population. Our results indicate that Sh1 genes for seed shattering were under parallel selection during sorghum, rice and maize domestication.

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.

Parent-of-Origin Effects on Gene Expression and DNA Methylation in the Maize Endosperm

This work uses deep sequencing of RNAs from maize endosperm and embryo to identify 54 maternally expressed genes and 46 paternally expressed genes and then examines genome-wide DNA methylation and gene expression, finding hypomethylation of the maternal allele and endosperm-specific expression for many of the imprinted genes. Imprinting describes the differential expression of alleles based on their parent of origin. Deep sequencing of RNAs from maize (Zea mays) endosperm and embryo tissue 14 d after pollination was used to identify imprinted genes among a set of ~12,000 genes that were expressed and contained sequence polymorphisms between the B73 and Mo17 genotypes. The analysis of parent-of-origin patterns of expression resulted in the identification of 100 putative imprinted genes in maize endosperm, including 54 maternally expressed genes (MEGs) and 46 paternally expressed genes (PEGs). Three of these genes have been previously identified as imprinted, while the remaining 97 genes represent novel imprinted maize genes. A genome-wide analysis of DNA methylation identified regions with reduced endosperm DNA methylation in, or near, 19 of the 100 imprinted genes. The reduced levels of DNA methylation in endosperm are caused by hypomethylation of the maternal allele for both MEGs and PEGs in all cases tested. Many of the imprinted genes with reduced DNA methylation levels also show endosperm-specific expression patterns. The imprinted maize genes were compared with imprinted genes identified in genome-wide screens of rice (Oryza sativa) and Arabidopsis thaliana, and at least 10 examples of conserved imprinting between maize and each of the other species were identified.

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.

Spreading of Heterochromatin Is Limited to Specific Families of Maize Retrotransposons

Transposable elements (TEs) have the potential to act as controlling elements to influence the expression of genes and are often subject to heterochromatic silencing. The current paradigm suggests that heterochromatic silencing can spread beyond the borders of TEs and influence the chromatin state of neighboring low-copy sequences. This would allow TEs to condition obligatory or facilitated epialleles and act as controlling elements. The maize genome contains numerous families of class I TEs (retrotransposons) that are present in moderate to high copy numbers, and many are found in regions near genes, which provides an opportunity to test whether the spreading of heterochromatin from retrotransposons is prevalent. We have investigated the extent of heterochromatin spreading into DNA flanking each family of retrotransposons by profiling DNA methylation and di-methylation of lysine 9 of histone 3 (H3K9me2) in low-copy regions of the maize genome. The effects of different retrotransposon families on local chromatin are highly variable. Some retrotransposon families exhibit enrichment of heterochromatic marks within 800–1,200 base pairs of insertion sites, while other families exhibit very little evidence for the spreading of heterochromatic marks. The analysis of chromatin state in genotypes that lack specific insertions suggests that the heterochromatin in low-copy DNA flanking retrotransposons often results from the spreading of silencing marks rather than insertion-site preferences. Genes located near TEs that exhibit spreading of heterochromatin tend to be expressed at lower levels than other genes. Our findings suggest that a subset of retrotransposon families may act as controlling elements influencing neighboring sequences, while the majority of retrotransposons have little effect on flanking sequences.

Genic and nongenic contributions to natural variation of quantitative traits in maize

The complex genomes of many economically important crops present tremendous challenges to understand the genetic control of many quantitative traits with great importance in crop production, adaptation, and evolution. Advances in genomic technology need to be integrated with strategic genetic design and novel perspectives to break new ground. Complementary to individual-gene-targeted research, which remains challenging, a global assessment of the genomic distribution of trait-associated SNPs (TASs) discovered from genome scans of quantitative traits can provide insights into the genetic architecture and contribute to the design of future studies. Here we report the first systematic tabulation of the relative contribution of different genomic regions to quantitative trait variation in maize. We found that TASs were enriched in the nongenic regions, particularly within a 5-kb window upstream of genes, which highlights the importance of polymorphisms regulating gene expression in shaping the natural variation. Consistent with these findings, TASs collectively explained 44%-59% of the total phenotypic variation across maize quantitative traits, and on average, 79% of the explained variation could be attributed to TASs located in genes or within 5 kb upstream of genes, which together comprise only 13% of the genome. Our findings suggest that efficient, cost-effective genome-wide association studies (GWAS) in species with complex genomes can focus on genic and promoter regions.

Repeat subtraction-mediated sequence capture from a complex genome

Sequence capture technologies, pioneered in mammalian genomes, enable the resequencing of targeted genomic regions. Most capture protocols require blocking DNA, the production of which in large quantities can prove challenging. A blocker-free, two-stage capture protocol was developed using NimbleGen arrays. The first capture depletes the library of repetitive sequences, while the second enriches for target loci. This strategy was used to resequence non-repetitive portions of an approximately 2.2 Mb chromosomal interval and a set of 43 genes dispersed in the 2.3 Gb maize genome. This approach achieved approximately 1800-3000-fold enrichment and 80-98% coverage of targeted bases. More than 2500 SNPs were identified in target genes. Low rates of false-positive SNP predictions were obtained, even in the presence of captured paralogous sequences. Importantly, it was possible to recover novel sequences from non-reference alleles. The ability to design novel repeat-subtraction and target capture arrays makes this technology accessible in any species.

Complementation contributes to transcriptome complexity in maize (Zea mays L.) hybrids relative to their inbred parents

Typically, F(1)-hybrids are more vigorous than their homozygous, genetically distinct parents, a phenomenon known as heterosis. In the present study, the transcriptomes of the reciprocal maize (Zea mays L.) hybrids B73×Mo17 and Mo17×B73 and their parental inbred lines B73 and Mo17 were surveyed in primary roots, early in the developmental manifestation of heterotic root traits. The application of statistical methods and a suitable experimental design established that 34,233 (i.e., 86%) of all high-confidence maize genes were expressed in at least one genotype. Nearly 70% of all expressed genes were differentially expressed between the two parents and 42%-55% of expressed genes were differentially expressed between one of the parents and one of the hybrids. In both hybrids, ∼10% of expressed genes exhibited nonadditive gene expression. Consistent with the dominance model (i.e., complementation) for heterosis, 1124 genes that were expressed in the hybrids were expressed in only one of the two parents. For 65 genes, it could be shown that this was a consequence of complementation of genomic presence/absence variation. For dozens of other genes, alleles from the inactive inbred were activated in the hybrid, presumably via interactions with regulatory factors from the active inbred. As a consequence of these types of complementation, both hybrids expressed more genes than did either parental inbred. Finally, in hybrids, ∼14% of expressed genes exhibited allele-specific expression (ASE) levels that differed significantly from the parental-inbred expression ratios, providing further evidence for interactions of regulatory factors from one parental genome with target genes from the other parental genome.