Improved maize reference genome with single-molecule technologies

The Emergence of Nanofluidics for Single-Biomolecule Manipulation and Sensing

Driven by recent advancements in nanofabrication techniques, single-molecule sensing and manipulations in nanofluidic devices are rapidly evolving. These sophisticated biosensors have already had significant impacts on basic research as well as on applications in molecular diagnostics. The nanoscale dimensions of these devices introduce new physical phenomena by confining the biomolecules in at least one dimension, creating effects such as biopolymer linearization, stretching, and separation by mass that are utilized to enhance the biomolecule sensing resolutions. At the same time, the suppressed diffusional motion allows for better single-molecule SNR (signal-to-noise ratio) sensing over time. In particular, nanofluidic devices based on nanochannels have been established as promising technologies for the linearization of ultralong genomic DNA molecules and for optical genome mapping, opening a window to directly observe and infer genome organization. More recently, nanochannels have shown promising capabilities for single-molecule protein sizing, separation, and identification. Consequently, this technology is attracting remarkable interest for applications in single-molecule proteomics. In this review, we discuss the recent advancements of nanochannel-based technologies, focusing on their applications for single-molecule sensing and the characterization of a wide range of biomolecules.

Unveiling the heterosis pattern of modern maize breeding in Southwest China through population structure and genetic diversity analysis

Maize (Zea mays L.) is an important food crop throughout the world and is also one of the earliest crops to use heterosis. In this study, we evaluated the genetic diversity, population structure, and selective sweep of 100 elite inbred maize lines collected from the current breeding program in Sichuan province, Southwest China, using 5,261,175 high-quality single nucleotide polymorphisms (SNPs). We discovered an abundance of genetic diversities and classified them into four groups. By combining kinship relationships, these groups were further divided into Tropic-local A, Improved-tropic, Tropic-local B, and Improved-local. Genomic differentiation was assessed using Fst values (0.21-0.44) as well as genetic diversity (π = 6.07 × 10-4 - 6.61 × 10-4). We generated 900 (90 × 10) hybrids using 90 and 10 inbred maize lines from 100 diverse maize germplasms. All hybrids were evaluated for 10 traits in three replicate tests across two locations. We found that the patterns of G1 × G3, G1 × G4, G2 × G3, and G3 × G4 exhibited significant heterosis in yield-related traits and have been used in commercial breeding. In addition, we also explored the relationship between 10 traits of hybrid offspring and the number of heterozygous SNP. Under most heterosis modes, the best linear unbiased estimation (BLUE) value of the trait was highly consistent with the trend of deleterious SNPs, but there was a deviation in the G1 × G3 mode. Taken together, the results provide insight into the utilization of the current maize germplasm in Sichuan province to improve hybrid breeding.

Integrated transcriptome and metabolome analysis reveals the impacts of prolonged light exposure on starch and protein content in maize kernels

BACKGROUND

The light environment significantly influences crop growth, development, quality, and yield, particularly in controlled-environment agriculture. Recent advances in artificial lighting technology have allowed growers to precisely control the light environment in terms of duration, spectrum, and intensity. Starch and protein are the most significant nutritional constituents of maize kernels. However, little is known about the effects of the light environment on starch and protein content in maize kernels. Therefore, we investigated the effects of natural light and supplemental exposure to blue (B), far-red (FR), and red (R) light on starch and protein content in kernels of the inbred maize line B73.

RESULTS

Exposure to supplemental B, FR, or R light resulted in significant increases in starch content but decreases in protein content. Notably, protein content was lowest under B light. Substantial proportions of genes (5.03-75.23%) and metabolites (46.89-85.64%) were regulated by different wavelengths of light. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, as well as weighted gene co-expression network analysis (WGCNA), revealed that differentially expressed genes (DEGs) under B, FR, and R light were involved in pathways related to starch and protein synthesis. KEGG metabolomic analysis showed that differentially abundant metabolites (DAMs) were primarily associated with histidine, D-amino acid, cysteine, and methionine metabolism. Nine DEGs related to starch synthesis were identified as potential candidates for investigating the effects of light quality on starch synthesis, and 14 DEGs related to protein synthesis provided evidence for the influence of light quality on protein synthesis in maize.

CONCLUSIONS

This study identified the regulatory network governing starch and protein content in B73 maize kernels under different light conditions, contributing to a deeper understanding of how light quality affects the nutritional components of maize kernels.

Genome-wide comparison reveals large structural variants in cassava landraces

BACKGROUND

Structural variants (SVs) are critical for plant genomic diversity and phenotypic variation. This study investigates a large, 9.7 Mbp highly repetitive segment on chromosome 12 of TMEB117, a region not previously characterized in cassava (Manihot esculenta Crantz). We aim to explore its presence and variability across multiple cassava landraces, providing insights into its genomic significance and potential implications.

RESULTS

We validated the presence of the 9.7 Mbp segment in the TMEB117 genome, distinguishing it from other published cassava genome assemblies. By mapping short-read sequencing data from 16 cassava landraces to TMEB117 chromosome 12, we observed variability in read mapping, suggesting that while all genotypes contain the insertion region, some exhibit missing segments or sequence differences. Further analysis revealed two unique genes associated with deacetylase activity, HDA14 and SRT2, within the insertion. Additionally, the MUDR-Mutator transposable element was significantly overrepresented in this region.

CONCLUSIONS

This study uncovers a large structural variant in the TMEB117 cassava genome, highlighting its variability among different genotypes. The enrichment of HDA14 and SRT2 genes and the MUDR-Mutator elements within the insertion suggests potential functional significance, though further research is needed to explore this. These findings provide important insights into the role of structural variations in shaping cassava genomic diversity.

Large-scale metabolomic landscape of edible maize reveals convergent changes in metabolite differentiation and facilitates its breeding improvement

Edible maize is an important food crop that provides energy and nutrients to meet human health and nutritional requirements. However, how environmental pressures and human activity have shaped the metabolome of edible maize remains unclear. In this study, we collected 452 diverse edible maize accessions worldwide, including waxy, sweet, and field maize. A total of 3020 non-redundant metabolites, including 802 annotated metabolites, were identified using a two-step optimized approach, which generated the most comprehensive annotated metabolite dataset in plants to date. Although specific metabolite differentiation was detected between field and sweet maize and between field and waxy maize, convergent metabolite differentiation was the dominant pattern. We identified hub genes in all metabolite classes by hotspot analysis in a metabolite genome-wide association study. Seventeen and 15 hub genes were selected as the key differentiation genes for flavonoids and lipids, respectively. Surprisingly, almost all of these genes were under diversifying selection, suggesting that diversifying selection was the main genetic mechanism of convergent metabolic differentiation. Further genetic and molecular studies revealed the roles and genetic diversifying selection mechanisms of ZmGPAT11 in convergent metabolite differentiation in the lipid pathway. On the basis of our research, we established the first edible maize metabolome database, EMMDB (https://www.maizemdb.site/home/). We successfully used EMMDB for precision improvement of nutritional and flavor traits and bred the elite inbred line 6644_2, with greatly increased contents of flavonoids, lysophosphatidylcholines, lysophosphatidylethanolamines, and vitamins. Collectively, our study sheds light on the genetic mechanisms of metabolite differentiation in edible maize and provides a database for breeding improvement of flavor and nutritional traits in edible maize by metabolome precision design.

ZmSSRP1 facilitates the progression of RNA polymerase II and is essential for kernel development in maize

Transcript elongation controlled by RNA polymerase II (RNAP II) represents a key regulatory event in numerous cellular processes. However, the precise mechanisms underlying the regulation of RNAP II distribution and progression in plants remain largely elusive. Here, we positionally cloned the causal mutation in the defective kernel 59 (dek59) maize (Zea mays) mutant and demonstrated that Dek59 encodes Structure-Specific Recognition Protein 1 (ZmSSRP1), a subunit of the FAcilitates Chromatin Transcription (FACT) complex that regulates RNAP II. Using genome-wide mapping assays, we determined that ZmSSRP1-binding sites co-localize with those of RNAP II phosphorylated at its serine 2 residue (Ser2P) and are highly enriched within actively transcribed genes. Mutation of ZmSSRP1 resulted in Ser2P accumulation around the +1 nucleosome of genes, affecting gene expression in a gene length-dependent manner. The reduced amount of RNAP II in the dek59 mutant was rescued to wild-type-like levels by inhibiting the proteasome, indicating that arrested RNAP II degradation is proteasome-dependent. These findings reveal the indispensable role of ZmSSRP1 in regulating RNAP II-mediated transcription, which is critical for the proper expression of thousands of genes during maize seed development.

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.

Creating a gene-indexed EMS mutation library of Zheng58 for improving maize genetics research

KEY MESSAGE

We created the Zheng58 EMS Mutation Library, which provides a valuable resource for future investigations into the functions and roles of specific genes in maize physiology and development. Understanding the genetic variations present within maize is fundamental to improving maize breeding programs and developing crops with desirable traits. In this work, the Gene-Indexed Ethyl Methanesulfonate (EMS) Mutation Library of Zheng58 in maize was created to accelerate maize genetics research. By chemically inducing mutations in the Zheng58 maize inbred line using (EMS), 426 M2 EMS lines were generated, of which 185 exhibited heritable phenotypic changes. Coupling with high throughput sequencing techniques, over two million mutations, encompassing single nucleotide polymorphisms (SNPs) and small insertions and deletions (InDels), were subsequently identified. Functional annotation of mutation sites further indicated that a significant number of mutations influences important cellular processes, including translation termination (8,279), splice site disruption leading to mis-splicing (9,504), and missense mutations affecting protein-coding sequences (52,494). Therefore, the Gene-Indexed EMS Mutation Library of maize Zheng58 provides a valuable resource for future investigations into the functions and roles of specific genes in maize physiology and development.

Identification of the maize drought-resistant gene Zinc-finger Inflorescence Meristem 23 through high-resolution temporal transcriptome analysis

Drought is a major abiotic stress that significantly limits maize productivity. However, previous transcriptomic studies with limited time-point sampling have hindered the construction of robust co-expression networks, making it challenging to identify reliable hub genes involved in drought tolerance. To overcome this limitation, we generated a high-temporal-resolution transcriptome dataset spanning 108 time points from maize seedlings subjected to two consecutive rounds of drought and re-watering treatments. A total of 8477 drought-responsive genes (DRGs) were identified by comparing drought-stressed and well-watered controls. Using weighted gene co-expression network analysis (WGCNA), we constructed 17 co-expression modules, of which 8 were strongly associated with drought stress responses and collectively contained 353 hub genes. Among them, we validated the drought resistance functions of ZmCPK35, a known drought-responsive gene, and Zinc-finger Inflorescence Meristem 23 (ZmZIM23), a newly identified drought-regulatory gene, within the M10 module. Functional analysis revealed that ZmZIM23 enhances drought tolerance by improving water-use efficiency, reducing transpiration rates, and promoting biomass accumulation. Furthermore, yeast one-hybrid (Y1H) and dual-luciferase (LUC) assays demonstrated that ZmWRKY40, another M10 module member, transcriptionally regulates both ZmZIM23 and ZmCPK35. By integrating high-resolution transcriptomic data with co-expression network analyses, this study unveils key drought-responsive regulatory networks in maize and identifies novel candidate genes for improving drought tolerance. These findings provide valuable insights into the genetic foundation of drought adaptation and offer potential targets for the development of drought-resistant maize cultivars.

Total tocopherol levels in maize grain depend on chlorophyll biosynthesis within the embryo

BACKGROUND

Tocopherols are a class of lipid-soluble compounds that have multiple functional roles in plants and exhibit vitamin E activity, an essential nutrient for human and animal health. The tocopherol biosynthetic pathway is conserved across the plant kingdom, but source of the key tocopherol pathway precursor, phytol, is unclear. Two protochlorophyllide reductases (POR1 and POR2) were previously identified as loci controlling the natural variation of total tocopherols in maize grain, a non-photosynthetic tissue. POR1 and POR2 are key genes in chlorophyll biosynthesis yet the contribution of the chlorophyll biosynthetic pathway to tocopherol biosynthesis is still not understood.

RESULTS

We took two approaches to alter the activity of these two POR genes within kernel tissue, physiological treatments and CRISPR/Cas9-mediated knockouts, to determine the role of chlorophyll biosynthesis for tocopherol content. Since light is required for POR enzymatic activity, we imposed a dark treatment on developing kernels, which reduced chlorophyll a and tocopherols levels in embryo tissue by 92-99% and 87-90%, respectively, compared to the light treatment. In CRISPR/Cas9-mediated knockouts, the levels of chlorophyll a and tocopherols in embryos of the por1 por2 double homozygous mutant were reduced by 98-100% and 76-83%, respectively, compared to WT.

CONCLUSION

These findings demonstrate that tocopherol synthesis in maize grain depends almost entirely on phytol derived from chlorophyll biosynthesis within the embryo. POR1 and POR2 activity play crucial roles in chlorophyll biosynthesis, underscoring the importance of POR alleles and their activity in the biofortification of vitamin E levels in non-photosynthetic grain of maize.

Regulatory variation controlling architectural pleiotropy in maize

An early event in plant organogenesis is establishment of a boundary between the stem cell containing meristem and differentiating lateral organ. In maize (Zea mays), evidence suggests a common gene network functions at boundaries of distinct organs and contributes to pleiotropy between leaf angle and tassel branch number, two agronomic traits. To uncover regulatory variation at the nexus of these two traits, we use regulatory network topologies derived from specific developmental contexts to guide multivariate genome-wide association analyses. In addition to defining network plasticity around core pleiotropic loci, we identify new transcription factors that contribute to phenotypic variation in canopy architecture, and structural variation that contributes to cis-regulatory control of pleiotropy between tassel branching and leaf angle across maize diversity. Results demonstrate the power of informing statistical genetics with context-specific developmental networks to pinpoint pleiotropic loci and their cis-regulatory components, which can be used to fine-tune plant architecture for crop improvement.

Comparative transcriptome analysis reveals potential regulatory genes involved in the development and strength formation of maize stalks

BACKGROUND

Stalk strength is a critical trait in maize that influences plant architecture, lodging resistance and grain yield. The developmental stage of maize, spanning from the vegetative stage to the reproductive stage, is critical for determining stalk strength. However, the dynamics of the genetic control of this trait remains unclear.

RESULTS

Here, we report a temporal resolution study of the maize stalk transcriptome in one tropical line and one non-stiff-stalk line using 53 transcriptomes collected covering V7 (seventh leaf stage) through silking stage. The time-course transcriptomes were categorized into four phases corresponding to stalk early development, stalk early elongation, stalk late elongation, and stalk maturation. Fuzzy c-means clustering and Gene Ontology (GO) analyses elucidated the chronological sequence of events that occur at four phases of stalk development. Gene Ontology analysis suggests that active cell division occurs in the stalk during Phase I. During Phase II, processes such as cell wall extension, lignin deposition, and vascular cell development are active. In Phase III, lignin metabolic process, secondary cell wall biogenesis, xylan biosynthesis process, cell wall biogenesis, and polysaccharide biosynthetic process contribute to cell wall strengthening. Defense responses, abiotic stresses, and transport of necessary nutrients or substances are active engaged during Phase IV. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the two maize lines presented significant gene expression differences in the phenylpropanoid biosynthesis pathway and the flavonoid biosynthesis pathway. Certain differentially expressed genes (DEGs) encoding transcription factors, especially those in the NAC and MYB families, may be involved in stalk development. In addition, six potential regulatory genes associated with stalk strength were identified through weighted gene co-expression network analysis (WGCNA).

CONCLUSION

The data set provides a high temporal-resolution atlas of gene expression during maize stalk development. These phase-specific genes, differentially expressed genes, and potential regulatory genes reported in this study provide important resources for further studies to elucidate the genetic control of stalk development and stalk strength formation in maize.

The maize mTERF18 regulates transcriptional termination of the mitochondrial nad6 gene and is essential for kernel development

Mitochondria are semi-autonomous organelles present in eukaryotic cells, containing their own genome and transcriptional machinery. However, their functions are intricately linked to proteins encoded by the nuclear genome. Mitochondrial transcription termination factors (mTERFs) are nucleic acid-binding proteins involved in RNA splicing and transcription termination within plant mitochondria and chloroplasts. Despite their recognized importance, the specific roles of mTERF proteins in maize remain largely unexplored. Here, we clone and functionally characterize the maize mTERF18 gene. Our findings reveal that mTERF18 mutations lead to severely undifferentiated embryos, resulting in abortive phenotypes. Early kernel exhibits abnormal basal endosperm transfer layer and a significant reduction in both starch and protein accumulation in mterf18. We identify the mTERF18 gene through mapping-based cloning and validate this gene through allelic tests. mTERF18 is widely expressed across various maize tissues and encodes a highly conserved mitochondrial protein. Transcriptome data reveal that mTERF18 mutations disrupt transcriptional termination of the nad6 gene, leading to undetectable levels of Nad6 protein and reduced complex I assembly and activity. Furthermore, transmission electron microscopy observation of mterf18 endosperm uncover severe mitochondrial defects. Collectively, these findings highlight the critical role of mTERF18 in mitochondrial gene transcription termination and its pivotal impact on maize kernel development.

Time-course transcriptomic analysis reveals transcription factors involved in modulating nitrogen sensibility in maize

Nitrogen (N) serves both as a vital macronutrient and a signaling molecule for plants. Unveiling key regulators involved in N metabolism helps dissect the mechanisms underlying N metabolism, which is essential for developing maize with high N use efficiency. Two maize lines, B73 and Ki11, show differential chlorate and low-N tolerance. Time-course transcriptomic analysis reveals that the expression of N utilization genes (NUGs) in B73 and Ki11 have distinct responsive patterns to nitrate variation. By the coexpression networks, significant differences in the number of N response modules and regulatory networks of transcription factors (TFs) are revealed between B73 and Ki11. There are 23 unique TFs in B73 and 41 unique TFs in Ki11. MADS26 is a unique TF in the B73 N response network, with different expression levels and N response patterns in B73 and Ki11. Overexpression of MADS26 enhances the sensitivity to chlorate and the utilization of nitrate in maize, at least partially explaining the differential chlorate tolerance and low-N sensitivity between B73 and Ki11. The findings in this work provide unique insights and promising candidates for maize breeding to reduce unnecessary N overuse.

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.

Time-resolved oxidative signal convergence across the algae-embryophyte divide

The earliest land plants faced a significant challenge in adapting to environmental stressors. Stress on land is unique in its dynamics, entailing swift and drastic changes in light and temperature. While we know that land plants share with their closest streptophyte algal relatives key components of the genetic makeup for dynamic stress responses, their concerted action is little understood. Here, we combine time-course stress profiling using photophysiology, transcriptomics on 2.7 Tbp of data, and metabolite profiling analyses on 270 distinct samples, to study stress kinetics across three 600-million-year-divergent streptophytes. Through co-expression analysis and Granger causal inference we predict a gene regulatory network that retraces a web of ancient signal convergences at ethylene signaling components, osmosensors, and chains of major kinases. These kinase hubs already integrated diverse environmental inputs since before the dawn of plants on land.

A candidate association study of transcription factors in maize revealed the ZmPLATZ15-ZmEREB200 module as a key regulator of waterlogging tolerance at the seedling stage

Soil waterlogging is a major abiotic stress that severely impairs maize growth and development by inducing hypoxic conditions that disrupt essential physiological processes. Transcription factors (TFs) play crucial roles in modulating plant responses to waterlogging stress by regulating the expression of stress-related genes that enhance or diminish stress tolerance. In this study, we conducted an association analysis to identify 11 TFs closely associated with waterlogging stress in maize. Notably, the PLATZ family emerged as a novel and significant contributor to waterlogging stress. Overexpression of ZmPLATZ15 resulted in increased sensitivity to waterlogging at maize seedlings. Conversely, ZmEREB200, a member of the maize Group VII ERF (ZmERFVII) family, was significantly downregulated in the ZmPLATZ15 overexpression lines under waterlogging stress. Promoter analysis revealed that ZmPLATZ15 regulates ZmEREB200 by binding to the A/T-rich motifs in the ZmEREB200 promoter. Interestingly, overexpression of ZmEREB200 was found to enhance waterlogging tolerance at maize seedlings. To further elucidate their roles, we analyzed the transcriptomic profiles of ZmPLATZ15 and ZmEREB200 overexpression lines under waterlogging stress. The overlapping differentially expressed genes in both ZmPLATZ15 and ZmEREB200 overexpression lines were significantly enriched in pathways associated with redox balance and salicylic acid metabolism, both of which are crucial for modulating waterlogging tolerance at maize seedlings. Metabolomic analysis revealed that antioxidant enzyme activity, salicylic acid, and glutathione levels were decreased in OE-ZmPLATZ15, while these metabolites were significantly increased in OE-ZmEREB200. These contrasting metabolic responses in overexpression lines may underlie their different tolerances to waterlogging stress. Our findings provide valuable insights into the regulatory mechanisms underlying maize's response to waterlogging stress and highlight the potential of TFs as tools for developing maize varieties with enhanced waterlogging tolerance.

Near-complete assembly and comprehensive annotation of the wheat Chinese Spring genome

A complete reference genome assembly is crucial for biological research and genetic improvement. Owing to its large size and highly repetitive nature, there are numerous gaps in the globally used wheat Chinese Spring (CS) genome assembly. In this study, we generated a 14.46 Gb near-complete assembly of the CS genome, with a contig N50 of over 266 Mb and an overall base accuracy of 99.9963%. Among the 290 gaps that remained (26, 257, and 7 gaps from the A, B, and D subgenomes, respectively), 278 were extremely high-copy tandem repeats, whereas the remaining 12 were transposable-element-associated gaps. Four chromosome assemblies were completely gap-free, including chr1D, chr3D, chr4D, and chr5D. Extensive annotation of the near-complete genome revealed 151 405 high-confidence genes, of which 59 180 were newly annotated, including 7602 newly assembled genes. Except for the centromere of chr1B, which has a gap associated with superlong GAA repeat arrays, the centromeric sequences of all of the remaining 20 chromosomes were completely assembled. Our near-complete assembly revealed that the extent of tandem repeats, such as simple-sequence repeats, was highly uneven among different subgenomes. Similarly, the repeat compositions of the centromeres also varied among the three subgenomes. With the genome sequences of all six types of seed storage proteins (SSPs) fully assembled, the expression of ω-gliadin was found to be contributed entirely by the B subgenome, whereas the expression of the other five types of SSPs was most abundant from the D subgenome. The near-complete CS genome will serve as a valuable resource for genomic and functional genomic research and breeding of wheat as well as its related species.

A unified VCF dataset from nearly 1,500 diverse maize accessions and resources to explore the genomic landscape of maize

Efforts to capture and analyze maize nucleotide diversity have ranged widely in scope, but differences in reference genome version and software algorithms used in these efforts inhibit comparison, and these data are generally not available in an easy-to-use visualization platform for quick access and analysis. To address these issues, The Maize Genetics and Genomics Database has collaborated with maize researchers to offer variant data from a diverse set of 1,498 inbred lines, traditional varieties, and teosintes through a standardized variant-calling pipeline against version 5 of the B73 reference genome. The output was filtered for mapping quality, completeness, and linkage disequilibrium, and annotated based on variant effects relative to the B73 RefGen_v5 gene annotations. MaizeGDB has also updated a web tool, SNPversity 2.0, to filter, visualize, and download genotype sets based on genomic locations and accessions of interest, and added external datasets to demonstrate SNPversity 2.0's broad usage. MaizeGDB plans to host annual updates of these resources as additional resequencing data become available, with plans to expand to all publicly available sequence data.

Genome-wide association mapping and genomic prediction analyses reveal the genetic architecture of grain yield and agronomic traits under drought and optimum conditions in maize

BACKGROUND

Drought is a major abiotic stress in sub-Saharan Africa, impacting maize growth and development leading to severe yield loss. Drought tolerance is a complex trait regulated by multiple genes, making direct grain yield selection ineffective. To dissect the genetic architecture of grain yield and flowering traits under drought stress, a genome-wide association study (GWAS) was conducted on a panel of 236 maize lines testcrossed and evaluated under managed drought and optimal growing conditions in multiple environments using seven multi-locus GWAS models (mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB, pKWmEB, ISIS EM-BLASSO, and FARMCPU) from mrMLM and GAPIT R packages. Genomic prediction with RR-BLUP model was applied on BLUEs across locations under optimum and drought conditions.

RESULTS

A total of 172 stable and reliable quantitative trait nucleotides (QTNs) were identified, of which 77 are associated with GY, AD, SD, ASI, PH, EH, EPO and EPP under drought and 95 are linked to GY, AD, SD, ASI, PH, EH, EPO and EPP under optimal conditions. Among these QTNs, 17 QTNs explained over 10% of the phenotypic variation (R2 ≥ 10%). Furthermore, 43 candidate genes were discovered and annotated. Two major candidate genes, Zm00001eb041070 closely associated with grain yield near peak QTN, qGY_DS1.1 (S1_216149215) and Zm00001eb364110 closely related to anthesis-silking interval near peak QTN, qASI_DS8.2 (S8_167256316) were identified, encoding AP2-EREBP transcription factor 60 and TCP-transcription factor 20, respectively under drought stress. Haplo-pheno analysis identified superior haplotypes for qGY_DS1.1 (S1_216149215) associated with the higher grain yield under drought stress. Genomic prediction revealed moderate to high prediction accuracies under optimum and drought conditions.

CONCLUSION

The lines carrying superior haplotypes can be used as potential donors in improving grain yield under drought stress. Integration of genomic selection with GWAS results leads not only to an increase in the prediction accuracy but also to validate the function of the identified candidate genes as well increase in the accumulation of favorable alleles with minor and major effects in elite breeding lines. This study provides valuable insight into the genetic architecture of grain yield and secondary traits under drought stress.

Developing pangenomes for large and complex plant genomes and their representation formats

BACKGROUND

The development of pangenomes has revolutionized genomic studies by capturing the complete genetic diversity within a species. Pangenome assembly integrates data from multiple individuals to construct a comprehensive genomic landscape, revealing both core and accessory genomic elements. This approach enables the identification of novel genes, structural variations, and gene presence-absence variations, providing insights into species evolution, adaptation, and trait variation. Representing pangenomes requires innovative visualization formats that effectively convey the complex genomic structures and variations.

AIM

This review delves into contemporary methodologies and recent advancements in constructing pangenomes, particularly in plant genomes. It examines the structure of pangenome representation, including format comparison, conversion, visualization techniques, and their implications for enhancing crop improvement strategies.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Earlier comparative studies have illuminated novel gene sequences, copy number variations, and presence-absence variations across diverse crop species. The concept of a pan-genome, which captures multiple genetic variations from a broad spectrum of genotypes, offers a holistic perspective of a species' genetic makeup. However, constructing a pan-genome for plants with larger genomes poses challenges, including managing vast genome sequence data and comprehending the genetic variations within the germplasm. To address these challenges, researchers have explored cost-effective alternatives to encapsulate species diversity in a single assembly known as a pangenome. This involves reducing the volume of genome sequences while focusing on genetic variations. With the growing prominence of the pan-genome concept in plant genomics, several software tools have emerged to facilitate pangenome construction. This review sheds light on developing and utilizing software tools tailored for constructing pan-genomes in plants. It also discusses representation formats suitable for downstream analyses, offering valuable insights into the genetic landscape and evolutionary dynamics of plant species. In summary, this review underscores the significance of pan-genome construction and representation formats in resolving the genetic architecture of plants, particularly those with complex genomes. It provides a comprehensive overview of recent advancements, aiding in exploring and understanding plant genetic diversity.

Transcriptome profiling of maize transcription factor mutants to probe gene regulatory network predictions

Transcription factors play important roles in regulation of gene expression and phenotype. A variety of approaches have been utilized to develop gene regulatory networks to predict the regulatory targets for each transcription factor, such as yeast-1-hybrid screens and gene co-expression network analysis. Here we identified potential transcription factor targets and used a reverse genetics approach to test the predictions of several gene regulatory networks in maize. Loss-of-function mutant alleles were isolated for 22 maize transcription factors. These mutants did not exhibit obvious morphological phenotypes. However, transcriptomic profiling identified differentially expressed genes in each of the mutant genotypes, and targeted metabolic profiling indicated variable phenolic accumulation in some mutants. An analysis of expression levels for predicted target genes based on yeast-1-hybrid screens identified a small subset of predicted targets that exhibit altered expression levels. The analysis of predicted targets from gene co-expression network-based methods found significant enrichments for prediction sets of some transcription factors, but most predicted targets did not exhibit altered expression. This could result from false-positive gene co-expression network predictions, a transcription factor with a secondary regulatory role resulting in minor effects on gene regulation, or redundant gene regulation by other transcription factors. Collectively, these findings suggest that loss-of-function for single uncharacterized transcription factors might have limited phenotypic impacts but can reveal subsets of gene regulatory network predicted targets with altered expression.

Identification of cell-type specificity, trans- and cis-acting functions of plant lincRNAs from single-cell transcriptomes

Long noncoding RNAs, including intergenic lncRNAs (lincRNAs), play a key role in various biological processes throughout the plant life cycle, and the advent of single-cell RNA sequencing (scRNA-seq) technology has opened up a valuable avenue for scrutinizing the intricate roles of lincRNAs in cellular processes. Here, we identified a new batch of lincRNAs using scRNA-seq data from diverse tissues of plants (rice, Arabidopsis, tomato, and maize). Based on well-annotated single-cell transcriptome atlases, plant lincRNAs were found to possess the same level of cell-type specificity as mRNAs and to be involved in the differentiation of certain cell types based on pseudo-time analysis. Many lincRNAs were predicted to play a hub role in the cell-type-specific co-expression networks of lincRNAs and mRNAs, suggesting their trans-acting abilities. Besides, plant lincRNAs were revealed to have potential cis-acting properties based on their genomic distances and expression correlations with the neighboring mRNAs. Furthermore, an online platform, PscLncRNA (http://ibi.zju.edu.cn/psclncrna/), was constructed for searching and visualizing all identified plant lincRNAs with annotated potential functions. Our work provides new insights into plant lincRNAs at single-cell resolution and an important resource for understanding and further investigation of plant lincRNAs.

Embracing native diversity to enhance the maximum quantum efficiency of photosystem II in maize

The sustainability of maize cultivation would benefit tremendously from early sowing, but is hampered by low temperatures during early development in temperate climates. We show that allelic variation within the gene encoding subunit M of the NADH-dehydrogenase-like (NDH) complex (ndhm1) in a European maize landrace affects several quantitative traits that are relevant during early development in cold climates through NDH-mediated cyclic electron transport around photosystem I, a process crucial for photosynthesis and photoprotection. Beginning with a genome-wide association study for maximum potential quantum yield of photosystem II in dark-adapted leaves (Fv/Fm), we capitalized on the large phenotypic effects of a hAT transposon insertion in ndhm1 on multiple quantitative traits (early plant height [EPH], Fv/Fm, chlorophyll content, and cold tolerance) caused by the reduced protein levels of NDHM and associated NDH components. Analysis of the ndhm1 native allelic series revealed a rare allele of ndhm1 that is associated with small albeit significant improvements of Fv/Fm, photosystem II efficiency in light-adapted leaves (ΦPSII), and EPH compared with common alleles. Our work showcases the extraction of favorable alleles from locally adapted landraces, offering an efficient strategy for broadening the genetic variation of elite germplasm by breeding or genome editing.

Ribosome binding of phasiRNA precursors accelerates the 24-nt phasiRNA burst in meiotic maize anthers

Reproductive phasiRNAs (phased, secondary, small interfering RNAs), produced from numerous PHAS loci, are essential for plant anther development. PHAS transcripts are enriched on endoplasmic reticulum-bound ribosomes in maize (Zea mays), but the impact of ribosome binding on phasiRNA biogenesis remains elusive. Through ribosome profiling of maize anthers at 10 developmental stages, we demonstrated that 24-PHAS transcripts are bound by ribosomes, with patterns corresponding to the timing and abundance of 24-PHAS transcripts. Ribosome binding to 24-PHAS transcripts is conserved among different maize inbred lines, with ribosomes enriched upstream of the miR2275 target sites. We detected short open reading frames (sORFs) in the ribosome-binding regions of some 24-PHAS transcripts and observed a 3-nt periodicity in most sORFs, but mass spectrometry failed to detect peptides corresponding to the sORFs. Deletion of the entire ribosome-binding region of 24PHAS_NO296 locus eliminated ribosome binding and decreased 24-nt phasiRNA production, without affecting 24PHAS_NO296 transcript levels. In contrast, disrupting only the sORFs in 24PHAS_NO296 did not substantially affect the generation of 24-nt phasiRNAs. A newly formed sORF in these mutants may have re-directed ribosome binding to its transcripts. Overall, these findings demonstrate that sORFs facilitate ribosome binding to 24-PHAS transcripts, thereby promoting phasiRNA biogenesis in meiotic anthers.

Dynamic Chromatin Accessibility and Gene Expression Regulation During Maize Leaf Development

BACKGROUND/OBJECTIVES

Chromatin accessibility is closely associated with transcriptional regulation during maize (Zea mays) leaf development. However, its precise role in controlling gene expression at different developmental stages remains poorly understood. This study aimed to investigate the dynamics of chromatin accessibility and its influence on genome-wide gene expression during the BBCH_11, BBCH_13, and BBCH_17 stages of maize leaf development.

METHODS

Maize leaves were collected at the BBCH_11, BBCH_13, and BBCH_17 developmental stages, and chromatin accessibility was assessed using ATAC-seq. RNA-seq was performed to profile gene expression. Integrated analysis of ATAC-seq and RNA-seq data was conducted to elucidate the relationship between chromatin accessibility and transcriptional regulation.

RESULTS

A total of 46,808, 38,242, and 41,084 accessible chromatin regions (ACRs) were identified at the BBCH_11, BBCH_13, and BBCH_17 stages, respectively, with 23.4%, 12.2%, and 21.9% of these regions located near transcription start sites (TSSs). Integrated analyses revealed that both the number and intensity of ACRs significantly influence gene expression levels. Motif analysis identified key transcription factors associated with leaf development and potential transcriptional repressors among genes, showing divergent regulation patterns in ATAC-seq and RNA-seq datasets.

CONCLUSIONS

These findings demonstrate that chromatin accessibility plays a crucial role in regulating the spatial and temporal expression of key genes during maize leaf development by modulating transcription factor binding. This study provides novel insights into the regulatory mechanisms underlying maize leaf development, contributing to a deeper understanding of chromatin-mediated gene expression.

Evolution of Plant Genome Size and Composition

The rapid development of sequencing technology has led to an explosion of plant genome data, opening up more opportunities for research in the field of comparative evolutionary analysis of plant genomes. In this review, we focus on changes in plant genome size and composition, examining the effects of polyploidy, whole-genome duplication, and alternations in transposable elements on plant genome architecture and evolution, respectively. In addition, to address gaps in the available information, we also collected and analyzed 234 representative plant genome data as a supplement. We aim to provide a comprehensive, up-to-date summary of information on plant genome architecture and evolution in this review.

A metabolic roadmap of waxy corn flavor

As well as being a popular vegetable crop worldwide, waxy corn represents an important amylopectin source, but little is known about its breeding history and flavor characteristics. In this study, through comparative-omic analyses between 318 diverse waxy corn and 507 representative field corn inbred lines we revealed that many metabolic pathways and genes exhibited selection characteristics during the breeding history of waxy corn, contributing to the divergence between waxy and field corn. We showed that waxy corn is not only altered in its glutinous property but also its sweetness, aroma, and palatability are all significantly affected. A substantial proportion (43%) of flavor-related metabolites have pleiotropic effects, affecting both flavor and yield characteristics, and 27% of these metabolites are related to antagonistic outcomes on yield and flavor. Furthermore, through multiple concrete examples, we demonstrated how yield and quality are coordinately or antagonistically regulated at the genetic level. In particular, some sweet molecules, such as DIMBOA and raffinose, which do not participate in the starch biosynthesis pathway, were identified as potential targets for breeding a new type of "sweet-waxy" corn. Taken together, our findings shed light on the historical selection of waxy corn and demonstrate the genetic and metabolic basis of waxy corn flavor, collectively providing valuable resources and knowledge for future crop breeding for improved nutritional quality.

ZmICE1a regulates the defence-storage trade-off in maize endosperm

The endosperm of cereal grains feeds the entire world as a major food supply; however, little is known about its defence response during endosperm development. The Inducer of CBF Expression 1 (ICE1) is a well-known regulator of cold tolerance in plants. ICE1 has a monocot-specific homologue that is preferentially expressed in cereal endosperms but with an unclear regulatory function. Here we characterized the function of monocot-specific ZmICE1a, which is expressed in the entire endosperm, with a predominant expression in its peripheral regions, including the aleurone layer, subaleurone layer and basal endosperm transfer layer in maize (Zea mays). Loss of function of ZmICE1a reduced starch content and kernel weight. RNA sequencing and CUT&Tag-seq analyses revealed that ZmICE1a positively regulates genes in starch synthesis while negatively regulating genes in aleurone layer-specific defence and the synthesis of indole-3-acetic acid and jasmonic acid (JA). Exogenous indole-3-acetic acid and JA both induce the expression of numerous defence genes, which show distinct spatial-specific expression in the basal endosperm transfer layer and subaleurone layer, respectively. Moreover, we dissected a JA-ZmJAZ9-ZmICE1a-MPI signalling axis involved in JA-mediated defence regulation. Overall, our study revealed ZmICE1a as a key regulator of endosperm defence response and a coordinator of the defence-storage trade-off in endosperm development.

The C2H2-type zinc finger transcription factor ZmDi19-7 regulates plant height and organ size by promoting cell size in maize

The drought-induced protein 19 (Di19) gene family encodes a Cys2/His2 zinc-finger protein implicated in responses to diverse plant stressors. To date, potential roles of these proteins as transcription factors remain largely elusive in maize. Here, we show that ZmDi19-7 gene exerts pivotal functions in regulation of plant height and organ growth by modulating the cell size in maize. ZmDi19-7 physically interacts with ubiquitin receptor protein ZmDAR1b, which is indispensable in ubiquitination of ZmDi19-7 and affects its protein stability. Further genetic analysis demonstrated that ZmDAR1b act in a common pathway with ZmDi19-7 to regulate cell size in maize. ZmDi19-7, severing as a transcriptional factor, is significantly enriched in conserved DiBS element in the promoter region of ZmHSP22, ZmHSP18c, ZmSAUR25, ZmSAUR55, ZmSAUR7 and ZmXTH23 and orchestrates the expression of these genes involving in auxin-mediated cell expansion and protein processing in the endoplasmic reticulum. Thus, our findings demonstrate that ZmDi19-7 is an important newfound component of the ubiquitin-proteasome pathway in regulation of plant height and organ size in maize. These discoveries highlight potential targets for the genetic improvement of maize in the future.

Expanding the BonnMu sequence-indexed repository of transposon induced maize (Zea mays L.) mutations in dent and flint germplasm

The BonnMu resource is a transposon tagged mutant collection designed for functional genomics studies in maize. To expand this resource, we crossed an active Mutator (Mu) stock with dent (B73, Co125) and flint (DK105, EP1, and F7) germplasm, resulting in the generation of 8064 mutagenized BonnMu F2-families. Sequencing of these Mu-tagged families revealed 425 924 presumptive heritable Mu insertions affecting 36 612 (83%) of the 44 303 high-confidence gene models of maize (B73v5). On average, we observed 12 Mu insertions per gene (425 924 total insertions/36 612 affected genes) and 53 insertions per BonnMu F2-family (425 924 total insertions/8064 families). Mu insertions and photos of seedling phenotypes from segregating BonnMu F2-families can be accessed through the Maize Genetics and Genomics Database (MaizeGDB). Downstream examination via the automated Mutant-seq Workflow Utility (MuWU) identified 94% of the presumptive germinal insertion sites in genic regions and only a small fraction of 6% inserting in non-coding intergenic sequences of the genome. Consistently, Mu insertions aligned with gene-dense chromosomal arms. In total, 42% of all BonnMu insertions were located in the 5' untranslated region of genes, corresponding to accessible chromatin. Furthermore, for 38% of the insertions (163 843 of 425 924 total insertions) Mu1, Mu8 and MuDR were confirmed to be the causal Mu elements. Our publicly accessible European BonnMu resource has archived insertions covering two major germplasm groups, thus facilitating both forward and reverse genetics studies.

Transcriptomic Characterization of Genes Harboring Markers Linked to Maize Yield

BACKGROUND

It is currently believed that breeding priorities, including maize breeding, should focus on introducing varieties with greater utility value, specifically higher yields, into production. Global modern maize breeding relies on various molecular genetics techniques. Using the above mentioned technologies, we can identify regions of the genome that are associated with various phenotypic traits, including yield, which is of fundamental importance for understanding and manipulating these regions.

OBJECTIVES

The aim of the study was to analyze the expression of candidate genes associated with maize yield. To better understand the function of the analyzed genes in increasing maize yield, their expression in different organs and tissues was also assessed using publicly available transcriptome data.

METHODS

RT-qPCR analyses were performed using iTaq Universal SYBR Green Supermix (Bio-Rad, Hercules, CA, USA) and CFX96 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). Each of the performed RT-qPCR experiments consisted of three biological replicates and three technical replicates, the results of which were averaged.

RESULTS

The research results allowed us to select three out of six candidate genes (cinnamoyl-CoA reductase 1-CCR1, aspartate aminotransferase-AAT and sucrose transporter 1-SUT1), which can significantly affect grain yield in maize. Not only our studies but also literature reports clearly indicate the participation of CCR1, AAT and SUT1 in the formation of yield. Identified molecular markers located within these genes can be used in breeding programs to select high yielding maize genotypes.

Combined Genome-Wide Association Study and Linkage Analysis for Mining Candidate Genes for the Kernel Row Number in Maize (Zea mays L.)

Kernel row number (KRN) is one of the key traits that significantly affect maize yield and productivity. Therefore, investigating the candidate genes and their functions in regulating KRN provides a theoretical basis and practical direction for genetic improvement in maize breeding, which is vital for increasing maize yield and understanding domestication. In this study, three recombinant inbred line (RIL) populations were developed using the parental lines AN20, YML1218, CM395, and Ye107, resulting in a multiparent population comprising a total of 490 F9 RILs. Phenotypic evaluation of the RILs for KRN was performed in three distinct environments. The heritability estimates of the RILs ranged from 81.40% to 84.16%. Genotyping-by-sequencing (GBS) of RILs identified 569,529 high-quality single nucleotide polymorphisms (SNPs). Combined genome-wide association study (GWAS) and linkage analyses revealed 120 SNPs and 22 quantitative trait loci (QTLs) which were significantly associated with KRN in maize. Furthermore, two novel candidate genes, Zm00001d042733 and Zm00001d042735, regulating KRN in maize were identified, which were located in close proximity to the significant SNP3-178,487,003 and overlapping the interval of QTL qKRN3-1. Zm00001d042733 encodes ubiquitin carboxyl-terminal hydrolase and Zm00001d042735 encodes the Arabidopsis Tóxicos en Levadura family of proteins. This study identified novel candidate loci and established a theoretical foundation for further functional validation of candidate genes. These findings deepen our comprehension of the genetic mechanisms that underpin KRN and offer potential applications of KRN-related strategies in developing maize varieties with higher yield.

Transcriptome profiling of eight Zea mays lines identifies genes responsible for the resistance to Fusarium verticillioides

BACKGROUND

The cultivation of maize (Zea mays L.), one of the most important crops worldwide for food, feed, biofuels, and industrial applications, faces significant constraints due to Fusarium verticillioides, a fungus responsible for severe diseases including seedling blights, stalk rot, and ear rot. Its impact is worsened by the fact that chemical and agronomic measures used to control the infection are often inefficient. Hence, genetic resistance is considered the most reliable resource to reduce the damage. This study aims to elucidate the genetic basis of F. verticillioides resistance in maize.

RESULTS

Young seedlings of eight divergent maize lines, founders of the MAGIC population, were artificially inoculated with a F. verticillioides strain. Phenotypic analysis and transcriptome sequencing of both control and treated samples identified several hundred differentially expressed genes enriched in metabolic processes associated with terpene synthesis. A WGCNA further refined the pool of genes with potential implications in disease response and found a limited set of hub genes, encoding bZIP and MYB transcription factors, or involved in carbohydrate metabolism, solute transport processes, calcium signaling, and lipid pathways. Finally, additional gene resources were provided by combining transcriptomic data with previous QTL mapping, thereby shedding light on the molecular mechanisms in the maize-F. verticillioides interaction.

CONCLUSIONS

The transcriptome profiling of eight divergent MAGIC maize founder lines with contrasting levels of Fusarium verticillioides resistance combined with phenotypic analysis, clarifies the molecular mechanisms underlying the maize-F. verticillioides interaction.

The establishment of the anther somatic niche with single-cell sequencing

The anther is the developmental housing of pollen and therefore the male gametes of flowering plants. The meiotic cells from which pollen are derived must differentiate de novo from somatic anther cells and synchronously develop with the rest of the anther. Anthropogenic control over another development has become crucial for global agriculture so as to maintain inbred lines and generate hybrid seeds of many crops. Understanding the genes that underlie the proper differentiation, developmental landmarks, and functions of each anther cell type is thus fundamental to both basic and applied plant sciences. We investigated the development of the somatic niche of the maize (Zea mays) anther using single-cell RNA-seq (scRNA-seq). Extensive background knowledge on the birth then pace and pattern of cell division of the maize anther cell types and published examples of cell-type gene expression from in situ hybridization allowed us to identify the primary cell types within the anther lobe, as well as the connective cells between the four lobes. We established the developmental trajectories of somatic cell types from pre-meiosis to post-meiosis, identified putative marker genes for the somatic cell types that previously lacked any known specific functions, and addressed the possibility that tapetal cells sequentially differentiate. This comprehensive scRNA-seq dataset of the somatic niche of the maize anther will serve as a baseline for future analyses investigating male-sterile genotypes and the impact of environmental conditions on male fertility in flowering plants.

QTL mapping and genome-wide association analysis reveal genetic loci and candidate gene for resistance to gray leaf spot in tropical and subtropical maize germplasm

KEY MESSAGE

Using QTL mapping and GWAS, two candidate genes (Zm00001d051039 and Zm00001d051147) were consistently identified across the three different environments and BLUP values. GWAS analysis identified the candidate gene, Zm00001d044845. These genes were subsequently validated to exhibit a significant association with maize gray leaf spot (GLS) resistance. Gray leaf spot (GLS) is a major foliar disease of maize (Zea mays L.) that causes significant yield losses worldwide. Understanding the genetic mechanisms underlying gray leaf spot resistance is crucial for breeding high-yielding and disease-resistant varieties. In this study, eight tropical and subtropical germplasms were crossed with the temperate germplasm Ye107 to develop a nested association mapping (NAM) population comprising 1,653 F2:8 RILs, consisting of eight recombinant inbred line (RIL) subpopulations, using the single-seed descent method. The NAM population was evaluated for GLS resistance in three different environments, and genotyping by sequencing of the NAM population generated 593,719 high-quality single-nucleotide polymorphisms (SNPs). Linkage analysis and genome-wide association studies (GWASs) were conducted to identify candidate genes regulating GLS resistance in maize. Both analyses identified 25 QTLs and 149 SNPs that were significantly associated with GLS resistance. Candidate genes were screened 20 Kb upstream and downstream of the significant SNPs, and three novel candidate genes (Zm00001d051039, Zm00001d051147, and Zm00001d044845) were identified. Zm00001d051039 and Zm00001d051147 were located on chromosome 4 and co-localized in both linkage (qGLS4-1 and qGLS4-2) and GWAS analyses. SNP-138,153,206 was located 0.499 kb downstream of the candidate gene Zm00001d051039, which encodes the protein IN2-1 homolog B, a homolog of glutathione S-transferase (GST). GSTs and protein IN2-1 homolog B scavenge reactive oxygen species under various stress conditions, and GSTs are believed to protect plants from a wide range of biotic and abiotic stresses by detoxifying reactive electrophilic compounds. Zm00001d051147 encodes a probable beta-1,4-xylosyltransferase involved in the biosynthesis of xylan in the cell wall, enhancing resistance. SNP-145,813,215 was located 2.69 kb downstream of the candidate gene. SNP-5,043,412 was consistently identified in three different environments and BLUP values and was located 8.788 kb downstream of the candidate gene Zm00001d044845 on chromosome 9. Zm00001d044845 encodes the U-box domain-containing protein 4 (PUB4), which is involved in regulating plant immunity. qRT-PCR analysis showed that the relative expression levels of the three candidate genes were significantly upregulated in the leaves of the TML139 (resistant) parent, indicating that these three candidate genes could be associated with resistance to GLS. The findings of this study are significant for marker-assisted breeding aimed at enhancing resistance to GLS in maize and lay the foundation for further elucidation of the genetic mechanisms underlying resistance to gray leaf spot in maize and breeding of new disease-resistant varieties.

Tip of the iceberg? Three novel TOPLESS-interacting effectors of the gall-inducing fungus Ustilago maydis

Ustilago maydis is a biotrophic pathogen causing smut disease in maize. It secretes a cocktail of effector proteins, which target different host proteins during its biotrophic stages in the host plant. One such class of proteins we identified previously is TOPLESS (TPL) and TOPLESS-RELATED (TPR) transcriptional corepressors. Here, we screened 297 U. maydis effector candidates for their ability to interact with maize TPL protein RAMOSA 1 ENHANCER LOCUS 2 LIKE 2 (RELK2) and their ability to induce auxin signaling and thereby identified three novel TPL-interacting protein effectors (Tip6, Tip7, and Tip8). Structural modeling and mutational analysis allowed the identification of TPL-interaction motifs of Tip6 and Tip7. In planta interaction between Tip6 and Tip7 with RELK2 occurs mainly in nuclear compartments, whereas Tip8 colocalizes with RELK2 in a compartment outside the nucleus. Overexpression of Tip8 in nonhost plants leads to cell death, indicating recognition of the effector or its activity. By performing infection assays with single and multideletion mutants of U. maydis, we demonstrate a positive role of Tip6 and Tip7 in U. maydis virulence. Transcriptional profiling of maize leaves infected with Tip effector mutants in comparison with SG200 strain suggests Tip effector activities are not merely redundant.

C4 cereal and biofuel crop microbiomes

Microbiomes provide multiple life-support functions for plants, including nutrient acquisition and tolerance to abiotic and biotic stresses. Considering the importance of C4 cereal and biofuel crops for food security under climate change conditions, more attention has been given recently to C4 plant microbiome assembly and functions. Here, we review the current status of C4 cereal and biofuel crop microbiome research with a focus on beneficial microbial traits for crop growth and health. We highlight the importance of environmental factors and plant genetics in C4 crop microbiome assembly and pinpoint current knowledge gaps. Finally, we discuss the potential of foxtail millet as a C4 model species and outline future perspectives of C4 plant microbiome research.

GPCR-like Protein ZmCOLD1 Regulate Plant Height in an ABA Manner

G protein-coupled receptors (GPCRs) are sensors for the G protein complex to sense changes in environmental factors and molecular switches for G protein complex signal transduction. In this study, the homologous gene of GPCR-like proteins was identified from maize and named as ZmCOLD1. Subcellular analysis showed that the ZmCOLD1 protein is localized to the cell membrane and endoplasmic reticulum. A CRISPR/Cas9 knock-out line of ZmCOLD1 was further created and its plant height was significantly lower than the wild-type maize at both the seedling and adult stages. Histological analysis showed that the increased cell number but significantly smaller cell size may result in dwarfing of zmcold1, indicating that the ZmCOLD1 gene could regulate plant height development by affecting the cell division process. Additionally, ZmCOLD1 was verified to interact with the maize Gα subunit, ZmCT2, though the central hydrophilic loop domain by in vivo and in vitro methods. Abscisic acid (ABA) sensitivity analysis by seed germination assays exhibited that zmcold1 were hypersensitive to ABA, indicating its important roles in ABA signaling. Finally, transcriptome analysis was performed to investigate the transcriptional change in zmcold1 mutant. Overall, ZmCOLD1 functions as a GPCR-like protein and an important regulator to plant height.

Natural Diversity of Heat-Induced Transcription of Retrotransposons in Arabidopsis thaliana

Transposable elements (TEs) are major components of plant genomes, profoundly impacting the fitness of their hosts. However, technical bottlenecks have long hindered our mechanistic understanding of TEs. Using RNA-Seq and long-read sequencing with Oxford Nanopore Technologies' (ONT) direct cDNA sequencing, we analyzed the heat-induced transcription of TEs in three natural accessions of Arabidopsis thaliana (Cvi-0, Col-0, and Ler-1). In addition to the well-studied ONSEN retrotransposon family, we confirmed Copia-35 as a second heat-responsive retrotransposon family with particularly high activity in the relict accession Cvi-0. Our analysis revealed distinct expression patterns of individual TE copies and suggest different mechanisms regulating the GAG protein production in the ONSEN versus Copia-35 families. In addition, analogously to ONSEN, Copia-35 activation led to the upregulation of flanking genes such as APUM9 and potentially to the quantitative modulation of flowering time. ONT data allowed us to test the extent to which read-through formation is important in the regulation of adjacent genes. Unexpectedly, our results indicate that for both families, the upregulation of flanking genes is not predominantly directly initiated by transcription from their 3' long terminal repeats. These findings highlight the intraspecific expressional diversity linked to retrotransposon activation under stress.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Cuticle development and the underlying transcriptome-metabolome associations during early seedling establishment

The plant cuticle is a complex extracellular lipid barrier that has multiple protective functions. This study investigated cuticle deposition by integrating metabolomics and transcriptomics data gathered from six different maize seedling organs of four genotypes, the inbred lines B73 and Mo17, and their reciprocal hybrids. These datasets captured the developmental transition of the seedling from heterotrophic skotomorphogenic growth to autotrophic photomorphogenic growth, a transition that is highly vulnerable to environmental stresses. Statistical interrogation of these data revealed that the predominant determinant of cuticle composition is seedling organ type, whereas the seedling genotype has a smaller effect on this phenotype. Gene-to-metabolite associations assessed by integrated statistical analyses identified three gene networks associated with the deposition of different elements of the cuticle: cuticular waxes; monomers of lipidized cell wall biopolymers, including cutin and suberin; and both of these elements. These gene networks reveal three metabolic programs that appear to support cuticle deposition, including processes of chloroplast biogenesis, lipid metabolism, and molecular regulation (e.g. transcription factors, post-translational regulators, and phytohormones). This study demonstrates the wider physiological metabolic context that can determine cuticle deposition and lays the groundwork for new targets for modulating the properties of this protective barrier.

Genomic prediction of cereal crop architectural traits using models informed by gene regulatory circuitries from maize

Plant architecture is a major determinant of planting density, which enhances productivity potential for crops per unit area. Genomic prediction is well positioned to expedite genetic gain of plant architectural traits since they are typically highly heritable. Additionally, the adaptation of genomic prediction models to query predictive abilities of markers tagging certain genomic regions could shed light on the genetic architecture of these traits. Here, we leveraged transcriptional networks from a prior study that contextually described developmental progression during tassel and leaf organogenesis in maize (Zea mays) to inform genomic prediction models for architectural traits. Since these developmental processes underlie tassel branching and leaf angle, 2 important agronomic architectural traits, we tested whether genes prioritized from these networks quantitatively contribute to the genetic architecture of these traits. We used genomic prediction models to evaluate the ability of markers in the vicinity of prioritized network genes to predict breeding values of tassel branching and leaf angle traits for 2 diversity panels in maize and diversity panels from sorghum (Sorghum bicolor) and rice (Oryza sativa). Predictive abilities of markers near these prioritized network genes were similar to those using whole-genome marker sets. Notably, markers near highly connected transcription factors from core network motifs in maize yielded predictive abilities that were significantly greater than expected by chance in not only maize but also closely related sorghum. We expect that these highly connected regulators are key drivers of architectural variation that are conserved across closely related cereal crop species.

Review of the technology used for structural characterization of the GMO genome using NGS data

The molecular characterization of genetically modified organisms (GMOs) is essential for ensuring safety and gaining regulatory approval for commercialization. According to CODEX standards, this characterization involves evaluating the presence of introduced genes, insertion sites, copy number, and nucleotide sequence structure. Advances in technology have led to the increased use of next-generation sequencing (NGS) over traditional methods such as Southern blotting. While both methods provide high reproducibility and accuracy, Southern blotting is labor-intensive and time-consuming due to the need for repetitive probe design and analyses for each target, resulting in low throughput. Conversely, NGS facilitates rapid and comprehensive analysis by mapping whole-genome sequencing (WGS) data to plasmid sequences, accurately identifying T-DNA insertion sites and flanking regions. This advantage allows for efficient detection of T-DNA presence, copy number, and unintended gene insertions without additional probe work. This paper reviews the current status of GMO genome characterization using NGS and proposes more efficient strategies for this purpose.

Tubulin participates in establishing protoxylem vessel reinforcement patterns and hydraulic conductivity in maize

Water transportation to developing tissues relies on the structure and function of plant xylem cells. Plant microtubules govern the direction of cellulose microfibrils and guide secondary cell wall formation and morphogenesis. However, the relevance of microtubule-determined xylem wall thickening patterns in plant hydraulic conductivity remains unclear. In the present study, we identified a maize (Zea mays) semi-dominant mutant, designated drought-overly-sensitive1 (ZmDos1), the upper leaves of which wilted even when exposed to well-watered conditions during growth; the wilting phenotype was aggravated by increased temperatures and decreased humidity. Protoxylem vessels in the stem and leaves of the mutant showed altered thickening patterns of the secondary cell wall (from annular to spiral), decreased inner diameters, and limited water transport efficiency. The causal mutation for this phenotype was found to be a G-to-A mutation in the maize gene α-tubulin4, resulting in a single amino acid substitution at position 196 (E196K). Ectopic expression of the mutant α-tubulin4 in Arabidopsis (Arabidopsis thaliana) changed the orientation of microtubule arrays, suggesting a determinant role of this gene in microtubule assembly and secondary cell wall thickening. Our findings suggest that the spiral wall thickenings triggered by the α-tubulin mutation are stretched during organ elongation, causing a smaller inner diameter of the protoxylem vessels and affecting water transport in maize. This study underscores the importance of tubulin-mediated protoxylem wall thickening in regulating plant hydraulics, improves our understanding of the relationships between protoxylem structural features and functions, and offers candidate genes for the genetic enhancement of maize.

Trait association and prediction through integrative k-mer analysis

Genome-wide association study (GWAS) with single nucleotide polymorphisms (SNPs) has been widely used to explore genetic controls of phenotypic traits. Alternatively, GWAS can use counts of substrings of length k from longer sequencing reads, k-mers, as genotyping data. Using maize cob and kernel color traits, we demonstrated that k-mer GWAS can effectively identify associated k-mers. Co-expression analysis of kernel color k-mers and genes directly found k-mers from known causal genes. Analyzing complex traits of kernel oil and leaf angle resulted in k-mers from both known and candidate genes. A gene encoding a MADS transcription factor was functionally validated by showing that ectopic expression of the gene led to less upright leaves. Evolution analysis revealed most k-mers positively correlated with kernel oil were strongly selected against in maize populations, while most k-mers for upright leaf angle were positively selected. In addition, genomic prediction of kernel oil, leaf angle, and flowering time using k-mer data resulted in a similarly high prediction accuracy to the standard SNP-based method. Collectively, we showed k-mer GWAS is a powerful approach for identifying trait-associated genetic elements. Further, our results demonstrated the bridging role of k-mers for data integration and functional gene discovery.

Potent pollen gene regulation by DNA glycosylases in maize

Although DNA methylation primarily represses TEs, it also represses select genes that are methylated in plant body tissues but demethylated by DNA glycosylases (DNGs) in endosperm or pollen. Either one of two DNGs, MATERNAL DEREPRESSION OF R1 (MDR1) or DNG102, is essential for pollen viability in maize. Using single-pollen mRNA sequencing on pollen-segregating mutations in both genes, we identify 58 candidate DNG target genes that account for 11.1% of the wild-type transcriptome but are silent or barely detectable in other tissues. They are unusual in their tendency to lack introns but even more so in their TE-like methylation (teM) in coding DNA. The majority have predicted functions in cell wall modification, and they likely support the rapid tip growth characteristic of pollen tubes. These results suggest a critical role for DNA methylation and demethylation in regulating maize genes with the potential for extremely high expression in pollen but constitutive silencing elsewhere.

Cortical parenchyma wall width regulates root metabolic cost and maize performance under suboptimal water availability

We describe how increased root cortical parenchyma wall width (CPW) can improve tolerance to drought stress in maize by reducing the metabolic costs of soil exploration. Significant variation (1.0-5.0 µm) for CPW was observed in maize germplasm. The functional-structural model RootSlice predicts that increasing CPW from 2 µm to 4 µm is associated with a ~15% reduction in root cortical cytoplasmic volume, respiration rate, and nitrogen content. Analysis of genotypes with contrasting CPW grown with and without water stress in the field confirms that increased CPW is correlated with an ~32-42% decrease in root respiration. Under water stress in the field, increased CPW is correlated with 125% increased stomatal conductance, 325% increased leaf CO2 assimilation rate, 73-78% increased shoot biomass, and 92-108% increased yield. CPW was correlated with leaf mesophyll midrib parenchyma wall width, indicating pleiotropy. Genome-wide association study analysis identified candidate genes underlying CPW. OpenSimRoot modeling predicts that a reduction in root respiration due to increased CPW would also benefit maize growth under suboptimal nitrogen, which requires empirical testing. We propose CPW as a new phene that has utility under edaphic stress meriting further investigation.

A dynamic regulome of shoot-apical-meristem-related homeobox transcription factors modulates plant architecture in maize

BACKGROUND

The shoot apical meristem (SAM), from which all above-ground tissues of plants are derived, is critical to plant morphology and development. In maize (Zea mays), loss-of-function mutant studies have identified several SAM-related genes, most encoding homeobox transcription factors (TFs), located upstream of hierarchical networks of hundreds of genes.

RESULTS

Here, we collect 46 transcriptome and 16 translatome datasets across 62 different tissues or stages from the maize inbred line B73. We construct a dynamic regulome for 27 members of three SAM-related homeobox subfamilies (KNOX, WOX, and ZF-HD) through machine-learning models for the detection of TF targets across different tissues and stages by combining tsCUT&Tag, ATAC-seq, and expression profiling. This dynamic regulome demonstrates the distinct binding specificity and co-factors for these homeobox subfamilies, indicative of functional divergence between and within them. Furthermore, we assemble a SAM dynamic regulome, illustrating potential functional mechanisms associated with plant architecture. Lastly, we generate a wox13a mutant that provides evidence that WOX13A directly regulates Gn1 expression to modulate plant height, validating the regulome of SAM-related homeobox genes.

CONCLUSIONS

The SAM-related homeobox transcription-factor regulome presents an unprecedented opportunity to dissect the molecular mechanisms governing SAM maintenance and development, thereby advancing our understanding of maize growth and shoot architecture.

Combined Metabolome and Transcriptome Analyses of Maize Leaves Reveal Global Effect of Biochar on Mechanisms Involved in Anti-Herbivory to Spodoptera frugiperda

Fall armyworm (FAW, Spodoptera frugiperda) has now spread to more than 26 Chinese provinces. The government is working with farmers and researchers to find ways to prevent and control this pest. The use of biochar is one of the economic and environmentally friendly strategies to increase plant growth and improve pest resistance. We tested four v/v combinations of bamboo charcoal with coconut bran [BC1 (10:1), BC2(30:1), BC3(50:1)] against a control (CK) in maize. We found that plant height, stem thickness, fresh weight and chlorophyll content were significantly higher in BC2, in addition to the lowest FAW survival %. We then compared the metabolome and transcriptome profiles of BC2 and CK maize plants under FAW herbivory. Our results show that the levels of flavonoids, amino acids and derivatives, nucleotides and derivatives and most phenolic acids decreased, while terpenoids, organic acids, lipids and defense-related hormones increased in BC-grown maize leaves. Transcriptome sequencing revealed consistent expression profiles of genes enriched in these pathways. We also observed the increased expression of genes related to abscisic acid, jasmonic acid, auxin and MAPK signaling. Based on these observations, we discussed the possible pathways involved in maize against FAW herbivory. We conclude that bamboo charcoal induces anti-herbivory responses in maize leaves.