High-resolution mapping of QTL for fatty acid composition in soybean using specific-locus amplified fragment sequencing

Identification of candidate genes and genomic prediction of soybean fatty acid components in two soybean populations

Soybean, a source of plant-derived lipids, contains an array of fatty acids essential for health. A comprehensive understanding of the fatty acid profiles in soybean is crucial for enhancing soybean cultivars and augmenting their qualitative attributes. Here, 180 F10 generation recombinant inbred lines (RILs), derived from the cross-breeding of the cultivated soybean variety 'Jidou 12' and the wild soybean 'Y9,' were used as primary experimental subjects. Using inclusive composite interval mapping (ICIM), this study undertook a quantitative trait locus (QTL) analysis on five distinct fatty acid components in the RIL population from 2019 to 2021. Concurrently, a genome-wide association study (GWAS) was conducted on 290 samples from a genetically diverse natural population to scrutinize the five fatty acid components during the same timeframe, thereby aiming to identify loci closely associated with fatty acid profiles. In addition, haplotype analysis and the Kyoto Encyclopedia of Genes and Genomes pathway analysis were performed to predict candidate genes. The QTL analysis elucidated 23 stable QTLs intricately associated with the five fatty acid components, exhibiting phenotypic contribution rates ranging from 2.78% to 25.37%. In addition, GWAS of the natural population unveiled 102 significant loci associated with these fatty acid components. The haplotype analysis of the colocalized loci revealed that Glyma.06G221400 on chromosome 6 exhibited a significant correlation with stearic acid content, with Hap1 showing a markedly elevated stearic acid level compared with Hap2 and Hap3. Similarly, Glyma.12G075100 on chromosome 12 was significantly associated with the contents of oleic, linoleic, and linolenic acids, suggesting its involvement in fatty acid biosynthesis. In the natural population, candidate genes associated with the contents of palmitic and linolenic acids were predominantly from the fatty acid metabolic pathway, indicating their potential role as pivotal genes in the critical steps of fatty acid metabolism. Furthermore, genomic selection (GS) for fatty acid components was conducted using ridge regression best linear unbiased prediction based on both random single nucleotide polymorphisms (SNPs) and SNPs significantly associated with fatty acid components identified by GWAS. GS accuracy was contingent upon the SNP set used. Notably, GS efficiency was enhanced when using SNPs derived from QTL mapping analysis and GWAS compared with random SNPs, and reached a plateau when the number of SNP markers exceeded 3,000. This study thus indicates that Glyma.06G221400 and Glyma.12G075100 are genes integral to the synthesis and regulatory mechanisms of fatty acids. It provides insights into the complex biosynthesis and regulation of fatty acids, with significant implications for the directed improvement of soybean oil quality and the selection of superior soybean varieties. The SNP markers delineated in this study can be instrumental in establishing an efficacious pipeline for marker-assisted selection and GS aimed at improving soybean fatty acid components.

The challenge of breeding for reduced off-flavor in faba bean ingredients

The growing interest in plant protein sources, such as pulses, is driven by the necessity for sustainable food production and climate change mitigation strategies. Faba bean (Vicia faba L.) is a promising protein crop for temperate climates, owing to its remarkable yield potential (up to 8 tonnes ha-1 in favourable growing conditions) and high protein content (~29% dry matter basis). Nevertheless, the adoption of faba bean protein in plant-based products that aim to resemble animal-derived counterparts is hindered by its distinctive taste and aroma, regarded as "off-flavors". In this review, we propose to introduce off-flavor as a trait in breeding programs by identifying molecules involved in sensory perception and defining key breeding targets. We discuss the role of lipid oxidation in producing volatile and non-volatile compounds responsible for the beany aroma and bitter taste, respectively. We further investigate the contribution of saponin, tannin, and other polyphenols to bitterness and astringency. To develop faba bean varieties with diminished off-flavors, we suggest targeting genes to reduce lipid oxidation, such as lipoxygenases (lox) and fatty acid desaturases (fad), and genes involved in phenylpropanoid and saponin biosynthesis, such as zero-tannin (zt), chalcone isomerase (chi), chalcone synthase (chs), β-amyrin (bas1). Additionally, we address potential challenges, including the need for high-throughput phenotyping and possible limitations that could arise during the genetic improvement process. The breeding approach can facilitate the use of faba bean protein in plant-based food such as meat and dairy analogues more extensively, fostering a transition toward more sustainable and climate-resilient diets.

Natural variation of GmFATA1B regulates seed oil content and composition in soybean

Soybean (Glycine max) produces seeds that are rich in unsaturated fatty acids and is an important oilseed crop worldwide. Seed oil content and composition largely determine the economic value of soybean. Due to natural genetic variation, seed oil content varies substantially across soybean cultivars. Although much progress has been made in elucidating the genetic trajectory underlying fatty acid metabolism and oil biosynthesis in plants, the causal genes for many quantitative trait loci (QTLs) regulating seed oil content in soybean remain to be revealed. In this study, we identified GmFATA1B as the gene underlying a QTL that regulates seed oil content and composition, as well as seed size in soybean. Nine extra amino acids in the conserved region of GmFATA1B impair its function as a fatty acyl-acyl carrier protein thioesterase, thereby affecting seed oil content and composition. Heterogeneously overexpressing the functional GmFATA1B allele in Arabidopsis thaliana increased both the total oil content and the oleic acid and linoleic acid contents of seeds. Our findings uncover a previously unknown locus underlying variation in seed oil content in soybean and lay the foundation for improving seed oil content and composition in soybean.

Discovery of common loci and candidate genes for controlling salt-alkali tolerance and yield-related traits in Brassica napus L

KEY MESSAGE

Common loci and candidate genes for controlling salt-alkali tolerance and yield-related traits were identified in Brassica napus combining QTL mapping with transcriptome under salt and alkaline stresses. The yield of rapeseed (Brassica napus L.) is determined by multiple yield-related traits, which are susceptible to environmental factors. Many yield-related quantitative trait loci (QTLs) have been reported in Brassica napus; however, no studies have been conducted to investigate both salt-alkali tolerance and yield-related traits simultaneously. Here, specific-locus amplified fragment sequencing (SLAF-seq) technologies were utilized to map the QTLs for salt-alkali tolerance and yield-related traits. A total of 65 QTLs were identified, including 30 QTLs for salt-alkali tolerance traits and 35 QTLs for yield-related traits, accounting for 7.61-27.84% of the total phenotypic variations. Among these QTLs, 18 unique QTLs controlling two to four traits were identified by meta-analysis. Six novel and unique QTLs were detected for salt-alkali tolerance traits. By comparing these unique QTLs for salt-alkali tolerance traits with those previously reported QTLs for yield-related traits, seven co-localized chromosomal regions were identified on A09 and A10. Combining QTL mapping with transcriptome of two parents under salt and alkaline stresses, thirteen genes were identified as the candidates controlling both salt-alkali tolerance and yield. These findings provide useful information for future breeding of high-yield cultivars resistant to alkaline and salt stresses.

Dual-function C2H2-type zinc-finger transcription factor GmZFP7 contributes to isoflavone accumulation in soybean

Isoflavones are a class of secondary metabolites produced by legumes and play important roles in human health and plant stress tolerance. The C2H2 zinc-finger transcription factor (TF) functions in plant stress tolerance, but little is known about its function in isoflavone regulation in soybean (Glycine max). Here, we report a C2H2 zinc-finger TF gene, GmZFP7, which regulates isoflavone accumulation in soybean. Overexpressing GmZFP7 increased the isoflavone concentration in both transgenic hairy roots and plants. By contrast, silencing GmZFP7 expression significantly reduced isoflavone levels. Metabolomic and qRT-PCR analysis revealed that GmZFP7 can increase the flux of the phenylpropanoid pathway. Furthermore, dual-luciferase and electrophoretic mobility shift assays showed that GmZFP7 regulates isoflavone accumulation by influencing the expression of Isoflavone synthase 2 (GmIFS2) and Flavanone 3 β-hydroxylase 1 (GmF3H1). In this study, we demonstrate that GmZFP7 contributes to isoflavone accumulation by regulating the expression of the gateway enzymes (GmIFS2 and GmF3H1) of competing phenylpropanoid pathway branches to direct the metabolic flux into isoflavone. A haplotype analysis indicated that important natural variations were present in GmZFP7 promoters, with P-Hap1 and P-Hap3 being the elite haplotypes. Our findings provide insight into how GmZFP7 regulates the phenylpropanoid pathway and enhances soybean isoflavone content.

MOTHER-OF-FT-AND-TFL1 regulates the seed oil and protein content in soybean

Soybean is a major crop that produces valuable seed oil and protein for global consumption. Seed oil and protein are regulated by complex quantitative trait loci (QTLs) and have undergone intensive selections during the domestication of soybean. It is essential to identify the major genetic components and understand their mechanism behind seed oil and protein in soybean. We report that MOTHER-OF-FT-AND-TFL1 (GmMFT) is the gene of a classical QTL that has been reported to regulate seed oil and protein content in many studies. Mutation of MFT decreased seeds oil content and weight in both Arabidopsis and soybean, whereas increased expression of GmMFT enhanced seeds oil content and weight. Haplotype analysis showed that GmMFT has undergone selection, which resulted in the extended haplotype homozygosity in the cultivated soybean and the enriching of the oil-favorable allele in modern soybean cultivars. This work unraveled the GmMFT-mediated mechanism regulating seed oil and protein content and seed weight, and revealed a previously unknown function of MFT that provides new insights into targeted soybean improvement and breeding.

Balancing functional and health benefits of food products formulated with palm oil as oil sources

Palm oil (PO) is widely utilised in the food industry and consumed in large quantities by humans. Owing to its bioactive components, such as fatty acids, carotenoids, vitamin E, and phenolic compounds, PO has been utilised for generations. However, public concern about their adverse effects on human health is growing. A literature search was conducted to identify fractionated palm oil processing techniques, proof of their health advantages, and potential food applications. Refined palm oil (RPO) is made from crude palm oil (CPO) and can be fractionated into palm olein (POl) and palm stearin (PS). Fractional crystallisation, dry fractionation, and solvent fractionation are the three basic fractionation procedures used in the PO industry. The composition of triacylglycerols and fatty acids in refined and fractionated palm oil and other vegetable oils is compared to elucidate the triacylglycerols and fatty acids that may be important in product development. It is well proven that RPO, POl, and PS extends the oil's shelf life in the food business. These oils have a more significant saturated fat content and antioxidant compounds than some vegetable oils, such as olive and coconut oils, making them more stable. Palm olein and stearin are also superior shortening agents and frying mediums for baking goods and meals. Furthermore, when ingested modestly daily, palm oils, especially RPO and POl, provide health benefits such as cardioprotective, antidiabetic, anti-inflammatory, and antithrombotic effects. Opportunities exist for fractionated palm oil to become a fat substitute; however, nutrition aspects need to be considered in further developing the market.

QTL Mapping of Palmitic Acid Content Using Specific-Locus Amplified Fragment Sequencing (SLAF-Seq) Genotyping in Soybeans (Glycine max L.)

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.

Genome wide association study to detect genetic regions related to isoflavone content in a mutant soybean population derived from radiation breeding

Isoflavones are major secondary metabolites that are exclusively produced by legumes, including soybean. Soy isoflavones play important roles in human health as well as in the plant defense system. The isoflavone content is influenced by minor-effect quantitative trait loci, which interact with polygenetic and environmental factors. It has been difficult to clarify the regulation of isoflavone biosynthesis because of its complex heritability and the influence of external factors. Here, using a genotype-by-sequencing-based genome-wide association mapping study, 189 mutant soybean genotypes (the mutant diversity pool, MDP) were genotyped on the basis of 25,646 high-quality single nucleotide polymorphisms (SNPs) with minor allele frequency of >0.01 except for missing data. All the accessions were phenotyped by determining the contents of 12 isoflavones in the soybean seeds in two consecutive years (2020 and 2021). Then, quantitative trait nucleotides (QTNs) related to isoflavone contents were identified and validated using multi-locus GWAS models. A total of 112 and 46 QTNs related to isoflavone contents were detected by multiple MLM-based models in 2020 and 2021, respectively. Of these, 12 and 5 QTNs were related to more than two types of isoflavones in 2020 and 2021, respectively. Forty-four QTNs were detected within the 441-Kb physical interval surrounding Gm05:38940662. Of them, four QTNs (Gm05:38936166, Gm05:38936167, Gm05:38940662, and Gm05:38940717) were located at Glyma.05g206900 and Glyma.05g207000, which encode glutathione S-transferase THETA 1 (GmGSTT1), as determined from previous quantitative trait loci annotations and the literature. We detected substantial differences in the transcript levels of GmGSTT1 and two other core genes (IFS1 and IFS2) in the isoflavone biosynthetic pathway between the original cultivar and its mutant. The results of this study provide new information about the factors affecting isoflavone contents in soybean seeds and will be useful for breeding soybean lines with high and stable concentrations of isoflavones.

Identification of superior haplotypes in a diverse natural population for breeding desirable plant height in soybean

KEY MESSAGE

Plant height of soybean is associated with a haplotype block on chromosome 19, which classified 211 soybean accessions into five distinct groups showing significant differences for the target trait. Genetic variation is pivotal for crop improvement. Natural populations are precious genetic resources. However, efficient strategies for the targeted utilization of these resources for quantitative traits, such as plant height (PH), are scarce. Being an important agronomic trait associated with soybean yield and quality, it is imperative to unravel the genetic mechanisms underlying PH in soybean. Here, a genome-wide association study (GWAS) was performed to identify single nucleotide polymorphisms (SNPs) significantly associated with PH in a natural population of 211 cultivated soybeans, which was genotyped with NJAU 355 K Soy SNP Array and evaluated across six environments. A total of 128 SNPs distributed across 17 chromosomes were found to be significantly associated with PH across six environments and a combined environment. Three significant SNPs were consistently identified in at least three environments on Chr.02 (AX-93958260), Chr.17 (AX-94154834), and Chr.19 (AX-93897200). Genomic regions of ~ 130 kb flanking these three consistent SNPs were considered as stable QTLs, which included 169 genes. Of these, 22 genes (including Dt1) were prioritized and defined as putative candidates controlling PH. The genomic region flanking 12 most significant SNPs was in strong linkage disequilibrium (LD). These SNPs formed a single haplotype block containing five haplotypes for PH, namely Hap-A, Hap-B, Hap-C, Hap-D, and Hap-E. Deployment of such superior haplotypes in breeding programs will enable development of improved soybean varieties with desirable plant height.

Identification of Genomic Loci and Candidate Genes Related to Seed Tocopherol Content in Soybean

Soybean seeds are primary sources of natural tocopherols used by the food and pharmaceutical industries, owing to their beneficial impacts on human health. Selection for higher tocopherol contents in seeds along with other desirable traits is an important goal in soybean breeding. In order to identify the genomic loci and candidate genes controlling tocopherol content in soybean seeds, the bulked-segregant analysis technique was performed using a natural population of soybean consisting of 1525 accessions. We constructed the bulked-segregant analysis based on 98 soybean accessions that showed extreme phenotypic variation for the target trait, consisting of 49 accessions with extremely-high and 49 accessions with extremely-low tocopherol content. A total of 144 variant sites and 109 predicted genes related to tocopherol content were identified, in which a total of 83 genes were annotated by the gene ontology functions. Furthermore, 13 enriched terms (p < 0.05) were detected, with four of them found to be highly enriched: response to lipid, response to abscisic acid, transition metal ion transmembrane transporter activity, and double-stranded DNA binding. Especially, six candidate genes were detected at 41.8−41.9 Mb genomic hotspots on chromosome 5 based on ANNOtate VARiation analysis. Among the genes, only Glyma.05G243400 carried a non-synonymous mutation that encodes a “translation elongation factor EF1A or initiation factor IF2gamma family protein” was identified. The haplotype analysis confirmed that Glyma.05G243400 exhibited highly significant variations in terms of tocopherol content across multiple experimental locations, suggesting that it can be the key candidate gene regulating soybean seed tocopherols. The present findings provide novel gene resources related to seed tocopherols for further validation by genome editing, functional characterization, and genetic improvement targeting enhanced tocopherol composition in soybean molecular breeding.

Identification and Validation of New Stable QTLs for Grain Weight and Size by Multiple Mapping Models in Common Wheat

Improvement of grain weight and size is an important objective for high-yield wheat breeding. In this study, 174 recombinant inbred lines (RILs) derived from the cross between Jing 411 and Hongmangchun 21 were used to construct a high-density genetic map by specific locus amplified fragment sequencing (SLAF-seq). Three mapping methods, including inclusive composite interval mapping (ICIM), genome-wide composite interval mapping (GCIM), and a mixed linear model performed with forward-backward stepwise (NWIM), were used to identify QTLs for thousand grain weight (TGW), grain width (GW), and grain length (GL). In total, we identified 30, 15, and 18 putative QTLs for TGW, GW, and GL that explain 1.1-33.9%, 3.1%-34.2%, and 1.7%-22.8% of the phenotypic variances, respectively. Among these, 19 (63.3%) QTLs for TGW, 10 (66.7%) for GW, and 7 (38.9%) for GL were consistent with those identified by genome-wide association analysis in 192 wheat varieties. Five new stable QTLs, including 3 for TGW (Qtgw.ahau-1B.1, Qtgw.ahau-4B.1, and Qtgw.ahau-4B.2) and 2 for GL (Qgl.ahau-2A.1 and Qgl.ahau-7A.2), were detected by the three aforementioned mapping methods across environments. Subsequently, five cleaved amplified polymorphic sequence (CAPS) markers corresponding to these QTLs were developed and validated in 180 Chinese mini-core wheat accessions. In addition, 19 potential candidate genes for Qtgw.ahau-4B.2 in a 0.31-Mb physical interval were further annotated, of which TraesCS4B02G376400 and TraesCS4B02G376800 encode a plasma membrane H+-ATPase and a serine/threonine-protein kinase, respectively. These new QTLs and CAPS markers will be useful for further marker-assisted selection and map-based cloning of target genes.

Development of introgression lines in high yielding, semi-dwarf genetic backgrounds to enable improvement of modern rice varieties for tolerance to multiple abiotic stresses free from undesirable linkage drag

Occurrence of multiple abiotic stresses in a single crop season has become more frequent than before. Most of the traditional donors possessing tolerance to abiotic stresses are tall, low-yielding with poor grain quality. To facilitate efficient use of complex polygenic traits in rice molecular breeding research, we undertook development of introgression lines in background of high-yielding, semi-dwarf varieties with good grain quality. The study reports the development and evaluations of over 25,000 introgression lines in eleven elite rice genetic backgrounds for improvement of yield under multiple abiotic-stresses such as drought, flood, high/low temperature. The developed introgression lines within each genetic background are near isogenic/recombinant inbred lines to their recipient recurrent parent with 50 to 98% background recovery and additionally carry QTLs/genes for abiotic stresses. The multiple-stress tolerant pyramided breeding lines combining high yield under normal situation and good yield under moderate to severe reproductive-stage drought, semi-dwarf plant type with good grain quality traits have been developed. The introgression lines in dwarf backgrounds open new opportunity to improve other varieties without any linkage drag as well as facilitate cloning of QTLs, identification and functional characterization of candidate genes, mechanisms associated with targeted QTLs and the genetic networks underlying complex polygenic traits.

Identification of candidate genes for key fibre-related QTLs and derivation of favourable alleles in Gossypium hirsutum recombinant inbred lines with G. barbadense introgressions

Fine mapping QTLs and identifying candidate genes for cotton fibre-quality and yield traits would be beneficial to cotton breeding. Here, we constructed a high-density genetic map by specific-locus amplified fragment sequencing (SLAF-seq) to identify QTLs associated with fibre-quality and yield traits using 239 recombinant inbred lines (RILs), which was developed from LMY22 (a high-yield Gossypium hirsutumL. cultivar) × LY343 (a superior fibre-quality germplasm with G. barbadenseL. introgressions). The genetic map spanned 3426.57 cM, including 3556 SLAF-based SNPs and 199 SSR marker loci. A total of 104 QTLs, including 67 QTLs for fibre quality and 37 QTLs for yield traits, were identified with phenotypic data collected from 7 environments. Among these, 66 QTLs were co-located in 19 QTL clusters on 12 chromosomes, and 24 QTLs were detected in three or more environments and determined to be stable. We also investigated the genomic components of LY343 and their contributions to fibre-related traits by deep sequencing the whole genome of LY343, and we found that genomic components from G. hirsutum races (which entered LY343 via its G. barbadense parent) contributed more favourable alleles than those from G. barbadense. We further identified six putative candidate genes for stable QTLs, including Gh_A03G1147 (GhPEL6), Gh_D07G1598 (GhCSLC6) and Gh_D13G1921 (GhTBL5) for fibre-length QTLs and Gh_D03G0919 (GhCOBL4), Gh_D09G1659 (GhMYB4) and Gh_D09G1690 (GhMYB85) for lint-percentage QTLs. Our results provide comprehensive insight into the genetic basis of the formation of fibre-related traits and would be helpful for cloning fibre-development-related genes as well as for marker-assisted genetic improvement in cotton.

Quantitative trait loci analysis of seed oil content and composition of wild and cultivated soybean

BACKGROUND

Soybean oil is a major source of edible oil, and the domestication of wild soybean has resulted in significant changes in oil content and composition. Extensive efforts have been made to identify genetic loci that are related to soybean oil traits. The objective of this study was to identify quantitative trait loci (QTLs) related to soybean seed oil and compare the fatty acid composition between wild and cultivated soybean.

RESULTS

Using the specific-locus amplified fragment sequencing (SLAF-seq) method, a total of 181 recombinant inbred lines (RILs) derived from a cross between wild soybean ZYD00463 (Glycine soja) and cultivated soybean WDD01514 (Glycine max) were genotyped. Finally, a high-density genetic linkage map comprising 11,398 single-nucleotide polymorphism (SNP) markers on 20 linkage groups (LGs) was constructed. Twenty-four stable QTLs for seed oil content and composition were identified by model-based composite interval mapping (CIM) across multiple environments. Among these QTLs, 23 overlapped with or were adjacent to previously reported QTLs. One QTL, qPA10_1 (5.94-9.98 Mb) on Chr. Ten is a novel locus for palmitic acid. In the intervals of stable QTLs, some interesting genes involved in lipid metabolism were detected.

CONCLUSIONS

We developed 181 RILs from a cross between wild soybean ZYD00463 and cultivated soybean WDD01514 and constructed a high-density genetic map using the SLAF-seq method. We identified 24 stable QTLs for seed oil content and compositions, which includes qPA10_1 on Chr. 10, a novel locus for palmitic acid. Some interesting genes in the QTL regions were also detected. Our study will provide useful information for scientists to learn about genetic variations in lipid metabolism between wild and cultivated soybean.

Quantitative Trait Locus (QTLs) Mapping for Quality Traits of Wheat Based on High Density Genetic Map Combined With Bulked Segregant Analysis RNA-seq (BSR-Seq) Indicates That the Basic 7S Globulin Gene Is Related to Falling Number

Numerous quantitative trait loci (QTLs) have been identified for wheat quality; however, most are confined to low-density genetic maps. In this study, based on specific-locus amplified fragment sequencing (SLAF-seq), a high-density genetic map was constructed with 193 recombinant inbred lines derived from Chuanmai 42 and Chuanmai 39. In total, 30 QTLs with phenotypic variance explained (PVE) up to 47.99% were identified for falling number (FN), grain protein content (GPC), grain hardness (GH), and starch pasting properties across three environments. Five NAM genes closely adjacent to QGPC.cib-4A probably have effects on GPC. QGH.cib-5D was the only one detected for GH with high PVE of 33.31-47.99% across the three environments and was assumed to be related to the nearest pina-D1 and pinb-D1genes. Three QTLs were identified for FN in at least two environments, of which QFN.cib-3D had relatively higher PVE of 16.58-25.74%. The positive effect of QFN.cib-3D for high FN was verified in a double-haploid population derived from Chuanmai 42 × Kechengmai 4. The combination of these QTLs has a considerable effect on increasing FN. The transcript levels of Basic 7S globulin and Basic 7S globulin 2 in QFN.cib-3D were significantly different between low FN and high FN bulks, as observed through bulk segregant RNA-seq (BSR). These QTLs and candidate genes based on the high-density genetic map would be beneficial for further understanding of the genetic mechanism of quality traits and molecular breeding of wheat.

Effect of marker segregation distortion on high density linkage map construction and QTL mapping in Soybean (Glycine max L.)

Marker segregation distortion is a natural phenomenon. Severely distorted markers are usually excluded in the construction of linkage maps. We investigated the effect of marker segregation distortion on linkage map construction and quantitative trait locus (QTL) mapping. A total of 519 recombinant inbred lines of soybean from orthogonal and reciprocal crosses between LSZZH and NN493-1 were genotyped by specific length amplified fragment markers and seed linoleic acid content was measured in three environments. As a result, twenty linkage groups were constructed with 11,846 markers, including 1513 (12.77%) significantly distorted markers, on 20 chromosomes, and the map length was 2475.86 cM with an average marker-interval of 0.21 cM. The inclusion of distorted markers in the analysis was shown to not only improve the grouping of the markers from the same chromosomes, and the consistency of linkage maps with genome, but also increase genome coverage by markers. Combining genotypic data from both orthogonal and reciprocal crosses decreased the proportion of distorted markers and then improved the quality of linkage maps. Validation of the linkage maps was confirmed by the high collinearity between positions of markers in the soybean reference genome and in linkage maps and by the high consistency of 24 QTL regions in this study compared with the previously reported QTLs and lipid metabolism related genes. Additionally, linkage maps that include distorted markers could add more information to the outputs from QTL mapping. These results provide important information for linkage mapping, gene cloning and marker-assisted selection in soybean.

Deciphering the Genetic Architecture of Plant Height in Soybean Using Two RIL Populations Sharing a Common M8206 Parent

Plant height (PH) is an important agronomic trait that is closely related to soybean yield and quality. However, it is a complex quantitative trait governed by multiple genes and is influenced by environment. Unraveling the genetic mechanism involved in PH, and developing soybean cultivars with desirable PH is an imperative goal for soybean breeding. In this regard, the present study used high-density linkage maps of two related recombinant inbred line (RIL) populations viz., MT and ZM evaluated in three different environments to detect additive and epistatic effect quantitative trait loci (QTLs) as well as their interaction with environments for PH in Chinese summer planting soybean. A total of eight and 12 QTLs were detected by combining the composite interval mapping (CIM) and mixed-model based composite interval mapping (MCIM) methods in MT and ZM populations, respectively. Among these QTLs, nine QTLs viz., QPH-2, qPH-6-2_MT, QPH-6, qPH-9-1_ZM, qPH-10-1_ZM, qPH-13-1_ZM, qPH-16-1_MT, QPH-17 and QPH-19 were consistently identified in multiple environments or populations, hence were regarded as stable QTLs. Furthermore, Out of these QTLs, three QTLs viz., qPH-4-2_ZM, qPH-15-1_MT and QPH-17 were novel. In particular, QPH-17 could detect in both populations, which was also considered as a stable and major QTL in Chinese summer planting soybean. Moreover, eleven QTLs revealed significant additive effects in both populations, and out of them only six showed additive by environment interaction effects, and the environment-independent QTLs showed higher additive effects. Finally, six digenic epistatic QTLs pairs were identified and only four additive effect QTLs viz., qPH-6-2_MT, qPH-19-1_MT/QPH-19, qPH-5-1_ZM and qPH-17-1_ZM showed epistatic effects. These results indicate that environment and epistatic interaction effects have significant influence in determining genetic basis of PH in soybean. These results would not only increase our understanding of the genetic control of plant height in summer planting soybean but also provide support for implementing marker assisted selection (MAS) in developing cultivars with ideal plant height as well as gene cloning to elucidate the mechanisms of plant height.

Genome-wide association and transcriptional studies reveal novel genes for unsaturated fatty acid synthesis in a panel of soybean accessions

BACKGROUND

The nutritional value of soybean oil is largely influenced by the proportions of unsaturated fatty acids (FAs), including oleic acid (OA, 18:1), linoleic acid (LLA, 18:2), and linolenic acid (LNA, 18:3). Genome-wide association (GWAS) studies along with gene expression studies in soybean [Glycine max (L.) Merr.] were leveraged to dissect the genetics of unsaturated FAs.

RESULTS

A association panel of 194 diverse soybean accessions were phenotyped in 2013, 2014 and 2015 to identify Single Nucleotide Polymorphisms (SNPs) associated with OA, LLA, and LNA content, and determine putative candidate genes responsible for regulating unsaturated FAs composition. 149 SNPs that represented 73 genomic regions were found to be associated with the unsaturated FA contents in soybean seeds according to the results of GWAS. Twelve novel genes were predicted to be involved in unsaturated FA synthesis in soybean. The relationship between expression pattern of the candidate genes and the accumulation of unsaturated FAs revealed that multiple genes might be involved in unsaturated FAs regulation simultaneously but work in very different ways: Glyma.07G046200 and Glyma.20G245500 promote the OA accumulation in soybean seed in all the tested accessions; Glyma.13G68600 and Glyma.16G200200 promote the OA accumulation only in high OA germplasms; Glyma.07G151300 promotes OA accumulation in higher OA germplasms and suppresses that in lower OA germplasms; Glyma.16G003500 has the effect of increasing LLA accumulation in higher LA germplasms; Glyma.07G254500 suppresses the accumulation of LNA in lower OA germplasms; Glyma.14G194300 might be involved in the accumulation of LNA content in lower LNA germplasms.

CONCLUSIONS

The beneficial alleles and candidate genes identified might be valuable for improving marker-assisted breeding efficiency and exploring the molecular mechanisms underlying unsaturated fatty acid of soybean.

Construction of a High-Density Genetic Map and Analysis of Seed-Related Traits Using Specific Length Amplified Fragment Sequencing for Cucurbita maxima

Seed traits are agronomically important for Cucurbita breeding, but the genes controlling seed size, seed weight and seed number have not been mapped in Cucurbita maxima (C. maxima). In this study, 100 F₂ individual derived from two parental lines, "2013-12" and "9-6", were applied to construct a 3,376.87-cM genetic map containing 20 linkage groups (LGs) with an average genetic distance of 0.47 cM using a total of 8,406 specific length amplified fragment (SLAF) markers in C. maxima. Ten quantitative trait loci (QTLs) of seed width (SW), seed length (SL) and hundred-seed weight (HSW) were identified using the composite interval mapping (CIM) method. The QTLs affecting SW, SL and HSW explained a maximum of 38.6%, 28.9% and 17.2% of the phenotypic variation and were detected in LG6, LG6 and LG17, respectively. To validate these results, an additional 150 F₂ individuals were used for QTL mapping of SW and SL with cleaved amplified polymorphic sequence (CAPS) markers. We found that two major QTLs, SL6-1 and SW6-1, could be detected in both SLAF-seq and CAPS markers in an overlapped region. Based on gene annotation and non-synonymous single-nucleotide polymorphisms (SNPs) in the major SWand SL-associated regions, we found that two genes encoding a VQ motif and an E3 ubiquitin-protein ligase may be candidate genes influencing SL, while an F-box and leucinerich repeat (LRR) domain-containing protein is the potential regulator for SW in C. maxima. This study provides the first high-density linkage map of C. maxima using SNPs developed by SLAF-seq technology, which is a powerful tool for associated mapping of important agronomic traits, map-based gene cloning and marker-assisted selection (MAS)-based breeding in C. maxima.

Construction of High-Density Genetic Map and Mapping Quantitative Trait Loci for Growth Habit-Related Traits of Peanut (Arachis hypogaea L.)

Plant growth habit is an important and complex agronomic trait and is associated with yield, disease resistance, and mechanized harvesting in peanuts. There are at least two distinct growth habits (erect and prostrate) and several intermediate forms existing in the peanut germplasm. A recombinant inbred line population containing 188 individuals was developed from a cross of "Jihua 5" and "M130" for genetically dissecting the architecture of the growth habit. A new high-density genetic linkage map was constructed by using specific locus amplified fragment sequencing technology. The map contains 2,808 single-nucleotide polymorphism markers distributed on 20 linkage groups with a total length of 1,308.20 cM and an average inter-marker distance of 0.47 cM. The quantitative trait locus (QTL) analysis of the growth habit-related traits was conducted based on phenotyping data from seven environments. A total of 39 QTLs for growth habit-related traits was detected on 10 chromosomes explaining 4.55-27.74% of the phenotypic variance, in which 6 QTLs were for lateral branch angle, 8 QTLs were for extent radius, 7 QTLs were for the index of plant type, 11 QTLs were for main stem height, and 7 QTLs were for lateral branch length. Among these QTLs, 12 were co-localized on chromosome B05 spanning an approximately 0.17 Mb physical interval in comparison with the allotetraploid reference genome of "Tifrunner." Analysis of the co-localized genome region has shown that the putative genes are involved in light and hormones and will facilitate peanut growth habit molecular breeding and study of peanut domestication.

QTL Mapping Combined With Bulked Segregant Analysis Identify SNP Markers Linked to Leaf Shape Traits in Pisum sativum Using SLAF Sequencing

Leaf shape is an important trait that influences the utilization rate of light, and affects quality and yield of pea (Pisum sativum). In the present study, a joint method of high-density genetic mapping using specific locus amplified fragment sequencing (SLAF-seq) and bulked segregant analysis (BSA) was applied to rapidly detect loci with leaf shape traits. A total of 7,146 polymorphic SLAFs containing 12,213 SNP markers were employed to construct a high-density genetic map for pea. We conducted quantitative trait locus (QTL) mapping on an F2 population to identify QTLs associated with leaf shape traits. Moreover, SLAF-BSA was conducted on the same F2 population to identify the single nucleotide polymorphism (SNP) markers linked to leaf shape in pea. Two QTLs (qLeaf_or-1, qLeaf_or-2) were mapped on linkage group 7 (LG7) for pea leaf shape. Through alignment of SLAF markers with Cicer arietinum, Medicago truncatula, and Glycine max, the pea LGs were assigned to their corresponding homologous chromosomal groups. The comparative genetic analysis showed that pea is more closely related to M. truncatula. Based on the sequencing results of two pools with different leaf shape, 179 associated markers were obtained after association analysis. The joint analysis of SLAF-seq and BSA showed that the QTLs obtained from mapping on a high-density genetic map are convincing due to the closely associated map region with the BSA results, which provided more potential markers related to leaf shape. Thus, the identified QTLs could be used in marker-assisted selection for pea breeding in the future. Our study revealed that joint analysis of QTL mapping on a high-density genetic map and BSA-seq is a cost-effective and accurate method to reveal genetic architecture of target traits in plant species without a reference genome.

An integrated omics analysis reveals molecular mechanisms that are associated with differences in seed oil content between Glycine max and Brassica napus

BACKGROUND

Rapeseed (Brassica napus L.) and soybean (Glycine max L.) seeds are rich in both protein and oil, which are major sources of biofuels and nutrition. Although the difference in seed oil content between soybean (~ 20%) and rapeseed (~ 40%) exists, little is known about its underlying molecular mechanism.

RESULTS

An integrated omics analysis was performed in soybean, rapeseed, Arabidopsis (Arabidopsis thaliana L. Heynh), and sesame (Sesamum indicum L.), based on Arabidopsis acyl-lipid metabolism- and carbon metabolism-related genes. As a result, candidate genes and their transcription factors and microRNAs, along with phylogenetic analysis and co-expression network analysis of the PEPC gene family, were found to be largely associated with the difference between the two species. First, three soybean genes (Glyma.13G148600, Glyma.13G207900 and Glyma.12G122900) co-expressed with GmPEPC1 are specifically enriched during seed storage protein accumulation stages, while the expression of BnPEPC1 is putatively inhibited by bna-miR169, and two genes BnSTKA and BnCKII are co-expressed with BnPEPC1 and are specifically associated with plant circadian rhythm, which are related to seed oil biosynthesis. Then, in de novo fatty acid synthesis there are rapeseed-specific genes encoding subunits β-CT (BnaC05g37990D) and BCCP1 (BnaA03g06000D) of heterogeneous ACCase, which could interfere with synthesis rate, and β-CT is positively regulated by four transcription factors (BnaA01g37250D, BnaA02g26190D, BnaC01g01040D and BnaC07g21470D). In triglyceride synthesis, GmLPAAT2 is putatively inhibited by three miRNAs (gma-miR171, gma-miR1516 and gma-miR5775). Finally, in rapeseed there was evidence for the expansion of gene families, CALO, OBO and STERO, related to lipid storage, and the contraction of gene families, LOX, LAH and HSI2, related to oil degradation.

CONCLUSIONS

The molecular mechanisms associated with differences in seed oil content provide the basis for future breeding efforts to improve seed oil content.

Construction of a high-density genetic linkage map and QTL mapping of oleic acid content and three agronomic traits in sunflower (Helianthus annuus L.) using specific-locus amplified fragment sequencing (SLAF-seq)

High-density genetic linkage maps are particularly important for quantitative trait loci (QTL) mapping, genome assembly, and marker-assisted selection (MAS) in plants. In this study, a high-density genetic linkage map of sunflower (Helianthus annuus L.) was constructed using an F₂ population generated from a cross between Helianthus annuus L. '86-1' and 'L-1-OL-1' via specific-locus amplified fragment sequencing (SLAF-seq). After sequence preprocessing, 530.50 M reads (105.60 Gb) were obtained that contained a total of 343,197 SLAFs, of which 39,589 were polymorphic. Of the polymorphic SLAFs, 6,136 were organized into a linkage map consisting of 17 linkage groups (LGs) spanning 2,221.86 cM, with an average genetic distance of 0.36 cM between SLAFs. Based on this high-density genetic map, QTL analysis was performed that focused on four sunflower phenotypic traits: oleic acid content (OAC), plant height (PH), head diameter (HD), and stem diameter (SD). Subsequently, for these four traits eight QTLs were detected that will likely be useful for increasing our understanding of genetic factors underlying these traits and for use in marker-assisted selection (MAS) for future sunflower breeding.

Genome-wide discovery of DNA polymorphisms among chickpea cultivars with contrasting seed size/weight and their functional relevance

Seed size/weight is a major agronomic trait which determine crop productivity in legumes. To understand the genetic basis of seed size determination, we sought to identify DNA polymorphisms between two small (Himchana 1 and Pusa 362) and two large-seeded (JGK 3 and PG 0515) chickpea cultivars via whole genome resequencing. We identified a total of 75535 single nucleotide polymorphisms (SNPs), 6486 insertions and deletions (InDels), 1938 multi-nucleotide polymorphisms (MNPs) and 5025 complex variants between the two small and two large-seeded chickpea cultivars. Our analysis revealed 814, 244 and 72 seed-specific genes harboring DNA polymorphisms in promoter or non-synonymous and large-effect DNA polymorphisms, respectively. Gene ontology analysis revealed enrichment of cell growth and division related terms in these genes. Among them, at least 22 genes associated with quantitative trait loci, and those involved in cell growth and division and encoding transcription factors harbored promoter and/or large-effect/non-synonymous DNA polymorphisms. These also showed higher expression at late-embryogenesis and/or mid-maturation stages of seed development in the large-seeded cultivar, suggesting their role in seed size/weight determination in chickpea. Altogether, this study provided a valuable resource for large-scale genotyping applications and a few putative candidate genes that might play crucial role in governing seed size/weight in chickpea.

Construction of a high-density genetic map and mapping of QTLs for soybean (Glycine max) agronomic and seed quality traits by specific length amplified fragment sequencing

BACKGROUND

Soybean is not only an important oil crop, but also an important source of edible protein and industrial raw material. Yield-traits and quality-traits are increasingly attracting the attention of breeders. Therefore, fine mapping the QTLs associated with yield-traits and quality-traits of soybean would be helpful for soybean breeders. In the present study, a high-density linkage map was constructed to identify the QTLs for the yield-traits and quality-traits, using specific length amplified fragment sequencing (SLAF-seq).

RESULTS

SLAF-seq was performed to screen SLAF markers with 149 F8:11 individuals from a cross between a semi wild soybean, 'Huapidou', and a cultivated soybean, 'Qihuang26', which generated 400.91 M paired-end reads. In total, 53,132 polymorphic SLAF markers were obtained. The genetic linkage map was constructed by 5111 SLAF markers with segregation type of aa×bb. The final map, containing 20 linkage groups (LGs), was 2909.46 cM in length with an average distance of 0.57 cM between adjacent markers. The average coverage for each SLAF marker on the map was 81.26-fold in the male parent, 45.79-fold in the female parent, and 19.84-fold average in each F8:11 individual. According to the high-density map, 35 QTLs for plant height (PH), 100-seeds weight (SW), oil content in seeds (Oil) and protein content in seeds (Protein) were found to be distributed on 17 chromosomes, and 14 novel QTLs were identified for the first time. The physical distance of 11 QTLs was shorter than 100 Kb, suggesting a direct opportunity to find candidate genes. Furthermore, three pairs of epistatic QTLs associated with Protein involving 6 loci on 5 chromosomes were identified. Moreover, 13, 14, 7 and 9 genes, which showed tissue-specific expression patterns, might be associated with PH, SW, Oil and Protein, respectively.

CONCLUSIONS

With SLAF-sequencing, some novel QTLs and important QTLs for both yield-related and quality traits were identified based on a new, high-density linkage map. Moreover, 43 genes with tissue-specific expression patterns were regarded as potential genes in further study. Our findings might be beneficial to molecular marker-assisted breeding, and could provide detailed information for accurate QTL localization.

Genome-wide Scan for Seed Composition Provides Insights into Soybean Quality Improvement and the Impacts of Domestication and Breeding

The complex genetic architecture of quality traits has hindered efforts to modify seed nutrients in soybean. Genome-wide association studies were conducted for seed composition, including protein, oil, fatty acids, and amino acids, using 313 diverse soybean germplasm accessions genotyped with a high-density SNP array. A total of 87 chromosomal regions were identified to be associated with seed composition, explaining 8%-89% of genetic variances. The candidate genes GmSAT1, AK-HSDH, SACPD-C, and FAD3A of known function, and putative MtN21 nodulin, FATB, and steroid-5-α-reductase involved in N2 fixation, amino acid biosynthesis, and fatty acid metabolism were found at the major-effect loci. Further analysis of additional germplasm accessions indicated that these major-effect loci had been subjected to domestication or modern breeding selection, and the allelic variants and distributions were relevant to geographic regions. We also revealed that amino acid concentrations related to seed weight and to total protein had a different genetic basis. This helps uncover the in-depth genetic mechanism of the intricate relationships among the seed compounds. Thus, our study not only provides valuable genes and markers for soybean nutrient improvement, both quantitatively and qualitatively, but also offers insights into the alteration of soybean quality during domestication and breeding.