Identification of two major QTLs for pod shell thickness in peanut (Arachis hypogaea L.) using BSA-seq analysis.

BACKGROUND: Pod shell thickness (PST) is an important agronomic trait of peanut because it affects the ability of shells to resist pest infestations and pathogen attacks, while also influencing the peanut shelling process. However, very few studies have explored the genetic basis of PST. RESULTS: An F2 segregating population derived from a cross between the thick-shelled cultivar Yueyou 18 (YY18) and the thin-shelled cultivar Weihua 8 (WH8) was used to identify the quantitative trait loci (QTLs) for PST. On the basis of a bulked segregant analysis sequencing (BSA-seq), four QTLs were preliminarily mapped to chromosomes 3, 8, 13, and 18. Using the genome resequencing data of YY18 and WH8, 22 kompetitive allele-specific PCR (KASP) markers were designed for the genotyping of the F2 population. Two major QTLs (qPSTA08 and qPSTA18) were identified and finely mapped, with qPSTA08 detected on chromosome 8 (0.69-Mb physical genomic region) and qPSTA18 detected on chromosome 18 (0.15-Mb physical genomic region). Moreover, qPSTA08 and qPSTA18 explained 31.1-32.3% and 16.7-16.8% of the phenotypic variation, respectively. Fifteen genes were detected in the two candidate regions, including three genes with nonsynonymous mutations in the exon region. Two molecular markers (Tif2_A08_31713024 and Tif2_A18_7198124) that were developed for the two major QTL regions effectively distinguished between thick-shelled and thin-shelled materials. Subsequently, the two markers were validated in four F2:3 lines selected. CONCLUSIONS: The QTLs identified and molecular markers developed in this study may lay the foundation for breeding cultivars with a shell thickness suitable for mechanized peanut shelling.

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 candidate genes controlling cold tolerance at the early seedling stage from Dongxiang wild rice by QTL mapping, BSA-Seq and RNA-Seq.

BACKGROUND: The cold tolerance of rice is closely related to its production and geographic distribution. The identification of cold tolerance-related genes is of important significance for developing cold-tolerant rice. Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) is well-adapted to the cold climate of northernmost-latitude habitats ever found in the world, and is one of the most valuable rice germplasms for cold tolerance improvement. RESULTS: Transcriptome analysis revealed genes differentially expressed between Xieqingzao B (XB; a cold sensitive variety) and 19H19 (derived from an interspecific cross between DXWR and XB) in the room temperature (RT), low temperature (LT), and recovery treatments. The results demonstrated that chloroplast genes might be involved in the regulation of cold tolerance in rice. A high-resolution SNP genetic map was constructed using 120 BC5F2 lines derived from a cross between 19H19 and XB based on the genotyping-by-sequencing (GBS) technique. Two quantitative trait loci (QTLs) for cold tolerance at the early seedling stage (CTS), qCTS12 and qCTS8, were detected. Moreover, a total of 112 candidate genes associated with cold tolerance were identified based on bulked segregant analysis sequencing (BSA-seq). These candidate genes were divided into eight functional categories, and the expression trend of candidate genes related to 'oxidation-reduction process' and 'response to stress' differed between XB and 19H19 in the RT, LT and recovery treatments. Among these candidate genes, the expression level of LOC_Os12g18729 in 19H19 (related to 'response to stress') decreased in the LT treatment but restored and enhanced during the recovery treatment whereas the expression level of LOC_Os12g18729 in XB declined during recovery treatment. Additionally, XB contained a 42-bp deletion in the third exon of LOC_Os12g18729, and the genotype of BC5F2 individuals with a survival percentage (SP) lower than 15% was consistent with that of XB. Weighted gene coexpression network analysis (WGCNA) and modular regulatory network learning with per gene information (MERLIN) algorithm revealed a gene interaction/coexpression network regulating cold tolerance in rice. In the network, differentially expressed genes (DEGs) related to 'oxidation-reduction process', 'response to stress' and 'protein phosphorylation' interacted with LOC_Os12g18729. Moreover, the knockout mutant of LOC_Os12g18729 decreased cold tolerance in early rice seedling stage signifcantly compared with that of wild type. CONCLUSIONS: In general, study of the genetic basis of cold tolerance of rice is important for the development of cold-tolerant rice varieties. In the present study, QTL mapping, BSA-seq and RNA-seq were integrated to identify two CTS QTLs qCTS8 and qCTS12. Furthermore, qRT-PCR, genotype sequencing and knockout analysis indicated that LOC_Os12g18729 could be the candidate gene of qCTS12. These results are expected to further exploration of the genetic mechanism of CTS in rice and improve cold tolerance of cultivated rice by introducing the cold tolerant genes from DXWR through marker-assisted selection.

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.

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.

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.

Combined genomics to discover genes associated with tolerance to soil carbonate.

Carbonate-rich soils limit plant performance and crop production. Previously, local adaptation to carbonated soils was detected in wild Arabidopsis thaliana accessions, allowing the selection of two demes with contrasting phenotypes: A1 (carbonate tolerant, c+) and T6 (carbonate sensitive, c-). Here, A1(c+) and T6(c - ) seedlings were grown hydroponically under control (pH 5.9) and bicarbonate conditions (10 mM NaHCO3 , pH 8.3) to obtain ionomic profiles and conduct transcriptomic analysis. In parallel, A1(c+) and T6(c - ) parental lines and their progeny were cultivated on carbonated soil to evaluate fitness and segregation patterns. To understand the genetic architecture beyond the contrasted phenotypes, a bulk segregant analysis sequencing (BSA-Seq) was performed. Transcriptomics revealed 208 root and 2503 leaf differentially expressed genes in A1(c+) versus T6(c - ) comparison under bicarbonate stress, mainly involved in iron, nitrogen and carbon metabolism, hormones and glycosylates biosynthesis. Based on A1(c+) and T6(c - ) genome contrasts and BSA-Seq analysis, 69 genes were associated with carbonate tolerance. Comparative analysis of genomics and transcriptomics discovered a final set of 18 genes involved in bicarbonate stress responses that may have relevant roles in soil carbonate tolerance.

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.

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.

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.

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.

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.

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 F2 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.

Evaluation of nine statistics to identify QTLs in bulk segregant analysis using next generation sequencing approaches.

BACKGROUND: Bulk segregant analysis (BSA) combined with next generation sequencing is a powerful tool to identify quantitative trait loci (QTL). The impact of the size of the study population and the percentage of extreme genotypes analysed have already been assessed. But a good comparison of statistical approaches designed to identify QTL regions using next generation sequencing (NGS) technologies for BSA is still lacking. RESULTS: We developed an R code to simulate QTLs in bulks of F2 contrasted lines. We simulated a range of recombination rates based on estimations using different crop species. The simulations were used to benchmark the ability of statistical methods identify the exact location of true QTLs. A single QTL led to a shift in allele frequency across a large fraction of the chromosome for plant species with low recombination rate. The smoothed version of all statistics performed best notably the smoothed Euclidean distance-based statistics was always found to be more accurate in identifying the location of QTLs. We propose a simulation approach to build confidence interval statistics for the detection of QTLs. CONCLUSION: We highlight the statistical methods best suited for BSA studies using NGS technologies in crops even when recombination rate is low. We also provide simulation codes to build confidence intervals and to assess the impact of recombination for application to other studies. This computational study will help select NGS-based BSA statistics that are useful to the broad scientific community.

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.

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.

Conuping BSA-Seq and RNA-Seq Reveal the Molecular Pathway and Genes Associated with the Plant Height of Foxtail Millet (Setaria italica).

Foxtail millet (Setaria italica) plays an important role in C4 crop research and agricultural development in arid areas due to its short growth period, drought tolerance, and barren tolerance. Exploration of the dwarfing mechanism and the dwarf genes of foxtail millet can provide a reference for dwarf breeding and dwarf research of other C4 crops. In this study, genetic analysis was performed using phenotypic data, candidate genes were screened by bulk segregant analysis sequencing (BSA-Seq); differentially expressed genes and metabolic pathways in different strains of high samples were analyzed by RNA sequencing (RNA-Seq). The association analysis of BSA-Seq and RNA-Seq further narrowed the candidate range. As a result, a total of three quantitative trait loci (QTLs) and nine candidate genes related to plant height were obtained on chromosomes I and IX. Based on the functional prediction of the candidate genes, we propose a hypothetical mechanism for the formation of millet dwarfing, in which, metabolism and MAPK signaling play important roles in the formation of foxtail millet plant height.

Conjunctive Analyses of BSA-Seq and BSR-Seq Unveil the Msß-GAL and MsJMT as Key Candidate Genes for Cytoplasmic Male Sterility in Alfalfa (Medicago sativa L.).

Knowing the molecular mechanism of male sterility in alfalfa is important to utilize the heterosis more effectively. However, the molecular mechanisms of male sterility in alfalfa are still unclear. In this study, the bulked segregant analysis (BSA) and bulked segregant RNA-seq (BSR) were performed with F2 separation progeny to study the molecular mechanism of male sterility in alfalfa. The BSA-seq analysis was located in a candidate region on chromosome 5 containing 626 candidate genes which were associated with male sterility in alfalfa, while the BSR-seq analysis filtered seven candidate DEGs related to male sterility, and these candidate genes including EF-Tu, ß-GAL, CESA, PHGDH, and JMT. The conjunctive analyses of BSR and BSA methods revealed that the genes of Msß-GAL and MsJMT are the common detected candidate genes involved in male sterility in alfalfa. Our research provides a theory basis for further study of the molecular mechanism of male sterility in alfalfa and significant information for the genetic breeding of Medicago sativa.

Combined BSA-Seq Based Mapping and RNA-Seq Profiling Reveal Candidate Genes Associated with Plant Architecture in Brassica napus.

Plant architecture involves important agronomic traits affecting crop yield, resistance to lodging, and fitness for mechanical harvesting in Brassica napus. Breeding high-yield varieties with plant architecture suitable for mechanical harvesting is the main goal of rapeseed breeders. Here, we report an accession of B. napus (4942C-5), which has a dwarf and compact plant architecture in contrast to cultivated varieties. A BC8 population was constructed by crossing a normal plant architecture line, 8008, with the recurrent parent 4942C-5. To investigate the molecular mechanisms underlying plant architecture, we performed phytohormone profiling, bulk segregant analysis sequencing (BSA-Seq), and RNA sequencing (RNA-Seq) in BC8 plants with contrasting plant architecture. Genetic analysis indicated the plant architecture traits of 4942C-5 were recessive traits controlled by multiple genes. The content of auxin (IAA), gibberellin (GA), and abscisic acid (ABA) differed significantly between plants with contrasting plant architecture in the BC8 population. Based on BSA-Seq analysis, we identified five candidate intervals on chromosome A01, namely those of 0 to 6.33 Mb, 6.45 to 6.48 Mb, 6.51 to 6.53 Mb, 6.77 to 6.79 Mb, and 7 to 7.01 Mb regions. The RNA-Seq analysis revealed a total of 4378 differentially expressed genes (DEGs), of which 2801 were up-regulated and 1577 were down-regulated. There, further analysis showed that genes involved in plant hormone biosynthesis and signal transduction, cell structure, and the phenylpropanoid pathway might play a pivotal role in the morphogenesis of plant architecture. Association analysis of BSA-Seq and RNA-Seq suggested that seven DEGs involved in plant hormone signal transduction and a WUSCHEL-related homeobox (WOX) gene (BnaA01g01910D) might be candidate genes responsible for the dwarf and compact phenotype in 4942C-5. These findings provide a foundation for elucidating the mechanisms underlying rapeseed plant architecture and should contribute to breed new varieties suitable for mechanization.