QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations

Genome-wide association study and BSR-seq identify nitrate reductase-related genes in rice landraces (Oryza sativa L.)

Nitrogen (N) is an essential nutrient for rice (Oryza sativa L.) growth and development. However, the lower nitrogen use efficiency (NUE) results in an N fertilizer surplus, which causes many environmental problems. In this study, genome-wide association studies were used to detect nitrate reductase (NR)-related loci in 419 rice landraces. Using the general linear model (GLM), mixed linear model (MLM), linear model (LM), and linear mixed model (LMM), we found six, nine, seven, and six significant single-nucleotide polymorphisms (SNPs) associated (p < 1 × 10-5) for three traits. Moreover, 98 significant SNPs were associated (logarithm of odds ≥ 3) with three traits through 3 V multi-locus random-SNP-effect mixed linear model. Interestingly, we found that Chr1_15896481 was significantly associated in the GLM, MLM, LM, and LMM models. Meanwhile, this significant locus overlapped with a candidate region in bulked segregant RNA sequencing. Through integrated analysis, we identified a most likely candidate genomic region 15,627,420-16,084,761 bp on chromosome 1. By performing functional annotation, RNA sequencing, and real-time quantitative polymerase chain reaction (RT-qPCR) analysis for the genes within this interval, we identified five candidate genes that may affect NR activity. Os01g0378400 exhibits a gene expression pattern highly similar to that of OsNR1.2. It belongs to the NAC transcription factor family, which is involved in plant N metabolism. Os01g0377700 is homologous to an ammonium transporter gene (Cre06g293051). Os01g0383700 encodes a WD40 domain protein, Os01g0379400 encodes an F-box protein, and Os01g0382800 encodes a DYW-type PPR domain protein. These findings will provide valuable genetic resources for NUE genetic improvement in rice breeding.

SKSR1 identified as key virulence factor in Cryptosporidium by genetic crossing

Cryptosporidium is a major cause of severe diarrhea. Although Cryptosporidium isolates exhibit significant differences in infectivity and virulence, the genetic determinants for these traits are not clear. In this study, we use classical genetics to cross two Cryptosporidium parvum isolates of different virulence and use bulk segregant analysis of whole-genome sequences from the progeny to identify quantitative trait loci (QTL) associated with Cryptosporidium infectivity and virulence. Of the 23 genes in three QTL, two have loss-of-function mutations in the low-virulence isolates, including the SKSR1 gene encoding a variant secretory protein. Deletion of the SKSR1 gene or expression of the frame-shifted sequence reduces the pathogenicity of the virulent isolate. SKSR1 is expressed in small granules and secreted into the parasite-host interface during invasion. These results demonstrate that SKSR1 is an important virulence factor in Cryptosporidium, and suggest that the extended SKSR protein family, encoded by clusters of subtelomeric genes, may contribute to pathogenesis.

The knockout of ClaCSLH1 induced dwarfing in watermelon

In agriculture, selecting ideal plant types with desirable traits, such as dwarfing and upright stem structures, significantly enhances crop yield and quality by optimizing light absorption, spatial efficiency, and nutrient utilization. Developing new varieties of dwarf watermelon is a crucial objective in watermelon breeding. In this study, we constructed an F2 population using the wild-type V063 as the paternal parent and the dwarf variety dw-n as the maternal parent. The dwarfing trait was found to be governed by a pair of recessive alleles. Through bulk segregant analysis sequencing (BSA-seq) and RNA sequencing (RNA-seq), we identified the gene Cla97C02G035450, which encodes cellulose synthase-like H1 (CSLH1), as a candidate gene associated with the dwarfing phenotype. ClaCLSH1 belongs to the ClaCESA/CSLs family, which is involved in the cell wall formation by regulating the synthesis of cellulose and hemicellulose. Microscopic analyses revealed that dw-n exhibited shorter internode cells, thicker cell walls, and elevated hemicellulose content compared to V063. Subcellular localization studies demonstrated that the CLACSLH1 protein is primarily localized in the nucleus and the cell membrane/wall. Notably, the overexpression of CLACSLH1 in the dw-n background rescued its dwarf phenotype. Furthermore, experiments indicated that knockdown of CLACSLH1 resulted in excessive hemicellulose synthesis, inhibited internode cell elongation, and ultimately led to the stunted phenotype observed in dw-n. This research provides innovative insights into the development of superior dwarf watermelon varieties and advances our understanding of the molecular mechanisms underlying watermelon dwarfism.

Genetic analysis and fine mapping reveal that AhRt3, which encodes an anthocyanin reductase, is responsible for red testa in cultivated peanuts

KEY MESSAGE

AhRt3, which governs the red testa of peanut, was narrowed down to a 125.30 kb region, and one gene encoding anthocyanin reductase was identified as the putative candidate gene. Testa color is a special characteristic of peanuts (Arachis hypogaea L.), and those with dark testa have been focused on recent years owing to their high-anthocyanin content and increased antioxidant nutritional value. However, the genetic mechanisms underlying this trait remain limited. To identify the gene responsible for the red testa color in peanuts, an F2 population was constructed by crossing YH91 (pink testa) with JHT1 (red testa). Genetic analysis revealed that the red testa was controlled by a single dominant gene named AhRt3 (Arachis hypogaea Red Testa 3). Through bulked segregant analysis sequencing, AhRt3 was preliminarily mapped to the chromosome Arahy.03 and subsequently narrowed to a 125.30 kb genomic region containing 12 potential candidate genes. RNA-seq analysis revealed that 4,880 genes were differentially expressed in the seed testa, with only the candidate gene Arahy.W8TDEC exhibiting higher expression levels in JHT1 than in YH91. Additionally, sequence variation, functional annotation, and expression profiling confirmed that Arahy.W8TDEC, which encodes an anthocyanin reductase, may be a candidate gene for AhRt3. The structural variation involving an inversion between the sixth exon and the 3'UTR of Arahy.W8TDEC resulted in altered amino acids closely associated with the red testa phenotype in peanuts. In conclusion, this study highlights the role of a novel gene in regulating red testa and contributes valuable insights into the genetic basis of seed testa in peanuts.

A feedback loop at the THERMOSENSITIVE PARTHENOCARPY 4 locus controls tomato fruit set under heat stress

High temperatures compromise crop productivity worldwide, but breeding bottlenecks slow the delivery of climate-resilient crops. By investigating tomato fruit set under high temperatures, we discover a module comprising two linked genes, THERMOSENSITIVE PARTHENOCARPY 4a (TSP4a) and TSP4b, which encode the transcriptional regulators IAA9 and AINTEGUMENTA (ANT), respectively, to control thermosensitive parthenocarpy. TSP4a and TSP4b form a positive feedback loop upon heat stress to repress auxin signaling in ovaries. Natural TSP4a and TSP4b alleles bear regulatory-region polymorphisms and are differentially expressed to overcome the trade-off between fruit set and wider plant development. Gene editing of the TSP4a promoter and TSP4b 3' UTR in open-chromatin regions results in expression down-regulation, increased parthenocarpy without yield penalties and maintenance of fruit-sugar levels without broad auxin-related pleiotropic defects in greenhouse-grown plants. These mechanistic insights into heat-induced parthenocarpy and auxin signaling in reproductive organs demonstrate breeding utility to safeguard tomato yield under warming scenarios.

Identification and Characterization of a Male Sterile Rapeseed (Brassica napus) Line for Hybrid Seed Production

A male sterile mutant, S201, was identified in Brassica napus. Genetic analysis revealed that the male sterility trait was controlled by a recessive nuclear gene, male sterility (MS), which was stably inherited. The results of microscopy showed that the main reason for male sterility was a defect in microspore development, resulting in the absence of typical exine and mature microspores. Bulked segregant analysis (BSA) and genotyping of an F2 population showed that the MS gene was located in a 1.4 Mb region. Sequence analysis showed that the CYP704B1 gene in this region contained two non-synonymous SNPs, leading to substitutions of two amino acids. A high-throughput KASP marker was characterized to detect the presence of the ms gene in the breeding population. The data presented here indicate that the male sterile mutant S201 can be applied in rapeseed breeding by producing the male sterile line and that the KASP marker developed for male sterility will be useful in marker-assisted selection of male sterile individuals in rapeseed-breeding programs.

Identification and characterization of the powdery mildew resistance in cultivated emmer wheat accession Lxd-682 via bulked segregant RNA sequencing

BACKGROUND

Common wheat (Triticum aestivum L.) is a vital source of nutrition for human consumption. However, wheat production is significantly threatened by various diseases, such as powdery mildew, a widespread fungal disease caused by Blumeria graminis f. sp. tritici (Bgt). Utilizing and identifying resistance genes and elucidating the molecular mechanisms underlying this resistance are the most effective and sustainable ways to fight this disease.

RESULTS

Lxd-682, a cultivated emmer wheat accession, exhibited resistance to 12 out of 13 tested Bgt isolates at the seedling stage. Genetic analysis revealed that this resistance is conferred by a single dominant gene, tentatively designated as PmLxd-682. Molecular mapping positioned PmLxd-682 between the markers WGRE77413 and WGRC1096, with the Pm4-diagnostic marker JS717/JS718 co-segregating. Homology-based cloning and sequence alignment further confirmed that PmLxd-682 is identical to Pm4a. qRT-PCR analysis showed that the alternative splicing PmLxd-682-V2 exhibited higher expression level than that of PmLxd-682-V1 post-Bgt invasion, suggesting its prominent role in fighting Bgt invasion. Additionally, four pathogenesis-related (PR) genes were significantly up-regulated in both Lxd-682 and susceptible parent Langdon upon infection, revealing possibly unimportant roles in resistance pathway. Furthermore, 1,567 differentially expressed genes (DEGs) between resistant and susceptible bulks were identified through BSR-Seq, with 490 ones located within the candidate interval on chromosome 2AL, and potential biological processes associated with resistance were enriched via gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathway analysis. To verify the potential regulatory genes, three key genes, TRITD2 Av1G294940, TRITD2 Av1G036490 and TRITD2 Av1G295220 all encoding disease resistance protein, were selected from six candidates via qRT-PCR following post-Bgt invasion. Molecular markers JS717/JS718 and WGRC1096 were confirmed to be available for marker-assisted selection (MAS) of PmLxd-682 in breeding practices.

CONCLUSIONS

The study identified key genetic intervals and genes involved in the resistance of a cultivated emmer wheat accession Lxd-682 to powdery mildew. These findings significantly advance our understanding of plant-pathogen interactions and establish a solid foundation for future genetic and functional studies aimed at improving disease resistance in crops.

Mapping of resistance genes to powdery mildew based on DNA re-sequencing and bulk segregant analysis in Capsicum

Powdery mildew caused by Leveillula taurica adversely affects the development and growth of pepper plants. However, there have been few reports on the fine mapping and quantitative trait locus (QTLs) gene cloning of resistance genes to powdery mildew in pepper. Herein, an F2 segregating population was constructed using the high resistance material "NuMex Suave Red" and the extremely susceptible material "c89" for bulked segregant analysis and DNA re-sequencing (BSA-seq). Molecular markers were used to achieve fine mapping, followed by expression verification. A major QTL located on chromosome 5 (Chr5, 7.20-11.75 Mb) that is associated with resistance to powdery mildew in pepper was mapped using BSA-seq. A narrow interval of 64.86 kb encompassing five genes was refined using InDel and KSAP molecular markers developed from the QTL region. Among them, the expression of the ubiquitin-conjugating enzyme E2 gene, Capana05g000392, was significantly upregulated in multiple resistant materials. In addition, there was a single nucleotide polymorphism (SNP) of A/G in the 241st position of the CDS sequence of Capana05g000392, which in turn leads to an amino acid polymorphism of M/V between susceptible parent and resistant parent. Overall, these results indicate that the Capana05g000392 gene may serve as a robust potential factor against powdery mildew in pepper. These findings further elucidate the genetic mechanism of resistance to powdery mildew in pepper and facilitate molecular marker-assisted breeding.

Manipulation of the brown glume and internode 1 gene leads to alterations in the colouration of lignified tissues, lignin content and pathogen resistance in wheat

Lignin is a crucial component of the cell wall, providing mechanical support and protection against biotic and abiotic stresses. However, little is known about wheat lignin-related mutants and their roles in pathogen defence. Here, we identified an ethyl methanesulfonate (EMS)-derived Aegilops tauschii mutant named brown glume and internode 1 (bgi1), which exhibits reddish-brown pigmentation in various tissues, including internodes, spikes and glumes. Using map-based cloning and single nucleotide polymorphism (SNP) analysis, we identified AET6Gv20438400 (BGI1) as the leading candidate gene, encoding the TaCAD1 protein. The mutation occurred in the splice acceptor site of the first intron, resulting in a premature stop codon in BGI1. We validated the function of BGI1 using loss-of-function EMS and gene editing knockout mutants, both of which displayed reddish-brown pigmentation in lignified tissues. BGI1 knockout mutants exhibited reduced lignin content and shearing force relative to wild type, while BGI1 overexpression transgenic plants showed increased lignin content and enhanced disease resistance against common root rot and Fusarium crown rot. We confirmed that BGI1 exhibits CAD activity both in vitro and in vivo, playing an important role in lignin biosynthesis. BGI1 was highly expressed in the stem and spike, with its localisation observed in the cytoplasm. Transcriptome analysis revealed the regulatory networks associated with BGI1. Finally, we demonstrated that BGI1 interacts with TaPYL-1D, potentially involved in the abscisic acid signalling pathway. The identification and functional characterisation of BGI1 significantly advance our understanding of CAD proteins in lignin biosynthesis and plant defence against pathogen infection in wheat.

The Genome of Lolium multiflorum Reveals the Genetic Architecture of Paraquat Resistance

Herbicide resistance in agricultural weeds has become one of the greatest challenges for sustainable crop production. The repeated evolution of herbicide resistance provides an excellent opportunity to study the genetic and physiological basis of the resistance phenotype and the evolutionary responses to human-mediated selection pressures. Lolium multiflorum is a ubiquitous weed that has evolved herbicide resistance repeatedly around the world in various cropping systems. We assembled and annotated a chromosome-scale genome for L. multiflorum and elucidated the genetic architecture of paraquat resistance by performing quantitative trait locus analysis, genome-wide association studies, genetic divergence analysis and transcriptome analyses from paraquat-resistant and -susceptible L. multiflorum plants. We identified two regions on chromosome 5 that were associated with paraquat resistance. These regions both showed evidence for positive selection among the resistant populations we sampled, but the effects of this selection on the genome differed, implying a complex evolutionary history. In addition, these regions contained candidate genes that encoded cellular transport functions, including a novel multidrug and toxin extrusion (MATE) protein and a cation transporter previously shown to interact with polyamines. Given that L. multiflorum is a weed and a cultivated crop species, the genomic resources generated will prove valuable to a wide spectrum of the plant science community. Our work contributes to a growing body of knowledge on the underlying evolutionary and ecological dynamics of rapid adaptation to strong anthropogenic selection pressure that could help initiate efforts to improve weed management practices in the long term for a more sustainable agriculture.

Regulation of maize kernel development via divergent activation of α-Zein genes by transcription factors O11, O2, and PBF1

α-Zeins, the major maize endosperm storage proteins, are transcriptionally regulated by Opaque 2 (O2) and PROLAMIN-BOX BINDING FACTOR1 (PBF1), with Opaque 11 (O11) functioning upstream of them. However, whether O11 directly binds to α-zein genes and its regulatory interactions with O2 and PBF1 remain unclear. Using the small-kernel mutant sw1, which exhibits decreased 19-kDa and increased 22-kDa α-zein, we positionally cloned O11 and found it directly binds to G-box/E-box motifs. O11 activates 19-kDa α-zein transcription, stronger than PBF1 but weaker than O2. Notably, PBF1 competitively binds to overlapping E-box/P-box motif, and represses O11-mediated transactivation. Although O11 does not physically interact with O2, it participates in the O2-centered hierarchical network to enhance α-zein expression. sw1 o2 and sw1 pbf1 double mutants exhibit smaller, more opaque kernels with further reduced 19-kDa and 22-kDa α-zeins compared to the single mutants, suggesting distinct regulatory effects of these transcription factors on 19-kDa and 22-kDa α-zein genes. Promoter motif analysis suggests that O11, PBF1, and O2 directly regulate 19-kDa α-zein genes, while O11 indirectly controls 22-kDa α-zein genes via O2 and PBF1 modulation. These findings identify the unique and coordinated roles of O11, O2, and PBF1 in regulating α-zein genes and kernel development.

A 294 kb deletion causes reduced leaflet size and biomass in pigeonpea

KEY MESSAGE

BSA-seq and fine mapping revealed a 294 kb deletion on chromosome 9 regulating leaflet size and biomass in pigeonpea. Leaf size critically influences photosynthetic capacity, impacting organic matter production and biomass yield. This study reports the identification and characterization of a small leaflet mutant (sl1) in pigeonpea (Cajanus cajan) generated via aerial mutagenesis. Compared to the wild-type Qiongzhong, sl1 displayed significantly reduced leaf area, plant height, stem diameter, and biomass, characteristic of a dwarf phenotype. Genetic analysis confirmed a single recessive locus controlling the sl1 phenotype. Bulked segregant analysis sequencing (BSA-seq) and fine mapping identified the causal mutation as a 294 kb deletion encompassing 21 genes on chromosome 9. Transcriptomic analysis identified 1,039 differentially expressed genes (DEGs), indicating disruptions in, among others, plant hormone signaling pathways. Analysis of 28 target plant hormone metabolites revealed significant shifts in sl1 mutant compared to wild-type, including increased levels of strigolactone, methyl indole-3-acetate, and trans-zeatin-riboside, and decreases in gibberellin A3, N6-isopentenyladenine, and methyl jasmonate. Cytological analysis revealed a decreased cell number in sl1 leaves, contributing to the reduced leaflet size. Three candidate genes, CC09g01700, CC09g01705, and CC09g01707, within the deleted region were prioritized based on their altered expression patterns and their putative roles in leaf development. These findings elucidate the genetic regulation of leaf morphology and biomass in pigeonpea, offering potential targets for marker-assisted selection to enhance pigeonpea yield.

Chromosomal structural variation loci HSS1 and HSS6 lead to hybrid sterility in rice

KEY MESSAGE

Chromosomal structural variation leads to rice hybrid sterility, and this study has established a new inheritance model and offers a novel insight into hybrid sterility in rice. The utilization of heterosis among different species significantly enhances rice yield. However, the direct application of this advantage is hindered by hybrid sterility (HS). Here we identify a stable semi-sterile inbred line (SSIL) in the hybridization between rice species O. sativa indica and wild rice O. longistaminata. Both pollen and spikelet fertility in heterozygous SSIL plants are semi-sterile, and the homozygous plants show normal fertility. Interestingly, unlike previously reported hybrid sterility loci, SSIL does not induce segregation distortion in progeny. The genotypes of chromosome 1 (Chr.1) exhibit tight pseudo-linkage on Chr.6 in SSIL population. HSS1 and HSS6 are named as the two loci of Chr.1 and Chr.6. Cytological observations revealed abnormalities in male and female gametes during meiosis, ultimately resulting in semi-sterile pollen and spikelet. HSS1 was narrowed down to a 190.6-kb interval and HSS6 to a 1391.6-kb region using an SSIL population of 10,393 plants. Through high-throughput sequencing and observation of chromosomal behavior during meiosis, a reciprocal translocation between the short arm of Chr.1 and the end of Chr.6 in the translocated RD23-type (RT-type) gamete in SSIL was discovered. The presence of RT-type gametes in SSIL is the direct causative factor for semi-sterility in both male and female gametes. This translocation led to abnormal SSIL-RT/L synapses and the formation of tetravalent ring structures in chromosomes during pachynema. Our findings have uncovered a pair of reciprocal chromosome translocations that control reproductive isolation in rice. These insights offer valuable guidance for optimizing hybrid breeding applications, ultimately enhancing the benefits of heterosis.

Fine mapping and candidate gene analysis of major QTLs for number of seeds per pod in Arachis hypogaea L

BACKGROUND

Peanut (Arachis hypogaea L., 2n = 2x = 20) is an important industrial and oil crop that is widely grown in more than 100 countries. In recent years, breeders have focused on increasing the seed number per pod to improve their yield in addition to other breeding for other key components of yield, including the pod number, seeds per pod, and 100-seed weight.

RESULTS

In this study, a secondary population of 1,114 BC1F2 lines was derived from a cross between the parents R45 and JNH3. Two stable major-effect quantitative trait loci of qRMPA09.1 and qRMPA09.2 were detected simultaneously and mapped within chromosomal intervals of approximately 400 Kb and 600 Kb on chromosome A09. Additionally, combined whole-genome and RNA sequencing analyses showed the differential expression of the Arahy.04JNDX gene that belongs to a MYB transcription factor (TF) between the two parents. The AhMYB51 gene was also inferred to influence the number of seeds per pod in peanuts. An examination of the backcross lines L2/L4 showed that AhMYB51 increases the rate of multiple pods per plant (RMSP) primarily by affecting brassinosteroids in the flowers, while its overexpression promotes the length of siliques in Arabidopsis thaliana.

CONCLUSIONS

Our findings provide valuable insights for the cloning of favorable alleles for RMSP in peanuts. The qRMSPA09.1 and qRMSPA09.2 are two novel QTL associated with the RMSP trait, with AhMYB51 predicted as its candidate gene. Moreover, the closely linked polymorphic SNP markers for loci of two significant QTLs may be useful in accelerating marker-assisted breeding in peanuts.

Fine mapping and identification of CqMYB62 as the subgynoecy gene in chieh-qua (Benincasa hispida Cogn. var. Chieh-qua How)

KEY MESSAGE

It was hypothesized that Bch08G003160 (CqMYB62), located in the 51.08 Kb region on chromosome 08, might be an important candidate gene for the subgynoecy trait in chieh-qua, based on BSA-seq and linkage mapping approaches. In cucurbit crops, the use of female lines can greatly increase the yield of a single plant and is especially important for the production of hybrid generation seeds, thus being of great interest to breeders. To identify genes regulating sex differentiation in chieh-qua, genetic analysis of the subgynoecy trait was conducted using a chieh-qua F2 population. Initial localization of the locus was done using BSA-seq, followed by fine mapping with a large F2 population (n = 2,741). The locus was ultimately narrowed down to a 51.08 Kb region on Chr08, revealing a single gene Bch08G003160 (CqMYB62) in this region. Further analysis revealed that the presence of two variant loci (SNP_416 and SNP_317) in the coding region resulted in premature termination of the codon and loss of function of the structural domain of the PLN03212 superfamily. Moreover, our research indicated that the subgynoecy trait mediated by CqMYB62 in chieh-qua is potentially regulated by gibberellic acid (GA). Two efficient dCAPS markers were developed to distinguish subgynoecy. In summary, these findings highlight the critical role of CqMYB62 in subgynoecy trait regulation, offering potential implications for chieh-qua breeding programs.

A genomic variation map provides insights into potato evolution and key agronomic traits

Hybrid potato breeding based on diploid inbred lines is transforming the way of genetic improvement of this staple food crop, which requires a deep understanding of potato domestication and differentiation. In the present study, we resequenced 314 diploid wild and landrace accessions to generate a variome map of 47,203,407 variants. Using the variome map, we discovered the reshaping of tuber transcriptome during potato domestication, characterized genome-wide differentiation between landrace groups Stenotomum and Phureja. We identified a jasmonic acid biosynthetic gene possibly affecting the tuber dormancy period. Genome-wide association studies revealed a UDP-glycosyltransferase gene for the biosynthesis of anti-nutritional steroidal glycoalkaloids (SGAs), and a Dehydration Responsive Element Binding (DREB) transcription factor conferring increased average tuber weight. In addition, genome similarity and group-specific SNP analyses indicated that tetraploid potatoes originated from the diploid Solanum tuberosum group Stenotomum. These findings shed light on the evolutionary trajectory of potato domestication and improvement, providing a solid foundation for advancing hybrid potato-breeding practices.

Mapping and molecular marker development for the BnaSBT gene controlling inflorescence and plant architectures in B. napus

Exploring the molecular mechanism underlying plant architecture and breeding new varieties suitable for mechanized harvesting are primary objectives for rapeseed breeders in China. However, few genes controlling plant architecture have been cloned in Brassica napus. In this study, SX3, a scattered-bud B. napus line with a dwarf and compact plant architecture, was characterized. To identify the genes underlying bud arrangement, plant height and branch angle, segregating populations were constructed by crossing SX3 with two clustered-bud lines with a tall and loose plant architecture. Genetic analysis revealed that the scattered-bud trait (SBT) was controlled by a single dominant gene, BnaSBT. BnaSBT is likely a pleiotropic gene that simultaneously controls plant height and branch angle. Using BSA-seq analysis, BnaSBT was mapped to a 4.15 Mb region on ChrA10. Owing to the lack of recombinants within this region, it was infeasible to finely map BnaSBT. RNA-seq analysis of BC2 plants with contrasting inflorescence and plant architectures revealed that the upregulation of genes involved in amino acid and lipid metabolism and genes encoding MADS-box transcription factors is related to the the phenotype of SX3. These findings together with comparative sequencing indicated that BnaA10.SEP1, BnaA10.AGL15, BnaA10.GLN1-4 and BnaA10.AGP15 are candidate genes for BnaSBT. Markers closely linked to the scattered-bud trait were developed for selecting dwarf and compact plants. These findings provide molecular markers and germplasms for breeding new varieties with ideal plant types and lay a theoretical foundation for cloning key genes and elucidating the genetic basis of inflorescence and plant architectures in B. napus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-025-01556-2.

A premature termination codon mutation in the onion AcCER2 gene is associated with both glossy leaves and thrip resistance

Plant epicuticular waxes (EW) play a critical role in defending against biotic and abiotic stresses. Notably, onions (Allium cepa L.) present a distinctive case where the mutant with defect in leaf and stalk EW showed resistance to thrips compared with the wild type with integral EW. We identified a premature stop codon mutation in the AcCER2 gene, an ortholog of CER2 gene in Arabidopsis thaliana that has been proved essential for the biosynthesis of very long-chain fatty acids (VLCFAs), in the onions with glossy leaf and stalks in our experiments. The data hinted at the possibility that this mutation might impede the elongation process of VLCFAs from C28 to C32, thereby hindering the production of 16-hentriacontanone, a primary constituent of onion EW. Transcriptomic analysis revealed substantial alterations in expression of genes in the pathways related not only to lipid synthesis and transport but also to signal transduction and cell wall modification in glossy mutants. Meanwhile, metabolomic profiling indicates a remarkable increase in flavonoid accumulation and a significant reduction in soluble sugar content in glossy mutants. These findings suggested that the enhanced resistance of glossy mutants to thrips might be a consequence of multiple physiological changes, and our integrated multiomics analysis highlighting the regulatory role of AcCER2 in these processes. Our study has yielded valuable insights into the biosynthesis of onion EW and has provided an initial hypothesis for the mechanisms underlying thrip resistance. These findings hold significant promise for the breeding programs of thrip-resistant onion.

Solanum pan-genetics reveals paralogues as contingencies in crop engineering

Pan-genomics and genome-editing technologies are revolutionizing breeding of global crops1,2. A transformative opportunity lies in exchanging genotype-to-phenotype knowledge between major crops (that is, those cultivated globally) and indigenous crops (that is, those locally cultivated within a circumscribed area)3-5 to enhance our food system. However, species-specific genetic variants and their interactions with desirable natural or engineered mutations pose barriers to achieving predictable phenotypic effects, even between related crops6,7. Here, by establishing a pan-genome of the crop-rich genus Solanum8 and integrating functional genomics and pan-genetics, we show that gene duplication and subsequent paralogue diversification are major obstacles to genotype-to-phenotype predictability. Despite broad conservation of gene macrosynteny among chromosome-scale references for 22 species, including 13 indigenous crops, thousands of gene duplications, particularly within key domestication gene families, exhibited dynamic trajectories in sequence, expression and function. By augmenting our pan-genome with African eggplant cultivars9 and applying quantitative genetics and genome editing, we dissected an intricate history of paralogue evolution affecting fruit size. The loss of a redundant paralogue of the classical fruit size regulator CLAVATA3 (CLV3)10,11 was compensated by a lineage-specific tandem duplication. Subsequent pseudogenization of the derived copy, followed by a large cultivar-specific deletion, created a single fused CLV3 allele that modulates fruit organ number alongside an enzymatic gene controlling the same trait. Our findings demonstrate that paralogue diversifications over short timescales are underexplored contingencies in trait evolvability. Exposing and navigating these contingencies is crucial for translating genotype-to-phenotype relationships across species.

Genetic, metabolomic and transcriptomic analyses of the cotton yellow anther trait

In the fiber industry, cotton (Gossypium hirsutum L.) is an important crop. One of the most important morphology traits of plants is the color of the anthers, is closely related to pollen fertility and stress resistance. Upland cotton anthers appear white, while island cotton and many wild cotton species have yellow anthers. Carotenoids are natural pigments in plants which involved in many metabolic processes, including photosynthesis, photoprotection, photomorphogenesis, growth and development. Here, we characterized the yellow anther trait of G. hirsutum. Carotenoid and flavonoid profiles in the yellow anthers were greatly altered compared to that in the white anthers, indicating that both carotenoids and flavonoids contribute to the yellow anther phenotype. Map-based cloning identified GhYA (GH_A05G4013) encoding a phytoene synthase to be the candidate gene responsible for anther coloration. GhYA is predominantly expressed in anthers, with its expression level gradually decreasing with the development of anthers. Haplotype analysis revealed that white anthers are associated with two haplotypes, with X74 belonging to HAP1. Through evolutionary analysis, it was found that although there are many white anther Germplasm in upland cotton, the two types of white anther haplotypes were mutated from yellow anthers respectively. Comparative transcriptome analysis between the yellow anther and white anther accessions revealed differentially expressed genes related to both the carotenoid and flavonoid biosynthesis pathways, in line with the changed profiles of the two types of metabolites in yellow anthers; meanwhile, it also indicates potential cross-talk between the flavonoid and carotenoid pathways. According to the results, the PSY gene is critical for the regulation of carotenoids accumulation in cotton anthers.

Transcriptional Reprogramming and Allelic Variation in Pleiotropic QTL Regulates Days to Flowering and Growth Habit in Pigeonpea

The present study investigated the linkage between days to flowering (DTF) and growth habit (GH) in pigeonpea using QTL mapping, QTL-seq, and GWAS approaches. The linkage map developed here is the largest to date, spanning 1825.56 cM with 7987 SNP markers. In total, eight and four QTLs were mapped for DTF and GH, respectively, harbouring 78 pigeonpea orthologs of Arabidopsis flowering time genes. Corroboratively, QTL-seq analysis identified a single linked QTL for both traits on chromosome 3, possessing 15 genes bearing genic variants. Together, these 91 genes were clustered primarily into autonomous, photoperiod, and epigenetic pathways. Further, we identified 39 associations for DTF and 111 associations for GH through GWAS in the QTL regions. Of these, nine associations were consistent and constituted nine haplotypes (five late, two early, one each for super-early and medium duration). The involvement of multiple genes explained the range of allelic effects and the presence of multiple LD blocks. Further, the linked QTL on chromosome 3 was fine-mapped to the 0.24-Mb region with an LOD score of 8.56, explaining 36.47% of the phenotypic variance. We identified a 10-bp deletion in the first exon of TFL1 gene of the ICPL 20338 variety, which may affect its interaction with the Apetala1 and Leafy genes, resulting in determinate GH and early flowering. Further, the genic marker developed for the deletion in the TFL1 gene could be utilized as a foreground marker in marker-assisted breeding programmes to develop early-flowering pigeonpea varieties.

Regreening mechanisms in cucumber: insights from a CsSIG2 mutation affecting chloroplast development

KEY MESSAGE

CsSIG2 is essential for cucumber chloroplast development, and mutations in CsSIG2 reveal mechanisms that restore chloroplast functionality and drive the regreening phenotype in the mutant. Chloroplast development and leaf color are essential traits that significantly influence plant photosynthesis and overall vigor. This study investigates a natural mutation in the cucumber that leads to a virescent leaf-color (Csvl-6) phenotype characterized by an initial yellow color in cotyledons and young leaves, which gradually transition to green as the plant matures. We utilized bulked segregant analysis and genetic linkage mapping to locate the best candidate gene sigma factor 2 (CsSIG2) on chromosome 6, identifying a single nonsynonymous SNP resulting in an arginine to glycine substitution in the CsSIG2 protein. Comparative transcriptome analysis highlighted that this mutation disrupts early chloroplast biogenesis and delays chlorophyll accumulation, but the chloroplasts can recover, leading to greening during later stages of leaf development. Our findings reveal that the recovery phenomenon involves upregulation of chloroplast-encoded genes responsible for thylakoid membrane formation and photosystem function, alongside altered expression of transcription factors linked to chlorophyll metabolism. This study elucidates the genetic and molecular basis of chloroplast development in cucumber, providing valuable insights into the mechanisms underlying leaf greening, which could inform future breeding efforts focused on manipulating leaf color traits for enhanced crop performance.

A special short-wing petal faba genome and genetic dissection of floral and yield-related traits accelerate breeding and improvement of faba bean

BACKGROUND

A comprehensive study of the genome and genetics of superior germplasms is fundamental for crop improvement. As a widely adapted protein crop with high yield potential, the improvement in breeding and development of the seeds industry of faba bean have been greatly hindered by its giant genome size and high outcrossing rate.

RESULTS

To fully explore the genomic diversity and genetic basis of important agronomic traits, we first generate a de novo genome assembly and perform annotation of a special short-wing petal faba bean germplasm (VF8137) exhibiting a low outcrossing rate. Comparative genome and pan-genome analyses reveal the genome evolution characteristics and unique pan-genes among the three different faba bean genomes. In addition, the genome diversity of 558 accessions of faba bean germplasm reveals three distinct genetic groups and remarkable genetic differences between the southern and northern germplasms. Genome-wide association analysis identifies several candidate genes associated with adaptation- and yield-related traits. We also identify one candidate gene related to short-wing petals by combining quantitative trait locus mapping and bulked segregant analysis. We further elucidate its function through multiple lines of evidence from functional annotation, sequence variation, expression differences, and protein structure variation.

CONCLUSIONS

Our study provides new insights into the genome evolution of Leguminosae and the genomic diversity of faba bean. It offers valuable genomic and genetic resources for breeding and improvement of faba bean.

A glycogen synthase kinase-3 gene enhances grain yield heterosis in semi-dwarf rapeseed

Optimizing plant height is a key breeding objective in Brassica napus to enhance lodging resistance and increase yield potential. In the present study, we identified a semi-dwarf gene in rapeseed, BnDWARF5 (BnDF5), which encodes a glycogen synthase kinase 3, BRASSINOSTEROID-INSENSITIVE 2 (BnaC03.BIN2), primarily controlling the elongation of basal internodes by inhibiting the elongation of internode cells. Genetic mapping and cloning revealed that BnDF5 is governed by a semi-dominant/dominant gene located on chromosome C03. Sequencing uncovered an SNP in BnaC03.BIN2 due to an amino acid substitution, which was confirmed via kompetitive allele-specific polymerase chain reaction marker analysis, and expressing the mutated BnaC03.BIN2 in the wild type resulted in decreased plant height. Practical breeding applications showed that heterozygous BnDF5 plants exhibited optimal intermediate height and strong yield heterosis, making the semi-dwarf mutant a valuable genetic resource for developing semi-dwarf rapeseed varieties with improved lodging resistance and yield.

Mapping and characterization of a novel powdery mildew resistance locus (PM2) in Cannabis sativa L

Introduction

Breeding genetic resistance to economically important crop diseases is the most sustainable strategy for disease management and enhancing agricultural and horticultural productivity, particularly where the application of synthetic pesticides is prohibited. Powdery mildew disease, caused by the biotrophic fungal pathogen Golovinomyces ambrosiae, is one of the most prevalent threats to the cannabis and hemp industry worldwide.

Methods

In this study, we used bulked-segregant analysis combined with high-throughput RNA sequencing (BSRSeq) to identify and map a novel single dominant resistance (R) locus (designated PM2), that strongly suppresses powdery mildew infection and sporulation in Cannabis sativa.

Results and discussion

BSA mapped PM2 to chromosome 9. Histochemical analysis revealed that PM2-induced resistance is mediated by a highly localized hypersensitive response mainly in the epidermal cells of the host. Importantly, genetic markers capable of tracking PM2 resistance in breeding populations were developed using associated SNPs identified in this study. The ability to track PM2 will allow for successful introgression of PM resistance into elite cannabis cultivars and help move towards a more sustainable cannabis industry.

Meta genetic analysis of melon sweetness

KEY MESSAGE

Through meta-genetic analysis of Cucumis melo sweetness, we expand the description of the complex genetic architecture of this trait. Integration of extensive new results with published QTL data provides an outline towards construction of a melon sweetness pan-QTLome. An ultimate objective in crop genetics is describing the complete repertoire of genes and alleles that shape the phenotypic variation of a quantitative trait within a species. Flesh sweetness is a primary determinant of fruit quality and consumer acceptance of melons. Cucumis melo is a diverse species that, among other traits, displays extensive variation in total soluble solids (TSS) content in fruit flesh, ranging from 20 Brix in non-sweet to 180 Brix in sweet accessions. We present here meta-genetic analysis of TSS and sugar variation in melon, using six different populations and fruit measurements collected from more than 30,000 open-field and greenhouse-grown plants, integrated with 15 published melon sweetness-related quantitative trait loci (QTL) studies. Starting with characterization of sugar composition variation across 180 diverse accessions that represent 3 subspecies and 12 of their cultivar-groups, we mapped TSS and sugar QTLs, and confirmed that sucrose accumulation is the key variable explaining TSS variation. All modes-of-inheritance for TSS were displayed by multi-season analysis of a broad half-diallel population derived from 20 diverse founders, with significant prevalence of the additive component. Through parallel genetic mapping in four advanced bi-parental populations, we identified common as well as unique TSS QTLs in 12 chromosomal regions. We demonstrate the cumulative less-than-additive nature of favorable TSS QTL alleles and the potential of a QTL-stacking approach. Using our broad dataset, we were additionally able to show that TSS variation displays weak genetic correlations with melon fruit size and ripening behavior, supporting effective breeding for sweetness per se. Our integrated analysis, combined with additional layers of published QTL data, broadens the perspective on the complex genetic landscape of melon sweetness and proposes a scheme towards future construction of a crop community-driven melon sweetness pan-QTLome.

Molecular mapping and development of SSR markers associated with Chilli leaf curl virus resistance in chilli (Capsicum annuum L.)

Chilli leaf curl virus (ChLCV), caused by whiteflies transmitted begomoviruses, poses a significant threat to chilli cultivation and production all over the world. Exploring host plant resistance and identification of associated molecular markers will accelerate resistance breeding. QTL-seq analysis was employed in the IHR4615(R) × IHR2451(S) F2 population to identify QTLs associated with ChLCV-Raichur isolate resistance. A significant QTL was mapped on chromosome 6 associated with ChLCV-Raichur isolate resistance. Defense-related genes were predicted within the identified locus. Further refining of the identified locus with simple sequence repeats (SSR) markers led to the identification of two SSR markers IHR-LCV-SSR-76 and IHR-LCV-SSR-165 associated with the ChLCV-Raichur isolate resistance with 89.50 and 72.50 % prediction efficacy, respectively in IHR4615(R) × IHR2451(S) F2 population. These markers are located at 7 cM and 17.65 cM genetic distances from the resistant gene. These markers were further validated in another resistant source-based F2 population of IHR4392(S) × IHR4597(R). The developed and validated molecular markers can be explored in marker-assisted breeding programs aiming at developing resistant cultivars/ F1 hybrids of chilli.

Assessment of potential candidate genes for partial resistance to Sclerotinia stem rot caused by Sclerotinia sclerotiorum using real-time quantitative PCR

Sclerotinia stem rot (SSR), caused by Sclerotinia sclerotiorum (Lib) de Bary (S. sclerotiorum), is one of the most important diseases that causes significant soybean [Glycine max (L.) Merr.] seed yield and quality losses in Canada and globally. Initiation of plant defense mechanisms is crucial for establishing partial resistance to the pathogenic fungus. To understand plant response to S. sclerotiorum, we conducted a temporal (1, 3, and 5 days post-inoculation [DPI]) assessment of gene expression changes in the stem of soybean genotypes with contrasting phenotypic response. We focused on four genes that have been previously reported as associated with SSR partial resistance and are known to be involved in defense-related functions such as cell wall modification, signaling, response to wounding, and response to fungus. The results showed a higher and earlier expression of the genes in partially resistant cultivars compared to the susceptible. Expression of some genes increased up to 11- (Glyma.02G059700) to 16-fold (Glyma.09G232100) by 3 DPI in the partially resistant cultivar, OAC Drayton, while the genes were generally downregulated in the susceptible cultivar, OAC Shire, at the same DPI. This study improves our understanding of expression patterns of genes involved in plant defense against fungal pathogens in soybean. More importantly, the knowledge of genes that are essential in defense against S. sclerotiorum can be used to fine-map the quantitative trait loci for SSR resistance and facilitate accelerated breeding of SSR-resistant cultivars through gene-based marker-assisted selection.

Integrated GWAS, BSA-seq, and RNA-seq analyses to identify candidate genes associated with male fertility trait in peach

Male sterility in peach (Prunus persica L.), characterized by the absence of fertile pollen grains in the anther, is determined by a recessive allele in homozygosis of the major gene located on chromosome 6. Developing tightly linked molecular markers can help identify appropriate peach parents or male-sterile plants for early culling in segregating progenies, thereby increasing breeding efficiency. In this study, we performed comprehensive research integrating genome-wide association study, bulked segregant analysis, and tissue-specific transcriptome sequencing for precisely characterizing the genes associated with male sterility and fertility in peach. We identified the candidate gene Prupe.6G027000, which encodes an ATP-binding cassette transporter G family member 26 (ABCG26), as a reliable candidate for controlling the targeted traits, as indicated by gene expression profiling and validated by quantitative real-time polymerase chain reaction, in situ hybridization, and virus-induced gene silencing. Prupe.6G027000 was transcribed preferentially on the tapetum and microspore surface, and its transient silencing caused severe pollen abortion in peach. The genotypes of nonsynonymous single-nucleotide variation (T > C) harbored in the coding region of Prupe.6G027000 exhibited approximately 96.2% consistency with male fertile or sterile phenotype in 579 peach accessions. These findings lay the foundation for dissecting the genetic basis of male fertility traits, and facilitating the establishment of a marker-assisted selection system in peaches.

The allelic mutation of NBS-LRR gene causes premature senescence in wheat

Premature senescence has a significant impact on the yield and quality of wheat crops. The process is controlled by multiple and intricate genetic pathways and regulatory elements, whereby the discovery of additional mutants provides important insights into the molecular basis of this important trait. Here, we developed a premature senescence wheat mutant je0874, its leaves started to show yellow before heading stage; with plant growth and development, the degree of yellowing worsened rapidly, and chlorophyll content in flag leaf was reduced by 93.8 % at 15 days after heading, all other leaves became dryness at the grain filling stage. In the mutant, the reactive oxygen species (ROS) and its metabolites increased up to 34.8-47.3 %, while activities of ROS scavenging enzymes were reduced by 62.7-96.7 %. Premature senescence resulted in a reduction of thousand grain weight by over 50 %. Genetic analysis showed the mutation of senescence was controlled by a single recessive gene, and target gene was finely mapped to a 338 kb region of the long arm of chromosome 2D. This region contained a total of 6 annotated genes, while only gene TraesFLD2D01G513900 carried a SNP mutation. The gene contained an NBS-LRR domain, we named it Taps1. Allelic mutants of Taps1 exhibited a lesion mimic phenotype, and the mutant allele resulted in cell death in tobacco, which represent a novel gene controlling wheat senescence. Two haplotypes were identified in 180 accessions, which did not lead to cell death. These results contribute to increase our understanding of the regulation of premature plant senescence.

Mapping resistance to Sclerotinia white mold in two pinto bean recombinant inbred line populations

White mold, caused by the fungus Sclerotinia sclerotiorum (Lib.) de Bary, is a devastating disease affecting common bean (Phaseolus vulgaris L.) production worldwide. Breeding for resistance to white mold is challenging due to its quantitative inheritance and intricate genetic mechanisms. This research aimed to validate and characterize physiological resistance in the pinto dry bean market class through greenhouse straw tests under controlled conditions and field assessments under natural environments. Classical quantitative trait locus (QTL) mapping and Khufu de novo QTL-seq were employed to detect and narrow QTL intervals and identify candidate genes associated with white mold resistance in two pinto bean recombinant inbred line populations, PT9-5-6/USPT-WM-12 (P2) and PT12-37/VCP-13 (P3). Eleven QTL, five in P2 and six in P3, conditioning white mold resistance were identified. New QTL were discovered including WM1.4 and WM11.5 in P2, and WM1.5 and WM7.7 in P3. Existing major-effect QTL were validated: WM5.4 (34%-phenotypic variation explained) and WM7.4 (20%) in straw tests, and WM2.2 (15%) and WM3.1 (27%) under field conditions. QTL for avoidance traits such as resistance to lodging and late maturity overlapped WM2.2 in P2 and WM1.5, WM3.1, WM5.4, and WM7.7 in P3. WM5.4 (Pv05: 7.0-38.7 Mb) was associated with a large Phaseolus coccineus L. genome introgression in the resistant parent VCP-13. These findings offer narrowed genomic intervals and putative candidate genes for marker-assisted selection targeting white mold resistance improvement in pinto beans.

Antagonistic systemin receptors integrate the activation and attenuation of systemic wound signaling in tomato

Pattern recognition receptor (PRR)-mediated perception of damage-associated molecular patterns (DAMPs) triggers the first line of inducible defenses in both plants and animals. Compared with animals, plants are sessile and regularly encounter physical damage by biotic and abiotic factors. A longstanding problem concerns how plants achieve a balance between wound defense response and normal growth, avoiding overcommitment to catastrophic defense. Here, we report that two antagonistic systemin receptors, SYR1 and SYR2, of the wound peptide hormone systemin in tomato act in a ligand-concentration-dependent manner to regulate immune homeostasis. Whereas SYR1 acts as a high-affinity receptor to initiate systemin signaling, SYR2 functions as a low-affinity receptor to attenuate systemin signaling. The expression of systemin and SYR2, but not SYR1, is upregulated upon SYR1 activation. Our findings provide a mechanistic explanation for how plants appropriately respond to tissue damage based on PRR-mediated perception of DAMP concentrations and have implications for uncoupling defense-growth trade-offs.

Molecular mapping of candidate genes in determining red color of perilla leaf

Perilla frutescens is a traditional medicinal plant and functional food in Asian communities, characterized by distinct red and green leaf types that have significant phenotypic and medicinal implications. However, the genetic mechanisms controlling anthocyanin synthesis in this species remain unclear. Genetic analysis serves as a powerful tool for investigating the pivotal genes and regulatory mechanisms governing anthocyanin accumulation in red and green perilla. In this study, an F2 segregation population was constructed from a hybrid of red and green perilla, and representative samples were subjected to mix-sequencing using BSA-seq and BSR-seq. A 6.0 Mb candidate region on chromosome 8 was identified, pinpointing PfMYB113b, PfC4H1, and PfF3H as key genes involved in anthocyanin biosynthesis. The insertion of a repeat sequence in the promoter of PfMYB113b leads to alterations in gene expression levels. Furthermore, PfMYB113b regulates the transcription of PfC4H1 and PfF3H, thereby influencing anthocyanin synthesis. These findings enhance our understanding of the genetic regulatory mechanisms underlying leaf coloration in perilla.

Identification and validation of genomic regions for pod shatter resistance in Brassica rapa using QTL-seq and traditional QTL mapping

BACKGROUND

Pod shatter resistance is an important trait in Brassica species, significantly impacting the yield and profitability of growers. Identifying genomic regions and understanding genes underlying shatter resistance is a major objective of breeding programs. Brassica rapa, commonly known as rape or field mustard, is an ancestral species of Brassica napus and Brassica juncea - the most widely oilseed crops grown worldwide. In this study, we performed diversity analysis of B. rapa accessions, bulked segregant analysis based quantitative trait locus-sequencing (QTL-seq), and traditional quantitative trait locus (QTL) mapping in an F2 population to identify genomic regions associated with pod shatter resistance in B. rapa.

RESULTS

A considerable genetic variation for pod shatter resistance, measured as rupture energy (RE), varied from 0.63 to 3.49 mJ(½) was revealed among 90 accessions of B. rapa. Cluster analysis based on 10,324 DArTseq markers showed that pod shatter-resistant accessions originated from diverse sources. We further investigated the genetic and anatomical bases of variation in pod shatter resistance from two contrasting parental lines, ATC90153 (maternal parent with high RE) and ATC91215 (paternal parent with low RE). Bulked segregant resequencing analysis of parental lines and two pooled samples, prepared from 10 resistant and 10 sensitive lines to pod shatter, identified three genomic regions for shatter resistance on chromosomes A06 and A09. Traditional QTL analysis validated marker-pod shatter resistance associations on chromosomes A06 and A09 in the same F2 population using a linkage map based on 23,274 DArTseq markers. Physical positions of significantly associated markers and the priori pod dehiscence genes on the B. rapa reference genome sequence suggested BEE1/PEROXIDASE/TCP8 on A06 and ADPG1/SHP1/MYB116 genes on A09 as potential candidates for pod shatter resistance. Sequence comparison of parental lines identified sequence variants (194 SNPs and 74 InDELs on A06, and two SNPs and two InDELs on A09) in the promoter and downstream regions of B. rapa genes within the QTL.

CONCLUSIONS

We identified QTLs and priori candidate genes associated with variation in pod shatter resistance on chromosomes A06 and A09 in B. rapa. This study provides potential gene targets to understand molecular mechanisms and improve pod shatter resistance in Brassica crops.

Genetic Analysis of the Awn Length Gene in the Rice Chromosome Segment Substitution Line CSSL29

Awn length is a significant agronomic trait in rice. To analyze the genetic mechanism of awn length in the chromosome segment substitution line 29 (CSSL29) derived from 9311 (recipient) into Nipponbare (NIP, donor), an F2 segregated population was constructed from 9311 (indica) and CSSL29. The population and candidate genes were analyzed using quantitative trait loci sequencing (QTL-seq), yeast two-hybrid assays, and 3 k and 10 k rice population databases. The results indicated that the awn length in the F2 segregating population followed a normal distribution, and the long-awn phenotype in CSSL29 was controlled by multiple genes. Through BSA sequencing data, a major QTL qAWN4 associated with rice awn length was identified on chromosome 4, containing the cloned gene An-2. Further investigation of the CSSL29 long-awn substitution segment revealed the presence of the awn length gene An-1, with both genes exhibiting an additive effect on the regulation of the long-awn phenotype. Yeast two-hybrid experiments confirmed no interaction between An-2 and An-1, suggesting that additive effect awn length regulation is not mediated through simple protein-to-protein binding. Population genetic analysis indicated that the An-2 allele was artificially selected during domestication but did not significantly differ between indica and japonica subspecies. These findings enhance our understanding of the genetic regulation of rice awn length and the domestication of long-awn rice, laying the groundwork for future research in this area.

The calcium-dependent protein kinase CmaCPK4 regulates sex determination in pumpkin (Cucurbita maxima D.)

Pumpkin (Cucurbita maxima D.) is typically monoecious with individual male and female flowers, and its yield is associated with the degree of femaleness, i.e. the ratio of female to male flowers produced by the plant. Subgynoecy represents a sex form with a high degree of femaleness, but the regulatory mechanisms in pumpkin remain poorly understood. In this study, using the F2 population crossed from the subgynoecious line 2013-12 and the monoecious line 9-6, we initially identified a recessive locus to control the subgynoecious trait and named it sg1. After bulked segregant analysis with whole-genome resequencing and molecular marker linkage analysis, the sg1 locus was mapped to pumpkin Chromosome 2. Genetic sequence analysis found a pumpkin calcium-dependent protein kinase (CPK) gene, CmaCPK4, in the mapping interval as the candidate gene. A retrotransposon insertion identified within the promoter elevated CmaCPK4 expression in 2013-12. Morphological characterization of near-isogenic lines containing the sg1 allele showed increases in the ratio of female flowers and high ethylene contents in terminal buds compared with the receptor parent. Heterologous overexpression of CmaCPK4 significantly increased the ratio of female flowers in cucumber (Cucumis sativus). Furthermore, CmaCPK4 directly interacts with and phosphorylates 1-aminocyclopropane-1-carboxylate synthase 5 (CmaACS5) and 1-aminocyclopropane-1-carboxylate synthase 7 (CmaACS7), resulting in increased ethylene content in 2013-12, which affected pumpkin sex determination. These findings provide insights into the role of the CmaCPK4-CmaACS5/CmaACS7 module in ethylene-induced sex determination in pumpkin.

Exploration of the molecular mechanism behind a novel natural genic male-sterile mutation of 1205A in Brassica napus

The use of a male sterility hybrid seed production system has resulted in a significant increase in rapeseed yields by over 20%. Nevertheless, the mechanisms underlying male sterility remain largely unexamined. This study presents a spontaneous recessive genic male-sterile (RGMS) mutant of 1205A, which was employed to establish two two-line hybrid production systems: 1205AB and NT7G132AB. Cytological investigations reveal that the mutation occurs at the early microspore stage, resulting in premature degradation of pollen. Through inheritance analysis, linkage mapping, and bulked-segregant analysis sequencing (BSA-Seq), a single gene locus, designated Bna1205ams1, was identified within the QTL region on chrC03 (15.36-18.90 Mb). The development of three newly co-segregated kompetitive allele-specific PCR (KASP) markers, in conjunction with two traditional co-segregated markers, allowed for the refinement of the QTL of Bna1205ams1 to a segment of 181.47 kb. This refinement facilitated the identification of a candidate gene, BnaC03g27700D, through functional and expression analyses. Furthermore, the subcellular localization of BnaC03g27700D was examined. Metabolic fluctuations associated with the fertility gene were observed, particularly in processes related to aborted tapetal programmed cell death (PCD), which may contribute to reduced pollen fertility with abnormal pollen exine. A strong correlation was also established between BnaC03g27700D and thirteen metabolites. This study not only offers valuable insights into the research and practical application of plant male sterility but also serves as a case study on the genetic regulatory mechanisms governing male sterility.

Two novel alleles of the MYB transcription factor BjA06.GL1 and BjB02.GL1 control leaf trichomes and enhance resistance to aphids in Brassica juncea

Leaf trichome formation is a very important agronomic trait as it confers resistance to biotic and abiotic stresses, but the causal genes involved in this process in Brassica juncea remain largely unexplored. In this study, we first characterized the haplotypes of BjB02.GL1 among different inbred lines with leaf trichomes or glabrous leaves. A comparative analysis of the number and density of leaf trichomes between the two mustard inbred lines was then performed. BSA analysis of leaves with trichomes and glabrous pools from the F2 segregating population mapped the candidate genes on Chr.A06 and Chr.B02. Two candidate genes, BjA06.GL1 and BjB02.GL1, were subsequently cloned. After sequence alignment of the BjGL1 genes, both single-nucleotide polymorphisms (SNPs) and indel were identified in the BjA06.GL1 and BjB02.GL1 genes. And quantitative real-time polymerase chain reaction (qRT-PCR) analysis further confirmed that both the BjA06.GL1 and BjB02.GL1 genes were more highly expressed in leaves with trichomes than in glabrous leaves. As the leaf size increased, the leaf trichome density decreased. Gene editing of both BjA06.GL1 and BjB02.GL1 changed the leaf trichome to a glabrous leaf phenotype in mustard. In addition, plants with leaf trichomes presented greater resistance to aphids. Taken together, our results revealed that both BjA06.GL1 and BjB02.GL1 positively regulate leaf trichome formation and help increase aphid resistance in mustard. This study provides valuable resources and helps to elucidate the molecular mechanism of leaf trichome formation in B. juncea.

Fine mapping of the Chilli veinal mottle virus resistance 4 (cvr4) gene in pepper (Capsicum annuum L.)

KEY MESSAGE

The single recessive Chilli veinal mottle virus resistance locus, cvr4, was fine-mapped in pepper through bulked segregant RNA sequencing combined with gene silencing analysis. Chilli veinal mottle virus (ChiVMV) is a widespread pathogen affecting the production of peppers (Capsicum annuum L.) in Asia and Africa. Few loci conferring resistance to ChiVMV have been identified, severely limiting the development of resistant cultivars. To identify ChiVMV resistance genes, we constructed an F2:3 segregating population derived from a cross between the ChiVMV-resistant cultivar 'CV9' and the susceptible cultivar 'Jeju'. The inheritance study of F2:3 populations showed a 1:3 ratio of resistant to susceptible individuals, demonstrating the existence of a single recessive ChiVMV resistance gene in CV9; we named this gene cvr4. To map the cvr4 locus, we employed bulked segregant analysis by RNA sequencing (BSR-seq) of pools from resistant and susceptible F2:3 individuals. We mapped cvr4 to the telomeric region of pepper chromosome 11. To narrow down the cvr4 locus, we developed additional molecular markers in the cvr4 target region, leading to a 2-Mb region of chromosome 11 showing complete co-segregation with the ChiVMV resistance phenotype. Using the polymorphisms identified during BSR-seq, we defined a list of 15 candidate genes for cvr4, which we tested through virus-induced gene silencing analysis for ChiVMV resistance. Of these, the silencing of several genes (DEM.v1.00021323, DEM.v1.00021336, and DEM.v1.00021337) restricted virus spread. Although DEM.v1.00021323 transcript levels were similar between the resistant and susceptible bulks, its alternative spliced isoforms differed in abundance, suggesting that the splicing variants of DEM.v1.00021323 might affect viral infection. These findings may facilitate the breeding of ChiVMV-resistant cultivars in pepper.

A favorable natural variation in CCD7 from orchardgrass confers enhanced tiller number

Tiller number is a crucial determinant that significantly influences the productivity and reproductive capacity of forage. The regeneration potential, biomass production, and seed yield of perennial forage species are highly reliant on the development of tillering. Strigolactones (SLs) are recently discovered carotenoid-derived phytohormones that play a crucial role in the regulation of tillering in annual crops. However, the modulation of tiller growth in perennial forage by SLs remains insufficiently investigated. In this study, we identified two alleles of the SLs biosynthesis gene, DgCCD7A and DgCCD7D, which encode CAROTENOID CLEAVAGE DIOXYGENASE 7 (CCD7), from two distinct subspecies of orchardgrass (Dactylis glomerata) exhibiting contrasting tillering phenotype and SLs content. The functionality of the DgCCD7A allele derived from high-tillering phenotypic orchardgrass was found to be diminished compared to that of DgCCD7D from the low-tillering type in rescuing the increased branching phenotype of CCD7-defective mutants in Arabidopsis and rice (Oryza sativa). Notably, the introduction of DgCCD7A in rice resulted in an increase in tiller number without significantly compromising grain yield. Moreover, we demonstrated that the L309P variation in DgCCD7A is a rare natural variant exclusively found in orchardgrass. Our findings revealed that DgCCD7A, a rare favorable natural variation of CCD7 in orchardgrass, holds significant potential for breeding application in improving the plant architecture of perennial forage and crops.

Bcwf regulates the white petal color in pak choi [Brassica campestris (syn. Brassica rapa) ssp. chinensis]

Flower color is important in determining the ornamental value of Brassica species. However, our knowledge about the regulation of flower color in pak choi [Brassica campestris (syn. Brassica rapa) ssp. chinensis] is limited. In this study, we investigated the molecular mechanism underlying white flower traits in pak choi by analyzing a genetic population with white and yellow flowers. Our genetic analysis revealed that the white trait is controlled by a single recessive gene called Bcwf. Through BSA-Seq and fine mapping, we identified a candidate gene, BraC02g039450.1, which is similar to Arabidopsis AtPES2 involved in carotenoid ester synthesis. Sequence analysis showed some mutations in the promoter region of Bcwf in white flowers. Tobacco transient assay confirmed that these mutations reduce the promoter's activity, leading to downregulation of Bcwf expression in white flowers. Furthermore, the silencing of Bcwf in pak choi resulted in lighter petal color and reduced carotenoid content. These findings provide new insights into the molecular regulation of white flower traits in pak choi and highlight the importance of Bcwf in petal coloring and carotenoid accumulation.

Natural variation in MdNAC5 contributes to fruit firmness and ripening divergence in apple

Fruit firmness is an important trait for characterizing the quality and value of apple. It also serves as an indicator of fruit maturity, as it is a complex trait regulated by multiple genes. Resequencing techniques can be employed to elucidate variations in such complex fruit traits. Here, the whole genomes of 294 F 1 hybrids of 'Fuji' and 'Cripp's Pink' were resequenced, and a high-density binmap was constructed using 5014 bin markers with a total map distance of 2213.23 cM and an average map distance of 0.44 cM. Quantitative trait loci (QTLs) of traits related to fruit were mapped, and an A-T allele variant identified in the coding region of MdNAC5 was found to potentially regulate fruit firmness and ripening. The overexpression of MdNAC5 A resulted in higher production of methionine and 1-aminocyclopropanecarboxylic acid compared to MdNAC5 T , leading to reduced fruit firmness and accelerated ripening in apples and tomatoes. Furthermore, the activities of MdNAC5 A and MdNAC5 T were enhanced through their differential binding to the promoter regions of MdACS1 and MdERF3. Spatial variations in MdNAC5 A and MdNAC5 T caused changes in MdACS1 expression following their interaction with MdERF3. Ultimately, utilizing different MdNAC5 alleles offers a strategy to manipulate fruit firmness in apple breeding.

Identification of the Mitf gene mutation causing congenital deafness and pigmentation disorders in porcupines using BSA-Seq

Worldwide, congenital deafness and pigmentation disorders impact millions with their diverse manifestations, and among these genetic conditions, mutations in the Microphthalmia-associated transcription factor (MITF: OMIM#156845) gene are notable for their profound effects on melanocyte development and auditory functions. This study reports a novel porcupine model exhibiting spontaneous deafness and pigmentation abnormalities reminiscent of human Waardenburg Syndrome Type 2 (WS2: OMIM#193510). Through phenotypic characterization, including coat color, skin, eye morphology, and auditory brainstem response (ABR) assessments, we identified hypopigmentation and complete deafness in mutant porcupines. To pinpoint the genetic basis, a breeding program was established, and Bulk Segregant Analysis (BSA) combined with RNA sequencing was conducted. Primers based on the identified candidate genes were designed for PCR amplification, followed by verification through Sanger sequencing. Through BSA analysis, we identified a total of 88 SNP and 336 InDel candidate sites. By annotating the Mitf gene, we obtained four unique transcript sequences. The SNP and InDel sites within the porcupine Mitf gene sequence, identified through BSA screening, were analyzed in conjunction with the gene's annotation results. This analysis revealed a specific mutation site, Mitf c.875_877delGAA p. (Arg217del), which was subsequently verified by Sanger sequencing. This naturally occurring Mitf mutation in porcupines provides a valuable model for studying the mechanisms underlying WS2 and exploring potential therapeutic strategies for deafness and pigmentation-related disorders.

Irregular green netting of eggplant fruit peel: a domestication trait controlled by SmGLK2 with potential for fruit colour diversification

The distribution of chlorophylls in eggplant (Solanum melongena) peel exhibits either a uniform pattern or an irregular green netting pattern. The latter, manifested as a gradient of dark green netting that is intensified in the proximal part of the fruit on a pale green background, is common in wild relatives and some eggplant landraces. Despite the selection of uniform chlorophylls during domestication, the netting pattern contributes to a greater diversity of fruit colours. Here, we used over 2300 individuals from different populations, including a multi-parent advanced generation inter-cross population for candidate genomic region identification, an F2 population for bulked segregant analysis by sequencing, and advanced backcrosses for edges-to-core fine-mapping, to identify SmGLK2 gene as responsible for the irregular netting in eggplant fruits. We also analysed the gene sequence of 178 S. melongena accessions and 22 wild relative species for tracing the evolutionary changes that the gene has undergone during domestication. Three different mutations were identified leading to the absence of netting. The main causative indel induces a premature stop codon disrupting the protein conformation and function, which was confirmed by western blot analysis and confocal microscopy observations. SmGLK2 has a major role in regulating chlorophyll biosynthesis in eggplant fruit peel.

MdNAC5: a key regulator of fructose accumulation in apple fruit

The sweetness of apple fruit is a key factor in the improvement of apple varieties, with fructose being the sweetest of the soluble sugars, playing a crucial role in determining the overall sweetness of the apple. Therefore, uncovering the key genes controlling fructose accumulation and deciphering the regulatory mechanisms of fructose are vitally important for the improvement of apple varieties. In this study, through BSA-seq and transcriptome analysis of the 'Changfu 2' × 'Golden Delicious' F1 hybrid population, MdNAC5 was identified as a key regulatory gene for fructose content. MdNAC5 was shown to significantly influence fructose accumulation in both apples and tomatoes. Furthermore, we conducted a detailed identification of sugar transporters and metabolic enzymes in apples, discovering that MdNAC5 can enhance fructose accumulation in vacuoles and the conversion of sucrose to fructose by binding to and activating the promoters of the vacuolar sugar transporter MdTST2 and the neutral invertase MdNINV6. Additionally, MdNAC5 regulated the MdEIN3.4-MdSWEET15a module, strengthening the unloading of sucrose in the phloem of the fruit. Our results reveal a new mechanism by which MdNAC5 regulates fructose accumulation in apples and provide theoretical foundations for improving apple sweetness through genetic modification.

Natural variations in the Cis-elements of GhRPRS1 contributing to petal colour diversity in cotton

The cotton genus comprises both diploid and allotetraploid species, and the diversity in petal colour within this genus offers valuable targets for studying orthologous gene function differentiation and evolution. However, the genetic basis for this diversity in petal colour remains largely unknown. The red petal colour primarily comes from C, G, K, and D genome species, and it is likely that the common ancestor of cotton had red petals. Here, by employing a clone mapping strategy, we mapped the red petal trait to a specific region on chromosome A07 in upland cotton. Genomic comparisons and phylogenetic analyses revealed that the red petal phenotype introgressed from G. bickii. Transcriptome analysis indicated that GhRPRS1, which encodes a glutathione S-transferase, was the causative gene for the red petal colour. Knocking out GhRPRS1 resulted in white petals and the absence of red spots, while overexpression of both genotypes of GhRPRS1 led to red petals. Further analysis suggested that GhRPRS1 played a role in transporting pelargonidin-3-O-glucoside and cyanidin-3-O-glucoside. Promoter activity analysis indicated that variations in the promoter, but not in the gene body of GhRPRS1, have led to different petal colours within the genus. Our findings provide new insights into orthologous gene evolution as well as new strategies for modifying promoters in cotton breeding.

Polyploid QTL-seq identified QTLs controlling potato flesh color and tuber starch phosphorus content in a plexity-dependent manner

The progenies of polyploid crops inherit multiple sets of homoeologous chromosomes through various combinations, which impedes the identification of the quantitative trait loci (QTL) governing agronomic traits and the implementation of DNA marker-assisted breeding. Previously, we developed a whole-genome sequencing-based polyploid QTL-seq method that utilizes comprehensively extracted simplex polymorphisms for QTL mapping. Here, we verified the detection of duplex QTLs by modifying the analytical settings to explore the QTLs governing tuber flesh color and starch phosphorus content using tetraploid potato (Solanum tuberosum L.). The F1 progenies were obtained from a cross between 'Touya' (TY) and 'Benimaru' (BM). A single TY-derived QTL responsible for yellow flesh color was identified around a β-carotene hydroxylase gene on chromosome 3 using simplex polymorphisms, and a BM-derived QTL associated with decreased starch phosphorus content near a starch synthase II gene on chromosome 2 was detected using duplex polymorphisms. Furthermore, linked DNA markers were developed at the QTL sites. For the latter QTL, plexity-distinguishable markers were developed using quantitative PCR, fragment analysis, and amplicon sequencing. These revealed the allele dosage-dependent effect of the reduced starch phosphorus content. Thus, the polyploid QTL-seq pipeline can explore versatile QTLs beyond simplex, facilitating DNA marker-assisted breeding in various polyploid crops.

A novel Pik allele confers extended resistance to rice blast

In the ongoing arms race between rice and Magnaporthe oryzae, the pathogen employs effectors to evade the immune response, while the host develops resistance genes to recognise these effectors and confer resistance. In this study, we identified a novel Pik allele, Pik-W25, from wild rice WR25 through bulked-segregant analysis, creating the Pik-W25 NIL (Near-isogenic Lines) named G9. Pik-W25 conferred resistance to isolates expressing AvrPik-C/D/E alleles. CRISPR-Cas9 editing was used to generate transgenic lines with a loss of function in Pik-W25-1 and Pik-W25-2, resulting in loss of resistance in G9 to isolates expressing the three alleles, confirming that Pik-W25-induced immunity required both Pik-W25-1 and Pik-W25-2. Yeast two-hybrid (Y2H) and split luciferase complementation assays showed interactions between Pik-W25-1 and the three alleles, while Pik-W25-2 could not interact with AvrPik-C, -D, and -E alleles with Y2H assay, indicating Pik-W25-1 acts as an adaptor and Pik-W25-2 transduces the signal to trigger resistance. The Pik-W25 NIL exhibited enhanced field resistance to leaf and panicle blast without significant changes in morphology or development compared to the parent variety CO39, suggesting its potential for resistance breeding. These findings advance our knowledge of rice blast resistance mechanisms and offer valuable resources for effective and sustainable control strategies.

Exploring the genomic landscape of gummy stem blight resistance in watermelon through QTL-Seq

BACKGROUND

Watermelon is a nutritionally and economically significant crop in the US and globally. Gummy Stem Blight (GSB), caused by three cryptic Stagonosporopsis species, is one of the most devastating diseases affecting watermelon in the US, impacting most of the plant's above-ground parts. This study aimed to identify key Quantitative Trait Variants (QTVs) that include SNPs and In/Dels associated with GSB resistance in selfed derivatives of advanced multicross interspecific derivatives population derived from intercrosses between the most resistant lines of Citrullus amarus and highly susceptible cultivars of Citrullus lanatus.

RESULTS

Resistant and susceptible bulks were created by combining equimolar DNA concentrations from 30 extremely resistant derivatives and 30 extremely susceptible lines. These bulks underwent whole-genome sequencing, generating over 1 billion reads per bulk to achieve comprehensive genome coverage. The mapping percentage of the bulks to the parental genomes ranged from 92 to 99%. More than 6 million SNPs and 1 million indels were identified from the resistant parental genome, compared to fewer than 2 million SNPs and 0.4 million indels from the susceptible parental genome. QTNs associated with GSB resistance were identified using single-nucleotide polymorphism-index and Gprime methods. Statistically significant variants/loci linked to GSB resistance were found on chromosomes 1, 2, 3, 5, 7, 10, and 11. Notably, the genes Lipase class 3 family protein, Ribosome hibernation promotion factor (CaU02G00010), Ubiquitin-like-specific protease 1D (CaU03G04260), and Zinc finger CCCH domain-containing 15 (CaU03G10970) harbored the highest delta SNPs. Several previously published genes, including Avr9/Cf-9 Rapidly Elicited Protein (CaU07G12990) on chromosome 7, were also identified.

CONCLUSIONS

Identifying significant loci associated with GSB resistance has facilitated the development of PACE assays, which will aid in breeding GSB-resistant watermelon cultivars. These findings provide critical insights into the genetic basis of GSB resistance and represent a significant step towards improving the resilience of watermelon crops against this devastating disease.

QTL-seq and QTL mapping identify a new locus for Cercospora leaf spot (Cercospora canescens) resistance in mungbean (Vigna radiata) and a cluster of Receptor-like protein 12 (RLP12) genes as candidate genes for the resistance

KEY MESSAGE

QTL-seq, linkage mapping, and whole-genome resequencing revealed a new locus (qCLS5.1) controlling Cercospora canescens resistance in mungbean and Receptor-like protein 12 (RLP12) genes as candidate genes for the resistance. Cercospora leaf spot (CLS) disease, caused by Cercospora canescens, is a common disease of mungbean (Vigna radiata). In this study, the genetics of CLS resistance was investigated in a new source of resistance (accession V2817) and the resistance was finely mapped to identify candidate genes. F2 and F2:3 populations of the cross V1197 (susceptible) × V2718 and a BC1F1 population of the cross V1197 × (V1197 × V2817) were used in this study. Segregation analysis suggested that the resistance is controlled by a single dominant gene. QTL-seq using F2 individuals revealed that a single QTL (designated qCLS5.1) on chromosome 5 controlled the resistance. The qCLS5.1 was confirmed in the F2:3 and BC1F1 populations by QTL analysis. Fine mapping using 978 F2 individuals localized qCLS5.1 to a 48.94 Kb region containing three tandemly duplicated Receptor-like protein 12 (RLP12) genes. Whole-genome resequencing and alignment of V1197 and V2817 revealed polymorphisms causing amino acid changes and premature stop codons in the three RLP12 genes. Collectively, these results show that qCLS5.1 is a new locus for CLS resistance in mungbean, and a cluster of RLP12 genes are candidate genes for the resistance. The new locus qCLS5.1 will be useful for molecular breeding of durable CLS-resistant mungbean cultivars.