Identification of candidate genes associated with double flowers via integrating BSA-seq and RNA-seq in Brassica napus.

As a Brassica crop, Brassica napus typically has single flowers that contain four petals. The double-flower phenotype of rapeseed has been a desirable trait in China because of its potential commercial value in ornamental tourism. However, few double-flowered germplasms have been documented in B. napus, and knowledge of the underlying genes is limited. Here, B. napus D376 was characterized as a double-flowered strain that presented an average of 10.92 ± 1.40 petals and other normal floral organs. F1, F2 and BC1 populations were constructed by crossing D376 with a single-flowered line reciprocally. Genetic analysis revealed that the double-flower trait was a recessive trait controlled by multiple genes. To identify the key genes controlling the double-flower trait, bulk segregant analysis sequencing (BSA-seq) and RNA-seq analyses were conducted on F2 individual bulks with opposite extreme phenotypes. Through BSA-seq, one candidate interval was mapped at the region of chromosome C05: 14.56-16.17 Mb. GO and KEGG enrichment analyses revealed that the DEGs were significantly enriched in carbohydrate metabolic processes, notably starch and sucrose metabolism. Interestingly, five and thirty-six DEGs associated with floral development were significantly up- and down-regulated, respectively, in the double-flowered plants. A combined analysis of BSA-seq and RNA-seq data revealed that five genes were candidates associated with the double flower trait, and BnaC05.ERS2 was the most promising gene. These findings provide novel insights into the breeding of double-flowered varieties and lay a theoretical foundation for unveiling the molecular mechanisms of floral development in B. napus.

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.

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.

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.

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.

A novel Pik allele confers extended resistance to rice blast.

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

A novel SNP within the Rsa10025320 gene is highly associated with hollowness in red-skinned radish fleshy roots.

Hollowness is a physiological disorder that frequently occurs during the growth and postharvest storage phases of fleshy radish roots, significantly diminishing quality, yield, and marketability. However, the molecular mechanism for hollowness remains elusive. To identify the QTLs and potential candidate genes for hollowness tolerance in radish, F2 and BC1 populations were constructed from hollowness-tolerant radish (C16) and hollowness-sensitive radish (C17) in the present study. Genetic analysis indicated that hollowness tolerance may be governed by two independent recessive genes. By employing bulked segregant analysis sequencing (BSA-seq), two significant candidate genomic intervals were pinpointed on chromosomes R04 (960 kb, 6.48-7.44 Mb) and R05 (600 kb, 31.44-32.04 Mb), which together harbor 107 annotated genes. Transcriptomic sequencing revealed that the downregulated differentially expressed genes (DEGs) were significantly enriched in biological processes related to cell death and the response to water stress, whereas the upregulated DEGs were significantly associated with the chitin catabolic process and the cell wall macromolecule metabolic process. A total of 46 intersecting genes were identified among these DEGs within the genomic intervals of interest. One gene with high expression (Rsa10025345) and two with low expression (Rsa10025320 and Rsa10018106) were detected in the tolerant variety C16. Furthermore, a SNP within Rsa10025320 resulting in an amino acid change (A188E) was characterized through sequence variation observed in both BSA-seq and RNA-seq data and further developed as a derived cleaved amplified polymorphic sequence (dCAPS) marker. Our study reveals potential target genes for tolerance to hollowness and paves the way for marker-assisted breeding of hollowness tolerance in red-skinned radishes.

Development and application of Key Allele-Specific PCR (KASP) molecular markers for assessing apple fruit crispness.

Crispness stands as a pivotal criterion in assessing apple texture, widely cherished by consumers. Yet, owing to its multifaceted nature, crispness remains a formidable challenge in artificial enhancement efforts. To expedite the early and precise evaluation of apple crispness, this study centered on a hybrid population derived from 'Fuji' and 'Pink Lady' cultivars, showcasing segregating crispness traits. We conducted measurements of flesh water content, cellular anatomical morphology, and employed a texture analyzer to assess mechanical properties of the offspring flesh. Integrating these three dimensions, we conducted a comprehensive analysis of quantitative characteristics of apple crispness, juxtaposed with sensory evaluation. Utilizing BSA-seq technology, we scrutinized extreme phenotypic individuals, revealing QTL loci intricately linked to the aforementioned dimensions, and subsequently developed Key Allele-Specific PCR (KASP) markers. These markers underwent validation in hybrid populations of 'Hanfu' x 'Pink Lady' and 'Hanfu' x 'Honey Crisp'. Our findings underscored significant correlations between mechanical properties, water content, and cell size with crispness. Higher mechanical properties and water content, alongside smaller cell size, correlated with firmer flesh texture; moderate mechanical properties, and elevated water content and cell size, with crisper texture; whereas lower mechanical properties, water content, and cell size implied softer flesh.The study yielded KASP markers effectively reflecting flesh mechanical properties (SNP_24399345), water content (SNP_8667563), and cell size (SNP_15566229). Comprehensive analysis of these markers identified CC-CC-TT as an effective identifier of soft flesh individuals; while GG-TC-TT and GG-CC-TT combinations better represented individuals with harder flesh. The Crunchy subclass could be discerned by combinations of GG-TC-TC, GG-TC-CC, GG-TT-TC, and GG-TT-CC. These findings furnish effective molecular markers for the genetic enhancement of apple crispness, bearing significant implications for the cultivation of novel apple varieties. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01509-1.

Combined BSA-Seq and RNA-Seq Analysis to Identify Candidate Genes Associated with Aluminum Toxicity in Rapeseed (Brassica napus L.).

Exchangeable aluminum (Al) ions released from acidic soils with pH < 5.5 inhibit root elongation of crops, ultimately leading to yield reduced. It is necessary to identify the quantitative trait locus (QTLs) and candidate genes that confer toxicity resistance to understand the mechanism and improve tolerance of rapeseed. In this study, an F2 segregating population was derived from a cross between Al-tolerance inbred line FDH188 (R178) and -sensitive inbred line FDH152 (S169), and the F2:3 were used as materials to map QTLs associated with the relative elongation of taproot (RET) under Al toxicity stress. Based on bulked segregant analysis sequencing (BSA-seq), three QTLs (qAT-A07-1, qAT-A07-2, and qAT-A09-1) were detected as significantly associated with RET, and 656 candidate genes were screened. By combined BSA and RNA-seq analysis, 55 candidate genes showed differentially expressed, including genes encoding ABC transporter G (ABCG), zinc finger protein, NAC, ethylene-responsive transcription factor (ERF), etc. These genes were probably positive factors in coping with Al toxicity stress in rapeseed. This study provides new insight into exploring the QTLs and candidate genes' response to Al toxicity stress by combined BSA-seq and RNA-seq and is helpful to further research on the mechanism of Al resistance in rapeseed.

Identification of Rapeseed (Brassica napus L.) Plant Height-Associated QTL Using BSA-seq and RNA-seq.

Plant height (PH) is a critical agronomic trait in Brassica napus, significantly impacting yield. Consequently, identifying genes associated with plant height is a pivotal objective in oilseed rape breeding. This study employed a combination of bulk segregant analysis sequencing (BSA-seq) and RNA sequencing (RNA-seq) for analysis. A novel quantitative trait locus (QTL), qPH_C02, was identified between 63,989,634 and 64,945,122 bp on chromosome C02, from which eight candidate genes were screened. The Gene Ontology (GO) analysis revealed enrichment in peroxisomes, while the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated enrichment in the oxidative phosphorylation (OP) pathway. It is hypothesized that the observed differences in plant height and silique length may be attributed to the regulation of peroxidase activity in the OP pathway, which in turn alters plant energy metabolism and controls nutrient uptake. Subsequently, we will further test this hypothesis. The results of this study will contribute to our understanding of the genetic basis for differences in plant height and provide a foundation for the selection and breeding of Brassica napus varieties with desired plant shapes.

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.

Function Analysis of a Maize Endo-1,4-ß-xylanase Gene ZmHSL in Response to High-Temperature Stress.

Rising temperature is a major threat to the normal growth and development of maize, resulting in low yield production and quality. The mechanism of maize in response to heat stress remains uncertain. In this study, a maize mutant Zmhsl-1 (heat sensitive leaves) with wilting and curling leaves under high temperatures was identified from maize Zheng 58 (Z58) mutant lines generated by ethyl methanesulfonate (EMS) mutagenesis. The Zmhsl-1 plants were more sensitive to increased temperature than Z58 in the field during growth season. The Zmhsl-1 plants had lower plant height, lower yield, and lower content of photosynthetic pigments. A bulked segregant analysis coupled with whole-genome sequencing (BSA-seq) enabled the identification of the corresponding gene, named ZmHSL, which encodes an endo-ß-1,4-xylanase from the GH10 family. The loss-of-function of ZmHSL resulted in reduced lignin content in Zmhsl-1 plants, leading to defects in water transport and more severe leaf wilting with the increase in temperature. RNA-seq analysis revealed that the differentially expressed genes identified between Z58 and Zmhsl-1 plants are mainly related to heat stress-responsive genes and unfolded protein response genes. All these data indicated that ZmHSL plays a key role in lignin synthesis, and its defective mutation causes changes in the cell wall structure and gene expression patterns, which impedes water transport and confers higher sensitivity to high-temperature stress.

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.

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.