Potential locus W and candidate gene McPRR2 associated with pericarp pigment accumulation in bitter gourd (Momordica charantia L.) revealed via BSA-seq analysis

Pericarp color is a prominent agronomic trait that exerts a significant impact on consumer and breeder preferences. Genetic analysis has revealed that the pericarp color of bitter gourd is a quantitative trait. However, the underlying mechanism for this trait in bitter gourd remains largely unknown. In the present study, we employed bulked segregant analysis (BSA) to identify the candidate genes responsible for bitter gourd pericarp color (specifically, dark green versus white) within F2 segregation populations resulting from the crossing of B07 (dark green pericarp) and A06 (white pericarp). Through genomic variation, genetic mapping, and expression analysis, we identified a candidate gene named McPRR2, which was a homolog of Arabidopsis pseudo response regulator 2 (APRR2) encoded by LOC111023472. Sequence alignment of the candidate gene between the two parental lines revealed a 15-bp nucleotide insertion in the coding region of LOC111023472, leading to a premature stop codon and potentially causing a loss-of-function mutation. qRT-PCR analysis demonstrated that the expression of McPRR2 was significantly higher in B07 compared to A06, and it was primarily expressed in the immature fruit pericarp. Moreover, overexpression of McPRR2 in tomato could enhance the green color of immature fruit pericarp by increasing the chlorophyll content. Consequently, McPRR2 emerged as a strong candidate gene regulating the bitter gourd pericarp color by influencing chlorophyll accumulation. Finally, we developed a molecular marker linked to pericarp color, enabling the identification of genotypes in breeding populations. These findings provided valuable insights into the genetic improvement of bitter gourd pericarp color.

A 5.2-kb insertion in the coding sequence of PavSCPL, a serine carboxypeptidase-like enhances fruit firmness in Prunus avium

Fruit firmness is an important trait in sweet cherry breeding because it directly positively influences fruit transportability, storage and shelf life. However, the underlying genes responsible and the molecular mechanisms that control fruit firmness remain unknown. In this study, we identified a candidate gene, PavSCPL, encoding a serine carboxypeptidase-like protein with natural allelic variation, that controls fruit firmness in sweet cherry using map-based cloning and functionally characterized PavSCPL during sweet cherry fruit softening. Genetic analysis revealed that fruit firmness in the 'Rainier' × 'Summit' F1 population was controlled by a single dominant gene. Bulked segregant analysis combined with fine mapping narrowed the candidate gene to a 473-kb region (7418778-7 891 914 bp) on chromosome 6 which included 72 genes. The candidate gene PavSCPL, and a null allele harbouring a 5244-bp insertion in the second exon that completely inactivated PavSCPL expression and resulted in the extra-hard-flesh phenotype, were identified by RNA-sequencing analysis and gene cloning. Quantitative RT-PCR analysis revealed that the PavSCPL expression level was increased with fruit softening. Virus-induced gene silencing of PavSCPL enhanced fruit firmness and suppressed the activities of certain pectin-degrading enzymes in the fruit. In addition, we developed functional molecular markers for PavSCPL and the Pavscpl5.2-k allele that co-segregated with the fruit firmness trait. Overall, this research identified a crucial functional gene for fruit firmness. The results provide insights into the genetic control and molecular mechanism of the fruit firmness trait and present useful molecular markers for molecular-assisted breeding for fruit firmness in sweet cherry.

Rapidly mining candidate cotton drought resistance genes based on key indicators of drought resistance

BACKGROUND

Focusing on key indicators of drought resistance is highly important for quickly mining candidate genes related to drought resistance in cotton.

RESULTS

In the present study, drought resistance was identified in drought resistance-related RIL populations during the flowering and boll stages, and multiple traits were evaluated; these traits included three key indicators: plant height (PH), single boll weight (SBW) and transpiration rate (Tr). Based on these three key indicators, three groups of extreme mixing pools were constructed for BSA-seq. Based on the mapping interval of each trait, a total of 6.27 Mb QTL intervals were selected on chromosomes A13 (3.2 Mb), A10 (2.45 Mb) and A07 (0.62 Mb) as the focus of this study. Based on the annotation information and qRT‒PCR analysis, three key genes that may be involved in the drought stress response of cotton were screened: GhF6'H1, Gh3AT1 and GhPER55. qRT‒PCR analysis of parental and extreme germplasm materials revealed that the expression of these genes changed significantly under drought stress. Cotton VIGS experiments verified the important impact of key genes on cotton drought resistance.

CONCLUSIONS

This study focused on the key indicators of drought resistance, laying the foundation for the rapid mining of drought-resistant candidate genes in cotton and providing genetic resources for directed molecular breeding of drought resistance in cotton.

Identification of a major QTL and candidate genes analysis for branch angle in rapeseed (Brassica napus L.) using QTL-seq and RNA-seq

Introduction

Branching angle is an essential trait in determining the planting density of rapeseed (Brassica napus L.) and hence the yield per unit area. However, the mechanism of branching angle formation in rapeseed is not well understood.

Methods

In this study, two rapeseed germplasm with extreme branching angles were used to construct an F2 segregating population; then bulked segregant analysis sequencing (BSA-seq) and quantitative trait loci (QTL) mapping were utilized to localize branching anglerelated loci and combined with transcriptome sequencing (RNA-seq) and quantitative real-time PCR (qPCR) for candidate gene mining.

Results and discussion

A branching angle-associated quantitative trait loci (QTL) was mapped on chromosome C3 (C3: 1.54-2.65 Mb) by combining BSA-seq as well as traditional QTL mapping. A total of 54 genes had SNP/Indel variants within the QTL interval were identified. Further, RNA-seq of the two parents revealed that 12 of the 54 genes were differentially expressed between the two parents. Finally, we further validated the differentially expressed genes using qPCR and found that six of them presented consistent differential expression in all small branching angle samples and large branching angles, and thus were considered as candidate genes related to branching angles in rapeseed. Our results introduce new candidate genes for the regulation of branching angle formation in rapeseed, and provide an important reference for the subsequent exploration of its formation mechanism.

The QTL and Candidate Genes Regulating the Early Tillering Vigor Traits of Late-Season Rice in Double-Cropping Systems

Rice effective panicle is a major trait for grain yield and is affected by both the genetic tiller numbers and the early tillering vigor (ETV) traits to survive environmental adversities. The mechanism behind tiller bud formation has been well described, while the genes and the molecular mechanism underlying rice-regulating ETV traits are unclear. In this study, the candidate genes in regulating ETV traits have been sought by quantitative trait locus (QTL) mapping and bulk-segregation analysis by resequencing method (BSA-seq) conjoint analysis using rice backcross inbred line (BIL) populations, which were cultivated as late-season rice of double-cropping rice systems. By QTL mapping, seven QTLs were detected on chromosomes 1, 3, 4, and 9, with the logarithm of the odds (LOD) values ranging from 3.52 to 7.57 and explained 3.23% to 12.98% of the observed phenotypic variance. By BSA-seq analysis, seven QTLs on chromosomes 1, 2, 4, 5, 7, and 9 were identified using single-nucleotide polymorphism (SNP) and insertions/deletions (InDel) index algorithm and Euclidean distance (ED) algorithm. The overlapping QTL resulting from QTL mapping and BSA-seq analysis was shown in a 1.39 Mb interval on chromosome 4. In the overlap interval, six genes, including the functional unknown genes Os04g0455650, Os04g0470901, Os04g0500600, and ethylene-insensitive 3 (Os04g0456900), sialyltransferase family domain containing protein (Os04g0506800), and ATOZI1 (Os04g0497300), showed the differential expression between ETV rice lines and late tillering vigor (LTV) rice lines and have a missense base mutation in the genomic DNA sequences of the parents. We speculate that the six genes are the candidate genes regulating the ETV trait in rice, which provides a research basis for revealing the molecular mechanism behind the ETV traits in rice.

Unveiling Vacuolar H+-ATPase Subunit a as the Primary Target of the Pyridinylmethyl-Benzamide Fungicide, Fluopicolide

An estimated 240 fungicides are presently in use, but the direct targets for the majority remain elusive, constraining fungicide development and efficient resistance monitoring. In this study, we found that Pcα-actinin knockout did not influence the sensitivity of Phytophthora capsici to fluopicolide, which is a notable oomycete inhibitor. Using a combination of Bulk Segregant Analysis Sequencing and Drug Affinity Responsive Target Stability (DARTS) assays, the vacuolar H+-ATPase subunit a (PcVHA-a) was pinpointed as the target protein of fluopicolide. We also confirmed four distinct point mutations in PcVHA-a responsible for fluopicolide resistance in P. capsici through site-directed mutagenesis. Molecular docking, ATPase activity assays, and a DARTS assay suggested a fluopicolide-PcVHA-a interaction. Sequence analysis and further molecular docking validated the specificity of fluopicolide for oomycetes or fish. These findings support the claim that PcVHA-a is the target of fluopicolide, proposing vacuolar H+-ATPase as a promising target for novel fungicide development.

Integrated BSA-seq and RNA-seq analysis to identify candidate genes associated with nitrogen utilization efficiency (NUtE) in rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) is one of the important oil crops, with a high demand for nitrogen (N). It is essential to explore the potential of rapeseed to improve nitrogen utilization efficiency (NUtE). Rapeseed is an allotetraploid crop with a relatively large and complex genome, and there are few studies on the mapping of genes related to NUtE regulation. In this study, we used the combination of bulk segregant analysis sequencing (BSA-Seq) and RNA sequencing (RNA-Seq) to analyze the N-efficient genotype 'Zheyou 18' and N-inefficient genotype 'Sollux', to identify the genetic regulatory mechanisms. Several candidate genes were screened, such as the high-affinity nitrate transporter gene NRT2.1 (BnaC08g43370D) and the abscisic acid (ABA) signal transduction-related genes (BnaC02g14540D, BnaA03g20760D, and BnaA05g01330D). BnaA05g01330D was annotated as ABA-INDUCIBLE bHLH-TYPE TRANSCRIPTION FACTOR (AIB/bHLH17), which was highly expressed in the root. The results showed that the primary root length of the ataib mutant was significantly longer than that of the wild type under low N conditions. Overexpression of BnaA5.AIB could reduce the NUtE under low N levels in Arabidopsis (Arabidopsis thaliana). Candidate genes identified in this study may be involved in the regulation of NUtE in rapeseed, and new functions of AIB in orchestrating N uptake and utilization have been revealed. It is indicated that BnaA5.AIB may be the key factor that links ABA to N signaling and a negative regulator of NUtE. It will provide a theoretical basis and application prospect for resource conservation, environmental protection, and sustainable agricultural development.

TBtools-II: A "one for all, all for one" bioinformatics platform for biological big-data mining

Since the official release of the stand-alone bioinformatics toolkit TBtools in 2020, its superior functionality in data analysis has been demonstrated by its widespread adoption by many thousands of users and references in more than 5000 academic articles. Now, TBtools is a commonly used tool in biological laboratories. Over the past 3 years, thanks to invaluable feedback and suggestions from numerous users, we have optimized and expanded the functionality of the toolkit, leading to the development of an upgraded version-TBtools-II. In this upgrade, we have incorporated over 100 new features, such as those for comparative genomics analysis, phylogenetic analysis, and data visualization. Meanwhile, to better meet the increasing needs of personalized data analysis, we have launched the plugin mode, which enables users to develop their own plugins and manage their selection, installation, and removal according to individual needs. To date, the plugin store has amassed over 50 plugins, with more than half of them being independently developed and contributed by TBtools users. These plugins offer a range of data analysis options including co-expression network analysis, single-cell data analysis, and bulked segregant analysis sequencing data analysis. Overall, TBtools is now transforming from a stand-alone software to a comprehensive bioinformatics platform of a vibrant and cooperative community in which users are also developers and contributors. By promoting the theme "one for all, all for one", we believe that TBtools-II will greatly benefit more biological researchers in this big-data era.

Physiological Analysis and Genetic Mapping of Short Hypocotyl Trait in Brassica napus L

Hypocotyl length is a botanical trait that affects the cold tolerance of Brassica napus L. (B. napus). In this study, we constructed an F2 segregating population using the cold-resistant short hypocotyl variety '16VHNTS158' and the cold-sensitive long hypocotyl variety 'Tianyou 2288' as the parents, and BSA-seq was employed to identify candidate genes for hypocotyl length in B. napus. The results of parental differences showed that the average hypocotyl lengths of '16VHNTS158' and 'Tianyou 2288' were 0.41 cm and 0.77 cm at the 5~6 leaf stage, respectively, after different low-temperature treatments, and '16VHNTS158' exhibited lower relative ion leakage rates compared to 'Tianyou 2288'. The contents of indole acetic acid (IAA), gibberellin (GA), and brassinosteroid (BR) in hypocotyls of '16VHNTS158' and 'Tianyou 2288' increased with decreasing temperatures, but the IAA and GA contents were significantly higher than those of 'Tianyou 2288', and the BR content was lower than that of 'Tianyou 2288'. The genetic analysis results indicate that the genetic model for hypocotyl length follows the 2MG-A model. By using SSR molecular markers, a QTL locus associated with hypocotyl length was identified on chromosome C04. The additive effect value of this locus was 0.025, and it accounted for 2.5% of the phenotypic variation. BSA-Seq further localized the major effect QTL locus on chromosome C04, associating it with 41 genomic regions. The total length of this region was 1.06 Mb. Within this region, a total of 20 non-synonymous mutation genes were identified between the parents, and 26 non-synonymous mutation genes were found within the pooled samples. In the reference genome of B. napus, this region was annotated with 24 candidate genes. These annotated genes are predominantly enriched in four pathways: DNA replication, nucleotide excision repair, plant hormone signal transduction, and mismatch repair. The findings of this study provide a theoretical basis for cloning genes related to hypocotyl length in winter rapeseed and their utilization in breeding.

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

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

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

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

Fine Mapping and Identification of a Candidate Gene for the Glossy Green Trait in Cabbage (Brassica oleracea var. capitata)

In higher plants, cuticular wax deposited on the surface of epidermal cells plays an important role in protecting the plant from biotic and abiotic stresses; however, the molecular mechanism of cuticular wax production is not completely understood. In this study, we identified a glossy green mutant (98-1030gl) from the glaucous cabbage inbred line 98-1030. Scanning electron microscopy indicated that the amount of leaf cuticular wax significantly decreased in 98-1030gl. Genetic analysis showed that the glossy green trait was controlled by a single recessive gene. Bulked segregant analysis coupled with whole genome sequencing revealed that the candidate gene for the glossy green trait was located at 13,860,000-25,070,000 bp (11.21 Mb) on Chromosome 5. Based on the resequencing data of two parents and the F2 population, insertion-deletion markers were developed and used to reduce the candidate mapping region. The candidate gene (Bol026949) was then mapped in a 50.97 kb interval. Bol026949 belongs to the Agenet/Tudor domain protein family, whose members are predicted to be involved in chromatin remodeling and RNA transcription. Sequence analysis showed that a single nucleotide polymorphism mutation (C → G) in the second exon of Bol026949 could result in the premature termination of its protein translation in 98-1030gl. Phylogenetic analysis showed that Bol026949 is relatively conserved in cruciferous plants. Transcriptome profiling indicated that Bol026949 might participate in cuticular wax production by regulating the transcript levels of genes involved in the post-translational cellular process and phytohormone signaling. Our findings provide an important clue for dissecting the regulatory mechanisms of cuticular wax production in cruciferous crops.

Genetic basis of the oil biosynthesis in ultra-high-oil maize grains with an oil content exceeding 20

Vegetable oil is an important part of the human diet and has multiple industrial uses. The rapid increase in vegetable oil consumption has necessitated the development of viable methods for optimizing the oil content of plants. The key genes regulating the biosynthesis of maize grain oil remain mostly uncharacterized. In this study, by analyzing oil contents and performing bulked segregant RNA sequencing and mapping analyses, we determined that su1 and sh2-R mediate the shrinkage of ultra-high-oil maize grains and contribute to the increase in the grain oil content. Functional kompetitive allele-specific PCR (KASP) markers developed for su1 and sh2-R detected su1su1Sh2Sh2, Su1Su1sh2sh2, and su1su1sh2sh2 mutants among 183 sweet maize inbred lines. An RNA sequencing (RNA-seq) analysis indicated that genes differentially expressed between two conventional sweet maize lines and two ultra-high-oil maize lines were significantly associated with linoleic acid metabolism, cyanoamino acid metabolism, glutathione metabolism, alanine, aspartate, and glutamate metabolism, and nitrogen metabolism. A bulk segregant analysis and sequencing (BSA-seq) analysis identified another 88 genomic intervals related to grain oil content, 16 of which overlapped previously reported maize grain oil-related QTLs. The combined analysis of BSA-seq and RNA-seq data enabled the identification of candidate genes. The KASP markers for GRMZM2G176998 (putative WD40-like beta propeller repeat family protein), GRMZM2G021339 (homeobox-transcription factor 115), and GRMZM2G167438 (3-ketoacyl-CoA synthase) were significantly related to maize grain oil content. Another candidate gene, GRMZM2G099802 (GDSL-like lipase/acylhydrolase), catalyzes the final step of the triacylglycerol synthesis pathway and was expressed at significantly higher levels in the two ultra-high-oil maize lines than in the two conventional sweet maize lines. These novel findings will help clarify the genetic basis of the increased oil production in ultra-high-oil maize lines with grain oil contents exceeding 20%. The KASP markers developed in this study may be useful for breeding new high-oil sweet maize varieties.

Map-based cloning of the APRR2 gene controlling green stigma in bitter gourd (Momordica charantia)

Bitter gourd is an economically important vegetable and medicinal crop distinguished by its bitter fruits. Its stigma color is widely used to assess the distinctiveness, uniformity, and stability of bitter gourd varieties. However, limited researches have been dedicated to genetic basis of its stigma color. In this study, we employed bulked segregant analysis (BSA) sequencing to identify a single dominant locus McSTC1 located on pseudochromosome 6 through genetic mapping of an F₂ population (n =241) derived from the cross between green and yellow stigma parental lines. An F₂-derived F₃ segregation population (n = 847) was further adopted for fine mapping, which delimited the McSTC1 locus to a 13.87 kb region containing one predicted gene McAPRR2 (Mc06g1638), a homolog of the Arabidopsis two-component response regulator-like gene AtAPRR2. Sequence alignment analysis of McAPRR2 revealed that a 15 bp insertion at exon 9 results in a truncated GLK domain of its encoded protein, which existed in 19 bitter gourd varieties with yellow stigma. A genome-wide synteny search of the bitter gourd McAPRR2 genes in Cucurbitaceae family revealed its close relationship with other cucurbits APRR2 genes that are corresponding to white or light green fruit skin. Our findings provide insights into the molecular marker-assisted breeding of bitter gourd stigma color and the mechanism of gene regulation for stigma color.

Genetic mapping of AhVt1, a novel genetic locus that confers the variegated testa color in cultivated peanut (Arachis hypogaea L.) and its utilization for marker-assisted selection

INTRODUCTION

Peanut (Arachis hypogaea L.) is an important cash crop worldwide. Compared with the ordinary peanut with pure pink testa, peanut with variegated testa color has attractive appearance and a higher market value. In addition, the variegated testa represents a distinct regulation pattern of anthocyanin accumulation in integument cells.

METHODS

In order to identify the genetic locus underlying variegated testa color in peanut, two populations were constructed from the crosses between Fuhua 8 (pure-pink testa) and Wucai (red on white variegated testa), Quanhonghua 1 (pure-red testa) and Wucai, respectively. Genetic analysis and bulked sergeant analysis sequencing were applied to detect and identify the genetic locus for variegated testa color. Marker-assisted selection was used to develop new variegated testa peanut lines.

RESULTS

As a result, all the seeds harvested from the F1 individuals of both populations showed the variegated testa type with white trace. Genetic analysis revealed that the pigmentation of colored region in red on white variegated testa was controlled by a previous reported gene AhRt1, while the formation of white region (un-pigmented region) in variegated testa was controlled by another single genetic locus. This locus, named as AhVt1 (Arachis hypogaea Variegated Testa 1), was preliminary mapped on chromosome 08 through bulked sergeant analysis sequencing. Using a secondary mapping population derived from the cross between Fuhua 8 and Wucai, AhVt1 was further mapped to a 1.89-Mb genomic interval by linkage analysis, and several potential genes associated with the uneven distribution of anthocyanin, such as MADS-box, MYB, and Chalcone synthase-like protein, were harbored in the region. Moreover, the molecular markers closely linked to the AhVt1 were developed, and the new variegated testa peanut lines were obtained with the help of marker-assisted selection.

CONCLUSION

Our findings will accelerate the breeding program for developing new peanut varieties with "colorful" testa colors and laid a foundation for map-based cloning of gene responsible for variegated testa.

New Genes Identified as Modulating Salt Tolerance in Maize Seedlings Using the Combination of Transcriptome Analysis and BSA

(1) Background: Salt stress is an abiotic factor that limits maize yield and quality. A highly salt-tolerance inbred AS5 and a salt-sensitive inbred NX420 collected from Ningxia Province, China, were used to identify new genes for modulating salt resistance in maize. (2) Methods: To understand the different molecular bases of salt tolerance in AS5 and NX420, we performed BSA-seq using an F2 population for two extreme bulks derived from the cross between AS5 and NX420. Transcriptomic analysis was also conducted for AS5 and NX420 at the seedling stage after treatment with 150 mM of NaCl for 14 days. (3) Results: AS5 had a higher biomass and lower Na⁺ content than NX420 in the seedling stage after treatment with 150 mM NaCl for 14 days. One hundred and six candidate regions for salt tolerance were mapped on all of the chromosomes through BSA-seq using F2 in an extreme population. Based on the polymorphisms identified between both parents, we detected 77 genes. A large number of differentially expressed genes (DEGs) at the seedling stage under salt stress between these two inbred lines were detected using transcriptome sequencing. GO analysis indicated that 925 and 686 genes were significantly enriched in the integral component of the membrane of AS5 and NX420, respectively. Among these results, two and four DEGs were identified as overlapping in these two inbred lines using BSA-seq and transcriptomic analysis, respectively. Two genes (Zm00001d053925 and Zm00001d037181) were detected in both AS5 and NX420; the transcription level of Zm00001d053925 was induced to be significantly higher in AS5 than in NX420 (41.99 times versus 6.06 times after 150 mM of NaCl treatment for 48 h), while the expression of Zm00001d037181 showed no significant difference upon salt treatment in both lines. The functional annotation of the new candidate genes showed that it was an unknown function protein. (4) Conclusions: Zm00001d053925 is a new functional gene responding to salt stress in the seedling stage, which provides an important genetic resource for salt-tolerant maize breeding.

Identification and characterization analysis of candidate genes controlling mushroom leaf development in Chinese kale by BSA-seq

Mushroom leaves (MLs) are malformed leaves that develop from the leaf veins in some of Chinese kale genotypes. To study the genetic model and molecular mechanism of ML development in Chinese kale, the F2 segregation population was constructed by two inbred lines, genotype Boc52 with ML and genotype Boc55 with normal leaves (NL). In the present study, we have identified for the first time that the development of mushroom leaves may be affected by the change of adaxial-abaxial polarity of leaves. Examination of the phenotypes of F1 and F2 segregation populations suggested that ML development is controlled by two dominant major genes inherited independently. BSA-seq analysis showed that a major quantitative trait locus (QTL) qML4.1 that controls ML development is located within 7.4 Mb on chromosome kC4. The candidate region was further narrowed to 255 kb by linkage analysis combined with insertion/deletion (InDel) markers, and 37 genes were predicted in this region. According to the expression and annotation analysis, a B3 domain-containing transcription factor NGA1-like gene, BocNGA1, was identified as a key candidate gene for controlling ML development in Chinese kale. Fifteen single nucleotide polymorphisms (SNPs) were found in coding sequences and 21 SNPs and 3 InDels found in the promoter sequences of BocNGA1 from the genotype Boc52 with ML. The expression levels of BocNGA1 in ML genotypes are significantly lower than in the NL genotypes, which suggests that BocNGA1 may act as a negative regulator for ML genesis in Chinese kale. This study provides a new foundation for Chinese kale breeding and for the study of the molecular mechanism of plant leaf differentiation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01364-6.

QTL mapping and BSA-seq map a major QTL for the node of the first fruiting branch in cotton

Understanding the genetic basis of the node of the first fruiting branch (NFFB) improves early-maturity cotton breeding. Here we report QTL mapping on 200 F₂ plants and derivative F_2:3 and F_2:4 populations by genotyping by sequencing (GBS). BC₁F₂ population was constructed by backcrossing one F_2:4 line with the maternal parent JF914 and used for BSA-seq for further QTL mapping. A total of 1,305,642 SNPs were developed between the parents by GBS, and 2,907,790 SNPs were detected by BSA-seq. A high-density genetic map was constructed containing 11,488 SNPs and spanning 4,202.12 cM in length. A total of 13 QTL were mapped in the 3 tested populations. JF914 conferred favorable alleles for 11 QTL, and JF173 conferred favorable alleles for the other 2 QTL. Two stable QTL were repeatedly mapped in F_2:3 and F_2:4, including qNFFB-D3-1 and qNFFB-D6-1. Only qNFFB-D3-1 contributed more than 10% of the phenotypic variation. This QTL covered about 24.7 Mb (17,130,008-41,839,226 bp) on chromosome D3. Two regions on D3 (41,779,195-41,836,120 bp, 41,836,768-41,872,287 bp) were found by BSA-seq and covered about 92.4 Kb. This 92.4 Kb region overlapped with the stable QTL qNFFB-D3-1 and contained 8 annotated genes. By qRT-PCR, Ghir_D03G012430 showed a lower expression level from the 1- to 2-leaf stage and a higher expression level from the 3- to 6-leaf stage in the buds of JF173 than that of JF914. Ghir_D03G012390 reached the highest level at the 3- and 5-leaf stages in the buds of JF173 and JF914, respectively. As JF173 has lower NFFB and more early maturity than JF914, these two genes might be important in cell division and differentiation during NFFB formation in the seedling stage. The results of this study will facilitate a better understanding of the genetic basis of NFFB and benefit cotton molecular breeding for improving earliness traits.

Integrating BSA-Seq with RNA-Seq Reveals a Novel Fasciated Ear5 Mutant in Maize

Increasing grain yield is required to meet the rapidly expanding demands for food, feed, and fuel. Inflorescence meristems are central to plant growth and development. However, the question concerning whether inflorescence development can be regulated to improve grain yield remains unclear. Here, we describe a naturally occurring single recessive mutation called fea5 that can increase grain yield in maize. Using bulk segregant analysis sequencing (BSA-seq), the candidate region was initially mapped to a large region on chromosome 4 (4.68 Mb-11.26 Mb). Transcriptome sequencing (RNA-seq) revealed a total of 1246 differentially expressed genes (DEGs), of which 835 were up-regulated and 411 were down-regulated. Further analysis revealed the enrichment of DEGs in phytohormone signal transduction. Consistently, phytohormone profiling indicated that auxin (IAA), jasmonic acid (JA), ethylene (ETH), and cytokinin (CK) levels increased significantly, whereas the gibberellin (GA) level decreased significantly in fea5. By integrating BSA-seq with RNA-seq, we identified Zm00001d048841 as the most likely candidate gene. Our results provide valuable insight into this new germplasm resource and the molecular mechanism underlying fasciated ears that produce a higher kernel row number in maize.

Mapping of a QTL associated with sucrose content in peanut kernels using BSA-seq

As an important factor affecting the edible quality of peanut kernels, sucrose content is a complex quantitative trait regulated by multiple factors. In this study, an F₂ segregating population and a recombinant inbred line (RIL) population, derived from a cross between the high sucrose content variety Jihuatian 1 and the low sucrose content line PI478819, were used as materials to map a quantitative trait locus (QTL) associated with sucrose content in peanut kernels. Four QTLs were initially located on chromosomes A03 and A06 based on BSA-seq technology, and multiple kompetitive allele-specific PCR markers were developed based on single-nucleotide polymorphisms (SNPs) in the intervals. The markers were genotyped in the RIL population and finely mapped to a stable QTL, qSUCA06, located on chromosome A06 within a 0.29-Mb physical genomic interval (112367085-112662675 bp), which accounted for 31.95%-41.05% of the phenotypic variance explained. SNP and insertion/deletion annotations were performed on genes in the candidate interval, and having screened out those genes with mutations in exons, candidate genes were verified by qRT-PCR. The results revealed that Arahy.Y2LWD9 may be the main gene regulating sucrose content. The QTL identified in this study will not only contribute to marker-assisted breeding for improvement of peanut sucrose content but also paves the way for identifying gene function.

The AP2 Transcription Factor BrSHINE3 Regulates Wax Accumulation in Nonheading Chinese Cabbage

Wax is an acellular structural substance attached to the surface of plant tissues. It forms a protective barrier on the epidermis of plants and plays an important role in resisting abiotic and biotic stresses. In this paper, nonheading Chinese cabbage varieties with and without wax powder were observed using scanning electron microscopy, and the surface of waxy plants was covered with a layer of densely arranged waxy crystals, thus differentiating them from the surface of waxless plants. A genetic analysis showed that wax powder formation in nonheading Chinese cabbage was controlled by a pair of dominant genes. A preliminary bulked segregant analysis sequencing (BSA-seq) assay showed that one gene was located at the end of chromosome A09. Within this interval, we identified BraA09000626, encoding an AP2 transcription factor homologous to Arabidopsis AtSHINE3, and we named it BrSHINE3. By comparing the CDS of the gene in the two parental plants, a 35 bp deletion in the BrSHINE3 gene of waxless plants resulted in a frameshift mutation. Tissue analysis showed that BrSHINE3 was expressed at significantly higher levels in waxy plant rosette stage petioles and bolting stage stems than in the tissues of waxless plants. We speculate that this deletion in BrSHINE3 bases in the waxless material may inhibit wax synthesis. The overexpression of BrSHINE3 in Arabidopsis induced the accumulation of wax on the stem surface, indicating that BrSHINE3 is a key gene that regulates the formation of wax powder in nonheading Chinese cabbage. The analysis of the subcellular localization showed that BrSHINE3 is mainly located in the nucleus and chloroplast of tobacco leaves, suggesting that the gene may function as a transcription factor. Subsequent transcriptome analysis of the homology of BrSHINE3 downstream genes in nonheading Chinese cabbage showed that these genes were downregulated in waxless materials. These findings provide a basis for a better understanding of the nonheading Chinese cabbage epidermal wax synthesis pathway and provide important information for the molecular-assisted breeding of nonheading Chinese cabbage.

Slym1 control the color etiolation of leaves by facilitating the decomposition of chlorophyll in tomato

Photosynthesis, as an important biological process of plants, produces organic substances for plant growth and development. Although the molecular mechanisms of photosynthesis had been well investigated, the relationship between chlorophyll synthesis and photosynthesis remains largely unknown. The leaf-color mutant was an ideal material for studying photosynthesis and chlorophyll synthesis, which had been seldom investigated in tomato. Here, we obtained a yellow leaf tomato mutant ym (The mutant plants from the line of zs4) in field. Transmission electron microscopy (TEM) and photosynthetic parameters results demonstrated that chloroplast's structure was obviously destroyed and photosynthetic capacity gets weak. The mutant was hybridized with the control to construct the F₂ segregation population for sequencing. Slym1 gene, controlling yellow mutant trait, was identified using Bulked Segregation Analysis. Slym1 was up-regulated in the mutant and Slym1 was located in the nucleus. The genes associated with photosynthesis and chlorophyll synthesis were down-regulated in Slym1-OE transgenic tomato plants. The results suggested that Slym1 negatively regulate photosynthesis. Photosynthetic pigment synthesis related genes HEMA, HEMB1, CHLG and CAO were up-regulated in Slym1 silencing plants. The redundant Slym1 binding the intermediate proteins MP resulting in hindering the interaction between MP and HY5 due to the Slym1 with a high expression level in ym mutant, lead to lots of the HY5 with unbound state accumulates in cells, that could accelerate the decomposition of chlorophyll. Therefore, the yellow leaf-color mutant ym could be used as an ideal material for further exploring the relationship between leaf color mutant and photosynthesis and the specific mechanism.

Identification and fine-mapping of a major QTL (PH1.1) conferring plant height in broomcorn millet (Panicum miliaceum)

The plant height of broomcorn millet (Panicum miliaceum) is a significant agronomic trait that is closely related to its plant architecture, lodging resistance, and final yield. However, the genes underlying the regulation of plant height in broomcorn millet are rarely reported. Here, an F₂ population derived from a cross between a normal variety, "Longmi12," and a dwarf mutant, "Zhang778," was constructed. Genetic analysis for the F₂ and F_2:3 populations revealed that the plant height was controlled by more than one locus. A major quantitative trait locus (QTL), PH1.1, was preliminarily identified in chromosome 1 using bulked segregant analysis sequencing (BSA-seq). PH1.1 was fine-mapped to a 109-kb genomic region with 15 genes using a high-density map. Among them, longmi011482 and longmi011489, containing nonsynonymous variations in their coding regions, and longmi011496, covering multiple insertion/deletion sequences in the promoter regions, may be possible candidate genes for PH1.1. Three diagnostic markers closely linked to PH1.1 were developed to validate the PH1.1 region in broomcorn millet germplasm. These findings laid the foundation for further understanding of the molecular mechanism of plant height regulation in broomcorn millet and are also beneficial to the breeding program for developing new varieties with optimal height.

Fine mapping of the gene controlling the weeping trait of Prunus persica and its uses for MAS in progenies

BACKGROUND

Fruit tree yield and fruit quality are affected by the tree's growth type, and branching angle is an important agronomic trait of fruit trees, which largely determines the crown structure. The weeping type of peach tree shows good ventilation and light transmission; therefore, it is commonly cultivated. However, there is no molecular marker closely linked with peach weeping traits for target gene screening and assisted breeding.

RESULTS

First, we confirmed that the peach weeping trait is a recessive trait controlled by a single gene by constructing segregating populations. Based on BSA-seq, we mapped the gene controlling this trait within 159 kb of physical distance on chromosome 3. We found a 35 bp deletion in the candidate area in standard type, which was not lacking in weeping type. For histological assessments, different types of branches were sliced and examined, showing fiber bundles in the secondary xylem of ordinary branches but not in weeping branches.

CONCLUSIONS

This study established a molecular marker that is firmly linked to weeping trait. This marker can be used for the selection of parents in the breeding process and the early screening of hybrid offspring to shorten the breeding cycle. Moreover, we preliminary explored histological differences between growth types. These results lay the groundwork for a better understanding of the weeping growth habit of peach trees.

Flavonoid Metabolic Profiles and Gene Mapping of Rice (Oryza sativa L.) Purple Gradient Grain Hulls

Rice (Oryza sativa L.) grain hull color is an easily observable trait and regarded as a crucial morphological marker in rice breeding. Here, a purple gradient grain hull mutant (pg) was found from natural mutations of a straw-white grain hull rice variety IARI 6184B (Orzya sativa L. subsp. indica). The color of the mutant grain hulls changed from straw-white to pink, then purple, and finally brownish-yellow. Ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) identified 217 flavonoids, including 18 anthocyanins, among which cyanidin O-syringic acid had the highest concentration in pink (66.2 × 10⁶) and purple (68.0 × 10⁶) grain hulls. The relative contents of hesperetin O-malonyl-hexoside, apigenin derivatives, genistein derivatives, and kaempferol 3-O derivatives were consistently downregulated during pg grain hull development. Conversely, 12 anthocyanins were upregulated in colored hulls, and cyanidin 3-O-malonylhexoside was abundant only in pink and purple grain hulls. Moreover, the candidate gene was mapped into a 1.38 Mb region on chromosome 4 through bulked segregant analysis based on deep sequencing (BSA-seq) and gene mapping approaches. These results increased our understanding of anthocyanin biosynthesis in rice grains, helping rice breeders to select new rice varieties with desirable grain traits.

Identification of candidate genes and clarification of the maintenance of the green pericarp of weedy rice grains

The weedy rice (Oryza sativa f. spontanea) pericarp has diverse colors (e.g., purple, red, light-red, and white). However, research on pericarp colors has focused on red and purple, but not green. Unlike many other common weedy rice resources, LM8 has a green pericarp at maturity. In this study, the coloration of the LM8 pericarp was evaluated at the cellular and genetic levels. First, an examination of their ultrastructure indicated that LM8 chloroplasts were normal regarding plastid development and they contained many plastoglobules from the early immature stage to maturity. Analyses of transcriptome profiles and differentially expressed genes revealed that most chlorophyll (Chl) degradation-related genes in LM8 were expressed at lower levels than Chl a/b cycle-related genes in mature pericarps, suggesting that the green LM8 pericarp was associated with inhibited Chl degradation in intact chloroplasts. Second, the F₂ generation derived from a cross between LM8 (green pericarp) and SLG (white pericarp) had a pericarp color segregation ratio of 9:3:4 (green:brown:white). The bulked segregant analysis of the F₂ populations resulted in the identification of 12 known genes in the chromosome 3 and 4 hotspot regions as candidate genes related to Chl metabolism in the rice pericarp. The RNA-seq and sqRT-PCR assays indicated that the expression of the Chl a/b cycle-related structural gene DVR (encoding divinyl reductase) was sharply up-regulated. Moreover, genes encoding magnesium-chelatase subunit D and the light-harvesting Chl a/b-binding protein were transcriptionally active in the fully ripened dry pericarp. Regarding the ethylene signal transduction pathway, the CTR (encoding an ethylene-responsive protein kinase) and ERF (encoding an ethylene-responsive factor) genes expression profiles were determined. The findings of this study highlight the regulatory roles of Chl biosynthesis- and degradation-related genes influencing Chl accumulation during the maturation of the LM8 pericarp.

GWAS and bulked segregant analysis reveal the Loci controlling growth habit-related traits in cultivated Peanut (Arachis hypogaea L.)

Background

Peanut (Arachis hypogaea L.) is a grain legume crop that originated from South America and is now grown around the world. Peanut growth habit affects the variety’s adaptability, planting patterns, mechanized harvesting, disease resistance, and yield. The objective of this study was to map the quantitative trait locus (QTL) associated with peanut growth habit-related traits by combining the genome-wide association analysis (GWAS) and bulked segregant analysis sequencing (BSA-seq) methods.

Results

GWAS was performed with 17,223 single nucleotide polymorphisms (SNPs) in 103 accessions of the U.S. mini core collection genotyped using an Affymetrix version 2.0 SNP array. With a total of 12,342 high-quality polymorphic SNPs, the 90 suggestive and significant SNPs associated with lateral branch angle (LBA), main stem height (MSH), lateral branch height (LBL), extent radius (ER), and the index of plant type (IOPT) were identified. These SNPs were distributed among 15 chromosomes. A total of 597 associated candidate genes may have important roles in biological processes, hormone signaling, growth, and development. BSA-seq coupled with specific length amplified fragment sequencing (SLAF-seq) method was used to find the association with LBA, an important trait of the peanut growth habit. A 4.08 Mb genomic region on B05 was associated with LBA. Based on the linkage disequilibrium (LD) decay distance, we narrowed down and confirmed the region within the 160 kb region (144,193,467–144,513,467) on B05. Four candidate genes in this region were involved in plant growth. The expression levels of Araip.E64SW detected by qRT-PCR showed significant difference between ‘Jihua 5’ and ‘M130’.

Conclusions

In this study, the SNP (AX-147,251,085 and AX-144,353,467) associated with LBA by GWAS was overlapped with the results in BSA-seq through combined analysis of GWAS and BSA-seq. Based on LD decay distance, the genome range related to LBA on B05 was shortened to 144,193,467–144,513,467. Three candidate genes related to F-box family proteins (Araip.E64SW, Araip.YG1LK, and Araip.JJ6RA) and one candidate gene related to PPP family proteins (Araip.YU281) may be involved in plant growth and development in this genome region. The expression analysis revealed that Araip.E64SW was involved in peanut growth habits. These candidate genes will provide molecular targets in marker-assisted selection for peanut growth habits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08640-3.

BSA-Seq and Fine Linkage Mapping for the Identification of a Novel Locus (qPH9) for Mature Plant Height in Rice (Oryza sativa)

BACKGROUND

Plant height is a key factor in the determination of rice yield since excessive height can easily cause lodging and reduce yield. Therefore, the identification and analysis of plant height-related genes to elucidate their physiological, biochemical, and molecular mechanisms have significant implications for rice breeding and production.

RESULTS

High-throughput quantitative trait locus (QTL) sequencing analysis of a 638-individual F2:3 mapping population resulted in the identification of a novel height-related QTL (qPH9), which was mapped to a 2.02-Mb region of Chromosome 9. Local QTL mapping, which was conducted using 13 single nucleotide polymorphism (SNP)-based Kompetitive allele-specific PCR (KASP) markers for the qPH9 region, and traditional linkage analysis, facilitated the localization of qPH9 to a 126-kb region that contained 15 genes. Subsequent haplotype and sequence analyses indicated that OsPH9 was the most probable candidate gene for plant height at this locus, and functional analysis of osph9 CRISPR/Cas9-generated OsPH9 knockout mutants supported this conclusion.

CONCLUSION

OsPH9 was identified as a novel regulatory gene associated with plant height in rice, along with a height-reducing allele in 'Dongfu-114' rice, thereby representing an important molecular target for rice improvement. The findings of the present study are expected to spur the investigation of genetic mechanisms underlying rice plant height and further the improvement of rice plant height through marker-assisted selection.

Genome-wide association analysis reveals a novel QTL CsPC1 for pericarp color in cucumber

BACKGROUND

Cucumber is an important melon crop in the world, with different pericarp colors. However, the candidate genes and the underlying genetic mechanism for such an important trait in cucumber are unknown. In this study, a locus controlling pericarp color was found on chromosome 3 of cucumber genome.

RESULTS

In this study, the light green inbred line G35 and the dark green inbred line Q51 were crossed to produce one F₂ population. Consequently, we identified a major locus CsPC1 (Pericarp color 1). Next, we mapped the CsPC1 locus to a 94-kb region chromosome 3 which contains 15 genes. Among these genes, Csa3G912920, which encodes a GATA transcription factor, was expressed at a higher level in the pericarp of the NIL-1334 line (with light-green pericarp) than in that of the NIL-1325 line (with dark-green pericarp). This study provides a new allele for the improvement of cucumber pericarp color.

CONCLUSION

A major QTL that controls pericarp color in cucumber, CsPC1, was identified in a 94-kb region that harbors the strong candidate gene CsGATA1.

Genomic and Transcriptomic Analyses Reveal Pathways and Genes Associated With Brittle Stalk Phenotype in Maize

The mechanical strength of the stalk affects the lodging resistance and digestibility of the stalk in maize. The molecular mechanisms regulating the brittleness of stalks in maize remain undefined. In this study, we constructed the maize brittle stalk mutant (bk5) by crossing the W22:Mu line with the Zheng 58 line. The brittle phenotype of the mutant bk5 existed in all of the plant organs after the five-leaf stage. Compared to wild-type (WT) plants, the sclerenchyma cells of bk5 stalks had a looser cell arrangement and thinner cell wall. Determination of cell wall composition showed that obvious differences in cellulose content, lignin content, starch content, and total soluble sugar were found between bk5 and WT stalks. Furthermore, we identified 226 differentially expressed genes (DEGs), with 164 genes significantly upregulated and 62 genes significantly downregulated in RNA-seq analysis. Some pathways related to cellulose and lignin synthesis, such as endocytosis and glycosylphosphatidylinositol (GPI)-anchored biosynthesis, were identified by the Kyoto Encyclopedia of Gene and Genomes (KEGG) and gene ontology (GO) analysis. In bulked-segregant sequence analysis (BSA-seq), we detected 2,931,692 high-quality Single Nucleotide Polymorphisms (SNPs) and identified five overlapped regions (11.2 Mb) containing 17 candidate genes with missense mutations or premature termination codons using the SNP-index methods. Some genes were involved in the cellulose synthesis-related genes such as ENTH/ANTH/VHS superfamily protein gene (endocytosis-related gene) and the lignin synthesis-related genes such as the cytochrome p450 gene. Some of these candidate genes identified from BSA-seq also existed with differential expression in RNA-seq analysis. These findings increase our understanding of the molecular mechanisms regulating the brittle stalk phenotype in maize.

Conjunctive Analyses of Bulk Segregant Analysis Sequencing and Bulk Segregant RNA Sequencing to Identify Candidate Genes Controlling Spikelet Sterility of Foxtail Millet

Foxtail millet has gradually become a model gramineous C4 crop owing to its short growth period and small genome. Research on the development of its spikelets is not only directly related to the yield and economic value of foxtail millet but also can provide a reference for studying the fertility of other C4 crops. In this study, a hybrid population containing 200 offspring was constructed from the Xinong8852 and An15 parental lines, and two extreme trait populations were constructed from the F2 generation for analysis of the spikelet sterility. The F2 population conformed to a 3:1 Mendelian segregation ratio, and it was thus concluded that this trait is likely controlled by a single recessive gene. Bulk segregant analysis sequencing (BSA-Seq) was used to determine the candidate regions and candidate genes related to the development of foxtail millet spikelets. Additionally, the functional analysis of differentially expressed genes in populations with different traits was conducted by bulk segregant RNA sequencing (BSR-Seq). Finally, conjunctive analysis of BSA-Seq and BSR-Seq results, combined with biological information analysis, revealed six genes on chromosome VII that were ultimately identified as candidate genes controlling foxtail millet spikelet development. This study provides a new reference for research on foxtail millet sterility and lays a solid foundation for the examination of fertility in other gramineous crops.

Identification of Candidate Chromosome Region Related to Melon (Cucumis melo L.) Fruit Surface Groove Trait Through Biparental Genetic Mapping and Genome-Wide Association Study

The melon fruit surface groove (fsg) not only affects peel structure and causes stress-induced fruit cracking but also fits consumers' requirements in different regions. In this study, genetic inheritance analysis of three F₂ populations derived from six parental lines revealed that the fsg trait is controlled by a simple recessive inherited gene. Through bulked segregant analysis sequencing (BSA-seq), the Cmfsg locus was detected in an 8.96 Mb interval on chromosome 11 and then initially mapped to a region of approximately 1.15 Mb. Further fine mapping with a large F₂ population including 1,200 plants narrowed this region to 207 kb containing 11 genes. A genome-wide association study (GWAS) with 187 melon accessions also produced the same chromosome region for the Cmfsg locus. Due to the rare molecular markers and lack of mutations in the coding and promoter regions of the 11 candidate genes in the fine-mapped interval, we conducted in silico BSA to explore the natural melon panel to predict candidate genes for the Cmfsg locus. A 1.07 kb segment upstream of MELO3C019694.2 (annotated as the AGAMOUS MADS-box transcription factor) exhibited a correlation with the grooved and non-grooved accessions among the F₂ individuals, and a natural panel consisted of 17 melon accessions. The expression level of MELO3C019694.2 in the pericarp was higher in grooved lines than in non-grooved lines and was specifically expressed in fruit compared with other tissues (female flower, male flower, root, and leaf). This work provides fundamental information for further research on melon fsg trait formation and molecular markers for melon breeding.

Identifying QTLs for Grain Size in a Colossal Grain Rice (Oryza sativa L.) Line, and Analysis of Additive Effects of QTLs

Grain size is an important component of quality and harvest traits in the field of rice breeding. Although numerous quantitative trait loci (QTLs) of grain size in rice have been reported, the molecular mechanisms of these QTLs remain poorly understood, and further research on QTL observation and candidate gene identification is warranted. In our research, we developed a suite of F2 intercross populations from a cross of 9311 and CG. These primary populations were used to map QTLs conferring grain size, evaluated across three environments, and then subjected to bulked-segregant analysis-seq (BSA-seq). In total, 4, 11, 12 and 14 QTLs for grain length (GL), grain width (GW), 1000-grain weight (TGW), and length/width ratio (LWR), respectively, were detected on the basis of a single-environment analysis. In particular, over 200 splicing-related sites were identified by whole-genome sequencing, including one splicing-site mutation with G>A at the beginning of intron 4 on Os03g0841800 (qGL3.3), producing a smaller open reading frame, without the third and fourth exons. A previous study revealed that the loss-of-function allele caused by this splicing site can negatively regulate rice grain length. Furthermore, qTGW2.1 and qGW2.3 were new QTLs for grain width. We used the near-isogenic lines (NILs) of these GW QTLs to study their genetic effects on individuals and pyramiding, and found that they have additive effects on GW. In summary, these discoveries provide a valuable genetic resource, which will facilitate further study of the genetic polymorphism of new rice varieties in rice breeding.

Fine Mapping and Functional Analysis of the Gene PcTYR, Involved in Control of Cap Color of Pleurotus cornucopiae

Oyster mushrooms have a high biological efficiency and are easy to cultivate, which is why they are produced all over the world. Cap color is an important commercial trait for oyster mushrooms. Little is known about the genetic mechanism of the cap color trait in oyster mushrooms, which limits molecular breeding for the improvement of cap color-type cultivars. In this study, a 0.8-Mb major quantitative trait locus (QTL) region controlling cap color in the oyster mushroom Pleurotus cornucopiae was mapped on chromosome 7 through bulked-segregant analysis sequencing (BSA-seq) and extreme-phenotype genome-wide association studies (XP-GWAS). Candidate genes were further selected by comparative transcriptome analysis, and a tyrosinase gene (PcTYR) was identified as the highest-confidence candidate gene. Overexpression of PcTYR resulted in a significantly darker cap color, while the cap color of RNA interference (RNAi) strains for this gene was significantly lighter than that of the wild-type (WT) strains, suggesting that PcTYR plays an essential role in cap color formation. This is the first report about fine mapping and functional verification of a gene controlling cap color in oyster mushrooms. This will enhance our understanding of the genetic basis for cap color formation in oyster mushrooms and will facilitate molecular breeding for cap color. IMPORTANCE Oyster mushrooms are widely cultivated and consumed over the world for their easy cultivation and high biological efficiency (mushroom fresh weight/substrate dry weight × 100%). Fruiting bodies with dark caps are more and more popular according to consumer preferences, but dark varieties are rarely seen on the market. Little is known about the genetic mechanism of the cap color trait in oyster mushrooms, which limits molecular breeding for the improvement of cap color-type cultivars. A major QTL of cap color in oyster mushroom P. cornucopiae was fine mapped by using bulked-segregant analysis (BSA) and extreme-phenotype genome-wide association study (XP-GWAS) analysis. A candidate gene PcTYR coding tyrosinase was further identified with the help of comparative transcriptome analysis. qPCR analysis and genetic transformation tests proved that PcTYR played an essential role in cap color formation. This study will contribute to revealing the genetic mechanism of cap color formation in mushrooms, thereby facilitating molecular breeding for cap color trait.

CmFSI8/CmOFP13 encoding an OVATE family protein controls fruit shape in melon

Fruit shape is an important quality and yield trait in melon (Cucumis melo). Although some quantitative trait loci for fruit shape have been reported in in this species, the genes responsible and the underlying mechanisms remain poorly understood. Here, we identified and characterized a gene controlling fruit shape from two melon inbred lines, B8 with long-horn fruit and HP22 with flat-round fruit. Genetic analysis suggested that the shape was controlled by a single and incompletely dominant locus, which we designate as CmFSI8/CmOFP13. This gene was finely mapped to a 53.7-kb interval on chromosome 8 based on bulked-segregant analysis sequencing and map-based cloning strategies. CmFSI8/CmOFP13 encodes an OVATE family protein (OFP) and is orthologous to AtOFP1 and SlOFP20. The transcription level of CmFSI8/CmOFP13 in the ovary of HP22 was significantly higher than that in B8, and sequence analysis showed that a 12.5-kb genomic variation with a retrotransposon insertion identified in the promoter was responsible for elevating the expression, and this ultimately caused the differences in fruit shape. Ectopic overexpression of CmFSI8/CmOFP13 in Arabidopsis led to multiple phenotypic changes, including kidney-shaped leaves and shortened siliques. Taken together, our results demonstrate the involvement of an OFP in regulating fruit shape in melon, and our improved understanding of the molecular mechanisms will enable us to better manipulate fruit shape in breeding.

Fine-mapping and candidate gene analysis of a major locus controlling leaf thickness in rice (Oryza sativa L.)

Leaf thickness is an important trait in rice (Oryza sativa L.). It affects both photosynthesis and sink-resource efficiency. However, compared to leaf length and length width, reports seldom focused on leaf thickness due to the complicated measurement and minor difference. To identify the quantitative trait loci (QTL) and explore the genetic mechanism regulating the natural variation of leaf thickness, we crossed a high leaf thickness variety Aixiuzhan (AXZ) to a thin leaf thickness variety Yangdao No.6 (YD 6) and evaluated 585 F₂ individuals. We further use bulked sergeant analysis with whole-genome resequencing (BSA-seq) to identify five genomic regions, including chromosomes 1, 6, 9, 10, and 12. These regions represented significant allele frequency differentiation between thick and thin leaf thickness among the mixed pool offspring. Moreover, we conducted a linkage mapping using 276 individuals derived from the F₂ population. We fine-mapped and confirmed that chromosome 9 contributed the primary explanation of phenotypic variance. We fine-mapped the candidate regions and confirmed that the chromosome 9 region contributed to flag leaf thickness in rice. We observed the virtual cellular slices and found that the bundle sheath cells in YD 6 flag leaf veins are fewer than AXZ. We analyzed the potential regions on chromosome 9 and narrowed the QTL candidate intervals in the 928-kb region. Candidate genes of this major QTL were listed as potentially controlled leaf thickness. These results provide promising evidence that cloning leaf thickness is associated with yield production in rice.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11032-022-01275-y.

Identification, Analysis, and Confirmation of Seed Storability-Related Loci in Dongxiang Wild Rice (Oryza rufipogon Griff.)

Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) has strong seed storability and identifying its elite gene resources may facilitate genetic improvements in rice seed storability. In this study, we developed two backcross inbred lines (BILs) populations, with DXWR as a common donor parent and two rice varieties (F6 and R974) as recipient parents. Bulked segregant analysis via whole genome sequencing (BSA-seq) was used to identify seed storability-related loci in the DXWR and F6 population. Two main genomic regions containing 18,550,000–20,870,000 bp on chromosome 4 and 7,860,000–9,780,000 bp on chromosome 9 were identified as candidate loci of DXWR seed storability; these overlapped partially with seed storability-related quantitative trait loci (QTLs) discovered in previous studies, suggesting that these loci may provide important regions for isolating the responsible genes. In total, 448 annotated genes were predicted within the identified regions, of which 274 and 82 had nonsynonymous and frameshift mutations, respectively. We detected extensive metabolic activities and cellular processes during seed storability and confirmed the effects of the seed storability-related candidate loci using four BILs from DXWR and R974. These results may facilitate the cloning of DXWR seed storability-related genes, thereby elucidating rice seed storability and its improvement potential.

Genomic mapping and identification of candidate genes encoding nulliplex-branch trait in sea-island cotton (Gossypium barbadense L.) by multi-omics analysis

Nulliplex branch is a key architectural trait in sea-island cotton (Gossypium barbadense L.), but its genetic basis is not well understood. Here we investigated the genetic basis of the nulliplex-branch trait in cotton by combining newly created bulked segregant analysis (BSA)-seq data, published RNA-seq data, and published whole-genome resequencing (WGR) data. We delimited the nulliplex-branch locus (qD07-NB) to D07, region 14.8-17.1 Mb, using various BSA methods and markers. We integrated our BSA data with WGR data of sea-island cotton varieties and detected a missense single-nucleotide polymorphism in the candidate gene (Gbar_D07G011870) of qD07-NB. This gene was under positive selection during sea-island cotton breeding in the Xinjiang Uygur Autonomous Region, China. Notably, the nulliplex-branch varieties possessed a better fiber quality than the long-branch varieties, and a set of high-quality molecular markers was identified for molecular breeding of the nulliplex-branch trait in cotton. We combined BSA-seq and RNA-seq data to compare gene expression profiles between two elite sea-island cotton varieties during three developmental stages. We identified eleven relevant candidate genes, five downregulated and six upregulated, in the qD07-NB locus. This research will expand our understanding of the genetic basis of the nulliplex-branch trait and provide guidance for architecture-focused breeding in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01229-w.

Role of MdERF3 and MdERF118 natural variations in apple flesh firmness/crispness retainability and development of QTL-based genomics-assisted prediction

Retention of flesh texture attributes during cold storage is critical for the long-term maintenance of fruit quality. The genetic variations determining flesh firmness and crispness retainability are not well understood. The objectives of this study are to identify gene markers based on quantitative trait loci (QTLs) and to develop genomics-assisted prediction (GAP) models for apple flesh firmness and crispness retainability. Phenotype data of 2664 hybrids derived from three Malus domestica cultivars and a M. asiatica cultivar were collected in 2016 and 2017. The phenotype segregated considerably with high broad-sense heritability of 83.85% and 83.64% for flesh firmness and crispness retainability, respectively. Fifty-six candidate genes were predicted from the 62 QTLs identified using bulked segregant analysis and RNA-seq. The genotype effects of the markers designed on each candidate gene were estimated. The genomics-predicted values were obtained using pyramiding marker genotype effects and overall mean phenotype values. Fivefold cross-validation revealed that the prediction accuracy was 0.5541 and 0.6018 for retainability of flesh firmness and crispness, respectively. An 8-bp deletion in the MdERF3 promoter disrupted MdDOF5.3 binding, reduced MdERF3 expression, relieved the inhibition on MdPGLR3, MdPME2, and MdACO4 expression, and ultimately decreased flesh firmness and crispness retainability. A 3-bp deletion in the MdERF118 promoter decreased its expression by disrupting the binding of MdRAVL1, which increased MdPGLR3 and MdACO4 expression and reduced flesh firmness and crispness retainability. These results provide insights regarding the genetic variation network regulating flesh firmness and crispness retainability, and the GAP models can assist in apple breeding.

QTL Mapping and Identification of Candidate Genes for Heat Tolerance at the Flowering Stage in Rice

High-temperature stress can cause serious abiotic damage that limits the yield and quality of rice. Heat tolerance (HT) during the flowering stage of rice is a key trait that can guarantee a high and stable yield under heat stress. HT is a complex trait that is regulated by multiple quantitative trait loci (QTLs); however, few underlying genes have been fine mapped and cloned. In this study, the F_2:3 population derived from a cross between Huanghuazhan (HHZ), a heat-tolerant cultivar, and 9311, a heat-sensitive variety, was used to map HT QTLs during the flowering stage in rice. A new major QTL, qHTT8, controlling HT was identified on chromosome 8 using the bulked-segregant analysis (BSA)-seq method. The QTL qHTT8 was mapped into the 3,555,000–4,520,000 bp, which had a size of 0.965 Mb. The candidate region of qHTT8 on chromosome 8 contained 65 predicted genes, and 10 putative predicted genes were found to be associated with abiotic stress tolerance. Furthermore, qRT-PCR was performed to analyze the differential expression of these 10 genes between HHZ and 9311 under high temperature conditions. LOC_Os08g07010 and LOC_Os08g07440 were highly induced in HHZ compared with 9311 under heat stress. Orthologous genes of LOC_Os08g07010 and LOC_Os08g07440 in plants played a role in abiotic stress, suggesting that they may be the candidate genes of qHTT8. Generally, the results of this study will prove useful for future efforts to clone qHTT8 and breed heat-tolerant varieties of rice using marker-assisted selection.

Cytological Observations and Bulked-Segregant Analysis Coupled Global Genome Sequencing Reveal Two Genes Associated with Pollen Fertility in Tetraploid Rice

Neo-tetraploid rice with high fertility is a useful germplasm for polyploid rice breeding, which was developed from the crossing of different autotetraploid rice lines. However, little information is available on the molecular mechanism underlying the fertility of neo-tetraploid rice. Here, two contrasting populations of tetraploid rice, including one with high fertility (hereafter referred to as JG) and another with low fertility (hereafter referred to as JD), were generated by crossing Huaduo 3 (H3), a high fertility neo-tetraploid rice that was developed by crossing Jackson-4x with 96025-4x, and Huajingxian74-4x (T452), a low fertility autotetraploid rice parent. Cytological, global genome sequencing-based bulked-segregant (BSA-seq) and CRISPR/Cas9 technology were employed to study the genes associated with pollen fertility in neo-tetraploid rice. The embryo sacs of JG and JD lines were normal; however, pollen fertility was low in JD, which led to scarce fertilization and low seed setting. Cytological observations displayed low pollen fertility (25.1%) and approximately 31.3 and 27.2% chromosome lagging at metaphase I and II, and 28.8 and 24.8% chromosome straggling at anaphase I and II in JD, respectively. BSA-seq of F_2–3 generations and RNA-seq of F₄ generation detected a common fragment, i.e., 18,915,234–19,500,000, at chromosome 7, which was comprised of 78 genes associated with fertility. Among 78 genes, 9 genes had been known to be involved in meiosis and pollen development. Two mutants ny1 (LOC_Os07g32406) and ny2 (LOC_Os07g32040) were generated by CRISPR/Cas9 knockout in neo-tetraploid rice, and which exhibited low pollen fertility and abnormal chromosome behavior. Our study revealed that two unknown genes, LOC_Os07g32406 (NY1) and LOC_Os07g32040 (NY2) play an important role in pollen development of neo-tetraploid rice and provides a new perspective about the genetic mechanisms of fertility in polyploid rice.

Rapid Mining of Candidate Genes for Verticillium Wilt Resistance in Cotton Based on BSA-Seq Analysis

Cotton is a globally important cash crop. Verticillium wilt (VW) is commonly known as "cancer" of cotton and causes serious loss of yield and fiber quality in cotton production around the world. Here, we performed a BSA-seq analysis using an F_2:3 segregation population to identify the candidate loci involved in VW resistance. Two QTLs (qvw-D05-1 and qvw-D05-2) related to VW resistance in cotton were identified using two resistant/susceptible bulks from the F₂ segregation population constructed by crossing the resistant cultivar ZZM2 with the susceptible cultivar J11. A total of 30stop-lost SNPs and 42 stop-gained SNPs, which included 17 genes, were screened in the qvw-D05-2 region by SnpEff analysis. Further analysis of the transcriptome data and qRT-PCR revealed that the expression level of Ghir_D05G037630 (designated as GhDRP) varied significantly at certain time points after infection with V. dahliae. The virus-induced gene silencing of GhDRP resulted in higher susceptibility of the plants to V. dahliae than the control, suggesting that GhDRP is involved in the resistance to V. dahlia infection. This study provides a method for rapid mining of quantitative trait loci and screening of candidate genes, as well as enriches the genomic information and gene resources for the molecular breeding of disease resistance in cotton.

Intricate genetic variation networks control the adventitious root growth angle in apple

Background

The root growth angle (RGA) typically determines plant rooting depth, which is significant for plant anchorage and abiotic stress tolerance. Several quantitative trait loci (QTLs) for RGA have been identified in crops. However, the underlying mechanisms of the RGA remain poorly understood, especially in apple rootstocks. The objective of this study was to identify QTLs, validate genetic variation networks, and develop molecular markers for the RGA in apple rootstock.

Results

Bulked segregant analysis by sequencing (BSA-seq) identified 25 QTLs for RGA using 1955 hybrids of the apple rootstock cultivars ‘Baleng Crab’ (Malus robusta Rehd., large RGA) and ‘M9’ (M. pumila Mill., small RGA). With RNA sequencing (RNA-seq) and parental resequencing, six major functional genes were identified and constituted two genetic variation networks for the RGA. Two single nucleotide polymorphisms (SNPs) of the MdLAZY1 promoter damaged the binding sites of MdDREB2A and MdHSFB3, while one SNP of MdDREB2A and MdIAA1 affected the interactions of MdDREB2A/MdHSFB3 and MdIAA1/MdLAZY1, respectively. A SNP within the MdNPR5 promoter damaged the interaction between MdNPR5 and MdLBD41, while one SNP of MdLBD41 interrupted the MdLBD41/MdbHLH48 interaction that affected the binding ability of MdLBD41 on the MdNPR5 promoter. Twenty six SNP markers were designed on candidate genes in each QTL interval, and the marker effects varied from 0.22°-26.11°.

Conclusions

Six diagnostic markers, SNP592, G122, b13, Z312, S1272, and S1288, were used to identify two intricate genetic variation networks that control the RGA and may provide new insights into the accuracy of the molecular markers. The QTLs and SNP markers can potentially be used to select deep-rooted apple rootstocks.

Characterization and Fine Mapping of a Yellow-Virescent Gene Regulating Chlorophyll Biosynthesis and Early Stage Chloroplast Development in Brassica napus

Chlorophyll biosynthesis and chloroplast development are crucial to photosynthesis and plant growth, but their regulatory mechanism remains elusive in many crop species. We isolated a Brassica napus yellow-virescent leaf (yvl) mutant, which exhibited yellow-younger-leaf and virescent-older-leaf with decreased chlorophyll accumulation and delayed chloroplast development. We mapped yvl locus to a 70-kb interval between molecular markers yvl-O10 and InDel-O6 on chromosome A03 in BC₂F₂ population using whole genome re-sequencing and bulked segregant analysis. The mutant had a ‘C’ to ‘T’ substitution in the coding sequence of BnaA03.CHLH, which encodes putative H subunit of Mg-protoporphyrin IX chelatase (CHLH). The mutation resulted in an imperfect protein structure and reduced activity of CHLH. It also hampered the plastid encoded RNA polymerase which transcribes regulatory genes of photosystem II and I. Consequently, the chlorophyll a/b and carotenoid contents were reduced and the chloroplast ultrastructure was degraded in yvl mutant. These results explain that a single nucleotide mutation in BnaA03.CHLH impairs PEP activity to disrupt chloroplast development and chlorophyll biosynthesis in B. napus.

Comparative Transcriptome Analysis Provides Insights Into Yellow Rind Formation and Preliminary Mapping of the Clyr (Yellow Rind) Gene in Watermelon

As an important appearance trait, the rind color of watermelon fruit affects the commodity value and further determines consumption choices. In this study, a comparative transcriptome analysis was conducted to elucidate the genes and pathways involved in the formation of yellow rind fruit in watermelon using a yellow rind inbred line WT4 and a green rind inbred line WM102. A total of 2,362 differentially expressed genes (DEGs) between WT4 and WM102 at three different stages (0, 7, and 14 DAP) were identified and 9,770 DEGs were obtained by comparing the expression level at 7 DAP and 14 DAP with the former stages of WT4. The function enrichment of DEGs revealed a number of pathways and terms in biological processes, cellular components, and molecular functions that were related to plant pigment metabolism, suggesting that there may be a group of common core genes regulating rind color formation. In addition, next-generation sequencing aided bulked-segregant analysis (BSA-seq) of the yellow rind pool and green rind pool selected from an F₂ population revealed that the yellow rind gene (Clyr) was mapped on the top end of chromosome 4. Based on the BSA-seq analysis result, Clyr was further confined to a region of 91.42 kb by linkage analysis using 1,106 F₂ plants. These results will aid in identifying the key genes and pathways associated with yellow rind formation and elucidating the molecular mechanism of rind color formation in watermelon.

Genomic and transcriptomic analysis reveal molecular basis of salinity tolerance in a novel strong salt-tolerant rice landrace Changmaogu

BACKGROUND

Salt stress is an important factor that limits rice yield. We identified a novel, strongly salt tolerant rice landrace called Changmaogu (CMG) collected from a coastal beach of Zhanjiang, Guangdong Province, China. The salt tolerance of CMG was much better than that of the international recognized salt tolerant rice cultivar Pokkali in the germination and seedling stages.

RESULTS

To understand the molecular basis of salt tolerance in CMG, we performed BSA-seq for two extreme bulks derived from the cross between CMG and a cultivar sensitive to salt, Zhefu802. Transcriptomic sequencing was conducted for CMG at the germination and young seedling stages. Six candidate regions for salt tolerance were mapped on Chromosome 1 by BSA-seq using the extreme populations. Based on the polymorphisms identified between both parents, we detected 32 genes containing nonsynonymous coding single nucleotide polymorphisms (SNPs) and frameshift mutations in the open reading frame (ORF) regions. With transcriptomic sequencing, we detected a large number of differentially expressed genes (DEGs) at the germination and seedling stages under salt stress. KEGG analysis indicated two of 69 DEGs shared at the germination and seedling stages were significantly enriched in the pathway of carotenoid biosynthesis. Of the 169 overlapping DEGs among three sample points at the seedling stage, 13 and six DEGs were clustered into the pathways of ABA signal transduction and carotenoid biosynthesis, respectively. Of the 32 genes carrying sequence variation, only OsPP2C8 (Os01g0656200) was differentially expressed in the young seedling stage under salt stress and also showed sequence polymorphism in the ORFs between CMG and Zhefu802.

CONCLUSION

OsPP2C8 was identified as the target candidate gene for salinity tolerance in the seedling stage. This provides an important genetic resource for the breeding of novel salt tolerant rice cultivars.

Whole-genome resequencing from bulked-segregant analysis reveals gene set based association analyses for the Vibrio anguillarum resistance of turbot (Scophthalmus maximus)

Many achievements have been made to develop quantitative trait loci (QTLs) and gene-associated single nucleotide polymorphisms (SNPs) to facilitate practical marker-assisted selection (MAS) in aquatic animals. However, the systematic studies of SNPs associated with extreme threshold traits were poor in populations lacking of parental genomic information. Coupling next generation sequencing with bulked segregant analysis (BSA) should allow identification of numerous associated SNPs with extreme phenotypes. In the present study, using combination of SNP frequency difference and Euclidean distance, we conducted linkage analysis of SNPs located in genes involved in immune responses, and identified markers associated with Vibrio anguillarum resistance in turbot (Scophthalmus maximus). A total of 221 SNPs was found as candidate SNPs between resistant and susceptible individuals. Among these SNPs, 35 loci located in immune related genes were genotyped in verification population and 7 of them showed significant association with V. anguillarum resistance in both alleles and genotypes (P < 0.05). Among these 7 genes, PIK3CA-like, CYLD, VCAM1, RhoB and RhoGEF are involved in PI3K/Akt/mTOR pathway and NF-κB pathway, which influence the efficiency of bacteria entering the host and inflammation. SNP-SNP interaction analysis was performed by generalized multifactor dimensionality reduction (GMDR). The combination of SNP loci in RhoB, PIK3CA-like and ADCY3 showed a significant effect on V. anguillarum resistance with the verification rate in the sequencing population up to 70.8%. Taken all, our findings demonstrated the feasibility of BSA-seq approach in identifying genes responsible for the extreme phenotypes and will aid in performing MAS in turbot.

A Single Nucleotide Deletion in an ABC Transporter Gene Leads to a Dwarf Phenotype in Watermelon

Dwarf habit is one of the most important traits in crop plant architecture, as it can increase plant density and improved land utilization, especially for protected cultivation, as well as increasing lodging resistance and economic yield. At least four dwarf genes have been identified in watermelon, but none of them has been cloned. In the current study, the Cldw-1 gene was primary-mapped onto watermelon chromosome 9 by next-generation sequencing-aided bulked-segregant analysis (BSA-seq) of F₂ plants derived from a cross between a normal-height line, WT4, and a dwarf line, WM102, in watermelon. The candidate region identified by BSA-seq was subsequently validated and confirmed by linkage analysis using 30 simple sequence repeat (SSR) markers in an F₂ population of 124 plants. The Cldw-1 gene was further fine-mapped by chromosome walking in a large F₂ population of 1,053 plants and was delimited into a candidate region of 107.00 kb. Six genes were predicted to be in the candidate region, and only one gene, Cla010337, was identified to have two single nucleotide polymorphisms (SNPs) and a single nucleotide deletion in the exons in the dwarf line, WM102. A derived cleaved amplified polymorphic sequence (dCAPS) marker was developed from the single nucleotide deletion, co-segregated with the dwarf trait in both the F₂ population and a germplasm collection of 165 accessions. Cla010337 encoded an ATP-binding cassette transporter (ABC transporter) protein, and the expression levels of Cla010337 were significantly reduced in all the tissues tested in the dwarf line, WM102. The results of this study will be useful in achieving a better understanding of the molecular mechanism of the dwarf plant trait in watermelon and for the development of marker-assisted selection (MAS) for new dwarf cultivars.

Mapping Cucumber Vein Yellowing Virus Resistance in Cucumber (Cucumis sativus L.) by Using BSA-seq Analysis

Cucumber vein yellowing virus (CVYV) causes severe yield losses in cucurbit crops across Mediterranean countries. The control of this virus is based on cultural practices to prevent the presence of its vector (Bemisia tabaci) and breeding for natural resistance, which requires the identification of the loci involved and the development of molecular markers for linkage analysis. In this work, we mapped a monogenic locus for resistance to CVYV in cucumber by using a Bulked Segregant Analysis (BSA) strategy coupled with whole-genome resequencing. We phenotyped 135 F₃ families from a segregating population between a susceptible pickling cucumber and a resistant Long Dutch type cucumber for CVYV resistance. Phenotypic analysis determined the monogenic and incomplete dominance inheritance of the resistance. We named the locus CsCvy-1. For mapping this locus, 15 resistant and 15 susceptible homozygous F₂ individuals were selected for whole genome resequencing. By using a customized bioinformatics pipeline, we identified a unique region in chromosome 5 associated to resistance to CVYV, explaining more than 80% of the variability. The resequencing data provided us with additional SNP markers to decrease the interval of CsCvy-1 to 625 kb, containing 24 annotated genes. Markers flanking CsCvy-1 in a 5.3 cM interval were developed for marker-assisted selection (MAS) in breeding programs and will be useful for the identification of the target gene in future studies.

Dissection of quantitative traits by bulk segregant mapping in a protoploid yeast species

Since more than a decade ago, Saccharomyces cerevisiae has been used as a model to dissect complex traits, revealing the genetic basis of a large number of traits in fine detail. However, to have a more global view of the genetic architecture of traits across species, the examination of the molecular basis of phenotypes within non-conventional species would undoubtedly be valuable. In this respect, the Saccharomycotina yeasts represent ideal and potential non-model organisms. Here we sought to assess the feasibility of genetic mapping by bulk segregant analysis in the protoploid Lachancea kluyveri (formerly S. kluyveri) yeast species, a distantly related species to S. cerevisiae For this purpose, we designed a fluorescent mating-type marker, compatible with any mating-competent strains representative of this species, to rapidly create a large population of haploid segregants (>10(5) cells). Quantitative trait loci can be mapped by selecting and sequencing an enriched pool of progeny with extreme phenotypic values. As a test bed, we applied this strategy and mapped the causal loci underlying halotolerance phenotypes in L. kluyveri Overall, this study demonstrates that bulk segregant mapping is a powerful way for investigating the genetic basis of natural variations in non-model yeast organisms and more precisely in L. kluyveri.