Natural variation at the soybean J locus improves adaptation to the tropics and enhances yield

Molecular mapping of the Cf-10 gene by combining SNP/InDel-index and linkage analysis in tomato (Solanum lycopersicum)

BACKGROUND

Leaf mold, one of the major diseases of tomato caused by Cladosporium fulvum (C. fulvum), can dramatically reduce the yield and cause multimillion dollar losses annually worldwide. Mapping the resistance genes (R genes) of C. fulvum and devising MAS based strategies for breeding new cultivars is an effective approach to improve the resistance in tomato. Up to now, many C. fulvum genes or QTLs have been mapped using different genetic materials, but few studies focused on Cf-10 gene positioning.

RESULTS

In this study, we investigated the genetic rules for Cf-10 and used a novel combinatorial strategy to rapidly map the Cf-10 gene. Initially, the performance of F1, F2 and BC1F1 individuals after infection, demonstrated that the resistance against C. fulvum was controlled by a single dominant gene. Two pools of resistant and susceptible individuals from F2 population were investigated, using mapping by sequencing approach and Cf-10 was found to be localized to 3.35 Mb and 3.74 Mb on chromosome 1, employing SNP/InDel index methods, respectively. After accounting for overlapping regions, these two algorithms yielded a total length of 3.29 Mb, narrowing down the target region. We further developed five serviceable KASP markers for this region based on sequencing data and conducted local QTL mapping using individuals from the F2 population, except for mapping by sequencing as mentioned above. Finally Cf-10 gene was mapped spanning a region of 790 kb, where only one gene (Solyc01g007130.3) was annotated as probable receptor protein kinase TMK1 with a LRR motif, a common R gene characteristic. The RT-qPCR analysis further confirmed the localization and the relative expression of Solyc01g007130.3 in Ontario 792 and was found to be significantly higher than that in Moneymaker at 9 dpi and 12 dpi, respectively.

CONCLUSION

This study proposed a novel combinatorial strategy by combining SNP-index, InDel-index analyses and local QTL mapping using KASP genotyping approach to rapidly map genes responsible for specific traits and provided a robust base for cloning the Cf-10 gene. Furthermore, these analyses suggest that Solyc01g007130.3 is a potential candidate to be regarded as Cf-10 gene.

Combining Simple Phenotyping and Photothermal Algorithm for the Prediction of Soybean Phenology: Application to a Range of Common Cultivars Grown in Europe

Developing new cropping strategies (very early sowing, crop expansion at higher latitudes, double cropping) to improve soybean production in Europe under climate change needs a good prediction of phenology under different temperature and photoperiod conditions. For that purpose, a simple phenology algorithm (SPA) was developed and parameterized for 10 contrasting soybean cultivars (maturity group 000 to II). Two experiments were carried out at INRA Toulouse (France) for parameterization: 1) Phenological monitoring of plants in pots on an outdoor platform with 6 planting dates. 2) Response of seed germination to temperature in controlled conditions. Multi-location field trials including 5 sites, 4 years, 2 sowing dates, and 10 cultivars were used to evaluate the SPA phenology predictions. Mean cardinal temperatures (minimum, optimum, and maximum) for germination were ca. 2, 30, and 40°C, respectively with significant differences among cultivars. The photoperiod sensitivity coefficient varied among cultivars when fixing Popt and Pcrt, optimal and critical photoperiods respectively, by maturity group. The parameterized algorithm showed an RMSE of less than 6 days for the prediction of crop cycle duration (i.e. cotyledons stage to physiological maturity) in the field trials including 75 data points. Flowering (R1 stage), and beginning of grain filling (R5 stage) dates were satisfactorily predicted with RMSEs of 8.2 and 9.4 days respectively. Because SPA can be also parameterized using data from field experiments, it can be useful as a plant selection tool across environments. The algorithm can be readily applied to species other than soybean, and its incorporation into cropping systems models would enhance the assessment of the performance of crop cultivars under climate change scenarios.

Positional Cloning of the Flowering Time QTL qFT12-1 Reveals the Link Between the Clock Related PRR Homolog With Photoperiodic Response in Soybeans

Flowering time and maturity are important agronomic traits for soybean cultivars to adapt to different latitudes and achieve maximal yield. Genetic studies on genes and quantitative trait loci (QTL) that control flowering time and maturity are extensive. In particular, the molecular bases of E1-E4, E6, E9, E10, and J have been deciphered. For a better understanding of regulation of flowering time gene networks, we need to understand if more molecular factors carrying different biological functions are also involved in the regulation of flowering time in soybeans. We developed a population derived from a cross between a landrace Jilincailihua (male) and a Chinese cultivar Chongnong16 (female). Both parents carry the same genotypes of E1e2E3HaE4 at E1, E2, E3, and E4 loci. Nighty-six individuals of the F₂ population were genotyped with Illumina SoySNP8k iSelect BeadChip. A total of 2,407 polymorphic single nucleotide polymorphism (SNP) markers were used to construct a genetic linkage map. One major QTL, qFT12-1, was mapped to an approximately 567-kB region on chromosome 12. Genotyping and phenotyping of recombinant plant whose recombination events were occurring within the QTL region allowed us to narrow down the QTL region to 56.4 kB, in which four genes were annotated. Allelism and association analysis indicated Glyma.12G073900, a PRR7 homolog, is the strongest candidate gene for qFT12-1. The findings of this study disclosed the possible involvement of circadian clock gene in flowering time regulation of soybeans.

Transcriptome-Enabled Network Inference Revealed the GmCOL1 Feed-Forward Loop and Its Roles in Photoperiodic Flowering of Soybean

Photoperiodic flowering, a plant response to seasonal photoperiod changes in the control of reproductive transition, is an important agronomic trait that has been a central target of crop domestication and modern breeding programs. However, our understanding about the molecular mechanisms of photoperiodic flowering regulation in crop species is lagging behind. To better understand the regulatory gene networks controlling photoperiodic flowering of soybeans, we elucidated global gene expression patterns under different photoperiod regimes using the near isogenic lines (NILs) of maturity loci (E loci). Transcriptome signatures identified the unique roles of the E loci in photoperiodic flowering and a set of genes controlled by these loci. To elucidate the regulatory gene networks underlying photoperiodic flowering regulation, we developed the network inference algorithmic package CausNet that integrates sparse linear regression and Granger causality heuristics, with Gaussian approximation of bootstrapping to provide reliability scores for predicted regulatory interactions. Using the transcriptome data, CausNet inferred regulatory interactions among soybean flowering genes. Published reports in the literature provided empirical verification for several of CausNet's inferred regulatory interactions. We further confirmed the inferred regulatory roles of the flowering suppressors GmCOL1a and GmCOL1b using GmCOL1 RNAi transgenic soybean plants. Combinations of the alleles of GmCOL1 and the major maturity locus E1 demonstrated positive interaction between these genes, leading to enhanced suppression of flowering transition. Our work provides novel insights and testable hypotheses in the complex molecular mechanisms of photoperiodic flowering control in soybean and lays a framework for de novo prediction of biological networks controlling important agronomic traits in crops.

Creation of Early Flowering Germplasm of Soybean by CRISPR/Cas9 Technology

Soybean is an important economic crop and a typical short-day crop, sensitive to photoperiod, and has narrow geographical adaptative region, which limit the creation of transgenic materials and reduce the breeding efficiency of new varieties. In addition, the genetic transformation efficiency of soybean is lower than that of many other crops, and the available receptor genotypes are limited. In this study, Agrobacterium-mediated transformation were used to introduce the CRISPR/Cas9 expression vector into soybean cultivar Jack and generated targeted mutants of E1 gene controlling soybean flowering. We obtained two novel types of mutations, 11 bp and 40 bp deletion at E1 coding region, respectively, and frameshift mutations produced premature translation termination codons and truncated E1 proteins, causing obvious early flowering under long day condition. In addition, no off-target effects were observed by predicting and analyzing the potential off-target sites of E1 targets. Significant decreased E1 gene expression of two novel mutants showed that the truncated E1 protein disinhibited GmFT2a/5a and increasing GmFT2a/5a gene expressions resulted obvious early flowering. Homozygous trans-clean mutants without T-DNA elements were also obtained and showed early flowering under long day condition. The photo-insensitive soybean transformation receptor we created laid a foundation for breeding excellent transgenic receptors suitable for high latitudes.

The Complex Interactions Between Flowering Behavior and Fiber Quality in Hemp

Hemp, Cannabis sativa L., is a sustainable multipurpose fiber crop with high nutrient and water use efficiency and with biomass of excellent quality for textile fibers and construction materials. The yield and quality of hemp biomass are largely determined by the genetic background of the hemp cultivar but are also strongly affected by environmental factors, such as temperature and photoperiod. Hemp is a facultative short-day plant, characterized by a strong adaptation to photoperiod and a great influence of environmental factors on important agronomic traits such as "flowering-time" and "sex determination." This sensitivity of hemp can cause a considerable degree of heterogeneity, leading to unforeseen yield reductions. Fiber quality for instance is influenced by the developmental stage of hemp at harvest. Also, male and female plants differ in stature and produce fibers with different properties and quality. Next to these causes, there is evidence for specific genotypic variation in fiber quality among hemp accessions. Before improved hemp cultivars can be developed, with specific flowering-times and fiber qualities, and adapted to different geographical regions, a better understanding of the molecular mechanisms controlling important phenological traits such as "flowering-time" and "sex determination" in relation to fiber quality in hemp is required. It is well known that genetic factors play a major role in the outcome of both phenological traits, but the major molecular factors involved in this mechanism are not characterized in hemp. Genome sequences and transcriptome data are available but their analysis mainly focused on the cannabinoid pathway for medical purposes. Herein, we review the current knowledge of phenotypic and genetic data available for "flowering-time," "sex determination," and "fiber quality" in short-day and dioecious crops, respectively, and compare them with the situation in hemp. A picture emerges for several controlling key genes, for which natural genetic variation may lead to desired flowering behavior, including examples of pleiotropic effects on yield quality and on carbon partitioning. Finally, we discuss the prospects for using this knowledge for the molecular breeding of this sustainable crop via a candidate gene approach.

Efficient QTL detection of flowering date in a soybean RIL population using the novel restricted two-stage multi-locus GWAS procedure

KEY MESSAGE

Eighty-six R1 QTLs accounting for 89.92% phenotypic variance in a soybean RIL population were identified using RTM-GWAS with SNPLDB marker which performed superior over CIM and MLM-GWAS with BIN/SNPLDB marker. A population (NJRIKY) composed of 427 recombinant inbred lines (RILs) derived from Kefeng-1 × NN1138-2 (MGII × MGV, MG maturity group) was applied for detecting flowering date (R1) quantitative trait locus (QTL) system in soybean. From a low-depth re-sequencing (~ 0.75 ×), 576,874 SNPs were detected and organized into 4737 BINs (recombination breakpoint determinations) and 3683 SNP linkage disequilibrium blocks (SNPLDBs), respectively. Using the association mapping procedures "Restricted Two-stage Multi-locus Genome-wide Association Study" (RTM-GWAS), "Mixed Linear Model Genome-wide Association Study" (MLM-GWAS) and the linkage mapping procedure "Composite Interval Mapping" (CIM), 67, 36 and 10 BIN-QTLs and 86, 14 and 23 SNPLDB-QTLs were detected with their phenotypic variance explained (PVE) 88.70-89.92% (within heritability 98.2%), 146.41-353.62% (overflowing) and 88.29-172.34% (overflowing), respectively. The RTM-GWAS with SNPLDBs which showed to be more efficient and reasonable than the others was used to identify the R1 QTL system in NJRIKY. The detected 86 SNPLDB-QTLs with their PVE from 0.02 to 30.66% in a total of 89.92% covered 51 out of 104 R1 QTLs in 18 crosses in SoyBase and 26 out of 139 QTLs in a nested association mapping population, while the rest 29 QTLs were novel ones. From the QTL system, 52 candidate genes were annotated, including the verified gene E1, E2, E9 and J, and grouped into 3 categories of biological processes, among which 24 genes were enriched into three protein-protein interaction networks, suggesting gene networks working together. Since NJRIKY involves only MGII and MGV, the QTL/gene system among MG000-MGX should be explored further.

Transcriptomic comparison reveals genetic variation potentially underlying seed developmental evolution of soybeans

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. However, the genetic variations underlying soybean domestication are not well known. Comparative transcriptomics revealed that a small portion of the orthologous genes might have been fast evolving. In contrast, three gene expression clusters were identified as divergent by their expression patterns, which occupied 37.44% of the total genes, hinting at an essential role for gene expression alteration in soybean domestication. Moreover, the most divergent stage in gene expression between wild and cultivated soybeans occurred during seed development around the cotyledon stage (15 d after fertilization, G15). A module in which the co-expressed genes were significantly down-regulated at G15 of wild soybeans was identified. The divergent clusters and modules included substantial differentially expressed genes (DEGs) between wild and cultivated soybeans related to cell division, storage compound accumulation, hormone response, and seed maturation processes. Chromosomal-linked DEGs, quantitative trait loci controlling seed weight and oil content, and selection sweeps revealed candidate DEGs at G15 in the fruit-related divergence of G. max and G. soja. Our work establishes a transcriptomic selection mechanism for altering gene expression during soybean domestication, thus shedding light on the molecular networks underlying soybean seed development and breeding strategy.

Natural variation in GmGBP1 promoter affects photoperiod control of flowering time and maturity in soybean

The present study screened for polymorphisms in coding and non-coding regions of the GmGBP1 gene in 278 soybean accessions with variable maturity and growth habit characteristics under natural field conditions in three different latitudes in China. The results showed that the promoter region was highly diversified compared with the coding sequence of GmGBP1. Five polymorphisms and four haplotypes were closely related to soybean flowering time and maturity through association and linkage disequilibrium analyses. Varieties with the polymorphisms SNP_-796G, SNP_-770G, SNP_-307T, InDel_-242normal, SNP_353A, or haplotypes Hap-3 and Hap-4 showed earlier flowering time and maturity in different environments. The shorter growth period might be largely due to higher GmGBP1 expression levels in soybean that were caused by the TCT-motif with SNP_-796G in the promoter. In contrast, the lower expression level of GmGBP1 in soybean caused by RNAi interference of GmGBP1 resulted in a longer growth period under different day lengths. Furthermore, the gene interference of GmGBP1 also caused a reduction in photoperiod response sensitivity (PRS) before flowering in soybean. RNA-seq analysis on GmGBP1 underexpression in soybean showed that 94 and 30 predicted genes were significantly upregulated and downregulated, respectively. Of these, the diurnal photoperiod-specific expression pattern of three significant flowering time genes GmFT2a, GmFT5a, and GmFULc also showed constantly lower mRNA levels in GmGBP1-i soybean than in wild type, especially under short day conditions. Together, the results showed that GmGBP1 functioned as a positive regulator upstream of GmFT2a and GmFT5a to activate the expression of GmFULc to promote flowering on short days.

Allelic Variation and Distribution of the Major Maturity Genes in Different Soybean Collections

Soybean time of flowering and maturity are genetically controlled by E genes. Different allelic combinations of these genes determine soybean adaptation to a specific latitude. The paper describes the first attempt to assess adaptation of soybean genotypes developed and realized at Institute of Field and Vegetable Crops, Novi Sad, Serbia [Novi Sad (NS) varieties and breeding lines] based on E gene variation, as well as to comparatively assess E gene variation in North-American (NA), Chinese, and European genotypes, as most of the studies published so far deal with North-American and Chinese cultivars and breeding material. Allelic variation and distribution of the major maturity genes (E1, E2, E3, and E4) has been determined in 445 genotypes from soybean collections of NA ancestral lines, Chinese germplasm, and European varieties, as well as NS varieties and breeding lines. The study showed that allelic combinations of E1-E4 genes significantly determined the adaptation of varieties to different geographical regions, although they have different impacts on maturity. In general, each collection had one major E genotype haplogroup, comprising over 50% of the lines. The exceptions were European varieties that had two predominant haplogroups and NA ancestral lines distributed almost evenly among several haplogroups. As e1-as/e2/E3/E4 was the most common genotype in NS population, present in the best-performing genotypes in terms of yield, this specific allele combination was proposed as the optimal combination for the environments of Central-Eastern Europe.

Genomic and molecular analysis of conserved and unique features of soybean PIF4

Phytochrome-interacting factor 4 (PIF4) participates in light signaling by interacting with photoreceptors, phytochromes, and cryptochromes. Although well characterized in Arabidopsis, PIF4's role in crop plants is unknown. Here we performed the first integrated genomics, transcriptomics, and molecular characterization of PIF4 in soybean (Glycine max) plants. Fifteen identified Glycine max PIFs (GmPIFs) grouped into PIF3, PIF4, and PIF8 subfamilies based on their phylogenetic relationships. The GmPIF4 subfamily formed two distinct clades (GmPIF4 I and GmPIF4 II) with different amino acid sequences in the conserved bHLH region. Quantitative transcriptional analysis of soybean plants exposed to different photoperiods and temperatures indicated that all PIF4 I clade GmPIF4s conserved PIF4-like expression. Three out of four GmPIF4 transcripts of the GmPIF4 I clade increased at 35 °C compared to 25 °C under short day conditions. RNA sequencing of soybeans undergoing floral transition showed differential regulation of GmPIF4b, and ectopic GmPIF4b expression in wild type Arabidopsis resulted in an early flowering phenotype. Complementation of GmPIF4b in Arabidopsis pif4-101 mutants partially rescued the mutant phenotype. PIF4 protein levels peaked before dawn, and a GmPIF4b protein variant was observed in soybean plants treated at high temperatures.

De novo assembly of a Chinese soybean genome

Soybean was domesticated in China and has become one of the most important oilseed crops. Due to bottlenecks in their introduction and dissemination, soybeans from different geographic areas exhibit extensive genetic diversity. Asia is the largest soybean market; therefore, a high-quality soybean reference genome from this area is critical for soybean research and breeding. Here, we report the de novo assembly and sequence analysis of a Chinese soybean genome for "Zhonghuang 13" by a combination of SMRT, Hi-C and optical mapping data. The assembled genome size is 1.025 Gb with a contig N50 of 3.46 Mb and a scaffold N50 of 51.87 Mb. Comparisons between this genome and the previously reported reference genome (cv. Williams 82) uncovered more than 250,000 structure variations. A total of 52,051 protein coding genes and 36,429 transposable elements were annotated for this genome, and a gene co-expression network including 39,967 genes was also established. This high quality Chinese soybean genome and its sequence analysis will provide valuable information for soybean improvement in the future.

Functional Genomics and Flowering Time in Medicago truncatula: An Overview

Flowering time is an important trait that influences adaptation and yield in many crop legumes. Both the inherent earliness of flowering and the degree to which it is responsive to environmental factors determine both the eco-geographic range across which crops can be successfully grown and the seasonal cycles most suitable for production. This chapter will provide a brief review of studies investigating the genetic control of flowering time in Medicago truncatula.

Quantitative Trait Locus Mapping of Flowering Time and Maturity in Soybean Using Next-Generation Sequencing-Based Analysis

Soybean (Glycine max L.) is a major legume crop that is mainly distributed in temperate regions. The adaptability of soybean to grow at relatively high latitudes is attributed to natural variations in major genes and quantitative trait loci (QTLs) that control flowering time and maturity. Identification of new QTLs and map-based cloning of candidate genes are the fundamental approaches in elucidating the mechanism underlying soybean flowering and adaptation. To identify novel QTLs/genes, we developed two F8:10 recombinant inbred lines (RILs) and evaluated the traits of time to flowering (R1), maturity (R8), and reproductive period (RP) in the field. To rapidly and efficiently identify QTLs that control these traits, next-generation sequencing (NGS)-based QTL analysis was performed. This study demonstrates that only one major QTL on chromosome 4 simultaneously controls R1, R8, and RP traits in the Dongnong 50 × Williams 82 (DW) RIL population. Furthermore, three QTLs were mapped to chromosomes 6, 11, and 16 in the Suinong 14 × Enrei (SE) RIL population. Two major pleiotropic QTLs on chromosomes 4 and 6 were shown to affect flowering time, maturity, and RP. A QTL influencing RP was identified on chromosome 11, and QTL on chromosome 16 was associated with time to flowering responses. All these QTLs contributed to soybean maturation. The QTLs identified in this study may be utilized in fine mapping and map-based cloning of candidate genes to elucidate the mechanisms underlying flowering and soybean adaptation to different latitudes and to breed novel soybean cultivars with optimal yield-related traits.

Developmental Pathways Are Blueprints for Designing Successful Crops

Genes controlling plant development have been studied in multiple plant systems. This has provided deep insights into conserved genetic pathways controlling core developmental processes including meristem identity, phase transitions, determinacy, stem elongation, and branching. These pathways control plant growth patterns and are fundamentally important to crop biology and agriculture. This review describes the conserved pathways that control plant development, using Arabidopsis as a model. Historical examples of how plant development has been altered through selection to improve crop performance are then presented. These examples, drawn from diverse crops, show how the genetic pathways controlling development have been modified to increase yield or tailor growth patterns to suit local growing environments or specialized crop management practices. Strategies to apply current progress in genomics and developmental biology to future crop improvement are then discussed within the broader context of emerging trends in plant breeding. The ways that knowledge of developmental processes and understanding of gene function can contribute to crop improvement, beyond what can be achieved by selection alone, are emphasized. These include using genome re-sequencing, mutagenesis, and gene editing to identify or generate novel variation in developmental genes. The expanding scope for comparative genomics, the possibility to engineer new developmental traits and new approaches to resolve gene-gene or gene-environment interactions are also discussed. Finally, opportunities to integrate fundamental research and crop breeding are highlighted.

Genomic adaptation of flowering-time genes during the expansion of rice cultivation area

The diversification of flowering time in response to natural environments is critical for the spread of crops to diverse geographic regions. In contrast with recent advances in understanding the molecular basis of photoperiodic flowering in rice (Oryza sativa), little is known about how flowering-time diversification is structured within rice subspecies. By analyzing genome sequencing data and a set of 429 chromosome segment substitution lines (CSSLs) originating from 10 diverse rice accessions with wide distributions, we revealed diverse effects of allelic variations for common flowering-time quantitative trait loci in the recipient's background. Although functional variations associated with a few loci corresponded to standing variations among subspecies, the identified functional nucleotide polymorphisms occurred recently after rice subgroup differentiation, indicating that the functional diversity of flowering-time gene sequences was not particularly associated with phylogenetic relationship between rice subspecies. Intensive analysis of the Hd1 genomic region identified the signature of an early introgression of the Hd1 with key mutation(s) in aus and temperate japonica accessions. Our data suggested that, after such key introgressions, new mutations were selected and accelerated the flowering-time diversity within subspecies during the expansion of rice cultivation area. This finding may imply that new genome-wide changes for flowering-time adaptation are one of the critical determinants for establishing genomic architecture of local rice subgroups. In-depth analyses of various rice genomes coupling with the genetically confirmed phenotypic changes in a large set of CSSLs enabled us to demonstrate how rice genome dynamics has coordinated with the adaptation of cultivated rice during the expansion of cultivation area.

Genotyping of Soybean Cultivars With Medium-Density Array Reveals the Population Structure and QTNs Underlying Maturity and Seed Traits

Soybean was domesticated about 5,000 to 6,000 years ago in China. Although genotyping technologies such as genotyping by sequencing (GBS) and high-density array are available, it is convenient and economical to genotype cultivars or populations using medium-density SNP array in genetic study as well as in molecular breeding. In this study, 235 cultivars, collected from China, Japan, USA, Canada and some other countries, were genotyped using SoySNP8k iSelect BeadChip with 7,189 single nucleotide polymorphisms (SNPs). In total, 4,471 polymorphic SNP markers were used to analyze population structure and perform genome-wide association study (GWAS). The most likely K value was 7, indicating this population can be divided into 7 subpopulations, which is well in accordance with the geographic origins of cultivars or accession studied. The LD decay rate was estimated at 184 kb, where r2 dropped to half of its maximum value (0.205). GWAS using FarmCPU detected a stable quantitative trait nucleotide (QTN) for hilum color and seed color, which is consistent with the known loci or genes. Although no universal QTNs for flowering time and maturity were identified across all environments, a total of 30 consistent QTNs were detected for flowering time (R1) or maturity (R7 and R8) on 16 chromosomes, most of them were corresponding to known E1 to E4 genes or QTL region reported in SoyBase (soybase.org). Of 16 consistent QTNs for protein and oil contents, 11 QTNs were detected having antagonistic effects on protein and oil content, while 4 QTNs soly for oil content, and one QTN soly for protein content. The information gained in this study demonstrated that the usefulness of the medium-density SNP array in genotyping for genetic study and molecular breeding.

Identification of novel loci associated with maturity and yield traits in early maturity soybean plant introduction lines

BACKGROUND

To continue to meet the increasing demands of soybean worldwide, it is crucial to identify key genes regulating flowering and maturity to expand the cultivated regions into short season areas. Although four soybean genes have been successfully utilized in early maturity breeding programs, new genes governing maturity are continuously being identified suggesting that there remains as yet undiscovered loci governing agronomic traits of interest. The objective of this study was to identify novel loci and genes involved in a diverse set of early soybean maturity using genome-wide association (GWA) analyses to identify loci governing days to maturity (DTM), flowering (DTF) and pod filling (DTPF), as well as yield and 100 seed weight in Canadian environments. To do so, soybean plant introduction lines varying significantly for maturity, but classified as early varieties, were used. Plants were phenotyped for the five agronomic traits for five site-years and GWA approaches used to identify candidate loci and genes affecting each trait.

RESULTS

Genotyping using genotyping-by-sequencing and microarray methods identified 67,594 single nucleotide polymorphisms, of which 31,283 had a linkage disequilibrium < 1 and minor allele frequency > 0.05 and were used for GWA analyses. A total of 9, 6, 4, 5 and 2 loci were detected for GWA analyses for DTM, DTF, DTPF, 100 seed weight and yield, respectively. Regions of interest, including a region surrounding the E1 gene for flowering and maturity, and several novel loci, were identified, with several loci having pleiotropic effects. Novel loci affecting maturity were identified on chromosomes five and 13 and reduced maturity by 7.2 and 3.3 days, respectively. Novel loci for maturity and flowering contained genes orthologous to known Arabidopsis flowering genes, while loci affecting yield and 100 seed weight contained genes known to cause dwarfism.

CONCLUSIONS

This study demonstrated substantial variation in soybean agronomic traits of interest, including maturity and flowering dates as well as yield, and the utility of GWA analyses in identifying novel genetic factors underlying important agronomic traits. The loci and candidate genes identified serve as promising targets for future studies examining the mechanisms underlying the related soybean traits.

Soybean MADS-box gene GmAGL1 promotes flowering via the photoperiod pathway

BACKGROUND

The MADS-box transcription factors are an ancient family of genes that regulate numerous physiological and biochemical processes in plants and facilitate the development of floral organs. However, the functions of most of these transcription factors in soybean remain unknown.

RESULTS

In this work, a MADS-box gene, GmAGL1, was overexpressed in soybean. Phenotypic analysis showed that GmAGL1 overexpression not only resulted in early maturation but also promoted flowering and affected petal development. Furthermore, the GmAGL1 was much more effective at promoting flowering under long-day conditions than under short-day conditions. Transcriptome sequencing analysis showed that before flowering, the photoperiod pathway photoreceptor CRY2 and several circadian rhythm genes, such as SPA1, were significantly down-regulated, while some other flowering-promoting circadian genes, such as GI and LHY, and downstream genes related to flower development, such as FT, LEAFY, SEP1, SEP3, FUL, and AP1, were up-regulated compared with the control. Other genes related to the flowering pathway were not noticeably affected.

CONCLUSIONS

The findings reported herein indicate that GmAGL1 may promote flowering mainly through the photoperiod pathway. Interestingly, while overexpression of GmAGL1 promoted plant maturity, no reduction in seed production or oil and protein contents was observed.

Large-Scale Investigation of Soybean Gene Functions by Overexpressing a Full-Length Soybean cDNA Library in Arabidopsis

Molecular breeding has become an important approach for crop improvement, and a prerequisite for molecular breeding is elucidation of the functions of genetic loci or genes. Soybean is one of the most important food and oil crops worldwide. However, due to the difficulty of genetic transformation in soybean, studies of its functional genomics lag far behind those of other crops such as rice, which severely impairs the progress of molecular improvement in soybean. Here, we describe an effective large-scale strategy to investigate the functions of soybean genes via overexpression of a full-length soybean cDNA library in Arabidopsis. The overexpression vector pJL12 was modified for use in the construction of a normalized full-length cDNA library. The constructed cDNA library showed good quality; repetitive clones represented approximately 4%, insertion fragments were approximately 2.2 kb, and the full-length rate was approximately 98%. This cDNA library was then overexpressed in Arabidopsis, and approximately 2000 transgenic lines were preliminarily obtained. Phenotypic analyses of the positive T₁ transgenic plants showed that more than 5% of the T₁ transgenic lines displayed abnormal developmental phenotypes, and approximately 1% of the transgenic lines exhibited potentially favorable traits. We randomly amplified 4 genes with obvious phenotypes (enlarged seeds, yellowish leaves, more branches, and dense siliques) and repeated the transgenic analyses in Arabidopsis. Subsequent phenotypic observation demonstrated that these phenotypes were indeed due to the overexpression of soybean genes. We believe our strategy represents an effective large-scale approach to investigate the functions of soybean genes and further reveal genes favorable for molecular improvement in soybean.

ZmCCT9 enhances maize adaptation to higher latitudes

Flowering time is a critical determinant of crop adaptation to local environments. As a result of natural and artificial selection, maize has evolved a reduced photoperiod sensitivity to adapt to regions over 90° of latitude in the Americas. Here we show that a distant Harbinger-like transposon acts as a cis-regulatory element to repress ZmCCT9 expression to promote flowering under the long days of higher latitudes. The transposon at ZmCCT9 and another functional transposon at a second flowering-time gene, ZmCCT10, arose sequentially following domestication and were targeted by selection as maize spread from the tropics to higher latitudes. Our results demonstrate that new functional variation created by transposon insertions helped maize to spread over a broad range of latitudes rapidly.

From its tropical origin in southwestern Mexico, maize spread over a wide latitudinal cline in the Americas. This feat defies the rule that crops are inhibited from spreading easily across latitudes. How the widespread latitudinal adaptation of maize was accomplished is largely unknown. Through positional cloning and association mapping, we resolved a flowering-time quantitative trait locus to a Harbinger-like transposable element positioned 57 kb upstream of a CCT transcription factor (ZmCCT9). The Harbinger-like element acts in cis to repress ZmCCT9 expression to promote flowering under long days. Knockout of ZmCCT9 by CRISPR/Cas9 causes early flowering under long days. ZmCCT9 is diurnally regulated and negatively regulates the expression of the florigen ZCN8, thereby resulting in late flowering under long days. Population genetics analyses revealed that the Harbinger-like transposon insertion at ZmCCT9 and the CACTA-like transposon insertion at another CCT paralog, ZmCCT10, arose sequentially following domestication and were targeted by selection for maize adaptation to higher latitudes. Our findings help explain how the dynamic maize genome with abundant transposon activity enabled maize to adapt over 90° of latitude during the pre-Columbian era.

Dissecting the Genetic Basis of Local Adaptation in Soybean

Soybean (Glycine max) is the most widely grown oilseed in the world and is an important source of protein for both humans and livestock. Soybean is widely adapted to both temperate and tropical regions, but a changing climate demands a better understanding of adaptation to specific environmental conditions. Here, we explore genetic variation in a collection of 3,012 georeferenced, locally adapted landraces from a broad geographical range to help elucidate the genetic basis of local adaptation. We used geographic origin, environmental data and dense genome-wide SNP data to perform an environmental association analysis and discover loci displaying steep gradients in allele frequency across geographical distance and between landrace and modern cultivars. Our combined application of methods in environmental association mapping and detection of selection targets provide a better understanding of how geography and selection may have shaped genetic variation among soybean landraces. Moreover, we identified several important candidate genes related to drought and heat stress, and revealed important genomic regions possibly involved in the geographic divergence of soybean.

Photoperiodism dynamics during the domestication and improvement of soybean

Soybean (Glycine max) is a facultative short-day plant with a sensitive photoperiod perception and reaction system, which allows it to adjust its physiological state and gene regulatory networks to seasonal and diurnal changes in environmental conditions. In the past few decades, soybean cultivation has spread from East Asia to areas throughout the world. Biologists and breeders must now confront the challenge of understanding the molecular mechanism of soybean photoperiodism and improving agronomic traits to enable this important crop to adapt to geographical and environmental changes. In this review, we summarize the genetic regulatory network underlying photoperiodic responses in soybean. Genomic and genetic studies have revealed that the circadian clock, in conjunction with the light perception pathways, regulates photoperiodic flowering. Here, we provide an annotated list of 844 candidate flowering genes in soybean, with their putative biological functions. Many photoperiod-related genes have been intensively selected during domestication and crop improvement. Finally, we describe recent progress in engineering photoperiod-responsive genes for improving agronomic traits to enhance geographic adaptation in soybean, as well as future prospects for research on soybean photoperiodic responses.

AtHKT1 drives adaptation of Arabidopsis thaliana to salinity by reducing floral sodium content

Arabidopsis thaliana high-affinity potassium transporter 1 (AtHKT1) limits the root-to-shoot sodium transportation and is believed to be essential for salt tolerance in A. thaliana. Nevertheless, natural accessions with 'weak allele' of AtHKT1, e.g. Tsu-1, are mainly distributed in saline areas and are more tolerant to salinity. These findings challenge the role of AtHKT1 in salt tolerance and call into question the involvement of AtHKT1 in salinity adaptation in A. thaliana. Here, we report that AtHKT1 indeed drives natural variation in the salt tolerance of A. thaliana and the coastal AtHKT1, so-called weak allele, is actually hyper-functional in reducing flowers sodium content upon salt stress. Our data showed that AtHKT1 positively contributes to saline adaptation in a linear manner. Forward and reverse genetics analysis established that the single AtHKT1 locus is responsible for the variation in the salinity adaptation between Col-0 and Tsu-1. Reciprocal grafting experiments revealed that shoot AtHKT1 determines the salt tolerance of Tsu-1, whereas root AtHKT1 primarily drives the salt tolerance of Col-0. Furthermore, evidence indicated that Tsu-1 AtHKT1 is highly expressed in stems and is more effective compared to Col-0 AtHKT1 at limiting sodium flow to the flowers. Such efficient retrieval of sodium to the reproductive organ endows Tsu-1 with stronger fertility compared to Col-0 upon salt stress, thus improving Tsu-1 adaptation to a coastal environment. To conclude, our data not only confirm the role of AtHKT1 in saline adaptation, but also sheds light on our understanding of the salt tolerance mechanisms in plants.

Cross-species complementation reveals conserved functions for EARLY FLOWERING 3 between monocots and dicots

Plant responses to the environment are shaped by external stimuli and internal signaling pathways. In both the model plant Arabidopsis thaliana (Arabidopsis) and crop species, circadian clock factors are critical for growth, flowering, and circadian rhythms. Outside of Arabidopsis, however, little is known about the molecular function of clock gene products. Therefore, we sought to compare the function of Brachypodium distachyon (Brachypodium) and Setaria viridis (Setaria) orthologs of EARLY FLOWERING 3, a key clock gene in Arabidopsis. To identify both cycling genes and putative ELF3 functional orthologs in Setaria, a circadian RNA-seq dataset and online query tool (Diel Explorer) were generated to explore expression profiles of Setaria genes under circadian conditions. The function of ELF3 orthologs from Arabidopsis, Brachypodium, and Setaria was tested for complementation of an elf3 mutation in Arabidopsis. We find that both monocot orthologs were capable of rescuing hypocotyl elongation, flowering time, and arrhythmic clock phenotypes. Using affinity purification and mass spectrometry, our data indicate that BdELF3 and SvELF3 could be integrated into similar complexes in vivo as AtELF3. Thus, we find that, despite 180 million years of separation, BdELF3 and SvELF3 can functionally complement loss of ELF3 at the molecular and physiological level.

Assessing and Exploiting Functional Diversity in Germplasm Pools to Enhance Abiotic Stress Adaptation and Yield in Cereals and Food Legumes

There is a need to accelerate crop improvement by introducing alleles conferring host plant resistance, abiotic stress adaptation, and high yield potential. Elite cultivars, landraces and wild relatives harbor useful genetic variation that needs to be more easily utilized in plant breeding. We review genome-wide approaches for assessing and identifying alleles associated with desirable agronomic traits in diverse germplasm pools of cereals and legumes. Major quantitative trait loci and single nucleotide polymorphisms (SNPs) associated with desirable agronomic traits have been deployed to enhance crop productivity and resilience. These include alleles associated with variation conferring enhanced photoperiod and flowering traits. Genetic variants in the florigen pathway can provide both environmental flexibility and improved yields. SNPs associated with length of growing season and tolerance to abiotic stresses (precipitation, high temperature) are valuable resources for accelerating breeding for drought-prone environments. Both genomic selection and genome editing can also harness allelic diversity and increase productivity by improving multiple traits, including phenology, plant architecture, yield potential and adaptation to abiotic stresses. Discovering rare alleles and useful haplotypes also provides opportunities to enhance abiotic stress adaptation, while epigenetic variation has potential to enhance abiotic stress adaptation and productivity in crops. By reviewing current knowledge on specific traits and their genetic basis, we highlight recent developments in the understanding of crop functional diversity and identify potential candidate genes for future use. The storage and integration of genetic, genomic and phenotypic information will play an important role in ensuring broad and rapid application of novel genetic discoveries by the plant breeding community. Exploiting alleles for yield-related traits would allow improvement of selection efficiency and overall genetic gain of multigenic traits. An integrated approach involving multiple stakeholders specializing in management and utilization of genetic resources, crop breeding, molecular biology and genomics, agronomy, stress tolerance, and reproductive/seed biology will help to address the global challenge of ensuring food security in the face of growing resource demands and climate change induced stresses.

Characterization of Genetic Basis on Synergistic Interactions between Root Architecture and Biological Nitrogen Fixation in Soybean

Soybean [Glycine max (L.) Merr] is an important legume crop and its yield largely depends on root architecture (RA) and biological nitrogen fixation (BNF). However, the relationship between RA and BNF, and its genetics behind remain unclear. Here, two soybean genotypes contrasting in RA and their 175 F_9:11 recombinant inbred lines (RILs) were evaluated in field. The shallow-root parent, JD12, had better nodulation and higher yield than the deep-root parent, NF58. Strong correlations between shoot dry weight (SDW) and RA or BNF traits existed in the RILs, and the shallow-root group had more and heavier nodules, as well as higher SDW. After inoculating with rhizobia, roots became shallower and bigger, showing strong synergistic interactions between RA and BNF. In total, 70 QTLs were identified for the 21 tested traits. Among them, qBNF-RA-C2, qBNF-RA-O, and qBNF-RA-B1, were newly identified QTLs for BNF and/or RA traits in soybean, which co-located with the QTLs for SDW detected presently, and with the QTLs for yield identified previously. The results together suggest that there are synergistic interactions between RA and BNF, and the QTLs identified here could be used for breeding new soybean varieties with higher yields through optimization of RA traits and BNF capacity.