Gene expression and expression quantitative trait loci analyses uncover natural variations underlying the improvement of important agronomic traits during modern maize breeding.

Maize (Zea mays L.) is a major staple crop worldwide, and during modern maize breeding, cultivars with increased tolerance to high-density planting and higher yield per plant have contributed significantly to the increased yield per unit land area. Systematically identifying key agronomic traits and their associated genomic changes during modern maize breeding remains a significant challenge because of the complexity of genetic regulation and the interactions of the various agronomic traits, with most of them being controlled by numerous small-effect quantitative trait loci (QTLs). Here, we performed phenotypic and gene expression analyses for a set of 137 elite inbred lines of maize from different breeding eras in China. We found four yield-related traits are significantly improved during modern maize breeding. Through gene-clustering analyses, we identified four groups of expressed genes with distinct trends of expression pattern change across the historical breeding eras. In combination with weighted gene co-expression network analysis, we identified several candidate genes regulating various plant architecture- and yield-related agronomic traits, such as ZmARF16, ZmARF34, ZmTCP40, ZmPIN7, ZmPYL10, ZmJMJ10, ZmARF1, ZmSWEET15b, ZmGLN6 and Zm00001d019150. Further, by combining expression quantitative trait loci (eQTLs) analyses, correlation coefficient analyses and population genetics, we identified a set of candidate genes that might have been under selection and contributed to the genetic improvement of various agronomic traits during modern maize breeding, including a number of known key regulators of plant architecture, flowering time and yield-related traits, such as ZmPIF3.3, ZAG1, ZFL2 and ZmBES1. Lastly, we validated the functional variations in GL15, ZmPHYB2 and ZmPYL10 that influence kernel row number, flowering time, plant height and ear height, respectively. Our results demonstrates the effectiveness of our combined approaches for uncovering key candidate regulatory genes and functional variation underlying the improvement of important agronomic traits during modern maize breeding, and provide a valuable genetic resource for the molecular breeding of maize cultivars with tolerance for high-density planting.

Genetic variation in the aquaporin TONOPLAST INTRINSIC PROTEIN 4;3 modulates maize cold tolerance.

Cold stress is a major abiotic stress that threatens maize (Zea mays L.) production worldwide. Understanding the molecular mechanisms underlying cold tolerance is crucial for breeding resilient maize varieties. Tonoplast intrinsic proteins (TIPs) are a subfamily of aquaporins in plants. Here, we report that TIP family proteins are involved in maize cold tolerance. The expression of most TIP genes was responsive to cold stress. Overexpressing TIP2;1, TIP3;2 or TIP4;3 reduced the cold tolerance of maize seedlings, while loss-of-function mutants of TIP4;3 exhibited enhanced cold tolerance. Candidate gene-based association analysis revealed that a 328-bp transposon insertion in the promoter region of TIP4;3 was strongly associated with maize cold tolerance. This transposon insertion conferred cold tolerance by repressing TIP4;3 expression through increased methylation of its promoter region. Moreover, TIP4;3 was found to suppress stomatal closure and facilitate reactive oxygen species (ROS) accumulation under cold stress, thereby inhibiting the expression of cold-responsive genes, including DEHYDRATION-RESPONSIVE ELEMENT BINDING FACTOR 1 (DREB1) genes and a subset of peroxidase genes, ultimately attenuating maize cold tolerance. This study thus elucidates the mechanism underlying TIP-mediated cold tolerance and identifies a favourable TIP4;3 allele as a potential genetic resource for breeding cold-tolerant maize varieties.

ZmBELL10 interacts with other ZmBELLs and recognizes specific motifs for transcriptional activation to modulate internode patterning in maize.

Plant height is an important agronomic trait that affects crop yield. Elucidating the molecular mechanism underlying plant height regulation is also an important question in developmental biology. Here, we report that a BELL transcription factor, ZmBELL10, positively regulates plant height in maize (Zea mays). Loss of ZmBELL10 function resulted in shorter internodes, fewer nodes, and smaller kernels, while ZmBELL10 overexpression increased plant height and hundred-kernel weight. Transcriptome analysis and chromatin immunoprecipitation followed by sequencing showed that ZmBELL10 recognizes specific sequences in the promoter of its target genes and activates cell division- and cell elongation-related gene expression, thereby influencing node number and internode length in maize. ZmBELL10 interacted with several other ZmBELL proteins via a spatial structure in its POX domain to form protein complexes involving ZmBELL10. All interacting proteins recognized the same DNA sequences, and their interaction with ZmBELL10 increased target gene expression. We identified the key residues in the POX domain of ZmBELL10 responsible for its protein-protein interactions, but these residues did not affect its transactivation activity. Collectively, our findings shed light on the functions of ZmBELL10 protein complexes and provide potential targets for improving plant architecture and yield in maize.

Genomic and Transcriptomic Analyses Reveal Pathways and Genes Associated With Brittle Stalk Phenotype in Maize.

The mechanical strength of the stalk affects the lodging resistance and digestibility of the stalk in maize. The molecular mechanisms regulating the brittleness of stalks in maize remain undefined. In this study, we constructed the maize brittle stalk mutant (bk5) by crossing the W22:Mu line with the Zheng 58 line. The brittle phenotype of the mutant bk5 existed in all of the plant organs after the five-leaf stage. Compared to wild-type (WT) plants, the sclerenchyma cells of bk5 stalks had a looser cell arrangement and thinner cell wall. Determination of cell wall composition showed that obvious differences in cellulose content, lignin content, starch content, and total soluble sugar were found between bk5 and WT stalks. Furthermore, we identified 226 differentially expressed genes (DEGs), with 164 genes significantly upregulated and 62 genes significantly downregulated in RNA-seq analysis. Some pathways related to cellulose and lignin synthesis, such as endocytosis and glycosylphosphatidylinositol (GPI)-anchored biosynthesis, were identified by the Kyoto Encyclopedia of Gene and Genomes (KEGG) and gene ontology (GO) analysis. In bulked-segregant sequence analysis (BSA-seq), we detected 2,931,692 high-quality Single Nucleotide Polymorphisms (SNPs) and identified five overlapped regions (11.2 Mb) containing 17 candidate genes with missense mutations or premature termination codons using the SNP-index methods. Some genes were involved in the cellulose synthesis-related genes such as ENTH/ANTH/VHS superfamily protein gene (endocytosis-related gene) and the lignin synthesis-related genes such as the cytochrome p450 gene. Some of these candidate genes identified from BSA-seq also existed with differential expression in RNA-seq analysis. These findings increase our understanding of the molecular mechanisms regulating the brittle stalk phenotype in maize.

Genome-wide selection and genetic improvement during modern maize breeding.

Since the development of single-hybrid maize breeding programs in the first half of the twentieth century1, maize yields have increased over sevenfold, and much of that increase can be attributed to tolerance of increased planting density2-4. To explore the genomic basis underlying the dramatic yield increase in maize, we conducted a comprehensive analysis of the genomic and phenotypic changes associated with modern maize breeding through chronological sampling of 350 elite inbred lines representing multiple eras of germplasm from both China and the United States. We document several convergent phenotypic changes in both countries. Using genome-wide association and selection scan methods, we identify 160 loci underlying adaptive agronomic phenotypes and more than 1,800 genomic regions representing the targets of selection during modern breeding. This work demonstrates the use of the breeding-era approach for identifying breeding signatures and lays the foundation for future genomics-enabled maize breeding.

Dynamic Changes in Genome-Wide Histone3 Lysine27 Trimethylation and Gene Expression of Soybean Roots in Response to Salt Stress.

Soybean is an important economic crop for human diet, animal feeds and biodiesel due to high protein and oil content. Its productivity is significantly hampered by salt stress, which impairs plant growth and development by affecting gene expression, in part, through epigenetic modification of chromatin status. However, little is known about epigenetic regulation of stress response in soybean roots. Here, we used RNA-seq and ChIP-seq technologies to study the dynamics of genome-wide transcription and histone methylation patterns in soybean roots under salt stress. Eight thousand seven hundred ninety eight soybean genes changed their expression under salt stress treatment. Whole-genome ChIP-seq study of an epigenetic repressive mark, histone H3 lysine 27 trimethylation (H3K27me3), revealed the changes in H3K27me3 deposition during the response to salt stress. Unexpectedly, we found that most of the inactivation of genes under salt stress is strongly correlated with the de novo establishment of H3K27me3 in various parts of the promoter or coding regions where there is no H3K27me3 in control plants. In addition, the soybean histone modifiers were identified which may contribute to de novo histone methylation and gene silencing under salt stress. Thus, dynamic chromatin regulation, switch between active and inactive modes, occur at target loci in order to respond to salt stress in soybean. Our analysis demonstrates histone methylation modifications are correlated with the activation or inactivation of salt-inducible genes in soybean roots.

[Integrating genetic and gene expression data: methods and applications of eQTL mapping].

The availability of high-throughput genotyping technologies and microarray assays has allowed researchers to investigate genetic variations that influence levels of gene expression. Expression Quantitative Trait Locus (eQTL) mapping methods have been used to identify the genetic basis of gene expression. Similar to traditional QTL studies, the main goal of eQTL is to identify the genomic locations to which the expression traits are linked. Although microarrays provide the expression data of thousands of transcripts, standard QTL mapping methods, which are able to handle at most tens of traits, cannot be applied directly. As a result, it is necessary to consider the statistical principles involved in the design and analysis of these experiments. In this paper, we reviewed individual selection, experimental design of microarray, normalization of gene expression data, mapping methods, and explaining of results and proposed potential methodological problems for such analyses. Finally, we discussed the applications of this integrative genomic approach to estimate heritability of transcripts, identify candidate genes, construct gene networks, and understand interactions between genes, genes and environments.

Association Analysis and Identification of ZmHKT1;5 Variation With Salt-Stress Tolerance.

The high-affinity potassium transporter (HKT) genes are essential for plant salt stress tolerance. However, there were limited studies on HKTs in maize (Zea mays), and it is basically unknown whether natural sequence variations in these genes are associated with the phenotypic variability of salt tolerance. Here, the characterization of ZmHKT1;5 was reported. Under salt stress, ZmHKT1;5 expression increased strongly in salt-tolerant inbred lines, which accompanied a better-balanced Na+/K+ ratio and preferable plant growth. The association between sequence variations in ZmHKT1;5 and salt tolerance was evaluated in a diverse population comprising 54 maize varieties from different maize production regions of China. Two SNPs (A134G and A511G) in the coding region of ZmHKT1;5 were significantly associated with different salt tolerance levels in maize varieties. In addition, the favorable allele of ZmHKT1; 5 identified in salt tolerant maize varieties effectively endowed plant salt tolerance. Transgenic tobacco plants of overexpressing the favorable allele displayed enhanced tolerance to salt stress better than overexpressing the wild type ZmHKT1;5. Our research showed that ZmHKT1;5 expression could effectively enhance salt tolerance by maintaining an optimal Na+/K+ balance and increasing the antioxidant activity that keeps reactive oxygen species (ROS) at a low accumulation level. Especially, the two SNPs in ZmHKT1;5 might be related with new amino acid residues to confer salt tolerance in maize. Key Message: Two SNPs of ZmHKT1;5 related with salt tolerance were identified by association analysis. Overexpressing ZmHKT1;5 in tobaccos showed that the SNPs might enhance its ability to regulating Na+/K+ homeostasis.