Genetic variation in ZmVPP1 contributes to drought tolerance in maize seedlings

Genome-Wide Analysis of TCP Family Genes in Zea mays L. Identified a Role for ZmTCP42 in Drought Tolerance

The Teosinte-branched 1/Cycloidea/Proliferating (TCP) plant-specific transcription factors (TFs) have been demonstrated to play a fundamental role in plant development and organ patterning. However, it remains unknown whether or not the TCP gene family plays a role in conferring a tolerance to drought stress in maize, which is a major constraint to maize production. In this study, we identified 46 ZmTCP genes in the maize genome and systematically analyzed their phylogenetic relationships and synteny with rice, sorghum, and ArabidopsisTCP genes. Expression analysis of the 46 ZmTCP genes in different tissues and under drought conditions, suggests their involvement in maize response to drought stress. Importantly, genetic variations in ZmTCP32 and ZmTCP42 are significantly associated with drought tolerance at the seedling stage. RT-qPCR results suggest that ZmTCP32 and ZmTCP42 RNA levels are both induced by ABA, drought, and polyethylene glycol treatments. Based on the significant association between the genetic variation of ZmTCP42 and drought tolerance, and the inducible expression of ZmTCP42 by drought stress, we selected ZmTCP42, to investigate its function in drought response. We found that overexpression of ZmTCP42 in Arabidopsis led to a hypersensitivity to ABA in seed germination and enhanced drought tolerance, validating its function in drought tolerance. These results suggested that ZmTCP42 functions as an important TCP TF in maize, which plays a positive role in drought tolerance.

Parallelism and convergence in post-domestication adaptation in cereal grasses

The selection of desirable traits in crops during domestication has been well studied. Many crops share a suite of modified phenotypic characteristics collectively known as the domestication syndrome. In this sense, crops have convergently evolved. Previous work has demonstrated that, at least in some instances, convergence for domestication traits has been achieved through parallel molecular means. However, both demography and selection during domestication may have placed limits on evolutionary potential and reduced opportunities for convergent adaptation during post-domestication migration to new environments. Here we review current knowledge regarding trait convergence in the cereal grasses and consider whether the complexity and dynamism of cereal genomes (e.g., transposable elements, polyploidy, genome size) helped these species overcome potential limitations owing to domestication and achieve broad subsequent adaptation, in many cases through parallel means. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Genome-wide association study dissects the genetic bases of salt tolerance in maize seedlings

Excess salinity is a natural stress that causes crop yield losses worldwide. The genetic bases of maize salt tolerance remain largely unknown. Here we investigated the survival rates of 445 maize natural accessions after salt treatments. A skewed distribution of the salt-tolerant phenotypes was observed in this population. Genome-wide association studies (GWAS) revealed 57 loci significantly associated with salt tolerance. Forty-nine candidate genes were detected from these loci. About 10% of these genes were co-localized with loci from QTL mapping. Forty four percent of the candidate genes were involved in stress responses, ABA signaling, stomata division, DNA binding/transcription regulation and auxin signaling, suggesting that they are key genetic mechanisms of maize salt tolerance. Transgenic studies showed that two genes, the salt-tolerance-associated-gene 4 (SAG4, GRMZM2G077295) and SAG6 (GRMZM2G106056), which encode a protein transport protein and the double-strand break repair protein MRE11, respectively, had positive roles in plant salt tolerance, and their salt-tolerant haplotypes were revealed. The genes we identified in this study provide a list of candidate targets for further study of maize salt tolerance, and of genetic markers and materials that may be used for breeding salt-tolerance in maize.

A Novel Candidate Gene Associated With Body Weight in the Pacific White Shrimp Litopenaeus vannamei

Improvements of growth traits are always the focus in selective breeding programs for the Pacific white shrimp Litopenaeus vannamei (L. vannamei). Identification of growth-related genes or markers can contribute to the application of modern breeding technologies, and thus accelerate the genetic improvement of growth traits. The aim of this study was to identify the genes and molecular markers associated with the growth traits of L. vannamei. A population of 200 individuals was genotyped using 2b-RAD techniques for genome-wide linkage disequilibrium (LD) analysis and genome-wide association study (GWAS). The results showed that the LD decayed fast in the studied population, which suggest that it is feasible to fine map the growth-related genes with GWAS in L. vannamei. One gene designated as LvSRC, encoding the class C scavenger receptor (SRC), was identified as a growth-related candidate gene by GWAS. Further targeted sequencing of the candidate gene in another population of 322 shrimps revealed that several non-synonymous mutations within LvSRC were significantly associated with the body weight (P < 0.01), and the most significant marker (SRC_24) located in the candidate gene could explain 13% of phenotypic variance. The current results provide not only molecular markers for genetic improvement in L. vannamei, but also new insights for understanding the growth regulation mechanism in penaeid shrimp.

A large-scale circular RNA profiling reveals universal molecular mechanisms responsive to drought stress in maize and Arabidopsis

Drought is a major abiotic stress that threatens global food security. Circular RNAs (circRNAs) are endogenous RNAs. How these molecules influence plant stress responses remains elusive. Here, a large-scale circRNA profiling identified 2174 and 1354 high-confidence circRNAs in maize and Arabidopsis, respectively, and most were differentially expressed in response to drought. A substantial number of drought-associated circRNA-hosting genes were involved in conserved or species-specific pathways in drought responses. Comparative analysis revealed that circRNA biogenesis was more complex in maize than in Arabidopsis. In most cases, maize circRNAs were negatively correlated with sRNA accumulation. In 368 maize inbred lines, the circRNA-hosting genes were enriched for single nucleotide polymorphisms (SNPs) associated with circRNA expression and drought tolerance, implying either important roles of circRNAs in maize drought responses or their potential use as biomarkers for breeding drought-tolerant maize. Additionally, the expression levels of circRNAs derived from drought-responsible genes encoding calcium-dependent protein kinase and cytokinin oxidase/dehydrogenase were significantly associated with drought tolerance of maize seedlings. Specifically, Arabidopsis plants overexpressing circGORK (Guard cell outward-rectifying K⁺ -channel) were hypersensitive to abscisic acid, but insensitive to drought, suggesting a positive role of circGORK in drought tolerance. We report the transcriptomic profiling and transgenic studies of circRNAs in plant drought responses, and provide evidence highlighting the universal molecular mechanisms involved in plant drought tolerance.

Genome-wide analysis of maize OSCA family members and their involvement in drought stress

BACKGROUND

Worldwide cultivation of maize is often impacted negatively by drought stress. Hyperosmolality-gated calcium-permeable channels (OSCA) have been characterized as osmosensors in Arabidopsis. However, the involvement of members of the maize OSCA (ZmOSCA) gene family in response to drought stress is unknown. It is furthermore unclear which ZmOSCA gene plays a major role in genetic improvement of drought tolerance in Maize.

METHODS

We predicted the protein domain structure and transmembrane regions by using the NCBI Conserved Domain Database database and TMHMM server separately. The phylogeny tree was built by Mega7. We used the mixed linear model in TASSEL to perform the family-based association analysis.

RESULTS

In this report, 12 ZmOSCA genes were uncovered in the maize genome by a genome-wide survey and analyzed systematically to reveal their synteny and phylogenetic relationship with the genomes of rice, maize, and sorghum. These analyses indicated a relatively conserved evolutionary history of the ZmOSCA gene family. Protein domain and transmembrane analysis indicated that most of the 12 ZmOSCAs shared similar structures with their homologs. The result of differential expression analysis under drought at various stages, as well as the expression profiles in 15 tissues, revealed a functional divergence of ZmOSCA genes. Notably, the expression level of ZmOSCA4.1 being up-regulated in both seedlings and adult leaves. Notably, the association analysis between genetic variations in these genes and drought tolerance was detected. Significant associations between genetic variation in ZmOSCA4.1 and drought tolerance were found at the seedling stage. Our report provides a detailed analysis of the ZmOSCAs in the maize genome. These findings will contribute to future studies on the functional characterization of ZmOSCA proteins in response to water deficit stress, as well as understanding the mechanism of genetic variation in drought tolerance in maize.

Excess Pyrophosphate within Guard Cells Delays Stomatal Closure

A variety of cellular metabolic reactions generate inorganic pyrophosphate (PPi) as an ATP hydrolysis byproduct. The vacuolar H+-translocating pyrophosphatase (H+-PPase) loss-of-function fugu5 mutant is susceptible to drought and displays pleotropic postgerminative growth defects due to excess PPi. It was recently reported that stomatal closure after abscisic acid (ABA) treatment is delayed in vhp1-1, a fugu5 allele. In contrast, we found that specific removal of PPi rescued all of the above fugu5 developmental and growth defects. Hence, we speculated that excess PPi itself, rather than vacuolar acidification, might delay stomatal closure. To test this hypothesis, we constructed transgenic plants expressing the yeast IPP1 gene (encoding a cytosolic pyrophosphatase) driven by a guard cell-specific promoter (pGC1::IPP1) in the fugu5 background. Our measurements confirmed stomatal closure defects in fugu5, further supporting a role for H+-PPase in stomatal functioning. Importantly, while pGC1::IPP1 transgenics morphologically mimicked fugu5, stomatal closure was restored in response to ABA and darkness. Quantification of water loss revealed that fugu5 stomata were almost completely insensitive to ABA. In addition, growth of pGC1::IPP1 plants was promoted compared to fugu5 throughout their life; however, it did not reach the wild type level. fugu5 also displayed an increased stomatal index, in violation of the one-cell-spacing rule, and phenotypes recovered upon removal of PPi by pAVP1::IPP1 (FUGU5, VHP1 and AVP1 are the same gene encoding H+-PPase), but not in the pGC1::IPP1 line. Taken together, these results clearly support our hypothesis that dysfunction in stomata is triggered by excess PPi within guard cells, probably via perturbed guard cell metabolism.

A Genome-Wide Association Study Identifies New Loci Involved in Wound-Induced Lateral Root Formation in Arabidopsis thaliana

Root systems can display variable architectures that contribute to nutrient foraging or to increase the tolerance of abiotic stress conditions. Root tip excision promotes the developmental progression of previously specified lateral root (LR) founder cells, which allows to easily measuring the branching capacity of a given root as regards its genotype and/or growth conditions. Here, we describe the natural variation among 120 Arabidopsis thaliana accessions in root system architecture (RSA) after root tip excision. Wound-induced changes in RSA were associated with 19 genomic loci using genome-wide association mapping. Three candidate loci associated with wound-induced LR formation were investigated. Sequence variation in the hypothetical protein encoded by the At4g01090 gene affected wound-induced LR development and its loss-of-function mutants displayed a reduced number of LRs after root tip excision. Changes in a histidine phosphotransfer protein putatively involved in cytokinin signaling were significantly associated with LR number variation after root tip excision. Our results provide a better understanding of some of the genetic components involved in LR capacity variation among accessions.

Connecting genome structural variation with complex traits in crop plants

Structural genome variation is a major determinant of useful trait diversity. We describe how genome analysis methods are enabling discovery of trait-associated structural variants and their potential impact on breeding. As our understanding of complex crop genomes continues to grow, there is growing evidence that structural genome variation plays a major role in determining traits important for breeding and agriculture. Identifying the extent and impact of structural variants in crop genomes is becoming increasingly feasible with ongoing advances in the sophistication of genome sequencing technologies, particularly as it becomes easier to generate accurate long sequence reads on a genome-wide scale. In this article, we discuss the origins of structural genome variation in crops from ancient and recent genome duplication and polyploidization events and review high-throughput methods to assay such variants in crop populations in order to find associations with phenotypic traits. There is increasing evidence from such studies that gene presence-absence and copy number variation resulting from segmental chromosome exchanges may be at the heart of adaptive variation of crops to counter abiotic and biotic stress factors. We present examples from major crops that demonstrate the potential of pangenomic diversity as a key resource for future plant breeding for resilience and sustainability.

Natural variation of ZmHKT1 affects root morphology in maize at the seedling stage

Eight variants in ZmHKT1 promoter were significantly associated with root diameter, four haplotypes based on these significant variants were found, and Hap2 has the largest root diameter. Roots play an important role in uptake of water, nutrients and plant anchorage. Identification of gene and corresponding SNPs associated with root traits would enable develop maize lines with better root traits that might help to improve capacity for absorbing nutrients and water acquisition. The genomic sequences of a salt tolerance gene ZmHKT1 was resequenced in 349 maize inbred lines, and the association between nucleotide polymorphisms and seedling root traits was detected. A total of 269 variants in ZmHKT1 were identified, including 226 single nucleotide polymorphisms and 43 insertions and deletions. The gene displayed high level of nucleotide diversity, especially in non-genic regions. A total of 19 variations in untranslated region of ZmHKT1 were found to be associated with six seedling traits. Eight variants in promoter region were significantly associated with average root diameter (ARD), four haplotypes were found based on these significant variants, and Hap2 has the largest ARD. Two SNPs in high-linkage disequilibrium (SNP-415 and SNP 2169) with pleiotropic effects were significantly associated with plant height, root surface area, root volume, and shoot dry weight. This result revealed that ZmHKT1 was an important contributor to the phenotypic variations of seedling root traits in maize, these significant variants could use to develop functional markers to improve root traits.

Spatio-Temporal Transcriptional Dynamics of Maize Long Non-Coding RNAs Responsive to Drought Stress

Long non-coding RNAs (lncRNAs) have emerged as important regulators in plant stress response. Here, we report a genome-wide lncRNA transcriptional analysis in response to drought stress using an expanded series of maize samples collected from three distinct tissues spanning four developmental stages. In total, 3488 high-confidence lncRNAs were identified, among which 1535 were characterized as drought responsive. By characterizing the genomic structure and expression pattern, we found that lncRNA structures were less complex than protein-coding genes, showing shorter transcripts and fewer exons. Moreover, drought-responsive lncRNAs exhibited higher tissue- and development-specificity than protein-coding genes. By exploring the temporal expression patterns of drought-responsive lncRNAs at different developmental stages, we discovered that the reproductive stage R1 was the most sensitive growth stage with more lncRNAs showing altered expression upon drought stress. Furthermore, lncRNA target prediction revealed 653 potential lncRNA-messenger RNA (mRNA) pairs, among which 124 pairs function in cis-acting mode and 529 in trans. Functional enrichment analysis showed that the targets were significantly enriched in molecular functions related to oxidoreductase activity, water binding, and electron carrier activity. Multiple promising targets of drought-responsive lncRNAs were discovered, including the V-ATPase encoding gene, vpp4. These findings extend our knowledge of lncRNAs as important regulators in maize drought response.

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.

Crop genome-wide association study: a harvest of biological relevance

With the advent of rapid genotyping and next-generation sequencing technologies, genome-wide association study (GWAS) has become a routine strategy for decoding genotype-phenotype associations in many species. More than 1000 such studies over the last decade have revealed substantial genotype-phenotype associations in crops and provided unparalleled opportunities to probe functional genomics. Beyond the many 'hits' obtained, this review summarizes recent efforts to increase our understanding of the genetic architecture of complex traits by focusing on non-main effects including epistasis, pleiotropy, and phenotypic plasticity. We also discuss how these achievements and the remaining gaps in our knowledge will guide future studies. Synthetic association is highlighted as leading to false causality, which is prevalent but largely underestimated. Furthermore, validation evidence is appealing for future GWAS, especially in the context of emerging genome-editing technologies.

PbrmiR397a regulates lignification during stone cell development in pear fruit

Lignified stone cells substantially reduce fruit quality. Therefore, it is desirable to inhibit stone cell development using genetic technologies. However, the molecular mechanisms regulating lignification are poorly understood in fruit stone cells. In this study, we have shown that microRNA (miR) miR397a regulates fruit cell lignification by inhibiting laccase (LAC) genes that encode key lignin biosynthesis enzymes. Transient overexpression of PbrmiR397a, which is the miR397a of Chinese pear (Pyrus bretschneideri), and simultaneous silencing of three LAC genes reduced the lignin content and stone cell number in pear fruit. A single nucleotide polymorphism (SNP) identified in the promoter of the PbrmiR397a gene was found to associate with low levels of fruit lignin, after analysis of the genome sequences of sixty pear varieties. This SNP created a TCA element that responded to salicylic acid to induce gene expression as confirmed using a cell-based assay system. Furthermore, stable overexpression of PbrmiR397a in transgenic tobacco plants reduced the expression of target LAC genes and decreased the content of lignin but did not change the ratio of syringyl- and guaiacyl-lignin monomers. Consistent with reduction in lignin content, the transgenic plants showed fewer numbers of vessel elements and thinner secondary walls in the remaining elements compared to wild-type control plants. This study has advanced our understanding of the regulation of lignin biosynthesis and provided useful molecular genetic information for improving pear fruit quality.

Improving Genomic Selection With Quantitative Trait Loci and Nonadditive Effects Revealed by Empirical Evidence in Maize

Genomic selection (GS), a tool developed for molecular breeding, is used by plant breeders to improve breeding efficacy by shortening the breeding cycle and to facilitate the selection of candidate lines for creating hybrids without phenotyping in various environments. Association and linkage mapping have been widely used to explore and detect candidate genes in order to understand the genetic mechanisms of quantitative traits. In the current study, phenotypic and genotypic data from three experimental populations, including data on six agronomic traits (e.g., plant height, ear height, ear length, ear diameter, grain yield per plant, and hundred-kernel weight), were used to evaluate the effect of trait-relevant markers (TRMs) on prediction accuracy estimation. Integrating information from mapping into a statistical model can efficiently improve prediction performance compared with using stochastically selected markers to perform GS. The prediction accuracy can reach plateau when a total of 500-1,000 TRMs are utilized in GS. The prediction accuracy can be significantly enhanced by including nonadditive effects and TRMs in the GS model when genotypic data with high proportions of heterozygous alleles and complex agronomic traits with high proportion of nonadditive variancein phenotypic variance are used to perform GS. In addition, taking information on population structure into account can slightly improve prediction performance when the genetic relationship between the training and testing sets is influenced by population stratification due to different allele frequencies. In conclusion, GS is a useful approach for prescreening candidate lines, and the empirical evidence provided by the current study for TRMs and nonadditive effects can inform plant breeding and in turn contribute to the improvement of selection efficiency in practical GS-assisted breeding programs.

Sodium-Related Adaptations to Drought: New Insights From the Xerophyte Plant Zygophyllum xanthoxylum

Understanding the unusual physiological mechanisms that enable drought tolerance in xerophytes will be of considerable benefit because of the potential to identify novel and key genetic elements for future crop improvements. These plants are interesting because they are well-adapted for life in arid zones; Zygophyllum xanthoxylum, for example, is a typical xerophytic shrub that inhabits central Asian deserts, accumulating substantial levels of sodium (Na+) in its succulent leaves while growing in soils that contain very low levels of this ion. The physiological importance of this unusual trait to drought adaptations remains poorly understood, however. Thus, 2-week-old Z. xanthoxylum plants were treated with 50 mM NaCl (Na) for 7 days in this study in order to investigate their drought tolerance, leaf osmotic potential (Ψs) related parameters, anatomical characteristics, and transpiration traits. The results demonstrated that NaCl treatment significantly enhanced both the survivability and durability of Z. xanthoxylum plants under extreme drought conditions. The bulk of the Na+ ions encapsulated in plants was overwhelmingly allocated to leaves rather than roots or stems under drought conditions; thus, compared to the control, significantly more Na+ compared to other solutes such as K+, Ca2+, Cl-, sugars, and proline accumulated in the leaves of NaCl-treated plants and led to a marked decrease (31%) in leaf Ψs. In addition, the accumulation of Na+ ions also resulted in mesophyll cell enlargement and leaf succulence, enabling the additional storage of water; Na+ ions also reduced the rate of water loss by decreasing stomatal density and down-regulating stomatal aperture size. The results of this study demonstrate that Z. xanthoxylum has evolved a notable ability to utilize Na+ ions to lower Ψs, swell its leaves, and decrease stomatal aperture sizes, in order to enable the additional uptake and storage of water and mitigate losses. These distinctive drought adaption characteristics mean that the xerophytic plant Z. xanthoxylum presents a fascinating case study for the potential identification of important and novel genetic elements that could improve crops. This report provides insights on the eco-physiological role of sodium accumulation in xerophytes adapted to extremely arid habitats.

Vacuolar Proton Pyrophosphatase Is Required for High Magnesium Tolerance in Arabidopsis

Magnesium (Mg²⁺) is an essential nutrient in all organisms. However, high levels of Mg²⁺ in the environment are toxic to plants. In this study, we identified the vacuolar-type H⁺-pyrophosphatase, AVP1, as a critical enzyme for optimal plant growth under high-Mg conditions. The Arabidopsis avp1 mutants displayed severe growth retardation, as compared to the wild-type plants upon excessive Mg²⁺. Unexpectedly, the avp1 mutant plants retained similar Mg content to wild-type plants under either normal or high Mg conditions, suggesting that AVP1 may not directly contribute to Mg²⁺ homeostasis in plant cells. Further analyses confirmed that the avp1 mutant plants contained a higher pyrophosphate (PPi) content than wild type, coupled with impaired vacuolar H⁺-pyrophosphatase activity. Interestingly, expression of the Saccharomyces cerevisiae cytosolic inorganic pyrophosphatase1 gene IPP1, which facilitates PPi hydrolysis but not proton translocation into vacuole, rescued the growth defects of avp1 mutants under high-Mg conditions. These results provide evidence that high-Mg sensitivity in avp1 mutants possibly resulted from elevated level of cytosolic PPi. Moreover, genetic analysis indicated that mutation of AVP1 was additive to the defects in mgt6 and cbl2 cbl3 mutants that are previously known to be impaired in Mg²⁺ homeostasis. Taken together, our results suggest AVP1 is required for cellular PPi homeostasis that in turn contributes to high-Mg tolerance in plant cells.

Dissecting the genetic architecture of waterlogging stress-related traits uncovers a key waterlogging tolerance gene in maize

A key candidate gene, GRMZM2G110141, which could be used in marker-assisted selection in maize breeding programs, was detected among the 16 genetic loci associated with waterlogging tolerance identified through genome-wide association study. Waterlogging stress seriously affects the growth and development of upland crops such as maize (Zea mays L.). However, the genetic basis of waterlogging tolerance in crop plants is largely unknown. Here, we identified genetic loci for waterlogging tolerance-related traits by conducting a genome-wide association study using maize phenotypes evaluated in the greenhouse under waterlogging stress and normal conditions. A total of 110 trait-single nucleotide polymorphism associations spanning 16 genomic regions were identified; single associations explained 2.88-10.67% of the phenotypic variance. Among the genomic regions identified, 14 co-localized with previously detected waterlogging tolerance-related quantitative trail loci. Furthermore, 33 candidate genes involved in a wide range of stress-response pathways were predicted. We resequenced a key candidate gene (GRMZM2G110141) in 138 randomly selected inbred lines and found that variations in the 5'-UTR and in the mRNA abundance of this gene under waterlogging conditions were significantly associated with leaf injury. Furthermore, we detected favorable alleles of this gene and validated the favorable alleles in two different recombinant inbred line populations. These alleles enhanced waterlogging tolerance in segregating populations, strongly suggesting that GRMZM2G110141 is a key waterlogging tolerance gene. The set of waterlogging tolerance-related genomic regions and associated markers identified here could be valuable for isolating waterlogging tolerance genes and improving this trait in maize.

Variation in the regulatory region of FZP causes increases in secondary inflorescence branching and grain yield in rice domestication

Inflorescence branching is a key agronomic trait determining rice yield. The primary branch of the ancestral wild rice (Oryza rufipogon Griff.) bears few grains, due to minimal secondary branching. By contrast, Oryza sativa cultivars have been selected to produce large panicles with more secondary branches. Here we showed that the CONTROL OF SECONDARY BRANCH 1 (COS1) gene, which is identical to FRIZZY PANICLE (FZP), plays an important role in the key transition from few secondary branches in wild rice to more secondary branches in domesticated rice cultivars. A 4-bp tandem repeat deletion approximately 2.7 kb upstream of FZP may affect the binding activities of auxin response factors to the FZP promoter, decrease the expression level of FZP and significantly enhance the number of secondary branches and grain yield in cultivated rice. Functional analyses showed that NARROW LEAF 1 (NAL1), a trypsin-like serine and cysteine protease, interacted with FZP and promoted its degradation. Consistently, downregulating FZP expression or upregulating NAL1 expression in the commercial cultivar Zhonghua 17 increased the number of secondary branches per panicle, grain number per panicle and grain yield per plant. Our findings not only provide insights into the molecular mechanism of increasing grain number and yield during rice domestication, but also offer favorable genes for improving the grain yield of rice.

Pyrophosphate inhibits gluconeogenesis by restricting UDP-glucose formation in vivo

Pyrophosphate (PPi) is produced by anabolic reactions and serves as an energy donor in the cytosol of plant cells; however, its accumulation to toxic levels disrupts several common biosynthetic pathways and is lethal. Before acquiring photosynthetic capacity, young seedlings must endure a short but critical heterotrophic period, during which they are nourished solely by sugar produced from seed reserves by the anabolic process of gluconeogenesis. Previously, we reported that excess PPi in H⁺-PPase-knockout fugu5 mutants of Arabidopsis thaliana severely compromised gluconeogenesis. However, the precise metabolic target of PPi inhibition in vivo remained elusive. Here, CE-TOF MS analyses of major metabolites characteristic of gluconeogenesis from seed lipids showed that the Glc6P;Fru6P level significantly increased and that Glc1P level was consistently somewhat higher in fugu5 compared to wild type. In contrast, the UDP-Glc level decreased significantly in the mutants. Importantly, specific removal of PPi in fugu5, and thus in AVP1_pro:IPP1 transgenic lines, restored the Glc1P and the Glc6P;Fru6P levels, increased the UDP-Glc level ~2.0-fold, and subsequently increased sucrose synthesis. Given the reversible nature of the Glc1P/UDP-Glc reaction, our results indicate that UGP-Glc pyrophosphorylase is the major target when excess PPi exerts inhibitory effects in vivo. To validate our findings, we analyzed metabolite responses using a mathematical theory called structural sensitivity analysis (SSA), in which the responses of concentrations in reaction systems to perturbations in enzyme activity are determined from the structure of the network alone. A comparison of our experimental data with the results of pure structural theory predicted the existence of unknown reactions as the necessary condition for the above metabolic profiles, and confirmed the above results. Our data support the notion that H⁺-PPase plays a pivotal role in cytosolic PPi homeostasis in plant cells. We propose that the combination of metabolomics and SSA is powerful when seeking to identify and predict metabolic targets in living cells.

Genome-Wide Association Studies of Photosynthetic Traits Related to Phosphorus Efficiency in Soybean

Photosynthesis is the basis of plant growth and development, and is seriously affected by low phosphorus (P) stress. However, few studies have reported for the genetic foundation of photosynthetic response to low P stress in soybean. To address this issue, 219 soybean accessions were genotyped by 292,035 high-quality single nucleotide polymorphisms (SNPs) and phenotyped under normal and low P conditions in 2015 and 2016. These datasets were used to identify quantitative trait nucleotides (QTNs) for photosynthesis-related traits using mrMLM, ISIS EM-BLASSO, pLARmEB, FASTmrMLM, FASTmrEMMA, and pKWmEB methods. As a result, 159 QTNs within 31 genomic regions were found to be associated with four photosynthesis-related traits under different P stress conditions. Among the 31 associated regions, five (q7-2, q8-1, q9, q13-1, and q20-2) were detected commonly under both normal and low P conditions, indicating the insensitivity of these candidate genes to low P stress; five were detected only under normal P condition, indicating the sensitivity of these candidate genes to low P stress; six were detected only under low P condition, indicating the tolerantness of these candidate genes to low P stress; 20 were reported in previous studies. Around the 159 QTNs, 52 candidate genes were mined. These results provide the important information for marker-assisted breeding in soybean and further reveal the basis for the application of P tolerance to photosynthetic capacity.

Genetic characterization of inbred lines from Shaan A and B groups for identifying loci associated with maize grain yield

BACKGROUND

Increasing grain yield is a primary objective of maize breeding. Dissecting the genetic architecture of grain yield furthers genetic improvements to increase yield. Presented here is an association panel composed of 126 maize inbreds (AM126), which were genotyped by the genotyping-by-sequencing (tGBS) method. We performed genetic characterization and association analysis related to grain yield in the association panel.

RESULTS

In total, 46,046 SNPs with a minor allele frequency (MAF) ≥0.01 were used to assess genetic diversity and kinship in AM126. The results showed that the average MAF and polymorphism information content (PIC) were 0.164 and 0.198, respectively. The Shaan B group, with 11,284 unique SNPs, exhibited greater genetic diversity than did the Shaan A group, with 2644 SNPs. The 61.82% kinship coefficient in AM126 was equal to 0, and only 0.15% of that percentage was greater than 0.7. A total of 31,983 SNPs with MAF ≥0.05 were used to characterize population structure, LD decay and association mapping. Population structure analysis suggested that AM126 can be divided into 6 subgroups, which is consistent with breeding experience and pedigree information. The LD decay distance in AM126 was 150 kb. A total of 51 significant SNPs associated with grain yield were identified at P < 1 × 10- 3 across two environments (Yangling and Yulin). Among those SNPs, two loci displayed overlapping regions in the two environments. Finally, 30 candidate genes were found to be associated with grain yield.

CONCLUSIONS

These results contribute to the genetic characterization of this breeding population, which serves as a reference for hybrid breeding and population improvement, and demonstrate the genetic architecture of maize grain yield, potentially facilitating genetic improvement.

Genomic variations and association study of agronomic traits in flax

BACKGROUND

Flax (Linum usitatissimum. L) is an ancient oilseed and natural fiber crop. It could be divided into three categories by use, namely oil flax, fiber flax and oil-fiber dual purpose (OF). Cultivated flax is widely used in the food and textile industry. It is of great significance to elucidate the genetic characteristics of flax collections for accelerating the process of breeding improvement in this dual purpose crop. With the development of next-generation sequencing, we can use new methods, such as SLAF-seq (specific-locus amplified fragment sequencing), to decode unknown genomes of species. In this study, a high-through sequencing of flax collections using SLAF-seq was conducted. The evolutionary tendency was defined and candidate genes associated with agronomic traits of flax species were identified by Genome-Wide Association Studying (GWAS).

RESULTS

A flax collection consisting of 224 varieties were sequenced by SLAF-seq. In total, 346,639 SLAF tags were developed from all accessions, with an average sequencing depth of 7.19 for each accession. A total of 584,987 SNPs (single nucleotide polymorphism) with an MAF > 0.05 were identified from these SLAFs. The population structure division and phylogenetic analysis indicated a strong divergence among three kinds of flax groups. The genome-wide variation uncovered that oil flax had the highest genetic diversity and was considered to be the ancestor of fiber flax and oil-fiber flax. Sixteen associated peak SNPs for six traits were obtained by GWAS of oil-related traits using EMMAX (efficient mixed-model association eXpedited). Candidate genes and their related pathway were evaluated. A new GWAS was developed for fiber properties using the GLM (General linear model) model and a number of loci were identified.

CONCLUSIONS

To our knowledge, this is the first study on discovery multiple loci for important agronomic traits of flax species using GWAS strategy. These results will provide the highest possibility of incorporating both high fiber and good oil traits in a single variety.

Genomic variations and association study of agronomic traits in flax

BACKGROUND

Flax (Linum usitatissimum. L) is an ancient oilseed and natural fiber crop. It could be divided into three categories by use, namely oil flax, fiber flax and oil-fiber dual purpose (OF). Cultivated flax is widely used in the food and textile industry. It is of great significance to elucidate the genetic characteristics of flax collections for accelerating the process of breeding improvement in this dual purpose crop. With the development of next-generation sequencing, we can use new methods, such as SLAF-seq (specific-locus amplified fragment sequencing), to decode unknown genomes of species. In this study, a high-through sequencing of flax collections using SLAF-seq was conducted. The evolutionary tendency was defined and candidate genes associated with agronomic traits of flax species were identified by Genome-Wide Association Studying (GWAS).

RESULTS

A flax collection consisting of 224 varieties were sequenced by SLAF-seq. In total, 346,639 SLAF tags were developed from all accessions, with an average sequencing depth of 7.19 for each accession. A total of 584,987 SNPs (single nucleotide polymorphism) with an MAF > 0.05 were identified from these SLAFs. The population structure division and phylogenetic analysis indicated a strong divergence among three kinds of flax groups. The genome-wide variation uncovered that oil flax had the highest genetic diversity and was considered to be the ancestor of fiber flax and oil-fiber flax. Sixteen associated peak SNPs for six traits were obtained by GWAS of oil-related traits using EMMAX (efficient mixed-model association eXpedited). Candidate genes and their related pathway were evaluated. A new GWAS was developed for fiber properties using the GLM (General linear model) model and a number of loci were identified.

CONCLUSIONS

To our knowledge, this is the first study on discovery multiple loci for important agronomic traits of flax species using GWAS strategy. These results will provide the highest possibility of incorporating both high fiber and good oil traits in a single variety.

Genome-Wide Association Mapping for Cold Tolerance in a Core Collection of Rice (Oryza sativa L.) Landraces by Using High-Density Single Nucleotide Polymorphism Markers From Specific-Locus Amplified Fragment Sequencing

Understanding the genetic mechanism of cold tolerance in rice is important to mine elite genes from rice landraces and breed excellent cultivars for this trait. In this study, a genome-wide association study (GWAS) was performed using high-density single nucleotide polymorphisms (SNPs) obtained using specific-locus amplified fragment sequencing (SLAF-seq) technology from a core collection of landraces of rice. A total of 67,511 SNPs obtained from 116,643 SLAF tags were used for genotyping the 150 accessions of rice landraces in the Ting's rice core collection. A compressed mixed liner model was used to perform GWAS by using the high-density SNPs for cold tolerance in rice landraces at the seedling stage. A total of 26 SNPs were found to be significantly (P < 1.48 × 10-7) associated with cold tolerance, which could explained phenotypic variations ranging from 26 to 33%. Among them, two quantitative trait loci (QTLs) were mapped closely to the previously cloned/mapped genes or QTLs for cold tolerance. A newly identified QTL for cold tolerance in rice was further characterized by sequencing, real time-polymerase chain reaction, and bioinformatics analyses. One candidate gene, i.e., Os01g0620100, showed different gene expression levels between the cold tolerant and sensitive landraces under cold stress. We found the difference of coding amino acid in Os01g0620100 between cold tolerant and sensitive landraces caused by polymorphism within the coding domain sequence. In addition, the prediction of Os01g0620100 protein revealed a WD40 domain that was frequently found in cold tolerant landraces. Therefore, we speculated that Os01g0620100 was highly important for the response to cold stress in rice. These results indicated that rice landraces are important sources for investigating rice cold tolerance, and the mapping results might provide important information to breed cold-tolerant rice cultivars by using marker-assisted selection.

Genome-Wide Association Studies of Image Traits Reveal Genetic Architecture of Drought Resistance in Rice

Understanding how plants respond to drought can benefit drought resistance (DR) breeding. Using a non-destructive phenotyping facility, 51 image-based traits (i-traits) for 507 rice accessions were extracted. These i-traits can be used to monitor drought responses and evaluate DR. High heritability and large variation of these traits was observed under drought stress in the natural population. A genome-wide association study (GWAS) of i-traits and traditional DR traits identified 470 association loci, some containing known DR-related genes. Of these 470 loci, 443 loci (94%) were identified using i-traits, 437 loci (93%) co-localized with previously reported DR-related quantitative trait loci, and 313 loci (66.6%) were reproducibly identified by GWAS in different years. Association networks, established based on GWAS results, revealed hub i-traits and hub loci. This demonstrates the feasibility and necessity of dissecting the complex DR trait into heritable and simple i-traits. As proof of principle, we illustrated the power of this integrated approach to identify previously unreported DR-related genes. OsPP15 was associated with a hub i-trait, and its role in DR was confirmed by genetic transformation experiments. Furthermore, i-traits can be used for DR linkage analyses, and 69 i-trait locus associations were identified by both GWAS and linkage analysis of a recombinant inbred line population. Finally, we confirmed the relevance of i-traits to DR in the field. Our study provides a promising novel approach for the genetic dissection and discovery of causal genes for DR.

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.

Genetic-based dissection of arsenic accumulation in maize using a genome-wide association analysis method

Understanding the mechanism of arsenic (As) accumulation in plants is important in reducing As's toxicity to plants and its potential risks to human health. Here, we performed a genome-wide association study to dissect the genetic basis of the As contents of different maize tissues in Xixian, which was irrigated with As-rich surface water, and Changge using an association population consisting of 230 representative maize inbred lines. Phenotypic data revealed a wide normal distribution and high repeatability for the As contents in maize tissues. The As concentrations in maize tissues followed the same trend in the two locations: kernels < axes < stems < bracts < leaves. In total, 15, 16 and 15 non-redundant quantitative trait loci (QTLs) associated with As concentrations were identified (P ≤ 2.04 × 10-6 ) in five tissues from Xixian, Changge, and the combination of the locations, respectively, explaining 9.70%-24.65% of the phenotypic variation for each QTL, on average. Additionally, four QTLs [involving 15 single nucleotide polymorphisms (SNPs)] were detected in the single and the combined locations, indicating that these loci/SNPs might be stable across different environments. The candidate genes associated with these four loci were predicted. In addition, four non-redundant QTLs (6 SNPs), including a QTL that was detected in multiple locations according to the genome-wide association study, were found to co-localize with four previously reported QTL intervals. These results are valuable to understand the genetic architecture of As mechanism in maize and facilitate the genetic improvement of varieties without As toxicity.

Maize Urb2 protein is required for kernel development and vegetative growth by affecting pre-ribosomal RNA processing

Ribosome biogenesis is a fundamental process in eukaryotic cells. Although Urb2 protein has been implicated in ribosome biogenesis in yeast, the Urb2 domain is loosely conserved between plants and yeast, and the function of Urb2 protein in plants remains unknown. Here, we isolated a maize mutant, designated as urb2, with defects in kernel development and vegetative growth. Positional cloning and transgenic analysis revealed that urb2 encodes an Urb2 domain-containing protein. Compared with the wild-type (WT), the urb2 mutant showed decreased ratios of 60S/40S and 80S/40S and increased ratios of polyribosomes. The pre-rRNA intermediates of 35/33S rRNA, P-A3 and 18S-A3 were significantly accumulated in the urb2 mutant. Transcriptome profiling of the urb2 mutant indicated that ZmUrb2 affects the expression of a number of ribosome-related genes. We further demonstrated that natural variations in ZmUrb2 are significantly associated with maize kernel length. The overall results indicate that, by affecting pre-rRNA processing, the Urb2 protein is required for ribosome biogenesis in maize.

Effects of drought stress and water recovery on physiological responses and gene expression in maize seedlings

BACKGROUND

Drought is one of the major factors limiting global maize production. Exposure to long-term drought conditions inhibits growth and leads to yield losses. Although several drought-responsive genes have been identified and functionally analyzed, the mechanisms underlying responses to drought and water recovery treatments have not been fully elucidated. To characterize how maize seedling respond to drought stress at the transcriptional level, we analyzed physiological responses and differentially expressed genes (DEGs) in the inbred line B73 under water deficit and recovery conditions.

RESULTS

The data for relative leaf water content, leaf size, and photosynthesis-related parameters indicated that drought stress significantly repressed maize seedling growth. Further RNA sequencing analysis revealed that 6107 DEGs were responsive to drought stress and water recovery, with more down-regulated than up-regulated genes. Among the DEGs, the photosynthesis- and hormone-related genes were enriched in responses to drought stress and re-watering. Additionally, transcription factor genes from 37 families were differentially expressed among the three analyzed time-points. Gene ontology enrichment analyses of the DEGs indicated that 50 GO terms, including those related to photosynthesis, carbohydrate metabolism, oxidoreductase activities, nutrient metabolism and other drought-responsive pathways, were over-represented in the drought-treated seedlings. The content of gibberellin in drought treatment seedlings was decreased compared to that of control seedlings, while abscisic acid showed accumulated in the drought treated plants. The deep analysis of DEGs related to cell wall development indicated that these genes were prone to be down-regulated at drought treatment stage.

CONCLUSIONS

Many genes that are differentially expressed in responses to drought stress and water recovery conditions affect photosynthetic systems and hormone biosynthesis. The identified DEGs, especially those encoding transcription factors, represent potential targets for developing drought-tolerant maize lines.

Genomics-Enabled Next-Generation Breeding Approaches for Developing System-Specific Drought Tolerant Hybrids in Maize

Breeding science has immensely contributed to the global food security. Several varieties and hybrids in different food crops including maize have been released through conventional breeding. The ever growing population, decreasing agricultural land, lowering water table, changing climate, and other variables pose tremendous challenge to the researchers to improve the production and productivity of food crops. Drought is one of the major problems to sustain and improve the productivity of food crops including maize in tropical and subtropical production systems. With advent of novel genomics and breeding tools, the way of doing breeding has been tremendously changed in the last two decades. Drought tolerance is a combination of several component traits with a quantitative mode of inheritance. Rapid DNA and RNA sequencing tools and high-throughput SNP genotyping techniques, trait mapping, functional characterization, genomic selection, rapid generation advancement, and other tools are now available to understand the genetics of drought tolerance and to accelerate the breeding cycle. Informatics play complementary role by managing the big-data generated from the large-scale genomics and breeding experiments. Genome editing is the latest technique to alter specific genes to improve the trait expression. Integration of novel genomics, next-generation breeding, and informatics tools will accelerate the stress breeding process and increase the genetic gain under different production systems.

GWAS Uncovers Differential Genetic Bases for Drought and Salt Tolerances in Sesame at the Germination Stage

Sesame has great potential as an industrial crop but its production is challenged by drought and salt stresses. To unravel the genetic variants leading to salinity and drought tolerances at the germination stage, genome-wide association studies of stress tolerance indexes related to NaCl-salt and polyethylene glycol-drought induced stresses were performed with a diversity panel of 490 sesame accessions. An extensive variation was observed for drought and salt responses in the population and most of the accessions were moderately tolerant to both stresses. A total of 132 and 120 significant Single Nucleotide Polymorphisms (SNPs) resolved to nine and 15 Quantitative trait loci (QTLs) were detected for drought and salt stresses, respectively. Only two common QTLs for drought and salt responses were found located on linkage groups 5 and 7, respectively. This indicates that the genetic bases for drought and salt responses in sesame are different. A total of 13 and 27 potential candidate genes were uncovered for drought and salt tolerance indexes, respectively, encoding transcription factors, antioxidative enzymes, osmoprotectants and involved in hormonal biosynthesis, signal transduction or ion sequestration. The identified SNPs and potential candidate genes represent valuable resources for future functional characterization towards the enhancement of sesame cultivars for drought and salt tolerances.

The genetic architecture of amylose biosynthesis in maize kernel

Starch is the most abundant storage carbohydrate in maize kernel. The content of amylose and amylopectin confers unique properties in food processing and industrial application. Thus, the resurgent interest has been switched to the study of individual amylose or amylopectin rather than total starch, whereas the enzymatic machinery for amylose synthesis remains elusive. We took advantage of the phenotype of amylose content and the genotype of 9,007,194 single nucleotide polymorphisms from 464 inbred maize lines. The genome-wide association study identified 27 associated loci involving 39 candidate genes that were linked to amylose content including transcription factors, glycosyltransferases, glycosidases, as well as hydrolases. Except the waxy gene that encodes the granule-bound starch synthase, the remaining candidate genes were located in the upstream pathway of amylose synthesis, while the downstream members were already known from prior studies. The linked candidate genes could be transferred to manipulate amylose content and thus add value to maize kernel in the breeding programme.

A retrotransposon in an HKT1 family sodium transporter causes variation of leaf Na+ exclusion and salt tolerance in maize

Soil salinity is one of several major abiotic stresses that constrain maize productivity worldwide. An improved understanding of salt-tolerance mechanisms will thus enhance the breeding of salt-tolerant maize and boost productivity. Previous studies have indicated that the maintenance of leaf Na⁺ concentration is essential for maize salt tolerance, and the difference in leaf Na⁺ exclusion has previously been associated with variation in salt tolerance between maize varieties. Here, we report the identification and functional characterization of a maize salt-tolerance quantitative trait locus (QTL), Zea mays Na⁺ Content1 (ZmNC1), which encodes an HKT-type transporter (designated as ZmHKT1). We show that a natural ZmHKT1 loss-of-function allele containing a retrotransposon insertion confers increased accumulation of Na⁺ in leaves, and salt hypersensitivity. We next show that ZmHKT1 encodes a plasma membrane-localized Na⁺ -selective transporter, and is preferentially expressed in root stele (including the parenchyma cells surrounding the xylem vessels). We also show that loss of ZmHKT1 function increases xylem sap Na⁺ concentration and causes increased root-to-shoot Na⁺ delivery, indicating that ZmHKT1 promotes leaf Na⁺ exclusion and salt tolerance by withdrawing Na⁺ from the xylem sap. We conclude that ZmHKT1 is a major salt-tolerance QTL and identifies an important new gene target in breeding for improved maize salt tolerance.

Circular RNAs mediated by transposons are associated with transcriptomic and phenotypic variation in maize

Circular RNAs (circRNAs) are covalently closed RNA molecules. Recent studies have shown that circRNAs can arise from the transcripts of transposons. Given the prevalence of transposons in the maize genome and dramatic genomic variation driven by transposons, we hypothesize that transposons in maize may be involved in the formation of circRNAs and further modulate phenotypic variation. We performed circRNA-Seq on B73 seedling leaves and uncovered 2804 high-confidence maize circRNAs, which show distinct genomic features. Comprehensive analyses demonstrated that sequences related to LINE1-like elements (LLEs) and their Reverse Complementary Pairs (LLERCPs) are significantly enriched in the flanking regions of circRNAs. Interestingly, as the number of LLERCPs increase, the accumulation of circRNAs varies, whereas that of linear transcripts decreases. Furthermore, genes with LLERCP-mediated circRNAs are enriched among loci that are associated with phenotypic variation. These results suggest that circRNAs are likely to be involved in the modulation of phenotypic variation by LLERCPs. Further, we showed that the presence/absence variation of LLERCPs was associated with expression variation of circRNA-circ1690 and was related to ear height, potentially through the interplay between circRNAs and functional linear transcripts. Our first study of maize circRNAs uncovers a potential new way for transposons to modulate transcriptomic and phenotypic variations.

A Genome-Wide Association Study of Wheat Spike Related Traits in China

Rapid detection of allelic variation and identification of advantage haplotypes responsible for spike related traits play a crucial role in wheat yield improvement. The released genome sequence of hexaploid wheat (Chinese Spring) provides an extraordinary opportunity for rapid detection of natural variation and promotes breeding application. Here, selection signals detection and genome-wide association study (GWAS) were conducted for spike related traits. Based on the genotyping results by 90K SNP chip, 192 common wheat samples from southwest China were analyzed. One hundred and forty-six selective windows and one hundred and eighty-four significant SNPs (51 for spike length, 28 for kernels per spike, 39 for spikelet number, 30 for thousand kernel weight, and 36 for spike number per plant) were detected. Furthermore, tightly linkage and environmental stability window clusters and SNP clusters were also obtained. As a result, four SNP clusters associated with spike length were detected on chromosome 2A, 2B, 2D, and 6A. Two SNP clusters correlated to kernels per spike were detected on 2A and 2B. One pleiotropy SNP cluster correlated to spikelet number and kernels per spike was detected on 7B. According to the genome sequence, these SNP clusters and their overlapped/flanking QTLs which have been reported previously were integrated to a physical map. The candidate genes responsible for spike length, kernels per spike and spikelet number were predicted. Based on the genotypes of cultivars in south China, two advantage haplotypes associated with spike length and one advantage haplotype associated with kernels per spike/spikelet number were detected which have not been effectively transited into cultivars. According to these haplotypes, KASP markers were developed and diagnosed across landraces and cultivars which were selected from south and north China. Consequently, KASP assay, consistent with the GWAS results, provides reliable haplotypes for MAS in wheat yield improvement.

The Maize ABA Receptors ZmPYL8, 9, and 12 Facilitate Plant Drought Resistance

Drought is one of the major abiotic stresses affecting world agriculture. Breeding drought-resistant crops is one of the most important challenges for plant biologists. PYR1/PYL/RCARs, which encode the abscisic acid (ABA) receptors, play pivotal roles in ABA signaling, but how these genes function in crop drought response remains largely unknown. Here we identified 13 PYL family members in maize (ZmPYL1-13). Changes in expression of these genes under different stresses indicated that ZmPYLs played important roles in responding to multiple abiotic stresses. Transgenic analyses of ZmPYL genes in Arabidopsis showed that overexpression of ZmPYL3, ZmPYL9, ZmPYL10, and ZmPYL13 significantly enhanced the sensitivity of transgenic plants to ABA. Additionally, transgenic lines overexpressing ZmPYL8, ZmPYL9, and ZmPYL12 were more resistant to drought. Accumulation of proline and enhanced expression of drought-related marker genes in transgenic lines further confirmed the positive roles of ZmPYL genes in plant drought resistance. Association analyses with a panel of 368 maize inbred lines identified natural variants in ZmPYL8 and ZmPYL12 that were significantly associated with maize drought resistance. Our results deepen the knowledge of the function of maize PYL genes in responses to abiotic stresses, and the natural variants identified in ZmPYL genes may serve as potential molecular markers for breeding drought-resistant maize cultivars.

Association Analysis and Identification of ZmHKT1;5 Variation With Salt-Stress Tolerance

The high-affinity potassium transporter (HKT) genes are essential for plant salt stress tolerance. However, there were limited studies on HKTs in maize (Zea mays), and it is basically unknown whether natural sequence variations in these genes are associated with the phenotypic variability of salt tolerance. Here, the characterization of ZmHKT1;5 was reported. Under salt stress, ZmHKT1;5 expression increased strongly in salt-tolerant inbred lines, which accompanied a better-balanced Na⁺/K⁺ ratio and preferable plant growth. The association between sequence variations in ZmHKT1;5 and salt tolerance was evaluated in a diverse population comprising 54 maize varieties from different maize production regions of China. Two SNPs (A134G and A511G) in the coding region of ZmHKT1;5 were significantly associated with different salt tolerance levels in maize varieties. In addition, the favorable allele of ZmHKT1; 5 identified in salt tolerant maize varieties effectively endowed plant salt tolerance. Transgenic tobacco plants of overexpressing the favorable allele displayed enhanced tolerance to salt stress better than overexpressing the wild type ZmHKT1;5. Our research showed that ZmHKT1;5 expression could effectively enhance salt tolerance by maintaining an optimal Na⁺/K⁺ balance and increasing the antioxidant activity that keeps reactive oxygen species (ROS) at a low accumulation level. Especially, the two SNPs in ZmHKT1;5 might be related with new amino acid residues to confer salt tolerance in maize. Key Message: Two SNPs of ZmHKT1;5 related with salt tolerance were identified by association analysis. Overexpressing ZmHKT1;5 in tobaccos showed that the SNPs might enhance its ability to regulating Na⁺/K⁺ homeostasis.

Finger Millet [Eleusine coracana (L.) Gaertn.] Improvement: Current Status and Future Interventions of Whole Genome Sequence

The whole genome sequence (WGS) of the much awaited, nutrient rich and climate resilient crop, finger millet (Eleusine coracana (L.) Gaertn.) has been released recently. While possessing superior mineral nutrients and excellent shelf life as compared to other major cereals, multiploidy nature of the genome and relatively small plantation acreage in less developed countries hampered the genome sequencing of finger millet, disposing it as one of the lastly sequenced genomes in cereals. The genomic information available for this crop is very little when compared to other major cereals like rice, maize and barley. As a result, only a limited number of genetic and genomic studies has been undertaken for the improvement of this crop. Finger millet is known especially for its superior calcium content, but the high-throughput studies are yet to be performed to understand the mechanisms behind calcium transport and grain filling. The WGS of finger millet is expected to help to understand this and other important molecular mechanisms in finger millet, which may be harnessed for the nutrient fortification of other cereals. In this review, we discuss various efforts made so far on the improvement of finger millet including genetic improvement, transcriptome analysis, mapping of quantitative trait loci (QTLs) for traits, etc. We also discuss the pitfalls of modern genetic studies and provide insights for accelerating the finger millet improvement with the interventions of WGS in near future. Advanced genetic and genomic studies aided by WGS may help to improve the finger millet, which will be helpful to strengthen the nutritional security in addition to food security in the developing countries of Asia and Africa.

Genetic Diversity and Genome-Wide Association Study of Major Ear Quantitative Traits Using High-Density SNPs in Maize

Kernel and ear traits are key components of grain yield in maize (Zea mays L.). Investigation of these traits would help to develop high-yield varieties in maize. Genome-wide association study (GWAS) uses the linkage disequilibrium (LD) in the whole genome to determine the genes affecting certain phenotype. In this study, five ear traits (kernel length and width, ear length and diameter, cob diameter) were investigated across multi-environments for 2 years. Combining with the genotype obtained from Maize SNP50 chip, genetic diversity and association mapping in a set of 292 inbred lines were performed. Results showed that maize lines were clustered into seven subgroups and a total of 20 SNPs were found to be associated with ear traits significantly (P < 3.95E-05). The candidate genes identified by GWAS mainly encoded ubiquitin-activation enzymes (GRMZM2G015287), carotenoid cleavage dioxygenase (GRMZM2G446858), MYB-CC type transfactor, and phosphate starvation response protein 3, and they were associated with kernel length (KL) and ear diameter (ED), respectively. Moreover, two novel genes corresponding to RNA processing and fructose metabolism were found. Further, the SNPs detected by GWAS were confirmed by meta-QTL analysis. These genes and SNPs identified in the study would offer essential information for yield-related genes clone and breeding program in maize.

Genome-wide association analysis of lead accumulation in maize

Large phenotypic variations in the lead (Pb) concentration were observed in grains and leaves of maize plants. A further understanding of inheritance of Pb accumulation may facilitate improvement of low-Pb-accumulating cultivars in maize. A genome-wide association study was conducted in a population of 269 maize accessions with 43,737 single-nucleotide polymorphisms (SNPs). The Pb concentrations in leaves and kernels of 269 accessions were collected in pot-culture and field experiments in years of 2015 and 2016. Significant differences in Pb accumulation were found among individuals under different environments. Using the structure and kinship model, a total of 21 SNPs significantly associated with the Pb accumulation were identified with P < 2.28 × 10^-5 and FDR < 0.05 in the pot-culture and field experiments across 2 years. Three SNPs on chromosome 4 had significant associations simultaneously with the Pb concentrations of kernels and leaves and were co-localized with the previously detected quantitative trait loci. Through ridge regression best linear unbiased prediction Pb accumulation in the association population, the prediction accuracies by cross validation were 0.18-0.59 and 0.17-0.64, depending on the k-fold and the size of the training population. The results are helpful for genetic improvement and genomic prediction of Pb accumulation in maize.

Genetic variants associated with the root system architecture of oilseed rape (Brassica napus L.) under contrasting phosphate supply

Breeding crops with ideal root system architecture for efficient absorption of phosphorus is an important strategy to reduce the use of phosphate fertilizers. To investigate genetic variants leading to changes in root system architecture, 405 oilseed rape cultivars were genotyped with a 60K Brassica Infinium SNP array in low and high P environments. A total of 285 single-nucleotide polymorphisms were associated with root system architecture traits at varying phosphorus levels. Nine single-nucleotide polymorphisms corroborate a previous linkage analysis of root system architecture quantitative trait loci in the BnaTNDH population. One peak single-nucleotide polymorphism region on A3 was associated with all root system architecture traits and co-localized with a quantitative trait locus for primary root length at low phosphorus. Two more single-nucleotide polymorphism peaks on A5 for root dry weight at low phosphorus were detected in both growth systems and co-localized with a quantitative trait locus for the same trait. The candidate genes identified on A3 form a haplotype ‘BnA3Hap’, that will be important for understanding the phosphorus/root system interaction and for the incorporation into Brassica napus breeding programs.

Genome-Wide Association Study Reveals Novel Loci for SC7 Resistance in a Soybean Mutant Panel

Soybean mosaic virus (SMV) is a member of Potyvirus genus that causes severe yield loss and destroys seed quality in soybean [Glycine max (L.) Merr.]. It is important to explore new resistance sources and discover new resistance loci to SMV, which will provide insights to improve breeding strategies for SMV resistance. Here, a genome-wide association study was conducted to accelerate molecular breeding for the improvement of resistance to SMV in soybean. A population of 165 soybean mutants derived from two soybean parents was used in this study. There were 104 SNPs identified significantly associated with resistance to SC7, some of which were located within previous reported quantitative trait loci. Three putative genes on chromosome 1, 9, and 12 were homologous to WRKY72, eEF1Bβ, and RLP9, which were involved in defense response to insect and disease in Arabidopsis. Moreover, the expression levels of these three genes changed in resistance and susceptible soybean accessions after SMV infection. These three putative genes may involve in the resistance to SC7 and be worthy to further research. Collectively, markers significantly associated with resistance to SC7 will be helpful in molecular marker-assisted selection for breeding resistant soybean accessions to SMV, and the candidate genes identified would advance the functional study of resistance to SMV in soybean.

Natural variation identifies genes affecting drought-induced abscisic acid accumulation in Arabidopsis thaliana

Accumulation of the stress hormone abscisic acid (ABA) in response to drought and low water-potential controls many downstream acclimation mechanisms. However, mechanisms controlling ABA accumulation itself are less known. There was a 10-fold range of variation in ABA levels among nearly 300 Arabidopsis thaliana accessions exposed to the same low water-potential severity. Genome-wide association analysis (GWAS) identified genomic regions containing clusters of ABA-associated SNPs. Candidate genes within these regions included few genes with known stress or ABA-related function. The GWAS data were used to guide reverse genetic analysis, which found effectors of ABA accumulation. These included plasma-membrane-localized signaling proteins such as receptor-like kinases, aspartic protease, a putative lipid-binding START domain protein, and other membrane proteins of unknown function as well as a RING U-box protein and possible effect of tonoplast transport on ABA accumulation. Putative loss-of-function polymorphisms within the START domain protein were associated with climate factors at accession sites of origin, indicating its potential involvement in drought adaptation. Overall, using ABA accumulation as a basis for a combined GWAS-reverse genetic strategy revealed the broad natural variation in low-water-potential-induced ABA accumulation and was successful in identifying genes that affect ABA levels and may act in upstream drought-related sensing and signaling mechanisms. ABA effector loci were identified even when each one was of incremental effect, consistent with control of ABA accumulation being distributed among the many branches of ABA metabolism or mediated by genes with partially redundant function.

A systems approach to a spatio-temporal understanding of the drought stress response in maize

Crops are often subjected to periods of drought stress during their life cycle. However, how stress response mechanisms contribute to the crosstalk between stress signaling pathways and developmental signaling pathways is still unknown. We built a gene co-expression network from a spatio-temporal transcriptomic map of the drought stress response in maize (Zea mays), profiled from three tissues and four developmental stages and characterized hub genes associated with duplication events, selection, and regulatory networks. Co-expression analysis grouped drought-response genes into ten modules, covering 844 highly connected genes (hub genes). Of these, 15.4% hub genes had diverged by whole-genome duplication events and 2.5% might then have been selected during natural domestication and artificial improvement processes, successively. We identified key transcription factor hubs in a transcriptional regulatory network, which may function as a crosstalk mechanism between drought stress and developmental signalling pathways in maize. Understanding the evolutionary biases that have evolved to enhance drought adaptation lays the foundation for further dissection of crosstalk between stress signalling pathways and developmental signalling pathways in maize, towards molecular design of new cultivars with desirable yield and greater stress tolerance.

Genome-Wide Association Study Reveals Natural Variations Contributing to Drought Resistance in Crops

Crops are often cultivated in regions where they will face environmental adversities; resulting in substantial yield loss which can ultimately lead to food and societal problems. Thus, significant efforts have been made to breed stress tolerant cultivars in an attempt to minimize these problems and to produce more stability with respect to crop yields across broad geographies. Since stress tolerance is a complex and multi-genic trait, advancements with classical breeding approaches have been challenging. On the other hand, molecular breeding, which is based on transgenics, marker-assisted selection and genome editing technologies; holds great promise to enable farmers to better cope with these challenges. However, identification of the key genetic components underlying the trait is critical and will serve as the foundation for future crop genetic improvement. Recently, genome-wide association studies have made significant contributions to facilitate the discovery of natural variation contributing to stress tolerance in crops. From these studies, the identified loci can serve as targets for genomic selection or editing to enable the molecular design of new cultivars. Here, we summarize research progress on this issue and focus on the genetic basis of drought tolerance as revealed by genome-wide association studies and quantitative trait loci mapping. Although many favorable loci have been identified, elucidation of their molecular mechanisms contributing to increased stress tolerance still remains a challenge. Thus, continuous efforts are still required to functionally dissect this complex trait through comprehensive approaches, such as system biological studies. It is expected that proper application of the acquired knowledge will enable the development of stress tolerant cultivars; allowing agricultural production to become more sustainable under dynamic environmental conditions.

Molecular diversity analysis, drought related marker-traits association mapping and discovery of excellent alleles for 100-day old plants by EST-SSRs in cassava germplasms (Manihot esculenta Cranz)

Cassava is the third largest food crop of the world and has strong ability of drought tolerance. In order to evaluate the molecular diversity and to discover novel alleles for drought tolerance in cassava germplasms, we examined a total of 107 abiotic stress related expressed sequence tags-simple sequence repeat (EST-SSR) markers in 134 cassava genotypes coming from planting regions worldwide and performed drought related marker-traits association mapping. As results, we successfully amplified 98 of 107 markers in 97 polymorphic loci and 279 alleles, with 2.87 alleles per locus, gene diversity of 0.48 and polymorphic information content (PIC) of 0.41 on average. The genetic coefficient between every two lines was 0.37 on average, ranging from 0.21 to 0.82. According to our population structure analysis, these samples could be divided into three sub-populations showing obvious gene flow between them. We also performed water stress experiments using 100-day old cassava plants in two years and calculated the drought tolerance coefficients (DTCs) and used them as phenotypes for marker-trait association mapping. We found that 53 markers were significantly associated with these drought-related traits, with a contribution rate for trait variation of 8.60% on average, ranging between 2.66 and 28.09%. Twenty-four of these 53 associated genes showed differential transcription or protein levels which were confirmed by qRT-PCR under drought stress when compared to the control conditions in cassava. Twelve of twenty-four genes were the same differential expression patterns in omics data and results of qRT-PCR. Out of 33 marker-traits combinations on 24 loci, 34 were positive and 53 negative alleles according to their phenotypic effects and we also obtained the typical materials which carried these elite alleles. We also found 23 positive average allele effects while 10 loci were negative according to their allele effects (AAEs). Our results on molecular diversity, locus association and differential expression under drought can prove beneficial to select excellent materials through marker assisted selection and for functional genes research in the future.

An R2R3-type MYB transcription factor, GmMYB29, regulates isoflavone biosynthesis in soybean

Isoflavones comprise a group of secondary metabolites produced almost exclusively by plants in the legume family, including soybean [Glycine max (L.) Merr.]. They play vital roles in plant defense and have many beneficial effects on human health. Isoflavone content is a complex quantitative trait controlled by multiple genes, and the genetic mechanisms underlying isoflavone biosynthesis remain largely unknown. Via a genome-wide association study (GWAS), we identified 28 single nucleotide polymorphisms (SNPs) that are significantly associated with isoflavone concentrations in soybean. One of these 28 SNPs was located in the 5'-untranslated region (5'-UTR) of an R2R3-type MYB transcription factor, GmMYB29, and this gene was thus selected as a candidate gene for further analyses. A subcellular localization study confirmed that GmMYB29 was located in the nucleus. Transient reporter gene assays demonstrated that GmMYB29 activated the IFS2 (isoflavone synthase 2) and CHS8 (chalcone synthase 8) gene promoters. Overexpression and RNAi-mediated silencing of GmMYB29 in soybean hairy roots resulted in increased and decreased isoflavone content, respectively. Moreover, a candidate-gene association analysis revealed that 11 natural GmMYB29 polymorphisms were significantly associated with isoflavone contents, and regulation of GmMYB29 expression could partially contribute to the observed phenotypic variation. Taken together, these results provide important genetic insights into the molecular mechanisms underlying isoflavone biosynthesis in soybean.