Genetic variation in ZmVPP1 contributes to drought tolerance in maize seedlings

Abiotic stress responses in plants

Plants cannot move, so they must endure abiotic stresses such as drought, salinity and extreme temperatures. These stressors greatly limit the distribution of plants, alter their growth and development, and reduce crop productivity. Recent progress in our understanding of the molecular mechanisms underlying the responses of plants to abiotic stresses emphasizes their multilevel nature; multiple processes are involved, including sensing, signalling, transcription, transcript processing, translation and post-translational protein modifications. This improved knowledge can be used to boost crop productivity and agricultural sustainability through genetic, chemical and microbial approaches.

Exploitation of Drought Tolerance-Related Genes for Crop Improvement

Drought has become a major threat to food security, because it affects crop growth and development. Drought tolerance is an important quantitative trait, which is regulated by hundreds of genes in crop plants. In recent decades, scientists have made considerable progress to uncover the genetic and molecular mechanisms of drought tolerance, especially in model plants. This review summarizes the evaluation criteria for drought tolerance, methods for gene mining, characterization of genes related to drought tolerance, and explores the approaches to enhance crop drought tolerance. Collectively, this review illustrates the application prospect of these genes in improving the drought tolerance breeding of crop plants.

Genome-wide association study reveals the genetic architecture of root hair length in maize

Background

Root hair, a special type of tubular-shaped cell, outgrows from root epidermal cell and plays important roles in the acquisition of nutrients and water, as well as interactions with biotic and abiotic stress. Although many genes involved in root hair development have been identified, genetic basis of natural variation in root hair growth has never been explored.

Results

Here, we utilized a maize association panel including 281 inbred lines with tropical, subtropical, and temperate origins to decipher the phenotypic diversity and genetic basis of root hair length. We demonstrated significant associations of root hair length with many metabolic pathways and other agronomic traits. Combining root hair phenotypes with 1.25 million single nucleotide polymorphisms (SNPs) via genome-wide association study (GWAS) revealed several candidate genes implicated in cellular signaling, polar growth, disease resistance and various metabolic pathways.

Conclusions

These results illustrate the genetic basis of root hair length in maize, offering a list of candidate genes predictably contributing to root hair growth, which are invaluable resource for the future functional investigation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07961-z.

Genomic basis underlying the metabolome-mediated drought adaptation of maize

Background

Drought is a major environmental disaster that causes crop yield loss worldwide. Metabolites are involved in various environmental stress responses of plants. However, the genetic control of metabolomes underlying crop environmental stress adaptation remains elusive.

Results

Here, we perform non-targeted metabolic profiling of leaves for 385 maize natural inbred lines grown under well-watered as well as drought-stressed conditions. A total of 3890 metabolites are identified and 1035 of these are differentially produced between well-watered and drought-stressed conditions, representing effective indicators of maize drought response and tolerance. Genetic dissections reveal the associations between these metabolites and thousands of single-nucleotide polymorphisms (SNPs), which represented 3415 metabolite quantitative trait loci (mQTLs) and 2589 candidate genes. 78.6% of mQTLs (2684/3415) are novel drought-responsive QTLs. The regulatory variants that control the expression of the candidate genes are revealed by expression QTL (eQTL) analysis of the transcriptomes of leaves from 197 maize natural inbred lines. Integrated metabolic and transcriptomic assays identify dozens of environment-specific hub genes and their gene-metabolite regulatory networks. Comprehensive genetic and molecular studies reveal the roles and mechanisms of two hub genes, Bx12 and ZmGLK44, in regulating maize metabolite biosynthesis and drought tolerance.

Conclusion

Our studies reveal the first population-level metabolomes in crop drought response and uncover the natural variations and genetic control of these metabolomes underlying crop drought adaptation, demonstrating that multi-omics is a powerful strategy to dissect the genetic mechanisms of crop complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02481-1.

Dissection of the Genetic Basis of Rice Panicle Architecture Using a Genome-wide Association Study

Panicle architecture is one of the major factors influencing productivity of rice crops. The regulatory mechanisms underlying this complex trait are still unclear and genetic resources for rice breeders to improve panicle architecture are limited. Here, we have performed a genome-wide association study (GWAS) to analyze and identify genetic determinants underlying three panicle architecture traits. A population of 340 rice accessions from the 3000 Rice Genomes Project was phenotyped for panicle length, primary panicle number and secondary branch number over two years; GWAS was performed across the whole panel, and also across the japonica and indica sub-panels. A total of 153 quantitative trait loci (QTLs) were detected, of which 5 were associated with multiple traits, 8 were unique to either indica or japonica sub-panels, while 37 QTLs were stable across both years. Using haplotype and expression analysis, we reveal that genetic variations in the OsSPL18 promoter significantly affect gene expression and correlate with panicle length phenotypes. Three new candidate genes with putative roles in determining panicle length were also identified. Haplotype analysis of OsGRRP and LOC_Os03g03480 revealed high association with panicle length variation. Gene expression of DSM2, involved in abscisic acid biosynthesis, was up-regulated in long panicle accessions. Our results provide valuable information and resources for further unravelling the genetic basis determining rice panicle architecture. Identified candidate genes and molecular markers can be used in marker-assisted selection to improve rice panicle architecture through molecular breeding.

Breeding custom-designed crops for improved drought adaptation

The current pace of crop improvement is inadequate to feed the burgeoning human population by 2050. Higher, more stable, and sustainable crop production is required against a backdrop of drought stress, which causes significant losses in crop yields. Tailoring crops for drought adaptation may hold the key to address these challenges and provide resilient production systems for future harvests. Understanding the genetic and molecular landscape of the functionality of alleles associated with adaptive traits will make designer crop breeding the prospective approach for crop improvement. Here, we highlight the potential of genomics technologies combined with crop physiology for high-throughput identification of the genetic architecture of key drought-adaptive traits and explore innovative genomic breeding strategies for designing future crops.

Genetic Dissection of Phosphorus Use Efficiency in a Maize Association Population under Two P Levels in the Field

Phosphorus (P) deficiency is an important challenge the world faces while having to increase crop yields. It is therefore necessary to select maize (Zea may L.) genotypes with high phosphorus use efficiency (PUE). Here, we extensively analyzed the biomass, grain yield, and PUE-related traits of 359 maize inbred lines grown under both low-P and normal-P conditions. A significant decrease in grain yield per plant and biomass, an increase in PUE under low-P condition, as well as significant correlations between the two treatments were observed. In a genome-wide association study, 49, 53, and 48 candidate genes were identified for eleven traits under low-P, normal-P conditions, and in low-P tolerance index (phenotype under low-P divided by phenotype under normal-P condition) datasets, respectively. Several gene ontology pathways were enriched for the genes identified under low-P condition. In addition, seven key genes related to phosphate transporter or stress response were molecularly characterized. Further analyses uncovered the favorable haplotype for several core genes, which is less prevalent in modern lines but often enriched in a specific subpopulation. Collectively, our research provides progress in the genetic dissection and molecular characterization of PUE in maize.

Genetics of Germination and Seedling Traits under Drought Stress in a MAGIC Population of Maize

Drought is one of the most detrimental abiotic stresses hampering seed germination, development, and productivity. Maize is more sensitive to drought than other cereals, especially at seedling stage. Our objective was to study genetic regulation of drought tolerance at germination and during seedling growth in maize. We evaluated 420 RIL with their parents from a multi-parent advanced generation inter-cross (MAGIC) population with PEG-induced drought at germination and seedling establishment. A genome-wide association study (GWAS) was carried out to identify genomic regions associated with drought tolerance. GWAS identified 28 and 16 SNPs significantly associated with germination and seedling traits under stress and well-watered conditions, respectively. Among the SNPs detected, two SNPs had significant associations with several traits with high positive correlations, suggesting a pleiotropic genetic control. Other SNPs were located in regions that harbored major QTLs in previous studies, and co-located with QTLs for cold tolerance previously published for this MAGIC population. The genomic regions comprised several candidate genes related to stresses and plant development. These included numerous drought-responsive genes and transcription factors implicated in germination, seedling traits, and drought tolerance. The current analyses provide information and tools for subsequent studies and breeding programs for improving drought tolerance.

Favorable pleiotropic loci for fiber yield and quality in upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum L.) is an important economic crop for renewable textile fibers. However, the simultaneous improvement of yield and fiber quality in cotton is difficult as the linkage drag. Compared with breaking the linkage drag, identification of the favorable pleiotropic loci on the genome level by genome-wide association study (GWAS) provides a new way to improve the yield and fiber quality simultaneously. In our study restriction-site-associated DNA sequencing (RAD-seq) was used to genotype 316 cotton accessions. Eight major traits in three categories including yield, fiber quality and maturation were investigated in nine environments (3 sites × 3 years). 231 SNPs associated with these eight traits (− log₁₀(P) > 5.27) were identified, located in 27 genomic regions respectively by linkage disequilibrium analysis. Further analysis showed that four genomic regions (the region 1, 6, 8 and 23) held favorable pleiotropic loci and 6 candidate genes were identified. Through genotyping, 14 elite accessions carrying the favorable loci on four pleiotropic regions were identified. These favorable pleiotropic loci and elite genotypes identified in this study will be utilized to improve the yield and fiber quality simultaneously in future cotton breeding.

A forward genetics approach integrating genome-wide association study and expression quantitative trait locus mapping to dissect leaf development in maize (Zea mays)

The characterization of the genetic basis of maize (Zea mays) leaf development may support breeding efforts to obtain plants with higher vigor and productivity. In this study, a mapping panel of 197 biparental and multiparental maize recombinant inbred lines (RILs) was analyzed for multiple leaf traits at the seedling stage. RNA sequencing was used to estimate the transcription levels of 29 573 gene models in RILs and to derive 373 769 single nucleotide polymorphisms (SNPs), and a forward genetics approach combining these data was used to pinpoint candidate genes involved in leaf development. First, leaf traits were correlated with gene expression levels to identify transcript-trait correlations. Then, leaf traits were associated with SNPs in a genome-wide association (GWA) study. An expression quantitative trait locus mapping approach was followed to associate SNPs with gene expression levels, prioritizing candidate genes identified based on transcript-trait correlations and GWAs. Finally, a network analysis was conducted to cluster all transcripts in 38 co-expression modules. By integrating forward genetics approaches, we identified 25 candidate genes highly enriched for specific functional categories, providing evidence supporting the role of vacuolar proton pumps, cell wall effectors, and vesicular traffic controllers in leaf growth. These results tackle the complexity of leaf trait determination and may support precision breeding in maize.

Genome-wide association studies reveal the coordinated regulatory networks underlying photosynthesis and wood formation in Populus

Photosynthesis and wood formation underlie the ability of trees to provide renewable resources and perform ecological functions; however, the genetic basis and regulatory pathways coordinating these two linked processes remain unclear. Here, we used a systems genetics strategy, integrating genome-wide association studies, transcriptomic analyses, and transgenic experiments, to investigate the genetic architecture of photosynthesis and wood properties among 435 unrelated individuals of Populus tomentosa, and unravel the coordinated regulatory networks resulting in two trait categories. We detected 222 significant single-nucleotide polymorphisms, annotated to 177 candidate genes, for 10 traits of photosynthesis and wood properties. Epistasis uncovered 74 epistatic interactions for phenotypes. Strikingly, we deciphered the coordinated regulation patterns of pleiotropic genes underlying phenotypic variations for two trait categories. Furthermore, expression quantitative trait nucleotide mapping and coexpression analysis were integrated to unravel the potential transcriptional regulatory networks of candidate genes coordinating photosynthesis and wood properties. Finally, heterologous expression of two pleiotropic genes, PtoMYB62 and PtoMYB80, in Arabidopsis thaliana demonstrated that they control regulatory networks balancing photosynthesis and stem secondary cell wall components, respectively. Our study provides insights into the regulatory mechanisms coordinating photosynthesis and wood formation in poplar, and should facilitate genetic breeding in trees via molecular design.

Genetic basis of maize kernel oil-related traits revealed by high-density SNP markers in a recombinant inbred line population

BACKGROUND

Maize (Zea mays ssp. mays) is the most abundantly cultivated and highly valued food commodity in the world. Oil from maize kernels is highly nutritious and important for the diet and health of humans, and it can be used as a source of bioenergy. A better understanding of genetic basis for maize kernel oil can help improve the oil content and quality when applied in breeding.

RESULTS

In this study, a KUI3/SC55 recombinant inbred line (RIL) population, consisting of 180 individuals was constructed from a cross between inbred lines KUI3 and SC55. We phenotyped 19 oil-related traits and subsequently dissected the genetic architecture of oil-related traits in maize kernels based on a high-density genetic map. In total, 62 quantitative trait loci (QTLs), with 2 to 5 QTLs per trait, were detected in the KUI3/SC55 RIL population. Each QTL accounted for 6.7% (qSTOL1) to 31.02% (qBELI6) of phenotypic variation and the total phenotypic variation explained (PVE) of all detected QTLs for each trait ranged from 12.5% (OIL) to 52.5% (C16:0/C16:1). Of all these identified QTLs, only 5 were major QTLs located in three genomic regions on chromosome 6 and 9. In addition, two pairs of epistatic QTLs with additive effects were detected and they explained 3.3 and 2.4% of the phenotypic variation, respectively. Colocalization with a previous GWAS on oil-related traits, identified 19 genes. Of these genes, two important candidate genes, GRMZM2G101515 and GRMZM2G022558, were further verified to be associated with C20:0/C22:0 and C18:0/C20:0, respectively, according to a gene-based association analysis. The first gene encodes a kinase-related protein with unknown function, while the second gene encodes fatty acid elongase 2 (fae2) and directly participates in the biosynthesis of very long chain fatty acids in Arabidopsis.

CONCLUSIONS

Our results provide insights on the genetic basis of oil-related traits and a theoretical basis for improving maize quality by marker-assisted selection.

Manipulating ZmEXPA4 expression ameliorates the drought-induced prolonged anthesis and silking interval in maize

Drought poses a major environmental threat to maize (Zea mays) production worldwide. Since maize is a monoecious plant, maize grain yield is dependent on the synchronous development of male and female inflorescences. When a drought episode occurs during flowering, however, an asynchronism occurs in the anthesis and silking interval (ASI) that results in significant yield losses. The underlying mechanism responsible for this asynchronism is still unclear. Here, we obtained a comprehensive development-drought transcriptome atlas of maize ears. Genes that function in cell expansion and growth were highly repressed by drought in 50 mm ears. Notably, an association study using a natural-variation population of maize revealed a significant relationship between the level of α-expansin4 (ZmEXPA4) expression and drought-induced increases in ASI. Furthermore, genetic manipulation of ZmEXPA4 expression using a drought-inducible promoter in developing maize ears reduced the ASI under drought conditions. These findings provide important insights into the molecular mechanism underlying the increase in ASI in maize ears subjected to drought and provide a promising strategy that can be used for trait improvement.

Genome-Wide Association Study of Vascular Bundle-Related Traits in Maize Stalk

The vascular bundle plays an important role in nutrient transportation in plants and exerts great influence on crop yield. Maize is widely used for food, feed, and fuel, producing the largest yield in the world. However, genes and molecular mechanism controlling vascular bundle-related traits in maize have largely remained undiscovered. In this study, a natural population containing 248 diverse maize inbred lines genotyped with high-throughput SNP markers was used for genome-wide association study. The results showed that broad variations existed for the vascular bundle-related traits which are subject to genetic structure and it was suitable for association analysis. In this study, we identified 15, 13, 2, 1, and 5 SNPs significantly associated with number of small vascular bundle, number of large vascular bundle, average area of single small vascular bundle, average area of single large vascular bundle, and cross-sectional area, respectively. The 210 candidate genes in the confidence interval can be classified into ten biological processes, three cellular components, and eight molecular functions. As for the Kyoto Encyclopedia of Genes and Genomes analysis of the candidate genes, a total of six pathways were identified. Finally, we found five genes related to vascular development, three genes related to cell wall, and two genes related to the mechanical strength of the stalk. Our results provide the further understanding of the genetic foundation of vascular bundle-related traits in maize stalk.

Natural variations in the non-coding region of ZmNAC080308 contributes maintaining grain yield under drought stress in maize

BACKGROUND

Natural variations derived from both evolutionary selection and genetic recombination, presume to have important functions to respond to various abiotic stresses, which could be used to improve drought tolerance via genomic selection.

RESULTS

In the present study, the NAC-encoding gene of ZmNAC080308 was cloned and sequenced in 199 inbred lines in maize. Phylogenetic analysis showed that ZmNAC080308 is closely clusteredinto the same group with other well-known NAC genes responding to improve drought tolerance. In total, 86 SNPs and 47 InDels were identified in the generic region of ZmNAC080308, 19 of these variations were associated with GY (grain yield) in different environments. Nine variations in the 5'-UTR region of ZmNAC080308 are closely linked, they might regulate the gene expression and respond to improve GY under drought condition via Sp1-mediated transactivation. Two haplotypes (Hap1 and Hap2) identified in the, 5'-UTR region using the nine variations, and Hap2 containing insertion variants, exhibited 15.47 % higher GY under drought stress condition. Further, a functional marker was developed to predict the drought stress tolerance in a US maize inbred line panel. Lines carrying Hap2 exhibited > 10 % higher GY than those carrying Hap1 under drought stress condition. In Arabidopsis, overexpression ZmNAC080308 enhanced drought tolerance.

CONCLUSIONS

ZmNAC080308 is an important gene responding to drought tolerance, a functional marker is developed for improving maize drought tolerance by selecting this gene.

Genome-Wide Association Study of the Genetic Basis of Effective Tiller Number in Rice

BACKGROUND

Effective tiller number (ETN) has a pivotal role in determination of rice (Oryza sativa L.) grain yield. ETN is a complex quantitative trait regulated by both genetic and environmental factors. Despite multiple tillering-related genes have been cloned previously, few of them have been utilized in practical breeding programs.

RESULTS

In this study, we conducted a genome-wide association study (GWAS) for ETN using a panel of 490 rice accessions derived from the 3 K rice genomes project. Thirty eight ETN-associated QTLs were identified, interestingly, four of which colocalized with the OsAAP1, DWL2, NAL1, and OsWRKY74 gene previously reported to be involved in rice tillering regulation. Haplotype (Hap) analysis revealed that Hap5 of OsAAP1, Hap3 and 6 of DWL2, Hap2 of NAL1, and Hap3 and 4 of OsWRKY74 are favorable alleles for ETN. Pyramiding favorable alleles of all these four genes had more enhancement in ETN than accessions harboring the favorable allele of only one gene. Moreover, we identified 25 novel candidate genes which might also affect ETN, and the positive association between expression levels of the OsPILS6b gene and ETN was validated by RT-qPCR. Furthermore, transcriptome analysis on data released on public database revealed that most ETN-associated genes showed a relatively high expression from 21 days after transplanting (DAT) to 49 DAT and decreased since then. This unique expression pattern of ETN-associated genes may contribute to the transition from vegetative to reproductive growth of tillers.

CONCLUSIONS

Our results revealed that GWAS is a feasible way to mine ETN-associated genes. The candidate genes and favorable alleles identified in this study have the potential application value in rice molecular breeding for high ETN and grain yield.

Using high-throughput multiple optical phenotyping to decipher the genetic architecture of maize drought tolerance

BACKGROUND

Drought threatens the food supply of the world population. Dissecting the dynamic responses of plants to drought will be beneficial for breeding drought-tolerant crops, as the genetic controls of these responses remain largely unknown.

RESULTS

Here we develop a high-throughput multiple optical phenotyping system to noninvasively phenotype 368 maize genotypes with or without drought stress over a course of 98 days, and collected multiple optical images, including color camera scanning, hyperspectral imaging, and X-ray computed tomography images. We develop high-throughput analysis pipelines to extract image-based traits (i-traits). Of these i-traits, 10,080 were effective and heritable indicators of maize external and internal drought responses. An i-trait-based genome-wide association study reveals 4322 significant locus-trait associations, representing 1529 quantitative trait loci (QTLs) and 2318 candidate genes, many that co-localize with previously reported maize drought responsive QTLs. Expression QTL (eQTL) analysis uncovers many local and distant regulatory variants that control the expression of the candidate genes. We use genetic mutation analysis to validate two new genes, ZmcPGM2 and ZmFAB1A, which regulate i-traits and drought tolerance. Moreover, the value of the candidate genes as drought-tolerant genetic markers is revealed by genome selection analysis, and 15 i-traits are identified as potential markers for maize drought tolerance breeding.

CONCLUSION

Our study demonstrates that combining high-throughput multiple optical phenotyping and GWAS is a novel and effective approach to dissect the genetic architecture of complex traits and clone drought-tolerance associated genes.

Natural Variation in Crops: Realized Understanding, Continuing Promise

Crops feed the world's population and shape human civilization. The improvement of crop productivity has been ongoing for almost 10,000 years and has evolved from an experience-based to a knowledge-driven practice over the past three decades. Natural alleles and their reshuffling are long-standing genetic changes that affect how crops respond to various environmental conditions and agricultural practices. Decoding the genetic basis of natural variation is central to understanding crop evolution and, in turn, improving crop breeding. Here, we review current advances in the approaches used to map the causal alleles of natural variation, provide refined insights into the genetics and evolution of natural variation, and outline how this knowledge promises to drive the development of sustainable agriculture under the dome of emerging technologies.

Multiple functions of the vacuole in plant growth and fruit quality

Vacuoles are organelles in plant cells that play pivotal roles in growth and developmental regulation. The main functions of vacuoles include maintaining cell acidity and turgor pressure, regulating the storage and transport of substances, controlling the transport and localization of key proteins through the endocytic and lysosomal-vacuolar transport pathways, and responding to biotic and abiotic stresses. Further, proteins localized either in the tonoplast (vacuolar membrane) or inside the vacuole lumen are critical for fruit quality. In this review, we summarize and discuss some of the emerging functions and regulatory mechanisms associated with plant vacuoles, including vacuole biogenesis, vacuole functions in plant growth and development, fruit quality, and plant-microbe interaction, as well as some innovative research technology that has driven advances in the field. Together, the functions of plant vacuoles are important for plant growth and fruit quality. The investigation of vacuole functions in plants is of great scientific significance and has potential applications in agriculture.

Metabolomics-driven gene mining and genetic improvement of tolerance to salt-induced osmotic stress in maize

The farmland of the world's main corn-producing area is increasingly affected by salt stress. Therefore, the breeding of salt-tolerant cultivars is necessary for the long-term sustainability of global corn production. Previous studies have shown that natural maize varieties display a large diversity of salt tolerance, yet the genetic variants underlying such diversity remain poorly discovered and applied, especially those mediating the tolerance to salt-induced osmotic stress (SIOS). Here we report a metabolomics-driven understanding and genetic improvement of maize SIOS tolerance. Using a LC-MS-based untargeted metabolomics approach, we profiled the metabolomes of 266 maize inbred lines under control and salt conditions, and then identified 37 metabolite biomarkers of SIOS tolerance (METO1-37). Follow-up metabolic GWAS (mGWAS) and genotype-to-phenotype modeling identified 10 candidate genes significantly associating with the SIOS tolerance and METO abundances. Furthermore, we validated that a citrate synthase, a glucosyltransferase and a cytochrome P450 underlie the genotype-METO-SIOS tolerance associations, and showed that their favorable alleles additively improve the SIOS tolerance of elite maize inbred lines. Our study provides a novel insight into the natural variation of maize SIOS tolerance, which boosts the genetic improvement of maize salt tolerance, and demonstrates a metabolomics-based approach for mining crop genes associated with this complex agronomic trait.

Sat-BSA: an NGS-based method using local de novo assembly of long reads for rapid identification of genomic structural variations associated with agronomic traits

Advances in next generation sequencing (NGS)-based methodologies have accelerated the identifications of simple genetic variants such as point mutations and small insertions/deletions (InDels). Structural variants (SVs) including large InDels and rearrangements provide vital sources of genetic diversity for plant breeding. However, their analysis remains a challenge due to their complex nature. Consequently, novel NGS-based approaches are needed to rapidly and accurately identify SVs. Here, we present an NGS-based bulked-segregant analysis (BSA) technique called Sat-BSA (SVs associated with traits) for identifying SVs controlling traits of interest in crops. Sat-BSA targets allele frequencies at all SNP positions to first identify candidate genomic regions associated with a trait, which is then reconstructed by long reads-based local de novo assembly. Finally, the association between SVs, RNA-seq-based gene expression patterns and trait is evaluated for multiple cultivars to narrow down the candidate genes. We applied Sat-BSA to segregating F₂ progeny obtained from crosses between turnip cultivars with different tuber colors and successfully isolated two genes harboring SVs that are responsible for tuber phenotypes. The current study demonstrates the utility of Sat-BSA for the identification of SVs associated with traits of interest in species with large and heterozygous genomes.

Large-fragment insertion activates gene GaFZ (Ga08G0121) and is associated with the fuzz and trichome reduction in cotton (Gossypium arboreum)

Cotton seeds are typically covered by lint and fuzz fibres. Natural 'fuzzless' mutants are an ideal model system for identifying genes that regulate cell initiation and elongation. Here, using a genome-wide association study (GWAS), we identified a ~ 6.2 kb insertion, larINDEL_FZ , located at the end of chromosome 8, composed of a ~ 5.0 kb repetitive sequence and a ~ 1.2 kb fragment translocated from chromosome 12 in fuzzless Gossypium arboreum. The presence of larINDEL_FZ was associated with a fuzzless seed and reduced trichome phenotypes in G. arboreum. This distant insertion was predicted to be an enhancer, located ~ 18 kb upstream of the dominant-repressor GaFZ (Ga08G0121). Ectopic overexpression of GaFZ in Arabidopsis thaliana and G. hirsutum suggested that GaFZ negatively modulates fuzz and trichome development. Co-expression and interaction analyses demonstrated that GaFZ might impact fuzz fibre/trichome development by repressing the expression of genes in the very-long-chain fatty acid elongation pathway. Thus, we identified a novel regulator of fibre/trichome development while providing insights into the importance of noncoding sequences in cotton.

The genetic architecture of the dynamic changes in grain moisture in maize

Low grain moisture at harvest is crucial for safe production, transport and storage, but the genetic architecture of this trait in maize (Zea mays) remains elusive. Here, we measured the dynamic changes in grain moisture content in an association-mapping panel of 513 diverse maize inbred lines at five successive stages across five geographical environments. Genome-wide association study (GWAS) revealed 71 quantitative trait loci (QTLs) that influence grain moisture in maize. Epistatic effects play vital roles in the variability in moisture levels, even outperforming main-effect QTLs during the early dry-down stages. Distinct QTL-environment interactions influence the spatio-temporal variability of maize grain moisture, which is primarily triggered at specific times. By combining genetic population analysis, transcriptomic profiling and gene editing, we identified GRMZM5G805627 and GRMZM2G137211 as candidate genes underlying major QTLs for grain moisture in maize. Our results provide insights into the genetic architecture of dynamic changes in grain moisture, which should facilitate maize breeding.

Transcription-associated metabolomic adjustments in maize occur during combined drought and cold stress

Although simultaneous drought and cold stress occurs, especially in northwestern and eastern regions of China, and is an important factor limiting agricultural productivity, there are few studies focusing on plant responses to a combination of drought and cold stress. Here, by partially overlapping drought and cold stresses, we characterized the acclimation of maize (Zea mays B73) to these two stresses using physiological measurements, as well as comparative transcriptomics combined with metabolomics and hormonal analyses during the stress treatments and recovery stages. The combined drought and cold stress and drought stress alone were accompanied by a decline in photosynthetic capacity and enhanced transcriptional response, and subsequent recovery of these following removal from stress, whereas cold stress alone was accompanied by irreversible damage to photosynthetic capacity and chloroplast structure. The stress combination induced transcription-associated metabolomic alterations, in which raffinose, trehalose-6-phosphate, and proline accumulated, and monosaccharide abundance increased. Concomitantly, the increased abscisic acid (ABA) content and upregulated ABA signaling pathway may have provided the transcriptional regulation for the metabolic changes. In a parallel experiment, ABA treatments prior to exposure of the plants to cold stress primed the plants to survive the cold stress, thus confirming a key role for the endogenous ABA activated by the drought pretreatment in acclimation of the plants to cold. We present a model showing that the plant response to the combined stress is multi-faceted and reveal an ABA-dependent maize acclimation mechanism to the stress combination.

Genome-wide association mapping in maize: status and prospects

Genome-wide association study (GWAS) provides a robust and potent tool to retrieve complex phenotypic traits back to their underlying genetics. Maize is an excellent crop for performing GWAS due to diverse genetic variability, rapid decay of linkage disequilibrium, availability of distinct sub-populations and abundant SNP information. The application of GWAS in maize has resulted in successful identification of thousands of genomic regions associated with many abiotic and biotic stresses. Many agronomic and quality traits of maize are severely affected by such stresses and, significantly affecting its growth and productivity. To improve productivity of maize crop in countries like India which contribute only 2% to the world's total production in 2019-2020, it is essential to understand genetic complexity of underlying traits. Various DNA markers and trait associations have been revealed using conventional linkage mapping methods. However, it has achieved limited success in improving polygenic complex traits due to lower resolution of trait mapping. The present review explores the prospects of GWAS in improving yield, quality and stress tolerance in maize besides, strengths and challenges of using GWAS for molecular breeding and genomic selection. The information gathered will facilitate elucidation of genetic mechanisms of complex traits and improve efficiency of marker-assisted selection in maize breeding.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02799-4.

A capture-based assay for detection and characterization of transposon polymorphisms in maize

Transposons can create allelic diversity that affects gene expression and phenotypic diversity. The detection of transposon polymorphisms at a genome-wide scale across a large population is difficult. Here we developed a targeted sequencing approach to monitor transposon polymorphisms of interest. This approach can interrogate the presence or absence of transposons reliably across various genotypes. Using this approach, we genotyped a set of 965 transposon-related presence/absence polymorphisms in a diverse panel of 16 maize (Zea mays L.) inbred lines that are representative of the major maize breeding groups. About 70% of the selected regions can be effectively assayed in each genotype. The consistency between the capture-based assay and PCR-based assay are 98.6% based on analysis of 24 randomly selected transposon polymorphisms. By integrating the transposon polymorphisms data with gene expression data, ∼18% of the assayed transposon polymorphisms were found to be associated with variable gene expression levels. A detailed analysis of 18 polymorphisms in a larger association panel confirmed the effects of 10 polymorphisms, with one of them having stronger association with expression than nearby SNP markers. The effects of seven polymorphisms were tested using a luciferase-based expression assay, and one was confirmed. Together, this study demonstrates that the targeted sequencing assay is an effective way to explore transposon function in a high-throughput manner.

Zea mays RNA-seq estimated transcript abundances are strongly affected by read mapping bias

Background

Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines.

Results

In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2–4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters.

Conclusions

Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07577-3.

Association of TaD14-4D, a Gene Involved in Strigolactone Signaling, with Yield Contributing Traits in Wheat

Tillering is a crucial agronomic trait of wheat; it determines yield and plant architecture. Strigolactones (SLs) have been reported to inhibit plant branching. D14, a receptor of SLs, has been described to affect tillering in rice, yet it has seldomly been studied in wheat. In this study, three TaD14 homoeologous genes, TaD14-4A, TaD14-4B, and TaD14-4D, were identified. TaD14-4A, TaD14-4B, and TaD14-4D were constitutively expressed, and TaD14-4D had a higher expression level in most tissues. TaD14 proteins were localized in both cytoplasm and nucleus. An SNP and a 22 bp insertion/deletion (Indel) at the exon regions of TaD14-4D were detected, forming three haplotypes, namely 4D-HapI, 4D-HapII, and 4D-HapIII. Due to the frameshift mutation in the coding region of 4D-HapII, the interaction of 4D-HapII with TaMAX2 and TaD53 was blocked, which led to the blocking of SL signal transduction. Based on the two variation sites, two molecular markers, namely dCAPS-250 and Indel-747, were developed. Association analysis suggested that haplotypes of TaD14-4D were associated with effective tillering number (ETN) and thousand kernel weight (TKW) simultaneously in four environments. The favorable haplotype 4D-HapIII underwent positive selection in global wheat breeding. This study provides insights into understanding the function of natural variations of TaD14-4D and develops two useful molecular markers for wheat breeding.

A novel deletion in FLOWERING LOCUS T modulates flowering time in winter oilseed rape

A novel structural variant was discovered in the FLOWERING LOCUS T orthologue BnaFT.A02 by long-read sequencing. Nested association mapping in an elite winter oilseed rape population revealed that this 288 bp deletion associates with early flowering, putatively by modification of binding-sites for important flowering regulation genes. Perfect timing of flowering is crucial for optimal pollination and high seed yield. Extensive previous studies of flowering behavior in Brassica napus (canola, rapeseed) identified mutations in key flowering regulators which differentiate winter, semi-winter and spring ecotypes. However, because these are generally fixed in locally adapted genotypes, they have only limited relevance for fine adjustment of flowering time in elite cultivar gene pools. In crosses between ecotypes, the ecotype-specific major-effect mutations mask minor-effect loci of interest for breeding. Here, we investigated flowering time in a multiparental mapping population derived from seven elite winter oilseed rape cultivars which are fixed for major-effect mutations separating winter-type rapeseed from other ecotypes. Association mapping revealed eight genomic regions on chromosomes A02, C02 and C03 associating with fine modulation of flowering time. Long-read genomic resequencing of the seven parental lines identified seven structural variants coinciding with candidate genes for flowering time within chromosome regions associated with flowering time. Segregation patterns for these variants in the elite multiparental population and a diversity set of winter types using locus-specific assays revealed significant associations with flowering time for three deletions on chromosome A02. One of these was a previously undescribed 288 bp deletion within the second intron of FLOWERING LOCUS T on chromosome A02, emphasizing the advantage of long-read sequencing for detection of structural variants in this size range. Detailed analysis revealed the impact of this specific deletion on flowering-time modulation under extreme environments and varying day lengths in elite, winter-type oilseed rape.

Genetic dissection of drought resistance based on root traits at the bud stage in common bean

A whole-genome resequencing-derived SNP dataset used for genome-wide association analysis revealed 196 loci significantly associated with drought stress based on root traits. Candidate genes identified in the regions of these loci include homologs of known drought resistance genes in A. thaliana. Drought is the main abiotic constraint of the production of common bean. Improved adaptation to drought environments has become a main goal of crop breeding due to the increasing scarcity of water that will occur in the future. The overall objective of our study was to identify genomic regions associated with drought resistance based on root traits using genome-wide association analysis. A natural population of 438 common bean accessions was evaluated for root traits: root surface area, root average diameter, root volume, total root length, taproot length, lateral root number, root dry weight, lateral root length, special root weight/length, using seed germination pouches under drought conditions and in well-watered environments. The coefficient of variation ranged from 11.24% (root average diameter) to 38.19% (root dry weight) in the well-watered environment and from 9.61% (root average diameter) to 39.05% (lateral root length) under drought stress. A whole-genome resequencing-derived SNP dataset revealed 196 loci containing 230 candidate SNPs associated with drought resistance. Seventeen candidate SNPs were simultaneously associated with more than two traits. Forty-one loci were simultaneously associated with more than two traits, and eleven loci were colocated with loci previously reported to be related to drought resistance. Candidate genes of the associated loci included the ABA-responsive element-binding protein family, MYB, NAC, the protein kinase superfamily, etc. These results revealed promising alleles linked to drought resistance or root traits, providing insights into the genetic basis of drought resistance and roots, which will be useful for common bean improvement.

Recent advances in Cannabis sativa genomics research

Cannabis (Cannabis sativa L.) is one of the oldest cultivated plants purported to have unique medicinal properties. However, scientific research of cannabis has been restricted by the Single Convention on Narcotic Drugs of 1961, an international treaty that prohibits the production and supply of narcotic drugs except under license. Legislation governing cannabis cultivation for research, medicinal and even recreational purposes has been relaxed recently in certain jurisdictions. As a result, there is now potential to accelerate cultivar development of this multi-use and potentially medically useful plant species by application of modern genomics technologies. Whilst genomics has been pivotal to our understanding of the basic biology and molecular mechanisms controlling key traits in several crop species, much work is needed for cannabis. In this review we provide a comprehensive summary of key cannabis genomics resources and their applications. We also discuss prospective applications of existing and emerging genomics technologies for accelerating the genetic improvement of cannabis.

Genome-wide Identification and Characterization of FCS-Like Zinc Finger (FLZ) Family Genes in Maize (Zea mays) and Functional Analysis of ZmFLZ25 in Plant Abscisic Acid Response

FCS-like zinc finger family proteins (FLZs), a class of plant-specific scaffold of SnRK1 complex, are involved in the regulation of various aspects of plant growth and stress responses. Most information of FLZ family genes was obtained from the studies in Arabidopsis thaliana, whereas little is known about the potential functions of FLZs in crop plants. In this study, 37 maize FLZ (ZmFLZ) genes were identified to be asymmetrically distributed on 10 chromosomes and can be divided into three subfamilies. Protein interaction and subcellular localization assays demonstrated that eight typical ZmFLZs interacted and partially co-localized with ZmKIN10, the catalytic α-subunit of the SnRK1 complex in maize leaf mesophyll cells. Expression profile analysis revealed that several ZmFLZs were differentially expressed across various tissues and actively responded to diverse abiotic stresses. In addition, ectopic overexpression of ZmFLZ25 in Arabidopsis conferred hypersensitivity to exogenous abscisic acid (ABA) and triggered higher expression of ABA-induced genes, pointing to the positive regulatory role of ZmFLZ25 in plant ABA signaling, a scenario further evidenced by the interactions between ZmFLZ25 and ABA receptors. In summary, these data provide the most comprehensive information on FLZ family genes in maize, and shed light on the biological function of ZmFLZ25 in plant ABA signaling.

Cotton seedling drought tolerance is improved via salt preconditioning

In this study, 12 upland cotton seedlings were used as the material, and four treatments were designed (15% PEG for 6 h, 250 mM NaCl for 3 h, 15% PEG for 6 h after 250 mM NaCl pretreatment, and blank control). Various physiological indicators, including the malondialdehyde (MDA) and proline (Pro) contents and superoxide dismutase (SOD) and peroxidase (POD) activities, and the relative electrolyte leakage (REL), were measured during exposure to the aforementioned stresses, and three stress-related transcription factors (GhHsfA, GhbZIP, and GhNAC) were used to assess the differences in the drought resistance of cotton during exposure to PEG stress and NaCl/PEG combined stress. The analyses of the physiological and biochemical indicators revealed that the cotton seedlings exposed to NaCl/PEG combined stress exhibited the highest relative changes in the SOD and POD enzyme activities, while the relative changes in the MDA content and REL were relatively small. The cluster analysis showed that the treatments could be ranked as follows based on degree of damage exhibited by the exposed cotton seedlings: PEG > NaCl > NaCl/PEG. The exposure of cotton to NaCl/PEG combined stress resulted in a lower degree of damage than that obtained after exposure to PEG alone, which indicated that an appropriate amount of NaCl could partially relieve the adverse effects of drought on cotton seedlings. In addition, the relative expression levels of GhHsfA, GhbZIP, and GhNAC were significantly correlated with multiple physiological and biochemical indicators under different stresses, and the principal component analysis identified these transcription factors as important indicators. Based on these findings, these three transcription factors can be used as molecular indicators for the identification of drought resistance. A comprehensive D value cluster analysis ranked the 12 cotton varieties based on their drought resistance, and the most drought-resistant variety was ND359-5. This study provides new methods and materials for research on drought resistance in cotton.

Status and prospects of genome-wide association studies in plants

Genome-wide association studies (GWAS) have developed into a powerful and ubiquitous tool for the investigation of complex traits. In large part, this was fueled by advances in genomic technology, enabling us to examine genome-wide genetic variants across diverse genetic materials. The development of the mixed model framework for GWAS dramatically reduced the number of false positives compared with naïve methods. Building on this foundation, many methods have since been developed to increase computational speed or improve statistical power in GWAS. These methods have allowed the detection of genomic variants associated with either traditional agronomic phenotypes or biochemical and molecular phenotypes. In turn, these associations enable applications in gene cloning and in accelerated crop breeding through marker assisted selection or genetic engineering. Current topics of investigation include rare-variant analysis, synthetic associations, optimizing the choice of GWAS model, and utilizing GWAS results to advance knowledge of biological processes. Ongoing research in these areas will facilitate further advances in GWAS methods and their applications.

Genetic and molecular control of grain yield in maize

Understanding the genetic and molecular basis of grain yield is important for maize improvement. Here, we identified 49 consensus quantitative trait loci (cQTL) controlling maize yield-related traits using QTL meta-analysis. Then, we collected yield-related traits associated SNPs detected by association mapping and identified 17 consensus significant loci. Comparing the physical positions of cQTL with those of significant SNPs revealed that 47 significant SNPs were located within 20 cQTL regions. Furthermore, intensive reviews of 31 genes regulating maize yield-related traits found that the functions of many genes were conservative in maize and other plant species. The functional conservation indicated that some of the 575 maize genes (orthologous to 247 genes controlling yield or seed traits in other plant species) might be functionally related to maize yield-related traits, especially the 49 maize orthologous genes in cQTL regions, and 41 orthologous genes close to the physical positions of significant SNPs. In the end, we prospected on the integration of the public sources for exploring the genetic and molecular mechanisms of maize yield-related traits, and on the utilization of genetic and molecular mechanisms for maize improvement.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01214-3.

Genome-Wide Association Study Dissects Resistance Loci against Bacterial Blight in a Diverse Rice Panel from the 3000 Rice Genomes Project

BACKGROUND

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo) is one of the most devastating bacterial diseases of rice in temperate and tropical regions. Breeding and deployment of resistant cultivars carrying major resistance (R) genes has been the most effective approach for BB management. However, because of specific interaction of each R gene with the product of the corresponding pathogen avirulence or effector gene, new pathogen strains that can overcome the deployed resistance often emerge rapidly. To deal with ever-evolving Xoo, it is necessary to identify novel R genes and resistance quantitative trait loci (QTL).

RESULTS

BB resistance of a diverse panel of 340 accessions from the 3000 Rice Genomes Project (3 K RGP) was evaluated by artificial inoculation with four representative Xoo strains, namely Z173 (C4), GD1358 (C5), V from China and PXO339 (P9a) from Philippines. Using the 3 K RG 4.8mio filtered SNP Dataset, a total of 11 QTL associated with BB resistance on chromosomes 4, 5, 11 and 12 were identified through a genome-wide association study (GWAS). Among them, eight resistance loci, which were narrowed down to relatively small genomic intervals, coincided with previously reported QTL or R genes, e.g. xa5, xa25, xa44(t). The other three QTL were putative novel loci associated with BB resistance. Linear regression analysis showed a dependence of BB lesion length on the number of favorable alleles, suggesting that pyramiding QTL using marker-assisted selection would be an effective approach for improving resistance. In addition, the Hap2 allele of LOC_Os11g46250 underlying qC5-11.1 was validated as positively regulating resistance against strain C5.

CONCLUSIONS

Our findings provide valuable information for the genetic improvement of BB resistance and application of germplasm resources in rice breeding programs.

Genetic insights into natural variation underlying salt tolerance in wheat

Developing salt-tolerant crop varieties is one of the important approaches to cope with increasing soil salinization worldwide. In this study, a diversity panel of 323 wheat accessions and 150 doubled haploid lines were phenotyped for salt-responsive morphological and physiological traits across two growth stages. The comprehensive salt tolerance of each wheat accession was evaluated based on principal component analysis. A total of 269 associated loci for salt-responsive traits and/or salt tolerance indices were identified by genome-wide association studies using 395 675 single nucleotide polymorphisms, among which 22 overlapping loci were simultaneously identified by biparental quantitative trait loci mapping. Two novel candidate genes ROOT NUMBER 1 (TaRN1) and ROOT NUMBER 2 (TaRN2) involved in root responses to salt stress fell within overlapping loci, showing different expression patterns and a frameshift mutation (in TaRN2) in contrasting salt-tolerant wheat genotypes. Moreover, the decline in salt tolerance of Chinese wheat varieties was observed from genetic and phenotypic data. We demonstrate that a haplotype controlling root responses to salt stress has been diminished by strong selection for grain yield, which highlights that linkage drag constrains the salt tolerance of Chinese wheat. This study will facilitate salt-tolerant wheat breeding in terms of elite germplasm, favorable alleles and selection strategies.

Transcriptome analysis and molecular mechanism of linseed (Linum usitatissimum L.) drought tolerance under repeated drought using single-molecule long-read sequencing

Background

Oil flax (linseed, Linum usitatissimum L.) is one of the most important oil crops., However, the increases in drought resulting from climate change have dramatically reduces linseed yield and quality, but very little is known about how linseed coordinates the expression of drought resistance gene in response to different level of drought stress (DS) on the genome-wide level.

Results

To explore the linseed transcriptional response of DS and repeated drought (RD) stress, we determined the drought tolerance of different linseed varieties. Then we performed full-length transcriptome sequencing of drought-resistant variety (Z141) and drought-sensitive variety (NY-17) under DS and RD stress at the seedling stage using single-molecule real-time sequencing and RNA-sequencing. Gene Ontology (GO) and reduce and visualize GO (REVIGO) enrichment analysis showed that upregulated genes of Z141 were enriched in more functional pathways related to plant drought tolerance than those of NY-17 were under DS. In addition, 4436 linseed transcription factors were identified, and 1190 were responsive to stress treatments. Moreover, protein-protein interaction (PPI) network analysis showed that the proline biosynthesis pathway interacts with stress response genes through RAD50 (DNA repair protein 50) interacting protein 1 (RIN-1). Finally, proline biosynthesis and DNA repair structural gene expression patterns were verified by RT- PCR.

Conclusions

The drought tolerance of Z141 may be related to its upregulation of drought tolerance genes under DS. Proline may play an important role in linseed drought tolerance by maintaining cell osmotic and protecting DNA from ROS damage. In summary, this study provides a new perspective to understand the drought adaptability of linseed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07416-5.

Genetic dissection of maize drought tolerance for trait improvement

Maize is one of the most important crops, but its production is threatened by drought stress worldwide. Thus, increased drought tolerance has been a major goal of maize breeding. Conventional breeding strategies have led to significantly increase of maize yields; however, these strategies often fail to meet the need for drought stress tolerance enhancement. Here, we focus on progress related to the genetic dissection of drought tolerance in maize at different developmental stages achieved through linkage mapping and association mapping. Moreover, recent molecular breeding systems, including transgenic, genome-wide marker-assisted selection, and genome editing technologies, have provided a more direct, efficient, and accurate approach for trait improvement. We also provide perspectives on future directions regarding multi-omics studies and maize improvement. Overall, the application of acquired knowledge will facilitate maize breeding to meet the challenges.

Uncovering the genetic mechanisms regulating panicle architecture in rice with GPWAS and GWAS

Background

The number of panicles per plant, number of grains per panicle, and 1000-grain weight are important factors contributing to the grain yield per plant in rice. The Rice Diversity Panel 1 (RDP1) contains a total of 421 purified, homozygous rice accessions representing diverse genetic variations within O. sativa. The release of High-Density Rice Array (HDRA, 700 k SNPs) dataset provides a new opportunity to discover the genetic variants of panicle architectures in rice.

Results

In this report, a new method genome-phenome wide association study (GPWAS) was performed with 391 individuals and 27 traits derived from RDP1 to scan the relationship between the genes and multi-traits. A total of 1985 gene models were linked to phenomic variation with a p-value cutoff of 4.49E-18. Besides, 406 accessions derived from RDP1 with 411,066 SNPs were used to identify QTLs associated with the total spikelets number per panicle (TSNP), grain number per panicle (GNP), empty grain number per panicle (EGNP), primary branch number (PBN), panicle length (PL), and panicle number per plant (PN) by GLM, MLM, FarmCPU, and BLINK models for genome-wide association study (GWAS) analyses. A total of 18, 21, 18, 17, 15, and 17 QTLs were identified tightly linked with TSNP, GNP, EGNP, PBN, PL, and PN, respectively. Then, a total of 23 candidate genes were mapped simultaneously using both GWAS and GPWAS methods, composed of 6, 4, 5, 4, and 4 for TSNP, GNP, EGNP, PBN, and PL. Notably, one overlapped gene (Os01g0140100) were further investigated based on the haplotype and gene expression profile, indicating this gene might regulate the TSNP or panicle architecture in rice.

Conclusions

Nearly 30 % (30/106) QTLs co-located with the previous published genes or QTLs, indicating the power of GWAS. Besides, GPWAS is a new method to discover the relationship between genes and traits, especially the pleiotropy genes. Through comparing the results from GWAS and GPWAS, we identified 23 candidate genes related to panicle architectures in rice. This comprehensive study provides new insights into the genetic basis controlling panicle architectures in rice, which lays a foundation in rice improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07391-x.

Genome-wide association study of six quality-related traits in common wheat (Triticum aestivum L.) under two sowing conditions

We identified genomic regions associated with six quality-related traits in wheat under two sowing conditions and analyzed the effects of multienvironment-significant SNPs on the stability of these traits. Grain quality affects the nutritional and commercial value of wheat (Triticum aestivum L.) and is a critical factor influencing consumer preferences for specific wheat varieties. Climate change is predicted to increase environmental stress and thereby reduce wheat quality. Here, we performed a genotyping assay involving the use of the wheat 90 K array in a genome-wide association study of six quality-related traits in 486 wheat accessions under two sowing conditions (normal and late sowing) over 4 years. We identified 64 stable quantitative trait loci (QTL), including 10 for grain protein content, 9 for wet gluten content, 4 for grain starch content, 14 for water absorption, 15 for dough stability time and 12 for grain hardness in wheat under two sowing conditions. These QTL harbored 175 single nucleotide polymorphisms (SNPs), explaining approximately 3-13% of the phenotypic variation in multiple environments. Some QTL on chromosomes 6A and 5D were associated with multiple traits simultaneously, and two (QNGPC.cau-6A, QNGH.cau-5D) harbored known genes, such as NAM-A1 for grain protein content and Pinb for grain hardness, whereas other QTL could facilitate gene discovery. Forty-three SNPs that were detected under late or both normal and late sowing conditions appear to be related to phenotypic stability. The effects of these SNP alleles were confirmed in the association population. The results of this study will be useful for further dissecting the genetic basis of quality-related traits in wheat and developing new wheat cultivars with desirable alleles to improve the stability of grain quality.

Evolution and Adaptation of Forest and Crop Pathogens in the Anthropocene

Anthropocene marks the era when human activity is making a significant impact on earth, its ecological and biogeographical systems. The domestication and intensification of agricultural and forest production systems have had a large impact on plant and tree health. Some pathogens benefitted from these human activities and have evolved and adapted in response to the expansion of crop and forest systems, resulting in global outbreaks. Global pathogen genomics data including population genomics and high-quality reference assemblies are crucial for understanding the evolution and adaptation of pathogens. Crops and forest trees have remarkably different characteristics, such as reproductive time and the level of domestication. They also have different production systems for disease management with more intensive management in crops than forest trees. By comparing and contrasting results from pathogen population genomic studies done on widely different agricultural and forest production systems, we can improve our understanding of pathogen evolution and adaptation to different selection pressures. We find that in spite of these differences, similar processes such as hybridization, host jumps, selection, specialization, and clonal expansion are shaping the pathogen populations in both crops and forest trees. We propose some solutions to reduce these impacts and lower the probability of global pathogen outbreaks so that we can envision better management strategies to sustain global food production as well as ecosystem services.

The FUSED LEAVES1-ADHERENT1 regulatory module is required for maize cuticle development and organ separation

All aerial epidermal cells in land plants are covered by the cuticle, an extracellular hydrophobic layer that provides protection against abiotic and biotic stresses and prevents organ fusion during development. Genetic and morphological analysis of the classic maize adherent1 (ad1) mutant was combined with genome-wide binding analysis of the maize MYB transcription factor FUSED LEAVES1 (FDL1), coupled with transcriptional profiling of fdl1 mutants. We show that AD1 encodes an epidermally-expressed 3-KETOACYL-CoA SYNTHASE (KCS) belonging to a functionally uncharacterized clade of KCS enzymes involved in cuticular wax biosynthesis. Wax analysis in ad1 mutants indicates that AD1 functions in the formation of very-long-chain wax components. We demonstrate that FDL1 directly binds to CCAACC core motifs present in AD1 regulatory regions to activate its expression. Over 2000 additional target genes of FDL1, including many involved in cuticle formation, drought response and cell wall organization, were also identified. Our results identify a regulatory module of cuticle biosynthesis in maize that is conserved across monocots and eudicots, and highlight previously undescribed factors in lipid metabolism, transport and signaling that coordinate organ development and cuticle formation.

QTL Mapping for Leaf Area of Tea Plants (Camellia sinensis) Based on a High-Quality Genetic Map Constructed by Whole Genome Resequencing

High-quality genetic maps play important roles in QTL mapping and molecular marker-assisted breeding. Tea leaves are not only important vegetative organs but are also the organ for harvest with important economic value. However, the key genes and genetic mechanism of regulating leaf area have not been clarified. In this study, we performed whole-genome resequencing on "Jinxuan," "Yuncha 1" and their 96 F1 hybrid offspring. From the 1.84 Tb of original sequencing data, abundant genetic variation loci were identified, including 28,144,625 SNPs and 2,780,380 indels. By integrating the markers of a previously reported genetic map, a high-density genetic map consisting of 15 linkage groups including 8,956 high-quality SNPs was constructed. The total length of the genetic map is 1,490.81 cM, which shows good collinearity with the genome. A total of 25 representative markers (potential QTLs) related to leaf area were identified, and there were genes differentially expressed in large and small leaf samples near these markers. GWAS analysis further verified the reliability of QTL mapping. Thirty-one pairs of newly developed indel markers located near these potential QTLs showed high polymorphism and had good discrimination between large and small leaf tea plant samples. Our research will provide necessary support and new insights for tea plant genetic breeding, quantitative trait mapping and yield improvement.

Loci harboring genes with important role in drought and related abiotic stress responses in flax revealed by multiple GWAS models

QTNs associated with drought tolerance traits and indices were identified in a flax mini-core collection through multiple GWAS models and phenotyping at multiple locations under irrigated and non-irrigated field conditions. Drought is a critical phenomenon challenging today's agricultural sector. Crop varieties adapted to moisture deficit are becoming vital. Flax can be greatly affected by limiting moisture conditions, especially during the early development and reproductive stages. Here, a mini-core collection comprising genotypes from more than 20 major growing countries was evaluated for 11 drought-related traits in irrigated and non-irrigated fields for 3 years. Heritability of the traits ranged from 44.7 to 86%. Six of the 11 traits showed significant phenotypic difference between irrigated and non-irrigated conditions. A genome-wide association study (GWAS) was performed for these six traits and their corresponding stress indices based on 106 genotypes and 12,316 single nucleotide polymorphisms (SNPs) using six multi-locus and one single-locus models. The SNPs were then assigned to 8050 linkage disequilibrium (LD) blocks to which a restricted two-stage multi-locus multi-allele GWAS was applied. A total of 144 quantitative trait nucleotides (QTNs) and 13 LD blocks were associated with at least one trait or stress index. Of these, 16 explained more than 15% of the genetic variance. Most large-effect QTN loci harbored gene(s) previously predicted to play role(s) in the associated traits. Genes mediating responses to abiotic stresses resided at loci associated with stress indices. Flax genes Lus10009480 and Lus10030150 that are predicted to encode WAX INDUCER1 and STRESS-ASSOCIATED PROTEIN (SAP), respectively, are among the important candidates detected. Accessions with multiple favorable alleles outperformed others for grain yield, thousand seed weight and fiber/biomass in non-irrigated conditions, suggesting their potential usefulness in breeding and genomic selection.

Maize microRNA166 Inactivation Confers Plant Development and Abiotic Stress Resistance

MicroRNAs are important regulators in plant developmental processes and stress responses. In this study, we generated a series of maize STTM166 transgenic plants. Knock-down of miR166 resulted in various morphological changes, including rolled leaves, enhanced abiotic stress resistance, inferior yield-related traits, vascular pattern and epidermis structures, tassel architecture, as well as abscisic acid (ABA) level elevation and indole acetic acid (IAA) level reduction in maize. To profile miR166 regulated genes, we performed RNA-seq and qRT-PCR analysis. A total of 178 differentially expressed genes (DEGs) were identified, including 118 up-regulated and 60 down-regulated genes. These DEGs were strongly enriched in cell and intercellular components, cell membrane system components, oxidoreductase activity, single organism metabolic process, carbohydrate metabolic process, and oxidation reduction process. These results indicated that miR166 plays important roles in auxin and ABA interaction in monocots, yet the specific mechanism may differ from dicots. The enhanced abiotic stress resistance is partly caused via rolling leaves, high ABA content, modulated vascular structure, and the potential changes of cell membrane structure. The inferior yield-related traits and late flowering are partly controlled by the decreased IAA content, the interplay of miR166 with other miRNAs and AGOs. Taken together, the present study uncovered novel functions of miR166 in maize, and provide insights on applying short tandem target mimics (STTM) technology in plant breeding.

Distinct nucleotide patterns among three subgenomes of bread wheat and their potential origins during domestication after allopolyploidization

BACKGROUND

The speciation and fast global domestication of bread wheat have made a great impact on three subgenomes of bread wheat. DNA base composition is an essential genome feature, which follows the individual-strand base equality rule and [AT]-increase pattern at the genome, chromosome, and polymorphic site levels among thousands of species. Systematic analyses on base compositions of bread wheat and its wild progenitors could facilitate further understanding of the evolutionary pattern of genome/subgenome-wide base composition of allopolyploid species and its potential causes.

RESULTS

Genome/subgenome-wide base-composition patterns were investigated by using the data of polymorphic site in 93 accessions from worldwide populations of bread wheat, its diploid and tetraploid progenitors, and their corresponding reference genome sequences. Individual-strand base equality rule and [AT]-increase pattern remain in recently formed hexaploid species bread wheat at the genome, subgenome, chromosome, and polymorphic site levels. However, D subgenome showed the fastest [AT]-increase across polymorphic site from Aegilops tauschii to bread wheat than that on A and B subgenomes from wild emmer to bread wheat. The fastest [AT]-increase could be detected almost all chromosome windows on D subgenome, suggesting different mechanisms between D and other two subgenomes. Interestingly, the [AT]-increase is mainly contributed by intergenic regions at non-selective sweeps, especially the fastest [AT]-increase of D subgenome. Further transition frequency and sequence context analysis indicated that three subgenomes shared same mutation type, but D subgenome owns the highest mutation rate on high-frequency mutation type. The highest mutation rate on D subgenome was further confirmed by using a bread-wheat-private SNP set. The exploration of loci/genes related to the [AT] value of D subgenome suggests the fastest [AT]-increase of D subgenome could be involved in DNA repair systems distributed on three subgenomes of bread wheat.

CONCLUSIONS

The highest mutation rate is detected on D subgenome of bread wheat during domestication after allopolyploidization, leading to the fastest [AT]-increase pattern of D subgenome. The phenomenon may come from the joint action of multiple repair systems inherited from its wild progenitors.

Phenomics-based GWAS analysis reveals the genetic architecture for drought resistance in cotton

Drought resistance (DR) is a complex trait that is regulated by a variety of genes. Without comprehensive profiling of DR-related traits, the knowledge of the genetic architecture for DR in cotton remains limited. Thus, there is a need to bridge the gap between genomics and phenomics. In this study, an automatic phenotyping platform (APP) was systematically applied to examine 119 image-based digital traits (i-traits) during drought stress at the seedling stage, across a natural population of 200 representative upland cotton accessions. Some novel i-traits, as well as some traditional i-traits, were used to evaluate the DR in cotton. The phenomics data allowed us to identify 390 genetic loci by genome-wide association study (GWAS) using 56 morphological and 63 texture i-traits. DR-related genes, including GhRD2, GhNAC4, GhHAT22 and GhDREB2, were identified as candidate genes by some digital traits. Further analysis of candidate genes showed that Gh_A04G0377 and Gh_A04G0378 functioned as negative regulators for cotton drought response. Based on the combined digital phenotyping, GWAS analysis and transcriptome data, we conclude that the phenomics dataset provides an excellent resource to characterize key genetic loci with an unprecedented resolution which can inform future genome-based breeding for improved DR in cotton.

Natural Genetic Resources from Diverse Plants to Improve Abiotic Stress Tolerance in Plants

The current agricultural system is biased for the yield increase at the cost of biodiversity. However, due to the loss of precious genetic diversity during domestication and artificial selection, modern cultivars have lost the adaptability to cope with unfavorable environments. There are many reports on variations such as single nucleotide polymorphisms (SNPs) and indels in the stress-tolerant gene alleles that are associated with higher stress tolerance in wild progenitors, natural accessions, and extremophiles in comparison with domesticated crops or model plants. Therefore, to gain a better understanding of stress-tolerant traits in naturally stress-resistant plants, more comparative studies between the modern crops/model plants and crop progenitors/natural accessions/extremophiles are required. In this review, we discussed and summarized recent progress on natural variations associated with enhanced abiotic stress tolerance in various plants. By applying the recent biotechniques such as the CRISPR/Cas9 gene editing tool, natural genetic resources (i.e., stress-tolerant gene alleles) from diverse plants could be introduced to the modern crop in a non-genetically modified way to improve stress-tolerant traits.

Mapping Regulatory Determinants in Plants

The domestication and improvement of many plant species have frequently involved modulation of transcriptional outputs and continue to offer much promise for targeted trait engineering. The cis-regulatory elements (CREs) controlling these trait-associated transcriptional variants however reside within non-coding regions that are currently poorly annotated in most plant species. This is particularly true in large crop genomes where regulatory regions constitute only a small fraction of the total genomic space. Furthermore, relatively little is known about how CREs function to modulate transcription in plants. Therefore understanding where regulatory regions are located within a genome, what genes they control, and how they are structured are important factors that could be used to guide both traditional and synthetic plant breeding efforts. Here, we describe classic examples of regulatory instances as well as recent advances in plant regulatory genomics. We highlight valuable molecular tools that are enabling large-scale identification of CREs and offering unprecedented insight into how genes are regulated in diverse plant species. We focus on chromatin environment, transcription factor (TF) binding, the role of transposable elements, and the association between regulatory regions and target genes.