Identification of markers linked to disease-resistance genes by bulked segregant analysis: a rapid method to detect markers in specific genomic regions by using segregating populations

Genomics and transcriptomics identify quantitative trait loci affecting growth-related traits in silver pomfret (Pampus argenteus)

Pampus argenteus, a species distributed throughout the Indo-West Pacific, plays a significant role in the yield of aquaculture species. However, cultured P. argenteus has always been characterised by unbalanced growth synchronisation among individuals, slow growth rate, and lack of excellent germplasm resources. Therefore, we conducted mass selection for fast-growing strain P. argenteus for several consecutive years. Various genetic improvement programs have modified its genome sequence through selective pressure, leaving nucleotide signals that can be detected at the genomic level. In the present study, we combined bulked segregant analysis and transcriptome sequencing to identify candidate single nucleotide polymorphisms (SNPs) and key genes for growth-related traits in P. argenteus. A total of 7,280,936 SNPs and 2,212,379 insertions/deletions were identified in the extreme phenotypes of the fast-growing and slow-growing groups. Based on the examination of SNP frequency differences and sliding-window analysis, 42 SNPs were identified as candidate markers. Moreover, 14 of the 42 SNPs linked to growth-related traits were confirmed to be credible SNPs, and eight growth-related genes were screened, namely myb-binding protein 1 A, insulin A/B chains, α-1B adrenoceptor, engulfment and cell motility protein 3, myosin light chain kinase family member 4, insulin receptor located, unconventional myosin-9b, and matrilin-1. An optimal three-factor model (SNP4&SNP12&SNP14) was constructed using the generalized multifactor dimensionality reduction method, and its accuracy was verified as 67.72 %. These results may benefit genetic studies and accelerate genetic improvement of fast-growing strains of P. argenteus.

Fine mapping and characterization of Fusarium wilt (Fusarium oxysporum f. sp. benincasae) resistance gene Fob1(t) in wax gourd (Benincasa hispida Cogn.)

Fusarium wilt (FW), caused by the plant fungus Fusarium oxysporum, causes severe economic losses in wax gourd (Benincasa hispida Cogn.). Developing disease-resistant varieties is an effective measure. Currently, no reports exist on gene localization or cloning of genes associated with FW resistance in wax gourd. However, our team has identified an FW resistance gene and plans to apply it in resistant varieties. In this study, we used bulked segregant analysis sequencing and quantitative trait locus detection on 1,304 inbred lines derived from resistant GD68 and susceptible HM25 parents to identify FW resistance genes. We successfully identified a resistance locus between the 3M13.385 and 3M16.869 markers on chromosome 3, named Fob1(t). Fine mapping between markers 3M15.904 and 3M16.373 (469 kb apart) identified 22 candidate genes. The transcriptome, sequencing comparison, and qRT-PCR analyses suggested that the endochitinase gene Bch03G006380 is the resistance gene, with significant expression differences between the parents and a one-base mutation in the first exon. The study also revealed the roles of Fob1(t) in plant hormones, transcription factors, phenylpropane metabolism, oxidation-reduction, disease progression, enzymes, and cell wall modification pathways, enhancing FW resistance in GD68. Finally, we identified two closely linked insertion-deletion markers that can assist in the transfer and utilization of the Fob1(t) gene, significantly improving the screening rate of positive individual plants and reducing breeding time.

Identification of candidate genes associated with resistance against Vibrio parahaemolyticus causing acute hepatopancreatic necrosis disease in Litopenaeus vannamei via bulked segregant analysis

Acute hepatopancreatic necrosis disease (AHPND), a predominant bacterial affliction in the cultivation of the Pacific white shrimp Litopenaeus vannamei, is triggered by Vibrio parahaemolyticus (VpAHPND). In the present study, two populations of L. vannamei were used to construct a full-sib family aimed at identifying SNPs and genes associated with resistance to VpAHPND. Through whole-genome resequencing and bulked segregant analysis, 57 scaffolds were monitored via SNP and INDEL data, respectively, including 215 annotated genes enriched in the MAPK and FoXO pathways and other immune-related pathways. Among these, two scaffolds (Scaffold_1573 and 1504) containing 22 annotated genes with nonsynonymous mutations were both screened out on the basis of the SNP and InDel data. A key immune-related gene, LvPLP1, was identified, and its antibacterial function was further confirmed through dsRNA interference under VpAHPND challenge. Consequently, the mortality of shrimp and Vibrio loads in the hepatopancreas significantly increased under dsRNA interference with LvPLP1 during VpAHPND challenge, suggesting that the LvPLP1 gene is a key functional gene involved in anti-VpAHPND infection in L. vannamei. This study identified functional genes underlying resistance to VpAHPND infection and laid a genetic foundation for marker-assisted selection of L. vannamei shrimp resistant to VpAHPND.

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.

Identification of core candidate genes responding to Verticillium wilt (Verticillium dahliae) in cotton via integrated methods

Cotton is a vital natural fiber and oil crop, yet it is severely affected by verticillium wilt (VW), known as the 'cancer' of cotton, hindering the industry's sustainable development. Upland cotton, which is widely cultivated, lacks effective resistance to VW, while most sea island cotton shows strong resistance. In this study, an F2:3 population was constructed by hybridizing the verticillium wilt-resistant island cotton variety 'Hai7124' with the susceptible variety 'Xinhai14'. Using Bulked Segregant Analysis (BSA-seq), we identified 10 genetic intervals significantly associated with resistance. Additionally, two pathogenic strains of Verticillium dahliae, Vd592 (a strong pathogenic type) and VdKT (a weak pathogenic type), were used to infect the 'Hai7124' and 'Xinhai14' for RNA-seq analysis, focusing on differentially expressed genes and signaling pathways in samples treated with different resistant and susceptible materials and infected with different pathogens. By integrating BSA-seq and RNA-seq association analyses, the candidate gene range was further refined. Five genes (GBMYB102, GBWRKY65, GBRDA2, GBSOT16, and GBCWINV1) were validated through virus-induced gene silencing (VIGS). The results revealed that reduced expression of these genes significantly decreases plant disease resistance and leads to a reduction in the activity of defense-related enzymes (such as SOD, CAT or PAL) and secondary metabolites (including lignin or flavonoids). Based on the preliminary functional analysis of these candidate genes, we speculate that redox metabolism and secondary metabolites play crucial roles in the resistance of island cotton to Verticillium wilt, and that the resistance of island cotton to verticillium wilt is the result of multiple genes working together.

Genomics assisted mapping of earliness in pea (Pisum sativum L.)

BACKGROUND

Garden pea (Pisum sativum L.), is a temperate crop belonging to the Leguminosae family. Early maturing pea varieties complete their growth cycle in ∼80-90 days and fits very well within the crop rotation of rice, wheat, and maize, thereby providing an extra source of income to the farmers. Identification of genes associated with the earliness is very important for developing early maturing pea varieties.

METHODS AND RESULTS

In the present study we investigated the genetics of earliness and identified the putative genomic regions associated with the earliness in F2 population derived from a cross between early-maturing (Matar Ageta-10) and late-maturing (Punjab-89) pea varieties using BSA-Seq approach. Genetic analysis revealed that earliness follows a monogenic recessive inheritance pattern. Two extreme phenotypic pools were constructed by identifying ten extreme early and ten extreme late plants from the F2 population, and QTL-seq analysis was performed to obtain major genomic region of 6.5 Mb located at 418.46 Mb to 424.97 Mb on chromosome 7 and has been designated as PsE7. Further, a total of 907 SNPs were identified within this 6.5 Mb genomic region of which seven SNPs were validated through KASP markers. Among these, one marker namely PS423028253 showed association with the earliness trait at distance of 1.7 cM.

CONCLUSION

This novel genomic region along with KASP marker (PS423028253) identified in this study could be used for marker-assisted selection in pea breeding programs and will aid in the identification of the candidate genes in future studies.

Epistasis and cryptic QTL identified using modified bulk segregant analysis of copper resistance in budding yeast

The contributions of genetic interactions to natural trait variation are challenging to estimate experimentally, as current approaches for detecting epistasis are often underpowered. Powerful mapping approaches such as bulk segregant analysis (BSA), wherein individuals with extreme phenotypes are pooled for genotyping, obscure epistasis by averaging over genotype combinations. To accurately characterize and quantify epistasis underlying natural trait variation, we have engineered strains of the budding yeast Saccharomyces cerevisiae to enable crosses where one parent's chromosome is fixed while the rest of the chromosomes segregate. These crosses allow us to use BSA to identify quantitative trait loci (QTL) whose effects depend on alleles on the fixed parental chromosome, indicating a genetic interaction with that chromosome. Our method, which we term epic-QTL (for epistatic-with-chromosome QTL) analysis, can thus identify interaction loci with high statistical power. Here, we perform epic-QTL analysis of copper resistance with chromosome I or VIII fixed in a cross between divergent naturally derived strains. We find 7 loci that interact significantly with chromosome VIII and none that interact with chromosome I, the smallest of the 16 budding yeast chromosomes. Each of the 7 interactions alters the magnitude, rather than the direction, of an additive QTL effect. We also show that fixation of one source of variation-in this case, chromosome VIII, which contains the large-effect QTL mapping to CUP1-increases power to detect the contributions of other loci to trait differences.

The Identification of a Single-Base Mutation in the Maize Dwarf 1 Gene Responsible for Reduced Plant Height in the Mutant 16N125

Maize (Zea mays L.) is a globally vital crop for food, feed, and biofuel production, with plant height (PH) being a key agronomic trait that significantly influences yield, lodging resistance, and stress tolerance. This study identified a single-base mutation in the D1 (Dwarf 1) gene responsible for the dwarf phenotype in the maize mutant 16N125. Through genetic analysis and fine mapping, the candidate region was localized to chromosome 3, narrowing it down to an interval containing three genes. Sequencing revealed a non-synonymous mutation in D1, which encodes a gibberellin 3-beta-dioxygenase, leading to amino acid substitutions at positions 61 and 123. Genetic analysis of F2 populations confirmed that the mutation at position 61 was responsible for the dwarf trait. Furthermore, the mutation was detected in several Chinese inbred lines, indicating its potential role in dwarfing under specific conditions. These findings provide critical insights into the genetic mechanisms regulating maize plant height, offering valuable information for breeding programs focused on improving crop architecture and yield to address the challenges of global food security and climate change.

Genetic and Genomic Analysis Identifies bcltf1 as the Transcription Factor Coding Gene Mutated in Field Isolate Bc116, Deficient in Light Responses, Differentiation and Pathogenicity in Botrytis cinerea

Natural populations provide valuable information and resources for addressing the genetic characterization of biological systems. Botrytis cinerea is a necrotrophic fungus that exhibits complex responses to light. Physiological analysis of B. cinerea populations from vineyards in Castilla y León (Spain) allowed for the identification of isolate Bc116. This field isolate shows a reduced pathogenicity that is conditioned by the light regime. Light also delays germination and accentuates the negative effect it exerts on the vegetative growth of B. cinerea. Bc116 also displays a marked hyperconidiation phenotype and a characteristic sclerotia production pattern. Genetic analysis demonstrates that the alternative phenotypes regarding pathogenicity, conidiation, and pattern of sclerotia production co-segregate in the progeny of crosses between isolate Bc116 and a wild-type field isolate, Bc448, showing that they are under the control of a single locus. By applying a strategy based on bulked segregant analysis, the mutation in Bc116 has been mapped to a 200 kb region on Chr14 and the analysis of this region identifies a 2 kb deletion affecting the bcltf1 gene, encoding the B. cinerea Light Responsive Transcription Factor 1 described in the reference isolate B05.10. Transformation of Bc116 with the B05.10 bcltf1 allele restored the wild-type phenotypes, providing functional evidence that the natural mutant Bc116 is altered in gene bcltf1. This study offers additional information, derived from the analysis of the genetic background of a natural mutant, on the physiological processes regulated by BcLTF1 and supports the key role of this TF in the pathogenicity and photobiology of B. cinerea.

Twisted Sister1: an agravitropic mutant of bread wheat (Triticum aestivum) with altered root and shoot architectures

We identified a mutant of hexaploid wheat (Triticum aestivum) with impaired responses to gravity. The mutant, named Twisted Sister1 (TS1), had agravitropic roots that were often twisted along with altered shoot phenotypes. Roots of TS1 were insensitive to externally applied auxin, with the genetics and physiology suggestive of a mutated AUX/IAA transcription factor gene. Hexaploid wheat possesses over 80 AUX/IAA genes, and sequence information did not identify an obvious candidate. Bulked segregant analysis of an F2 population mapped the mutation to chromosome 5A, and subsequent mapping located the mutation to a 41 Mbp region. RNA-seq identified the TraesCS5A03G0149800 gene encoding a TaAUX/IAA protein to be mutated in the highly conserved domain II motif. We confirmed TraesCS5A03G0149800 as underlying the mutant phenotype by generating transgenic Arabidopsis thaliana. Analysis of RNA-seq data suggested broad similarities between Arabidopsis and wheat for the role of AUX/IAA genes in gravity responses, although there were marked differences. Here we show that the sequenced wheat genome, along with previous knowledge of the physiology of gravity responses from other plant species, gene mapping, RNA-seq, and expression in Arabidopsis have enabled the cloning of a key wheat gene that defines plant architecture.

Differential methylation of a retrotransposon upstream of a MYB gene causes variegation of lettuce leaves, which is abolished by the presence of an (AT)5 repeat in the promoter

Variegation, a common phenomenon in plants, can be the result of several genetic, developmental, and physiological factors. Leaves of some lettuce cultivars exhibit dramatic red variegation; however, the genetic mechanisms underlying this variegation remain unknown. In this study, we cloned the causal gene for variegation on lettuce leaves and elucidated the underlying molecular mechanisms. Genetic analysis revealed that the polymorphism of variegated versus uniformly red leaves is caused by an "AT" repeat in the promoter of the RLL2A gene encoding a MYB transcription factor. Complementation tests demonstrated that the RLL2A allele (RLL2AV) with (AT)n repeat numbers other than five led to variegated leaves. RLL2AV was expressed in the red spots but not in neighboring green regions. This expression pattern was in concert with a relatively low level of methylation in a retrotransposon inserted in -761 bp of the gene in the red spots compared to high methylation of the retrotransposon in the green region. The presence of (AT)5 in the promoter region, however, stabilized the expression of RLL2A, resulting in uniformly red leaves. In summary, we identified a novel promoter mechanism controlling variegation through inconsistent levels of methylation and showed that the presence of a simple sequence repeat of specific size could stabilize gene expression.

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.

The CsTM alters multicellular trichome morphology and enhances resistance against aphid by interacting with CsTIP1;1 in cucumber

The multicellular trichomes of cucumber (Cucumis sativus L.) serve as the primary defense barrier against external factors, whose impact extends beyond plant growth and development to include commercial characteristics of fruits. The aphid (Aphis gossypii Glover) is one of prominent pests in cucumber cultivation. However, the relationship between physical properties of trichomes and the aphid resistance at molecular level remains largely unexplored. Here, a spontaneous mutant trichome morphology (tm) was characterized by increased susceptibility towards aphid. Further observations showed the tm exhibited a higher and narrower trichome base, which was significantly distinguishable from that in wild-type (WT). We conducted map-based cloning and identified the candidate, CsTM, encoding a C-lectin receptor-like kinase. The knockout mutant demonstrated the role of CsTM in trichome morphogenesis. The presence of SNP does not regulate the relative expression of CsTM, but diminishes the CsTM abundance of membrane proteins in tm. Interestingly, CsTM was found to interact with CsTIP1;1, which encodes an aquaporin with extensive reports in plant resistance and growth development. The subsequent aphid resistance experiments revealed that both CsTM and CsTIP1;1 regulated the development of trichomes and conferred resistance against aphid by affecting cytoplasmic H2O2 contents. Transcriptome analysis revealed a significant enrichment of genes associated with pathogenesis, calcium binding and cellulose synthase. Overall, our study elucidates an unidentified mechanism that CsTM-CsTIP1;1 alters multicellular trichome morphology and enhances resistance against aphid, thus providing a wholly new perspective for trichome morphogenesis in cucumber.

Positional mapping - constitutive purple trait locus (pl) in guava (Psidium guajava L.) in F2 and BC1F1 populations of Purple Local × Allahabad Safeda

Guava (Psidium guajava L.) is a popular fruit crop in Southeast Asia. Landrace Purple Local (PL), also known as Malaysian guava, is rich in anthocyanin content in all the plant parts but has poor yield. Genetic inheritance and physical location of the allele(s) controlling the constitutive purple trait in guava have not yet been reported. The F1 plants generated by cross hybridizing PL × green cv. Allahabad Safeda (AS) and analysed for 7 years did not exhibit purple trait. Evaluation of F1, F2 and BC1F1 populations derived from the cross between PL x AS revealed that purple color is a recessive trait in guava. Equally spaced (~10 Mb) co-dominant polymorphic markers developed by in silico analysis of AS and PL genome (mapped to AS genome assembly) into PCR-based assay mapped the purple color locus (pl) on pseudochromosome 11 (PC 11) of guava genome. Further mapping on PC 11 identified 2 InDel markers at 2.49 Mb (Pg11_INDL_2.49 M) and 4.99 Mb (Pg11_INDL_4.99 M) closely associated with pl. Also, QTLseqr for purple and non-purple bulks in F2 provided two co-localized significant peak ΔSNP-indices at positions 2489072 and 4978573 on PC 11. The mapped genomic interval harbours 85 coding genes, including the potential candidates MYB-like ETC1, anthocyanidin reductase, MYB41-like transcription factors and F-box protein SKIP27-like. Markers flanking pl would potentiate the marker-assisted introgression of anthocyanin trait in popular cultivars of guava.

Breeding Potential for Increasing Carbon Sequestration via Rhizomatous Grain Sorghum

Rhizomes, key carbon sequestration sinks in perennial crops, are hypothesized to exhibit a trade-off with grain yield. This study evaluated rhizomatous grain sorghum populations for increasing carbon sequestration potential. Twelve F3:4 heterogeneous inbred families (HIFs) from a Sorghum bicolor (L.) Moench × Sorghum propinquum (Kunth) Hitchc cross were tested in a greenhouse, and two F4:5 HIF progenies were field tested. Traits measured included rhizome biomass, root biomass, total belowground biomass, and grain yield. Rhizome biomass showed high heritability (0.723) and correlated strongly with belowground biomass (r1 = 0.95; r2 = 0.97) in both F4:5 HIFs, suggesting the potential of rhizomes to sequester carbon. Contrary to the hypothesized trade-off, a positive relationship between rhizome biomass and grain yield was observed, potentially via rhizome-derived shoots, and individual plants pyramiding high rhizome biomass, biomass yield, and grain yield were also identified. Using bulked segregant analysis (BSA), twenty simple sequence repeat (SSR) markers linked to eight genomic regions associated with rhizome presence were identified, with five regions potentially being novel. This study suggests that breeding rhizomatous grain sorghum with high rhizome biomass could enhance carbon sequestration while preserving agronomic yields, offering new insights for future breeding and mapping initiatives.

Exploring resistance mechanisms and identifying QTLs for brown planthopper in tropical and subtropical rice (Oryza sativa L.) germplasm

KEY MESSAGE

A total of 4006 tropical and subtropical rice germplasms were screened for brown planthopper resistance, and the resistance mechanisms of 63 highly resistant accessions were characterized. This led to the designation of three novel resistance QTLs: Bph47, Bph48, and Bph49. The brown planthopper (BPH) is a significant piercing-sucking pest of rice plants that causes widespread destruction globally. Discovering new germplasms and genes for BPH resistance is essential for enhancing genetic diversity in rice breeding. In this study, 4006 rice accessions from tropical and subtropical regions were screened for BPH resistance at the seedling stage, and 63 accessions with high-resistant were identified. Of these, 59 accessions exhibited high resistance to BPH at the adult stage. The 63 accessions displayed widespread variation in key agronomic traits, though most were generally unsatisfactory. Assessments of antixenosis, antibiosis, and tolerance indicated diverse resistance mechanisms in the 63 accessions, with the majority (39/63) demonstrating both antixenosis and antibiosis. Microscopic observations and physiological assessments revealed significant differences in vascular bundle structure, fiber content, and activity of defense-related enzymes between the 63 high-resistance and 27 susceptible ones. Furthermore, correlation analysis highlighted a substantial positive relationship between BPH resistance and parameters such as rice trypsin inhibitor (RTI) levels and width of the sclerenchyma layer (WSL). Genetic analysis of F2:3 segregating populations from four resistant accessions crossed with the susceptible rice variety 9311 identified three novel major-effect quantitative-trait loci (QTLs) located on chromosome 1L (690 kb and 1.84 Mb) and 5S (295 kb). This study significantly enriched the BPH-resistant germplasm sources and genes, highlighting the varied resistance mechanisms of rice against BPH.

Development and Application of Pik Locus-Specific Molecular Markers for Blast Resistance Genes in Yunnan Japonica Rice Cultivars

Rice blast, caused by the fungal pathogen Magnaporthe oryzae, is one of the most devastating diseases affecting rice production worldwide, resulting in significant yield losses and threatening global food security. The severity of rice blast, particularly in susceptible regions, underscores the urgent need for available effective resistance strategies. In this study, six sets of gene-specific molecular markers for the Pik locus associated with rice blast resistance were developed based on publicly available gene sequences. Experimental validation confirmed their high accuracy. During the marker development process, a novel haplotype of the Pik locus was identified. This haplotype is characterized by 14 bp mutations and a 9 bp insertion within the coding sequence region when compared to the Pikh allele. Subsequently, a molecular marker specific to this haplotype was developed and validated. The application of these seven sets of markers to analyze 163 japonica rice cultivars bred in Yunnan Province between 1980 and 2020 revealed that 38.65% of the cultivars carry the Piks allele, indicating a low resistance frequency against the rice blast fungus under field conditions. In contrast, only a small proportion of cultivars possess other Pik locus alleles, which exhibit higher resistance frequencies. These findings highlight the limited utilization of Pik locus genes in japonica rice breeding in Yunnan. Furthermore, 21.47% of the cultivars lack any of the aforementioned Pik locus alleles, indicating the genetic diversity and complexity of the rice genetic resources of Yunnan Province.

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.

Identification of Quantitative Trait Loci and Development of Intermediate Breeding Parent for Rice Sheath Blight Resistance

The soil-borne pathogen Rhizoctonia solani is one of the most devastating necrotrophic pathogens worldwide, responsible for causing rice sheath blight (RSB). This pathogen has a broad host range, affecting economically important monocots and dicots such as rice, wheat, potato, soybean, sugar beet, and cucumber. Despite extensive screening of rice germplasm, genes that confer full resistance to RSB have rarely been identified, leading to slow progress in breeding resistant varieties. To identify RSB-resistant rice cultivars in Korea, variations in quantitatively inherited resistance have been observed. We conducted a study to visually assess the RSB resistance phenotypes of 250 cultivated varieties under natural disease conditions in the field over several years. Notable candidates included P1401, which showed resistance, while Junam was susceptible. To identify the quantitative trait loci (QTLs) associated with resistance, we developed an F2 mapping population by crossing P1401 and Junam, followed by bulked segregant analysis. These QTLs were mapped to specific locations on seven of the 12 rice chromosomes. This mapping population and the resulting datasets provide valuable resources for advancing genomic research in rice, particularly for marker-assisted breeding strategies for enhancing resistance to R. solani and other important agronomic traits.

Unraveling the role of OsSCL26 in transcriptional regulation in rice: Insights into grain shape, heading date, and carbohydrates

Grain shape, heading date, and amylase content are pivotal traits influencing rice yield, quality, distribution, and regional adaptability. Through our investigation, we identified a mutant, characterized by slender grains, elevated amylose content, and early heading date. Histocytologic scrutiny unveiled heightened cell proliferation in the spikelet hull contributing to the slender grain morphology. The OsSCL26 gene, governing these significant traits, was meticulously cloned via fine-mapping. Phenotypic scrutiny of OsSCL26 knockout and overexpression lines validated its pivotal role in trait regulation. Further analysis disclosed a substitution in the OsSCL26 promoter region, creating a novel binding site for the transcript factor OsbZIP47, thereby modulating its expression in the osscl26 mutant. Functionally, OsSCL26, acting as a serine/arginine-rich SC35-like protein, interacted with U1-70K in vivo and in vitro. OsSCL26 exhibited direct binding to genes implicated in grain shape and carbohydrates, thereby regulating their splicing. Moreover, OsSCL26 showed direct and indirect associations with target RNAs involved in circadian rhythm. Overall, our findings elucidate the mechanism of OsSCL26, an RNA binding protein interacting with splicing factor, as a crucial member of the spliceosome, thereby impacting post-transcriptional splicing and regulating grain shape, heading date, and carbohydrates in rice.

The MAP kinase scaffold MORG1 shapes cell death in unresolved ER stress in Arabidopsis

Governed by the unfolded protein response (UPR), the ability to counteract endoplasmic reticulum (ER) stress is critical for maintaining cellular homeostasis under adverse conditions. Unresolved ER stress leads to cell death through mechanisms that are yet not completely known. To identify key UPR effectors involved in unresolved ER stress, we performed an ethyl methanesulfonate (EMS) suppressor screen on the Arabidopsis bzip28/60 mutant, which is impaired in activating cytoprotective UPR pathways. This screen identified MAP kinase organizer 1 (MORG1), a conserved MAP kinase scaffold protein, as a previously uncharacterized regulator of ER stress tolerance. The coffin1 mutant, which carries a mutation in MORG1 , exhibited enhanced resilience to ER stress by partially restoring UPR gene expression and promoting growth under stress conditions. Mechanistically, we found that MORG1 modulates MPK6-dependent phosphorylation of the stress-responsive transcription factor WRKY8. Loss of WRKY8 phenocopied the coffin1 mutant, highlighting WRKY8's role as a key repressor in the UPR. Together, these findings reveal a MORG1-MPK6-WRKY8 signaling axis that fine-tunes UPR gene expression, providing new insights into ER stress regulation and strategies for improving stress tolerance in multicellular eukaryotes.

Resistance genes against yellow rust pathogen in Triticum spelta: a possible new Yr resistance gene in accession IARI276 and Yr5 presence confirmed in PI348764

Yellow/stripe rust caused by Puccinia striiformis f. sp. tritici is a major biotic stress in global wheat production. Introgression lines derived from the Triticum spelta accessions PI348764 and IARI276 showed high levels of yellow rust resistance at seedling and adult plant stage. The Yr5 gene located on 2B chromosome was previously the only stripe rust resistance gene described in T. spelta gene pool. By genotyping parental and introgressed material with markers linked with the Yr5 gene, we demonstrate that PI348764 likely carries Yr5, and that it appears to be absent from IARI276. By employing a combination of methods, including screening for adult plant resistance and seedling resistance at multiple field trials, bulked segregant analysis (BSA) on F5 families, and genotyping using wheat Breeders' 35K array, we show that YrIARI276 is a novel stripe rust resistance gene with putative chromosomal locations on 1BL, 1DL, 5AL, or 7BL. Furthermore, genetic analysis revealed that YrIARI276 showed a goodness of fit to Mendelian ratios for a single dominant gene. As the gene is distinct from Yr5 and the chromosomal location is unique from earlier reported Yr genes, it will be useful in improving diversity of Yr gene repertoire in disease resistance breeding programmes.

Genetic diversity and population structure analysis of short-day onions using molecular markers in association with resistance to Fusarium basal rot

In this research, we analyzed Random Amplified Polymorphic DNA (RAPD), Inter Simple Sequence Repeats (ISSR) and Sequence-related amplified polymorphism (SRAP) markers to evaluate the genetic diversity of eighteen different onion genotypes with various resistant levels to FOC. The results showed that the polymorphism means between RAPD primers was 61.11 to 81.81%; ISSR primers, 62.50 to 81.81%; and SRAP primers, 56.25 to 76.25%. Overall, by assessing MI, PIC, I and H indices, indicating the best thrive in evaluating the genetic diversity of the related onion populations. There is a significant correlation between the generated dendrograms based on similarity matrices. The classification pattern in dendrograms shows a corresponding correlation with the FOC disease severity bunches. So in all three markers studied, 'Saba' and 'Saba - HS', the most resistant ones to FOC disease, were grouped in a branch, and the 'Sahar - HS' and 'Golden Eye', the most susceptible ones were also grouped in another branch separately. This finding indicates that predominant primers act as markers linked to resistance gene(s) against FOC, which can be used to select onions resistant to FOC disease in any breeding scheme.

Genetic Analysis of Stripe Rust Resistance in the Chinese Wheat Cultivar Luomai 163

Stripe or yellow rust (YR) caused by Puccinia striiformis tritici (Pst) is an important foliar disease affecting wheat production globally. Resistant varieties are the most economically and environmentally effective way to manage this disease. The common winter wheat (Triticum aestivum L.) cultivar Luomai 163 exhibited resistance to the Pst races CYR32 and CYR33 at the seedling stage and showed a high level of adult plant resistance in the field. To understand the genetic basis of YR resistance in this cultivar, 142 F5 recombinant inbred lines (RILs) derived from cross Apav#1 × LM163 and both parents were genotyped with the 16K SNP array and bulked segregant analysis sequencing. The analysis detected a major gene, YrLM163, at the seedling stage associated with the 1BL.1RS translocation. Additionally, three genes for resistance at the adult plant stage were detected on chromosome arms 1BL (Lr46/Yr29/Pm39/Sr58), 6BS, and 6BL in Luomai 163, whereas Apav#1 contributed resistance at a quantitative trait locus (QTL) on 2BL. These QTL explained YR disease severity variations ranging from 6.9 to 54.8%. The kompetitive allele-specific PCR (KASP) markers KASP-2BL, KASP-6BS, and KASP-6BL for the three novel loci QYr.hzau-2BL, QYr.hzau-6BS, and QYr.hzau-6BL were developed and validated. QYr.hzau-1BL, QYr.hzau-2BL, and QYr.hzau-6BS showed varying degrees of resistance to YR when present individually or in combination based on genotype and phenotype analysis of a panel of 570 wheat accessions. Six RILs combining resistance alleles of all QTL, showing higher resistance to YR in the field than Luomai 163 with disease severities of 10.7 to 16.0%, are important germplasm resources for breeding programs to develop YR-resistant wheat varieties with good agronomic traits.

A single-base deletion in exon 2 of Hd1 delineates monogenic recessive photoperiod insensitivity in aromatic Joha rice: a novel allele for seasonal adaptability

BACKGROUND

Assam's aromatic Joha rice is a unique rice class famous for its aroma, taste, and nutritional benefits, which fetch high market prices in domestic and international markets. Joha landraces are inherently poor yielders due to their strong aroma and predominantly photoperiod sensitivity. Hybridization involving non-aromatic HYVs improves yield with concomitant loss of quality. In this context, mutation breeding, a sustainable approach where genetic mutations are induced to create desirable traits, often provides useful allelic variation in specific morpho-agronomic traits. The present study delves into the genetic characterization of a photoperiod-insensitive mutant. As part of our mutation breeding programme, this mutant was isolated from a gamma ray-induced M2 population of a Joha rice landrace, Kon Joha.

RESULTS

The mutant was unique, and a single recessive gene conditions the induced photoperiod insensitivity. Mutant gene tagging involved 402 SSR and InDel markers, and later polymorphic markers were used for bulk segregant analysis (BSA) in the F2 population of 'mutant × Kalijeera (distant parent)'. BSA revealed an association between the SSR marker RM527 and this mutant trait. This marker is present on chromosome 6 of the rice genome. Using chromosome 6-specific SSR markers in polymorphic screening and BSA revealed another associated marker, RM19725, for the mutant trait. The genomic interval between RM527 and RM19725 harbors a photoperiod-insensitive gene, Hd1, on chromosome 6. Cloning and sequencing of Hd1 genomic fragments from the parents and mutants revealed a single-base deletion in exon 2, leading to a frameshift mutation in the Hd1 protein. This mutation in exon 2 leads to severe structural abnormalities in the CCT domain of the Hd1 protein that is critical for the interaction of the repressing complex with conserved response elements in the florigen gene under long-day conditions, thereby causing photoperiod insensitivity.

CONCLUSIONS

The mutant's pleasant aroma and other quality characteristics, comparable to those of the parent cultivar, hold significant promise. They expand its potential use in a structured breeding programme aimed at developing high-value aromatic Joha rice. This rice, resilient to winter- and summer-growing environments and with broad seasonal adaptability, could revolutionize the rice market. The practical value of our research is underscored by this exciting possibility.

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.

Genotype-by-environment interactions shape ubiquitin-proteasome system activity

In genotype-by-environment interactions (GxE), the effect of a genetic variant on a trait depends on the environment. GxE influences numerous organismal traits across eukaryotic life. However, we have a limited understanding of how GxE shapes the molecular processes that give rise to organismal traits. Here, we characterized how GxE shapes protein degradation, an essential molecular process that influences numerous aspects of cellular and organismal physiology. Using the yeast Saccharomyces cerevisiae, we characterized GxE in the activity of the ubiquitin-proteasome system (UPS), the primary protein degradation system in eukaryotes. By mapping genetic influences on the degradation of six substrates that engage multiple distinct UPS pathways across eight diverse environments, we discovered extensive GxE in the genetics of UPS activity. Hundreds of locus effects on UPS activity varied depending on the substrate, the environment, or both. Most of these cases corresponded to loci that were present in one environment but not another ("presence / absence" GxE), while a smaller number of loci had opposing effects in different environments ("sign change" GxE). The number of loci exhibiting GxE, their genomic location, and the type of GxE (presence / absence or sign change) varied across UPS substrates. Loci exhibiting GxE were clustered at genomic regions that contain core UPS genes and especially at regions containing variation that affects the expression of thousands of genes, suggesting indirect contributions to UPS activity. Our results reveal highly complex interactions at the level of substrates and environments in the genetics of protein degradation.

Tailoring Genomic Selection for Bos taurus indicus: A Comprehensive Review of SNP Arrays and Reference Genomes

BACKGROUND

Advances in SNP arrays and reference genome assemblies have significantly transformed cattle genomics, particularly for Bos taurus indicus (Zebu cattle). Many commercial SNP arrays were originally designed for Bos taurus taurus, leading to ascertainment bias and the exclusion of crucial SNPs specific to Zebu populations. This review assesses progress in SNP array and reference genome development, with a focus on efforts tailored to Zebu populations and their impact on genomic selection and breeding efficiency.

METHODS

We reviewed the relevant literature on the development of SNP arrays, reference genome assemblies, and SNP genotyping techniques used for Zebu cattle. Emphasis was placed on SNP arrays specifically designed for Zebu breeds, evaluating their contributions to genomic evaluations and identifying limitations in prediction accuracy.

RESULTS

Recent advancements, such as GeneSeek's low- and high-density SNP panels, have aimed to reduce ascertainment bias and include key SNPs for Zebu populations by providing breed-specific panels. These panels have been instrumental in identifying genomic regions associated with economically important traits in Nellore cattle. Studies show that tailored SNP arrays and breed-specific reference genomes can enhance genetic diversity assessment and improve genomic predictions, supporting more effective breeding programs for Zebu cattle.

CONCLUSIONS

Improved SNP arrays and breed-specific reference genomes are crucial for accurate genomic selection in Zebu cattle. Future efforts should prioritize expanding de novo genome assemblies, reducing ascertainment bias, and developing cost-effective genotyping solutions tailored to Zebu populations. Targeted genomic tools will ultimately enable more efficient breeding practices and enhance genomic selection for economically important traits in B. t. indicus cattle.

Epistasis and cryptic QTL identified using modified bulk segregant analysis of copper resistance in budding yeast

The contributions of genetic interactions to natural trait variation are challenging to estimate experimentally, as current approaches for detecting epistasis are often underpowered. Powerful mapping approaches such as bulk segregant analysis, wherein individuals with extreme phenotypes are pooled for genotyping, obscure epistasis by averaging over genotype combinations. To accurately characterize and quantify epistasis underlying natural trait variation, we have engineered strains of the budding yeast Saccharomyces cerevisiae to enable crosses where one parent's chromosome is fixed while the rest of the chromosomes segregate. These crosses allow us to use bulk segregant analysis to identify quantitative trait loci (QTL) whose effects depend on alleles on the fixed parental chromosome, indicating a genetic interaction with that chromosome. Our method, which we term epic-QTL (for epistatic-with-chromosome QTL) analysis, can thus identify interaction loci with high statistical power. Here we perform epic-QTL analysis of copper resistance with chromosome I or VIII fixed in a cross between divergent naturally derived strains. We find seven loci that interact significantly with chromosome VIII and none that interact with chromosome I, the smallest of the 16 budding yeast chromosomes. Each of the seven interactions alters the magnitude, rather than the direction, of an additive QTL effect. We also show that fixation of one source of variation - in this case chromosome VIII, which contains the large-effect QTL mapping to CUP1 - increases power to detect the contributions of other loci to trait differences.

Novel Allelic Mutations in Dw3 Gene That Affect the Height of Sorghum Plants

Breeding for dwarfing traits in sorghum is crucial. However, only three genes (Dw1-Dw3) that control plant height have been mapped. In this study, 634 sorghum cultivars were collected to investigate plant height and genotypes. Four were genotyped Dw1DW2Dw3 (wild type) but with different plant heights, and they were selected to construct two populations and map new dwarf genes. Bulked segregant analysis with whole-genome resequencing of the two populations identified the candidate gene in one same genomic region-on chromosome 7. Then, it was narrowed down to a region containing nine genes. Amino acid and DNA sequence analysis of the parent and offspring plants revealed that two novel allelic mutations in the Dw3 gene play a role in reducing the plant height-8R262 or 8R417, including 1 bp substitution and 2 bp deletions. Furthermore, we sequenced 19 cultivars that primarily exhibited a "one-dwarf" hybrid or wild-type and presumed another allelic mutation via the amino acid alignment of 8R019, 8R100, and 8R402, which was another one-base substitution. These results indicate that multiple types of allelic mutations in the Dw3 gene should be considered when identified or applied.

Introgression and Mapping of a Novel Bacterial Blight Resistance Gene xa49(t) from Oryza rufipogon acc. CR100098A into O. sativa

Bacterial blight (BB) caused by Xanthomonas oryzae pv. oryzae is one of the epidemic diseases in rice. Rapid changes in the pathogenicity of the X. oryzae pv. oryzae pathogen demand the identification and characterization of novel BB resistance genes. Here, we report the transfer and mapping of a new BB resistance gene from Oryza rufipogon acc. CR100098A. Inheritance studies on the BC2F2 population, BC2F3 progenies, and backcross-derived recombinant inbred lines derived from a cross between Pusa44/O. rufipogon acc. CR100098A//2*PR114 showed that a single recessive gene confers resistance in O. rufipogon acc. CR100098A. Bulked segregant analysis using 203 simple sequence repeat (SSR) markers localized the BB resistance gene on chromosome 11 bracketed between two SSR markers, RM27235 and RM2136. Using PR114 and O. rufipogon acc. CR100098A genotyping by sequencing data, 86 KASP markers within the bracketed region were designed and tested for bulked segregant analysis. Only five KASP markers showed polymorphism between parents, and three were associated with the target gene. Seventy-seven new SSR markers were designed from the same interval. A total of 33 polymorphic markers were analyzed on the whole population and mapped the BB gene in an interval of 2.8 cM flanked by SSR markers PAU11_65 and PAU11_44 within a physical distance of 376.3 kb. The BB resistance gene mapped in this study is putatively new and designated as xa49(t). Fourteen putative candidate genes were identified within the xa49(t) region having a role in biotic stress resistance. The linked markers to the xa49(t) gene were validated in other rice cultivars for its successful deployment in BB resistance breeding.

Identification of candidate genes associating with soybean cyst nematode in soybean (Glycine max L.) using BSA-seq

Soybean cyst nematode disease represents the major soil-borne disease of soybean. Identifying disease-resistant genes in soybean has a substantial impact on breeding of disease-resistant crops and genetic improvement. The present work created the F2 population with the disease-resistant line H-10 and disease-susceptible line Chidou4. 30 respective F2disease-resistant and disease-susceptible individuals for forming two DNA pools for whole-genome re-sequencing were selected. As a result, a total of 11,522,230 single nucleotide polymorphism (SNPs) markers from these two parental lines and two mixed pools were obtained. Accordng to SNP-index based association analysis, there were altogether 741 genes out of 99% confidence interval, which were mainly enriched into regions of 38,524,128∼39,849,988 bp with a total length of 1.33 Mb contain 111 genes on chromosome 2, 27,821,012∼29,612,574 bp with a total length of 1.79 Mb contain 92 genes on chromosome 3, 308∼348,214 bp with a total of length 0.35 Mb contain 34 genes on chromosome 10, and 53,867,581∼58,017, 852 bp with a total length of 4.15 Mb contain 504 genes on chromosome 18. Bulk segregant analysis in F2 generations (BSA-seq) was correlated with a disease resistance interval containing 15 genes. Then, using bioinformatics analysis and differential expression analysis, five candidate genes were identified: Glyma.02G211400, Glyma.18G252800, Glyma.18G285800, Glyma.18G287400 and Glyma.18G298200. Our results provides a key basis for analyzing the soybean resistance mechanism against soybean cyst nematode and cloning soybean resistance genes.

The auxin efflux carrier PIN1a regulates vascular patterning in cereal roots

Barley (Hordeum vulgare) is an important global cereal crop and a model in genetic studies. Despite advances in characterising barley genomic resources, few mutant studies have identified genes controlling root architecture and anatomy, which plays a critical role in capturing soil resources. Our phenotypic screening of a TILLING mutant collection identified line TM5992 exhibiting a short-root phenotype compared with wild-type (WT) Morex background. Outcrossing TM5992 with barley variety Proctor and subsequent SNP array-based bulk segregant analysis, fine mapped the mutation to a cM scale. Exome sequencing pinpointed a mutation in the candidate gene HvPIN1a, further confirming this by analysing independent mutant alleles. Detailed analysis of root growth and anatomy in Hvpin1a mutant alleles exhibited a slower growth rate, shorter apical meristem and striking vascular patterning defects compared to WT. Expression and mutant analyses of PIN1 members in the closely related cereal brachypodium (Brachypodium distachyon) revealed that BdPIN1a and BdPIN1b were redundantly expressed in root vascular tissues but only Bdpin1a mutant allele displayed root vascular defects similar to Hvpin1a. We conclude that barley PIN1 genes have sub-functionalised in cereals, compared to Arabidopsis (Arabidopsis thaliana), where PIN1a sequences control root vascular patterning.

Fine mapping and identification of candidate genes associated with powdery mildew resistance in melon (Cucumis melo L.)

Powdery mildew (PM), a common disease of many major crop species, including melon (Cucumis melo L.), affects plant growth and fruit quality and seriously reduces production. Using a combined morphological and molecular approach, we attribute the PM pathogen that naturally occurs in melon to Podosphaera xanthii, and specifically to physiological race 1. An investigation into the genetic basis of PM resistance in melon using the resistant accession 'PI 164637' and susceptible counterpart 'HDZ' reveals dominant inheritance of PM resistance at the seedling stage, supported by F2 and backcross population segregation ratios. Adult plant assessments indicate a major gene with an additive effect for PM resistance. Bulk segregant analysis coupled with high-throughput sequencing identified a significant quantitative trait locus on chromosome 6 that is associated with PM resistance. Genetic mapping narrowed down the candidate region to 63.5 kb using InDel molecular markers, harboring 12 candidate genes. The marker chr06_indel_5 047 127 demonstrated high accuracy in screening PM resistance in an F2 segregating population and 30 inbred lines as natural populations. Functional annotation and expression analysis of candidate genes revealed that MYB transcription factor MELO3C006700, GATA transcription factor MELO3C028829 and heparanase-like protein MELO3C006697 are promising candidate genes for PM resistance in melon. The genetic architecture underlying this resistance in melon offers valuable insights for breeding programs, and the identified markers, especially chr06_indel_5 047 127, may enable practical applications for marker-assisted selection in developing PM-resistant melon varieties.

CO2 response screen in grass Brachypodium reveals the key role of a MAP kinase in CO2-triggered stomatal closure

Plants respond to increased CO2 concentrations through stomatal closure, which can contribute to increased water use efficiency. Grasses display faster stomatal responses than eudicots due to dumbbell-shaped guard cells flanked by subsidiary cells working in opposition. However, forward genetic screening for stomatal CO2 signal transduction mutants in grasses has yet to be reported. The grass model Brachypodium distachyon is closely related to agronomically important cereal crops, sharing largely collinear genomes. To gain insights into CO2 control mechanisms of stomatal movements in grasses, we developed an unbiased forward genetic screen with an EMS-mutagenized B. distachyon M5 generation population using infrared imaging to identify plants with altered leaf temperatures at elevated CO2. Among isolated mutants, a "chill1" mutant exhibited cooler leaf temperatures than wild-type Bd21-3 parent control plants after exposure to increased CO2. chill1 plants showed strongly impaired high CO2-induced stomatal closure despite retaining a robust abscisic acid-induced stomatal closing response. Through bulked segregant whole-genome sequencing analyses followed by analyses of further backcrossed F4 generation plants and generation and characterization of sodium azide and CRISPR-cas9 mutants, chill1 was mapped to a protein kinase, Mitogen-Activated Protein Kinase 5 (BdMPK5). The chill1 mutation impaired BdMPK5 protein-mediated CO2/HCO3- sensing together with the High Temperature 1 (HT1) Raf-like kinase in vitro. Furthermore, AlphaFold2-directed structural modeling predicted that the identified BdMPK5-D90N chill1 mutant residue is located at the interface of BdMPK5 with the BdHT1 Raf-like kinase. BdMPK5 is a key signaling component that mediates CO2-induced stomatal movements and is proposed to function as a component of the primary CO2 sensor in grasses.

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.

Fine-mapping of LrN3B on wheat chromosome arm 3BS, one of the two complementary genes for adult-plant leaf rust resistance

The common wheat line 4N0461 showed adult-plant resistance to leaf rust. 4N0461 was crossed with susceptible cultivars Nongda4503 and Shi4185 to map the causal resistance gene(s). Segregation of leaf rust response in F2 populations from both crosses was 9 resistant:7 susceptible, indicative of two complementary dominant resistance genes. The genes were located on chromosome arms 3BS and 4BL and temporarily named LrN3B and LrN4B, respectively. Subpopulations from 4N0461 × Nongda4503 with LrN3B segregating as a single allele were used to fine-map LrN3B locus. LrN3B was delineated in a genetic interval of 0.07 cM, corresponding to 106 kb based on the Chinese Spring reference genome (IWGSC RefSeq v1.1). Four genes were annotated in this region, among which TraesCS3B02G014800 and TraesCS3B02G014900 differed between resistant and susceptible genotypes, and both were required for LrN3B resistance in virus-induced gene silencing experiments. Diagnostic markers developed for checking the polymorphism of each candidate gene, can be used for marker-assisted selection in wheat breeding programs.

Identification of stripe rust resistance gene YrBDT in Chinese landrace wheat Baidatou using BSE-seq and BSR-seq

KEY MESSAGE

A new stripe rust resistance gene YrBDT in Chinese landrace wheat Baidatou was mapped to a 943.6-kb interval on chromosome arm 6DS and co-segregated with a marker CAPS3 developed from candidate gene TraesCS6D03G0027300. Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is a devastating foliar disease of wheat. Chinese landrace wheat Baidatou has shown high resistance to a broad spectrum of Pst races at both the seedling and adult-plant stages for decades in the Longnan region of Gansu province, a hot spot for stripe rust epidemics. Here, we report fine mapping and candidate gene analysis of stripe rust resistance gene YrBDT in Baidatou. Analysis of F1, F2 plants and F2:3 lines indicated that resistance in Baidatou to Pst race CYR31 was conferred by a single dominant gene, temporarily designated YrBDT. Bulked segregant exome capture sequencing (BSE-seq) analysis revealed 61 high-confidence polymorphic SNPs concentrated in a 5.4-Mb interval at the distal of chromosome arm 6DS. Several SNPs and InDels were also identified by genome mining of DNA sampled from the parents and contrasting bulks. The YrBDT locus was mapped to a 943.6-kb (4,658,322-5,601,880 bp) genomic region spanned by markers STS2 and STS3 based on IWGSC RefSeq v2.1, including five putative disease resistance genes. There was high collinearity of the target interval among Chinese Spring RefSeq v2.1, Ae. tauschii AL8/78 and Fielder genomes. The expression level of TraesCS6D03G0027300 showed significant association with Pst infection, and a gene-specific marker CAPS3 developed from TraesCS6D03G0027300 co-segregated with YrBDT suggesting this gene as a candidate of YrBDT. The resistance gene and flanking markers can be used in marker-assisted selection for improvement of stripe rust resistance.

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.

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

INTRODUCTION

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSION

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

Identification of Rmg11 in Tetraploid Wheat as a New Blast Resistance Gene with Tolerance to High Temperature

Wheat blast caused by Pyricularia oryzae pathotype Triticum has spread to Asia (Bangladesh) and Africa (Zambia) from the endemic region of South America. Wheat varieties with durable resistance are needed, but very limited resistance resources are currently available. After screening tetraploid wheat accessions, we found an exceptional accession St19 (Triticum dicoccum, KU-114). Primary leaves of St19 were resistant not only to Brazilian isolate Br48 (a carrier of Type eI of AVR-Rmg8) but also to Br48ΔA8, an AVR-Rmg8 disruptant of Br48, even at 30°C, suggesting that the resistance of St19 is tolerant to high temperature and controlled by a gene or genes other than Rmg8. When an F2 population derived from a cross between St19 and St30 (a susceptible accession of T. paleocolchicum, KU-191) was inoculated with Br48, resistant and susceptible seedlings segregated in a 3:1 ratio, indicating that resistance of St19 is conferred by a single gene. We designated this gene Rmg11. Molecular mapping revealed that the RMG11 locus is located on the short arm of chromosome 7A. Rmg11 is effective not only against other two Brazilian isolates (Br5 and Br116.5) but also against Bangladeshi isolates (T-108 and T-109) at the seedling stage. At the heading stage, lines containing Rmg11 were highly susceptible to the Bangladeshi isolates but moderately resistant to the Brazilian isolates. Stacking of Rmg11 with Rmg8 and the 2NS segment is highly recommended to achieve durable wheat blast resistance.

A pipeline for identification of causal mutations in barley identifies Xantha-j as the chlorophyll synthase gene

Thousands of barley (Hordeum vulgare L.) mutants have been isolated over the last century, and many are stored in gene banks across various countries. In the present work, we developed a pipeline to efficiently identify causal mutations in barley. The pipeline is also efficient for mutations located in centromeric regions. Through bulked segregant analyses using whole genome sequencing of pooled F2 seedlings, we mapped 2 mutations and identified a limited number of candidate genes. We applied the pipeline on F2 mapping populations made from xan-j.59 (unknown mutation) and xan-l.82 (previously known). The Xantha-j (xan-j) gene was identified as encoding chlorophyll synthase, which catalyzes the last step in the chlorophyll biosynthetic pathway: the addition of a phytol moiety to the propionate side chain of chlorophyllide. Key amino acid residues in the active site, including the binding sites of the isoprenoid and chlorophyllide substrates, were analyzed in an AlphaFold2-generated structural model of the barley chlorophyll synthase. Three allelic mutants, xan-j.19, xan-j.59, and xan-j.64, were characterized. While xan-j.19 is a 1 base pair deletion and xan-j.59 is a nonsense mutation, xan-j.64 causes an S212F substitution in chlorophyll synthase. Our analyses of xan-j.64 and treatment of growing barley with clomazone, an inhibitor of chloroplastic isoprenoid biosynthesis, suggest that binding of the isoprenoid substrate is a prerequisite for the stable maintenance of chlorophyll synthase in the plastid. We further suggest that chlorophyll synthase is a sensor for coordinating chlorophyll and isoprenoid biosynthesis.

Advances in genome sequencing and artificially induced mutation provides new avenues for cotton breeding

Cotton production faces challenges in fluctuating environmental conditions due to limited genetic variation in cultivated cotton species. To enhance the genetic diversity crucial for this primary fiber crop, it is essential to augment current germplasm resources. High-throughput sequencing has significantly impacted cotton functional genomics, enabling the creation of diverse mutant libraries and the identification of mutant functional genes and new germplasm resources. Artificial mutation, established through physical or chemical methods, stands as a highly efficient strategy to enrich cotton germplasm resources, yielding stable and high-quality raw materials. In this paper, we discuss the good foundation laid by high-throughput sequencing of cotton genome for mutant identification and functional genome, and focus on the construction methods of mutant libraries and diverse sequencing strategies based on mutants. In addition, the important functional genes identified by the cotton mutant library have greatly enriched the germplasm resources and promoted the development of functional genomes. Finally, an innovative strategy for constructing a cotton CRISPR mutant library was proposed, and the possibility of high-throughput screening of cotton mutants based on a UAV phenotyping platform was discussed. The aim of this review was to expand cotton germplasm resources, mine functional genes, and develop adaptable materials in a variety of complex environments.

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.

Fine mapping and identification of a Fusarium wilt resistance gene FwS1 in pea

KEY MESSAGE

A Fusarium wilt resistance gene FwS1 on pea chromosome 6 was identified and mapped to a 91.4 kb region by a comprehensive genomic-based approach, and the gene Psat6g003960 harboring NB-ARC domain was identified as the putative candidate gene. Pea Fusarium wilt, incited by Fusarium oxysporum f. sp. pisi (Fop), has always been a devastating disease that causes severe yield losses and economic damage in pea-growing regions worldwide. The utilization of pea cultivars carrying resistance gene is the most efficient approach for managing this disease. In order to finely map resistance gene, F2 populations were established through the cross between Shijiadacaiwan 1 (resistant) and Y4 (susceptible). The resistance genetic analysis indicated that the Fop resistance in Shijiadacaiwan 1 was governed by a single dominant gene, named FwS1. Based on the bulked segregant analysis sequencing analyses, the gene FwS1 was initially detected on chromosome 6 (i.e., linking group II, chr6LG2), and subsequent linkage mapping with 589 F2 individuals fine-mapped the gene FwS1 into a 91.4 kb region. The further functional annotation and haplotype analysis confirmed that the gene Psat6g003960, characterized by a NB-ARC (nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4) domain, was considered as the most promising candidate gene. The encoding amino acids were altered by a "T/C" single-nucleotide polymorphism (SNP) in the first exon of the Psat6g003960, and based on this SNP locus, the molecular marker A016180 was determined to be a diagnostic marker for FwS1 by validating its specificity in both pea accessions and genetic populations with different genetic backgrounds. The FwS1 with diagnostic KASP marker A016180 could facilitate marker-assisted selection in resistance pea breeding in pea. In addition, a comparison of the candidate gene Psat6g003960 in 74SN3B and SJ1 revealed the same sequences. This finding indicated that 74SN3B carried the candidate gene for FwS1, suggesting that FwS1 and Fwf may be closely linked or an identical resistant gene against Fusarium wilt.

Fine mapping of leaf delayed virescence gene dv4 in Triticum aestivum

Wheat (Triticum aestivum L.) is one of the most important crops worldwide, and its yield affects national food security. Wheat leaves are key photosynthetic organs where carbohydrates are synthesized for grain yield. Leaf colour mutants are ideal germplasm resources for molecular genetic studies of wheat chloroplast development, chlorophyll synthesis and photosynthesis. We obtained a wheat mutant delayed virescence 4 (dv4) from cultivar Guomai 301. The leaves of mutant dv4 were pale yellow at the seedling stage, golden yellow at the turning green stage, and they started to turn green at the jointing stage. Genetic analysis demonstrated that the yellow-leaf phenotype was controlled by a single recessive gene named as dv4. Gene dv4 was fine mapped in a 1.46 Mb region on chromosome 7DS by SSR and dCAPS marker assays. Three putative candidate genes were identified in this region. Because no leaf colour genes have been reported on wheat chromosome arm 7DS previously, dv4 is a novel leaf colour gene. The result facilitates map-based cloning of dv4 and provides information for the construction of a high-photosynthetic efficiency ideotype for improving wheat yield.

High-resolution bulked segregant analysis enables candidate gene identification for bacterial wilt resistance in Italian ryegrass (Lolium multiflorum Lam.)

Bacterial wilt, caused by Xanthomonas translucens pv. graminis (Xtg), is a serious disease of economically important forage grasses, including Italian ryegrass (Lolium multiflorum Lam.). A major QTL for resistance to Xtg was previously identified, but the precise location as well as the genetic factors underlying the resistance are yet to be determined. To this end, we applied a bulked segregant analysis (BSA) approach, using whole-genome deep sequencing of pools of the most resistant and most susceptible individuals of a large (n = 7484) biparental F2 population segregating for resistance to Xtg. Using chromosome-level genome assemblies as references, we were able to define a ~300 kb region highly associated with resistance on pseudo-chromosome 4. Further investigation of this region revealed multiple genes with a known role in disease resistance, including genes encoding for Pik2-like disease resistance proteins, cysteine-rich kinases, and RGA4- and RGA5-like disease resistance proteins. Investigation of allele frequencies in the pools and comparative genome analysis in the grandparents of the F2 population revealed that some of these genes contain variants with allele frequencies that correspond to the expected heterozygosity in the resistant grandparent. This study emphasizes the efficacy of combining BSA studies in very large populations with whole genome deep sequencing and high-quality genome assemblies to pinpoint regions associated with a binary trait of interest and accurately define a small set of candidate genes. Furthermore, markers identified in this region hold significant potential for marker-assisted breeding strategies to breed resistance to Xtg in Italian ryegrass cultivars more efficiently.

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

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