Genome-Wide Association Study of Root and Shoot Related Traits in Spring Soybean (Glycine max L.) at Seedling Stages Using SLAF-Seq

SLAF-seq technology-based genome-wide association and population structure analyses of hot pepper and sweet pepper

BACKGROUND

Utilizing Single Nucleotide Polymorphism (SNP) marker technology, a phylogenetic and agronomic trait network analysis was conducted on the collected hot pepper and sweet pepper germplasm resources, providing a theoretical basis for parental selection and new varieties.

RESULTS

Specific-locus amplified fragment sequencing (SLAF-seq) technology was employed for a genome-wide association study (GWAS) on 197 hot pepper and sweet pepper germplasm resources, generating 1404.88 Mb clean reads data with an average Q30 of 91.5% and mean GC content of 37.96%. Through sequencing data analysis, a total of 639,815 SLAF tags were obtained with an average sequencing depth of 12.16x, among which 86,381 were polymorphic SLAF tags, leading to the development of 18,145,155 SNP markers. The identified SNP markers were used for cluster analysis of the genetic structure and phylogenetic relationships of hot pepper and sweet pepper germplasm resources, dividing the 197 hot pepper and sweet pepper germplasm resources into 9 clusters. Additionally, a genome-wide association analysis was conducted on 25 agronomic traits of the 197 hot pepper and sweet pepper materials, yielding a substantial number of significantly associated SNP loci with agronomic traits. A correlation network analysis diagram was drawn among the various agronomic traits, preliminarily determining the relationships between the 25 agronomic characteristics of hot pepper and sweet pepper and the positions of 15 agronomic traits (p < 1.707 × 10-8) on the chromosomes were annotated, forming multi-trait aggregation regions.

CONCLUSIONS

Our research reveals the genetic diversity, phylogenetic relationships, and population structure of 197 hot pepper and sweet pepper germplasm resources, providing a basis for germplasm identification, resource utilization, and breeding.

Characterization of genetic diversity and identification of genetic loci associated with carbon allocation in N2 fixing soybean

BACKGROUND

Efficient capture and use of resources is critical for optimal plant growth and productivity. Both shoot and root growth are essential for resource acquisition, namely light and CO2 by the shoot and water and mineral nutrients by roots. Soybean [Glycine max (L.) Merr.], one of the most valuable crops world-wide, uses an additional strategy, symbiotic N fixation (SNF), for N acquisition. SNF relies on development of specialized root organs known as nodules, which represent a distinct C sink. The genetic diversity of C partitioning in N fixing soybean to shoots, roots, and nodules has not been previously investigated but is valuable to better understand consequences of differential C allocation and to develop genetic resources, including identification of quantitative trait loci (QTLs).

RESULTS

A diversity panel of 402 soybean genotypes was phenotyped outdoors in a deep-tube system without addition of mineral N to measure allocation of biomass to the shoot, root, and nodules, as well as to determine nodule number, mean nodule biomass, and total shoot N accumulation. Wide ranges in phenotypes were observed for each of these traits, demonstrating extensive natural diversity in C partitioning and SNF in soybean. Using a set of 35,647 single nucleotide polymorphism (SNP) markers, we identified 121 SNPs tagging 103 QTLs that include both 84 novel and 19 previously identified QTLs for the eight examined traits. A candidate gene search identified 79 promising gene models in the vicinity of these QTLs. Favorable alleles of QTLs identified here may be used in breeding programs to develop elite cultivars with altered C partitioning.

CONCLUSIONS

This study provides novel insights into the diversity of biomass allocation in soybean and illustrates that the traits measured here are heritable and quantitative. QTLs identified in this study can be used in genomic prediction models as well as for further investigation of candidate genes and their roles in determining partitioning of fixed C. Enhancing our understanding of C partitioning in plants may lead to elite cultivars with optimized resource use efficiencies.

Genome-wide analysis for root and leaf architecture traits associated with drought tolerance at the seedling stage in a highly ecologically diverse wheat population

Drought stress occurred at early growth stages in wheat affecting the following growth stages. Therefore, selecting promising drought-tolerant genotypes with highly adapted traits at the seedling stage is an important task for wheat breeders and geneticists. Few research efforts were conducted on the genetic control for drought-adaptive traits at the seedling stage in wheat. In this study, a set of 146 highly diverse spring wheat core collections representing 28 different countries was evaluated under drought stress at the seedling stage. All genotypes were exposed to drought stress for 13 days by water withholding. Leaf traits including seedling length, leaf wilting, days to wilting, leaf area, and leaf rolling were scored. Moreover, root traits such as root length, maximum width, emergence angle, tip angle, and number of roots were scored. Considerable significant genetic variation was found among all genotypes tested in these experiments. The heritability estimates ranged from 0.74 (leaf witling) to 0.99 (root tip angle). A set of nine genotypes were selected and considered drought-tolerant genotypes. Among all leaf traits, shoot length had significant correlations with all root traits under drought stress. The 146 genotypes were genotyped using the Infinium Wheat 15 K single nucleotide polymorphism (SNP) array and diversity arrays technology (DArT) marker platform. The result of genotyping revealed 12,999 SNPs and 2150 DArT markers which were used to run a genome-wide association study (GWAS). The results of GWAS revealed 169 markers associated with leaf and root traits under drought stress. Out of the 169 markers, 82 were considered major quantitative trait loci (QTL). The GWAS revealed 95 candidate genes were identified with 53 genes showing evidence for drought tolerance in wheat, while the remaining candidate genes were considered novel. No shared markers were found between leaf and root traits. The results of the study provided mapping novel markers associated with new root traits at the seedling stage. Also, the selected genotypes from different countries could be employed in future wheat breeding programs not only for improving adaptive drought-tolerant traits but also for expanding genetic diversity.

Assessment of phylogenetic relationships and genetic diversity of Sagittaria trifolia using phenotypic traits and SNP markers

The aquatic perennial herb Sagittaria trifolia L. commonly known as arrowhead, has been utilized in China both as a culinary vegetable and in traditional medicines. Characterizing the phylogenetic relationships and genetic diversity of arrowheads is crucial for improved management, conservation, and efficient utilization of the germplasm resources associated with this species. Herein, we presented the phenotypic traits and genome-wide DNA marker-based analyses of 111 arrowhead accessions, most of which were from China. Cluster analysis revealed that arrowhead could be categorized into two clusters based on 11 phenotypic traits, with Cluster 1 comprising two subclusters. All accessions were clustered into three sub-clusters based primarily on leaf shape and tuber weight. A set of 759,237 high-quality single-nucleotide polymorphisms was identified and used to assess the phylogenetic relationships. Population structure and phylogenetic tree analyses suggested that the accessions could be classified into two major groups, Group I was further subdivided into two subgroups, aligning with the clusters identified through morphological classification. By employing Sagittaria lichuanensis as an outgroup, the rooted tree revealed that the evolutionary relationships within the three groups followed a progression from Group I-1 to Group I-2 and finally to Group II. The landraces were clustered into one group along with the remaining wild accessions. The level of genetic diversity for Group I (π = 0.26) was slightly lower than that which was estimated for Group II (π = 0.29). The lowest pairwise differentiation levels (Fst, 0.008) were obtained from the comparison between groups I-2 and II, indicating that the two groups were the most closely related. This study provides novel insights into germplasm classification, evolutionary relationships, genomics and arrowhead breeding.

Pilot-scale genome-wide association mapping in diverse sorghum germplasms identified novel genetic loci linked to major agronomic, root and stomatal traits

This genome-wide association studies (GWAS) used a subset of 96 diverse sorghum accessions, constructed from a large collection of 219 accessions for mining novel genetic loci linked to major agronomic, root morphological and physiological traits. The subset yielded 43,452 high quality single nucleotide polymorphic (SNP) markers exhibiting high allelic diversity. Population stratification showed distinct separation between caudatum and durra races. Linkage disequilibrium (LD) decay was rapidly declining with increasing physical distance across all chromosomes. The initial 50% LD decay was ~ 5 Kb and background level was within ~ 80 Kb. This study detected 42 significant quantitative trait nucleotide (QTNs) for different traits evaluated using FarmCPU, SUPER and 3VmrMLM which were in proximity with candidate genes related and were co-localized in already reported quantitative trait loci (QTL) and phenotypic variance (R2) of these QTNs ranged from 3 to 20%. Haplotype validation of the candidate genes from this study resulted nine genes showing significant phenotypic difference between different haplotypes. Three novel candidate genes associated with agronomic traits were validated including Sobic.001G499000, a potassium channel tetramerization domain protein for plant height, Sobic.010G186600, a nucleoporin-related gene for dry biomass, and Sobic.002G022600 encoding AP2-like ethylene-responsive transcription factor for plant yield. Several other candidate genes were validated and associated with different root and physiological traits including Sobic.005G104100, peroxidase 13-related gene with root length, Sobic.010G043300, homologous to Traes_5BL_8D494D60C, encoding inhibitor of apoptosis with iWUE, and Sobic.010G125500, encoding zinc finger, C3HC4 type domain with Abaxial stomatal density. In this study, 3VmrMLM was more powerful than FarmCPU and SUPER for detecting QTNs and having more breeding value indicating its reliable output for validation. This study justified that the constructed subset of diverse sorghums can be used as a panel for mapping other key traits to accelerate molecular breeding in sorghum.

Genome-wide association study and high-quality gene mining related to soybean protein and fat

BACKGROUND

Soybean is one of the most important oil crops in the world, and its protein and fat are the primary sources of edible oil and vegetable protein. The effective components in soybean protein and fat have positive effects on improving human immunity, anti-tumor, and regulating blood lipids and metabolism. Therefore, increasing the contents of protein and fat in soybeans is essential for improving the quality of soybeans.

RESULTS

This study selected 292 soybean lines from different regions as experimental materials, based on SLAF-seq sequencing technology, and performed genome-wide association study (GWAS) on the phenotype data from 2019-2021 Planted at the experimental base of Jilin Agricultural University, such as the contents of protein and fat of soybeans. Through the GLM model and MLM model, four SNP sites (Gm09_39012959, Gm12_35492373, Gm16_9297124, and Gm20_24678362) that were significantly related to soybean fat content were associated for three consecutive years, and two SNP sites (Gm09_39012959 and Gm20_24678362) that were significantly related to soybean protein content were associated. By the annotation and enrichment of genes within the 100 Kb region of SNP loci flanking, two genes (Glyma.09G158100 and Glyma.09G158200) related to soybean protein synthesis and one gene (Glyma.12G180200) related to lipid metabolism were selected. By the preliminary verification of expression levels of genes with qPCR, it is found that during the periods of R6 and R7 of the accumulation of soybean protein and fat, Glyma.09G158100 and Glyma.09G158200 are positive regulatory genes that promote protein synthesis and accumulation, while Glyma.12G180200 is the negative regulatory gene that inhibits fat accumulation.

CONCLUSIONS

These results lay the basis for further verifying the gene function and studying the molecular mechanisms regulating the accumulation of protein and fat in soybean seeds.

Ideotype breeding and genome engineering for legume crop improvement

Ideotype breeding is a strategy whereby traits are modeled a priori and then introduced into a model or crop species to assess their impact on yield. Thus, knowledge about the connection between genotype and phenotype is required for ideotype breeding to be deployed successfully. The growing understanding of the genetic basis of yield-related traits, combined with increasingly efficient genome engineering tools, improved transformation efficiency, and high-throughput genotyping of regenerants paves the way for the widespread adoption of ideotype breeding as a complement to conventional breeding. We briefly discuss how ideotype breeding, coupled with such state-of-the-art biotechnological tools, could contribute to knowledge-based legume breeding and accelerate yield gains to ensure food security in the coming decades.

Genetic diversity, population structure, and genome-wide association study for the flowering trait in a diverse panel of 428 moth bean (Vigna aconitifolia) accessions using genotyping by sequencing

BACKGROUND

Moth bean (Vigna aconitifolia) is an underutilized, protein-rich legume that is grown in arid and semi-arid areas of south Asia and is highly resistant to abiotic stresses such as heat and drought. Despite its economic importance, the crop remains unexplored at the genomic level for genetic diversity and trait mapping studies. To date, there is no report of SNP marker discovery and association mapping of any trait in this crop. Therefore, this study aimed to dissect the genetic diversity, population structure and marker-trait association for the flowering trait in a diversity panel of 428 moth bean accessions using genotyping by sequencing (GBS) approach.

RESULTS

A total of 9078 high-quality single nucleotide polymorphisms (SNPs) were discovered by genotyping of 428 moth bean accessions. Model-based structure analysis and PCA grouped the moth bean accessions into two subpopulations. Cluster analysis revealed accessions belonging to the Northwestern region of India had higher variability than accessions from the other regions suggesting that this region represents its center of diversity. AMOVA revealed more variations within individuals (74%) and among the individuals (24%) than among the populations (2%). Marker-trait association analysis using seven multi-locus models including mrMLM, FASTmrEMMA FASTmrEMMA, ISIS EM-BLASSO, MLMM, BLINK and FarmCPU revealed 29 potential genomic regions for the trait days to 50% flowering, which were consistently detected in three or more models. Analysis of the allelic effect of the major genomic regions explaining phenotypic variance of more than 10% and those detected in at least 2 environments showed 4 genomic regions with significant phenotypic effect on this trait. Further, we also analyzed genetic relationships among the Vigna species using SNP markers. The genomic localization of moth bean SNPs on genomes of closely related Vigna species demonstrated that maximum numbers of SNPs were getting localized on Vigna mungo. This suggested that the moth bean is most closely related to V. mungo.

CONCLUSION

Our study shows that the north-western regions of India represent the center of diversity of the moth bean. Further, the study revealed flowering-related genomic regions/candidate genes which can be potentially exploited in breeding programs to develop early-maturity moth bean varieties.

Genome-wide association study reveals genetic loci and candidate genes for meat quality traits in a four-way crossbred pig population

Meat quality traits (MQTs) have gained more attention from breeders due to their increasing economic value in the commercial pig industry. In this genome-wide association study (GWAS), 223 four-way intercross pigs were genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) and phenotyped for PH at 45 min post mortem (PH45), meat color score (MC), marbling score (MA), water loss rate (WL), drip loss (DL) in the longissimus muscle, and cooking loss (CL) in the psoas major muscle. A total of 227, 921 filtered single nucleotide polymorphisms (SNPs) evenly distributed across the entire genome were detected to perform GWAS. A total of 64 SNPs were identified for six meat quality traits using the mixed linear model (MLM), of which 24 SNPs were located in previously reported QTL regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43% to 16.32%. The genomic heritability estimates based on SNP for six meat-quality traits were low to moderate (0.07-0.47) being the lowest for CL and the highest for DL. A total of 30 genes located within 10 kb upstream or downstream of these significant SNPs were found. Furthermore, several candidate genes for MQTs were detected, including pH45 (GRM8), MC (ANKRD6), MA (MACROD2 and ABCG1), WL (TMEM50A), CL (PIP4K2A) and DL (CDYL2, CHL1, ABCA4, ZAG and SLC1A2). This study provided substantial new evidence for several candidate genes to participate in different pork quality traits. The identification of these SNPs and candidate genes provided a basis for molecular marker-assisted breeding and improvement of pork quality traits.

Construction of a high-resolution genetic map and identification of single nucleotide polymorphism markers relevant to flower stalk height in onion

INTRODUCTION

Onion (Allium cepa L., 2n=16) is an economically and nutritionally important vegetable crop worldwide. Construction of a high-resolution genetic map and map-based gene mining in onion have lagged behind other vegetable crops such as tomato and pepper.

METHODS

In this study, we constructed a high-resolution genetic map of onion using 321 F2 individuals from a cross between two double haploid lines DH-1×DH-17 and employing specific length amplified fragment (SLAF)-seq technology. The genetic map containing 10,584 polymorphic SLAFs with 21,250 single nucleotide polymorphism (SNP) markers and 8 linkage groups was developed for onion, which spanned 928.32 cM, with an average distance of 0.09 cM between adjacent markers.

RESULTS

Using this map, we carried out QTL mapping of Ms locus related to the male-fertile trait and reproduced previous mapping results, which proved that this map was of good quality. Then, four QTLs (located on LG2, LG5, and LG8) were detected for flower stalk height, explaining 26.60% of the phenotypic variance. Among them, we proposed that 20 SLAF markers (in three QTLs) of flower stalk height trait were effective favorable allelic variant markers associated with heterosis.

DISCUSSION

Overall, the genetic map was structured using SLAF-seq based on DH lines, and it is the highest-quality and highest-resolution linkage map of onion to date. It lays a foundation for the fine mapping and candidate gene identification of flower stalk height, and provides new insights into the developmental genetic mechanisms in onion breeding.

Genome-Wide Association Studies of Seven Root Traits in Soybean (Glycine max L.) Landraces

Soybean [Glycine max (L.) Merr.], an important oilseed crop, is a low-cost source of protein and oil. In Southeast Asia and Africa, soybeans are widely cultivated for use as traditional food and feed and industrial purposes. Given the ongoing changes in global climate, developing crops that are resistant to climatic extremes and produce viable yields under predicted climatic conditions will be essential in the coming decades. To develop such crops, it will be necessary to gain a thorough understanding of the genetic basis of agronomic and plant root traits. As plant roots generally lie beneath the soil surface, detailed observations and phenotyping throughout plant development present several challenges, and thus the associated traits have tended to be ignored in genomics studies. In this study, we phenotyped 357 soybean landraces at the early vegetative (V2) growth stages and used a 180 K single-nucleotide polymorphism (SNP) soybean array in a genome-wide association study (GWAS) conducted to determine the phenotypic relationships among root traits, elucidate the genetic bases, and identify significant SNPs associated with root trait-controlling genomic regions/loci. A total of 112 significant SNP loci/regions were detected for seven root traits, and we identified 55 putative candidate genes considered to be the most promising. Our findings in this study indicate that a combined approach based on SNP array and GWAS analyses can be applied to unravel the genetic basis of complex root traits in soybean, and may provide an alternative high-resolution marker strategy to traditional bi-parental mapping. In addition, the identified SNPs, candidate genes, and diverse variations in the root traits of soybean landraces will serve as a valuable basis for further application in genetic studies and the breeding of climate-resilient soybeans characterized by improved root traits.

Genome-Wide Association Studies of Salt Tolerance at the Seed Germination Stage and Yield-Related Traits in Brassica napus L

Salt stress severely affects crop growth and development and reduces the yield of Brassica napus. Exploring natural genetic variations for high salt tolerance in B. napus seedlings is an effective approach to improve productivity under salt stress. Using 10,658 high-quality single nucleotide polymorphic (SNP) markers developed by specific-locus amplified fragment sequencing (SLAF-seq) technology, genome-wide association studies (GWAS) were performed to investigate the genetic basis of salt tolerance and yield-related traits of B. napus. The results revealed that 77 and 497 SNPs were significantly associated with salt tolerance and yield-related traits, of which 40 and 58 SNPs were located in previously reported QTLs/SNPs, respectively. We identified nineteen candidate genes orthologous with Arabidopsis genes known to be associated with salt tolerance and seven potential candidates controlling both salt tolerance and yield. Our study provides a novel genetic resource for the breeding of high-yield cultivars resistant to salt stress.

Deciphering the Genetic Basis of Root and Biomass Traits in Rapeseed (Brassica napus L.) through the Integration of GWAS and RNA-Seq under Nitrogen Stress

An excellent root system is responsible for crops with high nitrogen-use efficiency (NUE). The current study evaluated the natural variations in 13 root- and biomass-related traits under a low nitrogen (LN) treatment in a rapeseed association panel. The studied traits exhibited significant phenotypic differences with heritabilities ranging from 0.53 to 0.66, and most of the traits showed significant correlations with each other. The genome-wide association study (GWAS) found 51 significant and 30 suggestive trait–SNP associations that integrated into 14 valid quantitative trait loci (QTL) clusters and explained 5.7–21.2% phenotypic variance. In addition, RNA sequencing was performed at two time points to examine the differential expression of genes (DEGs) between high and low NUE lines. In total, 245, 540, and 399 DEGs were identified as LN stress-specific, high nitrogen (HN) condition-specific, and HNLN common DEGs, respectively. An integrated analysis of GWAS, weighted gene co-expression network, and DEGs revealed 16 genes involved in rapeseed root development under LN stress. Previous studies have reported that the homologs of seven out of sixteen potential genes control root growth and NUE. These findings revealed the genetic basis underlying nitrogen stress and provided worthwhile SNPs/genes information for the genetic improvement of NUE in rapeseed.

Fine Mapping and Functional Research of Key Genes for Photoperiod Sensitivity in Maize

Maize is native to the tropics and is very sensitive to photoperiod. Planting in temperate regions with increased hours of daylight always leads to late flowering, sterility, leggy plants, and increased numbers of maize leaves. This phenomenon severely affects the utilization of tropical maize germplasm resources. The sensitivity to photoperiod is mainly reflected in differences in plant height (PH), ear height (EH), total leaf number (LN), leaf number under ear (LE), silking stage (SS), and anthesis stage (AT) in the same variety under different photoperiod conditions. These differences are more pronounced for varieties that are more sensitive to photoperiod. In the current study, a high-density genetic map was constructed from a recombinant inbred line (RIL) population containing 209 lines to map the quantitative trait loci (QTL) for photoperiod sensitivity of PH, EH, LN, LE, SS, and AT. A total of 39 QTL were identified, including three consistent major QTL. We identified candidate genes in the consensus major QTL region by combined analysis of transcriptome data, and after enrichment by GO and KEGG, we identified a total of four genes (Zm00001d006212, Zm00001d017241, Zm00001d047761, and Zm00001d047632) enriched in the plant circadian rhythm pathway (KEGG:04712). We analyzed the expression levels of these four genes, and the analysis results showed that there were significant differences in response under different photoperiod conditions for three of them (Zm00001d047761, Zm00001d006212 and Zm00001d017241). The results of functional verification showed that the expression patterns of genes rhythmically oscillated, which can affect the length of the hypocotyl and the development of the shoot apical meristem. We also found that the phenotypes of the positive plants were significantly different from the control plants when they overexpressed the objective gene or when it was knocked out, and the expression period, phase, and amplitude of the target gene also shifted. The objective gene changed its own rhythmic oscillation period, phase, and amplitude with the change in the photoperiod, thereby regulating the photoperiod sensitivity of maize. These results deepen our understanding of the genetic structure of photoperiod sensitivity and lay a foundation for further exploration of the regulatory mechanism of photoperiod sensitivity.

High-Altitude Genetic Selection and Genome-Wide Association Analysis of Yield-Related Traits in Elymus sibiricus L. Using SLAF Sequencing

The genetic adaptations to harsh climatic conditions in high altitudes and genetic basis of important agronomic traits are poorly understood in Elymus sibiricus L. In this study, an association population of 210 genotypes was used for population structure, selective sweep analysis, and genome-wide association study (GWAS) based on 88,506 single nucleotide polymorphisms (SNPs). We found 965 alleles under the natural selection of high altitude, which included 7 hub genes involved in the response to UV, and flavonoid and anthocyanin biosynthetic process based on the protein-protein interaction (PPI) analysis. Using a mixed linear model (MLM), the GWAS test identified a total of 1,825 significant loci associated with 12 agronomic traits. Based on the gene expression data of two wheat cultivars and the PPI analysis, we finally identified 12 hub genes. Especially, in plant height traits, the top hub gene (TOPLESS protein) encoding auxins and jasmonic acid signaling pathway, shoot apical meristem specification, and xylem and phloem pattern formation was highly overexpressed. These genes might play essential roles in controlling the growth and development of E. sibiricus. Therefore, this study provides fundamental insights relevant to hub genes and will benefit molecular breeding and improvement in E. sibiricus and other Elymus species.