1. Introduction
Soybeans are one of the world’s most essential crops for oil and vegetable protein (
Palmitic acid (16:0) is one of two primary saturated fatty acids in soybeans, accounting for about 11% of soybean oil [9]. Increased palmitic acid content can improve the stability of soybean oil, whereas lower palmitic acid content can reduce the risk of cardiovascular diseases in humans [10]. Traditional breeding and recent advances in genetic engineering have been used to obtain soybean cultivars with low polyunsaturated fatty acid content and improved oxidative stability without chemical hydrogenation [11,12]. Modification of genes related to palmitic acid content has been reported to develop a suitable germplasm with palmitic acid content (>4 to <40%) that can meet health requirements [10]. In recent years, near-infrared spectroscopy has been used to detect fatty acid components in seeds without damage, especially in agriculture, industrial production, and scientific research [13,14,15,16]. To obtain a more reasonable fatty acid composition of soybean oil for industrial use, and due to greater consumer awareness of health issues, QTL mapping, gene cloning, and identifying molecular markers for palmitic acid content have been used to develop new soybean cultivars with desirable fatty acid compositions through molecular marker-assisted breeding of soybeans.
However, variation in the palmitic acid content of soybeans is influenced by a complex network of genetic and environmental factors [3,17,18]. Palmitic acid slightly affects yield-related phenotypes [19]. The current population size and stability used for palmitic acid QTL mapping are small, and the marker density of the constructed genetic linkage map needs to be increased. In previous studies, only 37 palmitic acid-related QTLs of linkage disequilibrium (LD) and 47 QTLs of a genome-wide association study (GWAS) were mapped (
In addition to palmitic acid, overall fatty acid content is a quantitative trait that depends on the combined effects of several major and minor genes [5,17,24]. QTL mapping is an effective method used to reveal the genetic basis of fatty acid formation and, to date, many QTLs for fatty acid content have been detected. However, due to the relatively low density of genetic maps, these QTLs span a considerable genomic region. The relatively low accuracy of QTL mapping using these maps has limited, not only the identification of fatty acid biosynthesis and regulatory networks, but also the application of these QTLs in soybean marker-assisted selection (MAS) breeding efforts [24].
With the publication of the soybean reference genome and continuous development of molecular marker technology, single nucleotide polymorphisms (SNPs) have become a critical component of genetic map construction, map-based cloning, and genome-wide association analysis, and can be used to locate important soybean pathways for complex quantitative traits [25,26]. Researchers have recently used SLAF-seq, the Soybean SNP BeadChip, SNP genotyping and bin map, and other methods to identify genotypes and construct genetic linkage maps. The number of genetic map markers range from 2086 to 8691, and the total genetic lengths are 1478.86–3780.98 cM. The average map distance is 0.4–1.3 cM, and significant progress has been made in important traits of QTL mapping, such as yield [27,28,29,30,31,32,33], quality [24,28,34,35,36], and stress tolerance [37,38,39,40,41]. Specific-locus amplified fragment sequencing (SLAF-seq) is a rapid and large-scale SNP genotyping method based on next-generation sequencing (NGS) technology [42,43], with effective reads and high throughput. It is widely used in genetic map construction, quantitative trait locus (QTL) mapping, and other fields of population genetic research related to molecular breeding [37]. Compared with the difficulty and high costs of chemical detection, molecular marker-assisted breeding is a simple and reasonable approach to improving the palmitic acid content of soybean cultivars [24]. The continuous improvement of soybean genetic mapping based on molecular markers has dramatically improved the speed and accuracy of QTL mapping for many critical agronomic traits and is a powerful method for genetic map construction, QTL mapping, and molecular marker-assisted selection.
In this study, QTL and linkage markers associated with palmitic acid components were obtained and finely mapped based on the phenotypes over three years, and a soybean high-density genetic linkage map was constructed by SLAF-seq technology. The genome and transcriptome were combined to obtain the candidate genes. The high-density genetic map from our research could further explore more effective and novel QTLs in soybeans and identify candidate genes and associated markers that could contribute to the gene cloning, functional verification, and improvement of palmitic acid composition.
2. Results
2.1. Phenotypic Evaluation
The palmitic acid content of the female parent, Dongnong42, (DN42, Heilongjiang province of China, official registration number: ZDD22643) was significantly higher than that of the male parent, Hobbit, (OH, USA, Chinese official registration number: WDD00573). Over three years, the palmitic acid content of the offspring conformed to a normal distribution, and palmitic acid content within the population varied widely (Figure 1, Table 1). The phenotypes of palmitic acid over the three years—2017 Sanya, 2018 Harbin, and 2019 Suihua—were strongly correlated (Table S1).
2.2. SLAF Library Construction and Genotyping of the RIL Population
Before library construction, sequencing, and genotyping of SLAF-seq, the enzyme combination, RsaI and HaeIII, was used to digest the genomic DNA of the two parents and the 181 RIL5 population. The sequencing depths of ‘DN42’ and ‘Hobbit’ were 120.52× and 117.77×, and the average sequencing depth of the offspring was 22.06×. The enzyme digestion efficiency was 86.62%. A total of 460,888 sequencing reads (Figure 2A), of which 73,579 were polymorphic (Figure 2B), were compared with the soybean genome (
2.3. Construction of the Soybean Genetic Map
A total of 9980 out of 48,082 SLAF markers were used to map out 20 chromosomes. The high-density genetic map had a total genetic distance of 2602.58 cM and an average genetic distance between markers of 0.39 cM, respectively (Figure 3A, Table 2). The collinearity analysis of the position and genetic map of the SLAF markers on the reference genome evaluated the quality of the high-density genetic map (Figure 3B). The average Spearman coefficient was 94.8%, showing the reliable consistency of SLAF markers and the soybean genome (Table 2, Figure 2, Table S2).
2.4. QTL Mapping of the RIL Population
The ICIM method was used to map the QTL of palmitic acid content in three environments. Two QTL loci related to palmitic acid were co-located (Table 3) on chromosomes 10 and 15, including three loci, explaining 4.2–10.1% of the phenotypic variance for palmitic acid across environments (Table 3). In Sanya, 2017, the qPA10-1 was mapped between two markers, Marker7632755 and Marker7723217. The odds and the additive effect logarithm were 2.3 and 4.2%, respectively (Figure S2). From 2018 to 2019, we simultaneously located two QTLs, qPA15-1 and qPA15-2, between Marker5381863 and Marker5394023 on chromosome 15, explaining 10.066 and 9.526% of the phenotypic variance, respectively. The logarithms of odds were 6.528 and 3.999, and the additive effects were 0.259 and 0.261, respectively (Figures S3 and S4), which showed these QTLs were associated with a high palmitic acid content.
2.5. Verification of Candidate Region and Screening of Linkage Markers
Based on the QTL mapping results of SLAF-seq, we screened 10 polymorphic markers from 27 SSR markers (Table S3) related to the candidate region with SLAF-seq to verify the region and screen associated markers. The region between two markers, G420 and G483, was confirmed as the candidate region, and the genotypes of two markers, G420 and G483, were significantly correlated with palmitic acid content (Figure 4).
2.6. Analysis of Genome and Transcriptome
After Illumina sequencing and alignment with the reference genome, Glycine max Wm82.a2.v1, 48.84 G of raw data were generated by re-sequencing, of which the raw data yields of DN42 and Hobbit were 25.48 and 23.36 G, respectively, and the filtered clean data were 25.2 and 23.1 G, respectively. The error rate (%) was 0.03, Q30 values (percentage of bases with a Phred value greater than 30 in the total bases) were 94.17 and 94.23%, respectively, and GC contents were 37.5 and 36.99%, respectively. Based on Wm82.a2.v1, seven candidate protein-encoding genes (Glyma.15g119200–Glyma.15g119800) were present in qPA15-1 and qPA15-2. The annotations of candidate genes were all bifunctional inhibitors/lipid-transfer protein/seed storage 2 S albumin superfamily protein, whereas Glyma.15g119300 was not present in Glycine max Wm82.a4.v1 (
To further select candidate genes, the developmental seeds of the six stages (R5–R6) of DN42 and Hobbit were used for seed transcriptome analysis. After the Q-PCR accurately quantified the effective concentrations, the libraries were sequenced on the Illumina platform. After transcriptome analysis, a total of 289.52 G of clean data was obtained; the clean data of each sample reached 5.48 G and the percentage of Q30 exceeded 92.31%. The efficiencies of clean reads of each sample sequence with the reference genome ranged from 91.11 to 95.19%. Based on the transcriptome, we found that three genes, Glyma.15g119200, Glyma.15g119300, and Glyma.15g119400, were not expressed in the two parents, and the expression levels of Glyma.15g119500 in Hobbit and Glyma.15g119600 in DN42 were deficient. The expression level of Glyma.15g119700 in DN42 was significantly higher than that of Hobbit at the six periods from R5 to R6. From R5 to R5.4, the expression level of Glyma.15g119800 in Hobbit was higher than that of DN42, but the difference was not significant in the early stages, and became significant at R5.4. On the contrary, the expression level in Hobbit was significantly lower than that of DN42 at R6 (Figure 5).There was a high similarity in the coding sequence among four genes: Glyma.15g119500, Glyma.15g119600, Glyma.15g119700, and Glyma.15g119800 (Figure S5).
3. Discussion
3.1. Construction of the High-Density Genetic Map
Previous studies have shown that relatively low-density genetic maps are not sufficient for the quantitative analysis of trait genes controlled by multiple genes, especially genes with low effects and high susceptibility to the environment [3,44]. According to our statistics, since 2010, about 40 genetic maps have been developed from soybean research on the yield, quality, and resistance to biotic or abiotic stresses (references not shown), with 2086–8691 markers of genetic populations, and total lengths of 1478.86–3780.98 cM. Only six of these maps had higher densities, exceeding 6000 markers [18,37,38,40,45].
In this study, we exploited 181 RIL5 populations to develop 460,888 SLAF markers, and the effective parental polymorphism was 10.43%. The polymorphism SLAF tag of populations was 73,579, located in 20 linkage groups. Finally, 9980 SLAF markers were constructed on the genetic map, with the total genetic distance and average map distances being 2602.58 cM and 0.39 cM, respectively. Compared with previous studies, the marker density of the genetic map in this study was still very high. SoySNP50K was used to detect two populations to construct two genetic maps with 21,478 and 11,922 markers, and, with the updating of the Glyma1.01 build based on two high-resolution linkage maps, the genomic positions of the commonly used markers in the BARCSOYSSR_1.0 database and SoySNP50K BeadChip were updated, based on the Wm82.a2.v1 assembly [46]. The high-density genetic map from this study could facilitate the identification of genes or QTL controlling yield, seed quality, and tolerance of biotic or abiotic stress, as well as other factors.
3.2. QTL Mapping of Palmitic Acid Content
The fatty acid content of soybean has high variability, and the palmitic acid content of inbred lines is significantly higher than that of landraces, and is greatly affected by the environment [4,18]. Previous studies have shown that reduced palmitic acid content may lead to increased oil content [47]. Currently, numerous QTLs related to palmitic acid have been identified, with 39 QTLs identified on 15 of 20 chromosomes and 47 reported SNPs (
Based on our genetic map, we located a QTL on chromosome 10, explaining 4.151% of the phenotypic variance in 2017 Sanya. This QTL was not repeatedly located in the following two years and may be an environment-specific QTL. Across two environments, qPA15-1 and qPA15-2 on chromosome 15 were repeatedly identified, with LOD values of 6.5 and 4.0, and explained 10.1 and 9.5% of the phenotypic variance, respectively. This QTL may also be affected by the environment, which may be why it was not detected in 2017 Sanya. Based on the QTL mapping region of SLAF-seq, we screened two markers based on 27 SSR markers, G420 (Gm15:9,358,076) and G483 (Gm15:9,466,527), which were associated with palmitic acid based on the RIL population derived from DN42 × Hobbit. This mapping region of chromosome 15 has been reported previously [53], shown as Seed palmitic 7-3 (
3.3. Analysis of Candidate Genes
High palmitic acid levels could improve the oxidative stability of soybean oil, and reducing palmitic acid in oil (16:0) from 11 to <4% makes it possible to achieve a low content of saturated fatty acids, which is desirable for cardiovascular health [11]. In soybeans, significant efforts have been made to identify and characterize mutant alleles related to palmitic acid content to improve oil composition [57,58]. For example, mutations at GmFATB (Glyma05g08060 (Wm82.a1.v1) and Glyma.05g012300 (Wm82.a2.v1)) present low palmitic acid (as low as 5.6%) [10]. New alleles of GmFATB1A (Glyma.05g012300) elongate FATB with an acyl-ATP thioesterase, which could reduce palmitic acid levels [59]. Glyma.14g27990 (Glyma.14g121400), coding Stearoyl-ACP desaturase isoform C (SACPD-C), could influence palmitic and stearic acid [60]. A beta-ketoacyl-[acyl-carrier-protein] synthase III gene (KAS III / Glyma09g41380, Glyma. 09g 277400) was highly correlated with the reduced palmitic acid phenotype due to fap1 [20]. Candidate genes of Abscisic Acid Insensitive 3 (ABI3, Glyma18g38490, Glyma.18g176100) were significantly associated with palmitic acid concentration, explaining 5.42% of phenotypic variation [61]. However, the small number of genes that have been published are insufficient for the identification and breeding of ideal soybean varieties of specific palmitic acid content. We analyzed the expression levels of the four published palmitic acid-related genes between parents, DN42 and Hobbit. Except for Glyma.09g277400 in the R6 stage and Glyma.18g176100 in the R5 stage, the expression levels of the four genes in the other seed development stages were not significantly different(Figure S6). In the candidate region of our study, Glyma.15g119300 was not present in Glycine max Wm82.a4.v1, which may have been caused by genome splicing errors. Twelve nonsynonymous SNVs existed in the coding regions of the six genes. Interestingly, the identical gene annotation of seven genes was consistent: bifunctional inhibitors/lipid-transfer protein/seed storage 2 S albumin superfamily protein. After transcriptome analysis, we found that Glyma.15g119200, Glyma.15g119300, and Glyma.15g119400 were not expressed in the two parents; and the expression level of Glyma.15g119500 in Hobbit and Glyma.15g119600 in DN42 were deficient. Furthermore, through the genomic information of DN42 and Hobbit, and we found that there were two and five non-synonymous mutations in Glyma.15g119700 and Glyma.15g119800, respectively, so we reasoned that these two genes were the candidate genes for the palmitic acid content. Four genes with high similarity in coding sequence, Glyma.15g119500, Glyma.15g119600, Glyma.15g119700, and Glyma.15g119800 (Figure S5), may act as a repeating sequence. However, further research is required to confirm candidate genes and mechanisms. Identifying QTLs, linkage markers, and candidate genes may help improve soybean oil composition traits to further enhance soybean oil quality.
4. Materials and Methods
4.1. Plant Material and Phenotyping
Soybean RIL population lines and their parents, Dongnong42 (DN42, Heilongjiang province of China, official registration number: ZDD22643) and Hobbit (OH, USA, Chinese official registration number: WDD00573), were grown in three different environments (Evr.1–Evr.3), including Sanya (Evr.1, Sanya city, 18°30′ N 109°52′ E, China) in 2017, Harbin (Evr.2, Harbin city, 45°83′ N 126°84′ E, China) in 2018, and Suihua (Evr.3, Suihua city, 46°63′ N 126°98′ E, China) in 2019. These locations belonged to two soybean ecotopes (Northern and Southern eco-regions) in China, and the soils were Northeast black soil and loam-type soil. Sanya’s climate and soil fertility differed from the other two locations. The irrigation and fertilization regimes of the exact location in 2018 and 2019 were similar.
In addition, the other field management practices were the same as the general practices of the local region. Randomized complete block designs were used for all experiments with three rows, 4 m long and 0.60 m apart, with a space of 6.67 cM between plants. Three replications were performed.
The soybean seeds were matured and dried to a constant weight before analyzing palmitic acid content. Measurement of soybean seeds was conducted using a near-infrared spectrophotometer (NIR) seed analyzer (FOSS NIR DS2500, FOSS Instrument, Denmark). Briefly, approximately 50 g of dry seeds per sample was fitted into a 50 mm diameter cup, and the cup was rotated during NIRS scanning. The average value of three scans per sample was used in data analysis.
4.2. Genotype Analysis, Construction, and Genotyping of the SLAF Library, and Construction of the High-Density Genetic Map
The DNA of young, healthy, fresh leaves from both parent lines and all 181 RILs were extracted using the modified CTAB protocol [62]. The specific-locus amplified fragment sequencing library of the two parents and 181 RILs was constructed and sequenced [42]. Based on the soybean reference genome (
4.3. QTL Analysis for Palmitic Acid Content
According to the phenotypes of the three environments, QTLs for palmitic acid content were detected using inclusive composite interval mapping (ICIM) in the R/qtl package [63,64,65]. One thousand permutation tests at a 95% confidence level were used to set the LOD threshold. Based on 1000 permutations, LOD = 2.5 was used to determine the presence of a putative QTL associated with a target trait in a particular genomic region.
4.4. Validation of Candidate Region and Development of Linkage Marker
In Harbin, 2020, young leaves of 30 plants selected from all population lines were collected and mixed in equal amounts, 30 days after the F10 generation plants were sown at the Nangang experimental station (Harbin, China). Genomic DNA of the two parents and RIL population were extracted by the CTAB method for genotype determination.
Based on the QTL mapping results of SLAF-seq, we screened the polymorphic markers of 27 SSR markers (Table S3) [46] related to the region for verification and linkage markers, which were well-distributed inside and outside of the region. Genomic DNA from parents and 260 RILs were used as templates. The 20-μL PCR contained 100 ng of template DNA, 1.5 mmol/L primers (Tsingke Biotech.), 1 × PCR buffer, 2 mmol/L dNTPs, and 1 unit of Taq polymerase (TransGen Biotech.). The PCR amplification program was as follows: 95 °C for 5 min, 34 cycles of 95 °C denaturations for 30 s, 58 °C annealing for 30 s, 72 °C extensions for 30 s, and a final 72 °C extension for 5 min.
A total of 3 uL 6 × Loading Buffer (loading buffer containing bromophenol blue indicator) was added to a 20 uL sample of PCR amplification product and mixed, then 2 uL of prepared PCR amplification was detected using 6% polyacrylamide gel electrophoresis (PAGE) under 120 V, 400 mA, lasting 30 min. The gel was placed and shaken in 1% AgNO3 solution for 3 min and then washed with water; the genotype band was visualized with a developing solution (1000 mL ddH2O, 1.5 g NaOH, and 5 mL formaldehyde (35%), mixed). According to general genetic coding rules, the genotypes of the RIL population were assigned as A and B, which were consistent with DN42, Hobbit, and H (heterozygous genotype).
QTL IciMapping was used to validate the candidate region-based SSR marker.
4.5. Re-Sequencing and RNA-Seq Analysis of Parents
Re-sequencing and RNA-seq of the two parents were accomplished at Beijing Bio Marker Bioinformatics Technology Co., Ltd. (Beijing, China). For re-sequencing, the DNA of young, healthy, and fresh leaves from both parent lines were extracted using the modified CTAB protocol. Libraries for the four DNA pools were prepared according to Illumina’s library preparation protocol. The nucleic acid sequences of Wm82.a2.v1 (
For RNA-seq, developing seeds of the four stages were collected between R5 (beginning seed stage; seed is 1/8 inch long in one of the four uppermost nodes on the main stem; primary and lateral roots grow strongly until R5) and R6 (full seed stage; pods contain a green seed that fills the cavity in one of the four upper most nodes on the main stem; most nutrients have been taken up by the time the plant reaches the R6 stage) [66] from the two parents, DN42 and Hobbit. All tissues were immediately frozen in liquid nitrogen and then stored at −80 °C. Each sample was collected from at least five independent plants and pooled together. According to the supplier’s instructions, total RNA was isolated from the tissues using the RNA prep Pure Plant Kit (TIANGEN, Beijing, China). RNA concentration and quality were measured by spectrophotometer (NanoDrop 2000, Thermo Scientific, waltham, MA, USA). The integrity of RNA was assessed by electrophoresis on a 1% agarose gel and observed by Good View TM staining. First-strand complementary cDNA was synthesized from 2 µg total RNA using the Revert Aid First Strand cDNA Synthesis Kit (Thermo Scientific, Waltham, MA, USA). Each cDNA was diluted 1:10 by adding RNase-free ddH2O and stored at −80 °C. qRT-PCR was performed on a CFX96 Touch™ Cycler (Bio-Rad, Hercules, CA, USA) system using SYBR Premix Ex Taq TM II (Takara Biomedical Technology, Beijing, China). Three independent biological replicates were performed for each tissue. Sequencing libraries were generated using NEB Next UltraTM RNA Library Prep Kit for Illumina (San Diego, USA), following the manufacturer’s recommendations, and sequenced on an Illumina platform. Differential expression analysis of groups was performed using DESeq2. The p-values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate. Genes with an adjusted p-value < 0.01 found by DESeq2 were assigned as being differentially expressed. RNA-seq analysis was performed using BMK Cloud (
4.6. Statistical Analysis
Microsoft Excel 2016 and SPSS19.0 were used for statistical analysis of phenotypic data and correlation analysis. One-way variance (ANOVA) analysis of SSR markers and expression analysis were carried out using the aov function in R.
5. Conclusions
In this study, a high-density genetic map related to palmitic acid content was constructed, and two QTLs were mapped based on the genetic map. SSR markers that could be helpful for breeding low palmitic acid content cultivars were developed. Two candidate genes of palmitic acid content were screened.
Y.X., performed the expression analysis and RIL population plants; H.G. and X.L. (Xinlei Liu) performed the SLAF-Seq and phenotype observations; L.Z. and X.L. (Xiaoyan Luan) completed linkage map construction and candidate genes analysis and verification; X.T. and D.C. performed the qPCR and SSR marker PCR trial and data analysis; L.Z. and X.L. (Xinlei Liu) performed genotype analysis; L.Q., Y.X. and H.G. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
Sequence data from this article can be found in the Phytozome database under the following accession numbers: Glyma.15g119500, Glyma.15g119600, Glyma.15g119700, and Glyma.15g119800.
We thank Qingshan Chen (Northeast Agricultural University) for sharing cultivar Dongnong42.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure 1. The frequency distribution of palmitic acid content (x-axis, %) of the RIL population and parents in (A) 2017, (B) 2018, (C) and 2019. F: female parent (Dongnong42, DN42); M: male parent (Hobbit).
Figure 2. Distribution of SLAF and polymorphic SLAF markers in the genome. (A) The picture shows the distribution of SLAF labels. (B) Distribution of polymorphic SLAF labels. The abscissa is the length of the chromosome. Each yellow band represents a chromosome. The genome is divided according to the size of 1 M. The darker areas in the figure are the areas where the SLAF tags are concentrated.
Figure 3. High density genetic map and collinearity analysis of the genetic map and genome. Distribution of markers on 20 chromosomes. (A) Black bars on each chromosome represent mapped SLAF-seq markers. The chromosome number is shown on the x-axis, and genetic distance (cM) is shown on the y-axis. (B) The x-axis represents the genetic distance of each linkage group. The y-axis represents the physical length of each linkage group, with the collinearity of genomic markers and genetic maps represented in scatter. Different colors represent different chromosomes or linkage groups.
Figure 4. Validation of candidate region and analysis of linkage markers. (A) Ten polymorphic markers related to the candidate region with SLAF-seq used to verify the region and screen a QTL region between G420 and G483 markers. (B) The lines of the RIL population were divided into two groups based on G420 genotypes. The palmitic acid content of the two groups was significantly different. Student’s t-test, * p < 0.05. (C) The lines of the RIL population were divided into two groups based on G483 genotypes. The palmitic acid content of the two groups was significantly different. Student’s t-test, ** p < 0.01. The “A” and “B” represent the genotypes of DN42 and Hobbit, respectively.
Figure 5. The relative expression (y-axis) of the four genes in the candidate region between the growth period of R5–R6 (x-axis) in RNA-seq. (A) Expression level of Glyma.15g119500 at the seed development R5–R6 stage. (B) Expression level of Glyma.15g119600 at the seed developing R5–R6 stage. (C) Expression level of Glyma.15g119700 at the seed development R5–R6 stage. (D) Expression level of Glyma.15g119800 at the seed development R5–R6 stage. Bars mean SD. *, p < 0.05; **, p < 0.01.
Statistical analysis of palmitic acid content (%) in the parent and RIL population over 3 years.
Year | Parents | RILs | ||||||
---|---|---|---|---|---|---|---|---|
DN42 | Hobbit | Mean | SD | Range | Skewness | Kurtosis | CV (%) | |
2017 | 11.99 | 9.32 | 9.97 | 0.8 | 1.69–11.2 | −6.17 | 61.25 | 0.08 |
2018 | 11.99 | 8.12 | 9.52 | 0.76 | 1.57–11.03 | −6.17 | 63.94 | 0.08 |
2019 | 10.86 | 8.92 | 9.92 | 0.81 | 1.9–11.69 | −5.25 | 51.7 | 0.08 |
Characteristics of the high-density genetic map.
Chromosome | Number of Markers | Map Distance (cM) | Average Map Distance (cM) | Gap < 5 |
Max |
---|---|---|---|---|---|
1 | 539 | 184.4 | 0.34 | 100 | 4.53 |
2 | 170 | 102.71 | 0.6 | 99.41 | 6.15 |
3 | 754 | 106.6 | 0.14 | 100 | 2.93 |
4 | 136 | 156.16 | 1.15 | 100 | 3.24 |
5 | 319 | 117.18 | 0.37 | 99.06 | 12.36 |
6 | 356 | 137.53 | 0.39 | 99.72 | 5.46 |
7 | 401 | 131.11 | 0.33 | 99.75 | 5.23 |
8 | 190 | 128.56 | 0.68 | 100 | 3.26 |
9 | 420 | 96.8 | 0.23 | 100 | 3.24 |
10 | 269 | 117.89 | 0.44 | 98.88 | 12.47 |
11 | 145 | 116.45 | 0.8 | 99.31 | 5.2 |
12 | 134 | 98.46 | 0.73 | 96.99 | 7.57 |
13 | 684 | 146.24 | 0.21 | 99.85 | 6.49 |
14 | 699 | 112.55 | 0.16 | 99.71 | 8.86 |
15 | 414 | 149.5 | 0.36 | 100 | 4.95 |
16 | 845 | 123.57 | 0.15 | 100 | 2.93 |
17 | 969 | 136.54 | 0.14 | 100 | 4.53 |
18 | 1119 | 169.34 | 0.15 | 100 | 2.43 |
19 | 1093 | 158.84 | 0.15 | 99.73 | 11.84 |
20 | 324 | 112.15 | 0.35 | 100 | 4.53 |
Total | 9980 | 2602.58 | 0.39 | 19.62 | 5.91 |
QTLs of the palmitic acid content of the RIL population in the period 2017–2019.
QTL | Year | Chr. | Left Marker | Right Marker | Genetic Position (cM) | Physical Position (bp) | LOD | ADD | PVE |
||
---|---|---|---|---|---|---|---|---|---|---|---|
Start | End | Start | End | ||||||||
qPA10-1 | 2017 | 10 | Marker7632755 | Marker7723217 | 76.341 | 76.341 | 39,427,048 | 39,517,256 | 2.298 | 0.174 | 4.151 |
qPA15-1 | 2018 | 15 | Marker5381863 | Marker5394023 | 43.494 | 43.78 | 9,218,139 | 9,387,049 | 6.528 | 0.259 | 10.066 |
qPA15-2 | 2019 | 15 | Marker5381863 | Marker5394023 | 43.494 | 43.78 | 9,218,139 | 9,387,049 | 3.999 | 0.261 | 9.526 |
Chr.: chromosome, LOD: logarithm of odds, PVE: phenotypic variance explained, ADD: additive effect.
SNPs and Indels in the candidate regions of Dongnong42 and Hobbit.
Gene | SNP | Indel | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Intergenic | Upstream | Downstream | Nonsynonymous |
Synonymous |
Intronic | 3′UTR | Intergenic | Upstream | Downstream | Non-Frameshift |
Non-Frameshift |
5′ UTR | |
Glyma.15g119200 | 5 | 10 | 2 | 3 | |||||||||
Glyma.15g119300 | 3 | 8 | 1 | 2 | 1 | ||||||||
Glyma.15g119400 | 3 | 1 | 1 | 4 | 3 | 1 | |||||||
Glyma.15g119500 | 7 | 6 | 1 | 3 | 4 | 7 | 1 | 3 | |||||
Glyma.15g119600 | 14 | 23 | 2 | 5 | |||||||||
Glyma.15g119700 | 15 | 13 | 2 | 1 | 4 | 3 | 1 | ||||||
Glyma.15g119800 | 21 | 20 | 5 | 8 | 3 | ||||||||
Total | 319 | 68 | 80 | 12 | 1 | 1 | 11 | 77 | 17 | 23 | 1 | 1 | 5 |
SNP and Indel numbers were contained in 2000 bp region upstream and downstream from the gene.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Crapiste, G.H.; Brevedan, M.I.; Carelli, A.A. Oxidation of sunflower oil during storage. J. Am. Oil Chem. Soc.; 1999; 76, 1437. [DOI: https://dx.doi.org/10.1007/s11746-999-0181-5]
2. Sala-Vila1, A.; Ros, E. Mounting evidence that increased consumption of a. linolenic acid, the vegetablen. 3 fatty acid, may benefit cardiovascular health. Clin. Lipidol.; 2011; 6, pp. 365-369. [DOI: https://dx.doi.org/10.2217/clp.11.36]
3. Fan, S.; Li, B.; Yu, F.; Han, F.; Yan, S.; Wang, L.; Sun, J. Analysis of additive and epistatic quantitative trait loci underlying fatty acid concentrations in soybean seeds across multiple environments. Euphytica; 2015; 206, pp. 689-700. [DOI: https://dx.doi.org/10.1007/s10681-015-1491-3]
4. Abdelghany, A.M.; Zhang, S.; Azam, M.; Shaibu, A.S.; Feng, Y.; Li, Y.; Tian, Y.; Hong, H.; Li, B.; Sun, J. Profiling of seed fatty acid composition in 1025 Chinese soybean accessions from diverse ecoregions. Crop. J.; 2020; 8, pp. 635-644. [DOI: https://dx.doi.org/10.1016/j.cj.2019.11.002]
5. Panthee, D.; Pantalone, V.; Saxton, A. Modifier QTL for fatty acid composition in soybean oil. Euphytica; 2006; 152, pp. 67-73. [DOI: https://dx.doi.org/10.1007/s10681-006-9179-3]
6. Hu, F.B.; Stampfer, M.J.; Manson, J.E.; Rimm, E.; Colditz, G.A.; Rosner, B.A.; Hennekens, C.H.; Willett, W.C. Dietary fat intake and the risk of coronary heart disease in women. N. Engl. J. Med.; 1997; 337, pp. 1491-1499. [DOI: https://dx.doi.org/10.1056/NEJM199711203372102] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9366580]
7. Thelen, J.J.; Ohlrogge, J.B. Metabolic engineering of fatty acid biosynthesis in plants. Metab. Eng.; 2002; 4, pp. 12-21. [DOI: https://dx.doi.org/10.1006/mben.2001.0204]
8. Canakci, M.; Monyem, A.; Van Gerpen, J. Accelerated oxidation processes in biodiesel. Trans. ASAE; 1999; 42, 1565. [DOI: https://dx.doi.org/10.13031/2013.13321]
9. Graef, G.; LaVallee, B.J.; Tenopir, P.; Tat, M.; Schweiger, B.; Kinney, A.J.; Van Gerpen, J.H.; Clemente, T.E. A high-oleic-acid and low-palmitic-acid soybean: Agronomic performance and evaluation as a feedstock for biodiesel. Plant Biotechnol. J.; 2009; 7, pp. 411-421. [DOI: https://dx.doi.org/10.1111/j.1467-7652.2009.00408.x]
10. Zhou, Z.; Lakhssassi, N.; Knizia, D.; Cullen, M.A.; El Baz, A.; Embaby, M.G.; Liu, S.; Badad, O.; Vuong, T.D.; AbuGhazaleh, A. et al. Genome-wide identification and analysis of soybean acyl-ACP thioesterase gene family reveals the role of GmFAT to improve fatty acid composition in soybean seed. Theor. Appl. Genet.; 2021; 134, pp. 3611-3623. [DOI: https://dx.doi.org/10.1007/s00122-021-03917-9]
11. Fehr, W.R. Breeding for modified fatty acid composition in soybean. Crop Sci.; 2007; 47, pp. S72-S87. [DOI: https://dx.doi.org/10.2135/cropsci2007.04.0004IPBS]
12. Pham, A.-T.; Bilyeu, K.; Chen, P.; Boerma, H.R.; Li, Z. Characterization of the fan1 locus in soybean line A5 and development of molecular assays for high-throughput genotyping of FAD3 genes. Mol. Breed.; 2014; 33, pp. 895-907. [DOI: https://dx.doi.org/10.1007/s11032-013-0003-1]
13. Pazdernik, D.L.; Killam, A.S.; Orf, J.H. Analysis of amino and fatty acid composition in soybean seed, using near infrared reflectance spectroscopy. Agron. J.; 1997; 89, pp. 679-685. [DOI: https://dx.doi.org/10.2134/agronj1997.00021962008900040022x]
14. Naeve, S.L.; Proulx, R.A.; Hulke, B.S.; O’Neill, T.A. Sample size and heterogeneity effects on the analysis of whole soybean seed using near infrared spectroscopy. Agron. J.; 2008; 100, pp. 231-234. [DOI: https://dx.doi.org/10.2134/agronj2007.0230]
15. Ambrose, A.; Lohumi, S.; Lee, W.-H.; Cho, B.K. Comparative nondestructive measurement of corn seed viability using Fourier transform near-infrared (FT-NIR) and Raman spectroscopy. Sens. Actuators B; 2016; 224, pp. 500-506. [DOI: https://dx.doi.org/10.1016/j.snb.2015.10.082]
16. Karn, A.; Heim, C.; Flint-Garcia, S.; Bilyeu, K.; Gillman, J. Development of rigorous fatty acid near-infrared spectroscopy quantitation methods in support of soybean oil improvement. J. Am. Oil Chem. Soc.; 2017; 94, pp. 69-76. [DOI: https://dx.doi.org/10.1007/s11746-016-2916-4]
17. Rebetzke, G.; Pantalone, V.; Burton, J.; Carter, T.; Wilson, R. Genetic background and environment influence palmitate content of soybean seed oil. Crop Sci.; 2001; 41, pp. 1731-1736. [DOI: https://dx.doi.org/10.2135/cropsci2001.1731]
18. Song, W.; Yang, R.; Wu, T.; Wu, C.; Sun, S.; Zhang, S.; Jiang, B.; Tian, S.; Liu, X.; Han, T. Analyzing the effects of climate factors on soybean protein, oil contents, and composition by extensive and high-density sampling in China. J. Agric. Food Chem.; 2016; 64, pp. 4121-4130. [DOI: https://dx.doi.org/10.1021/acs.jafc.6b00008] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27022763]
19. Cardinal, A.J.; Dewey, R.E.; Burton, J.W. Estimating the individual effects of the reduced palmitic acid fapnc and fap1 alleles on agronomic traits in two soybean populations. Crop Sci.; 2008; 48, pp. 633-639. [DOI: https://dx.doi.org/10.2135/cropsci2007.05.0251]
20. Cardinal, A.J.; Whetten, R.; Wang, S.; Auclair, J.; Hyten, D.; Cregan, P.; Bachlava, E.; Gillman, J.; Ramirez, M.; Dewey, R. et al. Mapping the low palmitate fap1 mutation and validation of its effects in soybean oil and agronomic traits in three soybean populations. Theor. Appl. Genet.; 2014; 127, pp. 97-111. [DOI: https://dx.doi.org/10.1007/s00122-013-2204-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24132738]
21. De Vries, B.D.; Fehr, W.R.; Welke, G.A.; Dewey, R.E. Molecular characterization of the mutant fap3 (A22) allele for reduced palmitate concentration in soybean. Crop Sci.; 2011; 51, pp. 1611-1616. [DOI: https://dx.doi.org/10.2135/cropsci2010.10.0619]
22. Cardinal, A.J.; Burton, J.W.; Camacho-Roger, A.M.; Yang, J.H.; Wilson, R.F.; Dewey, R.E. Molecular analysis of soybean lines with low palmitic acid content in the seed oil. Crop Sci.; 2007; 47, pp. 304-310. [DOI: https://dx.doi.org/10.2135/cropsci2006.04.0272]
23. Buhr, T.; Sato, S.; Ebrahim, F.; Xing, A.; Zhou, Y.; Mathiesen, M.; Schweiger, B.; Kinney, A.; Staswick, P.; Clemente, T. Ribozyme termination of RNA transcripts down-regulate seed fatty acid genes in transgenic soybean. Plant J.; 2002; 30, pp. 155-163. [DOI: https://dx.doi.org/10.1046/j.1365-313X.2002.01283.x]
24. Li, B.; Fan, S.; Yu, F.; Chen, Y.; Zhang, S.; Han, F.; Yan, S.; Wang, L.; Sun, J. High-resolution mapping of QTL for fatty acid composition in soybean using specific-locus amplified fragment sequencing. Theor. Appl. Genet.; 2017; 130, pp. 1467-1479. [DOI: https://dx.doi.org/10.1007/s00122-017-2902-8]
25. Schmutz, J.; Cannon, S.B.; Schlueter, J.; Ma, J.; Mitros, T.; Nelson, W.; Hyten, D.L.; Song, Q.; Thelen, J.J.; Cheng, J. et al. Genome sequence of the palaeopolyploid soybean. Nature; 2010; 463, pp. 178-183. [DOI: https://dx.doi.org/10.1038/nature08670] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20075913]
26. Soranzo, N.; Spector, T.D.; Mangino, M.; Kühnel, B.; Rendon, A.; Teumer, A.; Willenborg, C.; Wright, B.; Chen, L.; Li, M. et al. A genome-wide meta-analysis identifies 22 loci associated with eight hematological parameters in the HaemGen consortium. Nat. Genet.; 2009; 41, pp. 1182-1190. [DOI: https://dx.doi.org/10.1038/ng.467] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19820697]
27. Wang, W.; Liu, M.; Wang, Y.; Li, X.; Cheng, S.; Shu, L.; Yu, Z.; Kong, J.; Zhao, T.; Gai, J. Characterizing two inter-specific bin maps for the exploration of the QTLs/genes that confer three soybean evolutionary traits. Front. Plant Sci.; 2016; 7, 1248. [DOI: https://dx.doi.org/10.3389/fpls.2016.01248]
28. Cao, Y.; Li, S.; Wang, Z.; Chang, F.; Kong, J.; Gai, J.; Zhao, T. Identification of major quantitative trait loci for seed oil content in soybeans by combining linkage and genome-wide association mapping. Front. Plant Sci.; 2017; 8, 1222. [DOI: https://dx.doi.org/10.3389/fpls.2017.01222]
29. Liu, N.; Li, M.; Hu, X.; Ma, Q.; Mu, Y.; Tan, Z.; Xia, Q.; Zhang, G.; Nian, H. Construction of high-density genetic map and QTL mapping of yield-related and two quality traits in soybean RILs population by RAD-sequencing. BMC Genom.; 2017; 18, pp. 1-13. [DOI: https://dx.doi.org/10.1186/s12864-017-3854-8]
30. Kong, L.; Lu, S.; Wang, Y.; Fang, C.; Wang, F.; Nan, H.; Su, T.; Li, S.; Zhang, F.; Li, X. et al. Quantitative trait locus mapping of flowering time and maturity in soybean using next-generation sequencing-based analysis. Front. Plant Sci.; 2018; 9, 995. [DOI: https://dx.doi.org/10.3389/fpls.2018.00995]
31. Watanabe, S.; Shimizu, T.; Machita, K.; Tsubokura, Y.; Xia, Z.; Yamada, T.; Hajika, M.; Ishimoto, M.; Katayose, Y.; Harada, K. et al. Development of a high-density linkage map and chromosome segment substitution lines for Japanese soybean cultivar Enrei. DNA Res.; 2018; 25, pp. 123-136. [DOI: https://dx.doi.org/10.1093/dnares/dsx043] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29186379]
32. Han, J.; Han, D.; Guo, Y.; Yan, H.; Wei, Z.; Tian, Y.; Qiu, L. QTL mapping pod dehiscence resistance in soybean (Glycine max L. Merr.) using specific-locus amplified fragment sequencing. Theor. Appl. Genet.; 2019; 132, pp. 2253-2272. [DOI: https://dx.doi.org/10.1007/s00122-019-03352-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31161230]
33. Li, Y.; Dong, Y.; Wu, H.; Hu, B.; Zhai, H.; Yang, J.; Xia, Z. Positional cloning of the flowering time QTL qFT12-1 reveals the link between the clock related PRR homolog with photoperiodic response in soybeans. Front. Plant Sci.; 2019; 10, 1303. [DOI: https://dx.doi.org/10.3389/fpls.2019.01303] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31681389]
34. Qi, Z.; Huang, L.; Zhu, R.; Xin, D.; Liu, C.; Han, X.; Jiang, H.; Hong, W.; Hu, G.; Zheng, H. et al. A high-density genetic map for soybean based on specific length amplified fragment sequencing. PLoS ONE; 2014; 9, e104871. [DOI: https://dx.doi.org/10.1371/journal.pone.0104871] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25118194]
35. Cai, Z.; Cheng, Y.; Ma, Z.; Liu, X.; Ma, Q.; Xia, Q.; Zhang, G.; Mu, Y.; Nian, H. Fine-mapping of QTLs for individual and total isoflavone content in soybean (Glycine max L.) using a high-density genetic map. Theor. Appl. Genet.; 2018; 131, pp. 555-568. [DOI: https://dx.doi.org/10.1007/s00122-017-3018-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29159422]
36. Patil, G.; Vuong, T.D.; Kale, S.; Valliyodan, B.; Deshmukh, R.; Zhu, C.; Wu, X.; Bai, Y.; Yungbluth, D.; Lu, F. Dissecting genomic hotspots underlying seed protein, oil, and sucrose content in an interspecific mapping population of soybean using high-density linkage mapping. Plant Biotechnol. J.; 2018; 16, pp. 1939-1953. [DOI: https://dx.doi.org/10.1111/pbi.12929] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29618164]
37. Zhang, D.; Li, H.; Wang, J.; Zhang, H.; Hu, Z.; Chu, S.; Lv, H.; Yu, D. High-density genetic mapping identifies new major loci for tolerance to low-phosphorus stress in soybean. Front. Plant Sci.; 2016; 7, 372. [DOI: https://dx.doi.org/10.3389/fpls.2016.00372]
38. Lee, J.S.; Kim, S.; Ha, B.-K.; Kang, S. Positional mapping and identification of novel quantitative trait locus responsible for UV-B radiation tolerance in soybean [Glycine max (L.) Merr.]. Mol. Breed.; 2016; 36, 50. [DOI: https://dx.doi.org/10.1007/s11032-016-0471-1]
39. Niu, J.; Guo, N.; Sun, J.; Li, L.; Cao, Y.; Li, S.; Huang, J.; Zhao, J.; Zhao, T.; Xing, H. Fine mapping of a resistance gene RpsHN that controls Phytophthorasojae using recombinant inbred lines and secondary populations. Front. Plant Sci.; 2017; 8, 538. [DOI: https://dx.doi.org/10.3389/fpls.2017.00538]
40. Ren, H.; Han, J.; Wang, X.; Zhang, B.; Yu, L.; Gao, H.; Hong, H.; Sun, R.; Tian, Y.; Qi, X. et al. QTL mapping of drought tolerance traits in soybean with SLAF sequencing. Crop J.; 2020; 8, pp. 977-989. [DOI: https://dx.doi.org/10.1016/j.cj.2020.04.004]
41. Wilkes, J.; Saski, C.; Klepadlo, M.; Fallen, B.; Agudelo, P. Quantitative trait loci associated with Rotylenchulusreniformis host suitability in soybean. Phytopathology; 2020; 110, pp. 1511-1521. [DOI: https://dx.doi.org/10.1094/PHYTO-02-20-0035-R] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32370659]
42. Sun, X.; Liu, D.; Zhang, X.; Li, W.; Liu, H.; Hong, W.; Jiang, C.; Guan, N.; Ma, C.; Zeng, H. et al. SLAF-seq: An efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS ONE; 2013; 8, e58700. [DOI: https://dx.doi.org/10.1371/journal.pone.0058700] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23527008]
43. Davey, J.W.; Hohenlohe, P.A.; Etter, P.D.; Boone, J.Q.; Catchen, J.M.; Blaxter, M.L. Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat. Rev. Genet.; 2011; 12, pp. 499-510. [DOI: https://dx.doi.org/10.1038/nrg3012] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21681211]
44. Li, Y.; Reif, J.C.; Ma, Y.; Hong, H.; Liu, Z.; Chang, R.; Qiu, L. Targeted association mapping demonstrating the complex molecular genetics of fatty acid formation in soybean. BMC Genom.; 2015; 16, 841. [DOI: https://dx.doi.org/10.1186/s12864-015-2049-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26494482]
45. Wang, Q.; Du, W.; Yu, W.; Zhang, W.; Huang, F.; Cheng, H.; Yu, D. Genome-wide association analysis discovered new loci and candidate genes associated with low-phosphorus tolerance based on shoot mineral elements concentrations in soybean. Mol. Genet. Genom.; 2022; 297, pp. 843-858. [DOI: https://dx.doi.org/10.1007/s00438-022-01895-1]
46. Song, Q.; Jenkins, J.; Jia, G.; Hyten, D.L.; Pantalone, V.; Jackson, S.A.; Schmutz, J.; Cregan, P.B. Construction of high resolution genetic linkage maps to improve the soybean genome sequence assembly Glyma1. 01. BMC Genom.; 2016; 17, 33. [DOI: https://dx.doi.org/10.1186/s12864-015-2344-0]
47. Rebetzke, G.J.; Burton, J.W.; Carter Jr, T.E.; Wilson, R.F. Changes in agronomic and seed characteristics with selection for reduced palmitic acid content in soybean. Crop Sci.; 1998; 38, pp. 297-302. [DOI: https://dx.doi.org/10.2135/cropsci1998.0011183X003800020003x]
48. Reinprecht, Y.; Poysa, V.W.; Yu, K.; Rajcan, I.; Ablett, G.R.; Pauls, K.P. Seed and agronomic QTL in low linolenic acid, lipoxygenase-free soybean (Glycine max (L.) Merrill) germplasm. Genome; 2006; 49, pp. 1510-1527. [DOI: https://dx.doi.org/10.1139/g06-112]
49. Hyten, D.L.; Pantalone, V.R.; Saxton, A.M.; Schmidt, M.E.; Sams, C.E. Molecular mapping and identification of soybean fatty acid modifier quantitative trait loci. J. Am. Oil Chem. Soc.; 2004; 81, pp. 1115-1118. [DOI: https://dx.doi.org/10.1007/s11746-004-1027-z]
50. Li, Z.; Wilson, R.F.; Rayford, W.E.; Boerma, H.R. Molecular mapping genes conditioning reduced palmitic acid content in N87-2122-4 soybean. Crop Sci.; 2002; 42, pp. 373-378.
51. Leamy, L.J.; Zhang, H.; Li, C.; Chen, C.Y.; Song, B.-H. A genome-wide association study of seed composition traits in wild soybean (Glycine soja). BMC Genom.; 2017; 18, 18. [DOI: https://dx.doi.org/10.1186/s12864-016-3397-4]
52. Priolli, R.; Carvalho, C.; Bajay, M.; Pinheiro, J.B.; Vello, N.A. Genome analysis to identify SNPs associated with oil content and fatty acid components in soybean. Euphytica; 2019; 215, 54. [DOI: https://dx.doi.org/10.1007/s10681-019-2378-5]
53. Priolli, R.H.G.; Campos, J.; Stabellini, N.; Pinheiro, J.B.; Vello, N.A. Association mapping of oil content and fatty acid components in soybean. Euphytica; 2015; 203, pp. 83-96. [DOI: https://dx.doi.org/10.1007/s10681-014-1264-4]
54. Huai, Y.Y.; Zhang, S.R.; Wu, T.T.; Muhammad, A.; Li, J.; Sun, S.; Han, T.F.; Li Bin Sun, J.M. Potential evaluation of molecular markers related to major nutritional quality traits in soybean breeding. Acta Agr. Sin.; 2022; 48, pp. 1957-1976.
55. Bachleda, N.; Pham, A.; Li, Z. Identifying FATB1a deletion that causes reduced palmitic acid content in soybean N87-2122-4 to develop a functional marker for marker-assisted selection. Mol. Breed.; 2016; 36, 45. [DOI: https://dx.doi.org/10.1007/s11032-016-0468-9]
56. Wang, X.; Jiang, G.-L.; Green, M.; Scott, R.A.; Hyten, D.L.; Cregan, P.B. Quantitative trait locus analysis of saturated fatty acids in a population of recombinant inbred lines of soybean. Mol. Breed.; 2012; 30, pp. 1163-1179. [DOI: https://dx.doi.org/10.1007/s11032-012-9704-0]
57. Wilson, R.F. The role of genomics and biotechnology in achieving global food security for high-oleic vegetable oil. J. Oleo Sci.; 2012; 61, pp. 357-367. [DOI: https://dx.doi.org/10.5650/jos.61.357]
58. Cahoon, E.B.; Clemente, T.E.; Damude, H.G.; Kinney, A.J. Modifying vegetable oils for food and non-food purposes. Oil Crops; Vollmann, J.; Rajcan, I. Springer: New York, NY, USA, 2009; pp. 31-56.
59. Thapa, R.; Carrero-Colón, M.; Hudson, K.A. New alleles of FATB1A to reduce palmitic acid levels in soybean. Crop Sci.; 2016; 56, pp. 1076-1080. [DOI: https://dx.doi.org/10.2135/cropsci2015.09.0597]
60. Gillman, J.D.; Stacey, M.G.; Cui, Y.; Berg, H.R.; Stacey, G. Deletions of the SACPD-C locus elevate seed stearic acid levels but also result in fatty acid and morphological alterations in nitrogen fixing nodules. BMC Plant Biol.; 2014; 14, 143. [DOI: https://dx.doi.org/10.1186/1471-2229-14-143] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24886084]
61. Bueno, R.D.; God, P.I.V.G.; Prata, I.O.; Pereira, P.H.S.; Teixeira, A.I.; Piovesan, N.D.; De Barros, E.G. Association of candidate genes for fatty acid content in soybean by temperature-switch PCR (TSP) genotyping. Crop Breed. Appl. Biotechnol.; 2018; 18, pp. 244-251. [DOI: https://dx.doi.org/10.1590/1984-70332018v18n3a37]
62. Saghaimaroof, M.A.; Soliman, K.M.; Jorgensen, R.A.; Allard, R.W. Ribosomal DNA spacer-length polymorphisms in barley: Mendelian inheritance, chromosomal location, and population dynamics. Proc. Natl. Acad. Sci. USA; 1984; 81, pp. 8014-8018. [DOI: https://dx.doi.org/10.1073/pnas.81.24.8014] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/6096873]
63. Li, H.; Ye, G.; Wang, J. A modified algorithm for the improvement of composite interval mapping. Genetics; 2007; 175, pp. 361-374. [DOI: https://dx.doi.org/10.1534/genetics.106.066811]
64. Li, H.; Ribaut, J.-M.; Li, Z.; Wang, J. Inclusive composite interval mapping (ICIM) for digenic epistasis of quantitative traits in biparental populations. Theor. Appl. Genet.; 2008; 116, pp. 243-260. [DOI: https://dx.doi.org/10.1007/s00122-007-0663-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17985112]
65. Broman, K.W.; Wu, H.; Sen, Ś.; Churchill, G.A. R/qtl: QTL mapping in experimental crosses. Bioinformatics; 2003; 19, pp. 889-890. [DOI: https://dx.doi.org/10.1093/bioinformatics/btg112]
66. Fehr, W.; Caviness, C. Stages of soybean development. Iowa Agricultural and Home Economics Experiment Station Special Report; Iowa State University: Ames, IA, USA, 1977; pp. 3-11.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Soybeans are essential crops that supply protein and oil. The composition and contents of soybean fatty acids are relevant to human health and have a significant relationship with soybean oil processing and applications. Identifying quantitative trait locus (QTL) genes related to palmitic acid could facilitate the development of a range of nutritive soybean cultivars using molecular marker-assisted selection. In this study, we used a cultivar with higher palmitic acid content, ‘Dongnong42’, and a lower palmitic acid content cultivar, ‘Hobbit’, to establish F2:6 recombinant inbred lines. A high-density genetic map containing 9980 SLAF markers was constructed and distributed across 20 soybean chromosomes. The genetic map contained a total genetic distance of 2602.58 cM and an average genetic distance of 0.39 cM between adjacent markers. Two QTLs related to palmitic acid content were mapped using inclusive composite interval mapping, explaining 4.2–10.1% of the phenotypic variance in three different years and environments, including the QTL included in seed palmitic 7-3, which was validated by developing SSR markers. Based on the SNP/Indel and significant differential expression analyses of Dongnong42 and Hobbit, two genes, Glyma.15g119700 and Glyma.15g119800, were selected as candidate genes. The high-density genetic map, QTLs, and molecular markers will be helpful for the map-based cloning of palmitic acid content genes. These could be used to accelerate breeding for high nutritive value cultivars via molecular marker-assisted breeding.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
1 Institute of Soybean Research, Heilongjiang Provincial Academy of Agricultural Sciences, Harbin 150086, China; Key Laboratory of Soybean Biology of Ministry of Education China, Northeast Agricultural University, Harbin 150030, China
2 National Key Facility for Crop Gene Resources and Genetic Improvemen, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China
3 Institute of Soybean Research, Heilongjiang Provincial Academy of Agricultural Sciences, Harbin 150086, China
4 Key Laboratory of Soybean Biology of Ministry of Education China, Northeast Agricultural University, Harbin 150030, China
5 Key Laboratory of Soybean Biology of Ministry of Education China, Northeast Agricultural University, Harbin 150030, China; National Key Facility for Crop Gene Resources and Genetic Improvemen, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China