Genomic, Transcriptomic, and Phenomic Variation Reveals the Complex Adaptation of Modern Maize Breeding

Genomic basis determining root system architecture in maize

KEY MESSAGE

A total of 389 and 344 QTLs were identified by GWAS and QTL mapping explaining accumulatively 32.2-65.0% and 23.7-63.4% of phenotypic variation for 14 shoot-borne root traits using more than 1300 individuals across multiple field trails. Efficient nutrient and water acquisition from soils depends on the root system architecture (RSA). However, the genetic determinants underlying RSA in maize remain largely unexplored. In this study, we conducted a comprehensive genetic analysis for 14 shoot-borne root traits using 513 inbred lines and 800 individuals from four recombinant inbred line (RIL) populations at the mature stage across multiple field trails. Our analysis revealed substantial phenotypic variation for these 14 root traits, with a total of 389 and 344 QTLs identified through genome-wide association analysis (GWAS) and linkage analysis, respectively. These QTLs collectively explained 32.2-65.0% and 23.7-63.4% of the trait variation within each population. Several a priori candidate genes involved in auxin and cytokinin signaling pathways, such as IAA26, ARF2, LBD37 and CKX3, were found to co-localize with these loci. In addition, a total of 69 transcription factors (TFs) from 27 TF families (MYB, NAC, bZIP, bHLH and WRKY) were found for shoot-borne root traits. A total of 19 genes including PIN3, LBD15, IAA32, IAA38 and ARR12 and 19 GWAS signals were overlapped with selective sweeps. Further, significant additive effects were found for root traits, and pyramiding the favorable alleles could enhance maize root development. These findings could contribute to understand the genetic basis of root development and evolution, and provided an important genetic resource for the genetic improvement of root traits in maize.

Identification of Novel QTL for Mercury Accumulation in Maize Using an Enlarged SNP Panel

Mercury (Hg) pollution not only poses a threat to the environment but also adversely affects the growth and development of plants, with potential repercussions for animals and humans through bioaccumulation in the food chain. Maize, a crucial source of food, industrial materials, and livestock feed, requires special attention in understanding the genetic factors influencing mercury accumulation. Developing maize varieties with low mercury accumulation is vital for both maize production and human health. In this study, a comprehensive genome-wide association study (GWAS) was conducted using an enlarged SNP panel comprising 1.25 million single nucleotide polymorphisms (SNPs) in 230 maize inbred lines across three environments. The analysis identified 111 significant SNPs within 78 quantitative trait loci (QTL), involving 169 candidate genes under the Q model. Compared to the previous study, the increased marker density and optimized statistical model led to the discovery of 74 additional QTL, demonstrating improved statistical power. Gene ontology (GO) enrichment analysis revealed that most genes participate in arsenate reduction and stress responses. Notably, GRMZM2G440968, which has been reported in previous studies, is associated with the significant SNP chr6.S_155668107 in axis tissue. It encodes a cysteine proteinase inhibitor, implying its potential role in mitigating mercury toxicity by inhibiting cysteine. Haplotype analyses provided further insights, indicating that lines carrying hap3 exhibited the lowest mercury content compared to other haplotypes. In summary, our study significantly enhances the statistical power of GWAS, identifying additional genes related to mercury accumulation and metabolism. These findings offer valuable insights into unraveling the genetic basis of mercury content in maize and contribute to the development of maize varieties with low mercury accumulation.

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

Background

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

Result

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

Conclusion

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

Supplementary information

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

Genetic basis of resistance to southern corn leaf blight in the maize multi-parent population and diversity panel

Southern corn leaf blight (SLB), caused by the necrotrophic pathogen Cochliobolus heterostrophus, is one of the maize foliar diseases and poses a great threat to corn production around the world. Identification of genetic variations underlying resistance to SLB is of paramount importance to maize yield and quality. Here, we used a random-open-parent association mapping population containing eight recombinant inbred line populations and one association mapping panel consisting of 513 diversity maize inbred lines with high-density genetic markers to dissect the genetic basis of SLB resistance. Overall, 109 quantitative trait loci (QTLs) with predominantly small or moderate additive effects, and little epistatic effects were identified. We found 35 (32.1%) novel loci in comparison with the reported QTLs. We revealed that resistant alleles were significantly enriched in tropical accessions and the frequency of about half of resistant alleles decreased during the adaptation process owing to the selection of agronomic traits. A large number of annotated genes located in the SLB-resistant QTLs were shown to be involved in plant defence pathways. Integrating genome-wide association study, transcriptomic profiling, resequencing and gene editing, we identified ZmFUT1 and MYBR92 as the putative genes responsible for the major QTLs for resistance to C. heterostrophus. Our results present a comprehensive insight into the genetic basis of SLB resistance and provide resistant loci or genes as direct targets for crop genetic improvement.

A teosinte-derived allele of an HKT1 family sodium transporter improves salt tolerance in maize

The sodium cation (Na⁺ ) is the predominant cation with deleterious effects on crops in salt-affected agricultural areas. Salt tolerance of crop can be improved by increasing shoot Na⁺ exclusion. Therefore, it is crucial to identify and use genetic variants of various crops that promote shoot Na⁺ exclusion. Here, we show that a HKT1 family gene ZmNC3 (Zea mays L. Na⁺ Content 3; designated ZmHKT1;2) confers natural variability in shoot-Na⁺ accumulation and salt tolerance in maize. ZmHKT1;2 encodes a Na⁺ -preferential transporter localized in the plasma membrane, which mediates shoot Na⁺ exclusion, likely by withdrawing Na⁺ from the root xylem flow. A naturally occurring nonsynonymous SNP (SNP947-G) increases the Na⁺ transport activity of ZmHKT1;2, promoting shoot Na⁺ exclusion and salt tolerance in maize. SNP947-G first occurred in the wild grass teosinte (at a allele frequency of 43%) and has become a minor allele in the maize population (allele frequency 6.1%), suggesting that SNP947-G is derived from teosinte and that the genomic region flanking SNP947 likely has undergone selection during domestication or post-domestication dispersal of maize. Moreover, we demonstrate that introgression of the SNP947-G ZmHKT1;2 allele into elite maize germplasms reduces shoot Na⁺ content by up to 80% and promotes salt tolerance. Taken together, ZmNC3/ZmHKT1;2 was identified as an important QTL promoting shoot Na⁺ exclusion, and its favourable allele provides an effective tool for developing salt-tolerant maize varieties.

De Novo Domestication in the Multi-Omics Era

Most cereal crops were domesticated within the last 12,000 years and subsequently spread around the world. These crops have been nourishing the world by supplying a primary energy and nutrient source, thereby playing a critical role in determining the status of human health and sustaining the global population. Here, we review the major challenges of future agriculture and emphasize the utilization of wild germplasm. De novo domestication is one of the most straightforward strategies to manipulate domestication-related and/or other genes with known function, and thereby introduce desired traits into wild plants. We also summarize known causal variations and their corresponding pathways in order to better understand the genetic basis of crop evolution, and how this knowledge could facilitate de novo domestication. Indeed knowledge-driven de novo domestication has great potential for the development of new sustainable crops that have climate-resilient high yield with low resource input and meet individual nutrient needs. Finally, we discuss current opportunities for and barriers to knowledge-driven de novo domestication.

Machine learning assisted dynamic phenotypes and genomic variants help understand the ecotype divergence in rapeseed

Three ecotypes of rapeseed, winter, spring, and semi-winter, have been formed to enable the plant to adapt to different geographic areas. Although several major loci had been found to contribute to the flowering divergence, the genomic footprints and associated dynamic plant architecture in the vegetative growth stage underlying the ecotype divergence remain largely unknown in rapeseed. Here, a set of 41 dynamic i-traits and 30 growth-related traits were obtained by high-throughput phenotyping of 171 diverse rapeseed accessions. Large phenotypic variation and high broad-sense heritability were observed for these i-traits across all developmental stages. Of these, 19 i-traits were identified to contribute to the divergence of three ecotypes using random forest model of machine learning approach, and could serve as biomarkers to predict the ecotype. Furthermore, we analyzed genomic variations of the population, QTL information of all dynamic i-traits, and genomic basis of the ecotype differentiation. It was found that 213, 237, and 184 QTLs responsible for the differentiated i-traits overlapped with the signals of ecotype divergence between winter and spring, winter and semi-winter, and spring and semi-winter, respectively. Of which, there were four common divergent regions between winter and spring/semi-winter and the strongest divergent regions between spring and semi-winter were found to overlap with the dynamic QTLs responsible for the differentiated i-traits at multiple growth stages. Our study provides important insights into the divergence of plant architecture in the vegetative growth stage among the three ecotypes, which was contributed to by the genetic differentiation, and might contribute to environmental adaption and yield improvement.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Network and Evolutionary Analysis Reveals Candidate Genes of Membrane Trafficking Involved in Maize Seed Development and Immune Response

The plant membrane-trafficking system plays a crucial role in maintaining proper cellular functions and responding to various developmental and environmental cues. Thus far, our knowledge of the maize membrane-trafficking system is still limited. In this study, we systematically identified 479 membrane-trafficking genes from the maize genome using orthology search and studied their functions by integrating transcriptome and evolution analyses. These genes encode the components of coated vesicles, AP complexes, autophagy, ESCRTs, retromers, Rab GTPases, tethering factors, and SNAREs. The maize genes exhibited diverse but coordinated expression patterns, with 249 genes showing elevated expression in reproductive tissues. Further WGCNA analysis revealed that five COPII components and four Rab GTPases had high connectivity with protein biosynthesis during endosperm development and that eight components of autophagy, ESCRT, Rab, and SNARE were strongly co-upregulated with defense-related genes and/or with secondary metabolic processes to confer basal resistance to Fusarium graminearum. In addition, we identified 39 membrane-trafficking genes with strong selection signals during maize domestication and/or improvement. Among them, ZmSec23a and ZmVPS37A were selected for kernel oil production during improvement and pathogen resistance during domestication, respectively. In summary, these findings will provide important hints for future appreciation of the functions of membrane-trafficking genes in maize.

High-Throughput Phenotyping Accelerates the Dissection of the Phenotypic Variation and Genetic Architecture of Shank Vascular Bundles in Maize (Zea mays L.)

The vascular bundle of the shank is an important ‘flow’ organ for transforming maize biological yield to grain yield, and its microscopic phenotypic characteristics and genetic analysis are of great significance for promoting the breeding of new varieties with high yield and good quality. In this study, shank CT images were obtained using the standard process for stem micro-CT data acquisition at resolutions up to 13.5 μm. Moreover, five categories and 36 phenotypic traits of the shank including related to the cross-section, epidermis zone, periphery zone, inner zone and vascular bundle were analyzed through an automatic CT image process pipeline based on the functional zones. Next, we analyzed the phenotypic variations in vascular bundles at the base of the shank among a group of 202 inbred lines based on comprehensive phenotypic information for two environments. It was found that the number of vascular bundles in the inner zone (IZ_VB_N) and the area of the inner zone (IZ_A) varied the most among the different subgroups. Combined with genome-wide association studies (GWAS), 806 significant single nucleotide polymorphisms (SNPs) were identified, and 1245 unique candidate genes for 30 key traits were detected, including the total area of vascular bundles (VB_A), the total number of vascular bundles (VB_N), the density of the vascular bundles (VB_D), etc. These candidate genes encode proteins involved in lignin, cellulose synthesis, transcription factors, material transportation and plant development. The results presented here will improve the understanding of the phenotypic traits of maize shank and provide an important phenotypic basis for high-throughput identification of vascular bundle functional genes of maize shank and promoting the breeding of new varieties with high yield and good quality.

Population genomics of Zea species identifies selection signatures during maize domestication and adaptation

BACKGROUND

Maize (Zea mays L. ssp. mays) was domesticated from teosinte (Zea mays ssp. parviglumis) about 9000 years ago in southwestern Mexico and adapted to a range of environments worldwide. Researchers have depicted the maize domestication and adaptation processes over the past two decades, but efforts have been limited either in sample size or genetic diversity. To better understand these processes, we conducted a genome-wide survey of 982 maize inbred lines and 190 teosinte accessions using over 40,000 single-nucleotide polymorphism markers.

RESULTS

Population structure, principal component analysis, and phylogenetic trees all confirmed the evolutionary relationship between maize and teosinte, and determined the evolutionary lineage of all species within teosinte. Shared haplotype analysis showed similar levels of ancestral alleles from Zea mays ssp. parviglumis and Zea mays ssp. mexicana in maize. Scans for selection signatures identified 394 domestication sweeps by comparing wild and cultivated maize and 360 adaptation sweeps by comparing tropical and temperate maize. Permutation tests revealed that the public association signals for flowering time were highly enriched in the domestication and adaptation sweeps. Genome-wide association study identified 125 loci significantly associated with flowering-time traits, ten of which identified candidate genes that have undergone selection during maize adaptation.

CONCLUSIONS

In this study, we characterized the history of maize domestication and adaptation at the population genomic level and identified hundreds of domestication and adaptation sweeps. This study extends the molecular mechanism of maize domestication and adaptation, and provides resources for basic research and genetic improvement in maize.

Natural variation of ZmCGT1 is responsible for isoorientin accumulation in maize silk

Maize (Zea mays L.) silk contains high levels of flavonoids and is widely used to promote human health. Isoorientin, a natural C-glycoside flavone abundant in maize silk, has attracted considerable attention due to its potential value. Although different classes of flavonoid have been well characterized in plants, the genes involved in the biosynthesis of isoorientin in maize are largely unknown. Here, we used targeted metabolic profiling of isoorientin on the silks in an association panel consisting of 294 maize inbred lines. We identified the gene ZmCGT1 by genome-wide association analysis. The ZmCGT1 protein was characterized as a 2-hydroxyflavanone C-glycosyltransferase that can C-glycosylate 2-hydroxyflavanone to form flavone-C-glycoside after dehydration. Moreover, ZmCGT1 overexpression increased isoorientin levels and RNA interference-mediated ZmCGT1 knockdown decreased accumulation of isoorientin in maize silk. Further, two nucleotide polymorphisms, A502C and A1022G, which led to amino acid changes I168L and E341G, respectively, were identified to be functional polymorphisms responsible for the natural variation in isoorientin levels. In summary, we identified the gene ZmCGT1, which plays an important role in isoorientin biosynthesis, providing insights into the genetic basis of the natural variation in isoorientin levels in maize silk. The identified favorable CG allele of ZmCGT1 may be further used for genetic improvement of nutritional quality in maize.

The arches and spandrels of maize domestication, adaptation, and improvement

People living in the Balsas River basin in southwest México domesticated maize from the bushy grass teosinte. Nine thousand years later, in 2021, Ms. Deb Haaland - a member of the Pueblo of Laguna tribe of New Mexico - wore a dress adorned with a cornstalk when she was sworn in as the Secretary of Interior of the United States of America. This choice of garment highlights the importance of the coevolution of maize and the farmers who, through careful selection over thousands of years, domesticated maize and adapted the physiology and shoot architecture of maize to fit local environments and growth habits. Some traits such as tillering were directly selected on (arches), and others such as tassel size are the by-products (spandrels) of maize evolution. Here, we review current knowledge of the underlying cellular, developmental, physiological, and metabolic processes that were selected by farmers and breeders, which have positioned maize as a top global staple crop.

Phenotypic Plasticity Contributes to Maize Adaptation and Heterosis

Plant phenotypic plasticity describes altered phenotypic performance of an individual when grown in different environments. Exploring genetic architecture underlying plant plasticity variation may help mitigate the detrimental effects of a rapidly changing climate on agriculture, but little research has been done in this area to date. In the present study, we established a population of 976 maize F₁ hybrids by crossing 488 diverse inbred lines with two elite testers. Genome-wide association study identified hundreds of quantitative trait loci associated with phenotypic plasticity variation across diverse F₁ hybrids, the majority of which contributed very little variance, in accordance with the polygenic nature of these traits. We identified several quantitative trait locus regions that may have been selected during the tropical-temperate adaptation process. We also observed heterosis in terms of phenotypic plasticity, in addition to the traditional genetic value differences measured between hybrid and inbred lines, and the pattern of which was affected by genetic background. Our results demonstrate a landscape of phenotypic plasticity in maize, which will aid in the understanding of its genetic architecture, its contribution to adaptation and heterosis, and how it may be exploited for future maize breeding in a rapidly changing environment.

High-depth resequencing of 312 accessions reveals the local adaptation of foxtail millet

Based on the high-density variation map, we identified genome-level evidence for local adaptation and demonstrated that Siprr37 with transposon insertion contributes to the fitness of foxtail millet in the northeastern ecoregion. Adaptation is a robust way through which plants are able to overcome environmental constraints. The mechanisms of adaptation in heterogeneous natural environments are largely unknown. Deciphering the genomic basis of local adaptation will contribute to further improvement in domesticated plants. To this end, we describe a high-depth (19.4 ×) haplotype map of 3.02 million single nucleotide polymorphisms in foxtail millet (Setaria italica) from whole-genome resequencing of 312 accessions. In the genome-wide scan, we identified a set of improvement signals (including the homologous gene of OsIPA1, a key gene controlling ideal plant architecture) related to the geographical adaptation to four ecoregions in China. In particular, based on the genome-wide association analysis results, we identified the contribution of a pseudo-response regulator gene, SiPRR37, to heading date adaptation in foxtail millet. We observed the expression changes of SiPRR37 resulted from a key Tc1-Mariner transposon insertion in the first intron. Positive selection analyses revealed that SiPRR37 mainly contributed to the adaptation of northeastern ecoregions. Taken together, foxtail millet adapted to the northeastern region by regulating the function of SiPRR37, which sheds lights on genome-level evidence for adaptive geographical divergence. Besides, our data provide a nearly complete catalog of genomic variation aiding the identification of functionally important variants.

Identifying quantitative trait loci for the general combining ability of yield-relevant traits in maize

Maize is the most important staple crop worldwide. Many of its agronomic traits present with a high level of heterosis. Combining ability was proposed to exploit the rule of heterosis, and general combining ability (GCA) is a crucial measure of parental performance. In this study, a recombinant inbred line population was used to construct testcross populations by crossing with four testers based on North Carolina design II. Six yield-relevant traits were investigated as phenotypic data. GCA effects were estimated for three scenarios based on the heterotic group and the number of tester lines. These estimates were then used to identify quantitative trait loci (QTL) and dissect genetic basis of GCA. A higher heritability of GCA was obtained for each trait. Thus, testing in early generation of breeding may effectively select candidate lines with relatively superior GCA performance. The GCA QTL detected in each scenario was slightly different according to the linkage mapping. Most of the GCA-relevant loci were simultaneously detected in all three datasets. Therefore, the genetic basis of GCA was nearly constant although discrepant inbred lines were appointed as testers. In addition, favorable alleles corresponding to GCA could be pyramided via marker-assisted selection and made available for maize hybrid breeding.

Dissecting the phenotypic components and genetic architecture of maize stem vascular bundles using high-throughput phenotypic analysis

High-throughput phenotyping is increasingly becoming an important tool for rapid advancement of genetic gain in breeding programmes. Manual phenotyping of vascular bundles is tedious and time-consuming, which lags behind the rapid development of functional genomics in maize. More robust and automated techniques of phenotyping vascular bundles traits at high-throughput are urgently needed for large crop populations. In this study, we developed a standard process for stem micro-CT data acquisition and an automatic CT image process pipeline to obtain vascular bundle traits of stems including geometry-related, morphology-related and distribution-related traits. Next, we analysed the phenotypic variation of stem vascular bundles between natural population subgroup (480 inbred lines) based on 48 comprehensively phenotypic information. Also, the first database for stem micro-phenotypes, MaizeSPD, was established, storing 554 pieces of basic information of maize inbred lines, 523 pieces of experimental information, 1008 pieces of CT scanning images and processed images, and 24 192 pieces of phenotypic data. Combined with genome-wide association studies (GWASs), a total of 1562 significant single nucleotide polymorphism (SNPs) were identified for 30 stem micro-phenotypic traits, and 84 unique genes of 20 traits such as VBNum, VBAvArea and PZVBDensity were detected. Candidate genes identified by GWAS mainly encode enzymes involved in cell wall metabolism, transcription factors, protein kinase and protein related to plant signal transduction and stress response. The results presented here will advance our knowledge about phenotypic trait components of stem vascular bundles and provide useful information for understanding the genetic controls of vascular bundle formation and development.

Using single-plant-omics in the field to link maize genes to functions and phenotypes

Most of our current knowledge on plant molecular biology is based on experiments in controlled laboratory environments. However, translating this knowledge from the laboratory to the field is often not straightforward, in part because field growth conditions are very different from laboratory conditions. Here, we test a new experimental design to unravel the molecular wiring of plants and study gene-phenotype relationships directly in the field. We molecularly profiled a set of individual maize plants of the same inbred background grown in the same field and used the resulting data to predict the phenotypes of individual plants and the function of maize genes. We show that the field transcriptomes of individual plants contain as much information on maize gene function as traditional laboratory-generated transcriptomes of pooled plant samples subject to controlled perturbations. Moreover, we show that field-generated transcriptome and metabolome data can be used to quantitatively predict individual plant phenotypes. Our results show that profiling individual plants in the field is a promising experimental design that could help narrow the lab-field gap.

Transcriptomic divergence between upland and lowland ecotypes contributes to rice adaptation to a drought-prone agroecosystem

INTRODUCTION

Transcriptomic divergence drives plant ecological adaptation. Upland rice is differentiated in drought tolerance from lowland rice during its adaptation to the drought-prone environment. They provide a good system to learn the evolution of drought tolerance in rice.

METHODS AND RESULTS

We estimate morphological differences between the two rice ecotypes under well-watered and drought conditions, as well as their genetic and transcriptomic divergences by the high-throughput sequencing. Upland rice possesses higher expression diversity than lowland rice does. Thousands of genes exhibit expression divergences between the two rice ecotypes, which contributes to their morphological differences in drought tolerance. These transcriptomic divergences contribute to drought adaptation of upland rice during its domestication. Mutations in transcriptional regulatory regions, which cause presence and absence of cis-elements, are the cause of expression divergence. About 15.3% transcriptionally selected genes also receive sequence-based selection in upland or lowland ecotype. Some highly differentiated genes promote the transcriptomic divergence between rice ecotypes via gene co-expression network. In addition, we also detected transcriptomic trade-offs between drought tolerance and productivity.

DISCUSSION

Many key genes, which promote transcriptomic adaptation to drought in upland rice, have great prospective in breeding water-saving and drought-resistant rice. Meanwhile, appropriate strategies are required in breeding to overcome the potential transcriptomic trade-off.

QTL mapping and GWAS for field kernel water content and kernel dehydration rate before physiological maturity in maize

Kernel water content (KWC) and kernel dehydration rate (KDR) are two main factors affecting maize seed quality and have a decisive influence on the mechanical harvest. It is of great importance to map and mine candidate genes related to KWCs and KDRs before physiological maturity in maize. 120 double-haploid (DH) lines constructed from Si287 with low KWC and JiA512 with high KWC were used as the mapping population. KWCs were measured every 5 days from 10 to 40 days after pollination, and KDRs were calculated. A total of 1702 SNP markers were used to construct a linkage map, with a total length of 1,309.02 cM and an average map distance of 0.77 cM. 10 quantitative trait loci (QTLs) and 27 quantitative trait nucleotides (QTNs) were detected by genome-wide composite interval mapping (GCIM) and multi-locus random-SNP-effect mixed linear model (mrMLM), respectively. One and two QTL hotspot regions were found on Chromosome 3 and 7, respectively. Analysis of the Gene Ontology showed that 2 GO terms of biological processes (BP) were significantly enriched (P ≤ 0.05) and 6 candidate genes were obtained. This study provides theoretical support for marker-assisted breeding of mechanical harvest variety in maize.

Metabolomics analysis reveals differences in evolution between maize and rice

Metabolites are the intermediate and final products of metabolism, which play essential roles in plant growth, evolution and adaptation to changing climates. However, it is unclear how evolution contributes to metabolic variation in plants. Here, we investigated the metabolomics data from leaf and seed tissues in maize and rice. Using principal components analysis based on leaf metabolites but not seed metabolites, metabolomics data could be clearly separated for rice Indica and Japonica accessions, while two maize subgroups, temperate and tropical, showed more visible admixture. Rice and maize seed exhibited significant interspecific differences in metabolic variation, while within rice, leaf and seed displayed similar metabolic variations. Among 10 metabolic categories, flavonoids had higher variation in maize than rice, indicating flavonoids are a key constituent of interspecific metabolic divergence. Interestingly, metabolic regulation was also found to be reshaped dramatically from positive to negative correlations, indicative of the differential evolutionary processes in maize and rice. Moreover, perhaps due to this divergence significantly more metabolic interactions were identified in rice than maize. Furthermore, in rice, the leaf was found to harbor much more intense metabolic interactions than the seed. Our result suggests that metabolomes are valuable for tracking evolutionary history, thereby complementing and extending genomic insights concerning which features are responsible for interspecific differentiation in maize and rice.

Genome-wide expression quantitative trait locus analysis in a recombinant inbred line population for trait dissection in peanut

The transcriptome connects genome to the gene function and ultimate phenome in biology. So far, transcriptomic approach was not used in peanut for performing trait mapping in bi-parental populations. In this research, we sequenced the whole transcriptome in immature seeds in a peanut recombinant inbred line (RIL) population and explored thoroughly the landscape of transcriptomic variations and its genetic basis. The comprehensive analysis identified total 49 691 genes in RIL population, of which 92 genes followed a paramutation-like expression pattern. Expression quantitative trait locus (eQTL) analysis identified 1207 local eQTLs and 15 837 distant eQTLs contributing to the whole-genome transcriptomic variation in peanut. There were 94 eQTL hot spot regions detected across the genome with the dominance of distant eQTL. By integrating transcriptomic profile and annotation analyses, we unveiled a putative candidate gene and developed a linked marker InDel02 underlying a major QTL responsible for purple testa colour in peanut. Our result provided a first understanding of genetic basis of whole-genome transcriptomic variation in peanut and illustrates the potential of the transcriptome-aid approach in dissecting important traits in non-model plants.

The Past, Present, and Future of Maize Improvement: Domestication, Genomics, and Functional Genomic Routes toward Crop Enhancement

After being domesticated from teosinte, cultivated maize (Zea mays ssp. mays) spread worldwide and now is one of the most important staple crops. Due to its tremendous phenotypic and genotypic diversity, maize also becomes to be one of the most widely used model plant species for fundamental research, with many important discoveries reported by maize researchers. Here, we provide an overview of the history of maize domestication and key genes controlling major domestication-related traits, review the currently available resources for functional genomics studies in maize, and discuss the functions of most of the maize genes that have been positionally cloned and can be used for crop improvement. Finally, we provide some perspectives on future directions regarding functional genomics research and the breeding of maize and other crops.

Divergent selection and drift shape the genomes of two avian sister species spanning a saline-freshwater ecotone

The role of species divergence due to ecologically based divergent selection-or ecological speciation-in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline-freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes-the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)-to evaluate the influence of selective and demographic processes in shaping genome-wide patterns of divergence. Genome-wide divergence between these two recently diverged sister species was notably high (genome-wide F ST = 0.32). Against a background of high genome-wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism-all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.

Investigation and development of maize fused network analysis with multi-omics

Maize is a critically important staple crop in the whole world, which has contributed to both economic security and food in planting areas. The main target for researchers and breeding is the improvement of maize quality and yield. The use of computational biology methods combined with multi-omics for selecting biomolecules of interest for maize breeding has been receiving more attention. Moreover, the rapid growth of high-throughput sequencing data provides the opportunity to explore biomolecules of interest at the molecular level in maize. Furthermore, we constructed weighted networks for each of the omics and then integrated them into a final fused weighted network based on a nonlinear combination method. We also analyzed the final fused network and mined the orphan nodes, some of which were shown to be transcription factors that played a key role in maize development. This study could help to improve maize production via insights at the multi-omics level and provide a new perspective for maize researchers. All related data have been released at http://lab.malab.cn/∼jj/maize.htm.

Large-scale metabolite quantitative trait locus analysis provides new insights for high-quality maize improvement

It is generally recognized that many favorable genes which were lost during domestication, including those related to both nutritional value and stress resistance, remain hidden in wild relatives. To uncover such genes in teosinte, an ancestor of maize, we conducted metabolite profiling in a BC₂ F₇ population generated from a cross between the maize wild relative (Zea mays ssp. mexicana) and maize inbred line Mo17. In total, 65 primary metabolites were quantified in four tissues (seedling-stage leaf, grouting-stage leaf, young kernel and mature kernel) with clear tissue-specific patterns emerging. Three hundred and fifty quantitative trait loci (QTLs) for these metabolites were obtained, which were distributed unevenly across the genome and included two QTL hotspots. Metabolite concentrations frequently increased in the presence of alleles from the teosinte genome while the opposite was observed for grain yield and shape trait QTLs. Combination of the multi-tissue transcriptome and metabolome data provided considerable insight into the metabolic variations between maize and its wild relatives. This study thus identifies favorable genes hidden in the wild relative which should allow us to balance high yield and quality in future modern crop breeding programs.

De Novo Domestication: An Alternative Route toward New Crops for the Future

Current global agricultural production must feed over 7 billion people. However, productivity varies greatly across the globe and is under threat from both increased competitions for land and climate change and associated environmental deterioration. Moreover, the increase in human population size and dietary changes are putting an ever greater burden on agriculture. The majority of this burden is met by the cultivation of a very small number of species, largely in locations that differ from their origin of domestication. Recent technological advances have raised the possibility of de novo domestication of wild plants as a viable solution for designing ideal crops while maintaining food security and a more sustainable low-input agriculture. Here we discuss how the discovery of multiple key domestication genes alongside the development of technologies for accurate manipulation of several target genes simultaneously renders de novo domestication a route toward crops for the future.

A domestication-associated reduction in K+ -preferring HKT transporter activity underlies maize shoot K+ accumulation and salt tolerance

Maize was domesticated from Balsas teosinte c. 10 000 yr ago. Previous studies have suggested that increased tolerance to environmental stress occurred during maize domestication. However, the underlying genetic basis remains largely unknown. We used a maize (W22)-teosinte recombinant inbred line (RIL) to investigate the salt wild-type tolerance aspects of maize domestication. We revealed that ZmHKT2 is a major QTL that regulates K⁺ homeostasis in saline soils. ZmHKT2 encodes a K⁺ -preferring HKT family transporter and probably reduces shoot K⁺ content by removing K⁺ ions from root-to-shoot flowing xylem sap, ZmHKT2 deficiency increases xylem sap and shoot K⁺ concentrations, and increases salt tolerance. A coding sequence polymorphism in the ZmHKT2^W22 allele (SNP389-G) confers an amino acid variant ZmHKT2 that increases xylem sap K⁺ concentration, thereby increasing shoot K⁺ content and salt tolerance. Additional analyses showed that SNP389-G first existed in teosinte (allele frequency 56% in assayed accessions), then swept through the maize population (allele frequency 98%), and that SNP389-G probably underwent positive selection during maize domestication. We conclude that a domestication-associated reduction in K⁺ transport activity in ZmHKT2 underlies maize shoot K⁺ content and salt tolerance, and propose that CRISPR-based editing of ZmHKT2 might provide a feasible strategy for improving maize salt tolerance.

Dynamic plant height QTL revealed in maize through remote sensing phenotyping using a high-throughput unmanned aerial vehicle (UAV)

Plant height (PH) is a key factor in maize (Zea mays L.) yield, biomass, and plant architecture. We investigated the PH of diverse maize inbred lines (117 temperate lines, 135 tropical lines) at four growth stages using unmanned aerial vehicle high-throughput phenotypic platforms (UAV-HTPPs). We extracted PH data using an automated pipeline based on crop surface models and orthomosaic model. The correlation between UAV and manually measured PH data reached 0.95. Under temperate field conditions, temperate maize lines grew faster than tropical maize lines at early growth stages, but tropical lines grew faster at later growth stages and ultimately became taller than temperate lines. A genome-wide association study identified 68 unique quantitative trait loci (QTLs) for seven PH-related traits, and 35% of the QTLs coincided with those previously reported to control PH. Generally, different QTLs controlled PH at different growth stages, but eight QTLs simultaneously controlled PH and growth rate at multiple growth stages. Based on gene annotations and expression profiles, we identified candidate genes controlling PH. The PH data collected by the UAV-HTPPs were credible and the genetic mapping power was high. Therefore, UAV-HTPPs have great potential for use in studies on PH.

Characterization of Proteome Variation During Modern Maize Breeding

The success of modern maize breeding has been demonstrated by remarkable increases in productivity with tremendous modification of agricultural phenotypes over the last century. Although the underlying genetic changes of the maize adaptation from tropical to temperate regions have been extensively studied, our knowledge is limited regarding the accordance of protein and mRNA expression levels accompanying such adaptation. Here we conducted an integrative analysis of proteomic and transcriptomic changes in a maize association panel. The minimum extent of correlation between protein and RNA levels suggests that variation in mRNA expression is often not indicative of protein expression at a population scale. This is corroborated by the observation that mRNA- and protein-based coexpression networks are relatively independent of each other, and many pQTLs arise without the presence of corresponding eQTLs. Importantly, compared with transcriptome, the subtypes categorized by the proteome show a markedly high accuracy to resemble the genomic subpopulation. These findings suggest that proteome evolved under a greater evolutionary constraint than transcriptome during maize adaptation from tropical to temperate regions. Overall, the integrated multi-omics analysis provides a functional context to interpret gene expression variation during modern maize breeding.

Genome-scale identification, classification, and tissue specific expression analysis of late embryogenesis abundant (LEA) genes under abiotic stress conditions in Sorghum bicolor L

Late embryogenesis abundant (LEA) proteins, the space fillers or molecular shields, are the hydrophilic protective proteins which play an important role during plant development and abiotic stress. The systematic survey and characterization revealed a total of 68 LEA genes, belonging to 8 families in Sorghum bicolor. The LEA-2, a typical hydrophobic family is the most abundant family. All of them are evenly distributed on all 10 chromosomes and chromosomes 1, 2, and 3 appear to be the hot spots. Majority of the S. bicolor LEA (SbLEA) genes are intron less or have fewer introns. A total of 22 paralogous events were observed and majority of them appear to be segmental duplications. Segmental duplication played an important role in SbLEA-2 family expansion. A total of 12 orthologs were observed with Arabidopsis and 13 with Oryza sativa. Majority of them are basic in nature, and targeted by chloroplast subcellular localization. Fifteen miRNAs targeted to 25 SbLEAs appear to participate in development, as well as in abiotic stress tolerance. Promoter analysis revealed the presence of abiotic stress-responsive DRE, MYB, MYC, and GT1, biotic stress-responsive W-Box, hormone-responsive ABA, ERE, and TGA, and development-responsive SKn cis-elements. This reveals that LEA proteins play a vital role during stress tolerance and developmental processes. Using microarray data, 65 SbLEA genes were analyzed in different tissues (roots, pith, rind, internode, shoot, and leaf) which show clear tissue specific expression. qRT-PCR analysis of 23 SbLEA genes revealed their abundant expression in various tissues like roots, stems and leaves. Higher expression was noticed in stems compared to roots and leaves. Majority of the SbLEA family members were up-regulated at least in one tissue under different stress conditions. The SbLEA3-2 is the regulator, which showed abundant expression under diverse stress conditions. Present study provides new insights into the formation of LEAs in S. bicolor and to understand their role in developmental processes under stress conditions, which may be a valuable source for future research.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Exploiting SPL genes to improve maize plant architecture tailored for high-density planting

Maize (Zea mays ssp. mays) is an agronomically important crop and also a classical genetic model for studying the regulation of plant architecture formation, which is a critical determinant of grain yield. Since the 1930s, increasing planting density has been a major contributing factor to the >7-fold increase in maize grain yield per unit land area in the USA, which is accompanied by breeding and utilization of cultivars characterized by high-density-tolerant plant architecture, including decreased ear height, lodging resistance, more upright leaves, reduced tassel branch number, and reduced anthesis-silking interval (ASI). Recent studies demonstrated that phytochrome-mediated red/far-red light signaling pathway and the miR156/SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) regulatory module co-ordinately regulate the shade avoidance response and diverse aspects of plant architecture in responding to shading in Arabidopsis. The maize genome contains 30 ZmSPL genes, and 18 of them are predicted as direct targets of zma-miR156s. Accumulating evidence indicates that ZmSPL genes play important roles in regulating maize flowering time, plant/ear height, tilling, leaf angle, tassel and ear architecture, and grain size and shape. Finally, we discuss ways to exploit maize SPL genes and downstream targets for improving maize plant architecture tailored for high-density planting.

ZmCOL3, a CCT gene represses flowering in maize by interfering with the circadian clock and activating expression of ZmCCT

Flowering time is a trait vital to the adaptation of flowering plants to different environments. Here, we report that CCT domain genes play an important role in flowering in maize (Zea mays L.). Among the 53 CCT family genes we identified in maize, 28 were located in flowering time quantitative trait locus regions and 15 were significantly associated with flowering time, based on candidate-gene association mapping analysis. Furthermore, a CCT gene named ZmCOL3 was shown to be a repressor of flowering. Overexpressing ZmCOL3 delayed flowering time by approximately 4 d, in either long-day or short-day conditions. The absence of one cytosine in the ZmCOL3 3'UTR and the presence of a 551 bp fragment in the promoter region are likely the causal polymorphisms contributing to the maize adaptation from tropical to temperate regions. We propose a modified model of the maize photoperiod pathway, wherein ZmCOL3 acts as an inhibitor of flowering either by transactivating transcription of ZmCCT, one of the key genes regulating maize flowering, or by interfering with the circadian clock.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

Background

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

Results

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

Conclusion

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0596-1) contains supplementary material, which is available to authorized users.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

BACKGROUND

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

RESULTS

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

CONCLUSION

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

RNA sequencing provides insights into the evolution of lettuce and the regulation of flavonoid biosynthesis

Different horticultural types of lettuce exhibit tremendous morphological variation. However, the molecular basis for domestication and divergence among the different horticultural types of lettuce remains unknown. Here, we report the RNA sequencing of 240 lettuce accessions sampled from the major horticultural types and wild relatives, generating 1.1 million single-nucleotide polymorphisms (SNPs). Demographic modeling indicates that there was a single domestication event for lettuce. We identify a list of regions as putative selective sweeps that occurred during domestication and divergence, respectively. Genome-wide association studies (GWAS) identify 5311 expression quantitative trait loci (eQTL) regulating the expression of 4105 genes, including nine eQTLs regulating genes associated with flavonoid biosynthesis. GWAS for leaf color detects six candidate loci responsible for the variation of anthocyanins in lettuce leaves. Our study provides a comprehensive understanding of the domestication and the accumulation of anthocyanins in lettuce, and will facilitate the breeding of cultivars with improved nutritional value.

Genome-wide characterization of non-reference transposable element insertion polymorphisms reveals genetic diversity in tropical and temperate maize

Background

Maize was originally domesticated in a tropical environment but is now widely cultivated at temperate latitudes. Temperate and tropical maize populations have diverged both genotypically and phenotypically. Tropical maize lines grown in temperate environments usually exhibit delayed flowering, pollination, and seed set, which reduces their grain yield relative to temperate adapted maize lines. One potential mechanism by which temperate maize may have adapted to a new environment is novel transposable element insertions, which can influence gene regulation. Recent advances in sequencing technology have made it possible to study variation in transposon content and insertion location in large sets of maize lines.

Results

In total, 274,408 non-redundant TEs (NRTEs) were identified using resequencing data generated from 83 maize inbred lines. The locations of DNA TEs and copia-superfamily retrotransposons showed significant positive correlations with gene density and genetic recombination rates, whereas gypsy-superfamily retrotransposons showed a negative correlation with these two parameters. Compared to tropical maize, temperate maize had fewer unique NRTEs but higher insertion frequency, lower background recombination rates, and higher linkage disequilibrium, with more NRTEs close to flowering and stress-related genes in the genome. Association mapping demonstrated that the presence/absence of 48 NRTEs was associated with flowering time and that expression of neighboring genes differed between haplotypes where a NRTE was present or absent.

Conclusions

This study suggests that NRTEs may have played an important role in creating the variation in gene regulation that enabled the rapid adaptation of maize to diverse environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4103-x) contains supplementary material, which is available to authorized users.

Fast diffusion of domesticated maize to temperate zones

Adaptation to a temperate climate was a prerequisite for the spread of maize across a broad geographical range. To explicitly explore the demographic process underlying maize adaptation, we used a diffusion-based method to model the differentiation between temperate and tropical populations using the Non-Stiff Stalk group as a proxy for temperate maize. Based on multiple sequential Markovian coalescent approaches, we estimate that tropical and temperate maize diverged approximately 3'000 to 5'000 years ago and the population size shrank after the split. Using composite likelihood approaches, we identified a distinct tropical-temperate divergence event initiated 4'958 years ago (95% confidence interval (CI): 4'877-5'039) from an ancestral population whose effective size was 24,162 (95% CI: 23,914-24,409). We found that continuous gene flow between tropical and temperate maize accompanied the differentiation of temperate maize. Long identical-by-descent tracts shared by tropical and temperate inbred lines have been identified, which might be the result of gene flow between tropical and temperate maize or artificial selection during domestication and crop improvement. Understanding the demographic history of maize diffusion not only provides evidence for population dynamics of maize, but will also assist the identification of regions under selection and the genetic basis of complex traits of agronomic importance.

The genomic basis of adaptation in plants

Plants are powerful models for the study of adaptive evolution. Since they are rooted in place, they must directly face environmental insults, making adaptation to local conditions vital. In addition to adaptation to natural conditions, some plant species have held a central role in human subsistence over the past several thousand years. In these species, humans exerted strong selective pressures on traits of agricultural importance. Recently, an increasing number of studies have aimed to identify the genomic basis of adaptation. These studies have provided insights into the mechanisms through which the raw materials of adaptation were introduced as well as the modes of adaptation in wild and domesticated species.

Distant eQTLs and Non-coding Sequences Play Critical Roles in Regulating Gene Expression and Quantitative Trait Variation in Maize

A detailed understanding of genetic architecture of mRNA expression by millions of genetic variants is important for studying quantitative trait variation. In this study, we identified 1.25M SNPs with a minor allele frequency greater than 0.05 by combining reduced genome sequencing (GBS), high-density array technologies (600K), and previous deep RNA-sequencing data from 368 diverse inbred lines of maize. The balanced allelic frequencies and distributions in a relatively large and diverse natural panel helped to identify expression quantitative trait loci (eQTLs) associated with more than 18 000 genes (63.4% of tested genes). We found that distant eQTLs were more frequent (∼75% of all eQTLs) across the whole genome. Thirteen novel associated loci affecting maize kernel oil concentration were identified using the new dataset, among which one intergenic locus affected the kernel oil variation by controlling expression of three other known oil-related genes. Altogether, this study provides resources for expanding our understanding of cellular regulatory mechanisms of transcriptome variation and the landscape of functional variants within the maize genome, thereby enhancing the understanding of quantitative variations.

Genome-wide Association Studies in Maize: Praise and Stargaze

Genome-wide association study (GWAS) has become a widely accepted strategy for decoding genotype-phenotype associations in many species thanks to advances in next-generation sequencing (NGS) technologies. Maize is an ideal crop for GWAS and significant progress has been made in the last decade. This review summarizes current GWAS efforts in maize functional genomics research and discusses future prospects in the omics era. The general goal of GWAS is to link genotypic variations to corresponding differences in phenotype using the most appropriate statistical model in a given population. The current review also presents perspectives for optimizing GWAS design and analysis. GWAS analysis of data from RNA, protein, and metabolite-based omics studies is discussed, along with new models and new population designs that will identify causes of phenotypic variation that have been hidden to date. The joint and continuous efforts of the whole community will enhance our understanding of maize quantitative traits and boost crop molecular breeding designs.

Joint genome-wide association and transcriptome sequencing reveals a complex polygenic network underlying hypocotyl elongation in rapeseed (Brassica napus L.)

Hypocotyl elongation is considered an important typical seedling trait contributing directly to an increase in and stabilization of the yield in Brassica napus, but its molecular genetic mechanism is poorly understood. In the present study, hypocotyl lengths of 210 lines were measured in an illuminated culture room. A genome-wide association study (GWAS) was performed with 23,435 single nucleotide polymorphisms (SNPs) for hypocotyl length. Three lines with long hypocotyl length and three lines with short hypocotyl length from one doubled haploid line (DH) population were used for transcriptome sequencing. A GWAS followed by transcriptome analysis identified 29 differentially expressed genes associated with significant SNPs in B. napus. These genes regulate hypocotyl elongation by mediating flowering morphogenesis, circadian clock, hormone biosynthesis, or important metabolic signaling pathways. Among these genes, BnaC07g46770D negatively regulates hypocotyl elongation directly, as well as flowering time. Our results indicate that a joint GWAS and transcriptome analysis has significant potential for identifying the genes responsible for hypocotyl elongation; The extension of hypocotyl is a complex biological process regulated by a polygenic network.

A Gene-Oriented Haplotype Comparison Reveals Recently Selected Genomic Regions in Temperate and Tropical Maize Germplasm

The extensive genetic variation present in maize (Zea mays) germplasm makes it possible to detect signatures of positive artificial selection that occurred during temperate and tropical maize improvement. Here we report an analysis of 532,815 polymorphisms from a maize association panel consisting of 368 diverse temperate and tropical inbred lines. We developed a gene-oriented approach adapting exonic polymorphisms to identify recently selected alleles by comparing haplotypes across the maize genome. This analysis revealed evidence of selection for more than 1100 genomic regions during recent improvement, and included regulatory genes and key genes with visible mutant phenotypes. We find that selected candidate target genes in temperate maize are enriched in biosynthetic processes, and further examination of these candidates highlights two cases, sucrose flux and oil storage, in which multiple genes in a common pathway can be cooperatively selected. Finally, based on available parallel gene expression data, we hypothesize that some genes were selected for regulatory variations, resulting in altered gene expression.

Exploring Identity-By-Descent Segments and Putative Functions Using Different Foundation Parents in Maize

Maize foundation parents (FPs) play no-alternative roles in hybrid breeding because they were widely used in the development of new lines and hybrids. The combination of different identity-by-descent (IBD) segments and genes could account for the formation patterns of different FPs, and knowledge of these IBD regions would provide an extensive foundation for the development of new candidate FP lines in future maize breeding. In this paper, a panel of 304 elite lines derived from FPs, i.e., B73, 207, Mo17, and Huangzaosi (HZS), was collected and analyzed using 43,252 single nucleotide polymorphism (SNP) markers. Most IBD segments specific to particular FP groups were identified, including 116 IBD segments in B73, 105 in Mo17, 111 in 207, and 190 in HZS. In these regions, 423 quantitative trait nucleotides (QTNs) associated with 15 agronomic traits and 804 candidate genes were identified. Some known adaptation-related genes, e.g., dwarf8 and vgt1 in HZS, zcn8 and epc in Mo17, and ZmCCT in 207, were validated as being tightly linked to particular IBD segments. In addition, numerous new candidate genes were also identified. For example, GRMZM2G154278 in HZS, which belongs to the cell cycle control family, was closely linked to a QTN of the ear height/plant height (EH/PH) trait; GRMZM2G051943 in 207, which encodes an endochitinase precursor (EP) chitinase, was closely linked to a QTN for kernel density; and GRMZM2G170586 in Mo17 was closely linked to a QTN for ear diameter. Complex correlations among these genes were also found. Many IBD segments and genes were included in the formation of FP lines, and complex regulatory networks exist among them. These results provide new insights on the genetic basis of complex traits and provide new candidate IBD regions or genes for the improvement of special traits in maize production.

Characterizing the population structure and genetic diversity of maize breeding germplasm in Southwest China using genome-wide SNP markers

Background

Maize breeding germplasm used in Southwest China has high complexity because of the diverse ecological features of this area. In this study, the population structure, genetic diversity, and linkage disequilibrium decay distance of 362 important inbred lines collected from the breeding program of Southwest China were characterized using the MaizeSNP50 BeadChip with 56,110 single nucleotide polymorphisms (SNPs).

Results

With respect to population structure, two (Tropical and Temperate), three (Tropical, Stiff Stalk and non-Stiff Stalk), four [Tropical, group A germplasm derived from modern U.S. hybrids (PA), group B germplasm derived from modern U.S. hybrids (PB) and Reid] and six (Tropical, PB, Reid, Iowa Stiff Stalk Synthetic, PA and North) subgroups were identified. With increasing K value, the Temperate group showed pronounced hierarchical structure with division into further subgroups. The Genetic Diversity of each group was also estimated, and the Tropical group was more diverse than the Temperate group. Seven low-genetic-diversity and one high-genetic-diversity regions were collectively identified in the Temperate, Tropical groups, and the entire panel. SNPs with significant variation in allele frequency between the Tropical and Temperate groups were also evaluated. Among them, a region located at 130 Mb on Chromosome 2 showed the highest genetic diversity, including both number of SNPs with significant variation and the ratio of significant SNPs to total SNPs. Linkage disequilibrium decay distance in the Temperate group was greater (2.5–3 Mb) than that in the entire panel (0.5–0.75 Mb) and the Tropical group (0.25–0.5 Mb). A large region at 30–120 Mb of Chromosome 7 was concluded to be a region conserved during the breeding process by comparison between S37, which was considered a representative tropical line in Southwest China, and its 30 most similar derived lines.

Conclusions

For the panel covered most of widely used inbred lines in Southwest China, this work representatively not only illustrates the foundation and evolution trend of maize breeding resource as a theoretical reference for the improvement of heterosis, but also provides plenty of information for genetic researches such as genome-wide association study and marker-assisted selection in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3041-3) contains supplementary material, which is available to authorized users.

Genome-wide identification, structural analysis and new insights into late embryogenesis abundant (LEA) gene family formation pattern in Brassica napus

Late embryogenesis abundant (LEA) proteins are a diverse and large group of polypeptides that play important roles in desiccation and freezing tolerance in plants. The LEA family has been systematically characterized in some plants but not Brassica napus. In this study, 108 BnLEA genes were identified in the B. napus genome and classified into eight families based on their conserved domains. Protein sequence alignments revealed an abundance of alanine, lysine and glutamic acid residues in BnLEA proteins. The BnLEA gene structure has few introns (<3), and they are distributed unevenly across all 19 chromosomes in B. napus, occurring as gene clusters in chromosomes A9, C2, C4 and C5. More than two-thirds of the BnLEA genes are associated with segmental duplication. Synteny analysis revealed that most LEA genes are conserved, although gene losses or gains were also identified. These results suggest that segmental duplication and whole-genome duplication played a major role in the expansion of the BnLEA gene family. Expression profiles analysis indicated that expression of most BnLEAs was increased in leaves and late stage seeds. This study presents a comprehensive overview of the LEA gene family in B. napus and provides new insights into the formation of this family.