Candidate genes for first flower node identified in pepper using combined SLAF-seq and BSA

Development of a 45K pepper GBTS liquid-phase gene chip and its application in genome-wide association studies

Introduction

Pepper (Capsicum spp.) is a vegetable that is cultivated globally and has undergone extensive domestication, leading to a significant diversification in its agronomic traits. With the advancement of genomics in pepper and the reduction in sequencing costs, the high-throughput detection of single nucleotide polymorphisms (SNPs) and small insertions-deletions (indels) has become increasingly critical for analyzing pepper germplasms and improving breeding programs. As a result, there is a pressing need for a cost-effective, high-throughput, and versatile technique suitable for both foreground and background selection in pepper breeding.

Methods

In the present study, Python-based web scraping scripts were utilized to systematically extract data from published literatures and relevant sequence databases focusing on pepper genomes. Subsequent to data extraction, SNPs and indels were meticulously identified and filtered. This process culminated in the delineation of core polymorphic sites, which were instrumental in the development of specific probes. Following this, comprehensive phenotypic and genotypic analyses were conducted on a diverse collection of 420 pepper germplasms. Concurrently, a genome-wide association study (GWAS) was conducted to elucidate the genetic determinants of helical fruit shape in peppers.

Results

In this study, a 45K pepper Genotyping-By-Target-Sequencing (GBTS) liquid-phase gene chip was developed on the GenoBaits platform. This chip is composed of 45,389 probes, of which 42,535 are derived from core polymorphic sites (CPS) in the background genetic landscape, while 2,854 are associated with foreground agronomic traits, spanning across 43 traits. The CPS probes are spaced at an average interval of 68 Kb. We have assessed the performance of this chip on 420 pepper germplasms, with successful capture of target DNA fragments by 45,387 probes. Furthermore, the probe capture ratio surpassed 70% in 410 of the 420 germplasms tested. Using this chip, we have efficiently genotyped 273 germplasms for spiciness levels and elucidated the genetic relationships among 410 pepper germplasms. Our results allowed for precise clustering of sister lines and C. chinense germplasms. In addition, through a GWAS for helical fruit shape, we identified three quantitative trait loci (QTLs): heli2.1, heli11.1, and heli11.2. Within the heli11.1 QTL, a gene encoding the tubulin alpha chain was identified, suggesting its potential role in the helical growth pattern of pepper fruits.

Discussion

In summary, the 45K pepper GBTS liquid-phase gene chip offers robust detection of polymorphic sites and is a promising tool for advancing research into pepper germplasm and the breeding of new pepper varieties.

Identification of molecular markers and candidate regions associated with grain number per spike in Pubing3228 using SLAF-BSA

Grain number per spike, a pivotal agronomic trait dictating wheat yield, lacks a comprehensive understanding of its underlying mechanism in Pubing3228, despite the identification of certain pertinent genes. Thus, our investigation sought to ascertain molecular markers and candidate regions associated with grain number per spike through a high-density genetic mapping approach that amalgamates site-specific amplified fragment sequencing (SLAF-seq) and bulked segregation analysis (BSA). To facilitate this, we conducted a comparative analysis of two wheat germplasms, Pubing3228 and Jing4839, known to exhibit marked discrepancies in spike shape. By leveraging this methodology, we successfully procured 2,810,474 SLAF tags, subsequently resulting in the identification of 187,489 single nucleotide polymorphisms (SNPs) between the parental strains. We subsequently employed the SNP-index association algorithm alongside the extended distribution (ED) association algorithm to detect regions associated with the trait. The former algorithm identified 24 trait-associated regions, whereas the latter yielded 70. Remarkably, the intersection of these two algorithms led to the identification of 25 trait-associated regions. Amongst these regions, we identified 399 annotated genes, including three genes harboring non-synonymous mutant SNP loci. Notably, the APETALA2 (AP2) transcription factor families, which exhibited a strong correlation with spike type, were also annotated. Given these findings, it is plausible to hypothesize that these genes play a critical role in determining spike shape. In summation, our study contributes significant insights into the genetic foundation of grain number per spike. The molecular markers and candidate regions we have identified can be readily employed for marker-assisted breeding endeavors, ultimately leading to the development of novel wheat cultivars possessing enhanced yield potential. Furthermore, conducting further functional analyses on the identified genes will undoubtedly facilitate a comprehensive elucidation of the underlying mechanisms governing spike development in wheat.

Aspartyl proteases identified as candidate genes of a fiber length QTL, qFLD05, that regulates fiber length in cotton (Gossypium hirsutum L.)

KEY MESSAGE

GhAP genes were identified as the candidates involved in cotton fiber length under the scope of fine mapping a stable fiber length QTL, qFLD05. Moreover, the transcription factor GhWRKY40 positively regulated GhAP3 to decrease fiber length. Fiber length (FL) is an economically important fiber quality trait. Although several genes controlling cotton fiber development have been identified, our understanding of this process remains limited. In this study, an FL QTL (qFLD05) was fine-mapped to a 216.9-kb interval using a secondary F2:3 population derived from the upland hybrid cultivar Ji1518. This mapped genomic segment included 15 coding genes, four of which were annotated as aspartyl proteases (GhAP1-GhAP4). GhAPs were identified as candidates for qFLD05 as the sequence variations in GhAPs were associated with FL deviations in the mapping population, and functional validation of GhAP3 and GhAP4 indicated a longer FL following decreases in their expression levels through virus-induced gene silencing (VIGS). Subsequently, the potential involvement of GhWRKY40 in the regulatory network was revealed: GhWRKY40 positively regulated GhAP3's expression according to transcriptional profiling, VIGS, yeast one-hybrid assays and dual-luciferase experiments. Furthermore, alterations in the expression of the eight previously reported cotton FL-responsive genes from the above three VIGS lines (GhAP3, GhAP4 and GhWRKY40) implied that MYB5_A12 was involved in the GhWRKY40-GhAP network. In short, we unveiled the unprecedented FL regulation roles of GhAPs in cotton, which was possibly further regulated by GhWRKY40. These findings will reveal the genetic basis of FL development associated with qFLD05 and be beneficial for the marker-assisted selection of long-staple cotton.

A point mutation in MC06g1112 encoding FLOWERING LOCUS T decreases the first flower node in bitter gourd (Momordica charantia L.)

In Cucurbitaceae crops, the first flower node (FFN) is an important agronomic trait which can impact the onset of maturity, the production of female flowers, and yield. However, the gene responsible for regulating FFN in bitter gourd is unknown. Here, we used a gynoecious line (S156G) with low FFN as the female parent and a monoecious line (K8-201) with high FFN as the male parent to obtain F1 and F2 generations. Genetic analysis indicated that the low FFN trait was incompletely dominant over the high FFN trait. A major quantitative trait locus (QTL)-Mcffn and four minor effect QTLs-Mcffn1.1, Mcffn1.2, Mcffn1.3, and Mcffn1.4 were detected by whole-genome re-sequencing-based QTL mapping in the S156G×K8-201 F2 population (n=234) cultivated in autumn 2019. The Mcffn locus was further supported by molecular marker-based QTL mapping in three S156G×K8-201 F2 populations planted in autumn 2019 (n=234), autumn 2020 (n=192), and spring 2022 (n=205). Then, the Mcffn locus was fine-mapped into a 77.98-kb physical region on pseudochromosome MC06 using a large S156G×K8-201 F2 population (n=2,402). MC06g1112, which is a homolog of FLOWERING LOCUS T (FT), was considered as the most likely Mcffn candidate gene according to both expression and sequence variation analyses between parental lines. A point mutation (C277T) in MC06g1112, which results in a P93S amino acid mutation between parental lines, may be responsible for decreasing FFN in bitter gourd. Our findings provide a helpful resource for the molecular marker-assisted selective breeding of bitter gourd.

Combined GWAS and Transcriptome Analyses Provide New Insights Into the Response Mechanisms of Sunflower Against Drought Stress

Sunflower is one of the most important oil crops in the world, and drought stress can severely limit its production and quality. To understand the underlying mechanism of drought tolerance, and identify candidate genes for drought tolerance breeding, we conducted a combined genome-wide association studies (GWAS) and RNA-seq analysis. A total of 226 sunflower inbred lines were collected from different regions of China and other countries. Eight phenotypic traits were evaluated under control and drought stress conditions. Genotyping was performed using a Specific-Locus Amplified Fragment Sequencing (SLAF-seq) approach. A total of 934.08 M paired-end reads were generated, with an average Q30 of 91.97%. Based on the 243,291 polymorphic SLAF tags, a total of 94,162 high-quality SNPs were identified. Subsequent analysis of linkage disequilibrium (LD) and population structure in the 226 accessions was carried out based on the 94,162 high-quality SNPs. The average LD decay across the genome was 20 kb. Admixture analysis indicated that the entire population most likely originated from 11 ancestors. GWAS was performed using three methods (MLM, FarmCPU, and BLINK) simultaneously. A total of 80 SNPs showed significant associations with the 8 traits (p < 1.062 × 10-6). Next, a total of 118 candidate genes were found. To obtain more reliable candidate genes, RNA-seq analysis was subsequently performed. An inbred line with the highest drought tolerance was selected according to phenotypic traits. RNA was extracted from leaves at 0, 7, and 14 days of drought treatment. A total of 18,922 differentially expressed genes were obtained. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis showed up-regulated genes were mainly enriched in the branched-chain amino acid catabolic process, while the down-regulated genes were mainly enriched in the photosynthesis-related process. Six DEGs were randomly selected from all DEGs for validation; these genes showed similar patterns in RNA-seq and RT-qPCR analysis, with a correlation coefficient of 0.8167. Through the integration of the genome-wide association study and the RNA-sequencing, 14 candidate genes were identified. Four of them (LOC110885273, LOC110872899, LOC110891369, LOC110920644) were abscisic acid related protein kinases and transcription factors. These genes may play an important role in sunflower drought response and will be used for further study. Our findings provide new insights into the response mechanisms of sunflowers against drought stress and contribute to further genetic breeding.

Broadening the horizon of crop research: a decade of advancements in plant molecular genetics to divulge phenotype governing genes

Advancements in sequencing, genotyping, and computational technologies during the last decade (2011-2020) enabled new forward-genetic approaches, which subdue the impediments of precise gene mapping in varied crops. The modern crop improvement programs rely heavily on two major steps-trait-associated QTL/gene/marker's identification and molecular breeding. Thus, it is vital for basic and translational crop research to identify genomic regions that govern the phenotype of interest. Until the advent of next-generation sequencing, the forward-genetic techniques were laborious and time-consuming. Over the last 10 years, advancements in the area of genome assembly, genotyping, large-scale data analysis, and statistical algorithms have led faster identification of genomic variations regulating the complex agronomic traits and pathogen resistance. In this review, we describe the latest developments in genome sequencing and genotyping along with a comprehensive evaluation of the last 10-year headways in forward-genetic techniques that have shifted the focus of plant research from model plants to diverse crops. We have classified the available molecular genetic methods under bulk-segregant analysis-based (QTL-seq, GradedPool-Seq, QTG-Seq, Exome QTL-seq, and RapMap), target sequence enrichment-based (RenSeq, AgRenSeq, and TACCA), and mutation-based groups (MutMap, NIKS algorithm, MutRenSeq, MutChromSeq), alongside improvements in classical mapping and genome-wide association analyses. Newer methods for outcrossing, heterozygous, and polyploid plant genetics have also been discussed. The use of k-mers has enriched the nature of genetic variants which can be utilized to identify the phenotype-causing genes, independent of reference genomes. We envisage that the recent methods discussed herein will expand the repertoire of useful alleles and help in developing high-yielding and climate-resilient crops.

Rapid Identification of a Stripe Rust Resistance Gene YrXK in Chinese Wheat Line Xike01015 Using Specific Locus Amplified Fragment (SLAF) Sequencing

Wheat stripe rust, an airborne fungal disease caused by Puccinia striiformis Westend. f. sp. tritici, is one of the most devastating diseases of wheat. Chinese wheat cultivar Xike01015 displays high levels of all-stage resistance (ASR) to the current predominant P. striiformis f. sp. tritici race CYR33. In this study, a single dominant gene, designated YrXk, was identified in Xike01015 conferring resistance to CYR33 with genetic analysis of F₂ and BC₁ populations from a cross of Mingxian169 (susceptible) and Xike01015. The specific length amplified fragment sequencing (SLAF-seq) strategy was used to construct a linkage map in the F₂ population. Quantitative trait loci (QTL) analysis mapped YrXk to a 12.4-Mb segment on chromosome1 BS, explaining >86.96% of the phenotypic variance. Gene annotation in the QTL region identified three differential expressed candidate genes, TraesCS1B02G168600.1, TraesCS1B02G170200.1, and TraesCS1B02G172400.1. The qRT-PCR results showed that TraesCS1B02G172400.1 and TraesCS1B02G168600.1 are upregulated and that TraesCS1B02G170200.1 is slightly downregulated after inoculation with CYR33 in the seedling stage, which indicates that these genes may function in wheat resistance to stripe rust. The results of this study can be used in wheat breeding for improving resistance to stripe rust.

Identification of starch candidate genes using SLAF-seq and BSA strategies and development of related SNP-CAPS markers in tetraploid potato

Potato starch is an essential nutrient for humans and is widely used worldwide. Locating relevant genomic regions, mining stable genes and developing candidate gene markers can promote the breeding of new high-starch potato varieties. A total of 106 F₁ individuals and their parents (YSP-4 × MIN-021) were used as test materials, from which 20 plants with high starch content and 20 with low starch content were selected to construct DNA pools for site-specific amplified fragment sequencing (SLAF-seq) and bulked segregation analysis (BSA). A genomic region related to the starch traits was first identified in the 0–5.62 Mb of chromosome 2 in tetraploid potato. In this section, a total of 41 non-synonymous genes, which were considered as candidate genes related to the starch trait, were annotated through a basic local alignment search tool (BLAST) search of multiple databases. Six candidate genes for starch (PGSC0003DMG400017793, PGSC0003DMG400035245, PGSC0003DMG400036713, PGSC0003DMG400040452, PGSC0003DMG400006636 and PGSC0003DMG400044547) were further explored. In addition, cleaved amplified polymorphic sequence (CAPS) markers were developed based on single nucleotide polymorphism (SNP) sites associated with the starch candidate genes. SNP-CAPS markers chr2-CAPS6 and chr2-CAPS21 were successfully developed and validated with the F₂ population and 24 tetraploid potato varieties (lines). Functional analysis and cloning of the candidate genes associated with potato starch will be performed in further research, and the SNP-CAPS markers chr2-CAPS6 and chr2-CAPS21 can be further used in marker-assisted selection breeding of tetraploid potato varieties with high starch content.

The APETALA2 homolog CaFFN regulates flowering time in pepper

Flowering time is an important agronomic trait that contributes to fitness in plants. However, the genetic basis of flowering time has not been extensively studied in pepper. To understand the genetics underlying flowering time, we constructed an F₂ population by crossing a spontaneous early flowering mutant and a late-flowering pepper line. Using bulked segregant RNA-seq, a major locus controlling flowering time in this population was mapped to the end of chromosome 2. An APETALA2 (AP2) homolog (CaFFN) cosegregated with flowering time in 297 individuals of the F₂ population. A comparison between the parents revealed a naturally occurring rare SNP (SNP2T > C) that resulted in the loss of a start codon in CaFFN in the early flowering mutant. Transgenic Nicotiana benthamiana plants with high CaFFN expression exhibited a delay in flowering time and floral patterning defects. On the other hand, pepper plants with CaFFN silencing flowered early. Therefore, the CaFFN gene acts as a flowering repressor in pepper. CaFFN may function as a transcriptional activator to activate the expression of CaAGL15 and miR156e and as a transcriptional repressor to repress the expression of CaAG, CaAP1, CaSEP3, CaSOC1, and miR172b based on a qRT-PCR assay. Direct activation of CaAGL15 by CaFFN was detected using yeast one-hybrid and dual-luciferase reporter assays, consistent with the hypothesis that CaFFN regulates flowering time. Moreover, the CaFFN gene association analysis revealed a significant association with flowering time in a natural pepper population, indicating that the CaFFN gene has a broad effect on flowering time in pepper. Finally, the phylogeny, evolutionary expansion and expression patterns of CaFFN/AP2 homologs were analyzed to provide valuable insight into CaFFN. This study increases our understanding of the involvement of CaFFN in controlling flowering time in pepper, thus making CaFFN a target gene for breeding early maturing pepper.

RNA-seq bulked segregant analysis combined with KASP genotyping rapidly identified PmCH7087 as responsible for powdery mildew resistance in wheat

Powdery mildew causes considerable yield losses in common wheat (Triticum aestivum L.) production. Mapping and cloning powdery mildew-resistant quantitative trait loci can benefit stable yield production by facilitating the breeding of resistant varieties. In this study, we used the powdery mildew resistance introgression line 'CH7087' (harboring the resistance gene PmCH7087) and developed a large F₂ population and a corresponding F_2:3 segregation population containing 2,000 family lines for molecular mapping of PmCH7087. Genetic analysis demonstrated that the resistance phenotype was controlled by a single dominant gene. According to the performance exhibited by the F_2:3 lines, 50 resistant lines and 50 susceptible lines without phenotype segregation were chosen for pooling and bulked segregant RNA sequencing (BSR-Seq) analysis. A region spanning 42.77 Mb was identified, and genotyping of an additional 183 F_2:3 lines with extreme phenotypes using 20 kompetitive allele-specific polymerase chain reaction (KASP) markers in the BSR-Seq mapping regions confirmed this region and narrowed it to 9.68 Mb, in which 45 genes were identified and annotated. Five of these transcripts harbored nonsynonymous single nucleotide polymorphisms between the two parents, with the transcripts of TraesCS2B01G302800 being involved in signal transduction. Furthermore, TraesCS2B01G302800.2 was annotated as the closest homologue of serine/threonine-protein kinase PBS1, a typical participant in the plant disease immune response, indicating that TraesCS2B01G302800 was the candidate gene of PmCH7087. Our results may facilitate future research attempting to improve powdery mildew resistance in wheat and to identify candidate genes for further verification and gene cloning.

Genomic Approaches to Identify Molecular Bases of Crop Resistance to Diseases and to Develop Future Breeding Strategies

Plant diseases are responsible for substantial crop losses each year and affect food security and agricultural sustainability. The improvement of crop resistance to pathogens through breeding represents an environmentally sound method for managing disease and minimizing these losses. The challenge is to breed varieties with a stable and broad-spectrum resistance. Different approaches, from markers to recent genomic and 'post-genomic era' technologies, will be reviewed in order to contribute to a better understanding of the complexity of host-pathogen interactions and genes, including those with small phenotypic effects and mechanisms that underlie resistance. An efficient combination of these approaches is herein proposed as the basis to develop a successful breeding strategy to obtain resistant crop varieties that yield higher in increasing disease scenarios.

Genetic Improvement of Cereals and Grain Legumes

The anticipated population growth by 2050 will be coupled with increased food demand. To achieve higher and sustainable food supplies in order to feed the global population by 2050, a 2.4% rise in the yield of major crops is required. The key to yield improvement is a better understanding of the genetic variation and identification of molecular markers, quantitative trait loci, genes, and pathways related to higher yields and increased tolerance to biotic and abiotic stresses. Advances in genetic technologies are enabling plant breeders and geneticists to breed crop plants with improved agronomic traits. This Special Issue is an effort to report the genetic improvements by adapting genomic techniques and genomic selection.

Morphological, cytological and genetic analyzes of the 'sango' mutant with the defects in basidiocarp development in edible mushroom Pleurotus pulmonarius

A spontaneous, morphological variation 'sango' was observed in the progeny of a Pleurotus pulmonarius (Fr.) Quél. wild-type basidiocarp (also known as fruiting body) collected from the field. This variant developed wart- and coral-like structures instead of normal basidiocarps. Microscopic analysis showed that the sango phenotype had defects in the differentiation of the pileus and hymenium. Basidiocarp phenotypic data analysis in the progenies revealed that the sango trait is a heritable mutation character controlled by a single recessive gene. This mutation locus was mapped on linkage group III of a previously constructed genetic linkage map by amplified fragment length polymorphism (AFLP) technique in P. pulmonarius. Four AFLP markers identified by bulked segregant analysis showed linkage to the sango mutation locus, with the genetic distance ranging from 0 to 2.1 cM. Of these markers, one marker was co-segregated with the sango mutation locus. This knowledge will be a useful foundation for practical breeding as well as for elucidating molecular mechanisms in basidiocarp development of main edible mushrooms.

Genome-based breeding approaches in major vegetable crops

Vegetable crops are major nutrient sources for humanity and have been well-cultivated since thousands of years of domestication. With the rapid development of next-generation sequencing and high-throughput genotyping technologies, the reference genome of more than 20 vegetables have been well-assembled and published. Resequencing approaches on large-scale germplasm resources have clarified the domestication and improvement of vegetable crops by human selection; its application on genetic mapping and quantitative trait locus analysis has led to the discovery of key genes and molecular markers linked to important traits in vegetables. Moreover, genome-based breeding has been utilized in many vegetable crops, including Solanaceae, Cucurbitaceae, Cruciferae, and other families, thereby promoting molecular breeding at a single-nucleotide level. Thus, genome-wide SNP markers have been widely used, and high-throughput genotyping techniques have become one of the most essential methods in vegetable breeding. With the popularization of gene editing technology research on vegetable crops, breeding efficiency can be rapidly increased, especially by combining the genomic and variomic information of vegetable crops. This review outlines the present genome-based breeding approaches used for major vegetable crops to provide insights into next-generation molecular breeding for the increasing global population.

Molecular-Assisted Distinctness and Uniformity Testing Using SLAF-Sequencing Approach in Soybean

Distinctness, uniformity and stability (DUS) testing of cultivars through morphological descriptors is an important and compulsory part of soybean breeding. Molecular markers are usually more effective and accurate in describing the genetic features for the identification and purity assessment of cultivars. In the present study, we assessed the distinctness and uniformity of five soybean cultivars using both single nucleotide polymorphism (SNP) markers developed by specific-locus amplified fragment sequencing (SLAF-seq) technology, and simple sequence repeat (SSR) markers. The phylogenetic tree and principal component analysis (PCA) from both the SLAF-seq and SSR methods showed a clear distinction among cultivars Zhonghuang 18, Zhonghuang 68 and Zhonghuang 35, while no clear distinction was observed between cultivars Zhonghuang 13 and Hedou 13. Using the SLAF-seq method, we determined the proportion of homozygous loci for the five soybean cultivars. The heterozygosity of each individual plant was estimated for the assessment of cultivar purity and the purity levels of the five soybean cultivars ranged from 91.89% to 93.96%. To further validate the applicability of the SLAF-seq approach for distinctness testing, we used the SNP information of 150 soybean cultivars with different origins. The cultivars were also distinguished clearly. Taken together, SLAF-seq can be used as an accurate and reliable method in the assessment of the distinctness and uniformity of soybean cultivars.

Identification of a Major QTL (qRRs-10.1) That Confers Resistance to Ralstonia solanacearum in Pepper (Capsicum annuum) Using SLAF-BSA and QTL Mapping

The soilborne pathogen Ralstonia solanacearum is the causal agent of bacterial wilt (BW), a major disease of pepper (Capsicum annuum). The genetic basis of resistance to this disease in pepper is not well known. This study aimed to identify BW resistance markers in pepper. Analysis of the dynamics of bioluminescent R. solanacearum colonization in reciprocal grafts of a resistant (BVRC 1) line and a susceptible (BVRC 25) line revealed that the resistant rootstock effectively suppressed the spreading of bacteria into the scion. The two clear-cut phenotypic distributions of the disease severity index in 440 F2 plants derived from BVRC 25 × BVRC 1 indicated that a major genetic factor as well as a few minor factors that control BW resistance. By specific-locus amplified fragment sequencing combined with bulked segregant analysis, two adjacent resistance-associated regions on chromosome 10 were identified. Quantitative trait (QTL) mapping revealed that these two regions belong to a single QTL, qRRs-10.1. The marker ID10-194305124, which reached a maximum log-likelihood value at 9.79 and accounted for 19.01% of the phenotypic variation, was located the closest to the QTL peak. A cluster of five predicted R genes and three defense-related genes, which are located in close proximity to the significant markers ID10-194305124 or ID10-196208712, are important candidate genes that may confer BW resistance in pepper.

BRM: a statistical method for QTL mapping based on bulked segregant analysis by deep sequencing

Motivation

Bulked segregant analysis by deep sequencing (BSA-seq) has been widely used for quantitative trait locus (QTL) mapping in recent years. A number of different statistical methods for BSA-seq have been proposed. However, determination of significance threshold, the key point for QTL identification, remains to be a problem that has not been well solved due to the difficulty of multiple testing correction. In addition, estimation of the confidence interval is also a problem to be solved.

Results

In this paper, we propose a new statistical method for BSA-seq, named Block Regression Mapping (BRM). BRM is robust to sequencing noise and is applicable to the case of low sequencing depth. Significance threshold can be reasonably determined by taking multiple testing correction into account. Meanwhile, the confidence interval of QTL position can also be estimated.

Availability and implementation

The R scripts of our method are open source under GPLv3 license at https://github.com/huanglikun/BRM.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-Wide Correlation of 36 Agronomic Traits in the 287 Pepper (Capsicum) Accessions Obtained from the SLAF-seq-Based GWAS

There are many agronomic traits of pepper (Capsicum L.) with abundant phenotypes that can benefit pepper growth. Using specific-locus amplified fragment sequencing (SLAF-seq), a genome-wide association study (GWAS) of 36 agronomic traits was carried out for 287 representative pepper accessions. To ensure the accuracy and reliability of the GWAS results, we analyzed the genetic diversity, distribution of labels (SLAF tags and single nucleotide polymorphisms (SNPs)) and population differentiation and determined the optimal statistical model. In our study, 1487 SNPs were highly significantly associated with 26 agronomic traits, and 2126 candidate genes were detected in the 100-kb region up- and down-stream near these SNPs. Furthermore, 13 major association peaks were identified for 11 key agronomic traits. Then we examined the correlations among the 36 agronomic traits and analyzed SNP distribution and found 37 SNP polymerization regions (total size: 264.69 Mbp) that could be selected areas in pepper breeding. We found that the stronger the correlation between the two traits, the greater the possibility of them being in more than one polymerization region, suggesting that they may be linked or that one pleiotropic gene controls them. These results provide a theoretical foundation for future multi-trait pyramid breeding of pepper. Finally, we found that the GWAS signals were highly consistent with those from the nuclear restorer-of-fertility (Rf) gene for cytoplasmic male sterility (CMS), verifying their reliability. We further identified Capana06g002967 and Capana06g002969 as Rf candidate genes by functional annotation and expression analysis, which provided a reference for the study of cytoplasmic male sterility in Capsicum.

High-density genetic map construction and QTL mapping of first flower node in pepper (Capsicum annuum L.)

BACKGROUND

First flower node (FFN) is an important trait for evaluating fruit earliness in pepper (Capsicum annuum L.). The trait is controlled by quantitative trait loci (QTL); however, studies have been limited on QTL mapping and genes contributing to the trait.

RESULTS

In this study, we developed a high density genetic map using specific-locus amplified fragment sequencing (SLAF-seq), a high-throughput strategy for de novo single nucleotide polymorphism discovery, based on 146 recombinant inbred lines (RILs) derived from an intraspecific cross between PM702 and FS871. The map contained 9328 SLAF markers on 12 linkage groups (LGs), and spanned a total genetic distance of 2009.69 centimorgan (cM) with an average distance of 0.22 cM. The sequencing depth for the map was 72.39-fold in the male parent, 57.04-fold in the female parent, and 15.65-fold in offspring. Using the genetic map, two major QTLs, named Ffn2.1 and Ffn2.2, identified on LG02 were strongly associated with FFN, with a phenotypic variance explanation of 28.62 and 19.56%, respectively. On the basis of the current annotation of C. annuum cv. Criollo de Morelos (CM334), 59 candidate genes were found within the Ffn2.1 and Ffn2.2 region, but only 3 of 59 genes were differentially expressed according to the RNA-seq results. Eventually we identified one gene associated with the FFN based on the function through GO, KEGG, and Swiss-prot analysis.

CONCLUSIONS

Our research showed that the construction of high-density genetic map using SLAF-seq is a valuable tool for fine QTL mapping. The map we constructed is by far the most saturated complete genetic map of pepper, and using it we conducted fine QTL mapping for the important trait, FFN. QTLs and candidate genes obtained in this study lay a good foundation for the further research on FFN-related genes and other genetic applications in pepper.

Construction of a high density genetic map of an interspecific cross of Capsicum chinense and Capsicum annuum and QTL analysis of floral traits

The yield of pepper plants (Capsicum spp.) is their most important trait and is affected by the flower number and flowering time. Capsicum annuum produces a single flower per node and has an early flowering habit. By contrast, Capsicum chinense yields multiple flowers per node and has a late flowering character. However, the genetic mechanism underlying the control of these floral traits remains largely unknown. In this study, 150 F₂ populations from an interspecific cross between the inbred lines 740 (C. chinense) and CA1 (C. annuum) and their parents were used to construct a molecular genetic linkage map using the specific length amplified fragment sequencing (SLAF-seq) technique. This linkage map, spanning 1,586.78 cM in length, contained 9,038 markers on 12 chromosomes, with a mean marker distance of 0.18 cM. Phenotypic data on the flowering time and flower number per node were collected over multiple years, and QTL analysis identified 6 QTLs for the flowering time and flower number per node by composite interval mapping (CIM) and genome-wide composite interval mapping (GCIM) methods at least in two environments. The candidate genes within the major QTL were predicted. In the major flowering time QTL, the candidate gene Capana02g000700, which encodes the homeotic protein APETALA2, was identified. Quantitative reverse-transcription PCR (qRT-PCR) analysis indicated that its expression level in 740 was higher than that in CA1. Gene expression analysis indicated that the expression of Capana02g000700 was significantly upregulated in flowers, and many floral development-related genes were found to be coexpressed with Capana02g000700, supporting the function of this gene in association with flowering time in C. chinense and C. annuum species.

Genetic mapping of anthocyanin accumulation-related genes in pepper fruits using a combination of SLAF-seq and BSA

Anthocyanins have significant functions in stress tolerance in pepper (Capsicum annuum L.) and also benefit human health. Nevertheless, the key structural genes and regulatory genes involved in anthocyanin accumulation in pepper fruits are still not well understood and fine mapped. For the present study, 383 F₂ plants from a cross between the green-fruited C. annuum line Z5 and the purple-fruited line Z6 was developed. Two separate bulked DNA pools were constructed with DNAs extracted from either 37 plants with high anthocyanin content or from 18 plants with no anthocyanin. A combination of specific-locus amplified fragment sequencing (SLAF-seq) and bulked segregant analysis (BSA) was used to identify candidates for regions associated with anthocyanin accumulation. We identified a total of 127,004 high-quality single nucleotide polymorphism (SNP) markers, and detected 1674 high-quality SNP markers associated with anthocyanin accumulation. Three candidate anthocyanin-associated regions including the intervals from 12.48 to 20.00 Mb, from 54.67 to 56.59 Mb, and from 192.17 to 196.82 Mb were identified within a 14.10-Mb interval on chromosome 10 containing 126 candidate genes. Based on their annotations, we identified 12 candidate genes that are predicted to be associated with anthocyanin expression. The present results provide an efficient strategy for genetic mapping of and valuable insights into the genetic mechanisms of anthocyanin accumulation in pepper fruit, and allow us to clone and functionally analyze the genes that influence anthocyanin accumulation in this species.