Legume Research

  • Chief EditorJ. S. Sandhu

  • Print ISSN 0250-5371

  • Online ISSN 0976-0571

  • NAAS Rating 6.80

  • SJR 0.391

  • Impact Factor 0.8 (2024)

Frequency :
Monthly (January, February, March, April, May, June, July, August, September, October, November and December)
Indexing Services :
BIOSIS Preview, ISI Citation Index, Biological Abstracts, Elsevier (Scopus and Embase), AGRICOLA, Google Scholar, CrossRef, CAB Abstracting Journals, Chemical Abstracts, Indian Science Abstracts, EBSCO Indexing Services, Index Copernicus
Legume Research, volume 44 issue 5 (may 2021) : 515-521

Field-based Screening and Haplotyping of J Locus for Long Juvenile Trait in Tropical Soybean Genotypes

Enoch Sapey1,2, Yudong Fang1, Dalia Mohamedkheir Khojely1, Wenwen Song1, Bingjun Jiang1, Shan Yuan1, Shi Sun1, Tianfu Han1,*
1MARA Key Laboratory of Soybean Biology (Beijing), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China.
2CSIR-Oil Palm Research Institute, P.O Box 74. Kade, Ghana.
  • Submitted08-12-2020|

  • Accepted02-03-2021|

  • First Online 26-03-2021|

  • doi 10.18805/LR-603

Cite article:- Sapey Enoch, Fang Yudong, Khojely Mohamedkheir Dalia, Song Wenwen, Jiang Bingjun, Yuan Shan, Sun Shi, Han Tianfu (2021). Field-based Screening and Haplotyping of J Locus for Long Juvenile Trait in Tropical Soybean Genotypes . Legume Research. 44(5): 515-521. doi: 10.18805/LR-603.
Background: The long juvenile trait influences flowering time of soybean under tropical conditions. The trait ensures sufficient vegetative growth prior to flowering. The present study aimed at identifying tropical soybean genotypes with the long juvenile trait and harboring the loss-of-function alleles for the J gene and verifying the effect of loss-of-function alleles on the long juvenile trait. 

Methods: A total of 159 soybean genotypes were evaluated on the field for two years in tropical city of Sanya, Hainan Islands, China for days to beginning bloom (VE-R1), days to physiological maturity (VE-R7), nodes/plant, pods/plant and plant height at maturity. The full length of the J gene was cloned in the 159 soybean genotypes. The sequence data was subjected to haplotype analysis to determine the genotypes with the functional and loss of function alleles. 

Result: Significant differences (p<0.05) were observed among genotypes for all phenotypic traits. 53 genotypes were identified with delayed flowering based on the phenotypic data. Sequence comparison of the 159 genotypes identified 10 polymorphisms comprising 7 SNPs and three deletions in the coding sequences, the three deletions resulted in the loss of function of the J gene. Six genotypes with delayed flowering had the loss-of-function alleles of the J gene.
Soybean is one of the key leguminous crops in the world that provides protein and oil for human use and also serves as animal feed. This rising demand for soybean is creating economic opportunities to expand soybean production in tropical regions. Soybean was domesticated ~5000 years ago in the upper reach of Yellow River in China at approximately 35°N latitude characterized by long days (Hymowitz, 1970). However, soybean as a short-day plant is sensitive to photoperiod and its flowering is induced when the day length is shorter than a critical threshold (Yang et al., 2019). These characteristic limits the expansion of soybean cultivation in latitudes less than 20° where day length deviates around 12 to 13 hours.
 
Production of soybean has conventionally centered around mid- and high-latitude temperate regions. Low-latitude areas with tropical and sub-tropical climates were previously considered unsuitable for soybean production (Destro et al., 2001). Hartwig and Kiihl (1979) identified a recessive trait in soybean germplasm that delayed flowering under short-day conditions and has been referred to the long juvenile trait. The long juvenile trait influences flowering time in low latitude areas, which characterize tropical climates and allows for a broader range of sowing dates, delay flowering and ensures sufficient vegetative growth prior to reproductive growth (Carpentieri-Pípoloet_al2000).
 
The use of long juvenile trait in soybean breeding programs in South America has led to a tremendous increase in soybean production, even though almost half of its total soybean production area is located below 24°S outside the traditional soybean adaptation zone (Neumaier and James 1993). In Africa, the long juvenile trait, together with the improved agronomic performances is encouraging an increase in soybean production.
 
To date, only two loci, E6 and J, have been reported to control the long juvenile trait (Ray et al., 1995; Bonato and Vello 1999). Studies by Yue et al., (2017) and Lu et al., (2017) identified the J gene as the Arabidopsis flowering ortholog ELF3. GmELF3 has 4 exons/3 introns and is a highly conserved protein that controls flowering time in multiple species (Lu et al., 2017). They suggested that J is the dominant functional allele of ELF3 that promotes flowering, while j is the recessive, loss-of-function allele that causes delayed flowering and confers the long juvenile trait that is associated with the adaptation of soybean to tropical regions.

The objectives of this study were (i) to identify tropical soybean genotypes carrying the long juvenile trait and the loss-of-function alleles that are responsible for long juvenile trait in soybean, (ii) to verify the relationship between the J mutations and the long juvenile trait. The newly identified genotypes and the data provided will useful for tropical soybean breeding programs.
A total of 159 tropical soybean genotypes from different tropical regions were used in the study (Table 1). The tropical soybean accessions were retrieved from the Gene bank of the Institute of Crop Sciences, Chinese Academy of Agricultural Sciences (CAAS) and Agricultural Research Corporation (ARC), Wad Medani, Sudan.
 
Field phenotyping
 
The 159 soybean genotypes were planted in 1.5 m single row plot; the spacing between rows was 30 cm with 5 cm spacing within plants with two replications. Seeds were sown at CAAS Experimental Station-Sanya (18.27.17° N, 109.11° E), Hainan, China, for two years on the 19th and 9th of November, 2018 and 2019, respectively. The daylength in Sanya during the field evaluation ranged from 11.4 to 11.9 hours and the temperature ranged from 19-26°C. The management of the field experiments was done according to recommended agronomic practices for soybean. Five plants were randomly tagged per row for data collection. Flowering time was recorded at the VE-R1 stage [days from emergence (VE) to the beginning bloom (R1)] and physiological maturity (R7) recorded at VE-R7 stage (days from emergence to when at least a pod had mature color), by the methods described by Fehr et al., (1971). Plant height from the ground to the terminal raceme apex, number of stem nodes per plant and number of pods per plant were all recorded at physiological maturity.
 
Statistical analysis of the phenotypic data  
 
The data collected were subjected to analysis of variance (ANOVA) and descriptive statistics computed using SPSS software. Boxplots were constructed using the ggplot2 package in R software. 
 
Genotyping for J gene
 
DNA was extracted from the 159 soybean accessions using the NuClean Plant Genomic DNA Kit (CoWin Biosciences) according to the manufacturer’s recommendation. The quality of the extracted DNA was evaluated by 1% agarose gel electrophoresis. The DNA concentration of each sample was determined using Nanodrop 2000 (Thermo Scientific). Six primers were designed covering the four exons of the J gene. Two were used for Polymerase Chain Reaction (PCR) and the remaining four for sequencing. Table 2 shows the list of the six primers used for genotyping.
 
PCR was done to amplify genomic regions using KOD-Plus-Neo Kit (Toyobo Life Science). The total volume of the PCR mix of 50 µL consisted of 5 µL of PCR Buffer, 5 µL of 2 mM dNTPs, 3 µL of 25 mM MgSO4, 1.5 µL of 10 mM forward and reverse primer, 5 µL of DNA template at a concentration about 100 ng, 1 µL KOD Plus Neo and 27 µL of double-distilled water. The PCR conditions were as follows: denaturation for 2 minutes at 94°C followed by 35 cycles of 10 sec at 98°C and a primer extension reaction of 68°C for 5 minutes. The PCR products were separated on 1% agarose gel, stained with ethidium and visualized with Gel Doc XR Molecular Imager System (Bio-Rad). The PCR products were directly Sanger sequenced by TSINGKE Biotech, China. DNA Sequence assembly and alignment were done using the Mega software.
 
Definition of J/j alleles
 
The definition of the J/j alleles was done by comparing the sequence data with sequences deposited in the National Centre for Biotechnology Information (NCBI) database by Yue et al., (2017) and Lu et al., (2017).
Phenotypic variations of days to flowering, maturity and other associated agronomic traits
 
Significant differences (p <0.05) were observed among the genotypes for all the phenotypic traits measured (Table 3). There was also a significant difference between years and the genotype by year interactions for these traits, except for VE-R7. This shows a high genetic variability of the tropical genotypes used in the study. Days to from emergence to beginning bloom (VE-R1) ranged from 21.67 to 41.45 days among the 159 genotypes. A box plot of VE-RI (Fig 1) showed two outliers (TGX 1987-62F and IC192) with flowering days of 35.4 and 41.8, respectively. Days from emergence to physiological maturity (VE-R7) ranged from 63.89 to 85.08 days. The plant height of the genotypes ranged from 15.3 to 68.4 cm, a box plot of plant height of genotypes showed two outliers (IC192 and CIGRAS-06). The plant heights of the two genotypes were 68.5 and 64.3 cm, respectively and the number of nodes and pods per plant ranged from 5.43 to 13.92 and 4 to 27.59 respectively.

Fig 1: Box plot of the five phenotypic traits observed on the 159 genotypes.



Table 3: Combined ANOVA for all the five phenotypic traits.



Significant positive correlations were observed among all phenotypic traits studied (Table 4). A significant (p<0.05) and positive correlation was observed between VE-R1 and VE-R7. The correlations of plant height at maturity with the numbers of nodes per plant and pods per plant were highly positive. Plant height is an important trait that affects soybean adaptation and yield. A positive significant correlation was observed among the phenotypic traits indicating that delayed flowering resulted in an increase in the other important agronomic traits such as plant height, which sometimes correlates with high yields in tropical regions (Cober and Morrison, 2010).

Table 4: Correlation coefficients among the five phenotypic traits.

 

Selection of the genotypes that delayed flowering based on the phenotypic data
 
Sanya-Hainan Islands, the experimental site used for phenotyping of the long juvenile trait, is the southernmost region of China with a typical tropical climate. The day length and temperatures provided an ideal environment for screening for the long juvenile trait in China. Two long juvenile varieties Huaxia-3 and Huaxia-10 reported by Yue et al., (2017) and Lu et al., (2017) were included in the study to serve as checks. The two check varieties flowered about 32-33 days after emergence. Since the determination of juvenility of a genotype is determined by the photoperiod of the given environment, genotypes that have same or similar flowering days as the check varieties or flowered later than the check varieties were considered to be long juvenile. Based on the phenotypic data we have identified 53 genotypes that flowered 30 days and above (Table 5).

Table 5: List of selected delayed flowering genotypes.

 

SNP haplotyping to distinguish different alleles of the J gene
 
Based on the phenotypic differences observed, we have genotyped the 159 soybean genotypes and carried out haplotype analysis to determine the genotypes carrying the functional or non-functional alleles of the gene. Sequence comparison identified 10 polymorphisms comprising of seven SNPs and three deletions in the coding sequence, three of which resulted in the loss-of-function of the J gene. Fig 2 shows the position of the base substitution in the coding sequence and the number of haplotypes identified. Seven haplotypes were defined of which three haplotypes (j1, j2 and j4) resulted in the loss-of-function of the J gene. The other functional haplotypes include Hap1, Hap2, Hap3 and Hap 4. Only one soybean germplasm from Brazil MG/BR-48 had the loss-of-function allele j2. Huaxia 3 and Huaxia 10, the two Chinese tropical check varieties, had the j4 allele. Two genotypes TGx 1835-10E and TGx 1987-62F from IITA-Nigeria had the j2 allele. I.C.192 from the US had the j1 allele (Table 5, Fig 3). The identification of these alleles confirms studies by Lu et al., (2017) that the multiple loss-of-function J alleles existing in soybean germplasm worldwide may be responsible for delayed flowering and maturity.

Fig 2: Haplotypes and polymorphisms in the coding sequence of the J gene (Glyma.04G050200) in the 159 soybean genotypes. The haplotype is shown as a linear combination of alleles.



Fig 3: Variation in the five phenotypic traits of the various haplotypes of the J gene.


 
Association of the haplotypes with the phenotypic traits
 
Fig 3 shows the variability of the phenotypic traits among the haplotypes defined. Haplotypes Hap1 and Hap3 are the most common among the genotypes with high variability in all the six phenotypic traits. Hap 4 is the least common in only two genotypes. Some genotypes with functional J haplotypes had similar flowering days and, in some cases, even higher (Table 3) than those harboring the loss-of-function alleles indicating that the J gene might not be the only gene responsible for the long juvenile trait. The same observations were recorded on the other associated phenotypic traits such as the plant height, number of nodes per plant and number of pods per plant. This confirms the polygenic nature of the long juvenile trait by previous studies (Ray et al., 1995, Carpentieri-Pipolo et al., 2000, Cober and Morrison 2010, Lu et al., 2017, Yue et al., 2017 and Fang et al., 2020).
 
The present study mainly focused on identifying new long juvenile genotypes and verifying whether the J mutations are associated with delayed flowering in soybean. The long juvenile genetic mechanism in most tropical genotypes especially Sub-Sahara Africa is not known (Miranda et al., 2020). Huaxia 3, the check variety, was introduced to some African countries and resulted in an increased productivity (Yue et al., 2017). Introgression of the loss-of-function alleles Hap j2 and Hap j4 to temperate cultivars has also resulted in the development of new long juvenile varieties in Brazil and southern China (Lu et al., 2017), indicating that the loss-of-function J alleles will be a good strategy to enhance the adaptation and yield increases in the tropics. Therefore, there exists a possibility that incorporation of the identified long juvenile genotypes with or without the loss of function alleles for the J gene in tropical soybean breeding programs will lead to yield increases. However, proper yield trials have to be conducted first, as past studies reported background effects could influence delayed flowering in soybean (Ray et al., 1995). The lack of Selection of the genotypes that delayed flowering based on the phenotypic data Sanya-Hainan Islands, the experimental site used for phenotyping of the long juvenile trait, is the southernmost region of China with a typical tropical climate. The day length and temperatures provided an ideal environment for screening for the long juvenile trait in China. Two long juvenile varieties Huaxia-3 and Huaxia-10 reported by Yue et al., (2017) and Lu et al., (2017) were included in the study to serve as checks. The two check varieties flowered about 32-33 days after emergence. Since the determination of juvenility of a genotype is determined by the photoperiod of the given environment, genotypes that have same or similar flowering days as the check varieties or flowered later than the check varieties were considered to be long juvenile. Based on the phenotypic data we have identified 53 genotypes that flowered 30 days and above (Table 5).
 
SNP haplotyping to distinguish different alleles of the J gene
 
Based on the phenotypic differences observed, we have genotyped the 159 soybean genotypes and carried out haplotype analysis to determine the genotypes carrying the functional or non-functional alleles of the gene. Sequence comparison identified 10 polymorphisms comprising of seven SNPs and three deletions in the coding sequence, three of which resulted in the loss-of-function of the J gene. Fig 2 shows the position of the base substitution in the coding sequence and the number of haplotypes identified. Seven haplotypes were defined of which three haplotypes (j1, j2 and j4) resulted in the loss-of-function of the J gene. The other functional haplotypes include Hap1, Hap2, Hap3 and Hap 4. Only one soybean germplasm from Brazil MG/BR-48 had the loss-of-function allele j2. Huaxia 3 and Huaxia 10, the two Chinese tropical check varieties, had the j4 allele. Two genotypes TGx 1835-10E and TGx 1987-62F from IITA-Nigeria had the j2 allele. I.C.192 from the US had the j1 allele (Table 5, Fig 3). The identification of these alleles confirms studies by Lu et al., (2017) that the multiple loss-of-function J alleles existing in soybean germplasm worldwide may be responsible for delayed flowering and maturity.
 
Association of the haplotypes with the phenotypic traits
 
Fig 3 shows the variability of the phenotypic traits among the haplotypes defined. Haplotypes Hap1 and Hap3 are the most common among the genotypes with high variability in all the six phenotypic traits. Hap 4 is the least common in only two genotypes. Some genotypes with functional J haplotypes had similar flowering days and, in some cases, even higher (Table 3) than those harboring the loss-of-function alleles indicating that the J gene might not be the only gene responsible for the long juvenile trait. The same observations were recorded on the other associated phenotypic traits such as the plant height, number of nodes per plant and number of pods per plant. This confirms the polygenic nature of the long juvenile trait by previous studies (Ray et al., 1995, Carpentieri-Pipolo et al., 2000, Cober and Morrison 2010, Lu et al., 2017, Yue et al., 2017 and Fang et al., 2020).
 
The present study mainly focused on identifying new long juvenile genotypes and verifying whether the J mutations are associated with delayed flowering in soybean. The long juvenile genetic mechanism in most tropical genotypes especially Sub-Sahara Africa is not known (Miranda et al., 2020). Huaxia 3, the check variety, was introduced to some African countries and resulted in an increased productivity (Yue et al., 2017). Introgression of the loss-of-function alleles Hap j2 and Hap j4 to temperate cultivars has also resulted in the development of new long juvenile varieties in Brazil and southern China (Lu et al., 2017), indicating that the loss-of-function J alleles will be a good strategy to enhance the adaptation and yield increases in the tropics. Therefore, there exists a possibility that incorporation of the identified long juvenile genotypes with or without the loss of function alleles for the J gene in tropical soybean breeding programs will lead to yield increases. However, proper yield trials have to be conducted first, as past studies reported background effects could influence delayed flowering in soybean (Ray et al., 1995). The lack of proper yield data is a key limitation to the study. Future studies will focus on comprehensive yield trials on the identified long juvenile genotypes.
In this study, we screened 159 tropical soybean genotypes from tropical regions worldwide. We have identified 53 soybean genotypes with the delayed flowering. Haplotype analysis revealed six of them carried the loss-of-function alleles of the J gene. We also observed not all the genotypes with delayed flowering had the loss of function alleles for the J gene. This is indicative that some other genes may be acting with the J gene. The information we provided and the new long juvenile germplasm resources identified will be useful for tropical soybean breeding programs. Future studies should be directed towards mining other genes that may also be responsible for the long juvenile trait.
The authors acknowledged the support provided by Major Science and Technology Projects of China (2016ZX08004-003), China Agriculture Research System (CARS-04).

  1. Bonato, E.R. and Vello, N.A.J.G. (1999). E6, a dominant gene conditioning early flowering and maturity in soybeans. Journal of Genetics and Molecular Biology. 22: 229-232.

  2. Carpentieri-Pipolo, V., de Almeida, L.A., Kiihl, R.A.D. and Rosolem, C.A. (2000). Inheritance of long juvenile period under short day conditions for the BR80-6778 soybean [Glycine max (L.) Merrill] line. Euphytica. 112: 203-209.

  3. Cober, E.R. and Morrison, M.J. (2010). Regulation of seed yield and agronomic characters by photoperiod sensitivity and growth habit genes in soybean. Theoretical and Applied Genetics. 120: 1005-1012.

  4. Destro, D., Carpentieri-Pípolo, V., Kiihl, R.d.S. and de Almeida, L.A. (2001). Photoperiodism and genetic control of the long juvenile period in soybean: A review. Crop Breeding and Applied Biotechnology. 1: 72-92.

  5. Fang, X., Han, Y., Liu, M., Jiang, J., Li, X., Lian, Q., Xie, X., Huang, Y., Ma, Q., Nian, H., et al. (2020). Modulation of evening complex activity enables north-to-south adaptation of soybean. Sci China Life Sci. 64: 179-195.

  6. Fehr, W., Caviness, C., Burmood, D. and Pennington, J. (1971). Stage of development descriptions for soybeans, Glycine max (L.) Merrill. Crop Science. 11: 929-931.

  7. Hartwig, E.E. and Kiihl, R.A. (1979). Identification and utilization of a delayed flowering character in soybeans for short-day conditions. Field Crops Research. 2: 145-151.

  8. Hymowitz, T. (1970). On the domestication of the soybean. Economic Botany. 24: 408-421.

  9. Lu, S.J., Zhao, X.H., Hu, Y.L., Liu, S.L., Nan, H.Y. et al. (2017). Natural variation at the soybean J locus improves adaptation to the tropics and enhances yield. Nature Genetics. 49: 773-779.

  10. Miranda, C., Scaboo, A., Cober, E., Denwar, N., Bilyeu, K. (2020). The effects and interaction of soybean maturity gene alleles controlling flowering time, maturity and adaptation in tropical environments. BMC Plant Biology. 20: 65.

  11. Neumaier, N. and James, A. (1993). Exploiting the long-juvenile trait to improve adaptation of soybeans to the tropics. Food Legume Newsletters. 8.

  12. Ray, J.D., Hinson, K., Mankono, J.E.B. and Malo, M.F.J.C.S. (1995). Genetic control of a long juvenile trait in soybean. Crop Science. 35: 1001-1006.

  13. Yang, W., Wu, T., Zhang, X., Song, W., Xu, C., Sun, S., et al (2019). Critical Photoperiod Measurement of Soybean Genotypes in Different Maturity Groups. Crop Science. 59: 2055-2061.

  14. Yue, Y.L., Liu, N.X., Jiang, B.J., Li, M., et al. (2017). A single nucleotide deletion in J encoding GmELF3 confers long juvenility and is associated with adaption of tropic soybean. Molecular Plant. 10: 656-658.

Editorial Board

View all (0)