Genome sequencing reveals evidence of adaptive variation in the genus Zea

Pan-genome analyses of 11 Fraxinus species provide insights into salt adaptation in ash trees

Ash trees (Fraxinus) exhibit rich genetic diversity and wide adaptation to various ecological environments, and several species are highly salt tolerant. Dissecting the genomic basis of salt adaptation in Fraxinus is vital for its resistance breeding. Here, we present 11 high-quality chromosome-level genome assemblies for Fraxinus species, which reveal two unequal subgenome compositions and two recent whole-genome triplication events in their evolutionary history. A Fraxinus pan-genome was constructed on the basis of structural variations and revealed that presence-absence variations (PAVs) of transmembrane transport genes have likely contributed to salt adaptation in Fraxinus. Through whole-genome resequencing of an F1 population from an interspecies cross of F. velutina 'Lula 3' (salt tolerant) with F. pennsylvanica 'Lula 5' (salt sensitive), we mapped salt-tolerance PAV-based quantitative trait loci (QTLs) and pinpointed two PAV-QTLs and candidate genes associated with Fraxinus salt tolerance. Mechanistically, FvbHLH85 enhances salt tolerance by mediating reactive oxygen species and Na+/K+ homeostasis, whereas FvSWEET5 enhances salt tolerance by mediating osmotic homeostasis. Collectively, these findings provide valuable genomic resources for Fraxinus salt-resistance breeding and the research community.

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.

MaizeCODE reveals bi-directionally expressed enhancers that harbor molecular signatures of maize domestication

Modern maize (Zea mays ssp. mays) was domesticated from Teosinte parviglumis (Zea mays ssp. parviglumis), with subsequent introgressions from Teosinte mexicana (Zea mays ssp. mexicana), yielding increased kernel row number, loss of the hard fruit case and dissociation from the cob upon maturity, as well as fewer tillers. Molecular approaches have identified transcription factors controlling these traits, yet revealed that a complex regulatory network is at play. MaizeCODE deploys ENCODE strategies to catalog regulatory regions in the maize genome, generating histone modification and transcription factor ChIP-seq in parallel with transcriptomics datasets in 5 tissues of 3 inbred lines which span the phenotypic diversity of maize, as well as the teosinte inbred TIL11. Transcriptomic analysis reveals that pollen grains share features with endosperm, and express dozens of "proto-miRNAs" potential vestiges of gene drive and hybrid incompatibility. Integrated analysis with chromatin modifications results in the identification of a comprehensive set of regulatory regions in each tissue of each inbred, and notably of distal enhancers expressing non-coding enhancer RNAs bi-directionally, reminiscent of "super enhancers" in animal genomes. Furthermore, the morphological traits selected during domestication are recapitulated, both in gene expression and within regulatory regions containing enhancer RNAs, while highlighting the conflict between enhancer activity and silencing of the neighboring transposable elements.

Genome-wide identification, phylogeny, evolutionary expansion, and expression analyses of ABC gene family in Castanea mollissima under temperature stress

The ATP-binding cassette (ABC) gene family comprises some of the most critical transporter proteins in plants, playing vital roles in maintaining cellular homeostasis and adapting to environmental changes. While ABC transporters have been extensively characterized in various plant species, their profile in C. mollissima remains less understood. In this study, 164 ABC genes were identified and characterized within the C. mollissima genome, and subsequently classified into eight subfamilies. Collinear analysis suggested that dispersed duplication was the primary mechanism driving the expansion of the CmABC gene family. The study also examined morphological and physiological changes in C. mollissima under temperature stress, highlighting significant decreases in photosynthetic indicators and SOD enzyme activity, while other indicators varied. Transcriptome analysis revealed distinct expression patterns of various CmABC genes under temperature stress, identifying CmABCG29c and CmABCB11e as key candidates for responding to temperature stress. This was based on their expression patterns, correlation with physiological indicators, and WGCNA analysis. The expression levels of CmABC genes measured in RT-qPCR experiments were consistent with those observed in RNA-seq analysis. This research provides a theoretical foundation for understanding the physiological and gene expression responses of C. mollissima to temperature stress.

Genetic dissection of ear-related trait divergence between maize and teosinte

Maize has undergone remarkable domestication and shows striking differences in architecture and ear morphology compared to its wild progenitor, called teosinte. However, our understanding of the genetic mechanisms underlying the ear morphology differences between teosinte and cultivated maize is still limited. In this study, we explored the genetic basis of ear-related traits at both early and mature stages by analyzing a population derived from a cross between Mo17 and a teosinte line, mexicana. We identified 31 quantitative trait loci (QTLs) associated with four IM-related and four ear-related traits, with 27 QTLs subjected to selection during the domestication process. Several key genes related to ear development were found under selection, including KN1 and RA1. Analysis of gene expression in the IM of developing ears from the population revealed the prominent roles of cis-variants in gene regulation. We also identified a large number of trans-eQTLs responsible for gene expression variation, and enrichment analysis on a trans-eQTL hotspot revealed the possible involvement of the sulfur metabolic pathway in controlling ear traits. Integrating the expression and phenotypic mapping data, we pinpointed several candidate genes potentially influencing ear development. Our findings advance the understanding of the genetic basis driving ear trait variation during maize domestication.

A genomic journey across the past, present, and future of South American maize

Maize (Zea mays ssp. mays) diverged from one of its wild relatives, the teosinte Zea mays ssp. parviglumis, in the lowlands of southwest Mexico ~9000 years ago. Following this divergence, maize rapidly expanded throughout the Americas becoming a staple food. This dispersal was accompanied by significant demographic and selective changes, leading to the development of numerous local varieties with a complex evolutionary history that remains incompletely understood. In recent years, genomic advances have challenged traditional models of maize domestication and spread to South America. At least three distinct genetic lineages associated with different migratory waves have been described: ancestral Andean, ancestral Lowland, and Pan-American. Additionally, the significant role of the teosinte Zea mays ssp. mexicana in the evolution of modern maize has been recently uncovered. Genomic studies have shed light on highland adaptation processes, revealing largely independent adaptation events in Meso and South America. As new evidence emerges, the regional complexity underlying maize diversity and the need for comprehensive, multiscale approaches have become evident. In the face of climate change and evolving agricultural landscapes, the conservation of native maize in South America is of growing interest, with genomics serving as an invaluable tool for identifying and preserving the genetic variability of locally adapted germplasm.

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.

Genome-Wide Diversity in Lowland and Highland Maize Landraces From Southern South America: Population Genetics Insights to Assist Conservation

Maize (Zea mays ssp. mays L.) landraces are traditional American crops with high genetic variability that conform a source of original alleles for conventional maize breeding. Northern Argentina, one the southernmost regions of traditional maize cultivation in the Americas, harbours around 57 races traditionally grown in two regions with contrasting environmental conditions, namely, the Andean mountains in the Northwest and the tropical grasslands and Atlantic Forest in the Northeast. These races encounter diverse threats to their genetic diversity and persistence in their regions of origin, with climate change standing out as one of the major challenges. In this work, we use genome-wide SNPs derived from ddRADseq to study the genetic diversity of individuals representing the five groups previously described for this area. This allowed us to distinguish two clearly differentiated gene pools, the highland northwestern maize (HNWA) and the floury northeastern maize (FNEA). Subsequently, we employed essential biodiversity variables at the genetic level, as proposed by the Group on Earth Observations Biodiversity Observation Network (GEO BON), to evaluate the conservation status of these two groups. This assessment encompassed genetic diversity (Pi), inbreeding coefficient (F) and effective population size (Ne). FNEA showed low Ne values and high F values, while HNWA showed low Ne values and low Pi values, indicating that further genetic erosion is imminent for these landraces. Outlier detection methods allowed identification of putative adaptive genomic regions, consistent with previously reported flowering-time loci and chromosomal regions displaying introgression from the teosinte Zea mays ssp. mexicana. Finally, species distribution models were obtained for two future climate scenarios, showing a notable reduction in the potential planting area of HNWA and a shift in the cultivation areas of FNEA. These results suggest that maize landraces from Northern Argentina may be unable to cope with climate change. Therefore, active conservation policies are advisable.

Genetic engineering, including genome editing, for enhancing broad-spectrum disease resistance in crops

Plant diseases, caused by a wide range of pathogens, severely reduce crop yield and quality, posing a significant threat to global food security. Developing broad-spectrum resistance (BSR) in crops is a key strategy for controlling crop diseases and ensuring sustainable crop production. Cloning disease-resistance (R) genes and understanding their underlying molecular mechanisms provide new genetic resources and strategies for crop breeding. Novel genetic engineering and genome editing tools have accelerated the study and engineering of BSR genes in crops, which is the primary focus of this review. We first summarize recent advances in understanding the plant immune system, followed by an examination of the molecular mechanisms underlying BSR in crops. Finally, we highlight diverse strategies employed to achieve BSR, including gene stacking to combine multiple R genes, multiplexed genome editing of susceptibility genes and promoter regions of executor R genes, editing cis-regulatory elements to fine-tune gene expression, RNA interference, saturation mutagenesis, and precise genomic insertions. The genetic studies and engineering of BSR are accelerating the breeding of disease-resistant cultivars, contributing to crop improvement and enhancing global food security.

Elymusmultiramosus (Poaceae), a new species from the north-western Qinghai-Tibetan Plateau, China

A new species from China, Elymusmultiramosus Y.C. Zhang, sp. nov. is described and illustrated herein, based on morphological characters and molecular phylogenetic analysis. The taxonomic descriptions of E.multiramosus and the comparison with related species are presented. The taxonomic distinctiveness of this new species was inferred by Maximum Likelihood (ML) analysis and Bayesian phylogenetic analysis, based on the complete chloroplast genome sequence. It is assigned to the Elymus section and bears similarity to Elymusnutans Griseb. However, it can be easily distinguished from other species by its compound spike, in contrast to the simple spike inflorescence typical of those species. The compound spike is characterised by rhachillas that are extended at the base of the main axis, giving rise to 3-6 mini-spike-like branches. Notably, these branches significantly increase in length from the top towards the bottom of the compound spike. In the molecular phylogeny, Elymusmultiramosus from Qinghai, north-western China, is phylogenetically positioned as a distinct lineage. The lineage comprising Elymussinosubmuticus from Sichuan, east of the Tibetan Plateau and Elymusnutans from the Himalayas forms a sister group to Elymusmultiramosus, suggesting that these three species share a common ancestor that is distinct from the lineage leading to Elymusatratus from Gansu, north of the Tibetan Plateau.

Genomic insights into adaptive evolution of the species-rich cosmopolitan plant genus Rhododendron

The species-rich cosmopolitan genus Rhododendron offers a good system for exploring the genomic mechanisms underlying adaptation to diverse habitats. Here, we report high-quality chromosomal-level genome assemblies of nine species, representing all five subgenera, different altitudinal distributions, and all flower color types of this genus. Further comprehensive genomic analyses indicate diverse adaptive strategies employed by Rhododendron, particularly adaptation to alpine and subalpine habitats by expansion/contraction of gene families involved in pathogen defense and oxidative phosphorylation, genomic convergent evolution, and gene copy-number variation. The convergent adaptation to high altitudes is further shown by population genomic analysis of R. nivale from the Himalaya-Hengduan Mountains. Moreover, we identify the genes involved in the biosynthesis of anthocyanins and carotenoids, which play a crucial role in shaping flower color diversity and environmental adaptation. Our study is significant for comprehending plant adaptive evolution and the uneven distribution of species diversity across different geographical regions.

Genomic variation of 363 diverse tea accessions unveils the genetic diversity, domestication, and structural variations associated with tea adaptation

Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

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.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Selfish genetic elements contribute to hybrid incompatibility and bias or 'drive' their own transmission1,2. Chromosomal drive typically functions in asymmetric female meiosis, whereas gene drive is normally post-meiotic and typically found in males. Here, using single-molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Z. mays ssp. mexicana) that depends on RNA interference (RNAi). 22-nucleotide small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1 and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas3, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize4. A survey of maize traditional varieties and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least four chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive probably had a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of 'self' small RNAs in the germ lines of plants and animals.

Linkage and association mapping in multi-parental populations reveal the genetic basis of carotenoid variation in maize kernels

Carotenoids are indispensable to plants and critical components of the human diet. The carotenoid metabolic pathway is conserved across plant species, but our understanding of the genetic basis of carotenoid variation remains limited for the seeds of most cereal crops. To address this issue, we systematically performed linkage and association mapping for eight carotenoid traits using six recombinant inbred line (RIL) populations. Single linkage mapping (SLM) and joint linkage mapping (JLM) identified 77 unique additive QTLs and 104 pairs of epistatic QTLs. Among these QTLs, we identified 22 overlapping hotspots of additive and epistatic loci, highlighting the important contributions of some QTLs to carotenoid levels through additive or epistatic mechanisms. A genome-wide association study based on all RILs detected 244 candidate genes significantly associated with carotenoid traits, 23 of which were annotated as carotenoid pathway genes. Effect comparisons suggested that a small number of loci linked to pathway genes have substantial effects on carotenoid variation in our tested populations, but many loci not associated with pathway genes also make important contributions to carotenoid variation. We identified ZmPTOX as the causal gene for a QTL hotspot (Q10/JLM10/GWAS019); this gene encodes a putative plastid terminal oxidase that produces plastoquinone-9 used by two enzymes in the carotenoid pathway. Natural variants in the promoter and second exon of ZmPTOX were found to alter carotenoid levels. This comprehensive assessment of the genetic mechanisms underlying carotenoid variation establishes a foundation for rewiring carotenoid metabolism and accumulation for efficient carotenoid biofortification.

Evolution and Genetic Differentiation of Pleurotus tuoliensis in Xinjiang, China, Based on Population Genomics

Pleurotus tuoliensis is a unique species discovered in Xinjiang, China, which is recognized for its significant edible, medicinal, and economic value. It has been successfully incorporated into industrial production. Controversy has emerged concerning the evolution and environmental adaptability of this species due to inadequate interspecific ecology and molecular data. This study examines the germplasm resources of P. tuoliensis in the Xinjiang region. A total of 225 wild and cultivated strains of P. tuoliensis were gathered from seven representative regions. Phylogenetic analysis revealed that seven populations were notably segregated into three distinct groups, primarily attributed to environmental factors as the underlying cause for this differentiation. Population historical size data indicate that P. tuoliensis underwent two expansion events, one between 2 and 0.9 Mya (Miocene) and the other between 15 and 4 Mya (Early Pleistocene). The ancient climate fluctuations in the Xinjiang region might have contributed to the comparatively modest population size during the Pliocene epoch. Moreover, through the integration of biogeography and ancestral state reconstruction, it was determined that group C of P. tuoliensis emerged initially and subsequently dispersed to groups D and B, in that order. Subsequently, group D underwent independent evolution, whereas group B continued to diversify into groups A and EFG. The primary factor influencing this mode of transmission route is related to the geographical conditions and prevailing wind direction of each group. Subsequent research endeavors focused on assessing the domestication adaptability of P. tuoliensis to different substrates. It was found that the metabolic processes adapted during the domestication process were mainly related to energy metabolism, DNA repair, and environmental adaptability. Processes adapted to the host adaptability include responses to the host (meiosis, cell cycle, etc.) and stress in the growth environment (cysteine and methionine metabolism, sulfur metabolism, etc.). This study analyzed the systematic evolution and genetic differentiation of P. tuoliensis in Xinjiang. The identified loci and genes provide a theoretical basis for the subsequent improvement of germplasm resources and conducting molecular breeding.

Population-level gene expression can repeatedly link genes to functions in maize

Transcriptome-wide association studies (TWAS) can provide single gene resolution for candidate genes in plants, complementing genome-wide association studies (GWAS) but efforts in plants have been met with, at best, mixed success. We generated expression data from 693 maize genotypes, measured in a common field experiment, sampled over a 2-h period to minimize diurnal and environmental effects, using full-length RNA-seq to maximize the accurate estimation of transcript abundance. TWAS could identify roughly 10 times as many genes likely to play a role in flowering time regulation as GWAS conducted data from the same experiment. TWAS using mature leaf tissue identified known true-positive flowering time genes known to act in the shoot apical meristem, and trait data from a new environment enabled the identification of additional flowering time genes without the need for new expression data. eQTL analysis of TWAS-tagged genes identified at least one additional known maize flowering time gene through trans-eQTL interactions. Collectively these results suggest the gene expression resource described here can link genes to functions across different plant phenotypes expressed in a range of tissues and scored in different experiments.

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.

Genomic variation in weedy and cultivated broomcorn millet accessions uncovers the genetic architecture of agronomic traits

Large-scale genomic variations are fundamental resources for crop genetics and breeding. Here we sequenced 1,904 genomes of broomcorn millet to an average of 40× sequencing depth and constructed a comprehensive variation map of weedy and cultivated accessions. Being one of the oldest cultivated crops, broomcorn millet has extremely low nucleotide diversity and remarkably rapid decay of linkage disequilibrium. Genome-wide association studies identified 186 loci for 12 agronomic traits. Many causative candidate genes, such as PmGW8 for grain size and PmLG1 for panicle shape, showed strong selection signatures during domestication. Weedy accessions contained many beneficial variations for the grain traits that are largely lost in cultivated accessions. Weedy and cultivated broomcorn millet have adopted different loci controlling flowering time for regional adaptation in parallel. Our study uncovers the unique population genomic features of broomcorn millet and provides an agronomically important resource for cereal crops.

Upland rice genomic signatures of adaptation to drought resistance and navigation to molecular design breeding

Upland rice is a distinctive drought-aerobic ecotype of cultivated rice highly resistant to drought stress. However, the genetic and genomic basis for the drought-aerobic adaptation of upland rice remains largely unclear due to the lack of genomic resources. In this study, we identified 25 typical upland rice accessions and assembled a high-quality genome of one of the typical upland rice varieties, IRAT109, comprising 384 Mb with a contig N50 of 19.6 Mb. Phylogenetic analysis revealed upland and lowland rice have distinct ecotype differentiation within the japonica subgroup. Comparative genomic analyses revealed that adaptive differentiation of lowland and upland rice is likely attributable to the natural variation of many genes in promoter regions, formation of specific genes in upland rice, and expansion of gene families. We revealed differentiated gene expression patterns in the leaves and roots of the two ecotypes and found that lignin synthesis mediated by the phenylpropane pathway plays an important role in the adaptive differentiation of upland and lowland rice. We identified 28 selective sweeps that occurred during domestication and validated that the qRT9 gene in selective regions can positively regulate drought resistance in rice. Eighty key genes closely associated with drought resistance were appraised for their appreciable potential in drought resistance breeding. Our study enhances the understanding of the adaptation of upland rice and provides a genome navigation map of drought resistance breeding, which will facilitate the breeding of drought-resistant rice and the "blue revolution" in agriculture.

Genomic insight into the origin, domestication, dispersal, diversification and human selection of Tartary buckwheat

BACKGROUND

Tartary buckwheat, Fagopyrum tataricum, is a pseudocereal crop with worldwide distribution and high nutritional value. However, the origin and domestication history of this crop remain to be elucidated.

RESULTS

Here, by analyzing the population genomics of 567 accessions collected worldwide and reviewing historical documents, we find that Tartary buckwheat originated in the Himalayan region and then spread southwest possibly along with the migration of the Yi people, a minority in Southwestern China that has a long history of planting Tartary buckwheat. Along with the expansion of the Mongol Empire, Tartary buckwheat dispersed to Europe and ultimately to the rest of the world. The different natural growth environments resulted in adaptation, especially significant differences in salt tolerance between northern and southern Chinese Tartary buckwheat populations. By scanning for selective sweeps and using a genome-wide association study, we identify genes responsible for Tartary buckwheat domestication and differentiation, which we then experimentally validate. Comparative genomics and QTL analysis further shed light on the genetic foundation of the easily dehulled trait in a particular variety that was artificially selected by the Wa people, a minority group in Southwestern China known for cultivating Tartary buckwheat specifically for steaming as a staple food to prevent lysine deficiency.

CONCLUSIONS

This study provides both comprehensive insights into the origin and domestication of, and a foundation for molecular breeding for, Tartary buckwheat.

Two teosintes made modern maize

The origins of maize were the topic of vigorous debate for nearly a century, but neither the current genetic model nor earlier archaeological models account for the totality of available data, and recent work has highlighted the potential contribution of a wild relative, Zea mays ssp. mexicana. Our population genetic analysis reveals that the origin of modern maize can be traced to an admixture between ancient maize and Zea mays ssp. mexicana in the highlands of Mexico some 4000 years after domestication began. We show that variation in admixture is a key component of maize diversity, both at individual loci and for additive genetic variation underlying agronomic traits. Our results clarify the origin of modern maize and raise new questions about the anthropogenic mechanisms underlying dispersal throughout the Americas.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Meiotic drivers subvert Mendelian expectations by manipulating reproductive development to bias their own transmission. Chromosomal drive typically functions in asymmetric female meiosis, while gene drive is normally postmeiotic and typically found in males. Using single molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Zea mays ssp. mexicana), that depends on RNA interference (RNAi). 22nt small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-Like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1, and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize. A survey of maize landraces and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least 4 chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive likely played a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of "self" small RNAs in the germlines of plants and animals.

The Ga1 locus of the genus Zea is associated with novel genome structures derived from multiple, independent nonhomologous recombination events

The Ga1 locus controls cross-incompatibility between field corn and popcorn. The Ga1-S haplotype contains 2 types of pectin methylesterase (PME) genes, ZmPme3 and several copies of ZmGa1P that are expressed in silk and pollen, respectively. The ga1 haplotype contains nonfunctional tandem repeat sequences related to ZmPme3 and ZmGa1P. This haplotype can cross-pollinate freely and is widely present in field corn. The primary objective of this study is to characterize the repeat sequences from a diverse collection of maize and teosinte lines and use this information to understand the evolution of the Ga1 locus. First, we characterized the complexity of the Ga1 genome region in high-quality maize genome assemblies that led to their categorization into 5 groups based on the number and type of PME-like sequences found at this region. Second, we studied duplication events that led to the ga1 and Ga1-S repeats using maximum likelihood phylogenetic reconstruction. Divergence estimates of the ga1 haplotype suggest that the duplication events occurred more than 600 KYA whereas those in Ga1-S occurred at 3 time points, i.e. >600, ∼260, and ∼100 KYA. These estimates suggest that the ga1 and Ga1-S tandem duplication events occurred independently. Finally, analysis of ZmPme3 and ZmGa1P homologs in Zea and Tripsacum genomes suggests that ga1 and Ga1-S repeats originated from an ancestral pair of PME genes that duplicated and diverged through 2 evolutionary branches prior to the domestication of maize.

ZmEREB92 plays a negative role in seed germination by regulating ethylene signaling and starch mobilization in maize

Rapid and uniform seed germination is required for modern cropping system. Thus, it is important to optimize germination performance through breeding strategies in maize, in which identification for key regulators is needed. Here, we characterized an AP2/ERF transcription factor, ZmEREB92, as a negative regulator of seed germination in maize. Enhanced germination in ereb92 mutants is contributed by elevated ethylene signaling and starch degradation. Consistently, an ethylene signaling gene ZmEIL7 and an α-amylase gene ZmAMYa2 are identified as direct targets repressed by ZmEREB92. OsERF74, the rice ortholog of ZmEREB92, shows conserved function in negatively regulating seed germination in rice. Importantly, this orthologous gene pair is likely experienced convergently selection during maize and rice domestication. Besides, mutation of ZmEREB92 and OsERF74 both lead to enhanced germination under cold condition, suggesting their regulation on seed germination might be coupled with temperature sensitivity. Collectively, our findings uncovered the ZmEREB92-mediated regulatory mechanism of seed germination in maize and provide breeding targets for maize and rice to optimize seed germination performance towards changing climates.

Genetic and molecular exploration of maize environmental stress resilience: Toward sustainable agriculture

Global climate change exacerbates the effects of environmental stressors, such as drought, flooding, extreme temperatures, salinity, and alkalinity, on crop growth and grain yield, threatening the sustainability of the food supply. Maize (Zea mays) is one of the most widely cultivated crops and the most abundant grain crop in production worldwide. However, the stability of maize yield is highly dependent on environmental conditions. Recently, great progress has been made in understanding the molecular mechanisms underlying maize responses to environmental stresses and in developing stress-resilient varieties due to advances in high-throughput sequencing technologies, multi-omics analysis platforms, and automated phenotyping facilities. In this review, we summarize recent advances in dissecting the genetic factors and networks that contribute to maize abiotic stress tolerance through diverse strategies. We also discuss future challenges and opportunities for the development of climate-resilient maize varieties.

Synergistic optimization of crops by combining early maturation with other agronomic traits

Many newly created early maturing varieties exhibit poor stress resistance and low yield, whereas stress-resistant varieties are typically late maturing. For this reason, the polymerization of early maturity and other desired agronomic qualities requires overcoming the negative connection between early maturity, multi-resistance, and yield, which presents a formidable challenge in current breeding techniques. We review the most salient constraints of early maturity breeding in current crop planting practices and the molecular mechanisms of different maturation timeframes in diverse crops from their origin center to production areas. We explore current breeding tactics and the future direction of crop breeding and the issues that must be resolved to accomplish the polymerization of desirable traits in light of the current obstacles and limitations.

The genetic architecture of prolificacy in maize revealed by association mapping and bulk segregant analysis

KEY MESSAGE

Here, we revealed maize prolificacy highly correlated with domestication and identified a causal gene ZmEN1 located in one novel QTL qGEN261 that regulating maize prolificacy by using multiple-mapping methods. The development of maize prolificacy (EN) is crucial for enhancing yield and breeding specialty varieties. To achieve this goal, we employed a genome-wide association study (GWAS) to analyze the genetic architecture of EN in maize. Using 492 inbred lines with a wide range of EN variability, our results demonstrated significant differences in genetic, environmental, and interaction effects. The broad-sense heritability (H2) of EN was 0.60. Through GWAS, we identified 527 significant single nucleotide polymorphisms (SNPs), involved 290 quantitative trait loci (QTL) and 806 genes. Of these SNPs, 18 and 509 were classified as major effect loci and minor loci, respectively. In addition, we performed a bulk segregant analysis (BSA) in an F2 population constructed by a few-ears line Zheng58 and a multi-ears line 647. Our BSA results identified one significant QTL, qBEN1. Importantly, combining the GWAS and BSA, four co-located QTL, involving six genes, were identified. Three of them were expressed in vegetative meristem, shoot tip, internode and tip of ear primordium, with ZmEN1, encodes an unknown auxin-like protein, having the highest expression level in these tissues. It suggested that ZmEN1 plays a crucial role in promoting axillary bud and tillering to encourage the formation of prolificacy. Haplotype analysis of ZmEN1 revealed significant differences between different haplotypes, with inbred lines carrying hap6 having more EN. Overall, this is the first report about using GWAS and BSA to dissect the genetic architecture of EN in maize, which can be valuable for breeding specialty maize varieties and improving maize yield.

Unveiling the Effects of Carbon-Based Nanomaterials on Crop Growth: From Benefits to Detriments

To systematically assess the impact of typical CNMs on the growth effects of cereal crops, we conducted a meta-analysis of 48 independent studies worldwide. The pooled results showed that shoot weight (13.39%) and antioxidant metabolite content (SOD: 106.32%, POD: 32.29%, CAT: 22.63%) of cereal crops exposed to the presence of CNMs were significantly increased, but phytohormones secretion (17.84%) was inhibited. The results of subgroup analysis showed that there were differences in the results of different CNM types with the same exposure concentration on growth effects. Short-term exposure adversely affected the root and photosynthetic capacity of the crop, but prolonged exposure instead showed a promoting effect. Multiple linear regression analysis showed that the concentration of CNMs and cereal variety variables were significantly associated with changes in multiple growth effect values. This work could offer references and fresh perspectives for investigating how nanoparticles and crops interact.

CIMBL55: a repository for maize drought resistance alleles

Droughts threaten crop yields worldwide. Compared to other major staple cereal crops, maize (Zea mays) is especially sensitive to drought, which can cause dramatic fluctuations in its yield potential. Natural maize populations contain many superior alleles that can enhance drought resistance through complex regulatory mechanisms. We recently de novo assembled the genome of a prominent drought-resistant maize germplasm, CIMBL55, and systematically dissected the genetic basis for its drought resistance on the genome, transcriptome, and epigenome levels. These analyses revealed 65 favorable drought resistance alleles in CIMBL55. Subsequently, we genetically verified the functions of the drought resistance genes ZmABF4, ZmNAC075, and ZmRtn16 and unraveled the function of ZmRtn16 on a molecular level.

Novel insights into maize (Zea mays) development and organogenesis for agricultural optimization

In maize, intrinsic hormone activities and sap fluxes facilitate organogenesis patterning and plant holistic development; these hormone movements should be a primary focus of developmental biology and agricultural optimization strategies. Maize (Zea mays) is an important crop plant with distinctive life history characteristics and structural features. Genetic studies have extended our knowledge of maize developmental processes, genetics, and molecular ecophysiology. In this review, the classical life cycle and life history strategies of maize are analyzed to identify spatiotemporal organogenesis properties and develop a definitive understanding of maize development. The actions of genes and hormones involved in maize organogenesis and sex determination, along with potential molecular mechanisms, are investigated, with findings suggesting central roles of auxin and cytokinins in regulating maize holistic development. Furthermore, investigation of morphological and structural characteristics of maize, particularly node ubiquity and the alternate attachment pattern of lateral organs, yields a novel regulatory model suggesting that maize organ initiation and subsequent development are derived from the stimulation and interaction of auxin and cytokinin fluxes. Propositions that hormone activities and sap flow pathways control organogenesis are thoroughly explored, and initiation and development processes of distinctive maize organs are discussed. Analysis of physiological factors driving hormone and sap movement implicates cues of whole-plant activity for hormone and sap fluxes to stimulate maize inflorescence initiation and organ identity determination. The physiological origins and biogenetic mechanisms underlying maize floral sex determination occurring at the tassel and ear spikelet are thoroughly investigated. The comprehensive outline of maize development and morphogenetic physiology developed in this review will enable farmers to optimize field management and will provide a reference for de novo crop domestication and germplasm improvement using genome editing biotechnologies, promoting agricultural optimization.

Large haplotypes highlight a complex age structure within the maize pan-genome

The genomes of maize and other eukaryotes contain stable haplotypes in regions of low recombination. These regions, including centromeres, long heterochromatic blocks, and rDNA arrays, have been difficult to analyze with respect to their diversity and origin. Greatly improved genome assemblies are now available that enable comparative genomics over these and other nongenic spaces. Using 26 complete maize genomes, we developed methods to align intergenic sequences while excluding genes and regulatory regions. The centromere haplotypes (cenhaps) extend for megabases on either side of the functional centromere regions and appear as evolutionary strata, with haplotype divergence/coalescence times dating as far back as 450 thousand years ago (kya). Application of the same methods to other low recombination regions (heterochromatic knobs and rDNA) and all intergenic spaces revealed that deep coalescence times are ubiquitous across the maize pan-genome. Divergence estimates vary over a broad timescale with peaks at ∼16 and 300 kya, reflecting a complex history of gene flow among diverging populations and changes in population size associated with domestication. Cenhaps and other long haplotypes provide vivid displays of this ancient diversity.

A common resequencing-based genetic marker data set for global maize diversity

Maize (Zea mays ssp. mays) populations exhibit vast ranges of genetic and phenotypic diversity. As sequencing costs have declined, an increasing number of projects have sought to measure genetic differences between and within maize populations using whole-genome resequencing strategies, identifying millions of segregating single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels). Unlike older genotyping strategies like microarrays and genotyping by sequencing, resequencing should, in principle, frequently identify and score common genetic variants. However, in practice, different projects frequently employ different analytical pipelines, often employ different reference genome assemblies and consistently filter for minor allele frequency within the study population. This constrains the potential to reuse and remix data on genetic diversity generated from different projects to address new biological questions in new ways. Here, we employ resequencing data from 1276 previously published maize samples and 239 newly resequenced maize samples to generate a single unified marker set of approximately 366 million segregating variants and approximately 46 million high-confidence variants scored across crop wild relatives, landraces as well as tropical and temperate lines from different breeding eras. We demonstrate that the new variant set provides increased power to identify known causal flowering-time genes using previously published trait data sets, as well as the potential to track changes in the frequency of functionally distinct alleles across the global distribution of modern maize.

Identification of a new QTL underlying seminal root number in a maize-teosinte population

Seminal roots play an important role in acquisition of water and nutrients by maize seedlings. Compared with its teosinte ancestor, maize underwent a change in seminal root number (SRN). Although several key genes controlling SRN have been cloned, identification and utilization of new genes from teosinte would be useful for improving maize root architecture. In this study, a maize-teosinte BC₂F₆ population containing 206 individuals genotyped by resequencing was used to conduct high-resolution quantitative trait locus (QTL) mapping of SRN. A new major QTL on chromosome 7 (qSRN7) was identified. Differentially expressed genes (DEGs) based on RNA-Seq were identified between two inbred lines with no SRN and multiple SRN at two periods of seminal roots primordia formation. A total of 116 DEGs detected in at least one period were identified within the qSRN7 interval. Three DEGs (Zm00001d021572, Zm00001d021579 and Zm00001d021861) associated with SRN were identified through regional association mapping. When compared with reported domestication-related selective sweeps, Zm00001d021572 was selected during maize domestication. Our findings provide important insights into the genetic basis of SRN and identify a promising candidate gene for further studies on SRN.

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

Integrated single-molecule real-time sequencing and RNA sequencing reveal the molecular mechanisms of salt tolerance in a novel synthesized polyploid genetic bridge between maize and its wild relatives

BACKGROUND

Tripsacum dactyloides (2n = 4x = 72) and Zea perennis (2n = 4x = 40) are tertiary gene pools of Zea mays L. and exhibit many abiotic adaptations absent in modern maize, especially salt tolerance. A previously reported allopolyploid (hereafter referred to as MTP, 2n = 74) synthesized using Zea mays, Tripsacum dactyloides, and Zea perennis has even stronger salt tolerance than Z. perennis and T. dactyloides. This allopolyploid will be a powerful genetic bridge for the genetic improvement of maize. However, the molecular mechanisms underlying its salt tolerance, as well as the key genes involved in regulating its salt tolerance, remain unclear.

RESULTS

Single-molecule real-time sequencing and RNA sequencing were used to identify the genes involved in salt tolerance and reveal the underlying molecular mechanisms. Based on the SMRT-seq results, we obtained 227,375 reference unigenes with an average length of 2300 bp; most of the unigenes were annotated to Z. mays sequences (76.5%) in the NR database. Moreover, a total of 484 and 1053 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Functional enrichment analysis of DEGs revealed that multiple pathways responded to salt stress, including "Flavonoid biosynthesis," "Oxidoreductase activity," and "Plant hormone signal transduction" in the leaves and roots, and "Iron ion binding," "Acetyl-CoA carboxylase activity," and "Serine-type carboxypeptidase activity" in the roots. Transcription factors, such as those in the WRKY, B3-ARF, and bHLH families, and cytokinin negatively regulators negatively regulated the salt stress response. According to the results of the short time series-expression miner analysis, proteins involved in "Spliceosome" and "MAPK signal pathway" dynamically responded to salt stress as salinity changed. Protein-protein interaction analysis revealed that heat shock proteins play a role in the large interaction network regulating salt tolerance.

CONCLUSIONS

Our results reveal the molecular mechanism underlying the regulation of MTP in the response to salt stress and abundant salt-tolerance-related unigenes. These findings will aid the retrieval of lost alleles in modern maize and provide a new approach for using T. dactyloides and Z. perennis to improve maize.

Root volatile profiles and herbivore preference are mediated by maize domestication, geographic spread, and modern breeding

Domestication affected the abundances and diversity of maize root volatiles more than northward spread and modern breeding, and herbivore preference for roots was correlated with volatile diversity and herbivore resistance. Studies show that herbivore defenses in crops are mediated by domestication, spread, and breeding, among other human-driven processes. They also show that those processes affected chemical communication between crop plants and herbivores. We hypothesized that (i) preference of the herbivore (Diabrotica virgifera virgifera) larvae for embryonic roots of maize (Zea mays mays) would increase and (ii) root volatile diversity would decrease with the crop's domestication, northward spread to present-day USA, and modern breeding. We used Balsas teosinte (Zea mays parviglumis), Mexican and USA landrace maizes, and US inbred maize lines to test these hypotheses. We found that herbivore preference and volatile diversity increased with maize domestication and northward spread but decreased with modern breeding. Additionally, we found that the abundances of single volatiles did not consistently increase or decrease with maize domestication, spread, and breeding; rather, volatiles grouped per their abundances were differentially affected by those processes, and domestication had the greatest effects. Altogether, our results suggested that: the herbivore's preference for maize roots is correlated with volatile diversity and herbivore resistance; changes in abundances of individual volatiles are evident at the level of volatile groups; and maize domestication, but not spread and breeding, affected the abundances of some green leaf volatiles and sesquiterpenes/sesquiterpenoids. In part, we discussed our results in the context of herbivore defense evolution when resources for plant growth and defense vary across environments. We suggested that variability in relative abundance of volatiles may be associated with their local, functional relevance across wild and agricultural environments.