Genetic Dissection of Seedling Root System Architectural Traits in a Diverse Panel of Hexaploid Wheat through Multi-Locus Genome-Wide Association Mapping for Improving Drought Tolerance

Unraveling the genetics of heat tolerance in chickpea landraces (Cicer arietinum L.) using genome-wide association studies

Chickpea, being an important grain legume crop, is often confronted with the adverse effects of high temperatures at the reproductive stage of crop growth, drastically affecting yield and overall productivity. The current study deals with an extensive evaluation of chickpea genotypes, focusing on the traits associated with yield and their response to heat stress. Notably, we observed significant variations for these traits under both normal and high-temperature conditions, forming a robust basis for genetic research and breeding initiatives. Furthermore, the study revealed that yield-related traits exhibited high heritability, suggesting their potential suitability for marker-assisted selection. We carried out single-nucleotide polymorphism (SNP) genotyping using the genotyping-by-sequencing (GBS) method for a genome-wide association study (GWAS). Overall, 27 marker-trait associations (MTAs) linked to yield-related traits, among which we identified five common MTAs displaying pleiotropic effects after applying a stringent Bonferroni-corrected p-value threshold of <0.05 [-log10(p) > 4.95] using the BLINK (Bayesian-information and linkage-disequilibrium iteratively nested keyway) model. Through an in-depth in silico analysis of these markers against the CDC Frontier v1 reference genome, we discovered that the majority of the SNPs were located at or in proximity to gene-coding regions. We further explored candidate genes situated near these MTAs, shedding light on the molecular mechanisms governing heat stress tolerance and yield enhancement in chickpeas such as indole-3-acetic acid-amido synthetase GH3.1 with GH3 auxin-responsive promoter and pentatricopeptide repeat-containing protein, etc. The harvest index (HI) trait was associated with marker Ca3:37444451 encoding aspartic proteinase ortholog sequence of Oryza sativa subsp. japonica and Medicago truncatula, which is known for contributing to heat stress tolerance. These identified MTAs and associated candidate genes may serve as valuable assets for breeding programs dedicated to tailoring chickpea varieties resilient to heat stress and climate change.

Genome-wide association study of drought tolerance in wheat (Triticum aestivum L.) identifies SNP markers and candidate genes

Drought stress poses a severe threat to global wheat production, necessitating an in-depth exploration of the genetic basis for drought tolerance associated traits. This study employed a 90 K SNP array to conduct a genome-wide association analysis, unravelling genetic determinants of key traits related to drought tolerance in wheat, namely plant height, root length, and root and shoot dry weight. Using the mixed linear model (MLM) method on 125 wheat accessions subjected to both well-watered and drought stress treatments, we identified 53 SNPs significantly associated with stress susceptibility (SSI) and tolerance indices (STI) for the targeted traits. Notably, chromosomes 2A and 3B stood out with ten and nine associated markers, respectively. Across 17 chromosomes, 44 unique candidate genes were pinpointed, predominantly located on the distal ends of 1A, 1B, 1D, 2A, 3A, 3B, 4A, 6A, 6B, 7A, 7B, and 7D chromosomes. These genes, implicated in diverse functions related to plant growth, development, and stress responses, offer a rich resource for future investigation. A clustering pattern emerged, notably with seven genes associated with SSI for plant height and four genes linked to both STI of plant height and shoot dry weight, converging on specific regions of chromosome arms of 2AS and 3BL. Additionally, shared genes encoding polygalacturonase, auxilin-related protein 1, peptide deformylase, and receptor-like kinase underscored the interconnectedness between plant height and shoot dry weight. In conclusion, our findings provide insights into the molecular mechanisms governing wheat drought tolerance, identifying promising genomic loci for further exploration and crop improvement strategies.

A Multi-Model Based Stability Analysis Employing Multi-Environmental Trials (METs) Data for Discerning Heat Tolerance in Chickpea (Cicer arietinum L.) Landraces

Identifying a congenially targeted production environment and understanding the effects of genotype by environmental interactions on the adaption of chickpea genotypes is essential for achieving an optimal yield stability. Different models like additive main effect and multiplicative interactions (AMMI 1, AMM2), weighted average absolute scores of BLUPs (WAASB), and genotype plus genotype-environment (GGE) interactions were used to understand their suitability in the precise estimation of variance and their interaction. Our experiment used genotypes that represent the West Asia-North Africa (WANA) region. This trial involved two different sowing dates, two distinct seasons, and three different locations, resulting in a total of 12 environments. Genotype IG 5871(G1) showed a lower heat susceptibility index (HSI) across environments under study. The first four interactions principal component axis (IPCA) explain 93.2% of variations with significant genotype-environment interactions. Considering the AMMI stability value (ASV), the genotypes IG5862(G7), IG5861(G6), ILC239(G40), IG6002(G26), and ILC1932(G39), showing ASV scores of 1.66, 1.80, 2.20, 2.60, and 2.84, respectively, were ranked as the most stable and are comparable to the weighted average absolute scores of BLUPs (WAASB) ranking of genotypes. The which-won-where pattern of genotype plus genotype-environment (GGE) interactions suggested that the target environment consists of one mega environment. IG5866(G10), IG5865(G9), IG5884(G14), and IG5862(G7) displayed higher stability, as they were nearer to the origin. The genotypes that exhibited a superior performance in the tested environments can serve as ideal parental lines for heat-stress tolerance breeding programs. The weighted average absolute scores of BLUPs (WAASB) serve as an ideal tool to discern the variations and identify the stable genotype among all methods.

Unravelling the Genetic Basis of Moisture Deficit Stress Tolerance in Wheat for Seedling Vigour-Related Traits and Root Traits Using Genome-Wide Association Study

A key abiotic stress that negatively affects seed germination, plant development, and crop yield is moisture deficit stress. Achieving higher vigour and uniform germination under stress conditions is essential for crop establishment and productivity and to enhance the yield. Hence, revealing wheat's capacity to withstand moisture deficit stress during seed germination and early growth stages is fundamental in improving its overall performance. However, the genetic regulation of moisture deficit stress tolerance during the seed germination phase remains largely unexplored. In this study, a total of 193 wheat genotypes were subjected to simulated moisture deficit stress using PEG-6000 (-0.4 MPa) during the seed germination stage. The induced moisture deficit stress significantly reduced various seedling-vigour-related traits. The genetic regions linked to these traits were found using a genome-wide association study (GWAS). The analysis identified 235 MTAs with a significance -log10(p) value of >4. After applying the Bonferroni correction, the study identified 47 unique single nucleotide polymorphisms (SNPs) that are linked to candidate genes important for the trait of interest. The current study emphasises the effectiveness of genome-wide association studies (GWAS) in identifying promising candidate genes, improving wheat seedling vigour and root traits, and offering essential information for the development of wheat cultivars tolerant to moisture deficit stress.

Evaluation of millets for physio-chemical and root morphological traits suitable for resilient farming and nutritional security in Eastern Himalayas

Introduction

Millets are nutritionally superior and climate-resilient short-duration crops and hold a prominent place in cropping sequences around the world. They have immense potential to grow in a marginal environment due to diverse adaptive mechanisms.

Methods

An experiment was conducted in an organic production system in the North Eastern Himalayan foothills of India for 3 consecutive years by evaluating high-yielding varieties (HYVs) of different millets, viz., finger millet, foxtail millet, little millet, barnyard millet, proso millet, and browntop millet, along with local landraces of finger millets (Sikkim-1 and Sikkim-2; Nagaland-1 and Nagaland-2) to identify stable, high-yielding, and nutritionally superior genotypes suited for the region.

Results

Among the various millets, finger millet, followed by little millet and foxtail millet, proved their superiority in terms of productivity (ranging between 1.16 and 1.43 Mg ha^-1) compared to other millets. Among different varieties of finger millets, cv. VL Mandua 352 recorded the highest average grain yield (1.43 Mg ha^-1) followed by local landraces, Nagaland-2 (1.31 Mg ha^-1) and Sikkim-1 (1.25 Mg ha^-1). Root traits such as total root length, root volume, average diameter of roots, and root surface area were significantly higher in finger millet landraces Nagaland-1, Nagaland-2, and Sikkim-1 compared to the rest of the millet genotypes. The different millets were found to be rich sources of protein as recorded in foxtail millet cv. SiA 3088 (12.3%), proso millet cv. TNAU 145 (11.5%), and finger millet landraces, Sikkim-1 and Nagaland-2 (8.7% each). Finger millet landrace Sikkim-2 recorded the highest omega-6 content (1.16%), followed by barnyard millet cv. VL 207 (1.09%). Barnyard millet cv. VL 207 recorded the highest polyunsaturated fatty acid (PUFA) content (1.23%), followed by foxtail millet cv. SiA 3088 (1.09%). The local finger millet landraces Sikkim-1 and Sikkim-2 recorded the highest levels of histidine (0.41%) and tryptophan (0.12%), respectively. Sikkim-1 and Nagaland-2 recorded the highest level of thiamine (0.32%) compared to the HYVs.

Conclusion

These findings indicate that finger millet has great potential in the organic production system of the North Eastern Himalayan Region (NEHR) of India, and apart from HYVs like VL Mandua 352, local landraces, viz., Nagaland-2 and Sikkim-1, should also be promoted for ensuring food and nutritional security in this fragile ecosystem.

Genome-wide association mapping of genomic regions associated with drought stress tolerance at seedling and reproductive stages in bread wheat

Understanding the genetic architecture of drought stress tolerance in bread wheat at seedling and reproductive stages is crucial for developing drought-tolerant varieties. In the present study, 192 diverse wheat genotypes, a subset from the Wheat Associated Mapping Initiative (WAMI) panel, were evaluated at the seedling stage in a hydroponics system for chlorophyll content (CL), shoot length (SLT), shoot weight (SWT), root length (RLT), and root weight (RWT) under both drought and optimum conditions. Following that, a genome-wide association study (GWAS) was carried out using the phenotypic data recorded during the hydroponics experiment as well as data available from previously conducted multi-location field trials under optimal and drought stress conditions. The panel had previously been genotyped using the Infinium iSelect 90K SNP array with 26,814 polymorphic markers. Using single as well as multi-locus models, GWAS identified 94 significant marker-trait associations (MTAs) or SNPs associated with traits recorded at the seedling stage and 451 for traits recorded at the reproductive stage. The significant SNPs included several novel, significant, and promising MTAs for different traits. The average LD decay distance for the whole genome was approximately 0.48 Mbp, ranging from 0.07 Mbp (chromosome 6D) to 4.14 Mbp (chromosome 2A). Furthermore, several promising SNPs revealed significant differences among haplotypes for traits such as RLT, RWT, SLT, SWT, and GY under drought stress. Functional annotation and in silico expression analysis revealed important putative candidate genes underlying the identified stable genomic regions such as protein kinases, O-methyltransferases, GroES-like superfamily proteins, NAD-dependent dehydratases, etc. The findings of the present study may be useful for improving yield potential, and stability under drought stress conditions.

Chromosome groups 5, 6 and 7 harbor major quantitative trait loci controlling root traits in bread wheat (Triticum aestivum L.)

Identifying genomic regions for root traits in bread wheat can help breeders develop climate-resilient and high-yielding wheat varieties with desirable root traits. This study used the recombinant inbred line (RIL) population of Synthetic W7984 × Opata M85 to identify quantitative trait loci (QTL) for different root traits such as rooting depth (RD), root dry mass (RM), total root length (RL), root diameter (Rdia) and root surface areas (RSA1 for coarse roots and RSA2 for fine roots) under controlled conditions in a semi-hydroponic system. We detected 14 QTL for eight root traits on nine wheat chromosomes; we discovered three QTL each for RD and RSA1, two QTL each for RM and RSA2, and one QTL each for RL, Rdia, specific root length and nodal root number per plant. The detected QTL were concentrated on chromosome groups 5, 6 and 7. The QTL for shallow RD (Q.rd.uwa.7BL: Xbarc50) and high RM (Q.rm.uwa.6AS: Xgwm334) were validated in two independent F₂ populations of Synthetic W7984 × Chara and Opata M85 × Cascade, respectively. Genotypes containing negative alleles for Q.rd.uwa.7BL had 52% shallower RD than other Synthetic W7984 × Chara population lines. Genotypes with the positive alleles for Q.rm.uwa.6AS had 31.58% higher RM than other Opata M85 × Cascade population lines. Further, we identified 21 putative candidate genes for RD (Q.rd.uwa.7BL) and 13 for RM (Q.rm.uwa.6AS); TraesCS6A01G020400, TraesCS6A01G024400 and TraesCS6A01G021000 identified from Q.rm.uwa.6AS, and TraesCS7B01G404000, TraesCS7B01G254900 and TraesCS7B01G446200 identified from Q.rd.uwa.7BL encoded important proteins for root traits. We found germin-like protein encoding genes in both Q.rd.uwa.7BL and Q.rm.uwa.6AS regions. These genes may play an important role in RM and RD improvement. The identified QTL, especially the validated QTL and putative candidate genes are valuable genetic resources for future root trait improvement in wheat.

Genome-wide association study of agronomical and root-related traits in spring barley collection grown under field conditions

The root system is a key component for plant survival and productivity. In particular, under stress conditions, developing plants with a better root architecture can ensure productivity. The objectives of this study were to investigate the phenotypic variation of selected root- and yield-related traits in a diverse panel of spring barley genotypes. By performing a genome-wide association study (GWAS), we identified several associations underlying the variations occurring in root- and yield-related traits in response to natural variations in soil moisture. Here, we report the results of the GWAS based on both individual single-nucleotide polymorphism markers and linkage disequilibrium (LD) blocks of markers for 11 phenotypic traits related to plant morphology, grain quality, and root system in a group of spring barley accessions grown under field conditions. We also evaluated the root structure of these accessions by using a nondestructive method based on electrical capacitance. The results showed the importance of two LD-based blocks on chromosomes 2H and 7H in the expression of root architecture and yield-related traits. Our results revealed the importance of the region on the short arm of chromosome 2H in the expression of root- and yield-related traits. This study emphasized the pleiotropic effect of this region with respect to heading time and other important agronomic traits, including root architecture. Furthermore, this investigation provides new insights into the roles played by root traits in the yield performance of barley plants grown under natural conditions with daily variations in soil moisture content.

Genetic dissection of marker trait associations for grain micro-nutrients and thousand grain weight under heat and drought stress conditions in wheat

Introduction

Wheat is grown and consumed worldwide, making it an important staple food crop for both its calorific and nutritional content. In places where wheat is used as a staple food, suboptimal micronutrient content levels, especially of grain iron (Fe) and zinc (Zn), can lead to malnutrition. Grain nutrient content is influenced by abiotic stresses, such as drought and heat stress. The best method for addressing micronutrient deficiencies is the biofortification of food crops. The prerequisites for marker-assisted varietal development are the identification of the genomic region responsible for high grain iron and zinc contents and an understanding of their genetics.

Methods

A total of 193 diverse wheat genotypes were evaluated under drought and heat stress conditions across the years at the Indian Agricultural Research Institute (IARI), New Delhi, under timely sown irrigated (IR), restricted irrigated (RI) and late sown (LS) conditions. Grain iron content (GFeC) and grain zinc content (GZnC) were estimated from both the control and treatment groups. Genotyping of all the lines under study was carried out with the single nucleotide polymorphisms (SNPs) from Breeder's 35K Axiom Array.

Result and Discussion

Three subgroups were observed in the association panel based on both principal component analysis (PCA) and dendrogram analysis. A large whole-genome linkage disequilibrium (LD) block size of 3.49 Mb was observed. A genome-wide association study identified 16 unique stringent marker trait associations for GFeC, GZnC, and 1000-grain weight (TGW). In silico analysis demonstrated the presence of 28 potential candidate genes in the flanking region of 16 linked SNPs, such as synaptotagmin-like mitochondrial-lipid-binding domain, HAUS augmin-like complex, di-copper center-containing domain, protein kinase, chaperonin Cpn60, zinc finger, NUDIX hydrolase, etc. Expression levels of these genes in vegetative tissues and grain were also found. Utilization of identified markers in marker-assisted breeding may lead to the rapid development of biofortified wheat genotypes to combat malnutrition.

Quantitative trait loci mapping reveals important genomic regions controlling root architecture and shoot biomass under nitrogen, phosphorus, and potassium stress in rapeseed (Brassica napus L.)

Plants rely on root systems for nutrient uptake from soils. Marker-assisted selection helps breeders to select desirable root traits for effective nutrient uptake. Here, 12 root and biomass traits were investigated at the seedling stage under low nitrogen (LN), low phosphorus (LP), and low potassium (LK) conditions, respectively, in a recombinant inbred line (RIL) population, which was generated from Brassica napus L. Zhongshuang11 and 4D122 with significant differences in root traits and nutrient efficiency. Significant differences for all the investigated traits were observed among RILs, with high heritabilities (0.43-0.74) and high correlations between the different treatments. Quantitative trait loci (QTL) mapping identified 57, 27, and 36 loci, explaining 4.1-10.9, 4.6-10.8, and 4.9-17.4% phenotypic variances under LN, LP, and LK, respectively. Through QTL-meta analysis, these loci were integrated into 18 significant QTL clusters. Four major QTL clusters involved 25 QTLs that could be repeatedly detected and explained more than 10% phenotypic variances, including two NPK-common and two specific QTL clusters (K and NK-specific), indicating their critical role in cooperative nutrients uptake of N, P, and K. Moreover, 264 genes within the four major QTL clusters having high expressions in roots and SNP/InDel variations between two parents were identified as potential candidate genes. Thirty-eight of them have been reported to be associated with root growth and development and/or nutrient stress tolerance. These key loci and candidate genes lay the foundation for deeper dissection of the NPK starvation response mechanisms in B. napus.

Genome-wide association mapping for component traits of drought and heat tolerance in wheat

Identification of marker trait association is a prerequisite for marker-assisted breeding. To find markers linked with traits under heat and drought stress in bread wheat (Triticum aestivum L.), we performed a genome-wide association study (GWAS). GWAS mapping panel used in this study consists of advanced breeding lines from the IARI stress breeding programme produced by pairwise and complex crosses. Phenotyping was done at multi locations namely New Delhi, Karnal, Indore, Jharkhand and Pune with augmented-RCBD design under different moisture and heat stress regimes, namely timely sown irrigated (IR), timely sown restricted irrigated (RI) and late sown (LS) conditions. Yield and its component traits, viz., Days to Heading (DH), Days to Maturity (DM), Normalized Difference Vegetation Index (NDVI), Chlorophyll Content (SPAD), Canopy temperature (CT), Plant Height (PH), Thousand grain weight (TGW), Grain weight per spike (GWPS), Plot Yield (PLTY) and Biomass (BMS) were phenotyped. Analysis of variance and descriptive statistics revealed significant differences among the studied traits. Genotyping was done using the 35k SNP Wheat Breeder's Genotyping Array. Population structure and diversity analysis using filtered 10,546 markers revealed two subpopulations with sufficient diversity. A large whole genome LD block size of 7.15 MB was obtained at half LD decay value. Genome-wide association search identified 57 unique markers associated with various traits across the locations. Twenty-three markers were identified to be stable, among them nine pleiotropic markers were also identified. In silico search of the identified markers against the IWGSC ref genome revealed the presence of a majority of the SNPs at or near the gene coding region. These SNPs can be used for marker-assisted transfer of genes/QTLs after validation to develop climate-resilient cultivars.

Genome-Wide Identification and Characterization of Receptor-Like Protein Kinase 1 (RPK1) Gene Family in Triticum aestivum Under Drought Stress

Receptor-like protein kinase1 (RPK1) genes play crucial roles in plant growth and development processes, root architecture, and abiotic stress regulation. A comprehensive study of the RPK1 gene family has not been reported in bread wheat (Triticum aestivum). Here, we reported the genome-wide identification, characterization, and expression patterns of the RPK1 gene family in wheat. Results confirmed 15 TaRPK1 genes, classified mainly into three sub-clades based on a phylogenetic tree. The TaRPK1 genes were mapped on chromosomes 1–3 in the respective A, B, and D genomes. Gene structure, motif conservation, collinearity prediction, and synteny analysis were carried out systematically. A Gene ontology study revealed that TaRPK1 genes play a vital role during molecular and biological processes. We also identified 18 putative miRNAs targeting TaRPK1 genes, suggesting their roles in growth, development, and stress responses. Cis-Regulatory elements interpreted the presence of light-related elements, hormone responsiveness, and abiotic stress-related motifs in the promoter regions. The SWISS_MODEL predicted the successful models of TaRPK1 proteins with at least 30% identity to the template, a widely accepted threshold for successful modeling. In silico expression analysis in different tissues and stages suggested that TaRPK1 genes exhibited the highest expression in root tissues. Moreover, qRT-PCR further validated the higher expression of TaRPK1 genes in roots of drought-tolerant varieties compared to the drought-susceptible variety. Collectively, the present study renders valuable information on the functioning of TaRPK1 genes in wheat that will be useful in further functional validation of these genes in future studies.

Multi-locus genome-wide association studies (ML-GWAS) reveal novel genomic regions associated with seedling and adult plant stage leaf rust resistance in bread wheat (Triticum aestivum L.)

Leaf rust is one of the important diseases limiting global wheat production and productivity. To identify quantitative trait nucleotides (QTNs) or genomic regions associated with seedling and adult plant leaf rust resistance, multilocus genome-wide association studies (ML-GWAS) were performed on a panel of 400 diverse wheat genotypes using 35 K single-nucleotide polymorphism (SNP) genotyping assays and trait data of leaf rust resistance. Association analyses using six multi-locus GWAS models revealed a set of 201 significantly associated QTNs for seedling and 65 QTNs for adult plant resistance (APR), explaining 1.98-31.72% of the phenotypic variation for leaf rust. Among these QTNs, 51 reliable QTNs for seedling and 15 QTNs for APR were consistently detected in at least two GWAS models and were considered reliable QTNs. Three genomic regions were pleiotropic, each controlling two to three pathotype-specific seedling resistances to leaf rust. We also identified candidate genes, such as leucine-rich repeat receptor-like (LRR) protein kinases, P-loop containing nucleoside triphosphate hydrolase and serine-threonine/tyrosine-protein kinases (STPK), which have a role in pathogen recognition and disease resistance linked to the significantly associated genomic regions. The QTNs identified in this study can prove useful in wheat molecular breeding programs aimed at enhancing resistance to leaf rust and developing next-generation leaf rust-resistant varieties.

The Adaptation and Tolerance of Major Cereals and Legumes to Important Abiotic Stresses

Abiotic stresses, including drought, extreme temperatures, salinity, and waterlogging, are the major constraints in crop production. These abiotic stresses are likely to be amplified by climate change with varying temporal and spatial dimensions across the globe. The knowledge about the effects of abiotic stressors on major cereal and legume crops is essential for effective management in unfavorable agro-ecologies. These crops are critical components of cropping systems and the daily diets of millions across the globe. Major cereals like rice, wheat, and maize are highly vulnerable to abiotic stresses, while many grain legumes are grown in abiotic stress-prone areas. Despite extensive investigations, abiotic stress tolerance in crop plants is not fully understood. Current insights into the abiotic stress responses of plants have shown the potential to improve crop tolerance to abiotic stresses. Studies aimed at stress tolerance mechanisms have resulted in the elucidation of traits associated with tolerance in plants, in addition to the molecular control of stress-responsive genes. Some of these studies have paved the way for new opportunities to address the molecular basis of stress responses in plants and identify novel traits and associated genes for the genetic improvement of crop plants. The present review examines the responses of crops under abiotic stresses in terms of changes in morphology, physiology, and biochemistry, focusing on major cereals and legume crops. It also explores emerging opportunities to accelerate our efforts to identify desired traits and genes associated with stress tolerance.

Identification of Novel Quantitative Trait Nucleotides and Candidate Genes for Bacterial Wilt Resistance in Tobacco (Nicotiana tabacum L.) Using Genotyping-by-Sequencing and Multi-Locus Genome-Wide Association Studies

Tobacco bacterial wilt (TBW) is a devastating soil-borne disease threatening the yield and quality of tobacco. However, its genetic foundations are not fully understood. In this study, we identified 126,602 high-quality single-nucleotide polymorphisms (SNPs) in 94 tobacco accessions using genotyping-by-sequencing (GBS) and a 94.56 KB linkage disequilibrium (LD) decay rate for candidate gene selection. The population structure analysis revealed two subpopulations with 37 and 57 tobacco accessions. Four multi-locus genome-wide association study (ML-GWAS) approaches identified 142 quantitative trait nucleotides (QTNs) in E1-E4 and the best linear unbiased prediction (BLUP), explaining 0.49-22.52% phenotypic variance. Of these, 38 novel stable QTNs were identified across at least two environments/methods, and their alleles showed significant TBW-DI differences. The number of superior alleles associated with TBW resistance for each accession ranged from 4 to 24; eight accessions had more than 18 superior alleles. Based on TBW-resistant alleles, the five best cross combinations were predicted, including MC133 × Ruyuan No. 1 and CO258 × ROX28. We identified 52 candidate genes around 38 QTNs related to TBW resistance based on homologous functional annotation and KEGG enrichment analysis, e.g., CYCD3;2, BSK1, Nitab4.5_0000641g0050, Nitab4.5_0000929g0030. To the best of our knowledge, this is the first comprehensive study to identify QTNs, superior alleles, and their candidate genes for breeding TBW-resistant tobacco varieties. The results provide further insight into the genetic architecture, marker-assisted selection, and functional genomics of TBW resistance, improving future breeding efforts to increase crop productivity.