High-density genetic map construction and QTL mapping of first flower node in pepper (Capsicum annuum L.)

High-density genetic map construction and QTL mapping of a zigzag-shaped stem trait in tea plant (Camellia sinensis)

The highly unique zigzag-shaped stem phenotype in tea plants boasts significant ornamental value and is exceptionally rare. To investigate the genetic mechanism behind this trait, we developed BC1 artificial hybrid populations. Our genetic analysis revealed the zigzag-shaped trait as a qualitative trait. Utilizing whole-genome resequencing, we constructed a high-density genetic map from the BC1 population, incorporating 5,250 SNP markers across 15 linkage groups, covering 3,328.51 cM with an average marker interval distance of 0.68 cM. A quantitative trait locus (QTL) for the zigzag-shaped trait was identified on chromosome 4, within a 61.2 to 97.2 Mb range, accounting for a phenotypic variation explained (PVE) value of 13.62%. Within this QTL, six candidate genes were pinpointed. To better understand their roles, we analyzed gene expression in various tissues and individuals with erect and zigzag-shaped stems. The results implicated CsXTH (CSS0035625) and CsCIPK14 (CSS0044366) as potential key contributors to the zigzag-shaped stem formation. These discoveries lay a robust foundation for future functional genetic mapping and tea plant genetic enhancement.

A SLAF-based high-density genetic map construction and genetic architecture of thermotolerant traits in maize (Zea mays L.)

The leaf scorching trait at flowering is a crucial thermosensitive phenotype in maize under high temperature stress (HS), yet the genetic basis of this trait remains poorly understood. In this study, we genotyped a 254 RIL-F2:8 population, derived from the leaf scorch-free parental inbred line Abe2 and the leaf scorching maternal inbred line B73, using the specific-locus amplified fragment sequencing (SLAF-seq) method. A total of 10,112 polymorphic SLAF markers were developed, and a high-density genetic map with a total length of 1,475.88 cM was constructed. The average sequencing depth of the parents was 55.23X, and that of the progeny was 12.53X. Then, we identified a total of 16 QTLs associated with thermotolerant traits at flowering, of which four QTLs of leaf scorching damage (LS) were distributed on chromosomes 1 (qLS1), 2 (qLS2.1, qLS2.2) and 3 (qLS3), which could explain 19.73% of phenotypic variation. Combining one qLS1 locus with QTL-seq results led to the identification of 6 candidate genes. Expression experiments and sequence variation indicated that Zm00001d033328, encoding N-acetyl-gamma-glutamyl-phosphate reductase, was the most likely candidate gene controlling thermotolerant traits at flowering. In summary, the high-density genetic map and genetic basis of thermotolerant traits lay a critical foundation for mapping other complex traits and identifying the genes associated with thermotolerant traits in maize.

QTL mapping: insights into genomic regions governing component traits of yield under combined heat and drought stress in wheat

Drought and heat frequently co-occur during crop growth leading to devastating yield loss. The knowledge of the genetic loci governing component traits of yield under combined drought and heat stress is essential for enhancing the climate resilience. The present study employed a mapping population of 180 recombinant inbred lines (RILs) derived from a cross between GW322 and KAUZ to identify quantitative trait loci (QTLs) governing the component traits of yield under heat and combined stress conditions. Phenotypic evaluation was conducted across two consecutive crop seasons (2021-2022 and 2022-2023) under late sown irrigation (LSIR) and late sown restricted irrigation (LSRI) conditions at the Indian Council of Agricultural Research Institute-Indian Agricultural Research Institute (ICAR-IARI), New Delhi. Various physiological and agronomic traits of importance were measured. Genotyping was carried out with 35K SNP Axiom breeder's genotyping array. The linkage map spanned a length of 6769.45 cM, ranging from 2.28 cM/marker in 1A to 14.21 cM/marker in 5D. A total of 35 QTLs were identified across 14 chromosomes with 6B containing the highest (seven) number of QTLs. Out of 35 QTLs, 16 were major QTLs explaining the phenotypic variance greater than 10%. The study identified eight stable QTLs along with two hotspots on chromosomes 6B and 5B. Five QTLs associated with traits thousand-grain weight (TGW), normalized difference vegetation index (NDVI), and plant height (PH) were successfully validated. Candidate genes encoding antioxidant enzymes, transcription factors, and growth-related proteins were identified in the QTL regions. In silico expression analysis highlighted higher expression of transcripts TraesCS2D02G021000.1, TraesCS2D02G031000, TraesCS6A02G247900, and TraesCS6B02G421700 under stress conditions. These findings contribute to a deeper understanding of the genetic architecture underlying combined heat and drought tolerance in wheat, providing valuable insights for wheat improvement strategies under changing climatic conditions.

A point mutation in MC06g1112 encoding FLOWERING LOCUS T decreases the first flower node in bitter gourd (Momordica charantia L.)

In Cucurbitaceae crops, the first flower node (FFN) is an important agronomic trait which can impact the onset of maturity, the production of female flowers, and yield. However, the gene responsible for regulating FFN in bitter gourd is unknown. Here, we used a gynoecious line (S156G) with low FFN as the female parent and a monoecious line (K8-201) with high FFN as the male parent to obtain F1 and F2 generations. Genetic analysis indicated that the low FFN trait was incompletely dominant over the high FFN trait. A major quantitative trait locus (QTL)-Mcffn and four minor effect QTLs-Mcffn1.1, Mcffn1.2, Mcffn1.3, and Mcffn1.4 were detected by whole-genome re-sequencing-based QTL mapping in the S156G×K8-201 F2 population (n=234) cultivated in autumn 2019. The Mcffn locus was further supported by molecular marker-based QTL mapping in three S156G×K8-201 F2 populations planted in autumn 2019 (n=234), autumn 2020 (n=192), and spring 2022 (n=205). Then, the Mcffn locus was fine-mapped into a 77.98-kb physical region on pseudochromosome MC06 using a large S156G×K8-201 F2 population (n=2,402). MC06g1112, which is a homolog of FLOWERING LOCUS T (FT), was considered as the most likely Mcffn candidate gene according to both expression and sequence variation analyses between parental lines. A point mutation (C277T) in MC06g1112, which results in a P93S amino acid mutation between parental lines, may be responsible for decreasing FFN in bitter gourd. Our findings provide a helpful resource for the molecular marker-assisted selective breeding of bitter gourd.

Construction of the first high-density genetic linkage map and QTL mapping of flavonoid and leaf-size related traits in Epimedium

BACKGROUND

Leaves are the main medicinal organ in Epimedium herbs, and leaf flavonoid content is an important criterion of Epimedium herbs. However, the underlying genes that regulate leaf size and flavonoid content are unclear, which limits the use of breeding for Epimedium development. This study focuses on QTL mapping of flavonoid and leaf-size related traits in Epimedium.

RESULTS

We constructed the first high-density genetic map (HDGM) using 109 F1 hybrids of Epimedium leptorrhizum and Epimedium sagittatum over three years (2019-2021). Using 5,271 single nucleotide polymorphism (SNP) markers, an HDGM with an overall distance of 2,366.07 cM and a mean gap of 0.612 cM was generated by utilizing genotyping by sequencing (GBS) technology. Every year for three years, 46 stable quantitative trait loci (QTLs) for leaf size and flavonoid contents were discovered, including 31 stable loci for Epimedin C (EC), one stable locus for total flavone content (TFC), 12 stable loci for leaf length (LL), and two stable loci for leaf area (LA). For flavonoid content and leaf size, the phenotypic variance explained for these loci varied between 4.00 and 16.80% and 14.95 and 17.34%, respectively.

CONCLUSIONS

Forty-six stable QTLs for leaf size and flavonoid content traits were repeatedly detected over three years. The HDGM and stable QTLs are laying the basis for breeding and gene investigation in Epimedium and will contribute to accelerating the identification of desirable genotypes for Epimedium breeding.

QTL mapping and BSA-seq map a major QTL for the node of the first fruiting branch in cotton

Understanding the genetic basis of the node of the first fruiting branch (NFFB) improves early-maturity cotton breeding. Here we report QTL mapping on 200 F₂ plants and derivative F_2:3 and F_2:4 populations by genotyping by sequencing (GBS). BC₁F₂ population was constructed by backcrossing one F_2:4 line with the maternal parent JF914 and used for BSA-seq for further QTL mapping. A total of 1,305,642 SNPs were developed between the parents by GBS, and 2,907,790 SNPs were detected by BSA-seq. A high-density genetic map was constructed containing 11,488 SNPs and spanning 4,202.12 cM in length. A total of 13 QTL were mapped in the 3 tested populations. JF914 conferred favorable alleles for 11 QTL, and JF173 conferred favorable alleles for the other 2 QTL. Two stable QTL were repeatedly mapped in F_2:3 and F_2:4, including qNFFB-D3-1 and qNFFB-D6-1. Only qNFFB-D3-1 contributed more than 10% of the phenotypic variation. This QTL covered about 24.7 Mb (17,130,008-41,839,226 bp) on chromosome D3. Two regions on D3 (41,779,195-41,836,120 bp, 41,836,768-41,872,287 bp) were found by BSA-seq and covered about 92.4 Kb. This 92.4 Kb region overlapped with the stable QTL qNFFB-D3-1 and contained 8 annotated genes. By qRT-PCR, Ghir_D03G012430 showed a lower expression level from the 1- to 2-leaf stage and a higher expression level from the 3- to 6-leaf stage in the buds of JF173 than that of JF914. Ghir_D03G012390 reached the highest level at the 3- and 5-leaf stages in the buds of JF173 and JF914, respectively. As JF173 has lower NFFB and more early maturity than JF914, these two genes might be important in cell division and differentiation during NFFB formation in the seedling stage. The results of this study will facilitate a better understanding of the genetic basis of NFFB and benefit cotton molecular breeding for improving earliness traits.

Fine mapping of the reduced height gene Rht22 in tetraploid wheat landrace Jianyangailanmai (Triticum turgidum L.)

Rht22 was fine mapped in the interval of 0.53-1.48 Mb on 7AS, which reduces cell number of internode to cause semi-dwarfism in Jianyangailanmai. As a valuable germplasm resource for wheat genetic improvement, tetraploid wheat has several reduced height (Rht) and enhanced harvest index genes. Rht22, discovered in Jianyangailanmai (JAM, Triticum turgidum L., 2n = 4x = 28, AABB), significantly increases the spikelet number per spike, but its accurate chromosomal position is still unknown. In this study, a high-density genetic map was constructed using specific-length amplified fragment sequencing in an F₇ RIL_DJ population, which was derived from a cross between dwarf Polish wheat (T. polonicum L., 2n = 4x = 28, AABB) and JAM. Two plant height loci, Qph.sicau-4B and Qph.sicau-7A, were mapped on chromosomes 4BS and 7AS, respectively. Qph.sicau-7A was mapped to the 0.33-4.46 Mb interval on 7AS and likely represents the candidate region of Rht22. Fine mapping confirmed and narrowed Rht22 on chromosome arm 7AS between Xbag295.s53 and Xb295.191 in three different populations. The physical region ranged from 0.53 to 1.48 Mb and included 18 candidate genes. Transcriptome analysis of two pairs of near-isogenic lines revealed that 135 differentially expressed genes (DEGs) were associated with semi-dwarfism. Of these, the expression of 83 annotated DEGs involved in hormones synthesis and signal transduction, cell wall composition, DNA replication, microtubule and phragmoplast arrays was significantly down-regulated in the semi-dwarf line. Therefore, Rht22 causes semi-dwarfism in JAM by disrupting these cellular processes, which impairs cell proliferation and reduces internode cell number.

Identification of Fruit Traits Related QTLs and a Candidate Gene, CaBRX, Controlling Locule Number in Pepper (Capsicum annuum L.)

Fruit traits are important in pepper (Capsicum annuum L.) and affect its quality and yield. These traits are controlled by quantitative trait loci (QTLs). In this study, we identified many major QTLs that control fruit length (Ftl), fruit diameter (Ftd), fruit shape (Fts), fruit weight (Ftw) and locule number (Lcn) in the F2 and F2:3 populations developed from the QTL mapping of GS6 (P1) and Qiemen (P2). A total of 111 simple sequence repeats and insertion/deletion markers were utilized to construct a linkage map with 12 linkage groups over a length of 1320.72 cM. An inclusive composite interval mapping analysis indicated that many QTLs were detected and included ftl2.1, ftd2.1, fts1.1, ftw2.1 and lcn1.1. As a novel QTL, lcn1.1 was located between HM1112 and EPMS709, and the genetic distance was 3.18 cM covering 60 predicted genes. Within the region, we identified Capana01g004285 as a candidate gene by functional annotation and expression analysis and found that it encodes the BREVIS RADIX (BRX) protein. Knockdown of CaBRX through the virus-induced gene silencing approach in GS6 reduced the number of locules and influenced the expressions of genes related to flower and locule development, suggesting that CaBRX plays an important function in the development of locules.

The APETALA2 homolog CaFFN regulates flowering time in pepper

Flowering time is an important agronomic trait that contributes to fitness in plants. However, the genetic basis of flowering time has not been extensively studied in pepper. To understand the genetics underlying flowering time, we constructed an F₂ population by crossing a spontaneous early flowering mutant and a late-flowering pepper line. Using bulked segregant RNA-seq, a major locus controlling flowering time in this population was mapped to the end of chromosome 2. An APETALA2 (AP2) homolog (CaFFN) cosegregated with flowering time in 297 individuals of the F₂ population. A comparison between the parents revealed a naturally occurring rare SNP (SNP2T > C) that resulted in the loss of a start codon in CaFFN in the early flowering mutant. Transgenic Nicotiana benthamiana plants with high CaFFN expression exhibited a delay in flowering time and floral patterning defects. On the other hand, pepper plants with CaFFN silencing flowered early. Therefore, the CaFFN gene acts as a flowering repressor in pepper. CaFFN may function as a transcriptional activator to activate the expression of CaAGL15 and miR156e and as a transcriptional repressor to repress the expression of CaAG, CaAP1, CaSEP3, CaSOC1, and miR172b based on a qRT-PCR assay. Direct activation of CaAGL15 by CaFFN was detected using yeast one-hybrid and dual-luciferase reporter assays, consistent with the hypothesis that CaFFN regulates flowering time. Moreover, the CaFFN gene association analysis revealed a significant association with flowering time in a natural pepper population, indicating that the CaFFN gene has a broad effect on flowering time in pepper. Finally, the phylogeny, evolutionary expansion and expression patterns of CaFFN/AP2 homologs were analyzed to provide valuable insight into CaFFN. This study increases our understanding of the involvement of CaFFN in controlling flowering time in pepper, thus making CaFFN a target gene for breeding early maturing pepper.

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.

Identification of fruit size associated quantitative trait loci featuring SLAF based high-density linkage map of goji berry (Lycium spp.)

BACKGROUND

Goji (Lycium spp., 2n = 24) is a fruit bearing woody plant popular as a superfood for extensive medicinal and nutritional advantages. Fruit size associated attributes are important for evaluating small-fruited goji berry and plant architecture. The domestication traits are regulated quantitatively in crop plants but few studies have attempted on genomic regions corresponding to fruit traits.

RESULTS

In this study, we established high-resolution map using specific locus amplified fragment (SLAF) sequencing for de novo SNPs detection based on 305 F₁ individuals derived from L. chinense and L. barbarum and performed quantitative trait loci (QTL) analysis of fruit size related traits in goji berry. The genetic map contained 3495 SLAF markers on 12 LGs, spanning 1649.03 cM with 0.47 cM average interval. Female and male parents and F₁ individuals` sequencing depth was 111.85-fold and 168.72-fold and 35.80-fold, respectively. The phenotype data were collected for 2 successive years (2018-2019); however, two-year mean data were combined in an extra year (1819). Total 117 QTLs were detected corresponding to multiple traits, of which 78 QTLs in 2 individual years and 36 QTLs in extra year. Six Promising QTLs (qFW10-6.1, qFL10-2.1, qLL10-2.1, qLD10-2.1, qLD12-4.1, qLA10-2.1) were discovered influencing fruit weight, fruit length and leaf related attributes covering an interval ranged from 27.32-71.59 cM on LG10 with peak LOD of 10.48 and 14.6% PVE. Three QTLs targeting fruit sweetness (qFS3-1, qFS5-2) and fruit firmness (qFF10-1) were also identified. Strikingly, various traits QTLs were overlapped on LG10, in particular, qFL10-2.1 was co-located with qLL10-2.1, qLD10-2.1 and qLA10-2.1 among stable QTLs, harbored tightly linked markers, while qLL10-1 was one major QTL with 14.21 highest LOD and 19.3% variance. As LG10 harbored important traits QTLs, we might speculate that it could be hotspot region regulating fruit size and plant architectures.

CONCLUSIONS

This report highlights the extremely saturated linkage map using SLAF-seq and novel loci contributing fruit size-related attributes in goji berry. Our results will shed light on domestication traits and further strengthen molecular and genetic underpinnings of goji berry; moreover, these findings would better facilitate to assemble the reference genome, determining potential candidate genes and marker-assisted breeding.

A New Intra-Specific and High-Resolution Genetic Map of Eggplant Based on a RIL Population, and Location of QTLs Related to Plant Anthocyanin Pigmentation and Seed Vigour

Eggplant is the second most important solanaceous berry-producing crop after tomato. Despite mapping studies based on bi-parental progenies and GWAS approaches having been performed, an eggplant intraspecific high-resolution map is still lacking. We developed a RIL population from the intraspecific cross ‘305E40’, (androgenetic introgressed line carrying the locus Rfo-Sa1 conferring Fusarium resistance) x ‘67/3’ (breeding line whose genome sequence was recently released). One hundred and sixty-three RILs were genotyped by a genotype-by-sequencing (GBS) approach, which allowed us to identify 10,361 polymorphic sites. Overall, 267 Gb of sequencing data were generated and ~773 M Illumina paired end (PE) reads were mapped against the reference sequence. A new linkage map was developed, including 7249 SNPs assigned to the 12 chromosomes and spanning 2169.23 cM, with iaci@liberoan average distance of 0.4 cM between adjacent markers. This was used to elucidate the genetic bases of seven traits related to anthocyanin content in different organs recorded in three locations as well as seed vigor. Overall, from 7 to 17 QTLs (at least one major QTL) were identified for each trait. These results demonstrate that our newly developed map supplies valuable information for QTL fine mapping, candidate gene identification, and the development of molecular markers for marker assisted selection (MAS) of favorable alleles.