Bulked sample analysis in genetics, genomics and crop improvement

Determination of lethal and mutation induction doses of gamma rays for gladiolus (Gladiolus grandifloras Hort.) genotypes

Gladiolus is a highly allogamous flower plant, but owing to the prolonged juvenile phase, asexual propagation is preferred, which acts as a barrier for the induction of natural genetic variability in gladiolus. Therefore, the induced mutagenesis could be utilized for the creation of desirable genotypes, without altering their basic agronomic features. An analysis of the optimum doses of γ radiation for the induction of fruitful mutations could be achieved in short period of time, compared with the conventional method of breeding. The objectives of this study were to perform radiosensitivity tests on various gladiolus genotypes using different doses of gamma rays and to determine the optimal dose of radiation dose for obtaining the greatest number of mutants. The present experiment was carried out during the winter-spring seasons, for the four consecutive years of 2017-18, 2018-19, 2019-20, and 2020-21. The seven genotypes of gladiolus were exposed to seven doses of gamma rays (60Cobalt). Plants irradiated with radiation doses lower than 4.5 Kr (G1) had greater plant survivability than the higher doses of gamma rays (≥5.0 Kr). The radiation of G0 (0 Kr) result in highest plant survivability, while radiation dose of G6 (6.5 Kr) resulted lowest survivability. LD25 and BD50 for all the genotypes were achieved except for V5 and V7, similarly the median lethal doses (LD50) for V3 and V4 genotypes had been achieved. The highest flower blindness percent and percent abnormal plants were observed at G5 and G6 and between the 4.0 Kr (G1) and 5.5 Kr (G4) gamma ray doses, respectively. The flower colour mutation frequency was recorded highest in genotypes Tiger Flame at 5.0 Kr (V7G3), while the Flower colour mutation spectrum was identified between 4.0 Kr (G1) to 5.5 Kr (G4) in all the genotypes except for genotypes V5 and V7. For the generation of higher phenotypic variations, radiation dose between 4.0 Kr (G1) and 5.5 Kr (G4) were found the most prominent. Specifically the gamma rays radiation dose of 5.5 Kr (G4) resulted in the highest flower colour mutation frequency. These isolated mutant lines will broaden the gladiolus gene pool and support future gladiolus breeding experiments.

Identification of Dwarfing Candidate Genes in Brassica napus L. LSW2018 through BSA-Seq and Genetic Mapping

Plant height, as a crucial component of plant architecture, exerts a significant influence on rapeseed (Brassica napus L.) lodging resistance, photosynthetic efficiency, yield, and mechanized harvest level. A previous study identified dwarf rapeseed LSW2018. In this study, LSW2018 (dwarf parent (PD)) was crossed with 389 (high parent (PH)) to establish the F2 population, and 30 extremely dwarf (bulk-D) and high (bulk-H) plants in the F2 population were respectively selected to construct two bulked DNA pools. Whole-genome sequencing and variation analysis (BSA-seq) were performed on these four DNA pools (PD, PH, bulk-D, and bulk-H). The BSA-seq results revealed that the genomic region responsible for the dwarf trait spanned from 19.30 to 22.19 Mb on chromosome A03, with a length of 2.89 Mb. After fine mapping with simple sequence repeat (SSR) markers, the gene was narrowed to a 0.71 Mb interval. Within this region, a total of 113 genes were identified, 42 of which contained large-effect variants. According to reference genome annotation and qRT-PCR analysis, there are 17 differentially expressed genes in this region between high and dwarf individuals. This study preliminarily reveals the genetic basis of LSW2018 dwarfing and provides a theoretical foundation for the molecular marker-assisted breeding of dwarf rapeseed.

Inferring the composition of a mixed culture of natural microbial isolates by deep sequencing

Next generation sequencing has unlocked a wealth of genotype information for microbial populations, but phenotyping remains a bottleneck for exploiting this information, particularly for pathogens that are difficult to manipulate. Here, we establish a method for high-throughput phenotyping of mixed cultures, in which the pattern of naturally occurring single-nucleotide polymorphisms in each isolate is used as intrinsic barcodes which can be read out by sequencing. We demonstrate that our method can correctly deconvolute strain proportions in simulated mixed-strain pools. As an experimental test of our method, we perform whole genome sequencing of 66 natural isolates of the thermally dimorphic pathogenic fungus Coccidioides posadasii and infer the strain compositions for large mixed pools of these strains after competition at 37°C and room temperature. We validate the results of these selection experiments by recapitulating the temperature-specific enrichment results in smaller pools. Additionally, we demonstrate that strain fitness estimated by our method can be used as a quantitative trait for genome-wide association studies. We anticipate that our method will be broadly applicable to natural populations of microbes and allow high-throughput phenotyping to match the rate of genomic data acquisition.

A complex receptor locus confers responsiveness to necrosis and ethylene-inducing like peptides in Brassica napus

Brassica crops are susceptible to diseases which can be mitigated by breeding for resistance. MAMPs (microbe-associated molecular patterns) are conserved molecules of pathogens that elicit host defences known as pattern-triggered immunity (PTI). Necrosis and Ethylene-inducing peptide 1-like proteins (NLPs) are MAMPs found in a wide range of phytopathogens. We studied the response to BcNEP2, a representative NLP from Botrytis cinerea, and showed that it contributes to disease resistance in Brassica napus. To map regions conferring NLP response, we used the production of reactive oxygen species (ROS) induced during PTI across a population of diverse B. napus accessions for associative transcriptomics (AT), and bulk segregant analysis (BSA) on DNA pools created from a cross of NLP-responsive and non-responsive lines. In silico mapping with AT identified two peaks for NLP responsiveness on chromosomes A04 and C05 whereas the BSA identified one peak on A04. BSA delimited the region for NLP-responsiveness to 3 Mbp, containing ~245 genes on the Darmor-bzh reference genome and four co-segregating KASP markers were identified. The same pipeline with the ZS11 genome confirmed the highest-associated region on chromosome A04. Comparative BLAST analysis revealed unannotated clusters of receptor-like protein (RLP) homologues on ZS11 chromosome A04. However, no specific RLP homologue conferring NLP response could be identified. Our results also suggest that BR-SIGNALLING KINASE1 may be involved with modulating the NLP response. Overall, we demonstrate that responsiveness to NLP contributes to disease resistance in B. napus and define the associated genomic location. These results can have practical application in crop improvement.

Construction of a high-density genetic map and QTL localization of body weight and wool production related traits in Alpine Merino sheep based on WGR

BACKGROUND

The Alpine Merino is a new breed of fine-wool sheep adapted to the cold and arid climate of the plateau in the world. It has been popularized in Northwest China due to its superior adaptability as well as excellent production performance. Those traits related to body weight, wool yield, and wool fiber characteristics, which are economically essential traits in Alpine Merino sheep, are controlled by QTL (Quantitative Trait Loci). Therefore, the identification of QTL and genetic markers for these key economic traits is a critical step in establishing a MAS (Marker-Assisted Selection) breeding program.

RESULTS

In this study, we constructed the high-density genetic linkage map of Alpine Merino sheep by sequencing 110 F1 generation individuals using WGR (Whole Genome Resequencing) technology. 14,942 SNPs (Single Nucleotide Polymorphism) were identified and genotyped. The map spanned 2,697.86 cM, with an average genetic marker interval of 1.44 cM. A total of 1,871 high-quality SNP markers were distributed across 27 linkage groups, with an average of 69 markers per LG (Linkage Group). Among them, the smallest genetic distance is 19.62 cM for LG2, while the largest is 237.19 cM for LG19. The average genetic distance between markers in LGs ranged from 0.24 cM (LG2) to 3.57 cM (LG17). The marker density in the LGs ranged from LG14 (39 markers) to LG1 (150 markers).

CONCLUSIONS

The first genetic map of Alpine Merino sheep we constructed included 14,942 SNPs, while 46 QTLs associated with body weight, wool yield and wool fiber traits were identified, laying the foundation for genetic studies and molecular marker-assisted breeding. Notably, there were QTL intervals for overlapping traits on LG4 and LG8, providing potential opportunities for multi-trait co-breeding and further theoretical support for selection and breeding of ultra-fine and meaty Alpine Merino sheep.

Pooled DNA sequencing in hairy vetch (Vicia villosa Roth) reveals QTL for seed dormancy but not pod dehiscence

Introduction

Hairy vetch (Vicia villosa Roth) is a promising legume cover crop, but its use is limited by high rates of pod dehiscence and seed dormancy.

Methods

We used phenotypically contrasting pooled DNA samples (n=24 with 29-74 individuals per sample) from an ongoing cover crop breeding program across four environments (site-year combinations: Maryland 2020, Maryland 2022, Wisconsin 2021, Wisconsin 2022) to find genetic associations and genomic prediction accuracies for pod dehiscence and seed dormancy. We also combined pooled DNA sample genetic association results with the results of a prior genome-wide association study.

Results and discussion

Genomic prediction resulted in positive predictive abilities for both traits between environments and with an independent dataset (0.34-0.50), but reduced predictive ability for DNA pools with divergent seed dormancy in the Maryland environments (0.07-0.15). The pooled DNA samples found six significant (false discovery rate q-value<0.01) quantitative trait loci (QTL) for seed dormancy and four significant QTL for pod dehiscence. Unfortunately, the minor alleles of the pod dehiscence QTL increased the rate of pod dehiscence and are not useful for marker-assisted selection. When combined with a prior association study, sixteen seed dormancy QTL and zero pod dehiscence QTL were significant. Combining the association studies did not increase the detection of useful QTL.

OcBSA: An NGS-based bulk segregant analysis tool for outcross populations

Constructing inbred lines for self-incompatible species and species with long generation times is challenging, making the use of F1 outcross/segregating populations the main strategy for genetic studies of such species. However, there is a lack of dedicated algorithms/tools for rapid quantitative trait locus (QTL) mapping using the F1 populations. To this end, we have designed and developed an algorithm/tool called OcBSA specifically for QTL mapping of F1 populations. OcBSA transforms the four-haplotype inheritance problem from the two heterozygous diploid parents of the F1 population into the two-haplotype inheritance problem common in current genetic studies by removing the two haplotypes from the heterozygous parent that do not contribute to phenotype segregation in the F1 population. Testing of OcBSA on 1800 simulated F1 populations demonstrated its advantages over other currently available tools in terms of sensitivity and accuracy. In addition, the broad applicability of OcBSA was validated by QTL mapping using seven reported F1 populations of apple, pear, peach, citrus, grape, tea, and rice. We also used OcBSA to map the QTL for flower color in a newly constructed F1 population of potato generated in this study. The OcBSA mapping result was verified by the insertion or deletion markers to be consistent with a previously reported locus harboring the ANTHOCYANIN 2 gene, which regulates potato flower color. Taken together, these results highlight the power and broad utility of OcBSA for QTL mapping using F1 populations and thus a great potential for functional gene mining in outcrossing species. For ease of use, we have developed both Windows and Linux versions of OcBSA, which are freely available at: https://gitee.com/Bioinformaticslab/OcBSA.

Application of bulk segregant RNA-Seq (BSR-Seq) and allele-specific primers to study soybean powdery mildew resistance

BACKGROUND

Powdery mildew (PM) is one of the important soybean diseases, and host resistance could practically contribute to soybean PM management. To date, only the Rmd locus on chromosome (Chr) 16 was identified through traditional QTL mapping and GWAS, and it remains unclear if the bulk segregant RNA-Seq (BSR-Seq) methodology is feasible to explore additional PM resistance that might exist in other varieties.

RESULTS

BSR-Seq was applied to contrast genotypes and gene expressions between the resistant bulk (R bulk) and the susceptible bulk (S bulk), as well as the parents. The ∆(SNP-index) and G' value identified several QTL and significant SNPs/Indels on Chr06, Chr15, and Chr16. Differentially expressed genes (DEGs) located within these QTL were identified using HISAT2 and Kallisto, and allele-specific primers (AS-primers) were designed to validate the accuracy of phenotypic prediction. While the AS-primers on Chr06 or Chr15 cannot distinguish the resistant and susceptible phenotypes, AS-primers on Chr16 exhibited 82% accuracy prediction with an additive effect, similar to the SSR marker Satt431.

CONCLUSIONS

Evaluation of additional AS-primers in the linkage disequilibrium (LD) block on Chr16 further confirmed the resistant locus, derived from the resistant parental variety 'Kaohsiung 11' ('KS11'), not only overlaps with the Rmd locus with unique up-regulated LRR genes (Glyma.16G213700 and Glyma.16G215100), but also harbors a down-regulated MLO gene (Glyma.16G145600). Accordingly, this study exemplified the feasibility of BSR-Seq in studying biotrophic disease resistance in soybean, and showed the genetic makeup of soybean variety 'KS11' comprising the Rmd locus and one MLO gene.

Genetics, genomics and mechanisms responsible for high levels of pyrethroid resistance in Musca domestica

Insecticide resistance is both economically important and evolutionarily interesting phenomenon. Identification of the mutations responsible for resistance allows for highly sensitive resistance monitoring and allows tools to study the forces (population genetics, fitness costs, etc.) that shape the evolution of resistance. Genes coding for insecticide targets have many well-characterized mutations, but the mutations responsible for enhanced detoxification have proven difficult to identify. We employed multiple strategies to identify the mutations responsible for the extraordinarily high permethrin resistance in the KS17-R strain of house fly (Musca domestica): insecticide synergist assays, linkage analysis, bulk segregant analyses (BSA), transcriptomics and long read DNA (Nanopore) sequencing. The >85,100-fold resistance in KS17-R was partially suppressed by the insecticide synergists piperonyl butoxide and S,S,S-tributylphosphorothionate, but not by diethyl maleate nor by injection. This suggests the involvement of target site insensitivity, CYP-mediated resistance, possibly hydrolase mediated resistance and potentially other unknown factors. Linkage analysis identified chromosomes 1, 2, 3 and 5 as having a role in resistance. BSA mapped resistance loci on chromosomes 3 and 5. The locus on chromosome 3 was centered on the voltage sensitive sodium channel. The locus on chromosome 5 was associated with a duplication of multiple detoxification genes. Transcriptomic analyses and long read DNA sequencing revealed overexpressed CYPs and esterases and identified a complex set of structural variants at the chromosome 5 locus.

Genome-wide family prediction unveils molecular mechanisms underlying the regulation of agronomic traits in Urochloa ruziziensis

Tropical forage grasses, particularly those belonging to the Urochloa genus, play a crucial role in cattle production and serve as the main food source for animals in tropical and subtropical regions. The majority of these species are apomictic and tetraploid, highlighting the significance of U. ruziziensis, a sexual diploid species that can be tetraploidized for use in interspecific crosses with apomictic species. As a means to support breeding programs, our study investigates the feasibility of genome-wide family prediction in U. ruziziensis families to predict agronomic traits. Fifty half-sibling families were assessed for green matter yield, dry matter yield, regrowth capacity, leaf dry matter, and stem dry matter across different clippings established in contrasting seasons with varying available water capacity. Genotyping was performed using a genotyping-by-sequencing approach based on DNA samples from family pools. In addition to conventional genomic prediction methods, machine learning and feature selection algorithms were employed to reduce the necessary number of markers for prediction and enhance predictive accuracy across phenotypes. To explore the regulation of agronomic traits, our study evaluated the significance of selected markers for prediction using a tree-based approach, potentially linking these regions to quantitative trait loci (QTLs). In a multiomic approach, genes from the species transcriptome were mapped and correlated to those markers. A gene coexpression network was modeled with gene expression estimates from a diverse set of U. ruziziensis genotypes, enabling a comprehensive investigation of molecular mechanisms associated with these regions. The heritabilities of the evaluated traits ranged from 0.44 to 0.92. A total of 28,106 filtered SNPs were used to predict phenotypic measurements, achieving a mean predictive ability of 0.762. By employing feature selection techniques, we could reduce the dimensionality of SNP datasets, revealing potential genotype-phenotype associations. The functional annotation of genes near these markers revealed associations with auxin transport and biosynthesis of lignin, flavonol, and folic acid. Further exploration with the gene coexpression network uncovered associations with DNA metabolism, stress response, and circadian rhythm. These genes and regions represent important targets for expanding our understanding of the metabolic regulation of agronomic traits and offer valuable insights applicable to species breeding. Our work represents an innovative contribution to molecular breeding techniques for tropical forages, presenting a viable marker-assisted breeding approach and identifying target regions for future molecular studies on these agronomic traits.

A papain-like cysteine protease-released small signal peptide confers wheat resistance to wheat yellow mosaic virus

Wheat yellow mosaic virus (WYMV), a soil-borne pathogen, poses a serious threat to global wheat production. Here, we identify a WYMV resistance gene, TaRD21A, that belongs to the papain-like cysteine protease family. Through genetic manipulation of TaRD21A expression, we establish its positive role in the regulation of wheat to WYMV resistance. Furthermore, our investigation shows that the TaRD21A-mediated plant antiviral response relies on the release of a small peptide catalyzed by TaRD21A protease activity. To counteract wheat resistance, WYMV-encoded nuclear inclusion protease-a (NIa) suppress TaRD21A activity to promote virus infection. In resistant cultivars, a natural variant of TaRD21A features a glycine-to-threonine substitution and this substitution enables the phosphorylation of threonine, thereby weakening the interaction between NIa and TaRD21A, reinforcing wheat resistance against WYMV. Our study not only unveils a WYMV resistance gene but also offers insights into the intricate mechanisms underpinning resistance against WYMV.

Identification of QTL associated with resistance to Phytophthora fruit rot in cucumber (Cucumis sativus L.)

Phytophthora fruit rot (PFR) caused by the soilborne oomycete pathogen, Phytophthora capsici, can cause severe yield loss in cucumber. With no resistant variety available, genetic resources are needed to develop resistant varieties. The goal of this work was to identify quantitative trait loci (QTL) associated with resistance to PFR using multiple genomic approaches and populations. Two types of resistances have been identified: age-related resistance (ARR) and young fruit resistance. ARR occurs at 12-16 days post pollination (dpp), coinciding with the end of exponential fruit growth. A major QTL for ARR was discovered on chromosome 3 and a candidate gene identified based on comparative transcriptomic analysis. Young fruit resistance, which is observed during the state of rapid fruit growth prior to commercial harvest, is a quantitative trait for which multiple QTL were identified. The largest effect QTL, qPFR5.1, located on chromosome 5 was fine mapped to a 1-Mb region. Genome-wide association studies (GWAS) and extreme-phenotype genome-wide association study (XP-GWAS) for young fruit resistance were also performed on a cucumber core collection representing > 96% of the genetic diversity of the USDA cucumber germplasm. Several SNPs overlapped with the QTL identified from QTL-seq analysis on biparental populations. In addition, novel SNPs associated with the resistance were identified from the germplasm. The resistant alleles were found mostly in accessions from India and South Asia, the center of diversity for cucumber. The results from this work can be applied to future disease resistance studies and marker-assisted selection in breeding programs.

Seed color in lettuce is determined by the LsTT2, LsCHS, and Ls2OGD genes from the flavonoid biosynthesis pathway

KEY MESSAGE

The mutated LsTT2 and Ls2OGD genes are responsible for white seeds and yellow seeds in lettuce, respectively. Three LsCHS genes are involved in the biosynthesis of flavonoid in seed coats. Lettuce seeds have several different colors, including black, yellow, and white. The genetic mechanisms underlying color variations of lettuce seeds remain unknown. We used genome-wide association studies (GWAS) and map-based cloning approaches to clone genes controlling the color of lettuce seeds. LsTT2, which encodes an R2R3-MYB transcription factor and is homologous to the TT2 gene in Arabidopsis, was shown to be the causal gene for the variation of black and white seeds in lettuce. A point mutation leads to the lack of stop codon in the LsTT2 transcript, resulting in white seeds. Knockout of the LsTT2 gene converted black seeds to white seeds. The locus controlling yellow seeds was mapped to Chromosome 2. Knockout of two 2-oxoglutarate-dependent dioxygenases (2OGD) genes from the candidate region converted black seeds to yellow seeds, suggesting that these two 2OGD proteins catalyze the conversion of yellow metabolites to black metabolites. We also showed that three LsCHS genes from the candidate region are associated with flavonoid biosynthesis in seeds. Knockout mutants of the three LsCHS genes decreased color intensity. This study provides new insights into the regulation of flavonoid biosynthesis in plants.

Data-driven approaches to improve water-use efficiency and drought resistance in crop plants

With the increasing population, there lies a pressing demand for food, feed and fibre, while the changing climatic conditions pose severe challenges for agricultural production worldwide. Water is the lifeline for crop production; thus, enhancing crop water-use efficiency (WUE) and improving drought resistance in crop varieties are crucial for overcoming these challenges. Genetically-driven improvements in yield, WUE and drought tolerance traits can buffer the worst effects of climate change on crop production in dry areas. While traditional crop breeding approaches have delivered impressive results in increasing yield, the methods remain time-consuming and are often limited by the existing allelic variation present in the germplasm. Significant advances in breeding and high-throughput omics technologies in parallel with smart agriculture practices have created avenues to dramatically speed up the process of trait improvement by leveraging the vast volumes of genomic and phenotypic data. For example, individual genome and pan-genome assemblies, along with transcriptomic, metabolomic and proteomic data from germplasm collections, characterised at phenotypic levels, could be utilised to identify marker-trait associations and superior haplotypes for crop genetic improvement. In addition, these omics approaches enable the identification of genes involved in pathways leading to the expression of a trait, thereby providing an understanding of the genetic, physiological and biochemical basis of trait variation. These data-driven gene discoveries and validation approaches are essential for crop improvement pipelines, including genomic breeding, speed breeding and gene editing. Herein, we provide an overview of prospects presented using big data-driven approaches (including artificial intelligence and machine learning) to harness new genetic gains for breeding programs and develop drought-tolerant crop varieties with favourable WUE and high-yield potential traits.

Heritable epigenetic modification of BpPIN1 is associated with leaf shapes in Betula pendula

The new variety Betula pendula 'Dalecarlica', selected from Betula pendula, shows high ornamental value owing to its lobed leaf shape. In this study, to identify the genetic components of leaf shape formation, we performed bulked segregant analysis and molecular marker-based fine mapping to identify the causal gene responsible for lobed leaves in B. pendula 'Dalecarlica'. The most significant variations associated with leaf shape were identified within the gene BpPIN1 encoding a member of the PIN-FORMED family, responsible for the auxin efflux carrier. We further confirmed the hypomethylation at the promoter region promoting the expression level of BpPIN1, which causes stronger and longer veins and lobed leaf shape in B. pendula 'Dalecarlica'. These results indicated that DNA methylation at the BpPIN1 promoter region is associated with leaf shapes in B. pendula. Our findings revealed an epigenetic mechanism of BpPIN1 in the regulation of leaf shape in Betula  Linn. (birch), which could help in the molecular breeding of ornamental traits.

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.

Next-generation bulked segregant analysis for Breeding 4.0

Functional cloning and manipulation of genes controlling various agronomic traits are important for boosting crop production. Although bulked segregant analysis (BSA) is an efficient method for functional cloning, its low throughput cannot satisfy the current need for crop breeding and food security. Here, we review the rationale and development of conventional BSA and discuss its strengths and drawbacks. We then propose next-generation BSA (NG-BSA) integrating multiple cutting-edge technologies, including high-throughput phenotyping, biological big data, and the use of machine learning. NG-BSA increases the resolution of genetic mapping and throughput for cloning quantitative trait genes (QTGs) and optimizes candidate gene selection while providing a means to elucidate the interaction network of QTGs. The ability of NG-BSA to efficiently batch-clone QTGs makes it an important tool for dissecting molecular mechanisms underlying various traits, as well as for the improvement of Breeding 4.0 strategy, especially in targeted improvement and population improvement of crops.

A major QTL simultaneously increases the number of spikelets per spike and thousand-kernel weight in a wheat line

KEY MESSAGE

A novel and stably expressed QTL QSNS.sicau-SSY-7A for spikelet number per spike in wheat without negative effects on thousand-kernel weight was identified and validated in different genetic backgrounds. Spikelet number per spike (SNS) is an important determinant of yield in wheat. In the present study, we combined bulked segregant analysis (BSA) and the wheat 660 K single-nucleotide polymorphism (SNP) array to rapidly identify genomic regions associated with SNS from a recombinant inbred line (RIL) population derived from a cross between the wheat lines S849-8 and SY95-71. A genetic map was constructed using Kompetitive Allele Specific PCR markers in the SNP-enriched region on the long arm of chromosome 7A. A major and stably expressed QTL, QSNS.sicau-SSY-7A, was detected in multiple environments. It was located in a 1.6 cM interval on chromosome arm 7AL flanked by the markers AX-109983514 and AX-109820548. This QTL explained 6.86-15.72% of the phenotypic variance, with LOD values ranging from 3.66 to 8.66. Several genes associated with plant growth and development were identified in the interval where QSNS.sicau-SSY-7A was located on the 'Chinese Spring' wheat and wild emmer reference genomes. Furthermore, the effects of QSNS.sicau-SSY-7A and WHEAT ORTHOLOG OFAPO1(WAPO1) on SNS were analyzed. Interestingly, QSNS.sicau-SSY-7A significantly increased SNS without negative effects on thousand-kernel weight, anthesis date and plant height, demonstrating its great potential for breeding aimed at improving grain yield. Taken together, these results indicate that QSNS.sicau-SSY-7A is a promising locus for yield improvement, and its linkage markers are helpful for fine mapping and molecular breeding.

QTL and PACE analyses identify candidate genes for anthracnose resistance in tomato

Anthracnose, caused by the fungal pathogen Colletotrichum spp., is one of the most significant tomato diseases in the United States and worldwide. No commercial cultivars with anthracnose resistance are available, limiting resistant breeding. Cultivars with genetic resistance would significantly reduce crop losses, reduce the use of fungicides, and lessen the risks associated with chemical application. A recombinant inbred line (RIL) mapping population (N=243) has been made from a cross between the susceptible US28 cultivar and the resistant but semiwild and small-fruited 95L368 to identify quantitative trait loci (QTLs) associated with anthracnose resistance. The RIL population was phenotyped for resistance by inoculating ripe field-harvested tomato fruits with Colletotrichum coccodes for two seasons. In this study, we identified twenty QTLs underlying resistance, with a range of phenotypic variance of 4.5 to 17.2% using a skeletal linkage map and a GWAS. In addition, a QTLseq analysis was performed using deep sequencing of extreme bulks that validated QTL positions identified using traditional mapping and resolved candidate genes underlying various QTLs. We further validated AP2-like ethylene-responsive transcription factor, N-alpha-acetyltransferase (NatA), cytochrome P450, amidase family protein, tetratricopeptide repeat, bHLH transcription factor, and disease resistance protein RGA2-like using PCR allelic competitive extension (PACE) genotyping. PACE assays developed in this study will enable high-throughput screening for use in anthracnose resistance breeding in tomato.

Mapping of the Cladosporium fulvum resistance gene Cf-16, a major gene involved in leaf mold disease in tomato

Tomato (Solanum lycopersicum) is widely cultivated and consumed worldwide. Tomato leaf mold, caused by Cladosporium fulvum, is one of the most devastating diseases in tomato production. At present, some tomato leaf mold resistance (Cf series) genes used in production gradually lose resistance due to the continuous and rapid differentiation of C. fulvum physiological races. The Cf-16 gene derived from the "Ontario7816" tomato cultivar has shown effective resistance in field trials for many years, but few studies have reported on the mapping of the Cf-16 gene, which has not been cloned, limiting its utilization in tomato breeding. Here, we mapped Cf-16 using a novel comprehensive strategy including bulk segregation analysis (BSA), genome resequencing and SSR molecular markers. A genetic analysis revealed that Cf-16 resistance in "Ontario7816" is controlled by one major dominant locus. The Cf-16 gene was mapped in a region of 2.6 cM at chromosome 6 between two markers, namely, TGS447 and TES312, by using an F2 population from a cross between the resistant cultivar "Ontario7816" and susceptible line "Moneymaker." Two nucleotide-binding-site-leucine-rich repeat (NBS-LRR) resistance genes, namely, XM_004240667.3 and XM_010323727.1, were identified in this interval. They are strong candidates for the Cf-16 gene. The mapping of Cf-16 may speed up its utilization for breeding resistant tomato varieties and represents an important step forward in our understanding of the mechanism underlying resistance to tomato leaf mold.

A BrLINE1-RUP insertion in BrCER2 alters cuticular wax biosynthesis in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Glossiness is an important quality-related trait of Chinese cabbage, which is a leafy vegetable crop in the family Brassicaceae. The glossy trait is caused by abnormal cuticular wax accumulation. In this study, on the basis of a bulked segregant analysis coupled with next-generation sequencing (BSA-seq) and fine-mapping, the most likely candidate gene responsible for the glossy phenotype of Chinese cabbage was identified. It was subsequently named Brcer2 because it is homologous to AtCER2 (At4g24510). A bioinformatics analysis indicated a long interspersed nuclear element 1 (LINE-1) transposable element (named BrLINE1-RUP) was inserted into the first exon of Brcer2 in HN19-G via an insertion-mediated deletion mechanism, which introduced a premature termination codon. Gene expression analysis showed that the InDel mutation of BrCER2 reduced the transcriptional expression levels of Brcer2 in HN19-G. An analysis of cuticular waxes suggested that a loss-of-function mutation to BrCER2 in Chinese cabbage leads to a severe decrease in the abundance of very-long-chain-fatty-acids (> C28), resulting in the production of a cauline leaf, inflorescence stem, flower, and pistil with a glossy phenotype. These findings imply the insertion of the LINE-1 transposable element BrLINE1-RUP into BrCER2 can modulate the waxy traits of Chinese cabbage plants.

Fine mapping and genetic analysis identified a C2H2-type zinc finger as a candidate gene for heading date regulation in wheat

KEY MESSAGE

A minor-effect QTL, Qhd.2AS, that affects heading date in wheat was mapped to a genomic interval of 1.70-Mb on 2AS, and gene analysis indicated that the C2H2-type zinc finger protein gene TraesCS2A02G181200 is the best candidate for Qhd.2AS. Heading date (HD) is a complex quantitative trait that determines the regional adaptability of cereal crops, and identifying the underlying genetic elements with minor effects on HD is important for improving wheat production in diverse environments. In this study, a minor QTL for HD that we named Qhd.2AS was detected on the short arm of chromosome 2A by Bulked Segregant Analysis and validated in a recombinant inbred population. Using a segregating population of 4894 individuals, Qhd.2AS was further delimited to an interval of 0.41 cM, corresponding to a genomic region spanning 1.70 Mb (from 138.87 to 140.57 Mb) that contains 16 high-confidence genes based on IWGSC RefSeq v1.0. Analyses of sequence variations and gene transcription indicated that TraesCS2A02G181200, which encodes a C2H2-type zinc finger protein, is the best candidate gene for Qhd.2AS that influences HD. Screening a TILLING mutant library identified two mutants with premature stop codons in TraesCS2A02G181200, both of which exhibited a delay in HD of 2-4 days. Additionally, variations in its putative regulatory sites were widely present in natural accession, and we also identified the allele which was positively selected during wheat breeding. Epistatic analyses indicated that Qhd.2AS-mediated HD variation is independent of VRN-B1 and environmental factors. Phenotypic investigation of homozygous recombinant inbred lines (RILs) and F2:3 families showed that Qhd.2AS has no negative effect on yield-related traits. These results provide important cues for refining HD and therefore improving yield in wheat breeding programs and will deepen our understanding of the genetic regulation of HD in cereal plants.

Development of SNP markers linked to purple blotch resistance for marker-assisted selection in onion (Allium cepa L.) breeding

Purple blotch (PB), caused by Alternaria porri (Ellis) Cifferi, is one of the most destructive diseases of onion worldwide. Rapid development and deployment of resistant onion varieties is the most effective approach to control this disease. A single dominant gene, ApR1 was previously linked to PB resistance in onion cultivar 'Arka Kalyan'. In this study, an advanced RIL population derived from a cross between the resistant (Arka Kalyan) and susceptible (Agrifound Rose) cultivar of onion was used to fine map the resistant locus with SNP markers. Twenty plants from the RIL population, ten each with disease resistance and susceptibility trait, were subjected to restriction site-associated DNA sequencing (RAD-Seq) and generated 7388 single nucleotide polymorphisms (SNPs). Correlation analysis between marker genotypes and PB disease phenotype on the 20 plants identified 27 SNPs as candidate markers linked to ApR1 gene for PB resistance. Six candidate SNPs were converted to Kompetitive Allele-Specific PCR (KASP) markers designated as ApRsnip5, ApRsnip8, ApRsnip14, ApRsnip21, ApRsnip23 and ApRsnip25. Marker-trait association based on disease phenotyping and KASP genotyping data on 153 RILs confirmed that all six KASP markers were tightly associated with ApR1 gene within the genetic distance of 1.3 CentiMorgan (cM). ApRsnip14 co-segregated with the ApR1 locus. Further, the six KASP markers were tested on 27 onion lines with different genetic backgrounds. ApRsnip14, ApRsnip21, ApRsnip5 and ApRsnip23 not only showed the correct resistance allele in 3 resistance genotypes, but also clustered together in the remaining 24 susceptible lines. Alternatively, ApRsnip8 and ApRsnip25 exhibited false positives in two onion lines which do not have the R-gene. Overall, our results suggest that ApRsnip14 and ApRsnip23 with their close linkage to ApR1 locus and greater applicability on breeding germplasm are recommended in marker-assisted selection for PB resistance in onion breeding program.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03562-7.

Transcript profiling of hairy vetch (Vicia villosa Roth) identified interesting genes for seed dormancy

Hairy vetch, a diploid annual legume species, has a robust growth habit, high biomass yield, and winter hardy characteristics. Seed hardness is a major constraint for growing hairy vetch commercially. Hard seeded cultivars are valuable as forages, whereas soft seeded and shatter resistant cultivars have advantages for their use as a cover crop. Transcript analysis of hairy vetch was performed to understand the genetic mechanisms associated with important hairy vetch traits. RNA was extracted from leaves, flowers, immature pods, seed coats, and cotyledons of contrasting soft and hard seeded "AU Merit" plants. A range of 31.22-79.18 Gb RNA sequence data per tissue sample were generated with estimated coverage of 1040-2639×. RNA sequence assembly and mapping of the contigs against the Medicago truncatula (V4.0) genome identified 76,422 gene transcripts. A total of 24,254 transcripts were constitutively expressed in hairy vetch tissues. Key genes, such as KNOX4 (a class II KNOTTED-like homeobox KNOXII gene), qHs1 (endo-1,4-β-glucanase), GmHs1-1 (calcineurin-like metallophosphoesterase), chitinase, shatterproof 1 and 2 (SHP1, SHP2), shatter resistant 1-5 (SHAT1-5)(NAC transcription factor), PDH1 (prephenate dehydrogenase 1), and pectin methylesterases with a potential role in seed hardness and pod shattering, were further explored based on genes involved in seed hardness from other species to query the hairy vetch transcriptome data. Identification of interesting candidate genes in hairy vetch can facilitate the development of improved cultivars with desirable seed characteristics for use as a forage and as a cover crop.

Rapid identification of a major locus qPRL-C06 affecting primary root length in Brassica napus by QTL-seq

BACKGROUND AND AIMS

Brassica napus is one of the most important oilseed crops worldwide. Seed yield of B. napus significantly correlates with the primary root length (PRL). The aims of this study were to identify quantitative trait loci (QTLs) for PRL in B. napus.

METHODS

QTL-seq and conventional QTL mapping were jointly used to detect QTLs associated with PRL in a B. napus double haploid (DH) population derived from a cross between 'Tapidor' and 'Ningyou 7'. The identified major locus was confirmed and resolved by an association panel of B. napus and an advanced backcross population. RNA-seq analysis of two long-PRL lines (Tapidor and TN20) and two short-PRL lines (Ningyou 7 and TN77) was performed to identify differentially expressed genes in the primary root underlying the target QTLs.

KEY RESULTS

A total of 20 QTLs impacting PRL in B. napus grown at a low phosphorus (P) supply were found by QTL-seq. Eight out of ten QTLs affecting PRL at a low P supply discovered by conventional QTL mapping could be detected by QTL-seq. The locus qPRL-C06 identified by QTL-seq was repeatedly detected at both an optimal P supply and a low P supply by conventional QTL mapping. This major constitutive QTL was further confirmed by regional association mapping. qPRL-C06 was delimited to a 0.77 Mb genomic region on chromosome C06 using an advanced backcross population. A total of 36 candidate genes within qPRL-C06 were identified that showed variations in coding sequences and/or exhibited significant differences in mRNA abundances in primary root between the long-PRL and short-PRL lines, including five genes involved in phytohormone biosynthesis and signaling.

CONCLUSIONS

These results both demonstrate the power of the QTL-seq in rapid QTL detection for root traits and will contribute to marker-assisted selective breeding of B. napus cultivars with increased PRL.

Identification of simple sequence repeat markers linked to heat tolerance in rice using bulked segregant analysis in F2 population of NERICA-L 44 × Uma

The damage caused by high temperature is one of the most important abiotic stress affecting rice production. Reproductive stage of rice is highly susceptible to high temperature. The present investigation was undertaken to identify polymorphic microsatellite markers (SSR) associated with heat tolerance. The rice cultivars NERICA- L 44 (heat tolerant) and Uma (heat susceptible) were crossed to generate F¹ and F² populations. The F² population was subjected to heat stress at >38°C and the 144 F² plants were evaluated for their tolerance. The results note that the mean of the F² population was influenced by the tolerant parent with regards to the traits of plant height, membrane stability index, photosynthetic rate, stomatal conductance, evapotranspiration rate, pollen viability, spikelet fertility and 1000 grain weight. Ten each of the extremely susceptible and tolerant plants were selected based on the spikelet fertility percentage. Their DNA was pooled into tolerant and susceptible bulks and Bulked Segregant Analysis (BSA) was carried out using 100 SSR markers to check for polymorphism. The survey revealed a polymorphism of 18% between the parents. RM337, RM10793, RM242, RM5749, RM6100, RM490, RM470, RM473, RM222 and RM556 are some of the prominent markers that were found to be polymorphic between the parents and the bulks. We performed gene annotation and enrichment analysis of identified polymorphic markers. Result revealed that the sequence specific site of that chromosome mostly enriched with biological processes like metabolic pathway, molecular mechanism, and subcellular function. Among that RM337 was newly reported marker for heat tolerance. Expression analysis of two genes corresponds to RM337 revealed that LOP1 (LOC_Os08g01330) was linked to high temperature tolerance in rice. The results demonstrate that BSA using SSR markers is useful in identifying genomic regions that contribute to thermotolerance.

Identification of a stable major sucrose-related QTL and diagnostic marker for flavor improvement in peanut

An InDel marker closely linked with a major and stable quantitative trait locus (QTL) on chromosome A08, qSUCA08.2, controlling sucrose content will benefit peanut flavor improvement. Sucrose is the main soluble sugar in mature peanut kernel, and its content is a key determinant of flavor. However, the genetic basis of sucrose content in peanut remains poorly understood, which limits the progress of flavor improvement. In the present study, two genomic regions (qSUCA08a and qSUCB06a) for sucrose content on chromosomes A08 and B06 were identified by QTL-seq in a RIL population derived from a cross between Zhonghua 10 and ICG 12625. In the interval of qSUCB06a, QTL qSUCB06.2 was detected through QTL mapping in a single environment. The qSUCA08a was further dissected into 3 adjacent genomic regions using linkage analysis including a major QTL qSUCA08.2 explaining 5.43-17.84% phenotypic variation across five environments. A 61-bp insertion at position 35,099,320 in the higher sucrose parent ICG 12625 was found in qSUCA08.2. An InDel marker SUC.InDel.A08 based on the insertion/deletion polymorphism was developed and validated within a natural population containing 172 peanut cultivars in two environments. The mean sucrose content of 93 cultivars with ICG 12625 allele was significantly higher than that of 79 cultivars with Zhonghua 10 allele. The qSUCA08.2 corresponding to a 2.11 Mb interval harbored 110 genes. Among these genes, a total of 19 genes were considered as candidate genes including 5 non-synonymous mutation genes and 14 differentially expressed genes during seed development. These results provide new insights into the genetic basis of sucrose regulation in peanut and benefit the breeding program for developing new varieties with excellent flavor.

Heritability and gene functions associated with sclerotia formation of Rhizoctonia solani AG-7 using whole genome sequencing and genome-wide association study

Sclerotia are specialized fungal structures formed by pigmented and aggregated hyphae, which can survive under unfavourable environmental conditions and serve as the primary inocula for several phytopathogenic fungi including Rhizoctonia solani. Among 154 R. solani anastomosis group 7 (AG-7) isolates collected in fields, the sclerotia-forming capability regarding sclerotia number and sclerotia size varied in the fungal population, but the genetic makeup of these phenotypes remained unclear. As limited studies have focused on the genomics of R. solani AG-7 and the population genetics of sclerotia formation, this study completed the whole genome sequencing and gene prediction of R. solani AG-7 using the Oxford NanoPore and Illumina RNA sequencing. Meanwhile, a high-throughput image-based method was established to quantify the sclerotia-forming capability, and the phenotypic correlation between sclerotia number and sclerotia size was low. A genome-wide association study identified three and five significant SNPs associated with sclerotia number and size in distinct genomic regions, respectively. Of these significant SNPs, two and four showed significant differences in the phenotypic mean separation for sclerotia number and sclerotia size, respectively. Gene ontology enrichment analysis focusing on the linkage disequilibrium blocks of significant SNPs identified more categories related to oxidative stress for sclerotia number, and more categories related to cell development, signalling and metabolism for sclerotia size. These results indicated that different genetic mechanisms may underlie these two phenotypes. Moreover, the heritability of sclerotia number and sclerotia size were estimated for the first time to be 0.92 and 0.31, respectively. This study provides new insights into the heritability and gene functions related to the development of sclerotia number and sclerotia size, which could provide additional knowledge to reduce fungal residues in fields and achieve sustainable disease management.

Identification of a Major QTL-Controlling Resistance to the Subtropical Race 4 of Fusarium oxysporum f. sp. cubense in Musa acuminata ssp. malaccensis

Vascular wilt caused by the ascomycete fungal pathogen Fusarium oxysporum f. sp. cubense (Foc) is a major constraint of banana production around the world. The virulent race, namely Tropical Race 4, can infect all Cavendish-type banana plants and is now widespread across the globe, causing devastating losses to global banana production. In this study, we characterized Foc Subtropical Race 4 (STR4) resistance in a wild banana relative which, through estimated genome size and ancestry analysis, was confirmed to be Musa acuminata ssp. malaccensis. Using a self-derived F₂ population segregating for STR4 resistance, quantitative trait loci sequencing (QTL-seq) was performed on bulks consisting of resistant and susceptible individuals. Changes in SNP index between the bulks revealed a major QTL located on the distal end of the long arm of chromosome 3. Multiple resistance genes are present in this region. Identification of chromosome regions conferring resistance to Foc can facilitate marker assisted selection in breeding programs and paves the way towards identifying genes underpinning resistance.

Integrating linkage mapping and comparative transcriptome analysis for discovering candidate genes associated with salt tolerance in rice

Salinity is one of the most widespread abiotic stresses affecting rice productivity worldwide. Understanding the genetic basis of salt tolerance is key for breeding salt-tolerant rice varieties. Numerous QTLs have been identified to help dissect rice salt-tolerance genetic mechanisms, yet only rare genes located in significant QTLs have been thoroughly studied or fine-mapped. Here, a combination of linkage mapping and transcriptome profiling analysis was used to identify salt tolerance-related functional candidate genes underlying stable QTLs. A recombinant inbred line (RIL) population derived from a cross between Jileng 1 (salt-sensitive) and Milyang 23 (salt-tolerant) was constructed. Subsequently, a high-density genetic map was constructed by using 2921 recombination bin markers developed from whole genome resequencing. A total of twelve QTLs controlling the standard evaluation score under salt stress were identified by linkage analysis and distributed on chromosomes 2, 3, 4, 6, 8 and 11. Notably, five QTL intervals were detected as environmentally stable QTLs in this study, and their functions were verified by comparative transcriptome analysis. By comparing the transcriptome profiles of the two parents and two bulks, we found 551 salt stress-specific differentially expressed genes. Among them, fifteen DEGs located in stable QTL intervals were considered promising candidate genes for salt tolerance. According to gene annotations, the gene OsRCI2-8(Os06g0184800) was the most promising, as it is known to be associated with salt stress, and its differential expression between the tolerant and sensitive RIL bulks highlights its important role in salt stress response pathways. Our findings provide five stable salt tolerance-related QTLs and one promising candidate gene, which will facilitate breeding for improved salt tolerance in rice varieties and promote the exploration of salt stress tolerance mechanisms in rice.

Fine mapping and candidate gene analysis of proportion of four-seed pods by soybean CSSLs

Soybean yield, as one of the most important and consistent breeding goals, can be greatly affected by the proportion of four-seed pods (PoFSP). In this study, QTL mapping was performed by PoFSP data and BLUE (Best Linear Unbiased Estimator) value of the chromosome segment substitution line population (CSSLs) constructed previously by the laboratory from 2016 to 2018, and phenotype-based bulked segregant analysis (BSA) was performed using the plant lines with PoFSP extreme phenotype. Totally, 5 ICIM QTLs were repeatedly detected, and 6 BSA QTLs were identified in CSSLs. For QTL (qPoFSP13-1) repeated in ICIM and BSA results, the secondary segregation populations were constructed for fine mapping and the interval was reduced to 100Kb. The mapping results showed that the QTL had an additive effect of gain from wild parents. A total of 14 genes were annotated in the delimited interval by fine mapping. Sequence analysis showed that all 14 genes had genetic variation in promoter region or CDS region. The qRT-PCR results showed that a total of 5 candidate genes were differentially expressed between the plant lines having antagonistic extreme phenotype (High PoFSP > 35.92%, low PoFSP< 17.56%). The results of haplotype analysis showed that all five genes had two or more major haplotypes in the resource population. Significant analysis of phenotypic differences between major haplotypes showed all five candidate genes had haplotype differences. And the genotypes of the major haplotypes with relatively high PoFSP of each gene were similar to those of wild soybean. The results of this study were of great significance to the study of candidate genes affecting soybean PoFSP, and provided a basis for the study of molecular marker-assisted selection (MAS) breeding and four-seed pods domestication.

Glycoalkaloid Composition and Flavonoid Content as Driving Forces of Phytotoxicity in Diploid Potato

Despite their advantages, biotechnological and omic techniques have not been applied often to characterize phytotoxicity in depth. Here, we show the distribution of phytotoxicity and glycoalkaloid content in a diploid potato population and try to clarify the source of variability of phytotoxicity among plants whose leaf extracts have a high glycoalkaloid content against the test plant species, mustard. Six glycoalkaloids were recognized in the potato leaf extracts: solasonine, solamargine, α-solanine, α-chaconine, leptinine I, and leptine II. The glycoalkaloid profiles of the progeny of the group with high phytotoxicity differed from those of the progeny of the group with low phytotoxicity, which stimulated mustard growth. RNA sequencing analysis revealed that the upregulated flavonol synthase/flavonone 3-hydroxylase-like gene was expressed in the progeny of the low phytotoxicity group, stimulating plant growth. We concluded that the metabolic shift among potato progeny may be a source of different physiological responses in mustard. The composition of glycoalkaloids, rather than the total glycoalkaloid content itself, in potato leaf extracts, may be a driving force of phytotoxicity. We suggest that, in addition to glycoalkaloids, other metabolites may shape phytotoxicity, and we assume that these metabolites may be flavonoids.

Recognition of Pep-13/25 MAMPs of Phytophthora localizes to an RLK locus in Solanum microdontum

Pattern-triggered immunity (PTI) in plants is mediated by cell surface-localized pattern recognition receptors (PRRs) upon perception of microbe-associated molecular pattern (MAMPs). MAMPs are conserved molecules across microbe species, or even kingdoms, and PRRs can confer broad-spectrum disease resistance. Pep-13/25 are well-characterized MAMPs in Phytophthora species, which are renowned devastating oomycete pathogens of potato and other plants, and for which genetic resistance is highly wanted. Pep-13/25 are derived from a 42 kDa transglutaminase GP42, but their cognate PRR has remained unknown. Here, we genetically mapped a novel surface immune receptor that recognizes Pep-25. By using effectoromics screening, we characterized the recognition spectrum of Pep-13/25 in diverse Solanaceae species. Response to Pep-13/25 was predominantly found in potato and related wild tuber-bearing Solanum species. Bulk-segregant RNA sequencing (BSR-Seq) and genetic mapping the response to Pep-25 led to a 0.081 cM region on the top of chromosome 3 in the wild potato species Solanum microdontum subsp. gigantophyllum. Some BAC clones in this region were isolated and sequenced, and we found the Pep-25 receptor locates in a complex receptor-like kinase (RLK) locus. This study is an important step toward the identification of the Pep-13/25 receptor, which can potentially lead to broad application in potato and various other hosts of Phytophthora species.

Mapping and characterization of a novel adult-plant leaf rust resistance gene LrYang16G216 via bulked segregant analysis and conventional linkage method

A novel adult-plant leaf rust resistance gene LrYang16G216 on wheat chromosome 6BL was identified and mapped to a 0.59 cM genetic interval by BSA and conventional linkage method. Leaf rust (Puccinia triticina) is one of the most devastating fungal diseases of wheat (Triticum aestivum L.). Discovery and identification of new resistance genes is essential to develop disease-resistant cultivars. An advanced breeding line Yang16G216 was previously identified to confer adult-plant resistance (APR) to leaf rust. In this research, a recombinant inbred line (RIL) population was constructed from the cross between Yang16G216 and a highly susceptible line Yang16M6393, and genotyped with exome capture sequencing and 55 K SNP array. Through bulked segregant analysis (BSA) and genetic linkage mapping, a stable APR gene, designated as LrYang16G216, was detected and mapped to the distal region of chromosome arm 6BL with a genetic interval of 2.8 cM. For further verification, another RIL population derived from the cross between Yang16G216 and a susceptible wheat variety Yangmai 29 was analyzed using the enriched markers in the target interval, and LrYang16G216 was further narrowed to a 0.59 cM genetic interval flanked by the KASP markers Ax109403980 and Ax95083494, corresponding to the physical position 712.34-713.94 Mb in the Chinese Spring reference genome, in which twenty-six disease resistance-related genes were annotated. Based on leaf rust resistance spectrum, mapping data and physical location, LrYang16G216 was identified to be a novel and effective APR gene. The LrYang16G216 with linked markers will be useful for marker-assisted selection in wheat resistance breeding.

A HST1-like gene controls tiller angle through regulating endogenous auxin in common wheat

Tiller angle is one of the most important agronomic traits and one key factor for wheat ideal plant architecture, which can both increase photosynthetic efficiency and greatly enhance grain yield. Here, a deacetylase HST1-like (TaHST1L) gene controlling wheat tiller angle was identified by the combination of a genome-wide association study (GWAS) and bulked segregant analysis (BSA). Ethyl methane sulfonate (EMS)-mutagenized tetraploid wheat lines with the premature stop codon of TaHST1L exhibited significantly smaller tiller angles than the wild type. TaHST1L-overexpressing (OE) plants exhibited significantly larger tiller angles and increased tiller numbers in both winter and spring wheat, while TaHST1L-silenced RNAi plants displayed significantly smaller tiller angles and decreased tiller numbers. Moreover, TaHST1L strongly interacted with TaIAA17 and inhibited its expression at the protein level, and thus possibly improved the content of endogenous auxin in the basal tissue of tillers. The transcriptomics and metabolomics results indicated that TaHST1L might change plant architecture by mediating auxin signal transduction and regulating endogenous auxin levels. In addition, a 242-bp insertion/deletion (InDel) in the TaHST1L-A1 promoter altered transcriptional activity and TaHST1L-A1b allele with the 242-bp insertion widened the tiller angle of TaHST1L-OE transgenic rice plants. Wheat varieties with TaHST1L-A1b allele possessed the increased tiller angle and grain yield. Further analysis in wheat and its progenitors indicated that the 242-bp InDel possibly originated from wild emmer and was strongly domesticated in the current varieties. Therefore, TaHST1L involved in the auxin signalling pathway showed the big potential to improve wheat yield by controlling plant architecture.

Identification of genes for seed isoflavones based on bulk segregant analysis sequencing in soybean natural population

We identified four hub genes for isoflavone biosynthesis based on BSA-seq and WGCNA methods and validated that GmIE3-1 positively contribute to isoflavone accumulation in soybean. Soybean isoflavones are secondary metabolites of great interest owing to their beneficial impact on human health. Herein, we profiled the seed isoflavone content by HPLC in 1551 soybean accessions grown in two locations for two years and constructed two extreme pools with high (4065.1 µg g^-1) and low (1427.23 µg g^-1) isoflavone contents to identify candidate genes involved in isoflavone biosynthesis pathways using bulk segregant analysis sequencing (BSA-seq) approach. The results showed that the average sequencing depths were 50.3× and 65.7× in high and low pools, respectively. A total of 23,626 polymorphic SNPs and 5299 InDels were detected between both pools and 1492 genes with different variations were identified. Based on differential genes in BSA-seq and weighted gene co-expression network analysis (WGCNA), four hub genes, Glyma.06G290400 (designated as GmIE3-1), Glyma.01G239200, Glyma.01G241500, Glyma.13G256100 were identified, encoding E3 ubiquitin-protein ligase, arm repeat protein interacting with ABF2, zinc metallopeptidase EGY3, and dynamin-related protein 3A, respectively. The allelic variation in GmIE3-1 showed a significant influence on isoflavone accumulation. The virus-induced gene silencing (VIGS) and RNAi hairy root transformation of GmIE3-1 revealed partial suppression of this gene could cause a significant decrease (P < 0.0001) of total isoflavone content, suggesting GmIE3-1 is a positive regulator for isoflavones. The present study demonstrated that the BSA-seq approach combined with WGCNA, VIGS and hairy root transformation can efficiently identify isoflavone candidate genes in soybean natural population.

Identification of two tandem genes associated with primary rosette branching in flowering Chinese cabbage

Branching is an important agronomic trait determining plant architecture and yield; however, the molecular mechanisms underlying branching in the stalk vegetable, flowering Chinese cabbage, remain unclear. The present study identified two tandem genes responsible for primary rosette branching in flowering Chinese cabbage by GradedPool-Seq (GPS) combined with Kompetitive Allele Specific PCR (KASP) genotyping. A 900 kb candidate region was mapped in the 28.0-28.9 Mb interval of chromosome A07 through whole-genome sequencing of three graded-pool samples from the F₂ population derived by crossing the branching and non-branching lines. KASP genotyping narrowed the candidate region to 24.6 kb. Two tandem genes, BraA07g041560.3C and BraA07g041570.3C, homologous to AT1G78440 encoding GA2ox1 oxidase, were identified as the candidate genes. The BraA07g041560.3C sequence was identical between the branching and non-branching lines, but BraA07g041570.3C had a synonymous single nucleotide polymorphic (SNP) mutation in the first exon (290^th bp, A to G). In addition, an ERE cis-regulatory element was absent in the promoter of BraA07g041560.3C, and an MYB cis-regulatory element in the promoter of BraA07g041570.3C in the branching line. Gibberellic acid (GA₃) treatment decreased the primary rosette branch number in the branching line, indicating the significant role of GA in regulating branching in flowering Chinese cabbage. These results provide valuable information for revealing the regulatory mechanisms of branching and contributing to the breeding programs of developing high-yielding species in flowering Chinese cabbage.

Bulk segregant linkage mapping for rodent and human malaria parasites

In 2005 Richard Carter's group surprised the malaria genetics community with an elegant approach to rapidly mapping the genetic basis of phenotypic traits in rodent malaria parasites. This approach, which he termed "linkage group selection", utilized bulk pools of progeny, rather than individual clones, and exploited simple selection schemes to identify genome regions underlying resistance to drug treatment (or other phenotypes). This work was the first application of "bulk segregant" methodologies for genetic mapping in microbes: this approach is now widely used in yeast, and across multiple recombining pathogens ranging from Aspergillus fungi to Schistosome parasites. Genetic crosses of human malaria parasites (for which Richard Carter was also a pioneer) can now be conducted in humanized mice, providing new opportunities for exploiting bulk segregant approaches for a wide variety of malaria parasite traits. We review the application of bulk segregant approaches to mapping malaria parasite traits and suggest additional developments that may further expand the utility of this powerful approach.

Fine Mapping of Stripe-Rust-Resistance Gene YrJ22 in Common Wheat by BSR-Seq and MutMap-Based Sequencing

Identification and accurate mapping of new resistance genes are essential for gene pyramiding in wheat breeding. The YrJ22 gene is a dominant stripe-rust-resistance gene located at the distal end of chromosome 2AL, which was identified in a leading Chinese-wheat variety, Jimai 22, showing high resistance to CYR32, a prevalent race of Puccinia striiformis tritici (Pst) in China. In the current study, 15 F₁ and 2273 F₂ plants derived from the cross of Jimai 22/Avocet S were used for the fine-mapping of YrJ22. The RNA-Seq of resistant and susceptible bulks of F₂ plants (designated BSR-Seq) identified 10 single-nucleotide polymorphisms (SNP) in a 12.09 Mb physical interval on chromosome 2AL. A total of 1022 EMS-induced M₃ lines of Jimai 22 were screened, to identify susceptible mutants for MutMap analysis. Four CAPS markers were developed from SNPs identified using BSR-Seq and MutMap. A linkage map for YrJ22 was constructed with 11 CAPS/STS and three SSR markers. YrJ22 was located at a 0.9 cM genetic interval flanked by markers H736 and H400, corresponding to a 340.46 kb physical region (768.7-769.0 Mb), including 13 high-confidence genes based on the Chinese Spring reference genome. TraesCS2A01G573200 is a potential candidate-gene, according to linkage and quantitative real-time PCR (qPCR) analyses. The CAPS marker H732 designed from an SNP in TraesCS2A01G573200 co-segregated with YrJ22. These results provide a useful stripe-rust-resistance gene and molecular markers for marker-assisted selection in wheat breeding and for further cloning of the gene.

Genome-wide association study and genetic mapping of BhWAX conferring mature fruit cuticular wax in wax gourd

BACKGROUND

Wax gourd [Benincasa hispida (Thunb) Cogn. (2n = 2x = 24)] is an economically important vegetable crop of genus Benincasa in the Cucurbitaceae family. Fruit is the main consumption organ of wax gourd. The mature fruit cuticular wax (MFCW) is an important trait in breeding programs, which is also of evolutionary significance in wax gourd. However, the genetic architecture of this valuable trait remains unrevealed.

RESULTS

In this study, genetic analysis revealed that the inheritance of MFCW was controlled by a single gene, with MFCW dominant over non-MFCW, and the gene was primarily named as BhWAX. Genome-wide association study (GWAS) highlighted a 1.1 Mb interval on chromosome 9 associated with MFCW in wax gourd germplasm resources. Traditional fine genetic mapping delimited BhWAX to a 0.5 Mb region containing 12 genes. Based on the gene annotation, expression analysis and co-segregation analysis, Bhi09G001428 that encodes a membrane bound O-acyltransferase (MBOAT) was proposed as the candidate gene for BhWAX. Moreover, it was demonstrated that the efficiency of a cleaved amplified polymorphic sequences (CAPS) marker in the determination of MFCW in wax gourd reached 80%.

CONCLUSIONS

In closing, the study identified the candidate gene controlling MFCW and provided an efficient molecular marker for the trait in wax gourd for the first time, which will not only be beneficial for functional validation of the gene and marker-assisted breeding of wax gourd, but also lay a foundation for analysis of its evolutionary meaning among cucurbits.

A reference-anchored oat linkage map reveals quantitative trait loci conferring adult plant resistance to crown rust (Puccinia coronata f. sp. avenae)

We mapped three adult plant resistance (APR) loci on oat chromosomes 4D and 6C and developed flanking KASP/PACE markers for marker-assisted selection and gene pyramiding. Using sequence orthology search and the available oat genomic and transcriptomic data, we surveyed these genomic regions for genes that may control disease resistance. Sources of durable disease resistance are needed to minimize yield losses in cultivated oat caused by crown rust (Puccinia coronata f. sp. avenae). In this study, we developed five oat recombinant inbred line mapping populations to identify sources of adult plant resistance from crosses between five APR donors and Otana, a susceptible variety. The preliminary bulk segregant mapping based on allele frequencies showed two regions in linkage group Mrg21 (Chr4D) that are associated with the APR phenotype in all five populations. Six markers from these regions in Chr4D were converted to high-throughput allele specific PCR assays and were used to genotype all individuals in each population. Simple interval mapping showed two peaks in Chr4D, named QPc.APR-4D.1 and QPc.APR-4D.2, which were detected in the OtanaA/CI4706-2 and OtanaA/CI9416-2 and in the Otana/PI189733, OtanaD/PI260616, and OtanaA/CI8000-4 populations, respectively. These results were validated by mapping two entire populations, Otana/PI189733 and OtanaA/CI9416, genotyped using Illumina HiSeq, in which polymorphisms were called against the OT3098 oat reference genome. Composite interval mapping results confirmed the presence of the two quantitative trait loci (QTL) located on oat chromosome 4D and an additional QTL with a smaller effect located on chromosome 6C. This mapping approach also narrowed down the physical intervals to between 5 and 19 Mb, and indicated that QPc.APR-4D.1, QPc.APR-4D.2, and QPc.APR-6C explained 43.4%, 38.5%, and 21.5% of the phenotypic variation, respectively. In a survey of the gene content of each QTL, several clusters of disease resistance genes that may contribute to APR were found. The allele specific PCR markers developed for these QTL regions would be beneficial for marker-assisted breeding, gene pyramiding, and future cloning of resistance genes from oat.

DeepBSA: A deep-learning algorithm improves bulked segregant analysis for dissecting complex traits

Bulked segregant analysis (BSA) is a rapid, cost-effective method for mapping mutations and quantitative trait loci (QTLs) in animals and plants based on high-throughput sequencing. However, the algorithms currently used for BSA have not been systematically evaluated and are complex and fallible to operate. We developed a BSA method driven by deep learning, DeepBSA, for QTL mapping and functional gene cloning. DeepBSA is compatible with a variable number of bulked pools and performed well with various simulated and real datasets in both animals and plants. DeepBSA outperformed all other algorithms when comparing absolute bias and signal-to-noise ratio. Moreover, we applied DeepBSA to an F₂ segregating maize population of 7160 individuals and uncovered five candidate QTLs, including three well-known plant-height genes. Finally, we developed a user-friendly graphical user interface for DeepBSA, by integrating five widely used BSA algorithms and our two newly developed algorithms, that is easy to operate and can quickly map QTLs and functional genes. The DeepBSA software is freely available to non-commercial users at http://zeasystemsbio.hzau.edu.cn/tools.html and https://github.com/lizhao007/DeepBSA.

A Genomic BSAseq Approach for the Characterization of QTLs Underlying Resistance to Fusarium oxysporum in Eggplant

Eggplant (Solanum melongena L.), similar to many other crops, suffers from soil-borne diseases, including Fusarium oxysporum f. sp. melongenae (Fom), causing wilting and heavy yield loss. To date, the genetic factors underlying plant responses to Fom are not well known. We previously developed a Recombinant Inbred Lines (RILs) population using as a female parent the fully resistant line '305E40' and as a male parent the partially resistant line '67/3'. The fully resistant trait to Fom was introgressed from the allied species S. aethiopicum. In this work, the RIL population was assessed for the responses to Fom and by using a genomic mapping approach, two major QTLs on chromosomes CH02 and CH11 were identified, associated with the full and partial resistance trait to Fom, respectively. A targeted BSAseq procedure in which Illumina reads bulks of RILs grouped according to their resistance score was aligned to the appropriate reference genomes highlighted differentially enriched regions between resistant/susceptible progeny in the genomic regions underlying both QTLs. The characterization of such regions allowed us to identify the most reliable candidate genes for the two resistance traits. With the aim of revealing exclusive species-specific contigs and scaffolds inherited from the allied species and thus associated with the full resistance trait, a draft de-novo assembly of available Illumina sequences of the '305E40' parent was developed to better resolve the non-recombining genomic region on its CH02 carrying the introgressed Fom resistance locus from S. aethiopicum.

Quick mapping and characterization of a co-located kernel length and thousand-kernel weight-related QTL in wheat

A co-located KL and TKW-related QTL with no negative effect on PH and AD was rapidly identified using BSA and wheat 660 K SNP array. Its effect was validated in a panel of 218 wheat accessions. Kernel length (KL) and thousand-kernel weight (TKW) of wheat (Triticum aestivum L.) contribute significantly to kernel yield. In the present study, a recombinant inbred line (RIL) population derived from the cross between the wheat line S849-8 with larger kernels and more spikelets per spike and the line SY95-71 was developed. Further, of both the bulked segregant analysis (BSA) and the wheat 660 K single nucleotide polymorphism (SNP) array were used to rapidly identify genomic regions for kernel-related traits from this RIL population. Kompetitive Allele Specific PCR markers were further developed in the SNP-enriched region on the 2D chromosome to construct a genetic map. Both QKL.sicau-SSY-2D for KL and QTKW.sicau-SSY-2D for TKW were identified at multiple environments on chromosome arm 2DL. These two QTLs explained 9.68-23.02% and 6.73-18.32% of the phenotypic variation, respectively. The effects of this co-located QTL were successfully verified in a natural population consisting of 218 Sichuan wheat accessions. Interestingly, the major QTL was significantly and positively correlated with spike length, but did not negatively affect spikelet number per spike (SNS), plant height, or anthesis date. These results indicated that it is possible to synchronously improve kernel weight and SNS by using this QTL. Additionally, several genes associated with kernel development and filling rate were predicted and sequenced in the QTL-containing physical intervals of reference genomes of 'Chinese spring' and Aegilops tauschii. Collectively, these results provide a QTL with great breeding potential and its linked markers which should be helpful for fine mapping and molecular breeding.

Utility of Triti-Map for bulk-segregated mapping of causal genes and regulatory elements in Triticeae

Triticeae species, including wheat, barley, and rye, are critical for global food security. Mapping agronomically important genes is crucial for elucidating molecular mechanisms and improving crops. However, Triticeae includes many wild relatives with desirable agronomic traits, and frequent introgressions occurred during Triticeae evolution and domestication. Thus, Triticeae genomes are generally large and complex, making the localization of genes or functional elements that control agronomic traits challenging. Here, we developed Triti-Map, which contains a suite of user-friendly computational packages specifically designed and optimized to overcome the obstacles of gene mapping in Triticeae, as well as a web interface integrating multi-omics data from Triticeae for the efficient mining of genes or functional elements that control particular traits. The Triti-Map pipeline accepts both DNA and RNA bulk-segregated sequencing data as well as traditional QTL data as inputs for locating genes and elucidating their functions. We illustrate the usage of Triti-Map with a combination of bulk-segregated ChIP-seq data to detect a wheat disease-resistance gene with its promoter sequence that is absent from the reference genome and clarify its evolutionary process. We hope that Triti-Map will facilitate gene isolation and accelerate Triticeae breeding.

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.

Capturing Wheat Phenotypes at the Genome Level

Recent technological advances in next-generation sequencing (NGS) technologies have dramatically reduced the cost of DNA sequencing, allowing species with large and complex genomes to be sequenced. Although bread wheat (Triticum aestivum L.) is one of the world's most important food crops, efficient exploitation of molecular marker-assisted breeding approaches has lagged behind that achieved in other crop species, due to its large polyploid genome. However, an international public-private effort spanning 9 years reported over 65% draft genome of bread wheat in 2014, and finally, after more than a decade culminated in the release of a gold-standard, fully annotated reference wheat-genome assembly in 2018. Shortly thereafter, in 2020, the genome of assemblies of additional 15 global wheat accessions was released. As a result, wheat has now entered into the pan-genomic era, where basic resources can be efficiently exploited. Wheat genotyping with a few hundred markers has been replaced by genotyping arrays, capable of characterizing hundreds of wheat lines, using thousands of markers, providing fast, relatively inexpensive, and reliable data for exploitation in wheat breeding. These advances have opened up new opportunities for marker-assisted selection (MAS) and genomic selection (GS) in wheat. Herein, we review the advances and perspectives in wheat genetics and genomics, with a focus on key traits, including grain yield, yield-related traits, end-use quality, and resistance to biotic and abiotic stresses. We also focus on reported candidate genes cloned and linked to traits of interest. Furthermore, we report on the improvement in the aforementioned quantitative traits, through the use of (i) clustered regularly interspaced short-palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated gene-editing and (ii) positional cloning methods, and of genomic selection. Finally, we examine the utilization of genomics for the next-generation wheat breeding, providing a practical example of using in silico bioinformatics tools that are based on the wheat reference-genome sequence.

The fruit glossiness locus, dull fruit (D), encodes a C2H2-type zinc finger transcription factor, CsDULL, in cucumber (Cucumis sativus L.)

Fruit glossiness is an important external fruit quality trait for fresh-consumed cucumber fruit, affecting its marketability. Dull fruit appearance is mainly controlled by a single gene, D (for dull fruit) that is dominant to glossy fruit (dd), but the molecular mechanism controlling fruit glossiness is unknown. In the present study, we conducted map-based cloning of the D locus in cucumber and identified a candidate gene (Csa5G577350) that encodes a C₂H₂-type zinc finger transcription factor, CsDULL. A 4895-bp deletion including the complete loss of CsDULL resulted in glossy fruit. CsDULL is highly expressed in the peel of cucumber fruit, and its expression level is positively correlated with the accumulation of cutin and wax in the peel. Through transcriptome analysis, yeast one-hybrid and dual-luciferase assays, we identified two genes potentially targeted by CsDULL for regulation of cutin and wax biosynthesis/transportation that included CsGPAT4 and CsLTPG1. The possibility that CsDULL controls both fruit glossiness and wart development in cucumber is discussed. The present work advances our understanding of regulatory mechanisms of fruit epidermal traits, and provides a useful tool for molecular breeding to improve external fruit quality in cucumber.

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

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