Diurnal Expression Pattern,Allelic Variation,and Association Analysis Reveal Functional Features of the E1 Gene in Control of Photoperiodic Flowering in Soybean

Although four maturity genes, E1 to E4, in soybean have been successfully cloned, their functional mechanisms and the regulatory network of photoperiodic flowering remain to be elucidated. In this study, we investigated how the diurnal expression pattern of the E1 gene is related to photoperiodic length; and to what extent allelic variation in the B3-like domain of the E1 gene is associated with flowering time phenotype. The bimodal expression of the E1 gene peaked first at around 2 hours after dawn in long-day condition. The basal expression level of E1 was enhanced by the long light phase, and decreased by duration of dark. We identified a 5bp (3 SNP and 2-bp deletion) mutation, referred to an e1-b3a, which occurs in the middle of B3 domain of the E1 gene in the early flowering cultivar Yanhuang 3. Subcellular localization analysis showed that the putative truncated e1-b3a protein was predominately distributed in nuclei, indicating the distribution pattern of e1-b3a was similar to that of E1, but not to that of e1-as. Furthermore, genetic analysis demonstrated allelic variations at the E1 locus significantly underlay flowering time in three F₂ populations. Taken together, we can conclude the legume specific E1 gene confers some special features in photoperiodic control of flowering in soybean. Further characterization of the E1 gene will extend our understanding of the soybean flowering pathway in soybean.     

Nucleotide diversity of the upstream region of the putative MADS-box gene controlling soybean maturity

MADS-box genes are involved in plant reproductive development. However, the role of gene nucleotide diversity in soybean flowering and maturity remains unknown. Therefore, in this study, the distribution of DNA polymorphisms in the putative MADS-box gene located near the quantitative trait loci (QTL) for flowering time and maturity was targeted for association analysis using Glycine max (cultivated soybean) and Glycine soja (wild soybean). Sixteen single nucleotide polymorphisms identified in the upstream region of the putative MADS-box gene around QTL Pod mat 13-7 and Fflr 4-2 on chromosome 7 were found to be highly associated with maturity in soybean. The genetic diversity between cultivated soybeans and the wild relative was comparable, although the early maturity group (EMG) was less diverse than the late maturity group (LMG) of the cultivated soybean. Population size changes of the MADS-box gene in this soybean germplasm appeared to result from non-random selection. A selective pressure seemed to act on this gene in the EMG, while the LMG and G. soja were in genetic equilibrium. Neutrality tests and the constructed neighbor-joining tree indicate that the EMG of G. max has experienced strong artificial selection for its domestication and genetic improvement.     

Mapping and identification of a potential candidate gene for a novel maturity locus, <Emphasis Type="Italic">E10</Emphasis>, in soybean

Key message

E10 is a new maturity locus in soybean and FT4 is the predicted/potential functional gene underlying the locus.

Abstract

Flowering and maturity time traits play crucial roles in economic soybean production. Early maturity is critical for north and west expansion of soybean in Canada. To date, 11 genes/loci have been identified which control time to flowering and maturity; however, the molecular bases of almost half of them are not yet clear. We have identified a new maturity locus called “E10” located at the end of chromosome Gm08. The gene symbol E10e10 has been approved by the Soybean Genetics Committee. The e10e10 genotype results in 5–10 days earlier maturity than E10E10. A set of presumed E10E10 and e10e10 genotypes was used to identify contrasting SSR and SNP haplotypes. These haplotypes, and their association with maturity, were maintained through five backcross generations. A functional genomics approach using a predicted protein–protein interaction (PPI) approach (Protein–protein Interaction Prediction Engine, PIPE) was used to investigate approximately 75 genes located in the genomic region that SSR and SNP analyses identified as the location of the E10 locus. The PPI analysis identified FT4 as the most likely candidate gene underlying the E10 locus. Sequence analysis of the two FT4 alleles identified three SNPs, in the 5′UTR, 3′UTR and fourth exon in the coding region, which result in differential mRNA structures. Allele-specific markers were developed for this locus and are available for soybean breeders to efficiently develop earlier maturing cultivars using molecular marker assisted breeding.

     

Soybean adaption to high-latitude regions is associated with natural variations of GmFT2b,an ortholog of FLOWERING LOCUS T

Day length has an important influence on flowering and growth habit in many plant species. In crops such as soybean, photoperiod sensitivity determines the geographical range over which a given cultivar can grow and flower. The soybean genome contains ~10 genes homologous to FT, a central regulator of flowering from Arabidopsis thaliana. However, the precise roles of these soybean FTs are not clearly. Here we show that one such gene, GmFT2b, promotes flowering under long-days (LDs). Overexpression of GmFT2b upregulates expression of flowering-related genes which are important in regulating flowering time. We propose a ‘weight’ model for soybean flowering under short-day (SD) and LD conditions. Furthermore, we examine GmFT2b sequences in 195 soybean cultivars, as well as flowering phenotypes, geographical distributions and maturity groups. We found that Hap3, a major GmFT2b haplotype, is associated with significantly earlier flowering at higher latitudes. We anticipate our assay to provide important resources for the genetic improvement of soybean, including new germplasm for soybean breeding, and also increase our understanding of functional diversity in the soybean FT gene family.     

Allelic Combinations of Soybean Maturity Loci E1, E2, E3 and E4 Result in Diversity of Maturity and Adaptation to Different Latitudes

Soybean cultivars are extremely diverse in time to flowering and maturation as a result of various photoperiod sensitivities. The underlying molecular genetic mechanism is not fully clear, however, four maturity loci E1, E2, E3 and E4 have been molecularly identified. In this report, cultivars were selected with various photoperiod sensitivities from different ecological zones, which covered almost all maturity groups (MG) from MG 000 to MG VIII and MG X adapted from latitude N 18° to N 53°. They were planted in the field under natural daylength condition (ND) in Beijing, China or in pots under different photoperiod treatments. Maturity-related traits were then investigated. The four E maturity loci were genotyped at the molecular level. Our results suggested that these four E genes have different impacts on maturity and their allelic variations and combinations determine the diversification of soybean maturity and adaptation to different latitudes. The genetic mechanisms underlying photoperiod sensitivity and adaptation in wild soybean seemed unique from those in cultivated soybean. The allelic combinations and functional molecular markers for the four E loci will significantly assist molecular breeding towards high productivity.     

Identification of mega‐environments in Europe and effect of allelic variation at maturity E loci on adaptation of European soybean

Soybean cultivation holds great potential for a sustainable agriculture in Europe, but adaptation remains a central issue. In this large mega‐environment (MEV) study, 75 European cultivars from five early maturity groups (MGs 000–II) were evaluated for maturity‐related traits at 22 locations in 10 countries across Europe. Clustering of the locations based on phenotypic similarity revealed six MEVs in latitudinal direction and suggested several more. Analysis of maturity identified several groups of cultivars with phenotypic similarity that are optimally adapted to the different growing regions in Europe. We identified several haplotypes for the allelic variants at the E1, E2, E3 and E4 genes, with each E haplotype comprising cultivars from different MGs. Cultivars with the same E haplotype can exhibit different flowering and maturity characteristics, suggesting that the genetic control of these traits is more complex and that adaptation involves additional genetic pathways, for example temperature requirement. Taken together, our study allowed the first unified assessment of soybean‐growing regions in Europe and illustrates the strong effect of photoperiod on soybean adaptation and MEV classification, as well as the effects of the E maturity loci for soybean adaptation in Europe.     

Genome-Wide Association Studies Identifies Seven Major Regions Responsible for Iron Deficiency Chlorosis in Soybean (Glycine max)

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.     

基于大豆胞囊线虫病抗性候选基因的SNP位点遗传变异分析

以我国363份栽培和野生大豆资源为材料, 对大豆胞囊线虫抗性候选基因(rhg1和Rhg4)的SNP位点(8个)进行遗传变异分析, 以期阐明野生和栽培大豆间遗传多样性及连锁不平衡水平差异。结果表明, 与野生大豆相比, 代表我国栽培大豆总体资源多样性的微核心种质及其补充材料的连锁不平衡水平较高(R²值为0.216)。在栽培大豆群体内, 基因内和基因间分别有100％和16.6％的SNP位点对连锁不平衡显著, 形成两个基因特异的连锁不平衡区间(Block)。在所有供试材料中共检测到单倍型46个, 野生大豆的单倍型数目(27)少于栽培大豆(31), 但单倍型多样性(0.916)稍高于栽培大豆(0.816)。单倍型大多数(67.4％)为群体所特有(31个), 其中15个为野生大豆特有单倍型。野生大豆的两个主要优势单倍型(Hap_10和Hap_11)在栽培大豆中的发生频率也明显下降, 推测野生大豆向栽培大豆进化过程中, 一方面形成了新的单倍型, 另一方面因为瓶颈效应部分单倍型的频率降低甚至消失。     

Molecular markers for the E2 and E3 genes controlling flowering and maturity in soybean

Natural variation in flowering time may play a role in the adaptation of plants to various environments, and understanding the genetic basis of flowering and maturity would facilitate the development of early maturing cultivars. Molecular markers for the E2 and E3 loci, which control the time of flowering and maturity in soybean (Glycine max), were developed in this study. Single nucleotide-amplified polymorphism (SNAP) markers were developed from the nonsense mutation in E2 (GmGIa), which is a circadian clock-controlled gene. The E2- and e2-specific SNAP markers were validated using six E2 isolines. The soybean E3 gene is a photoreceptor phytochrome A (GmPhyA3) gene, and a co-dominant marker was designed based on sequence deletions within the E3 allele. A multiplex PCR assay using three primers for the E3 gene allowed allelic discrimination based on the sizes of PCR products. Furthermore, this E3 marker successfully detected two alleles in a single reaction when two types of DNA were pooled. These markers determined the genotypes of our mapping population previously reported to detect flowering quantitative trait loci close to the E2 and E3 loci, confirming that the mutations are responsible for the early flowering phenotype. The use of SNAP markers for E2 and a co-dominant marker for E3 is a simple, fast, and reproducible method, requiring only PCR and agarose gel electrophoresis. The molecular resources developed in this study could accelerate marker-assisted selection and cultivar development for short-season areas in a soybean breeding program.     

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.     

Natural variation in the genes responsible for maturity loci E1, E2, E3 and E4 in soybean

Background and Aims

The timing of flowering has a direct impact on successful seed production in plants. Flowering of soybean (Glycine max) is controlled by several E loci, and previous studies identified the genes responsible for the flowering loci E1, E2, E3 and E4. However, natural variation in these genes has not been fully elucidated. The aims of this study were the identification of new alleles, establishment of allele diagnoses, examination of allelic combinations for adaptability, and analysis of the integrated effect of these loci on flowering.

Methods

The sequences of these genes and their flanking regions were determined for 39 accessions by primer walking. Systematic discrimination among alleles was performed using DNA markers. Genotypes at the E1–E4 loci were determined for 63 accessions covering several ecological types using DNA markers and sequencing, and flowering times of these accessions at three sowing times were recorded.

Key Results

A new allele with an insertion of a long interspersed nuclear element (LINE) at the promoter of the E1 locus (e1-re) was identified. Insertion and deletion of 36 bases in the eighth intron (E2-in and E2-dl) were observed at the E2 locus. Systematic discrimination among the alleles at the E1–E3 loci was achieved using PCR-based markers. Allelic combinations at the E1–E4 loci were found to be associated with ecological types, and about 62–66 % of variation of flowering time could be attributed to these loci.

Conclusions

The study advances understanding of the combined roles of the E1–E4 loci in flowering and geographic adaptation, and suggests the existence of unidentified genes for flowering in soybean,     

Phylogeny and genetic diversity of D-genome species of Aegilops and Triticum (Triticeae, Poaceae) from Iran based on microsatellites, ITS, and trnL-F

Cereal species of the grass tribe Triticeae are economically important and provide staple food for large parts of the human population. The Fertile Crescent of Southwest Asia harbors high genetic and morphological diversity of these species. In this study, we analyzed genetic diversity and phylogenetic relationships among D genome-bearing species of the wheat relatives of the genus Aegilops from Iran and adjacent areas using allelic diversity at 25 nuclear microsatellite loci, nuclear rDNA ITS, and chloroplast trnL-F sequences. Our analyses revealed high microsatellite diversity in Aegilops tauschii and the D genomes of Triticum aestivum and Ae. ventricosa, low genetic diversity in Ae. cylindrica, two different Ae. tauschii gene pools, and a close relationship among Ae. crassa, Ae. juvenalis, and Ae. vavilovii. In the latter species group, cloned sequences revealed high diversity at the ITS region, while in most other polyploids, homogenization of the ITS region towards one parental type seems to have taken place. The chloroplast genealogy of the trnL-F haplotypes showed close relationships within the D genome Aegilops species and T. aestivum, the presence of shared haplotypes in up to three species, and up to three different haplotypes within single species, and indicates chloroplast capture from an unidentified species in Ae. markgrafii. The ITS phylogeny revealed Triticum as monophyletic and Aegilops as monophyletic when Amblyopyrum muticum is included.     

Resequencing 250 Soybean Accessions: New Insights into Genes Associated with Agronomic Traits and Genetic Networks

The limited knowledge of genomic diversity and functional genes associated with the traits of soybean varieties has resulted in slow progress in breeding. In this study, we sequenced the genomes of 250 soybean landraces and cultivars from China, America, and Europe, and investigated their population structure, genetic diversity and architecture, and the selective sweep regions of these accessions. Five novel agronomically important genes were identified, and the effects of functional mutations in respective genes were examined. The candidate genes GSTT1, GL3, and GSTL3 associated with the isoflavone content, CKX3 associated with yield traits, and CYP85A2 associated with both architecture and yield traits were found. The phenotype–gene network analysis revealed that hub nodes play a crucial role in complex phenotypic associations. This study describes novel agronomic trait-associated genes and a complex genetic network, providing a valuable resource for future soybean molecular breeding.