QTL-seq identifies an early flowering QTL located near Flowering Locus T in cucumber

LATA1, a RING E3 ligase, modulates the tiller angle by affecting auxin asymmetric distribution and content in rice

The tiller angle is an important agronomic trait that determines plant architecture and grain yield in rice (Oryza sativa L.). However, the molecular regulation mechanism of the rice tiller angle remains unclear. Here, we identified a rice tiller angle gene, LARGE TILLER ANGLE 1 (LATA1), using the MutMap approach. LATA1 encodes a C3H2C3-type RING zinc finger E3 ligase and the conserved region of the RING zinc finger is essential for its E3 activity. LATA1 was highly expressed in the root and tiller base and LATA1-GFP fusion protein was specifically localized to the nucleus. The mutation of LATA1 significantly reduced indole-3-acetic acid content and attenuated lateral auxin transport, thereby resulting in defective shoot gravitropism and spreading plant architecture in rice. Further investigations found that LATA1 may indirectly affect gravity perception by modulating the sedimentation rate of gravity-sensing amyloplasts upon gravistimulation. Our findings provide new insights into the molecular mechanism underlying the rice tiller angle and new genetic resource for the improvement of plant architecture in rice.

Identification of Genomic Regions Associated with Differences in Flowering Time and Inflorescence Architecture between Melastoma candidum and M. normale

Understanding the genetic basis of species differences in flowering time and inflorescence architecture can shed light on speciation and molecular breeding. Melastoma shows rapid speciation, with about 100 species formed in the past few million years, and, meanwhile, possesses high ornamental values. Two largely sympatric and closely related species of this genus, M. candidum and M. normale, differ markedly in flowering time and flower number per inflorescence. Here, we constructed an F2 population between M. candidum and M. normale, and used extreme bulks for flowering time and flower number per inflorescence in this population to identify genomic regions underlying the two traits. We found high differentiation on nearly the whole chromosome 7 plus a few regions on other chromosomes between the two extreme bulks for flowering time. Large chromosomal inversions on chromosome 7 between the two species, which contain flowering-related genes, can explain recombinational suppression on the chromosome. We identified 1872 genes with one or more highly differentiated SNPs between the two bulks for flowering time, including CSTF77, FY, SPA3, CDF3, AGL8, AGL15, FHY1, COL9, CIB1, FKF1 and FAR1, known to be related to flowering. We also identified 680 genes with one or more highly differentiated SNPs between the two bulks for flower number per inflorescence, including PNF, FIL and LAS, knows to play important roles in inflorescence development. These large inversions on chromosome 7 prevent us from narrowing down the genomic region(s) associated with flowering time differences between the two species. Flower number per inflorescence in Melastoma appears to be controlled by multiple genes, without any gene of major effect. Our study indicates that large chromosomal inversions can hamper the identification of the genetic basis of important traits, and the inflorescence architecture of Melastoma species may have a complex genetic basis.

Identification and application of a candidate gene AhAftr1 for aflatoxin production resistance in peanut seed (Arachis hypogaea L.)

INTRODUCTION

Peanut is susceptible to infection of Aspergillus fungi and conducive to aflatoxin contamination, hence developing aflatoxin-resistant variety is highly meaningful. Identifying functional genes or loci conferring aflatoxin resistance and molecular diagnostic marker are crucial for peanut breeding.

OBJECTIVES

This work aims to (1) identify candidate gene for aflatoxin production resistance, (2) reveal the related resistance mechanism, and (3) develop diagnostic marker for resistance breeding program.

METHODS

Resistance to aflatoxin production in a recombined inbred line (RIL) population derived from a high-yielding variety Xuhua13 crossed with an aflatoxin-resistant genotype Zhonghua 6 was evaluated under artificial inoculation for three consecutive years. Both genetic linkage analysis and QTL-seq were conducted for QTL mapping. The candidate gene was further fine-mapped using a secondary segregation mapping population and validated by transgenic experiments. RNA-Seq analysis among resistant and susceptible RILs was used to reveal the resistance pathway for the candidate genes.

RESULTS

The major effect QTL qAFTRA07.1 for aflatoxin production resistance was mapped to a 1.98 Mbp interval. A gene, AhAftr1 (Arachis hypogaea Aflatoxin resistance 1), was detected structure variation (SV) in leucine rich repeat (LRR) domain of its production, and involved in disease resistance response through the effector-triggered immunity (ETI) pathway. Transgenic plants with overexpression of AhAftr1(ZH6) exhibited 57.3% aflatoxin reduction compared to that of AhAftr1(XH13). A molecular diagnostic marker AFTR.Del.A07 was developed based on the SV. Thirty-six lines, with aflatoxin content decrease by over 77.67% compared to the susceptible control Zhonghua12 (ZH12), were identified from a panel of peanut germplasm accessions and breeding lines through using AFTR.Del.A07.

CONCLUSION

Our findings would provide insights of aflatoxin production resistance mechanisms and laid meaningful foundation for further breeding programs.

Delineation of loci governing an extra-earliness trait in lentil (Lens culinaris Medik.) using the QTL-Seq approach

Developing early maturing lentil has the potential to minimize yield losses, mainly during terminal drought. Whole-genome resequencing (WGRS) based QTL-seq identified the loci governing earliness in lentil. The genetic analysis for maturity duration provided a good fit to 3:1 segregation (F2), indicating earliness as a recessive trait. WGRS of Globe Mutant (late parent), late-flowering, and early-flowering bulks (from RILs) has generated 1124.57, 1052.24 million raw and clean reads, respectively. The QTL-Seq identified three QTLs (LcqDTF3.1, LcqDTF3.2, and LcqDTF3.3) on chromosome 3 having 246244 SNPs and 15577 insertions/deletions (InDels) and 13 flowering pathway genes. Of these, 11 exhibited sequence variations between bulks and validation (qPCR) revealed a significant difference in the expression of nine candidate genes (LcGA20oxG, LcFRI, LcLFY, LcSPL13a, Lcu.2RBY.3g060720, Lcu.2RBY.3g062540, Lcu.2RBY.3g062760, LcELF3a, and LcEMF1). Interestingly, the LcELF3a gene showed significantly higher expression in late-flowering genotype and exhibited substantial involvement in promoting lateness. Subsequently, an InDel marker (I-SP-383.9; LcELF3a gene) developed from LcqDTF3.2 QTL region showed 82.35% PVE (phenotypic variation explained) for earliness. The cloning, sequencing, and comparative analysis of the LcELF3a gene from both parents revealed 23 SNPs and InDels. Interestingly, a 52 bp deletion was recorded in the LcELF3a gene of L4775, predicted to cause premature termination of protein synthesis after 4 missense amino acids beyond the 351st amino acid due to the frameshift during translation. The identified InDel marker holds significant potential for breeding early maturing lentil varieties.

Chromosome-level genome assembly and population genomics reveals crucial selection for subgynoecy development in chieh-qua

Chieh-qua is an important cucurbit crop and very popular in South China and Southeast Asia. Despite its significance, its genetic basis and domestication history are unclear. In this study, we have successfully generated a chromosome-level reference genome assembly for the chieh-qua 'A36' using a hybrid assembly strategy that combines PacBio long reads and Illumina short reads. The assembled genome of chieh-qua is approximately 953.3 Mb in size and is organized into 12 chromosomes, with contig N50 of 6.9 Mb and scaffold N50 of 68.2 Mb. Notably, the chieh-qua genome is comparable in size to the wax gourd genome. Through gene prediction analysis, we have identified a total of 24 593 protein-coding genes in the A36 genome. Additionally, approximately 56.6% (539.3 Mb) of the chieh-qua genome consists of repetitive sequences. Comparative genome analysis revealed that chieh-qua and wax gourd are closely related, indicating a close evolutionary relationship between the two species. Population genomic analysis, employing 129 chieh-qua accessions and 146 wax gourd accessions, demonstrated that chieh-qua exhibits greater genetic diversity compared to wax gourd. We also employed the GWAS method to identify related QTLs associated with subgynoecy, an interested and important trait in chieh-qua. The MYB59 (BhiCQ0880026447) exhibited relatively high expression levels in the shoot apex of four subgynoecious varieties compared with monoecious varieties. Overall, this research provides insights into the domestication history of chieh-qua and offers valuable genomic resources for further molecular research.

The Dark Side of the pollen: BSA-seq identified genomic regions linked to male sterility in globe artichoke

Globe artichoke (Cynara cardunculus var. scolymus; 2n = 2x = 34) is a food crop consumed for its immature flower heads. Traditionally, globe artichoke varietal types are vegetatively propagated. However, seed propagation makes it possible to treat the crop as annual, increasing field uniformity and reducing farmers costs, as well as pathogens diffusion. Despite globe artichoke's significant agricultural value and the critical role of heterosis in the development of superior varieties, the production of hybrids remains challenging without a reliable system for large-scale industrial seed production. Male sterility (MS) presents a promising avenue for overcoming these challenges by simplifying the hybridization process and enabling cost-effective seed production. However, within the Cynara genus, genic male sterility has been linked to three recessive loci in globe artichoke, with no definitive genetic mechanism elucidated to date. A 250 offsprings F2 population, derived from a cross between a MS globe artichoke and a male fertile (MF) cultivated cardoon (C. cardunculus var. altilis) and fitting a monogenic segregation model (3:1), was analyzed through BSA-seq, aiming at the identification of genomic regions/genes affecting male sterility. Four QTL regions were identified on chromosomes 4, 12, and 14. By analyzing the sequence around the highest pick on chromosome 14, a cytochrome P450 (CYP703A2) was identified, carrying a deleterious substitution (R/Q) fixed in the male sterile parent. A single dCAPS marker was developed around this SNP, allowing the discrimination between MS and MF genotypes within the population, suitable for applications in plant breeding programs. A 3D model of the protein was generated by homology modeling, revealing that the mutated amino acid is part of a highly conserved motif crucial for protein folding.

QTLs and Genes for Salt Stress Tolerance: A Journey from Seed to Seed Continued

Rice (Oryza sativa L.) is a crucial crop contributing to global food security; however, its production is susceptible to salinity, a significant abiotic stressor that negatively impacts plant germination, vigour, and yield, degrading crop production. Due to the presence of exchangeable sodium ions (Na+), the affected plants sustain two-way damage resulting in initial osmotic stress and subsequent ion toxicity in the plants, which alters the cell's ionic homeostasis and physiological status. To adapt to salt stress, plants sense and transfer osmotic and ionic signals into their respective cells, which results in alterations of their cellular properties. No specific Na+ sensor or receptor has been identified in plants for salt stress other than the SOS pathway. Increasing productivity under salt-affected soils necessitates conventional breeding supplemented with biotechnological interventions. However, knowledge of the genetic basis of salinity stress tolerance in the breeding pool is somewhat limited because of the complicated architecture of salinity stress tolerance, which needs to be expanded to create salt-tolerant variants with better adaptability. A comprehensive study that emphasizes the QTLs, genes and governing mechanisms for salt stress tolerance is discussed in the present study for future research in crop improvement.

Screening candidate genes for fruit size based on QTL-seq in Chinese jujube

Introduction

Fruit size is an important economic trait affecting jujube fruit quality, which has always been the focus of marker-assisted breeding of jujube traits. However, despite a large number of studies have been carried out, the mechanism and key genes regulating jujube fruit size are mostly unknown.

Methods

In this study, we used a new analysis method Quantitative Trait Loci sequencing (QTL-seq) (bulked segregant analysis) to screen the parents 'Yuhong' and 'Jiaocheng 5' with significant phenotypic differences and mixed offspring group with extreme traits of large fruit and small fruit, respectively, and, then, DNA mixed pool sequencing was carried out to further shortening the QTL candidate interval for fruit size trait and excavated candidate genes for controlling fruit size.

Results

The candidate intervals related to jujube fruit size were mainly located on chromosomes 1, 5, and 10, and the frequency of chromosome 1 was the highest. Based on the QTL-seq results, the annotation results of ANNOVAR were extracted from 424 SNPs (single-nucleotide polymorphisms) and 164 InDels (insertion-deletion), from which 40 candidate genes were selected, and 37 annotated candidate genes were found in the jujube genome. Four genes (LOC107428904, LOC107415626, LOC125420708, and LOC107418290) that are associated with fruit size growth and development were identified by functional annotation of the genes in NCBI (National Center for Biotechnology Information). The genes can provide a basis for further exploration and identification on genes regulating jujube fruit size.

Discussion

In summary, the data obtained in this study revealed that QTL intervals and candidate genes for fruit size at the genomic level provide valuable resources for future functional studies and jujube breeding.

Comparative Genomic Analysis of PEBP Genes in Cucurbits Explores the Interactors of Cucumber CsPEBPs Related to Flowering Time

The family of phosphatidylethanolamine-binding proteins (PEBPs) participates in various plant biological processes, mainly flowering regulation and seed germination. In cucurbit crops, several PEBP genes have been recognized to be responsible for flowering time. However, the investigation of PEBP family members across the genomes of cucurbit species has not been reported, and their conservation and divergence in structure and function remain largely unclear. Herein, PEBP genes were identified from seven cucurbit crops and were used to perform a comparative genomics analysis. The cucurbit PEBP proteins could be classified into MFT, FT, TFL, and PEBP clades, and further, the TFL clade was divided into BFT-like, CEN-like, and TFL1-like subclades. The MFT-like, FT-like, and TFL-like proteins were clearly distinguished by a critical amino acid residue at the 85th position of the Arabidopsis FT protein. In gene expression analysis, CsaPEBP1 was highly expressed in flowers, and its expression levels in females and males were 70.5 and 89.2 times higher, respectively, than those in leaves. CsaPEBP5, CsaPEBP6, and CsaPEBP7 were specifically expressed in male flowers, with expression levels 58.1, 17.3, and 15.7 times higher, respectively, than those of leaves. At least five CsaPEBP genes exhibited the highest expression during the later stages of corolla opening. Through clustering of time-series-based RNA-seq data, several potential transcription factors (TFs) interacting with four CsaPEBPs were identified during cucumber corolla opening. Because of the tandem repeats of binding sites in promoters, NF-YB (Csa4G037610) and GATA (Csa7G64580) TFs appeared to be better able to regulate the CsaPEBP2 and CsaPEBP5 genes, respectively. This study would provide helpful information for further investigating the roles of PEBP genes and their interacting TFs in growth and development processes, such as flowering time regulation in cucurbit crops.

Delineation of genes for a major QTL governing heat stress tolerance in chickpea

Chickpea (Cicer arietinum) is a cool season grain legume experiencing severe yield loss during heat stress due to the intensifying climate changes and its associated gradual increase of mean temperature. Hence, understanding the genetic architecture regulating heat stress tolerance has emerged as an important trait to be addressed for enhancing yield and productivity of chickpea under heat stress. The present study is intended to identify the major genomic region(s) governing heat stress tolerance in chickpea. For this, an integrated genomics-assisted breeding strategy involving NGS-based high-resolution QTL-seq assay, QTL region-specific association analysis and molecular haplotyping was deployed in a population of 206 mapping individuals and a diversity panel of 217 germplasm accessions of chickpea. This combinatorial strategy delineated a major 156.8 kb QTL genomic region, which was subsequently narrowed-down to a functional candidate gene CaHSFA5 and its natural alleles associated strongly with heat stress tolerance in chickpea. Superior natural alleles and haplotypes delineated from the CaHSFA5 gene have functional significance in regulating heat stress tolerance in chickpea. Histochemical staining, interaction studies along with differential expression profiling of CaHSFA5 and ROS scavenging genes suggest a cross talk between CaHSFA5 with ROS homeostasis pertaining to heat stress tolerance in chickpea. Heterologous gene expression followed by heat stress screening further validated the functional significance of CaHSFA5 for heat stress tolerance. The salient outcomes obtained here can have potential to accelerate multiple translational genomic analysis including marker-assisted breeding and gene editing in order to develop high-yielding heat stress tolerant chickpea varieties.

Map-based cloning and functional analysis of a major quantitative trait locus, BolC.Pb9.1, controlling clubroot resistance in a wild Brassica relative (Brassica macrocarpa)

KEY MESSAGE

A causal gene BoUGT76C2, conferring clubroot resistance in wild Brassica oleracea, was identified and functionally characterized. Clubroot is a devastating soil-borne disease caused by the obligate biotrophic pathogen Plasmodiophora brassica (P. brassicae), which poses a great threat to Brassica oleracea (B. oleracea) production. Although several QTLs associated with clubroot resistance (CR) have been mapped in cultivated B. oleracea, none have been cloned in B. oleracea. Previously, we found that the wild B. oleracea B2013 showed high resistance to clubroot. In this study, we constructed populations using B2013 and broccoli line 90196. CR in B2013 is quantitatively inherited, and a major QTL, BolC.Pb9.1, was identified on C09 using QTL-seq and linkage analysis. The BolC.Pb9.1 was finely mapped to a 56 kb genomic region using F2:3 populations. From the target region, the candidate BoUGT76C2 showed nucleotide variations between the parents, and was inducible in response to P. brassicae infection. We generated BoUGT76C2 overexpression lines in the 90196 background, which showed significantly enhanced resistance to P. brassicae compared to the WT line, suggesting that BoUGT76C2 corresponds to the resistance gene BolC.Pb.9.1. This is the first report on the CR gene map-based cloning and functional analysis from wild relatives, which provides a theoretical basis to the understanding of the molecular mechanism of CR, and lays a foundation to improve the CR of cultivated B. oleracea.

The chickpea WIP2 gene underlying a major QTL contributes to lateral root development

Lateral roots are a major component of root system architecture, and lateral root count (LRC) positively contributes to yield under drought in chickpea. To understand the genetic regulation of LRC, a biparental mapping population derived from two chickpea accessions having contrasting LRCs was genotyped by sequencing, and phenotyped to map four major quantitative trait loci (QTLs) contributing to 13-32% of the LRC trait variation. A single- nucleotide polymorphism tightly linked to the locus contributing to highest trait variation was located on the coding region of a gene (CaWIP2), orthologous to NO TRANSMITTING TRACT/WIP domain protein 2 (NTT/WIP2) gene of Arabidopsis thaliana. A polymorphic simple sequence repeat (SSR) in the CaWIP2 promoter showed differentiation between low versus high LRC parents and mapping individuals, suggesting its utility for marker-assisted selection. CaWIP2 promoter showed strong expression in chickpea apical root meristem and lateral root primordia. Expression of CaWIP2 under its native promoter in the Arabidopsis wip2wip4wip5 mutant rescued its rootless phenotype to produce more lateral roots than the wild-type plants, and led to formation of amyloplasts in the columella. CaWIP2 expression also induced the expression of genes that regulate lateral root emergence. Our study identified a gene-based marker for LRC which will be useful for developing drought-tolerant, high-yielding chickpea varieties.

Deploying QTL-seq rapid identification and separation of the major QTLs of tassel branch number for fine-mapping in advanced maize populations

The tassel competes with the ear for nutrients and shields the upper leaves, thereby reducing the yield of grain. The tassel branch number (TBN) is a pivotal determinant of tassel size, wherein the reduced TBN has the potential to enhance the transmission of light and reduce the consumption of nutrients, which should ultimately result in increased yield. Consequently, the TBN has emerged as a vital target trait in contemporary breeding programs that focus on compact maize varieties. In this study, QTL-seq technology and advanced population mapping were used to rapidly identify and dissect the major effects of the TBN on QTL. Advanced mapping populations (BC4F2 and BC4F3) were derived from the inbred lines 18-599 (8-11 TBN) and 3237 (0-1 TBN) through phenotypic recurrent selection. First, 13 genomic regions associated with the TBN were detected using quantitative trait locus (QTL)-seq and were located on chromosomes 2 and 5. Subsequently, validated loci within these regions were identified by QTL-seq. Three QTLs for TBN were identified in the BC4F2 populations by traditional QTL mapping, with each QTL explaining the phenotypic variation of 6.13-18.17%. In addition, for the major QTL (qTBN2-2 and qTBN5-1), residual heterozygous lines (RHLs) were developed from the BC4F2 population. These two major QTLs were verified in the RHLs by QTL mapping, with the phenotypic variation explained (PVE) of 21.57% and 30.75%, respectively. Near-isogenic lines (NILs) of qTBN2-2 and qTBN5-1 were constructed. There were significant differences between the NILs in TBN. These results will enhance our understanding of the genetic basis of TBN and provide a solid foundation for the fine-mapping of TBN.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01431-y.

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.

BSA-Seq for the Identification of Major Genes for EPN in Rice

Improving rice yield is one of the most important food issues internationally. It is an undeniable goal of rice breeding, and the effective panicle number (EPN) is a key factor determining rice yield. Increasing the EPN in rice is a major way to increase rice yield. Currently, the main quantitative trait locus (QTL) for EPN in rice is limited, and there is also limited research on the gene for EPN in rice. Therefore, the excavation and analysis of major genes related to EPN in rice is of great significance for molecular breeding and yield improvement. This study used japonica rice varieties Dongfu 114 and Longyang 11 to construct an F5 population consisting of 309 individual plants. Two extreme phenotypic pools were constructed by identifying the EPN of the population, and QTL-seq analysis was performed to obtain three main effective QTL intervals for EPN. This analysis also helped to screen out 34 candidate genes. Then, EPN time expression pattern analysis was performed on these 34 genes to screen out six candidate genes with higher expression levels. Using a 3K database to perform haplotype analysis on these six genes, we selected haplotypes with significant differences in EPN. Finally, five candidate genes related to EPN were obtained.

WEAK SEED DORMANCY 1, an aminotransferase protein, regulates seed dormancy in rice through the GA and ABA pathways

Seed dormancy is a critical trait that enhances plant survival by preventing seed germination at the wrong time or under unsuitable conditions. Lack of seed dormancy in rice can lead to pre-harvest sprouting on mother plants leading to reduced yield and seed quality. Although some genes have been identified, knowledge of regulation of seed dormancy is limited. Here, we characterized a weak seed dormancy mutant named weak seed dormancy 1 (wsd1) that showed a higher seed germination percentage than the wild-type following the harvest ripeness. We cloned the WSD1 encoding an aminotransferase protein using a MutMap approach. WSD1 was stably expressed after imbibition and its protein was localized in the endoplasm reticulum. A widely targeted metabolomics assay and amino acid analysis showed that WSD1 had a role in regulating homeostasis of amino acids. PAC treatment and RNA-seq analysis showed that WSD1 regulates seed dormancy by involvement in the GA biosynthesis pathway. GA1 content and expression of GA biosynthesis-related genes were increased in the wsd1 mutant compared with the wild-type. The wsd1 mutant had reduced sensitivity to ABA. Our overall results indicated that WSD1 regulates seed dormancy by balancing the ABA and GA pathways.

Mitochondrial Inner Membrane ABC Transporter Bcmdl1 Is Involved in Conidial Germination, Virulence, and Resistance to Anilinopyrimidine Fungicides in Botrytis cinerea

Botrytis cinerea causes gray mold on thousands of plants, leading to huge losses in production. Anilinopyrimidine (AP) fungicides have been applied to control B. cinerea since the 1990s. Although resistance to AP fungicides was detected soon after their application, the mechanism of AP resistance remains to be elucidated. In this study, a sexual cross between resistant and sensitive isolates was performed, and the genomes of parental isolates and progenies were sequenced to identify resistance-related single nucleotide polymorphisms (SNPs). After screening and verification, mutation E407K in the Bcmdl1 gene was identified and confirmed to confer resistance to AP fungicides in B. cinerea. Bcmdl1 was predicted to encode a mitochondrial protein that belonged to a half-type ATP-binding cassette (ABC) transporter. Although Bcmdl1 was a transporter, it did not mediate resistance to multiple fungicides but mediated resistance specifically to AP fungicides. On the other hand, reductions in conidial germination and virulence were observed in Bcmdl1 knockout transformants compared to the parental isolate and complemented transformants, illustrating the biological functions of Bcmdl1. Subcellular localization analysis indicated that Bcmdl1 was localized in mitochondria. Interestingly, the production of ATP was reduced after cyprodinil treatment in Bcmdl1 knockout transformants, suggesting that Bcmdl1 was involved in ATP synthesis. Since Mdl1 could interact with ATP synthase in yeast, we hypothesize that Bcmdl1 forms a complex with ATP synthase, which AP fungicides might target, thereby interfering with the metabolism of energy. IMPORTANCE Gray mold, caused by B. cinerea, causes huge losses in the production of many fruits and vegetables. AP fungicides have been largely adopted to control this disease since the 1990s, and the development of resistance to AP fungicides initiates new problems for disease control. Due to the unknown mode of action, information on the mechanism of AP resistance is also limited. Recently, mutations in mitochondrial genes were reported to be related to AP resistance. However, the mitochondrial process of these genes remains to be elucidated. In this study, we identified several AP resistance-related mutations by quantitative trait locus sequencing (QTL-seq) and confirmed that mutation E407K in Bcmdl1 conferred AP resistance. We further characterized the expression patterns, biological functions, subcellular localization, and mitochondrial processes of the Bcmdl1 gene. This study deepens our understanding of the mechanism of resistance to and mode of action of AP fungicides.

QTL Mapping for Ovary- and Fruit-Related Traits in Cucumis sativus-C. hystrix Introgression Line IL52

IL52 is a valuable introgression line obtained from interspecific hybridization between cultivated cucumber (Cucumis sativus L., 2n = 14) and the wild relative species C. hystrix Chakr. (2n = 24). IL52 exhibits high resistance to a number of diseases, including downy mildew, powdery mildew, and angular leaf spot. However, the ovary- and fruit-related traits of IL52 have not been thoroughly investigated. Here, we conducted quantitative trait loci (QTL) mapping for 11 traits related to ovary size, fruit size, and flowering time using a previously developed 155 F_7:8 RIL population derived from a cross between CCMC and IL52. In total, 27 QTL associated with the 11 traits were detected, distributed on seven chromosomes. These QTL explained 3.61% to 43.98% of the phenotypic variance. Notably, we identified a major-effect QTL (qOHN4.1) on chromosome 4 associated with the ovary hypanthium neck width and further delimited it into a 114-kb candidate region harboring 13 candidate genes. Furthermore, the QTL qOHN4.1 is co-localized with the QTL detected for ovary length, mature fruit length, and fruit neck length, all residing within the consensus QTL FS4.1, suggesting a plausible pleiotropic effect.

Combined QTL-sequencing, linkage mapping, and RNA-sequencing identify candidate genes and KASP markers for low-temperature germination in Oryza sativa L. ssp. Japonica

MAIN CONCLUSION

Through QTL-seq, QTL mapping and RNA-seq, six candidate genes of qLTG9 can be used as targets for cold tolerance functional characterization, and six KASP markers can be used for marker-assisted breeding to improve the germination ability of japonica rice at low temperature. The development of direct-seeded rice at high latitudes and altitudes depends on the seed germination ability of rice under a low-temperature environment. However, the lack of regulatory genes for low-temperature germination has severely limited the application of genetics in improving the breeds. Here, we used cultivars DN430 and DF104 with significantly different low-temperature germination (LTG) and 460 F_2:3 progeny derived from them to identify LTG regulators by combining QTL-sequencing, linkage mapping, and RNA-sequencing. The QTL-sequencing mapped qLTG9 within a physical interval of 3.4 Mb. In addition, we used 10 Kompetitive allele-specific PCR (KASP) markers provided by the two parents, and qLTG9 was optimized from 3.4 Mb to a physical interval of 397.9 kb and accounted for 20.4% of the phenotypic variation. RNA-sequencing identified qLTG9 as eight candidate genes with significantly different expression within the 397.9 kb interval, six of which possessed SNPs on the promoter and coding regions. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) completely validated the results of these six genes in RNA-sequencing. Subsequently, six non-synonymous SNPs were designed using variants in the coding region of these six candidates. Genotypic analysis of these SNPs in 60 individuals with extreme phenotypes indicated these SNPs determined the differences in cold tolerance between parents. The six candidate genes of qLTG9 and the six KASP markers could be used together for marker-assisted breeding to improve LTG.

Identification of a suppressor for the wheat stripe rust resistance gene Yr81 in Chinese wheat landrace Dahongpao

Combined with BSE-Seq analysis and multiple genetic populations, three genes involved in stripe rust resistance were identified in Chinese wheat landrace Dahongpao, including a novel suppressor on 2BS. Dahongpao (DHP), a landrace of hexaploid wheat in China, exhibits a high degree of stripe rust resistance in the field for many years. In this study, bulked segregant analysis coupled with exome capture sequencing (BSE-Seq) was used to identify genes encoding stripe rust resistance in multiple genetic populations from the cross between DHP and a susceptible hexaploid Australian cultivar, Avocet S (AvS). The most effective QTL in DHP was Yr18, explaining up to 53.08% of phenotypic variance in the F_2:3 families. To identify additional genes, secondary mapping populations SP1 and SP2 were produced by crossing AvS with two resistant lines derived from F_2:3 families lacking Yr18. An all-stage resistance gene, Yr.DHP-6AS, was identified via BSE-Seq analysis of SP1. Combined the recombinant plants from both SP1 and SP2, Yr.DHP-6AS was located between KP6A_1.66 and KP6A_8.18, corresponding to the same region as Yr81. In addition, secondary mapping populations SP3 and SP4 were developed by selfing a segregating line from F_2:3 families lacking Yr18. A novel suppressor gene on chromosome 2BS was identified from DHP for effectively suppressing the resistance of Yr.DHP-6AS in the SP3 and SP4. As a result, the wheat lines carrying both Yr18 and Yr.DHP-6AS show higher level of stripe rust resistance than DHP, providing an effective and simple combination for developing new wheat cultivars with ASR and APR genes. Further, the newly developed KASP markers, KP6A_1.99 and KP6A_5.22, will facilitate the application of Yr.DHP-6AS in wheat breeding via marker-assisted selection.

Fine mapping and causal gene identification of a novel QTL for early flowering by QTL-seq, Target-seq and RNA-seq in spring oilseed rape

A novel quantitative trait locus for early flowering in spring oilseed rape, BnaC08cqDTF, was mapped to an 86-kb region on chromosome C08, and its causal gene, CRY2, was uncovered. Days to flowering is a very important agronomic and adaptive trait of Brassica napus oilseed rape (AACC, 2n = 38). We previously identified BnaC08cqDTF as a novel candidate quantitative trait locus (QTL) for early flowering in spring oilseed rape. Here, we present fine mapping of the locus and a study of its causal gene. Initial mapping was performed by QTL sequencing of DNA pools of BC₃F₂ plants with extreme flowering times derived from crosses between the spring-type cv. No. 4512 (early flowering) and cv. No. 5246 (late flowering), along with fine mapping by target sequencing of the BC₃F₂ and BC₄F₂ populations. Fine mapping narrowed down BnaC08cqDTF to an 86-kb region on chromosome C08. The region harbored fifteen genes. After comparative analyses of the DNA sequences for mutation between A and C syntenic regions and detected by RNA-seq and qRT-PCR between the two parents, we found that BnaC08G0010400ZS harbors an A/G nonsynonymous mutation in exon 3. This single nucleotide polymorphism (SNP) haplotype was also correlated with early flowering in a 256 accession panel. BnaC08G0010400ZS is a homolog of the AT1G04400 gene (CRY2) in Arabidopsis. The analyses of transgenic Arabidopsis verified that BnaC08G0010400ZS is responsible for early flowering. Our results contribute to a better understanding of the genetic control mechanism of early flowering in spring Brassica napus and will promote the breeding for early mature varieties.

A 5.5-kb LTR-retrotransposon insertion inside phytochrome B gene (CsPHYB) results in long hypocotyl and early flowering in cucumber (Cucumis sativus L.)

The novel spontaneous long hypocotyl and early flowering (lhef) mutation in cucumber is due to a 5551-bp LTR-retrotransposon insertion in CsPHYB gene encoding PHYTOCHROME B, which plays a major role in regulating photomorphogenic hypocotyl growth and flowering. Hypocotyl length and flowering time are important for establishing high-quality seedlings in modern cucumber production, but little is known for the underlying molecular mechanisms of these two traits. In this study, a spontaneous cucumber long hypocotyl and early flowering mutant was identified and characterized. Based on multiple lines of evidence, we show that cucumber phytochrome B (CsPHYB) is the candidate gene for this mutation, and a 5551-bp LTR-retrotransposon insertion in the first exon of CsPHYB was responsible for the mutant phenotypes. Uniqueness of the mutant allele at CsPHYB was verified in 114 natural cucumber lines. Ectopic expression of the CsPHYB in Arabidopsis phyB mutant rescued the long hypocotyl and early flowering phenotype of phyB-9 mutant. The wild-type CsPHYB protein was localized on the membrane and cytoplasm under white light condition, whereas in the nucleus under red light, it is consistent with its roles as a red-light photoreceptor in Arabidopsis. However, the mutant csphyb protein was localized on the membrane and cytoplasm under both white and red-light conditions. Expression dynamics of CsPHYB and several cell elongation-related genes were positively correlated with hypocotyl elongation; the transcription levels of key positive and negative regulators for flowering time were also consistent with the anthesis dates in the mutant and wild-type plants. Yeast two hybrid and bimolecular fluorescence complementation assays identified physical interactions between CsPHYB and phytochrome interacting factor 3/4 (CsPIF3/4). These findings will provide new insights into the roles of the CsPHYB in cucumber hypocotyl growth and flowering time.

Molecular mapping of genomic regions and identification of possible candidate genes associated with gynoecious sex expression in bitter gourd

Bitter gourd is an important vegetable crop grown throughout the tropics mainly because of its high nutritional value. Sex expression and identification of gynoecious trait in cucurbitaceous vegetable crops has facilitated the hybrid breeding programme in a great way to improve productivity. In bitter gourd, gynoecious sex expression is poorly reported and detailed molecular pathways involve yet to be studied. The present experiment was conducted to study the inheritance, identify the genomic regions associated with gynoecious sex expression and to reveal possible candidate genes through QTL-seq. Segregation for the gynoecious and monoecious sex forms in the F₂ progenies indicated single recessive gene controlling gynoecious sex expression in the genotype, PVGy-201. Gynoecious parent, PVGy-201, Monoecious parent, Pusa Do Mausami (PDM), and two contrasting bulks were constituted for deep-sequencing. A total of 10.56, 23.11, 15.07, and 19.38 Gb of clean reads from PVGy-201, PDM, gynoecious bulk and monoecious bulks were generated. Based on the ΔSNP index, 1.31 Mb regions on the chromosome 1 was identified to be associated with gynoecious sex expression in bitter gourd. In the QTL region 293,467 PVGy-201 unique variants, including SNPs and indels, were identified. In the identified QTL region, a total of 1019 homozygous variants were identified between PVGy1 and PDM genomes and 71 among them were non-synonymous variants (SNPS and INDELs), out of which 11 variants (7 INDELs, 4 SNPs) were classified as high impact variants with frame shift/stop gain effect. In total twelve genes associated with male and female gametophyte development were identified in the QTL-region. Ethylene-responsive transcription factor 12, Auxin response factor 6, Copper-transporting ATPase RAN1, CBL-interacting serine/threonine-protein kinase 23, ABC transporter C family member 2, DEAD-box ATP-dependent RNA helicase 1 isoform X2, Polygalacturonase QRT3-like isoform X2, Protein CHROMATIN REMODELING 4 were identified with possible role in gynoecious sex expression. Promoter region variation in 8 among the 12 genes indicated their role in determining gynoecious sex expression in bitter gourd genotype, DBGy-1. The findings in the study provides insight about sex expression in bitter gourd and will facilitate fine mapping and more precise identification of candidate genes through their functional validation.

Functional Bph14 from Rathu Heenati promotes resistance to BPH at the early seedling stage of rice (Oryza sativa L.) as revealed by QTL-seq

A QTL associated with BPH resistance at the early seedling stage was identified on chromosome 3. Functional Bph14 in Rathu Heenati was associated with BPH resistance at the early seedling stage. Brown planthopper (BPH; Nilaparvata lugens Stål) is considered the most important rice pest in many Asian countries. Several BPH resistance genes have previously been identified. However, there are few reports of genes specific for BPH resistance at the early seedling stage, a crucial stage for direct-seeding cultivation. In this study, we performed a QTL-seq analysis using two bulks (20 F₂ lines in each bulk) of the F₂ population (n = 300) derived from a cross of Rathu Heenati (RH) × HCS-1 to identify QTL/genes associated with BPH resistance at the early seedling stage. An important QTL was identified on chromosome 3 and Bph14 was identified as a potential candidate gene based on the differences in gene expression and sequence variation when compared with the two parents. All plants in the resistant bulks possessed the functional Bph14 from RH and all plants in the susceptible bulk and HCS-1 contained a large deletion (2703 bp) in Bph14. The functional Bph14 gene of RH appears to be important for BPH resistance at the early seedling stage of rice and could be used in conjunction with other BPH resistance genes in rice breeding programs that confer resistance to BPH at the early and later growth stages.

Quantitative trait locus mapping and identification of candidate genes for resistance to Verticillium wilt in four recombinant inbred line populations of Gossypium hirsutum

Improving resistance to Verticillium wilt is of great significance for achieving high and stable yields of Upland cotton (Gossypium hirsutum). To deeply understand the genetic basis of cotton resistance to Verticillium wilt, Verticillium wilt-resistant Upland Lumianyan 28 and four Verticillium wilt-susceptible Acala cotton cultivars were used to create four recombinant inbred line (RIL) populations of 469 families through nested hybridization. Phenotypic data collected in five stressful environments were used to select resistant and sensitive lines and create a mixed pool of extreme phenotypes for BSA-seq. A total of 8 QTLs associated with Verticillium wilt resistance were identified on 4 chromosomes, of which qVW-A12-5 was detected simultaneously in the RIL populations and in one of the RIL populations and was identified for the first time. According to the sequence comparison and transcriptome analysis of candidate genes in the QTL interval between parents and pools, 4 genes were identified in the qVW-A12-5 interval. qRT-PCR of parental and phenotypically extreme lines revealed that Gh_CPR30 was induced by and may be a candidate gene for resistance to Verticillium wilt in G. hirsutum. Furthermore, VIGS technology revealed that the disease severity index (DSI) of the Gh_CPR30-silenced plants was significantly higher than that of the control. These results indicate that the Gh_CPR30 gene plays an important role in the resistance of G. hirsutum to Verticillium wilt, and the study provides a molecular basis for analyzing the molecular mechanism underlying G. hirsutum resistance to Verticillium wilt.

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 F2 plants and derivative F2:3 and F2:4 populations by genotyping by sequencing (GBS). BC1F2 population was constructed by backcrossing one F2: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 F2:3 and F2: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.

Whole genome resequencing identifies candidate genes and allelic diagnostic markers for resistance to Ralstonia solanacearum infection in cultivated peanut (Arachis hypogaea L.)

Bacterial wilt disease (BWD), caused by Ralstonia solanacearum is a major challenge for peanut production in China and significantly affects global peanut field productivity. It is imperative to identify genetic loci and putative genes controlling resistance to R. solanacearum (RRS). Therefore, a sequencing-based trait mapping approach termed "QTL-seq" was applied to a recombination inbred line population of 581 individuals from the cross of Yueyou 92 (resistant) and Xinhuixiaoli (susceptible). A total of 381,642 homozygous single nucleotide polymorphisms (SNPs) and 98,918 InDels were identified through whole genome resequencing of resistant and susceptible parents for RRS. Using QTL-seq analysis, a candidate genomic region comprising of 7.2 Mb (1.8-9.0 Mb) was identified on chromosome 12 which was found to be significantly associated with RRS based on combined Euclidean Distance (ED) and SNP-index methods. This candidate genomic region had 180 nonsynonymous SNPs and 14 InDels that affected 75 and 11 putative candidate genes, respectively. Finally, eight nucleotide binding site leucine rich repeat (NBS-LRR) putative resistant genes were identified as the important candidate genes with high confidence. Two diagnostic SNP markers were validated and revealed high phenotypic variation in the different resistant and susceptible RIL lines. These findings advocate the expediency of the QTL-seq approach for precise and rapid identification of candidate genomic regions, and the development of diagnostic markers that are applicable in breeding disease-resistant peanut varieties.

Fine mapping and candidate gene analysis of gynoecy trait in chieh-qua (Benincasa hispida Cogn. var. chieh-qua How)

Gynoecy demonstrates an earlier production of hybrids and a higher yield and improves the efficiency of hybrid seed production. Therefore, the utilization of gynoecy is beneficial for the genetic breeding of chieh-qua. However, little knowledge of gynoecious-related genes in chieh-qua has been reported until now. Here, we used an F₂ population from the cross between the gynoecious line 'A36' and the monoecious line 'SX' for genetic mapping and revealed that chieh-qua gynoecy was regulated by a single recessive gene. We fine-mapped it into a 530-kb region flanked by the markers Indel-3 and KASP145 on Chr.8, which harbors eight candidate genes. One of the candidate genes, Bhi08G000345, encoding networked protein 4 (CqNET4), contained a non-synonymous SNP resulting in the amino acid substitution of isoleucine (ATA; I) to methionine (ATG; M). CqNET4 was prominently expressed in the female flower, and only three genes related to ethylene synthesis were significantly expressed between 'A36' and 'SX.' The results presented here provide support for the CqNET4 as the most likely candidate gene for chieh-qua gynoecy, which differed from the reported gynoecious genes.

Identification of a major QTL, Parth6.1 associated with parthenocarpic fruit development in slicing cucumber genotype, Pusa Parthenocarpic Cucumber-6

Parthenocarpy is an extremely important trait that revolutionized the worldwide cultivation of cucumber under protected conditions. Pusa Parthenocarpic Cucumber-6 (PPC-6) is one of the important commercially cultivated varieties under protected conditions in India. Understanding the genetics of parthenocarpy, molecular mapping and the development of molecular markers closely associated with the trait will facilitate the introgression of parthenocarpic traits into non-conventional germplasm and elite varieties. The F1, F2 and back-crosses progenies with a non-parthenocarpic genotype, Pusa Uday indicated a single incomplete dominant gene controlling parthenocarpy in PPC-6. QTL-seq comprising of the early parthenocarpy and non-parthenocarpic bulks along with the parental lines identified two major genomic regions, one each in chromosome 3 and chromosome 6 spanning over a region of 2.7 Mb and 7.8 Mb, respectively. Conventional mapping using F2:3 population also identified two QTLs, Parth6.1 and Parth6.2 in chromosome 6 which indicated the presence of a major effect QTL in chromosome 6 determining parthenocarpy in PPC-6. The flanking markers, SSR01148 and SSR 01012 for Parth6.1 locus and SSR10476 and SSR 19174 for Parth6.2 locus were identified and can be used for introgression of parthenocarpy through the marker-assisted back-crossing programme. Functional annotation of the QTL-region identified two major genes, Csa_6G396640 and Csa_6G405890 designated as probable indole-3-pyruvate monooxygenase YUCCA11 and Auxin response factor 16, respectively associated with auxin biosynthesis as potential candidate genes. Csa_6G396640 showed only one insertion at position 2179 in the non-parthenocarpic parent. In the case of Csa_6G405890, more variations were observed between the two parents in the form of SNPs and InDels. The study provides insight about genomic regions, closely associated markers and possible candidate genes associated with parthenocarpy in PPC-6 which will be instrumental for functional genomics study and better understanding of parthenocarpy in cucumber.

Analysis of QTLs and Candidate Genes for Tassel Symptoms in Maize Infected with Sporisorium reilianum

Heat smut is a fungal soil-borne disease caused by Sporisorium reilianum, and affects the development of male and female tassels. Our previous research found that the tassel symptoms in maize infected with Sporisorium reilianum significantly differed in inbred lines with Sipingtou blood, and exhibited stable heredity over time at multiple locations. In this study, cytological analysis demonstrated that the cellular organization structures of three typical inbred lines (Huangzao4, Jing7, and Chang7-2) showed significant discrepancies at the VT stage. QTLs that control the different symptoms of maize tassels infected with Sporisorium reilianum were located in two F₂ populations, which were constructed using three typical inbred lines. The BSA (bulked segregation analysis) method was used to construct mixed gene pools based on typical tassel symptoms. The QTLs of different symptoms of maize tassels infected with Sporisorium reilianum were detected with 869 SSR markers covering the whole maize genome. The mixed gene pools were screened with polymorphic markers between the parents. Additional SSR markers were added near the above marker to detect genotypes in partially single plants in F₂ populations. The QTL controlling tassel symptoms in the Huangzao4 and Jing7 lines was located on the bin 1.06 region, between the markers of umc1590 and bnlg1598, and explained 21.12% of the phenotypic variation with an additive effect of 0.6524. The QTL controlling the tassel symptoms of the Jing7 and Chang7-2 lines was located on the bin 2.07 region, between the markers of umc1042 and bnlg1335, and explained 11.26% phenotypic variation with an additive effect of 0.4355. Two candidate genes (ZmABP2 and Zm00001D006403) were identified by a conjoint analysis of label-free quantification proteome sequencings.

QTL analysis of tuber shape in a diploid potato population

Tuber shape is one of the most important traits for potato breeding. Since poor or irregular shape increases the difficulty of handling and processing, researching the inheritance of potato tuber shape for potato breeding is highly important. To efficiently identify QTL for tuber shape, a diploid potato population (PM7) was generated by self-pollinated M6 (S. chacoense). A QTL TScha6 for tuber shape was identified by the QTL-seq approach at 50.91-59.93 Mb on chromosome 6 in the potato DM reference genome. To confirm TScha6, four SSR and twenty CAPS markers around the QTL were developed and the TScha6 was narrowed down to an interval of ~ 1.85 Mb. The CAPS marker C6-58.27_665 linked to TScha6 was then used to screen 86 potato cultivars and advanced breeding lines. The tuber length/width (LW) ratio was significantly correlated with the presence/absence of C6-58.27_665, and the correlation coefficient was r = 0.55 (p < 0.01). These results showed that C6-58.27_665 could be applied in marker-assisted selection (MAS) for tuber shape breeding in the future. Our research sets the important stage for the future cloning of the tuber shape gene and utilities of the marker in the breeding program.

Ribonuclease H-like gene SMALL GRAIN2 regulates grain size in rice through brassinosteroid signaling pathway

Grain size is a key agronomic trait that determines the yield in plants. Regulation of grain size by brassinosteroids (BRs) in rice has been widely reported. However, the relationship between the BR signaling pathway and grain size still requires further study. Here, we isolated a rice mutant, named small grain2 (sg2), which displayed smaller grain and a semi-dwarf phenotype. The decreased grain size was caused by repressed cell expansion in spikelet hulls of the sg2 mutant. Using map-based cloning combined with a MutMap approach, we cloned SG2, which encodes a plant-specific protein with a ribonuclease H-like domain. SG2 is a positive regulator downstream of GLYCOGEN SYNTHASE KINASE2 (GSK2) in response to BR signaling, and its mutation causes insensitivity to exogenous BR treatment. Genetical and biochemical analysis showed that GSK2 interacts with and phosphorylates SG2. We further found that BRs enhance the accumulation of SG2 in the nucleus, and subcellular distribution of SG2 is regulated by GSK2 kinase activity. In addition, Oryza sativa OVATE family protein 19 (OsOFP19), a negative regulator of grain shape, interacts with SG2 and plays an antagonistic role with SG2 in controlling gene expression and grain size. Our results indicated that SG2 is a new component of GSK2-related BR signaling response and regulates grain size by interacting with OsOFP19.

Nitrate signaling and use efficiency in crops

Nitrate (NO₃^-) is not only an essential nutrient but also an important signaling molecule for plant growth. Low nitrogen use efficiency (NUE) of crops is causing increasingly serious environmental and ecological problems. Understanding the molecular mechanisms of NO₃^- regulation in crops is crucial for NUE improvement in agriculture. During the last several years, significant progress has been made in understanding the regulation of NO₃^- signaling in crops, and some key NO₃^- signaling factors have been shown to play important roles in NO₃^- utilization. However, no detailed reviews have yet summarized these advances. Here, we focus mainly on recent advances in crop NO₃^- signaling, including short-term signaling, long-term signaling, and the impact of environmental factors. We also review the regulation of crop NUE by crucial genes involved in NO₃^- signaling. This review provides useful information for further research on NO₃^- signaling in crops and a theoretical basis for breeding new crop varieties with high NUE, which has great significance for sustainable agriculture.

Discovery of stripe rust resistance with incomplete dominance in wild emmer wheat using bulked segregant analysis sequencing

Durable crop disease resistance is an essential component of global food security. Continuous pathogen evolution leads to a breakdown of resistance and there is a pressing need to characterize new resistance genes for use in plant breeding. Here we identified an accession of wild emmer wheat (Triticum turgidum ssp. dicoccoides), PI 487260, that is highly resistant to multiple stripe rust isolates. Genetic analysis revealed resistance was conferred by a single, incompletely dominant gene designated as Yr84. Through bulked segregant analysis sequencing (BSA-Seq) we identified a 52.7 Mb resistance-associated interval on chromosome 1BS. Detected variants were used to design genetic markers for recombinant screening, further refining the interval of Yr84 to a 2.3-3.3 Mb in tetraploid wheat genomes. This interval contains 34 candidate genes encoding for protein domains involved in disease resistance responses. Furthermore, KASP markers closely-linked to Yr84 were developed to facilitate marker-assisted selection for rust resistance breeding.

Elucidation of the Flavor Aspects and Flavor-Associated Genomic Regions in Bottle Gourd (Lagenaria siceraria) by Metabolomic Analysis and QTL-seq

Bottle gourd (Lagenaria siceraria) is a commercially important cucurbitaceous vegetable with health-promoting properties whose collections and cultivars differ considerably in their flavor aspects. However, the metabolomic profile related to flavor has not yet been elucidated. In the present study, a comprehensive metabolite analysis revealed the metabolite profile of the strong-flavor collection “J120” and weak-flavor collection “G32”. The major differentially expressed metabolites included carboxylic acids, their derivatives, and organooxygen compounds, which are involved in amino acid biosynthesis and metabolism. QTL-seq was used to identify candidate genomic regions controlling flavor in a MAGIC population comprising 377 elite lines. Three significant genomic regions were identified, and candidate genes likely associated with flavor were screened. Our study provides useful information for understanding the metabolic causes of flavor variation among bottle gourd collections and cultivars. Furthermore, the identified candidate genomic regions may facilitate rational breeding programs to improve bottle gourd quality.

Identification of single major QTL and candidate gene(s) governing hull-less seed trait in pumpkin

The hull-less pumpkin (Cucurbita pepo) seed does not require de-hulling before use for human consumption, as a result highly preferred by the oil, nut, and baking industries. In hull-less seeds, a single recessive gene is responsible for the absence of outer thick seed coat layers; however, the genomic region and gene(s) controlling the trait are unclear to date. In this study, four crosses attempted to derive F₂ and backcross populations confirmed the single recessive gene inheritance of hull-less seed trait in populations adapted to the sub-tropical climate. The candidate genomic region for hull-less seed trait was identified through the BSA-QTLseq approach using bulks of F_2:3 progenies from a cross of HP111 (hulled) and HLP36 (hull-less). A novel genomic region on chromosome 12 ranging from 1.80 to 3.86 Mb was associated with the hull-less seed trait. The re-sequencing data identified a total of 396 SNPs within this region and eight were successfully converted into polymorphic KASP assays. The genotyping of segregating F₂ (n = 160) with polymorphic KASP assays resulted in a 40.3 cM partial linkage map and identified Cp_3430407 (10 cM) and Cp_3498687 (16.1 cM) as flanking markers for hull-less locus (Cphl-1). These flanking markers correspond to the 68.28 kb region in the reference genome, and the marker, Cp_3430407 successfully predicted the genotype in 93.33% of the C. pepo hull-less germplasm lines, thus can be used for marker-assisted selection in parents polymorphic for the hull-less seed trait. The Cphl-1-linked genomic region (2.06 Mb) encompasses a total of 182 genes, including secondary cell wall and lignin biosynthesis-related transcriptional factors viz., "NAC" (Cp4.1LG12g04350) and "MYB" (Cp4.1LG12g03120). These genes were differentially expressed in the seeds of hulled and hull-less genotypes, and therefore could be the potential candidate genes governing the hull-less seed trait in pumpkin.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

Modified photoperiod response of CsFT promotes day neutrality and early flowering in cultivated cucumber

Map-based cloning and photoperiod response detection suggested that CsFT is the critical gene for cucumber photoperiod domestication. Photoperiod sensitivity is important for sensing seasonal changes and local adaptation. However, day-length sensitivity limits crop geographical adaptation and it should be modified during domestication. Cucumber was domesticated in southern Asia and is currently cultivated worldwide across a wide range of latitudes, but its photoperiod sensitivity and its change during cucumber domestication are unknown. Here, we confirmed wild cucumber (Hardwickii) was a short-day plant, and its flowering depends on short-day (SD) conditions, while the cultivated cucumber (9930) is a day-neutral plant that flowers independently of day length. A photoperiod sensitivity locus (ps-1) was identified by the 9930 × Hardwickii F₂ segregating populations, which span a ~ 970 kb region and contain 60 predicted genes. RNA-seq analysis showed that the critical photoperiod pathway gene FLOWERING LOCUS T (CsFT) within the ps-1 locus exhibits differential expression between 9930 and Hardwickii, which was confirmed by qRT-PCR detection. CsFT in Hardwickii was sensitive to day length and could be significantly induced by SD conditions, whereas CsFT was highly expressed in 9930 and was insensitive to day length. Moreover, the role of CsFT in promoting flowering was verified by overexpression of CsFT in Arabidopsis. We also identified the genetic variations existing in the promoter of CsFT among the different geographic cucumbers and suggest they have possible roles in photoperiod domestication. The results of this study suggest that a variation in photoperiod sensitivity of CsFT is associated with day neutrality and early flowering in cultivated cucumber and could contribute to cucumber cultivation in diverse regions throughout the world.

Evaluation of nine statistics to identify QTLs in bulk segregant analysis using next generation sequencing approaches

Background

Bulk segregant analysis (BSA) combined with next generation sequencing is a powerful tool to identify quantitative trait loci (QTL). The impact of the size of the study population and the percentage of extreme genotypes analysed have already been assessed. But a good comparison of statistical approaches designed to identify QTL regions using next generation sequencing (NGS) technologies for BSA is still lacking.

Results

We developed an R code to simulate QTLs in bulks of F2 contrasted lines. We simulated a range of recombination rates based on estimations using different crop species. The simulations were used to benchmark the ability of statistical methods identify the exact location of true QTLs. A single QTL led to a shift in allele frequency across a large fraction of the chromosome for plant species with low recombination rate. The smoothed version of all statistics performed best notably the smoothed Euclidean distance-based statistics was always found to be more accurate in identifying the location of QTLs. We propose a simulation approach to build confidence interval statistics for the detection of QTLs.

Conclusion

We highlight the statistical methods best suited for BSA studies using NGS technologies in crops even when recombination rate is low. We also provide simulation codes to build confidence intervals and to assess the impact of recombination for application to other studies. This computational study will help select NGS-based BSA statistics that are useful to the broad scientific community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08718-y.

Fine mapping and analysis of candidate genes for qFT7.1, a major quantitative trait locus controlling flowering time in Brassica rapa L

qFT7.1, a major QTL for flowering time in Brassica rapa was fine-mapped to chromosome A07 in a 56.4-kb interval, in which the most likely candidate gene is BraA07g018240.3C. In Brassica rapa, flowering time (FT) is an important agronomic trait that affects the yield, quality, and adaption. FT is a complicated trait that is regulated by many genes and is affected greatly by the environment. In this study, a chromosome segment substitution line (CSSL), CSSL16, was selected that showed later flowering than the recurrent parent, a rapid-cycling inbred line of B. rapa (RcBr). Using Bulked Segregant RNA sequencing, we identified a late flowering quantitative trait locus (QTL), designated as qFT7.1, on chromosome A07, based on a secondary-F₂ population derived from the cross between CSSL16 and RcBr. qFT7.1 was further validated by conventional QTL mapping. This QTL explained 39.9% (logarithm of odds = 32.2) of the phenotypic variations and was fine mapped to a 56.4-kb interval using recombinant analysis. Expression analysis suggested that BraA07g018240.3C, which is homologous to ATC (encoding Arabidopsis thaliana CENTRORADIALIS homologue), a gene for delayed flowering in Arabidopsis, as the most promising candidate gene. Sequence analysis demonstrated that two synonymous mutations existed in the coding region and numerous bases replacements existed in promoter region between BraA07g018240.3C from CSSL16 and RcBr. The results will increase our knowledge related to the molecular mechanism of late flowering in B. rapa and lays a solid foundation for the breeding of late bolting B. rapa.

OsWRKY115 on qCT7 links to cold tolerance in rice

qCT7, a novel QTL for increasing seedling cold tolerance in rice, was fine-mapped to a 70.9-kb region on chromosome 7, and key OsWRKY115 was identified in transgenic plants. Cold stress caused by underground cold-water irrigation seriously limits rice productivity. We systemically measured the cold-responsive traits of 2,570 F2 individuals derived from two widely cultivated rice cultivars, Kong-Yu-131 and Dong-Nong-422, to identify the major genomic regions associated with cold tolerance. A novel major QTL, qCT7, was mapped on chromosome 7 associated with the cold tolerance and survival, using whole-genome re-sequencing with bulked segregant analysis. Local QTL linkage analysis with F₂ and fine mapping with recombinant plant revealed a 70.9-kb core region on qCT7 encoding 13 protein-coding genes. Only the LOC_Os07g27670 expression level encoding the OsWRKY115 transcription factor on the locus was specifically induced by cold stress in the cold-tolerant cultivar. Moreover, haplotype analysis and the KASP8 marker indicated that OsWRKY115 was significantly associated with cold tolerance. Overexpression and knockout of OsWRKY115 significantly affected cold tolerance in seedlings. Our experiments identified OsWRKY115 as a novel regulatory gene associated with cold response in rice, and the Kong-Yu-131 allele with specific cold-induced expression may be an important molecular variant.

The New Variation in the Promoter Region of FLOWERING LOCUS T Is Involved in Flowering in Brassica rapa

Flowering time is an important agronomic trait in Brassica rapa and has a wide range of variation. The change from vegetative to reproductive development is a major transition period, especially in flowering vegetable crops. In this study, two non-heading Chinese cabbage varieties with significantly different flowering times, Pak-choi (B. rapa var. communis Tesn et Lee) and Caitai (B. rapa var. tsaitai Hort.), were used to construct segregated F₂ populations. The bulk-segregant approach coupled with whole genome re-sequencing was used for QTL sequencing (QTL-seq) analysis to map flowering time traits. The candidate genes controlling flowering time in B. rapa were predicted by homologous gene alignment and function annotation. The major-effect QTL ft7.1 was detected on chromosome A07 of B. rapa, and the FT family gene BrFT was predicted as the candidate gene. Moreover, a new promoter regional difference of 1577 bp was revealed by analyzing the sequence of the BrFT gene. The promoter region activity analysis and divergent gene expression levels indicated that the difference in the promoter region may contribute to different flowering times. These findings provide insights into the mechanisms underlying the flowering time in Brassica and the candidate genes regulating flowering in production.

Single trait versus principal component based association analysis for flowering related traits in pigeonpea

Pigeonpea, a tropical photosensitive crop, harbors significant diversity for days to flowering, but little is known about the genes that govern these differences. Our goal in the current study was to use genome wide association strategy to discover the loci that regulate days to flowering in pigeonpea. A single trait as well as a principal component based association study was conducted on a diverse collection of 142 pigeonpea lines for days to first and fifty percent of flowering over 3 years, besides plant height and number of seeds per pod. The analysis used seven association mapping models (GLM, MLM, MLMM, CMLM, EMLM, FarmCPU and SUPER) and further comparison revealed that FarmCPU is more robust in controlling both false positives and negatives as it incorporates multiple markers as covariates to eliminate confounding between testing marker and kinship. Cumulatively, a set of 22 SNPs were found to be associated with either days to first flowering (DOF), days to fifty percent flowering (DFF) or both, of which 15 were unique to trait based, 4 to PC based GWAS while 3 were shared by both. Because PC1 represents DOF, DFF and plant height (PH), four SNPs found associated to PC1 can be inferred as pleiotropic. A window of ± 2 kb of associated SNPs was aligned with available transcriptome data generated for transition from vegetative to reproductive phase in pigeonpea. Annotation analysis of these regions revealed presence of genes which might be involved in floral induction like Cytochrome p450 like Tata box binding protein, Auxin response factors, Pin like genes, F box protein, U box domain protein, chromatin remodelling complex protein, RNA methyltransferase. In summary, it appears that auxin responsive genes could be involved in regulating DOF and DFF as majority of the associated loci contained genes which are component of auxin signaling pathways in their vicinity. Overall, our findings indicates that the use of principal component analysis in GWAS is statistically more robust in terms of identifying genes and FarmCPU is a better choice compared to the other aforementioned models in dealing with both false positive and negative associations and thus can be used for traits with complex inheritance.

Genetic basis of the lobedness degree in tomato fruit morphology

Fruit shape is a key trait in tomato (Solanum lycopersicum L.). Since most studies focused on proximo-distal fruit morphology, we hypothesized that unknown QTLs for medio-lateral direction ones could be found analysing segregating populations where major shape genes are fixed. We examined the diversity of fruit morphology in medio-lateral direction; defined divergent traits in cultivars carrying identical genetic constitution at LC and FAS genes; and identified QTLs for lobedness degree (LD) by a QTL-seq approach. We found that LC and FAS genes were not enough to explain LD variability in a large tomato collection. Then, we derived F₂ populations crossing cultivars divergent for LD where LC and FAS were fixed (Yellow Stuffer x Heinz 1439 [F₂YSxH] and Voyage x Old Brooks [F₂VxOB]). By QTL-seq we identified a QTL for LD on chromosome 8 in both F₂, which was validated in F₂YSxH by interval mapping accounting for ~ 17% of the variability. Other two QTLs located on chromosomes 6 and 11 with epistasis explained ~ 61% of the variability in the F₂VxOB. In conclusion, three novel QTLs with major effect for LD (ld6, ld8, and ld11) were identified through the study of diversity and genetic segregation in intraspecific tomato crosses.

Recent status of Genotyping by Sequencing (GBS) Technology in cucumber (Cucumis sativus L.): a review

Current and advanced breeding tools are being used to improve economically important horticultural crops to meet the consumers' needs and preferences. Genotyping-by-sequencing (GBS) is an extremely useful tool in the investigation and analysis of the genetic diversity of different cultivars. Based on a broad range of genetic backgrounds like single nucleotide polymorphism (SNPs), GBS is known as a novel technique to facilitate the detection of quantitative trait loci (QTL) regions robustly linked with interested traits compared to genome-wide association study (GWAS) and QTL. GBS has gained popularity among breeders in recent years and it is also employed in cucumber breeding programs. Cucumbers (C. sativus L.) are monoecious, gynoecious and some of them are parthenocarpic species. Cucumber is one of the most economical and essential crops in the Cucurbitaceae family. For time immemorial, cucumber has been produced and consumed all over the world like other cucurbits. To a large extent, cultivated cucurbits are beneficial to human health for providing necessary minerals and fibers.Therefore, this review portrays the current status of advances made by using GBS and its combination with other tools in various studies of cucumber such as the use of GBS and single nucleotide polymorphism (SNP) markers, GBS and GWAS, also with QTL and marker-assisted selection (MAS) are applied to display and detect explicit genetic architecture complex traits in crops and chromosome rearrangements.Cucumber breeding programs have undoubtedly benefited from genotyping-by-sequencing. Using the GBS method, research discovered lots of new candidate genes that control various traits including spine color, fruit stalk-end color, and disease resistance in cucumber lines.

Genetic dissection of heterotic loci associated with plant weight by Graded pool-seq in heading Chinese cabbage (Brassica rapa)

Four heterotic QTL and a heterozygous segment for plant weight were identified by Graded Pool-Seq, QTL-seq and traditional genetic linkage analysis in heading Chinese cabbage. Heading Chinese cabbage (Brassica rapa L. spp. pekinensis) is a cross-pollinated leafy vegetable with significant heterosis. The use of heterosis is important for breeding high-yield Chinese cabbage hybrids. However, the formation and mechanism of heterosis have not been studied. We dissected the molecular mechanism of heterosis of yield-related traits in Chinese cabbage. An F₁ hybrid with high-parent heterosis of yield-related traits was selected and self-pollinated to generate segregating F₂ populations. QTL-seq, Graded Pool-seq (GPS), and traditional genetic linkage analysis were used to identify four heterotic quantitative trait loci (QTL) for plant weight: qPW1.1, qPW5.1, qPW7.1, and qPW8.1. Traditional genetic linkage analysis over two years showed that qPW8.1, located in marker A08_S45 (18,172,719) and A08_S85 (18,196,752), was mapped to a 23.5 kb genomic region. QTL qPW8.1 explained 8.6% and 23.6% of the phenotypic variation in plant weight and the total numbers of head leaves, respectively, and contained a heterozygous segment that might control the heterosis of plant weight. The qPW1.1 made an 11.7% phenotypic contribution to plant weight. The qPW7.1 was sensitive to environmental influence and explained 10.7% of the phenotypic variance. QTL qPW5.1 had a significant signal and was located in a genetic region near the centromere showing high heterozygosity. The "pseudo-overdominance" and "synergistic allelic" effects from parent line "XJD4" appear to play an important role in heterosis for plant weight in Chinese cabbage. These results provide a basis for an improved understanding of the molecular mechanism of yield-related traits and their heterosis.

SNP Discovery Using BSR-Seq Approach for Spot Blotch Resistance in Wheat (Triticum aestivum L.), an Essential Crop for Food Security

The pathogenic fungus, Bipolaris sorokiniana, that causes spot blotch (SB) disease of wheat, is a major production constraint in the Eastern Gangetic Plains of South Asia and other warm, humid regions of the world. A recombinant inbred line population was developed and phenotyped at three SB-prone locations in India. The single nucleotide polymorphism (SNP) for SB resistance was identified using a bulked segregant RNA-Seq-based approach, referred to as “BSR-Seq.” Transcriptome sequencing of the resistant parent (YS#24), the susceptible parent (YS#58), and their resistant and susceptible bulks yielded a total of 429.67 million raw reads. The bulk frequency ratio (BFR) of SNPs between the resistant and susceptible bulks was estimated, and selection of SNPs linked to resistance was done using sixfold enrichments in the corresponding bulks (BFR >6). With additional filtering criteria, the number of transcripts was further reduced to 506 with 1055 putative polymorphic SNPs distributed on 21 chromosomes of wheat. Based on SNP enrichment on chromosomal loci, five transcripts were found to be associated with SB resistance. Among the five SB resistance-associated transcripts, four were distributed on the 5B chromosome with putative 52 SNPs, whereas one transcript with eight SNPs was present on chromosome 3B. The SNPs linked to the trait were exposed to a tetra-primer ARMS-PCR assay, and an SNP-based allele-specific marker was identified for SB resistance. The in silico study of these five transcripts showed homology with pathogenesis-related genes; the metabolic pathway also exhibits similar results, suggesting their role in the plant defense mechanism.

Exome Sequencing from Bulked Segregant Analysis Identifies a Gene for All-Stage Resistance to Stripe Rust on Chromosome 1AL in Chinese Wheat Landrace 'Xiaohemai'

Stripe rust caused by Puccinia striiformis f. sp. tritici is one of the most destructive diseases of wheat. Identifying novel resistance genes applicable for developing disease-resistant cultivars is important for the sustainable control of wheat stripe rust. Chinese wheat landrace 'Xiaohemai' ('XHM') is an elite germplasm line with all-stage resistance (ASR) effective against predominant Chinese P. striiformis f. sp. tritici races. In this study, we performed a bulked segregant analysis coupled with exome capture sequencing (BSE-seq) to identify a candidate genomic region strongly associated with stripe rust resistance on chromosome 1AL in 173 F_2:3 lines derived from the cross 'XHM' × 'Avocet S'. The gene, designated as YrXH-1AL, was validated by a conventional quantitative trait locus analysis using newly developed Kompetitive allele-specific PCR (KASP) markers, explaining up to 48.50% of the phenotypic variance. By testing a secondary mapping population comprising 144 lines from the same cross at the seedling stage with prevalent P. striiformis f. sp. tritici race CYR34, YrXH-1AL was identified as a single Mendelian factor in a 1.5-cM interval flanked by KASP markers KP1A_484.33 and KP1A_490.09. This region corresponded to a 5.76-Mb genomic interval on 'Chinese Spring' chromosome 1AL. Furthermore, two cosegregating KASP markers showed high polymorphisms among 130 Chinese wheat cultivars and could be used for marker-assisted selection. Because no other Yr genes for ASR that originated from common wheat have been detected on chromosome 1AL, YrXH-1AL is likely a novel gene that can be incorporated into modern breeding materials to develop wheat cultivars with enhanced stripe rust resistance.

An autopolyploid-suitable polyBSA-seq strategy for screening candidate genetic markers linked to leaf blight resistance in sugarcane

An autopolyploid-suitable polyBSA-seq strategy was developed for screening candidate genetic markers linked to leaf blight resistance in sugarcane. Due to the complex genome architecture, the quantitative trait loci mappings and linkage marker selections for agronomic traits of autopolyploid crops were mainly limited to the time-consuming and cost intensive construction of genetic maps. To map resistance-linked markers for sugarcane leaf blight (SLB) caused by Stagonospora tainanensis, the autopolyploid-suitable bulk-segregant analysis based on the sequencing (polyBSA-seq) strategy was successfully applied for the first time. Resistant- and susceptible-bulks (R- and S-bulks) constructed from the extreme-phenotypic sugarcane F₁ lines of YT93-159 × ROC22 were deep sequenced with 195.0 × for bulks and 74.4 × for parents. Informative single-dose variants (ISDVs) present as one copy in one parent and null in the other parent were detected based on the genome sequence of LA Purple, an autooctoploid Saccharum officinarum, to screen candidate linkage markers (CLMs). The proportion of the number of short reads harboring ISDVs in the total short reads covering a given genomic position was defined as ISDV index and the ISDVs with indices met the threshold set in this study (0.04-0.14) were selected as CLMs. In total, three resistance- and one susceptibility-related CLMs for SLB resistance were identified by the polyBSA-seq. Among them, two markers on chromosome 10 were less than 300 Kb apart. Furthermore, the RNA-seq was used to calculate the expression level of genes within 1.0 Mb from the aforementioned four CLMs, which demonstrated that twelve genes were differentially expressed between resistant and susceptible clones, including a receptor-like kinase and an ethylene-responsive transcription factor. This is the first reported polyBSA-seq in autopolyploid sugarcane, which specifically tailored for the fast selection of the CLMs and causal genes associated with important agronomic traits.

Broadening the horizon of crop research: a decade of advancements in plant molecular genetics to divulge phenotype governing genes

Advancements in sequencing, genotyping, and computational technologies during the last decade (2011-2020) enabled new forward-genetic approaches, which subdue the impediments of precise gene mapping in varied crops. The modern crop improvement programs rely heavily on two major steps-trait-associated QTL/gene/marker's identification and molecular breeding. Thus, it is vital for basic and translational crop research to identify genomic regions that govern the phenotype of interest. Until the advent of next-generation sequencing, the forward-genetic techniques were laborious and time-consuming. Over the last 10 years, advancements in the area of genome assembly, genotyping, large-scale data analysis, and statistical algorithms have led faster identification of genomic variations regulating the complex agronomic traits and pathogen resistance. In this review, we describe the latest developments in genome sequencing and genotyping along with a comprehensive evaluation of the last 10-year headways in forward-genetic techniques that have shifted the focus of plant research from model plants to diverse crops. We have classified the available molecular genetic methods under bulk-segregant analysis-based (QTL-seq, GradedPool-Seq, QTG-Seq, Exome QTL-seq, and RapMap), target sequence enrichment-based (RenSeq, AgRenSeq, and TACCA), and mutation-based groups (MutMap, NIKS algorithm, MutRenSeq, MutChromSeq), alongside improvements in classical mapping and genome-wide association analyses. Newer methods for outcrossing, heterozygous, and polyploid plant genetics have also been discussed. The use of k-mers has enriched the nature of genetic variants which can be utilized to identify the phenotype-causing genes, independent of reference genomes. We envisage that the recent methods discussed herein will expand the repertoire of useful alleles and help in developing high-yielding and climate-resilient crops.