Weibin Song

Genome-wide patterns of genetic variation among elite maize inbred lines

We have resequenced a group of six elite maize inbred lines, including the parents of the most productive commercial hybrid in China. This effort uncovered more than 1,000,000 SNPs, 30,000 indel polymorphisms and 101 low-sequence-diversity chromosomal intervals in the maize genome. We also identified several hundred complete genes that show presence/absence variation among these resequenced lines. We discuss the potential roles of complementation of presence/absence variations and other deleterious mutations in contributing to heterosis. High-density SNP and indel polymorphism markers reported here are expected to be a valuable resource for future genetic studies and the molecular breeding of this important crop.

Genome-wide genetic changes during modern breeding of maize

The success of modern maize breeding has been demonstrated by remarkable increases in productivity over the last four decades. However, the underlying genetic changes correlated with these gains remain largely unknown. We report here the sequencing of 278 temperate maize inbred lines from different stages of breeding history, including deep resequencing of 4 lines with known pedigree information. The results show that modern breeding has introduced highly dynamic genetic changes into the maize genome. Artificial selection has affected thousands of targets, including genes and non-genic regions, leading to a reduction in nucleotide diversity and an increase in the proportion of rare alleles. Genetic changes during breeding happen rapidly, with extensive variation (SNPs, indels and copy-number variants (CNVs)) occurring, even within identity-by-descent regions. Our genome-wide assessment of genetic changes during modern maize breeding provides new strategies as well as practical targets for future crop breeding and biotechnology.

Extensive, clustered parental imprinting of protein-coding and noncoding RNAs in developing maize endosperm

Although genetic imprinting was discovered in maize 40 years ago, its exact extent in the triploid endosperm remains unknown. Here, we have analyzed global patterns of allelic gene expression in developing maize endosperms from reciprocal crosses between inbreds B73 and Mo17. We have defined an imprinted gene as one in which the relative expression of the maternal and paternal alleles differ at least fivefold in both hybrids of the reciprocal crosses. We found that at least 179 genes (1.6% of protein-coding genes) expressed in the endosperm are imprinted, with 68 of them showing maternal preferential expression and 111 paternal preferential expression. Additionally, 38 long noncoding RNAs were imprinted. The latter are transcribed in either sense or antisense orientation from intronic regions of normal protein-coding genes or from intergenic regions. Imprinted genes show a clear pattern of clustering around the genome, with a number of imprinted genes being adjacent to each other. Analysis of allele-specific methylation patterns of imprinted loci in the hybrid endosperm identified 21 differentially methylated regions (DMRs) of several hundred base pairs in length, corresponding to both imprinted genes and noncoding transcripts. All DMRs identified are uniformly hypomethylated in maternal alleles and hypermethylated in paternal alleles, regardless of the imprinting direction of their corresponding loci. Our study indicates highly extensive and complex regulation of genetic imprinting in maize endosperm, a mechanism that can potentially function in the balancing of the gene dosage of this triploid tissue.

An ultra-high density bin-map for rapid QTL mapping for tassel and ear architecture in a large F2 maize population

BackgroundUnderstanding genetic control of tassel and ear architecture in maize (Zea mays L. ssp. mays) is important due to their relationship with grain yield. High resolution QTL mapping is critical for understanding the underlying molecular basis of phenotypic variation. Advanced populations, such as recombinant inbred lines, have been broadly adopted for QTL mapping; however, construction of large advanced generation crop populations is time-consuming and costly. The rapidly declining cost of genotyping due to recent advances in next-generation sequencing technologies has generated new possibilities for QTL mapping using large early generation populations.ResultsA set of 708 F2 progeny derived from inbreds Chang7-2 and 787 were generated and genotyped by whole genome low-coverage genotyping-by-sequencing method (average 0.04×). A genetic map containing 6,533 bin-markers was constructed based on the parental SNPs and a sliding-window method, spanning a total genetic distance of 1,396 cM. The high quality and accuracy of this map was validated by the identification of two well-studied genes, r1, a qualitative trait locus for color of silk (chromosome 10) and ba1 for tassel branch number (chromosome 3). Three traits of tassel and ear architecture were evaluated in this population, a total of 10 QTL were detected using a permutation-based-significance threshold, seven of which overlapped with reported QTL. Three genes (GRMZM2G316366, GRMZM2G492156 and GRMZM5G805008) encoding MADS-box domain proteins and a BTB/POZ domain protein were located in the small intervals of qTBN5 and qTBN7 (~800 Kb and 1.6 Mb in length, respectively) and may be involved in patterning of tassel architecture. The small physical intervals of most QTL indicate high-resolution mapping is obtainable with this method.ConclusionsWe constructed an ultra-high-dentisy linkage map for the large early generation population in maize. Our study provides an efficient approach for fast detection of quantitative loci responsible for complex trait variation with high accuracy, thus helping to dissect the underlying molecular basis of phenotypic variation and accelerate improvement of crop breeding in a cost-effective fashion.

Efficiency and Inheritance of Targeted Mutagenesis in Maize Using CRISPR-Cas9.

CRISPR/Cas (clustered regularly interspaced short palindromic repeats/CRISPR-associated proteins) is an adaptive immune system in bacteria and archaea to defend against invasion from foreign DNA fragments. Recently, it has been developed as a powerful targeted genome editing tool for a wide variety of species. However, its application in maize has only been tested with transiently expressed somatic cells or with a limited number of stable transgenic T0 plants. The exact efficiency and specificity of the CRISPR/Cas system in the highly complex maize genome has not been documented yet. Here we report an extensive study of the well-studied type II CRISPR-Cas9 system for targeted genome editing in maize, with the codon-optimized Cas9 protein and the short non-coding guide RNA generated through a functional maize U6 snRNA promoter. Targeted gene mutagenesis was detected for 90 loci by maize protoplast assay, with an average cleavage efficiency of 10.67%. Stable knockout transformants for maize phytoene synthase gene (PSY1) were obtained. Mutations occurred in germ cells can be stably inherited to the next generation. Moreover, no off-target effect was detected at the computationally predicted putative off-target loci. No significant difference between the transcriptomes of the Cas9 expressed and non-expressed lines was detected. Our results confirmed that the CRISPR-Cas9 could be successfully applied as a robust targeted genome editing system in maize.

Maize LAZY1 Mediates Shoot Gravitropism and Inflorescence Development through Regulating Auxin Transport, Auxin Signaling, and Light Response

Map-based cloning of a maize prostrate-growth mutant links gravitropism with inflorescence development through auxin transport, auxin signaling, and response to light. Auxin is a plant hormone that plays key roles in both shoot gravitropism and inflorescence development. However, these two processes appear to be parallel and to be regulated by distinct players. Here, we report that the maize (Zea mays) prostrate stem1 mutant, which is allelic to the classic mutant lazy plant1 (la1), displays prostrate growth with reduced shoot gravitropism and defective inflorescence development. Map-based cloning identified maize ZmLA1 as the functional ortholog of LAZY1 in rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana). It has a unique role in inflorescence development and displays enriched expression in reproductive organs such as tassels and ears. Transcription of ZmLA1 responds to auxin and is repressed by light. Furthermore, ZmLA1 physically interacts with a putative auxin transport regulator in the plasma membrane and a putative auxin signaling protein in the nucleus. RNA-SEQ data showed that dozens of auxin transport, auxin response, and light signaling genes were differentially expressed in la1 mutant stems. Therefore, ZmLA1 might mediate the cross talk between shoot gravitropism and inflorescence development by regulating auxin transport, auxin signaling, and probably light response in maize.

Ribosome profiling reveals dynamic translational landscape in maize seedlings under drought stress

Plants can respond to environmental changes with various mechanisms occurred at transcriptional and translational levels. Thus far, there have been relatively extensive understandings of stress responses of plants on transcriptional level, while little information is known about that on translational level. To uncover the landscape of translation in plants in response to drought stress, we performed the recently developed ribosome profiling assay with maize seedlings growing under normal and drought conditions. Comparative analysis of the ribosome profiling data and the RNA-seq data showed that the fold changes of gene expression at transcriptional level were moderately correlated with that of translational level globally (R(2) = 0.69). However, less than half of the responsive genes were shared by transcription and translation under drought condition, suggesting that drought stress can introduce transcriptional and translational responses independently. We found that the translational efficiencies of 931 genes were changed significantly in response to drought stress. Further analysis revealed that the translational efficiencies of genes were highly influenced by their sequence features including GC content, length of coding sequences and normalized minimal free energy. In addition, we detected potential translation of 3063 upstream open reading frames (uORFs) on 2558 genes and these uORFs may affect the translational efficiency of downstream main open reading frames (ORFs). Our study indicates that plant can respond to drought stress with highly dynamic translational mechanism, that acting synergistically with that of transcription.

Isolation and Analysis of Rice Rf1-Orthologus PPR Genes Co-segregating with Rf3 in Maize

Using an in silico cloning approach, five putative maize pentatricopeptide repeat (PPR)-containing protein genes (PPR-814a, PPR-814b, PPR-814c, PPR-816, PPR-817) with complete open reading frames were identified in the inbred line S-Mo17Rf3Rf3. The amino acid sequence indicated that these genes encoded mitochondrially targeted proteins containing repeats of a 35-aa PPR motif. The genes were mapped into the interval umc1525–bnlg1520 on chromosome 2. In a non-restoring genotype, we identified three homologous genes that contained deletions or nucleotide substitutions in the coding region. Sequence analysis revealed that one of the three genes (PPR-814a, PPR-814b, PPR-814c) could be considered a candidate restorer gene for S male sterility cytoplasm, and linkage analysis demonstrated that the genes co-segregated with the fertility restorer gene Rf3.

Identification of genetic factors affecting plant density response through QTL mapping of yield component traits in maize (Zea mays L.)

It is generally believed that grain yield per unit area of modern maize hybrids is related to their adaptability to high plant population density. In this study, the effects of two different plant densities (52,500 and 90,000 plants/hm2) on 12 traits associated with yield were evaluated using a set of 231 F2:3 families derived from two elite inbred lines, Zheng58 and Chang7-2. Evaluation of the phenotypes expressed under the two plant density conditions showed that high plant density condition could decrease the value of 10 measured yield component traits, while the final grain yield per hectare and the rate of kernel production were increased. Twenty-seven quantitative trait loci (QTLs) for 10 traits were detected in both high and low plant density conditions; among them, some QTLs were shown to locate in five clusters. Thirty QTLs were only detected under high plant density. These results suggest that some of the yield component traits perhaps were controlled by a common set of genes, and that kernel number per row, ear length, row number per ear, cob diameter, cob weight, and ear diameter may be influenced by additional genetic mechanisms when grown under high plant density. The QTLs identified in this study provide useful information for marker-assisted selection of varieties targeting increased plant density.

Map-Based Cloning of zb7 Encoding an IPP and DMAPP Synthase in the MEP Pathway of Maize

IspH is a key enzyme in the last step of the methyl-D-erythritol-4-phosphate (MEP) pathway. Loss of function of IspH can often result in complete yellow or albino phenotype in many plants. Here, we report the characterization of a recessive mutant of maize, zebra7 (zb7), showing transverse green/yellow striped leaves in young plants. The yellow bands of the mutant have decreased levels of chlorophylls and carotenoids with delayed chloroplast development. Low temperature suppressed mutant phenotype, while alternate light/dark cycle or high temperature enlarged the yellow section. Map-based cloning demonstrated that zb7 encodes the IspH protein with a mis-sense mutation in a conserved region. Transgenic silencing of Zb7 in maize resulted in complete albino plantlets that are aborted in a few weeks, confirming that Zb7 is important in the early stages of maize chloroplast development. Zb7 is constitutively expressed and its expression subject to a 16-h light/8-h dark cycle regulation. Our results suggest that the less effective or unstable IspH in zb7 mutant, together with its diurnal expression, are mechanistically accounted for the zebra phenotype. The increased IspH mRNA in the leaves of zb7 at the late development stage may explain the restoration of mutant phenotype in mature stages.