Conserved noncoding genomic sequences associated with a flowering-time quantitative trait locus in maize

Exploring the multifaceted dynamics of flowering time regulation in field crops: Insight and intervention approaches

The flowering time (FTi) plays a critical role in the reproductive success and yield of various crop species by directly impacting both the quality and quantity of grain yield. Achieving optimal FTi is crucial for maximizing reproductive success and ensuring overall agricultural productivity. While genetic factors undoubtedly influence FTi, photoperiodism and vernalization are recognized as key contributors to the complex physiological processes governing flowering in plants. Identifying candidate genes and pathways associated with FTi is essential for developing genomic interventions and plant breeding to enhance adaptability to diverse environmental conditions. This review highlights the intricate nature of the regulatory mechanisms of flowering and emphasizes the vital importance of precisely regulating FTi to ensure plant adaptability and reproductive success. Special attention is given to essential genes, pathways, and genomic interventions geared toward promoting early flowering, particularly under challenging environmental conditions such as drought, heat, and cold stress as well as other abiotic stresses that occur during the critical flowering stage of major field crops. Moreover, this review explores the significant progress achieved in omics technologies, offering valuable insights and tools for deciphering and regulating FTi. In summary, this review aims to provide a comprehensive understanding of the mechanisms governing FTi, with a particular focus on their crucial role in bolstering yields under adverse environmental conditions to safeguard food security.

A maize landrace introgression library reveals a negative effect of root-to-shoot ratio on water-use efficiency

Novel sources of genetic variability for water-use efficiency (WUE) are needed in order to breed varieties more suitable to sustainable cropping systems. Here, a maize (Zea mays L.) introgression library of the landrace Gaspé Flint into the reference line B73 was characterized in high-throughput phenotyping platforms, both in well-watered and moderate water-deficit conditions, for water use, WUE, and root and shoot growth. Traits heritability ranged from 0.77 to 0.93. The introgression of Gaspé Flint chromosome segments into the B73 genome significantly altered several traits. Some introgression lines exhibited a faster shoot biomass accumulation than B73, resulting in higher WUE at the expense of root growth. Quantitative trait loci (QTL) mapping identified seven major QTL clusters affecting shoot growth and WUE, two of which overlapped, with opposite effects, with QTLs for root biomass known to include root developmental genes. These results support the non-intuitive hypothesis that reduced root-to-shoot ratio positively affects maize WUE.

A copy number variation in the ZmMADS1 promoter enhances maize adaptation to high altitudes

A 178-bp copy number variation in the ZmMADS1 promoter represses ZmMADS1 expression and results in delayed flowering, which has been a target of selection during maize spread into higher altitudes. ZmMADS1 promotes maize flowering via the autonomous pathway by directly activating the expression of ZmMADS69 and ZCN8, while repressing ZmRap2.7.

Enhancers in Plant Development, Adaptation and Evolution

Understanding plant responses to developmental and environmental cues is crucial for studying morphological divergence and local adaptation. Gene expression changes, governed by cis-regulatory modules (CRMs) including enhancers, are a major source of plant phenotypic variation. However, while genome-wide approaches have revealed thousands of putative enhancers in mammals, far fewer have been identified and functionally characterized in plants. This review provides an overview of how enhancers function to control gene regulation, methods to predict DNA sequences that may have enhancer activity, methods utilized to functionally validate enhancers and the current knowledge of enhancers in plants, including how they impact plant development, response to environment and evolutionary adaptation.

Genetic basis of Fusarium ear rot resistance and productivity traits in a heterozygous multi-parent recombinant inbred intercross (RIX) maize population

Maize (Zea mays L.) is one of the most productive crops worldwide. As a heterotic crop predominantly grown as F1 hybrid, maize exhibits challenges for genetic studies of complex traits, since homozygous genotypes, which are largely used in these studies, may not accurately reflect what happens in cultivated conditions. To map Fusarium Ear Rot (FER) resistance to Fusarium verticillioides and traits with potential impact on yield, including phenology, we constructed a recombinant intercross (RIX) population. This was achived by crossing pairs of recombinant inbred lines (RILs) derived from a multi-parent maize population. We characterized the RIX population over two growing seasons, employing artificial F. verticillioides inoculation. The heterozygous background of the material enabled the identification of QTL and candidate genes through in silico reconstruction of RIX genotype probabilities. A total of 37 loci were identified using single-year BLUPs while 29 with joint-year BLUPs. These, included several known QTL associated with days to tasseling, kernel row number and a QTL on the chromosome 9 associated with FER resistance. In this region, we could identify candidates based on their predicted functions and potential roles in plant-pathogen interactions and/or resistance mechanisms. These QTL represent potential breeding targets to FER resistance and yield components in commercial maize varieties.

Diverse roles of ethylene in maize growth and development, and its importance in shaping plant architecture

The gaseous plant hormone ethylene is a key developmental and growth regulator, and a pivotal endogenous response signal to abiotic and biotic interactions, including stress. Much of what is known about ethylene biosynthesis, perception, and signaling comes from decades of research primarily in Arabidopsis thaliana and other eudicot model systems. In contrast, detailed knowledge on the ethylene pathway and response to the hormone is markedly limited in maize (Zea mays L.), a global cereal crop that is a major source of calories for humans and livestock, as well as a key industrial biofeedstock. Recent reports of forward screens and targeted reverse genetics have provided important insight into conserved and unique differences of the ethylene pathway and downstream responses. Natural and edited allelic variation in the promoter regions and coding sequences of ethylene biosynthesis and signaling genes alters maize shoot and root architectures, and plays a crucial role in biomass and grain yields. This review discusses recent advances in ethylene research in maize, with an emphasis on the role of ethylene in regulating growth and development of the shoot and root systems, and ultimately how this crucial hormone impacts plant architecture and grain yield.

The genetic architecture of cell type-specific cis regulation in maize

Gene expression and complex phenotypes are determined by the activity of cis-regulatory elements. However, an understanding of how extant genetic variants affect cis regulation remains limited. Here, we investigated the consequences of cis-regulatory diversity using single-cell genomics of more than 0.7 million nuclei across 172 Zea mays (maize) inbreds. Our analyses pinpointed cis-regulatory elements distinct to domesticated maize and revealed how historical transposon activity has shaped the cis-regulatory landscape. Leveraging population genetics principles, we fine-mapped about 22,000 chromatin accessibility-associated genetic variants with widespread cell type-specific effects. Variants in TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR-binding sites were the most prevalent determinants of chromatin accessibility. Finally, integrating chromatin accessibility-associated variants, organismal trait variation, and population differentiation revealed how local adaptation has rewired regulatory networks in unique cellular contexts to alter maize flowering.

Understanding the Regulation Activities of Transposons in Driving the Variation and Evolution of Polyploid Plant Genome

Transposon is the main component of the eukaryotic genome, and more and more plant genome data show that transposons are diverse in regulating genome structure, variation, function and evolution, with different transposition mechanisms in the genome. Hybridization and polyploidy play an important role in promoting plant speciation and evolution, and recent studies have shown that polyploidy is usually accompanied by the expansion of transposons, which affect the genome size and structure of polyploid plants. Transposons can insert into genes and intergenic regions, resulting in great differences in the overall genome structure of closely related plant species, and it can also capture gene segments in the genome to increase the copy number of genes. In addition, transposons influence the epigenetic modification state of the genome and regulate the expression of the gene, while plant phenotype, biological and abiotic stress response are also regulated by transposons. Overall, transposons play an important role in the plant genome, especially polyploid plant genome, adaptation and evolution.

Improving yield-related traits by editing the promoter and distal regulatory region of heading date genes Ghd7 and PRR37 in elite rice variety Mei Xiang Zhan 2

KEY MESSAGE

We revealed that editing the promoter and distal regulatory region of the pleiotropic genes Ghd7 and PRR37 reduces their ability to delay heading date while improving their capacity to boost crop yield, offering valuable resources for rice breeding. Heading date is a crucial agronomic characteristic in rice that governs the adaptability to different latitudes and the yield of various varieties. Optimizing the heading date of superior cultivars in breeding practice can significantly broaden their potential planting areas. Ghd7 and PRR37 are pivotal genes that control heading date and enhance agronomic traits. In the elite indica rice variety Mei Xiang Zhan 2 (MXZ2), we used CRISPR/Cas9 technology to effectively generate homozygous mutant lines with a gradient change in heading date by multi-target editing the promoter and distal regulatory region of Ghd7 and PRR37. Various degrees of down-regulation of Ghd7 or PRR37 expression, impaired gene functions, and advancement of the heading date were observed in the mutant lines. Certain mutant lines exhibited an early heading date and increased yield while preserving the exceptional quality of MXZ2. Our study revealed that editing the promoter and distal regulatory region of the pleiotropic genes Ghd7 and PRR37 reduces their ability to delay heading date while improving their capacity to boost crop yield, offering valuable resources for rice breeding.

A LTR retrotransposon insertion leads to leafy phenotype in maize by elevating ZmOM66 expression

Leafy (Lfy1) is a classical dominant mutant showing more leaf number above primary ear and later flowering time in maize, but the causal gene together with its underlying genetic mechanism are unknown. Here, we report the cloning of Lfy1 mutant, and find that a retrotransposon insertion leads to leafy phenotype by elevating expression of its neighboring gene ZmOM66. ZmOM66 encodes an AAA+ ATPase that locate in mitochondria and interacts with itself. ZmOM66 overexpression affects the starch degradation, as well as contents of glucose, pyruvic acid, trehalose-6-phosphate, and TCA cycle related amino acids, and influences expression patterns of circadian clock genes. Moreover, expressions of floral related genes, including photoperiod regulated gene ZmPHYB1, integrator genes ZCN7, ZNC8 and ZCN12, and floral meristem identity genes ZMM4, ZMM15, and MASD67, are also significantly decreased by ZmOM66 overexpression. These results deepen our understanding of the regulatory mechanism of floral transition and leaf number in plant.

Natural variation in ZmNRT2.5 modulates husk leaf width and promotes seed protein content in maize

The husk leaf of maize (Zea mays) encases the ear as a modified leaf and plays pivotal roles in protecting the ear from pathogen infection, translocating nutrition for grains and warranting grain yield. However, the natural genetic basis for variation in husk leaf width remains largely unexplored. Here, we performed a genome-wide association study for maize husk leaf width and identified a 3-bp InDel (insertion/deletion) in the coding region of the nitrate transporter gene ZmNRT2.5. This polymorphism altered the interaction strength of ZmNRT2.5 with another transporter, ZmNPF5, thereby contributing to variation in husk leaf width. We also isolated loss-of-function mutants in ZmNRT2.5, which exhibited a substantial decrease in husk leaf width relative to their controls. We demonstrate that ZmNRT2.5 facilitates the transport of nitrate from husk leaves to maize kernels in plants grown under low-nitrogen conditions, contributing to the accumulation of proteins in maize seeds. Together, our findings uncovered a key gene controlling maize husk leaf width and nitrate transport from husk leaves to kernels. Identification of the ZmNRT2.5 loci offers direct targets for improving the protein content of maize seeds via molecular-assisted maize breeding.

Precise engineering of gene expression by editing plasticity

BACKGROUND

Identifying transcriptional cis-regulatory elements (CREs) and understanding their role in gene expression are essential for the precise manipulation of gene expression and associated phenotypes. This knowledge is fundamental for advancing genetic engineering and improving crop traits.

RESULTS

We here demonstrate that CREs can be accurately predicted and utilized to precisely regulate gene expression beyond the range of natural variation. We firstly build two sequence-to-expression deep learning models to respectively identify distal and proximal CREs by combining them with interpretability methods in multiple crops. A large number of distal CREs are verified for enhancer activity in vitro using UMI-STARR-seq on 12,000 synthesized sequences. These comprehensively characterized CREs and their precisely predicted effects further contribute to the design of in silico editing schemes for precise engineering of gene expression. We introduce a novel concept of "editingplasticity" to evaluate the potential of promoter editing to alter expression of each gene. As a proof of concept, both exhaustive prediction and random knockout mutants are analyzed within the promoter region of ZmVTE4, a key gene affecting α-tocopherol content in maize. A high degree of agreement between predicted and observed expression is observed, extending the range of natural variation and thereby allowing the creation of an optimal phenotype.

CONCLUSIONS

Our study provides a robust computational framework that advances knowledge-guided gene editing for precise regulation of gene expression and crop improvement. By reliably predicting and validating CREs, we offer a tool for targeted genetic modifications, enhancing desirable traits in crops.

The chromatin accessibility landscape during early maize seed development

Cis-regulatory elements (CREs) are enriched in accessible chromatin regions (ACRs) of eukaryotes. Despite extensive research on genome-wide ACRs in various plant tissues, the global impact of these changes on developmental processes in maize seeds remains poorly understood. In this study, we employed the assay for transposase-accessible chromatin sequencing (ATAC-seq) to reveal the chromatin accessibility profile throughout the genome during the early stages of maize seed development. We identified a total of 37 952 to 59 887 high-quality ACRs in maize seeds at 0 to 8 days after pollination (DAP). Furthermore, we examined the correlation between the identified ACRs and gene expression. We observed a positive correlation between the open degree of promoter-ACRs and the expression of most genes. Moreover, we identified binding footprints of numerous transcription factors (TFs) within chromatin accessibility regions and revealed key TF families involved in different stages. Through the footprints of accessible chromatin regions, we predicted transcription factor regulatory networks during early maize embryo development. Additionally, we discovered that DNA sequence diversity was notably reduced at ACRs, yet trait-associated SNPs were more likely to be located within ACRs. We edited the ACR containing the trait-associated SNP of NKD1. Both NKD1pro-1 and NKD1pro-2 showed phenotypes corresponding to the trait-associated SNP. Our results suggest that alterations in chromatin accessibility play a crucial role in maize seed development and highlight the potential contribution of open chromatin regions to advancements in maize breeding.

Evolving best practices for transcriptome-wide association studies accelerate discovery of gene-phenotype links

Transcriptome-wide association studies (TWAS) complement genome-wide association studies (GWAS) by using gene expression data to link specific genes to phenotypes. This review examines 37 TWAS studies across eight plant species, evaluating the impact of methodological choices on outcomes using maize and soybean datasets. Large sample sizes and synchronized sample collection for gene expression measurement appear to significantly increase power for discovering gene-phenotype linkages, while matching tissue, stage, and environment may matter much less than previously believed, making it feasible to reuse large and well-collected expression datasets across multiple studies. The development of statistical approaches and computational tools specifically optimized for plant TWAS data will ultimately be needed, but further potential remains to adapt advances developed in GWAS to TWAS contexts.

Characterization and functional analysis of conserved non-coding sequences among poaceae: insights into gene regulation and phenotypic variation in maize

BACKGROUND

Conserved non-coding sequences (CNS) are islands of non-coding sequences conserved across species and play an important role in regulating the spatiotemporal expression of genes. Identification of CNS provides valuable information about potentially functional genomic elements, regulatory regions, and helps to gain insights into the genetic basis of crop agronomic traits.

RESULTS

Here, we comprehensively analyze CNS in maize, by comparing the genomes of maize inbred line B73 (Zea mays ssp. mays), its close wild relative Zea mays spp. mexicana, and other grasses in Poaceae, including sorghum (Sorghum bicolor), foxtail millet (Setaria italica) and two adlay (Coix lacryma) cultivars. There were 289,931 CNS found in two syntenic gene pairs, while 51,701 CNS were conserved within at least three species. To explore the regulatory characteristics of the CNS identified, the flanking regions of CNS were compared with the peaks called using both transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and chromatin immunoprecipitation with high-throughput sequencing (ChIP-Seq) data of histone modifications. It was found that CNS in maize were enriched in open chromatin regions compared with randomly selected non-coding regions of similar length. A significant enrichment of transcription factor binding sites was found within CNS sequences, including different transcription factors involved in abiotic stress response, such as OBP (OBF-BINDING PROTEIN) family and Adof1 (Encodes dof zinc finger protein). To investigate the epigenetic modification patterns in CNS, ChIP-Seq data for histone modifications H3K9ac, H3K4me3, H3K36me3, H3K9me3, and H3K27ac were further analyzed to depict the changes along CNS. Our findings revealed significantly elevated levels of transcription-promoting histone modifications in the CNS regions compared to randomly selected non-coding sequences with an equal number and similar length. Notably, CNS were also identified on both Vgt1 (Vegetative to generative transition 1) and ZmCCT10. In addition, CNS with potential functions were identified based on SNPs within CNS significantly associated with various agronomic traits in maize, which holds potential utility in molecular breeding for maize.

CONCLUSIONS

In summary, we identified and characterized CNS in maize through genomic comparative analysis, which provides valuable insights into their potential regulatory effects on gene expression and phenotypic variation.

A Zea genus-specific micropeptide controls kernel dehydration in maize

Kernel dehydration rate (KDR) is a crucial production trait that affects mechanized harvesting and kernel quality in maize; however, the underlying mechanisms remain unclear. Here, we identified a quantitative trait locus (QTL), qKDR1, as a non-coding sequence that regulates the expression of qKDR1 REGULATED PEPTIDE GENE (RPG). RPG encodes a 31 amino acid micropeptide, microRPG1, which controls KDR by precisely modulating the expression of two genes, ZmETHYLENE-INSENSITIVE3-like 1 and 3, in the ethylene signaling pathway in the kernels after filling. microRPG1 is a Zea genus-specific micropeptide and originated de novo from a non-coding sequence. Knockouts of microRPG1 result in faster KDR in maize. By contrast, overexpression or exogenous application of the micropeptide shows the opposite effect both in maize and Arabidopsis. Our findings reveal the molecular mechanism of microRPG1 in kernel dehydration and provide an important tool for future crop breeding.

LOGOWheat: deep learning-based prediction of regulatory effects for noncoding variants in wheats

Identifying the regulatory effects of noncoding variants presents a significant challenge. Recently, the accumulation of epigenomic profiling data in wheat has provided an opportunity to model the functional impacts of these variants. In this study, we introduce Language of Genome for Wheat (LOGOWheat), a deep learning-based tool designed to predict the regulatory effects of noncoding variants in wheat. LOGOWheat initially employs a self-attention-based, contextualized pretrained language model to acquire bidirectional representations of the unlabeled wheat reference genome. Epigenomic profiling data are also collected and utilized to fine-tune the model, enabling it to discern the regulatory code inherent in genomic sequences. The test results suggest that LOGOWheat is highly effective in predicting multiple chromatin features, achieving an average area under the receiver operating characteristic (AUROC) of 0.8531 and an average area under the precision-recall curve (AUPRC) of 0.7633. Two case studies illustrate and demonstrate the main functions provided by LOGOWheat: assigning scores and prioritizing causal variants within a given variant set and constructing a saturated mutagenesis map in silico to discover high-impact sites or functional motifs in a given sequence. Finally, we propose the concept of extracting potential functional variations from the wheat population by integrating evolutionary conservation information. LOGOWheat is available at http://logowheat.cn/.

Transgenerational epigenetic inheritance during plant evolution and breeding

Plants can program and reprogram their genomes to create genetic variation and epigenetic modifications, leading to phenotypic plasticity. Although consequences of genetic changes are comprehensible, the basis for transgenerational inheritance of epigenetic variation is elusive. This review addresses contributions of external (environmental) and internal (genomic) factors to the establishment and maintenance of epigenetic memory during plant evolution, crop domestication, and modern breeding. Dynamic and pervasive changes in DNA methylation and chromatin modifications provide a diverse repertoire of epigenetic variation potentially for transgenerational inheritance. Elucidating and harnessing epigenetic inheritance will help us develop innovative breeding strategies and biotechnological tools to improve crop yield and resilience in the face of environmental challenges. Beyond plants, epigenetic principles are shared across sexually reproducing organisms including humans with relevance to medicine and public health.

Genetic variability of aquaporin expression in maize: From eQTLs to a MITE insertion regulating PIP2;5 expression

Plant aquaporins are involved in numerous physiological processes, such as cellular homeostasis, tissue hydraulics, transpiration, and nutrient supply, and are key players of the response to environmental cues. While varying expression patterns of aquaporin genes have been described across organs, developmental stages, and stress conditions, the underlying regulation mechanisms remain elusive. Hence, this work aimed to shed light on the expression variability of 4 plasma membrane intrinsic protein (PIP) genes in maize (Zea mays) leaves, and its genetic causes, through expression quantitative trait locus (eQTL) mapping across a 252-hybrid diversity panel. Significant genetic variability in PIP transcript abundance was observed to different extents depending on the isoforms. The genome-wide association study mapped numerous eQTLs, both local and distant, thus emphasizing the existing natural diversity of PIP gene expression across the studied panel and the potential to reveal regulatory actors and mechanisms. One eQTL associated with PIP2;5 expression variation was characterized. Genomic sequence comparison and in vivo reporter assay attributed, at least partly, the local eQTL to a transposon-containing polymorphism in the PIP2;5 promoter. This work paves the way to the molecular understanding of PIP gene regulation and its possible integration into larger networks regulating physiological and stress adaptation processes.

GenomicLinks: deep learning predictions of 3D chromatin interactions in the maize genome

Gene regulation in eukaryotes is partly shaped by the 3D organization of chromatin within the cell nucleus. Distal interactions between cis-regulatory elements and their target genes are widespread, and many causal loci underlying heritable agricultural traits have been mapped to distal non-coding elements. The biology underlying chromatin loop formation in plants is poorly understood. Dissecting the sequence features that mediate distal interactions is an important step toward identifying putative molecular mechanisms. Here, we trained GenomicLinks, a deep learning model, to identify DNA sequence features predictive of 3D chromatin interactions in maize. We found that the presence of binding motifs of specific transcription factor classes, especially bHLH, is predictive of chromatin interaction specificities. Using an in silico mutagenesis approach we show the removal of these motifs from loop anchors leads to reduced interaction probabilities. We were able to validate these predictions with single-cell co-accessibility data from different maize genotypes that harbor natural substitutions in these TF binding motifs. GenomicLinks is currently implemented as an open-source web tool, which should facilitate its wider use in the plant research community.

A genome-wide association study uncovers a ZmRap2.7-ZCN9/ZCN10 module to regulate ABA signalling and seed vigour in maize

Seed vigour, including rapid, uniform germination and robust seedling establishment under various field conditions, is becoming an increasingly essential agronomic trait for achieving high yield in crops. However, little is known about this important seed quality trait. In this study, we performed a genome-wide association study to identify a key transcription factor ZmRap2.7, which regulates seed vigour through transcriptionally repressing expressions of three ABA signalling genes ZmPYL3, ZmPP2C and ZmABI5 and two phosphatidylethanolamine-binding genes ZCN9 and ZCN10. In addition, ZCN9 and ZCN10 proteins could interact with ZmPYL3, ZmPP2C and ZmABI5 proteins, and loss-of-function of ZmRap2.7 and overexpression of ZCN9 and ZCN10 reduced ABA sensitivity and seed vigour, suggesting a complex regulatory network for regulation of ABA signalling mediated seed vigour. Finally, we showed that four SNPs in ZmRap2.7 coding region influenced its transcriptionally binding activity to the downstream gene promoters. Together with previously identified functional variants within and surrounding ZmRap2.7, we concluded that the distinct allelic variations of ZmRap2.7 were obtained independently during maize domestication and improvement, and responded separately for the diversities of seed vigour, flowering time and brace root development. These results provide novel genes, a new regulatory network and an evolutional mechanism for understanding the molecular mechanism of seed vigour.

Enhancing genomic association studies in slash pine through close-range UAV-based morphological phenotyping

In forestry genetics and industry, tree morphological traits such as height, crown size, and shape are critical for understanding growth dynamics and productivity. Traditional methods for measuring these traits are limited in efficiency, scalability, and accuracy, posing challenges for large-scale forest assessments. This study focuses on integrating unmanned aerial vehicle (UAV) technology with GWAS to improve genomic association studies in slash pine (Pinus elliottii). Seven key morphological traits have been identified (canopy area (CA), crown base height (CBH), crown length (CL), canopy volume (CV), crown width (CW), crown width height (CWH), and tree height (H)) through advanced UAV-based phenotyping. These associations account for a remarkable range of heritability in slash pine, with traits such as CBH, CL, CV, and H showing relatively high heritability across both Single nucleotide polymorphisms (SNP) and pedigree methods, indicating strong genetic influence, while traits such as CWH show lower heritability, suggesting greater environmental influence or non-additive genetic variance. The GWAS identified 28 associations, including 22 different SNPs localized to 16 candidate genes, that were significantly associated with the morphological traits of Slash Pine. Notably, two of these candidate genes, annotated as putative DEAD-like helicase and ethylene-responsive element binding factor (ERF), were present at different mutation sites and were significantly associated with CW and CA traits, respectively. These results demonstrate that the UAV imaging enables a comprehensive analysis of the Morphological growth response of slash pine and can facilitate the discovery of informative alleles to elucidate the genetic structure underlying complex phenotypic variation in conifers.

Lateral root enriched Massilia associated with plant flowering in maize

BACKGROUND

Beneficial associations between plants and soil microorganisms are critical for crop fitness and resilience. However, it remains obscure how microorganisms are assembled across different root compartments and to what extent such recruited microbiomes determine crop performance. Here, we surveyed the root transcriptome and the root and rhizosphere microbiome via RNA sequencing and full-length (V1-V9) 16S rRNA gene sequencing from genetically distinct monogenic root mutants of maize (Zea mays L.) under different nutrient-limiting conditions.

RESULTS

Overall transcriptome and microbiome display a clear assembly pattern across the compartments, i.e., from the soil through the rhizosphere to the root tissues. Co-variation analysis identified that genotype dominated the effect on the microbial community and gene expression over the nutrient stress conditions. Integrated transcriptomic and microbial analyses demonstrated that mutations affecting lateral root development had the largest effect on host gene expression and microbiome assembly, as compared to mutations affecting other root types. Cooccurrence and trans-kingdom network association analysis demonstrated that the keystone bacterial taxon Massilia (Oxalobacteraceae) is associated with root functional genes involved in flowering time and overall plant biomass. We further observed that the developmental stage drives the differentiation of the rhizosphere microbial assembly, especially the associations of the keystone bacteria Massilia with functional genes in reproduction. Taking advantage of microbial inoculation experiments using a maize early flowering mutant, we confirmed that Massilia-driven maize growth promotion indeed depends on flowering time.

CONCLUSION

We conclude that specific microbiota supporting lateral root formation could enhance crop performance by mediating functional gene expression underlying plant flowering time in maize. Video Abstract.

Transcription factor binding site divergence across maize inbred lines drives transcriptional and phenotypic variation

Regulatory elements are important constituents of plant genomes that have shaped ancient and modern crops. Their identification, function, and diversity in crop genomes however are poorly characterized, thus limiting our ability to harness their power for further agricultural advances using induced or natural variation. Here, we use DNA affinity purification-sequencing (DAP-seq) to map transcription factor (TF) binding events for 200 maize TFs belonging to 30 distinct families and heterodimer pairs in two distinct inbred lines historically used for maize hybrid plant production, providing empirical binding site annotation for 5.3% of the maize genome. TF binding site comparison in B73 and Mo17 inbreds reveals widespread differences, driven largely by structural variation, that correlate with gene expression changes. TF binding site presence-absence variation helps clarify complex QTL such as vgt1, an important determinant of maize flowering time, and DICE, a distal enhancer involved in herbivore resistance. Modification of TF binding regions via CRISPR-Cas9 mediated editing alters target gene expression and phenotype. Our functional catalog of maize TF binding events enables collective and comparative TF binding analysis, and highlights its value for agricultural improvement.

FKF1b controls reproductive transition associated with adaptation to geographical distribution in maize

Maize (Zea mays subspecies mays) is an important commercial crop across the world, and its flowering time is closely related to grain yield, plant cycle and latitude adaptation. FKF1 is an essential clock-regulated blue-light receptor with distinct functions on flowering time in plants, and its function in maize remains unclear. In this study, we identified two FKF1 homologs in the maize genome, named ZmFKF1a and ZmFKF1b, and indicated that ZmFKF1a and ZmFKF1b independently regulate reproductive transition through interacting with ZmCONZ1 and ZmGI1 to increase the transcription levels of ZmCONZ1 and ZCN8. We demonstrated that ZmFKF1b underwent artificial selection during modern breeding in China probably due to its role in geographical adaptation. Furthermore, our data suggested that ZmFKF1bHap_C7 may be an elite allele, which increases the abundance of ZmCONZ1 mRNA more efficiently and adapt to a wider range of temperature zone than that of ZmFKF1bHap_Z58 to promote maize floral transition. It extends our understanding of the genetic diversity of maize flowering. This allele is expected to be introduced into tropical maize germplasm to enrich breeding resources and may improve the adaptability of maize at different climate zones, especially at temperate region.

Breeding maize of ideal plant architecture for high-density planting tolerance through modulating shade avoidance response and beyond

Maize is a major staple crop widely used as food, animal feed, and raw materials in industrial production. High-density planting is a major factor contributing to the continuous increase of maize yield. However, high planting density usually triggers a shade avoidance response and causes increased plant height and ear height, resulting in lodging and yield loss. Reduced plant height and ear height, more erect leaf angle, reduced tassel branch number, earlier flowering, and strong root system architecture are five key morphological traits required for maize adaption to high-density planting. In this review, we summarize recent advances in deciphering the genetic and molecular mechanisms of maize involved in response to high-density planting. We also discuss some strategies for breeding advanced maize cultivars with superior performance under high-density planting conditions.

Studies on Candidate Genes Related to Flowering Time in a Multiparent Population of Maize Derived from Tropical and Temperate Germplasm

A comprehensive study on maize flowering traits, focusing on the regulation of flowering time and the elucidation of molecular mechanisms underlying the genes controlling flowering, holds the potential to significantly enhance our understanding of the associated regulatory gene network. In this study, three tropical maize inbreds, CML384, CML171, and CML444, were used, along with a temperate maize variety, Shen137, as parental lines to cross with Ye107. The resulting F1s underwent seven consecutive generations of self-pollination through the single-seed descent (SSD) method to develop a multiparent population. To investigate the regulation of maize flowering time-related traits and to identify loci and candidate genes, a genome-wide association study (GWAS) was conducted. GWAS analysis identified 556 SNPs and 12 candidate genes that were significantly associated with flowering time-related traits. Additionally, an analysis of the effect of the estimated breeding values of the subpopulations on flowering time was conducted to further validate the findings of the present study. Collectively, this study offers valuable insights into novel candidate genes, contributing to an improved understanding of maize flowering time-related traits. This information holds practical significance for future maize breeding programs aimed at developing high-yielding hybrids.

QTL mapping for plant height and ear height using bi-parental immortalized heterozygous populations in maize

Introduction

Plant height (PH) and ear height (EH) are key plant architectural traits in maize, which will affect the photosynthetic efficiency, high plant density tolerance, suitability for mechanical harvesting.

Methods

QTL mapping were conducted for PH and EH using a recombinant inbred line (RIL) population and two corresponding immortalized backcross (IB) populations obtained from crosses between the RIL population and the two parental lines.

Results

A total of 17 and 15 QTL were detected in the RIL and IB populations, respectively. Two QTL, qPH1-1 (qEH1-1) and qPH1-2 (qEH1-4) in the RIL, were simultaneously identified for PH and EH. Combing reported genome-wide association and cloned PH-related genes, co-expression network analyses were constructed, then five candidate genes with high confidence in major QTL were identified including Zm00001d011117 and Zm00001d011108, whose homologs have been confirmed to play a role in determining PH in maize and soybean.

Discussion

QTL mapping used a immortalized backcross population is a new strategy. These identified genes in this study can provide new insights for improving the plant architecture in maize.

Striking a growth-defense balance: Stress regulators that function in maize development

Maize (Zea mays) cultivation is strongly affected by both abiotic and biotic stress, leading to reduced growth and productivity. It has recently become clear that regulators of plant stress responses, including the phytohormones abscisic acid (ABA), ethylene (ET), and jasmonic acid (JA), together with reactive oxygen species (ROS), shape plant growth and development. Beyond their well established functions in stress responses, these molecules play crucial roles in balancing growth and defense, which must be finely tuned to achieve high yields in crops while maintaining some level of defense. In this review, we provide an in-depth analysis of recent research on the developmental functions of stress regulators, focusing specifically on maize. By unraveling the contributions of these regulators to maize development, we present new avenues for enhancing maize cultivation and growth while highlighting the potential risks associated with manipulating stress regulators to enhance grain yields in the face of environmental challenges.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

A 49-bp deletion of PmAP2L results in a double flower phenotype in Prunus mume

The double flower is an important trait with substantial ornamental value. While mutations in PETALOSA TOE-type or AG (AGAMOUS) genes play a crucial role in enhancing petal number in ornamental plants, the complete mechanism underlying the formation of double flowers remains to be fully elucidated. Through the application of bulked segregant analysis (BSA), we identified a novel gene, APETALA2-like (PmAP2L), characterized by a 49-bp deletion in double-flowered Prunus mume. β-Glucuronidase (GUS) staining and luciferase reporter assays confirmed that the 49-bp deletion in PmAP2L reduced its binding with Pmu-miRNA172a. Phylogenetic analysis and microsynteny analysis suggested that PmAP2L was not a PETALOSA TOE-type gene, and it might be a new gene controlling the formation of double flower in P. mume. Subsequently, overexpression of PmAP2L-D in tobacco led to a significant rise in the number of stamens and the conversion of stamens to petals. Furthermore, silencing of the homologue of RC5G0530900 in rose significantly reduced the number of petals. Using transient gene expression in P. mume flower buds, we determined the functional differences between PmAP2L-D and PmAP2-S in controlling flower development. Meanwhile, DNA-affinity purification sequencing (DAP-seq), yeast hybrid assays and luciferase reporter assays indicated that PmAP2L negatively regulated the floral organ identity genes by forming a repressor complex with PmTPL and PmHDA6/19. Overall, these findings indicate that the variation in PmAP2L is associated with differences in the regulation of genes responsible for floral organ identity, providing new insights into the double-flower trait and double-flower breeding in plants.

Application of Single-Cell Assay for Transposase-Accessible Chromatin with High Throughput Sequencing in Plant Science: Advances, Technical Challenges, and Prospects

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Genome-wide DNase I-hypersensitive site assay reveals distinct genomic distributions and functional features of open chromatin in autopolyploid sugarcane

The characterization of cis-regulatory DNA elements (CREs) is essential for deciphering the regulation of gene expression in eukaryotes. Although there have been endeavors to identify CREs in plants, the properties of CREs in polyploid genomes are still largely unknown. Here, we conducted the genome-wide identification of DNase I-hypersensitive sites (DHSs) in leaf and stem tissues of the auto-octoploid species Saccharum officinarum. We revealed that DHSs showed highly similar distributions in the genomes of these two S. officinarum tissues. Notably, we observed that approximately 74% of DHSs were located in distal intergenic regions, suggesting considerable differences in the abundance of distal CREs between S. officinarum and other plants. Leaf- and stem-dependent transcriptional regulatory networks were also developed by mining the binding motifs of transcription factors (TFs) from tissue-specific DHSs. Four TEOSINTE BRANCHED 1, CYCLOIDEA, and PCF1 (TCP) TFs (TCP2, TCP4, TCP7, and TCP14) and two ethylene-responsive factors (ERFs) (ERF109 and ERF03) showed strong causal connections with short binding distances from each other, pointing to their possible roles in the regulatory networks of leaf and stem development. Through functional validation in transiently transgenic protoplasts, we isolate a set of tissue-specific promoters. Overall, the DHS maps presented here offer a global view of the potential transcriptional regulatory elements in polyploid sugarcane and can be expected to serve as a valuable resource for both transcriptional network elucidation and genome editing in sugarcane breeding.

Unlocking Nature's Clock: CRISPR Technology in Flowering Time Engineering

Flowering is a crucial process in the life cycle of most plants as it is essential for the reproductive success and genetic diversity of the species. There are situations in which breeders want to expedite, delay, or prevent flowering, for example, to shorten or prolong vegetative growth, to prevent unwanted pollination, to reduce the risk of diseases or pests, or to modify the plant's phenotypes. This review aims to provide an overview of the current state of knowledge to use CRISPR/Cas9, a powerful genome-editing technology to modify specific DNA sequences related to flowering induction. We discuss the underlying molecular mechanisms governing the regulation of the photoperiod, autonomous, vernalization, hormonal, sugar, aging, and temperature signal pathways regulating the flowering time. In addition, we are investigating the most effective strategies for nominating target genes. Furthermore, we have collected a dataset showing successful applications of CRISPR technology to accelerate flowering in several plant species from 2015 up to date. Finally, we explore the opportunities and challenges of using the potential of CRISPR technology in flowering time engineering.

A potential candidate gene associated with the angles of the ear leaf and the second leaf above the ear leaf in maize

BACKGROUND

Leaf angle is a key trait for maize plant architecture that plays a significant role in its morphological development, and ultimately impacting maize grain yield. Although many studies have been conducted on the association and localization of genes regulating leaf angle in maize, most of the candidate genes identified are associated with the regulation of ligule-ear development and phytohormone pathways, and only a few candidate genes have been reported to enhance the mechanical strength of leaf midrib and vascular tissues.

RESULTS

To address this gap, we conducted a genome-wide association study (GWAS) using the leaf angle phenotype and genotyping-by-sequencing data generated from three recombinant inbred line (RIL) populations of maize. Through GWAS analysis, we identified 156 SNPs significantly associated with the leaf angle trait and detected a total of 68 candidate genes located within 10 kb upstream and downstream of these individual SNPs. Among these candidate genes, Zm00001d045408, located on chromosome 9 emerged as a key gene controlling the angles of both the ear leaf and the second leaf above the ear leaf. Notably, this new gene's homolog in Arabidopsis promotes cell division and vascular tissue development. Further analysis revealed that a SNP transversion (G/T) at 7.536 kb downstream of the candidate gene Zm00001d045408 may have caused a reduction in leaf angles of the ear and the second leaf above the ear leaf. Our analysis of the 10 kb region downstream of this candidate gene revealed a 4.337 kb solo long-terminal reverse transcription transposon (solo LTR), located 3.112 kb downstream of Zm00001d045408, with the SNP located 87 bp upstream of the solo LTR.

CONCLUSIONS

In summary, we have identified a novel candidate gene, Zm00001d045408 and a solo LTR that are associated with the angles of both the ear leaf and the second leaf above the ear leaf. The future research holds great potential in exploring the precise role of newly identified candidate gene in leaf angle regulation. Functional characterization of this gene can help in gaining deeper insights into the complex genetic pathways underlying maize plant architecture.

Genetic analysis of phenotypic plasticity identifies BBX6 as the candidate gene for maize adaptation to temperate regions

Introduction

Climate changes pose a significant threat to crop adaptation and production. Dissecting the genetic basis of phenotypic plasticity and uncovering the responsiveness of regulatory genes to environmental factors can significantly contribute to the improvement of climate- resilience in crops.

Methods

We established a BC1F3:4 population using the elite inbred lines Zheng58 and PH4CV and evaluated plant height (PH) across four environments characterized by substantial variations in environmental factors. Then, we quantified the correlation between the environmental mean of PH (the mean performance in each environment) and the environmental parameters within a specific growth window. Furthermore, we performed GWAS analysis of phenotypic plasticity, and identified QTLs and candidate gene that respond to key environment index. After that, we constructed the coexpression network involving the candidate gene, and performed selective sweep analysis of the candidate gene.

Results

We found that the environmental parameters demonstrated substantial variation across the environments, and genotype by environment interaction contributed to the variations of PH. Then, we identified PTT(35-48) (PTT is the abbreviation for photothermal units), the mean PTT from 35 to 48 days after planting, as the pivotal environmental index that closely correlated with environmental mean of PH. Leveraging the slopes of the response of PH to both the environmental mean and PTT(35-48), we successfully pinpointed QTLs for phenotypic plasticity on chromosomes 1 and 2. Notably, the PH4CV genotypes at these two QTLs exhibited positive contributions to phenotypic plasticity. Furthermore, our analysis demonstrated a direct correlation between the additive effects of each QTL and PTT(35-48). By analyzing transcriptome data of the parental lines in two environments, we found that the 1009 genes responding to PTT(35-48) were enriched in the biological processes related to environmental sensitivity. BBX6 was the prime candidate gene among the 13 genes in the two QTL regions. The coexpression network of BBX6 contained other genes related to flowering time and photoperiod sensitivity. Our investigation, including selective sweep analysis and genetic differentiation analysis, suggested that BBX6 underwent selection during maize domestication.

Discussion

Th is research substantially advances our understanding of critical environmental factors influencing maize adaptation while simultaneously provides an invaluable gene resource for the development of climate-resilient maize hybrid varieties.

Genomic variations combined with epigenetic modifications rewire open chromatin in rice

Cis-regulatory elements (CREs) fine-tune gene transcription in eukaryotes. CREs with sequence variations play vital roles in driving plant or crop domestication. However, how global sequence and structural variations (SVs) are responsible for multilevel changes between indica and japonica rice (Oryza sativa) is still not fully elucidated. To address this, we conducted multiomic studies using MNase hypersensitivity sequencing (MH-seq) in combination with RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and bisulfite sequencing (BS-seq) between the japonica rice variety Nipponbare (NIP) and indica rice variety 93-11. We found that differential MNase hypersensitive sites (MHSs) exhibited some distinct intrinsic genomic sequence features between NIP and 93-11. Notably, through MHS-genome-wide association studies (GWAS) integration, we found that key sequence variations may be associated with differences of agronomic traits between NIP and 93-11, which is partly achieved by MHSs harboring CREs. In addition, SV-derived differential MHSs caused by transposable element (TE) insertion, especially by noncommon TEs among rice varieties, were associated with genes with distinct functions, indicating that TE-driven gene neo- or subfunctionalization is mediated by changes of chromatin openness. This study thus provides insights into how sequence and genomic SVs control agronomic traits of NIP and 93-11; it also provides genome-editing targets for molecular breeding aiming at improving favorable agronomic properties.

Identifying QTLs involved in hybrid performance and heterotic group complementarity: new GWAS models applied to factorial and admixed diallel maize hybrid panels

KEY MESSAGE

An original GWAS model integrating the ancestry of alleles was proposed and allowed the detection of background specific additive and dominance QTLs involved in heterotic group complementarity and hybrid performance. Maize genetic diversity is structured into genetic groups selected and improved relative to each other. This process increases group complementarity and differentiation over time and ensures that the hybrids produced from inter-group crosses exhibit high performances and heterosis. To identify loci involved in hybrid performance and heterotic group complementarity, we introduced an original association study model that disentangles allelic effects from the heterotic group origin of the alleles and compared it with a conventional additive/dominance model. This new model was applied on a factorial between Dent and Flint lines and a diallel between Dent-Flint admixed lines with two different layers of analysis: within each environment and in a multiple-environment context. We identified several strong additive QTLs for all traits, including some well-known additive QTLs for flowering time (in the region of Vgt1/2 on chromosome 8). Yield trait displayed significant non-additive effects in the diallel panel. Most of the detected Yield QTLs exhibited overdominance or, more likely, pseudo-overdominance effects. Apparent overdominance at these QTLs contributed to a part of the genetic group complementarity. The comparison between environments revealed a higher stability of additive QTL effects than non-additive ones. Several QTLs showed variations of effects according to the local heterotic group origin. We also revealed large chromosomic regions that display genetic group origin effects. Altogether, our results illustrate how admixed panels combined with dedicated GWAS modeling allow the identification of new QTLs that could not be revealed by a classical hybrid panel analyzed with traditional modeling.

AP2/EREBP Pathway Plays an Important Role in Chaling Wild Rice Tolerance to Cold Stress

Cold stress is the main factor limiting rice production and distribution. Chaling wild rice can survive in cold winters. AP2/EREBP is a known transcription factor family associated with abiotic stress. We identified the members of the AP2/EREBP transcription factor family in rice, maize, and Arabidopsis, and conducted collinearity analysis and gene family analysis. We used Affymetrix array technology to analyze the expression of AP2/EREBP family genes in Chaling wild rice and cultivated rice cultivar Pei'ai64S, which is sensitive to cold. According to the GeneChip results, the expression levels of AP2/EREBP genes in Chaling wild rice were different from those in Pei'ai64S; and the increase rate of 36 AP2/EREBP genes in Chaling wild rice was higher than that in Pei'ai64S. Meanwhile, the MYC elements in cultivated rice and Chaling wild rice for the Os01g49830, Os03g08470, and Os03g64260 genes had different promoter sequences, resulting in the high expression of these genes in Chaling wild rice under low-temperature conditions. Furthermore, we analyzed the upstream and downstream genes of the AP2/EREBP transcription factor family and studied the conservation of these genes. We found that the upstream transcription factors were more conserved, indicating that these upstream transcription factors may be more important in regulating cold stress. Meanwhile, we found the expression of AP2/EREBP pathway genes was significantly increased in recombinant inbred lines from Nipponbare crossing with Chaling wild rice, These results suggest that the AP2/EREBP signaling pathway plays an important role in Chaling wild rice tolerance to cold stress.

Phenomic data-driven biological prediction of maize through field-based high-throughput phenotyping integration with genomic data

High-throughput phenotyping (HTP) has expanded the dimensionality of data in plant research; however, HTP has resulted in few novel biological discoveries to date. Field-based HTP (FHTP), using small unoccupied aerial vehicles (UAVs) equipped with imaging sensors, can be deployed routinely to monitor segregating plant population interactions with the environment under biologically meaningful conditions. Here, flowering dates and plant height, important phenological fitness traits, were collected on 520 segregating maize recombinant inbred lines (RILs) in both irrigated and drought stress trials in 2018. Using UAV phenomic, single nucleotide polymorphism (SNP) genomic, as well as combined data, flowering times were predicted using several scenarios. Untested genotypes were predicted with 0.58, 0.59, and 0.41 prediction ability for anthesis, silking, and terminal plant height, respectively, using genomic data, but prediction ability increased to 0.77, 0.76, and 0.58 when phenomic and genomic data were used together. Using the phenomic data in a genome-wide association study, a heat-related candidate gene (GRMZM2G083810; hsp18f) was discovered using temporal reflectance phenotypes belonging to flowering times (both irrigated and drought) trials where heat stress also peaked. Thus, a relationship between plants and abiotic stresses belonging to a specific time of growth was revealed only through use of temporal phenomic data. Overall, this study showed that (i) it is possible to predict complex traits using high dimensional phenomic data between different environments, and (ii) temporal phenomic data can reveal a time-dependent association between genotypes and abiotic stresses, which can help understand mechanisms to develop resilient plants.

Enhancing adaptation of tropical maize to temperate environments using genomic selection

Tropical maize can be used to diversify the genetic base of temperate germplasm and help create climate-adapted cultivars. However, tropical maize is unadapted to temperate environments, in which sensitivities to long photoperiods and cooler temperatures result in severely delayed flowering times, developmental defects, and little to no yield. Overcoming this maladaptive syndrome can require a decade of phenotypic selection in a targeted, temperate environment. To accelerate the incorporation of tropical diversity in temperate breeding pools, we tested if an additional generation of genomic selection can be used in an off-season nursery where phenotypic selection is not very effective. Prediction models were trained using flowering time recorded on random individuals in separate lineages of a heterogenous population grown at two northern U.S. latitudes. Direct phenotypic selection and genomic prediction model training was performed within each target environment and lineage, followed by genomic prediction of random intermated progenies in the off-season nursery. Performance of genomic prediction models was evaluated on self-fertilized progenies of prediction candidates grown in both target locations in the following summer season. Prediction abilities ranged from 0.30 to 0.40 among populations and evaluation environments. Prediction models with varying marker effect distributions or spatial field effects had similar accuracies. Our results suggest that genomic selection in a single off-season generation could increase genetic gains for flowering time by more than 50% compared to direct selection in summer seasons only, reducing the time required to change the population mean to an acceptably adapted flowering time by about one-third to one-half.

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.

cis-Regulatory Elements in Plant Development, Adaptation, and Evolution

cis-Regulatory elements encode the genomic blueprints that ensure the proper spatiotemporal patterning of gene expression necessary for appropriate development and responses to the environment. Accumulating evidence implicates changes to gene expression as a major source of phenotypic novelty in eukaryotes, including acute phenotypes such as disease and cancer in mammals. Moreover, genetic and epigenetic variation affecting cis-regulatory sequences over longer evolutionary timescales has become a recurring theme in studies of morphological divergence and local adaptation. Here, we discuss the functions of and methods used to identify various classes of cis-regulatory elements, as well as their role in plant development and response to the environment. We highlight opportunities to exploit cis-regulatory variants underlying plant development and environmental responses for crop improvement efforts. Although a comprehensive understanding of cis-regulatory mechanisms in plants has lagged behind that in animals, we showcase several breakthrough findings that have profoundly influenced plant biology and shaped the overall understanding of transcriptional regulation in eukaryotes.

Effects of the quantitative trait locus qPss3 on inhibition of photoperiod sensitivity and resistance to stalk rot disease in maize

KEY MESSAGE

We identified a quantitative trait locus, qPss3, and fine-mapped the causal locus to a 120-kb interval in maize. This locus inhibits the photoperiod sensitivity caused by ZmCCT9 and ZmCCT10, resulting in earlier flowering by 2 ~ 4 days without reduction in stalk-rot resistance in certain genotypes. Photoperiod sensitivity is a key factor affecting the adaptation of maize (Zea mays L.) to high-latitude growing areas. Although many genes associated with flowering time have been identified in maize, no gene that inhibits photoperiod sensitivity has been reported. In our previous study, we detected large differences in photoperiod sensitivity among maize inbred lines with the same photoperiod-sensitive allele at the ZmCCT10 locus. Here, we used two segregating populations with the same genetic backgrounds but different ZmCCT10 alleles to perform quantitative trait locus (QTL) analysis. We identified a unique QTL, qPss3, on chromosome 3 in the population carrying the sensitive ZmCCT10 allele. After sequential fine-mapping, we eventually delimited qPss3 to an interval of ~ 120 kb. qPss3 behaved as a dominant locus and caused earlier flowering by 2-4 days via inhibiting ZmCCT10-induced photoperiod sensitivity under long-day conditions. qPss3 also inhibited the photoperiod sensitivity induced by another flowering-related gene, ZmCCT9. For application in agriculture, an F₁ hybrid heterozygous at both qPss3 and ZmCCT10 loci constitutes an optimal allele combination, showing high resistance to stalk rot without a significant delay in flowering time. Moreover, qPss3 is of great value in regulating the flowering time of tropical maize grown at high-latitude regions.

Tomato APETALA2 family member SlTOE1 regulates inflorescence branching by repressing SISTER OF TM3

Inflorescence architecture directly impacts yield potential in most crops. As a model of sympodial plants, tomato (Solanum lycopersicum) inflorescence exhibits highly structural plasticity. However, the genetic regulatory network of inflorescence architecture in tomato remains unclear. Here, we investigated a modulator of inflorescence branching in tomato, TARGET OF EAT1 (SlTOE1), an APETALA2 (AP2) family member found to be predominantly expressed in the floral meristem (FM) of tomato. sltoe1 knockout mutants displayed highly branched inflorescences and defective floral organs. Transcriptome analysis revealed that SISTER OF TM3 (STM3) and certain floral development-related genes were upregulated in the flower meristem of sltoe1. SlTOE1 could directly bind the promoters of STM3 and Tomato MADS-box gene 3 (TM3) to repress their transcription. Simultaneous mutation of STM3 and TM3 partially restored the inflorescence branching of the sltoe1cr mutants, suggesting that SlTOE1 regulates inflorescence development, at least in part through an SlTOE1STM3/TM3 module. Genetic analysis showed that SlTOE1 and ENHANCER OF JOINTLESS 2 (EJ2) additively regulate tomato inflorescence branching; their double mutants showed more extensive inflorescence branching. Our findings uncover a pathway controlling tomato inflorescence branching and offer deeper insight into the functions of AP2 subfamily members.

Genome-Wide Association Study Identified Novel SNPs Associated with Chlorophyll Content in Maize

Chlorophyll is an essential component that captures light energy to drive photosynthesis. Chlorophyll content can affect photosynthetic activity and thus yield. Therefore, mining candidate genes of chlorophyll content will help increase maize production. Here, we performed a genome-wide association study (GWAS) on chlorophyll content and its dynamic changes in 378 maize inbred lines with extensive natural variation. Our phenotypic assessment showed that chlorophyll content and its dynamic changes were natural variations with a moderate genetic level of 0.66/0.67. A total of 19 single-nucleotide polymorphisms (SNPs) were found associated with 76 candidate genes, of which one SNP, 2376873-7-G, co-localized in chlorophyll content and area under the chlorophyll content curve (AUCCC). Zm00001d026568 and Zm00001d026569 were highly associated with SNP 2376873-7-G and encoded pentatricopeptide repeat-containing protein and chloroplastic palmitoyl-acyl carrier protein thioesterase, respectively. As expected, higher expression levels of these two genes are associated with higher chlorophyll contents. These results provide a certain experimental basis for discovering the candidate genes of chlorophyll content and finally provide new insights for cultivating high-yield and excellent maize suitable for planting environment.

Comprehensive analysis and expression profiles of the AP2/ERF gene family during spring bud break in tea plant (Camellia sinensis)

BACKGROUND

AP2/ERF transcription factors (AP2/ERFs) are important regulators of plant physiological and biochemical metabolism. Evidence suggests that AP2/ERFs may be involved in the regulation of bud break in woody perennials. Green tea is economically vital in China, and its production value is significantly affected by the time of spring bud break of tea plant. However, the relationship between AP2/ERFs in tea plant and spring bud break remains largely unknown.

RESULTS

A total of 178 AP2/ERF genes (CsAP2/ERFs) were identified in the genome of tea plant. Based on the phylogenetic analysis, these genes could be classified into five subfamilies. The analysis of gene duplication events demonstrated that whole genome duplication (WGD) or segmental duplication was the primary way of CsAP2/ERFs amplification. According to the result of the Ka/Ks value calculation, purification selection dominated the evolution of CsAP2/ERFs. Furthermore, gene composition and structure analyses of CsAP2/ERFs indicated that different subfamilies contained a variety of gene structures and conserved motifs, potentially resulting in functional differences among five subfamilies. The promoters of CsAP2/ERFs also contained various signal-sensing elements, such as abscisic acid responsive elements, light responsive elements and low temperature responsive elements. The evidence presented here offers a theoretical foundation for the diverse functions of CsAP2/ERFs. Additionally, the expressions of CsAP2/ERFs during spring bud break of tea plant were analyzed by RNA-seq and grouped into clusters A-F according to their expression patterns. The gene expression changes in clusters A and B were more synchronized with the spring bud break of tea plant. Moreover, several potential correlation genes, such as D-type cyclin genes, were screened out through weighted correlation network analysis (WGCNA). Temperature and light treatment experiments individually identified nine candidate CsAP2/ERFs that may be related to the spring bud break of tea plant.

CONCLUSIONS

This study provides new evidence for role of the CsAP2/ERFs in the spring bud break of tea plant, establishes a theoretical foundation for analyzing the molecular mechanism of the spring bud break of tea plant, and contributes to the improvement of tea cultivars.

Genomic-regions associated with cold stress tolerance in Asia-adapted tropical maize germplasm

Maize is gaining impetus in non-traditional and non-conventional seasons such as off-season, primarily due to higher demand and economic returns. Maize varieties directed for growing in the winter season of South Asia must have cold resilience as an important trait due to the low prevailing temperatures and frequent cold snaps observed during this season in most parts of the lowland tropics of Asia. The current study involved screening of a panel of advanced tropically adapted maize lines to cold stress during vegetative and flowering stage under field conditions. A suite of significant genomic loci (28) associated with grain yield along and agronomic traits such as flowering (15) and plant height (6) under cold stress environments. The haplotype regression revealed 6 significant haplotype blocks for grain yield under cold stress across the test environments. Haplotype blocks particularly on chromosomes 5 (bin5.07), 6 (bin6.02), and 9 (9.03) co-located to regions/bins that have been identified to contain candidate genes involved in membrane transport system that would provide essential tolerance to the plant. The regions on chromosome 1 (bin1.04), 2 (bin 2.07), 3 (bin 3.05-3.06), 5 (bin5.03), 8 (bin8.05-8.06) also harboured significant SNPs for the other agronomic traits. In addition, the study also looked at the plausibility of identifying tropically adapted maize lines from the working germplasm with cold resilience across growth stages and identified four lines that could be used as breeding starts in the tropical maize breeding pipelines.

Environment-specific selection alters flowering-time plasticity and results in pervasive pleiotropic responses in maize

Crop genetic diversity for climate adaptations is globally partitioned. We performed experimental evolution in maize to understand the response to selection and how plant germplasm can be moved across geographical zones. Initialized with a common population of tropical origin, artificial selection on flowering time was performed for two generations at eight field sites spanning 25° latitude, a 2800 km transect. We then jointly tested all selection lineages across the original sites of selection, for the target trait and 23 other traits. Modeling intergenerational shifts in a physiological reaction norm revealed separate components for flowering-time plasticity. Generalized and local modes of selection altered the plasticity of each lineage, leading to a latitudinal pattern in the responses to selection that were strongly driven by photoperiod. This transformation led to widespread changes in developmental, architectural, and yield traits, expressed collectively in an environment-dependent manner. Furthermore, selection for flowering time alone alleviated a maladaptive syndrome and improved yields for tropical maize in the temperate zone. Our findings show how phenotypic selection can rapidly shift the flowering phenology and plasticity of maize. They also demonstrate that selecting crops to local conditions can accelerate adaptation to climate change.

A role for heritable transcriptomic variation in maize adaptation to temperate environments

Background

Transcription bridges genetic information and phenotypes. Here, we evaluated how changes in transcriptional regulation enable maize (Zea mays), a crop originally domesticated in the tropics, to adapt to temperate environments.

Result

We generated 572 unique RNA-seq datasets from the roots of 340 maize genotypes. Genes involved in core processes such as cell division, chromosome organization and cytoskeleton organization showed lower heritability of gene expression, while genes involved in anti-oxidation activity exhibited higher expression heritability. An expression genome-wide association study (eGWAS) identified 19,602 expression quantitative trait loci (eQTLs) associated with the expression of 11,444 genes. A GWAS for alternative splicing identified 49,897 splicing QTLs (sQTLs) for 7614 genes. Genes harboring both cis-eQTLs and cis-sQTLs in linkage disequilibrium were disproportionately likely to encode transcription factors or were annotated as responding to one or more stresses. Independent component analysis of gene expression data identified loci regulating co-expression modules involved in oxidation reduction, response to water deprivation, plastid biogenesis, protein biogenesis, and plant-pathogen interaction. Several genes involved in cell proliferation, flower development, DNA replication, and gene silencing showed lower gene expression variation explained by genetic factors between temperate and tropical maize lines. A GWAS of 27 previously published phenotypes identified several candidate genes overlapping with genomic intervals showing signatures of selection during adaptation to temperate environments.

Conclusion

Our results illustrate how maize transcriptional regulatory networks enable changes in transcriptional regulation to adapt to temperate regions.

Supplementary information

The online version contains supplementary material available at 10.1186/s13059-023-02891-3.