QTL mapping and GWAS reveal candidate genes controlling capsaicinoid content in Capsicum

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study

BACKGROUND

Phytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, 'PWB-185', 'PWB-186', and '6347', at the 4-8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype.

RESULTS

Most of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers.

CONCLUSIONS

Results will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

Genetic characterization of a locus responsible for low pungency using EMS-induced mutants in Capsicum annuum L

KEY MESSAGE

The pepper mutants ('221-2-1a' and '1559-1-2h') with very low pungency were genetically characterized. The Pun4 locus, responsible for the reduced pungency of the mutant fruits, was localized to a 208 Mb region on chromosome 6. DEMF06G16460, encoding 3-ketoacyl-CoA synthase, was proposed as a strong candidate gene based on the genetic analyses of bulked segregants, DEG, and expression analyses. Capsaicinoids are unique alkaloids present in pepper (Capsicum spp.), synthesized through the condensation of by-products from the phenylpropanoid and branched-chain fatty acid pathways, and accumulating in the placenta. In this study, we characterized two allelic ethyl methanesulfonate-induced mutant lines with extremely low pungency ('221-2-1a' and '1559-1-2h'). These mutants, derived from the pungent Korean landrace 'Yuwolcho,' exhibited lower capsaicinoid content than Yuwolcho but still contained a small amount of capsaicinoid with functional capsaicinoid biosynthetic genes. Genetic crosses between the mutants and Yuwolcho or pungent lines indicated that a single recessive mutation was responsible for the low-pungency phenotype of mutant 221-2-1a; we named the causal locus Pungency 4 (Pun4). To identify Pun4, we combined genome-wide polymorphism analysis and transcriptome analysis with bulked-segregant analysis. We narrowed down the location of Pun4 to a 208-Mb region on chromosome 6 containing five candidate genes, of which DEMF06G16460, encoding a 3-ketoacyl-CoA synthase associated with branched-chain fatty acid biosynthesis, is the most likely candidate for Pun4. The expression of capsaicinoid biosynthetic genes in placental tissues in Yuwolcho and the mutant was consistent with the branched-chain fatty acid pathway playing a pivotal role in the lower pungency observed in the mutant. We also obtained a list of differentially expressed genes in placental tissues between the mutant and Yuwolcho, from which we selected candidate genes using gene co-expression analysis. In summary, we characterized the capsaicinoid biosynthesis-related locus Pun4 through integrated of genetic, genomic, and transcriptome analyses. These findings will contribute to our understanding of capsaicinoid biosynthesis in pepper.

Genomic Analysis of Romanian Lycium Genotypes: Exploring BODYGUARD Genes for Stress Resistance Breeding

Goji berries, long valued in Traditional Chinese Medicine and Asian cuisine for their wide range of medicinal benefits, are now considered a 'superfruit' and functional food worldwide. Because of growing demand, Europe and North America are increasing their goji berry production, using goji berry varieties that are not originally from these regions. European breeding programs are focusing on producing Lycium varieties adapted to local conditions and market demands. By 2023, seven varieties of goji berries were successfully registered in Romania, developed using germplasm that originated from sources outside the country. A broader project focused on goji berry breeding was initiated in 2014 at USAMV Bucharest. In the present research, five cultivated and three wild L. barbarum genotypes were compared to analyse genetic variation at the whole genome level. In addition, a case study presents the differences in the genomic coding sequences of BODYGUARD (BDG) 3 and 4 genes from chromosomes 4, 8, and 9, which are involved in cuticle-related resistance. All three BDG genes show distinctive differences between the cultivated and wild-type genotypes at the SNP level. In the BDG 4 gene located on chromosome 8, 69% of SNPs differentiate the wild from the cultivated genotypes, while in BDG 3 on chromosome 4, 64% of SNPs could tell the difference between the wild and cultivated goji berry. The research also uncovered significant SNP and InDel differences between cultivated and wild genotypes, in the entire genome, providing crucial insights for goji berry breeders to support the development of goji berry cultivation in Romania.

Design principles for synthetic control systems to engineer plants

KEY MESSAGE

Synthetic control systems have led to significant advancement in the study and engineering of unicellular organisms, but it has been challenging to apply these tools to multicellular organisms like plants. The ability to predictably engineer plants will enable the development of novel traits capable of alleviating global problems, such as climate change and food insecurity. Engineering predictable multicellular phenotypes will require the development of synthetic control systems that can precisely regulate how the information encoded in genomes is translated into phenotypes. Many efficient control systems have been developed for unicellular organisms. However, it remains challenging to use such tools to study or engineer multicellular organisms. Plants are a good chassis within which to develop strategies to overcome these challenges, thanks to their capacity to withstand large-scale reprogramming without lethality. Additionally, engineered plants have great potential for solving major societal problems. Here we briefly review the progress of control system development in unicellular organisms, and how that information can be leveraged to characterize control systems in plants. Further, we discuss strategies for developing control systems designed to regulate the expression of transgenes or endogenous loci and generate dosage-dependent or discrete traits. Finally, we discuss the utility that mathematical models of biological processes have for control system deployment.

Complex traits and candidate genes: estimation of genetic variance components across multiple genetic architectures

Large-effect loci-those statistically significant loci discovered by genome-wide association studies or linkage mapping-associated with key traits segregate amidst a background of minor, often undetectable, genetic effects in wild and domesticated plants and animals. Accurately attributing mean differences and variance explained to the correct components in the linear mixed model analysis is vital for selecting superior progeny and parents in plant and animal breeding, gene therapy, and medical genetics in humans. Marker-assisted prediction and its successor, genomic prediction, have many advantages for selecting superior individuals and understanding disease risk. However, these two approaches are less often integrated to study complex traits with different genetic architectures. This simulation study demonstrates that the average semivariance can be applied to models incorporating Mendelian, oligogenic, and polygenic terms simultaneously and yields accurate estimates of the variance explained for all relevant variables. Our previous research focused on large-effect loci and polygenic variance separately. This work aims to synthesize and expand the average semivariance framework to various genetic architectures and the corresponding mixed models. This framework independently accounts for the effects of large-effect loci and the polygenic genetic background and is universally applicable to genetics studies in humans, plants, animals, and microbes.

Genome-wide association study of soybean (Glycine max [L.] Merr.) germplasm for dissecting the quantitative trait nucleotides and candidate genes underlying yield-related traits

Soybean (Glycine max [L.] Merr.) is one of the most significant crops in the world in terms of oil and protein. Owing to the rising demand for soybean products, there is an increasing need for improved varieties for more productive farming. However, complex correlation patterns among quantitative traits along with genetic interactions pose a challenge for soybean breeding. Association studies play an important role in the identification of accession with useful alleles by locating genomic sites associated with the phenotype in germplasm collections. In the present study, a genome-wide association study was carried out for seven agronomic and yield-related traits. A field experiment was conducted in 2015/2016 at two locations that include 155 diverse soybean germplasm. These germplasms were genotyped using SoySNP50K Illumina Infinium Bead-Chip. A total of 51 markers were identified for node number, plant height, pods per plant, seeds per plant, seed weight per plant, hundred-grain weight, and total yield using a multi-locus linear mixed model (MLMM) in FarmCPU. Among these significant SNPs, 18 were putative novel QTNs, while 33 co-localized with previously reported QTLs. A total of 2,356 genes were found in 250 kb upstream and downstream of significant SNPs, of which 17 genes were functional and the rest were hypothetical proteins. These 17 candidate genes were located in the region of 14 QTNs, of which ss715580365, ss715608427, ss715632502, and ss715620131 are novel QTNs for PH, PPP, SDPP, and TY respectively. Four candidate genes, Glyma.01g199200, Glyma.10g065700, Glyma.18g297900, and Glyma.14g009900, were identified in the vicinity of these novel QTNs, which encode lsd one like 1, Ergosterol biosynthesis ERG4/ERG24 family, HEAT repeat-containing protein, and RbcX2, respectively. Although further experimental validation of these candidate genes is required, several appear to be involved in growth and developmental processes related to the respective agronomic traits when compared with their homologs in Arabidopsis thaliana. This study supports the usefulness of association studies and provides valuable data for functional markers and investigating candidate genes within a diverse germplasm collection in future breeding programs.

Genetic Regulation, Environmental Cues, and Extraction Methods for Higher Yield of Secondary Metabolites in Capsicum

Capsicum (chili pepper) is a widely popular and highly consumed fruit crop with beneficial secondary metabolites such as capsaicinoids, carotenoids, flavonoids, and polyphenols, among others. Interestingly, the secondary metabolite profile is a dynamic function of biosynthetic enzymes, regulatory transcription factors, developmental stage, abiotic and biotic environment, and extraction methods. We propose active manipulable genetic, environmental, and extraction controls for the modulation of quality and quantity of desired secondary metabolites in Capsicum species. Specific biosynthetic genes such as Pun (AT3) and AMT in the capsaicinoids pathway and PSY, LCY, and CCS in the carotenoid pathway can be genetically engineered for enhanced production of capsaicinoids and carotenoids, respectively. Generally, secondary metabolites increase with the ripening of the fruit; however, transcriptional regulators such as MYB, bHLH, and ERF control the extent of accumulation in specific tissues. The precise tuning of biotic and abiotic factors such as light, temperature, and chemical elicitors can maximize the accumulation and retention of secondary metabolites in pre- and postharvest settings. Finally, optimized extraction methods such as ultrasonication and supercritical fluid method can lead to a higher yield of secondary metabolites. Together, the integrated understanding of the genetic regulation of biosynthesis, elicitation treatments, and optimization of extraction methods can maximize the industrial production of secondary metabolites in Capsicum.

Fine mapping of a QTL and identification of candidate genes associated with cold tolerance during germination in peanut (Arachis hypogaea L.) on chromosome B09 using whole genome re-sequencing

Low temperatures significantly affect the growth and yield of peanuts. Temperatures lower than 12 °C are generally detrimental for the germination of peanuts. To date, there has been no report on precise information on the quantitative trait loci (QTL) for cold tolerance during the germination in peanuts. In this study, we developed a recombinant inbred line (RIL) population comprising 807 RILs by tolerant and sensitive parents. Phenotypic frequencies of germination rate low-temperature conditions among RIL population showed normally distributed in five environments. Then, we constructed a high density SNP-based genetic linkage map through whole genome re-sequencing (WGRS) technique and identified a major quantitative trait locus (QTL), qRGRB09, on chromosome B09. The cold tolerance-related QTLs were repeatedly detected in all five environments, and the genetic distance was 6.01 cM (46.74 cM - 61.75 cM) after taking a union set. To further confirm that qRGRB09 was located on chromosome B09, we developed Kompetitive Allele Specific PCR (KASP) markers for the corresponding QTL regions. A regional QTL mapping analysis, which was conducted after taking the intersection of QTL intervals of all environments into account, confirmed that qRGRB09 was between the KASP markers, G22096 and G220967 (chrB09:155637831-155854093), and this region was 216.26 kb in size, wherein a total of 15 annotated genes were detected. This study illustrates the relevance of WGRS-based genetic maps for QTL mapping and KASP genotyping that facilitated QTL fine mapping of peanuts. The results of our study also provided useful information on the genetic architecture underlying cold tolerance during germination in peanuts, which in turn may be useful for those engaged in molecular studies as well as crop improvement in the cold-stressed environment.

Large-scale population structure and genetic architecture of agronomic traits of garlic

Garlic, an asexually propagated crop, is the second important bulb crop after the onion and is used as a vegetable and medicinal plant. Abundant and diverse garlic resources have been formed over thousands of years of cultivation. However, genome variation, population structure and genetic architecture of garlic agronomic traits were still not well elucidated. Here, 1 100 258 single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing in 606 garlic accessions collected from 43 countries. Population structure, principal component and phylogenetic analysis showed that these accessions were divided into five subpopulations. Twenty agronomic traits, including above-ground growth traits, bulb-related and bolt-related traits in two consecutive years were implemented in a genome-wide association study. In total, 542 SNPs were associated with these agronomic traits, among which 188 SNPs were repeatedly associated with more than two traits. One SNP (chr6: 1896135972) was repeatedly associated with ten traits. These associated SNPs were located within or near 858 genes, 56 of which were transcription factors. Interestingly, one non-synonymous SNP (Chr4: 166524085) in ribosomal protein S5 was repeatedly associated with above-ground growth and bulb-related traits. Additionally, gene ontology enrichment analysis of candidate genes for genomic selection regions between complete-bolting and non-bolting accessions showed that these genes were significantly enriched in 'vegetative to reproductive phase transition of meristem', 'shoot system development', 'reproductive process', etc. These results provide valuable information for the reliable and efficient selection of candidate genes to achieve garlic genetic improvement and superior varieties.

MYB24 Negatively Regulates the Biosynthesis of Lignin and Capsaicin by Affecting the Expression of Key Genes in the Phenylpropanoid Metabolism Pathway in Capsicum chinense

The wide application of pepper is mostly related to the content of capsaicin, and phenylpropanoid metabolism and its branch pathways may play an important role in the biosynthesis of capsaicin. The expression level of MYB24, a transcription factor screened from the transcriptome data of the pepper fruit development stage, was closely related to the spicy taste. In this experiment, CcMYB24 was cloned from Hainan Huangdenglong pepper, a hot aromatic pepper variety popular in the world for processing, and its function was confirmed by tissue expression characteristics, heterologous transformation in Arabidopsis thaliana, and VIGS technology. The results showed that the relative expression level of CcMYB24 was stable in the early stage of pepper fruit development, and increased significantly from 30 to 50 days after flowering. Heterologous expression led to a significant increase in the expression of CcMYB24 and decrease in lignin content in transgenic Arabidopsis thaliana plants. CcMYB24 silencing led to a significant increase in the expression of phenylpropanoid metabolism pathway genes PAL, 4CL, and pAMT; lignin branch CCR1 and CAD; and capsaicin pathway CS, AT3, and COMT genes in the placenta of pepper, with capsaicin content increased by more than 31.72% and lignin content increased by 20.78%. However, the expression of PAL, pAMT, AT3, COMT, etc., in the corresponding pericarps did not change significantly. Although CS, CCR1, and CAD increased significantly, the relative expression amount was smaller than that in placental tissue, and the lignin content did not change significantly. As indicated above, CcMYB24 may negatively regulate the formation of capsaicin and lignin by regulating the expression of genes from phenylpropanoid metabolism and its branch pathways.

Genetic analysis of pungency deficiency in Japanese chili pepper 'Shishito' (Capsicum annuum) revealed its unique heredity and brought the discovery of two genetic loci involved with the reduction of pungency

The sensation of pungency generated by capsaicinoids is a characteristic trait of chili peppers (Capsicum spp.), and the presence or absence of pungency is central in determining its usage as a spice or a vegetable. In the present study, we aimed to clarify the heredity and genetic factors involved in the deficiency of pungency (quite low pungency) that is uniquely observed in the Japanese chili pepper 'Shishito' (Capsicum annuum). First, the F₂ population ('Shishito' × pungent variety 'Takanotsume') was used for segregation analysis, and pungency level was investigated using capsaicinoid quantification with high-performance liquid chromatography. Also, restriction site associated DNA sequencing of the F₂ population was performed, and genetic map construction and quantitative trait locus (QTL) mapping were implemented. The results indicated that the F₂ population showed varying capsaicinoid content and two major QTLs were detected, Shql3 and Shql7, which explained 39.8 and 19.7% of the genetic variance, respectively. According to these results, the quite low pungency of 'Shishito' was a quantitative trait that involved at least the two loci. Further, this trait was completely separate from general non-pungent traits controlled by individual recessive genes, as described in previous studies. The present study is the first report to investigate the genetic mechanism of pungency deficiency in Japanese chili peppers, and our results provide new insights into the genetic regulation of pungency in chili pepper.

Integrating GWAS, linkage mapping and gene expression analyses reveal the genetic control of first branch height in Brassica napus L

Rapeseed (Brassica napus L.) is a crucial oil crop cultivated worldwide. First branch height, an essential component of rapeseed plant architecture, has an important effect on yield and mechanized harvesting; however, the underlying genetic mechanism remains unclear. In this study, based on the 60K single nucleotide polymorphism array and a recombinant inbred lines population derived from M083 and 888-5, a total of 19 QTLs were detected in five environments, distributed on linkage groups A02, A09, A10, C06, and C07, which explained phenotypic variation ranging from 4.87 to 29.87%. Furthermore, 26 significant SNPs were discovered on Chr.A02 by genome-wide association study in a diversity panel of 324 re-sequencing accessions. The major QTL of the first branch height trait was co-located on Chr.A02 by integrating linkage mapping and association mapping. Eleven candidate genes were screened via allelic variation analysis, inter-subgenomic synteny analysis, and differential expression of genes in parental shoot apical meristem tissues. Among these genes, BnaA02g13010D, which encodes a TCP transcription factor, was confirmed as the target gene according to gene function annotation, haplotype analysis, and full-length gene sequencing, which revealed that TATA insertion/deletion in the promoter region was closely linked to significantly phenotypic differences BnaA02.TCP1 ^M083 overexpression resulted in decreased branch height and increased branch number in Arabidopsis. These results provide a genetic basis for first branch height and the ideal architecture of B. napus.

Pepper variome reveals the history and key loci associated with fruit domestication and diversification

Pepper (Capsicum spp.) is an important vegetable crop that provides a unique pungent sensation when eaten. Through construction of a pepper variome map, we examined the main groups that emerged during domestication and breeding of C. annuum, their relationships and temporal succession, and the molecular events underlying the main transitions. The results showed that the initial differentiation in fruit shape and pungency, increase in fruit weight, and transition from erect to pendent fruits, as well as the recent appearance of large, blocky, sweet fruits (bell peppers), were accompanied by strong selection/fixation of key alleles and introgressions in two large genomic regions. Furthermore, we identified Up, which encodes a BIG GRAIN protein involved in auxin transport, as a key domestication gene that controls erect vs pendent fruit orientation. The up mutation gained increased expression especially in the fruit pedicel through a 579-bp sequence deletion in its 5' upstream region, resulting in the phenotype of pendent fruit. The function of Up was confirmed by virus-induced gene silencing. Taken together, these findings constitute a cornerstone for understanding the domestication and differentiation of a key horticultural crop.

Genetic mapping revealed that the Pun2 gene in Capsicum chacoense encodes a putative aminotransferase

Several genes regulating capsaicinoid biosynthesis including Pun1 (also known as CS), Pun3, pAMT, and CaKR1 have been studied. However, the gene encoded by Pun2 in the non-pungent Capsicum chacoense is unknown. This study aimed to identify the Pun2 gene by genetic mapping using interspecific (C. chacoense × Capsicum annuum) and intraspecific (C. chacoense × C. chacoense) populations. QTL mapping using the interspecific F₂ population revealed two major QTLs on chromosomes 3 and 9. Two bin markers within the QTL regions on two chromosomes were highly correlated with the capsaicinoid content in the interspecific population. The major QTL, Pun2_PJ_Gibbs_3.11 on chromosome 3, contained the pAMT gene, indicating that the non-pungency of C. chacoense may be attributed to a mutation in the pAMT gene. Sequence analysis revealed a 7 bp nucleotide insertion in the 8^th exon of pAMT of the non-pungent C. chacoense. This mutation resulted in the generation of an early stop codon, resulting in a truncated mutant lacking the PLP binding site, which is critical for pAMT enzymatic activity. This insertion co-segregated with the pungency phenotype in the intraspecific F₂ population. We named this novel pAMT allele pamt¹¹ . Taken together, these data indicate that the non-pungency of C. chacoense is due to the non-functional pAMT allele, and Pun2 encodes the pAMT gene.

Mapping of Candidate Genes in Response to Low Nitrogen in Rice Seedlings

Nitrogen is not only a macronutrient essential for crop growth and development, but also one of the most critical nutrients in farmland ecosystem. Insufficient nitrogen supply will lead to crop yield reduction, while excessive application of nitrogen fertilizer will cause agricultural and eco-environment damage. Therefore, mining low-nitrogen tolerant rice genes and improving nitrogen use efficiency are of great significance to the sustainable development of agriculture. This study was conducted by Genome-wide association study on a basis of two root morphological traits (root length and root diameter) and 788,396 SNPs of a natural population of 295 rice varieties. The transcriptome of low-nitrogen tolerant variety (Longjing 31) and low-nitrogen sensitive variety (Songjing 10) were sequenced between low and high nitrogen treatments. A total of 35 QTLs containing 493 genes were mapped. 3085 differential expressed genes were identified. Among these 493 genes, 174 genes showed different haplotype patterns. There were significant phenotype differences among different haplotypes of 58 genes with haplotype differences. These 58 genes were hypothesized as candidate genes for low nitrogen tolerance related to root morphology. Finally, six genes (Os07g0471300, Os11g0230400, Os11g0229300, Os11g0229400, Os11g0618300 and Os11g0229333) which expressed differentially in Longjing 31 were defined as more valuable candidate genes for low-nitrogen tolerance. The results revealed the response characteristics of rice to low-nitrogen, and provided insights into regulatory mechanisms of rice to nitrogen deficiency.

Investigation of genetic factors regulating chlorophyll and carotenoid biosynthesis in red pepper fruit

Chlorophylls and carotenoids are synthesized in the chloroplast and chromoplast, respectively. Even though the two pigments are generated from the same precursor, the genetic correlation between chlorophyll and carotenoid biosynthesis has not yet been fully understood. We investigated the genetic correlation of chlorophyll and carotenoid biosynthesis during fruit ripening. Two recombinant inbred lines populations, "Long Sweet" × "AC2212" ("LA") RILs derived from a cross between Capsicum annuum "Long Sweet" with light-green and light-red fruit and C. annuum "AC2212" with dark-green and brown-fruit and "3501 (F)" × "3509 (C)" ("FC") RILs from C. annuum "3501" with dark-green and dark-red fruit and C. annuum "3509" with intermediate green and light-red fruit, were used. As the fruit ripened, three accessions produced high levels of xanthophyll. The dark-green immature fruit accumulated more total carotenoids than the light-green fruit. This trend corresponded to the expression pattern of 1-deoxy-d-xylulose 5-phosphate synthase (DXS) and CaGLK2 genes during fruit development. The expression levels of DXS and CaGLK2 in the dark-green accession "3501" were significantly higher than those of "3509" and "Long Sweet" during the early stages of fruit development. Furthermore, the genotype analysis of the transcription factor controlling chloroplast development (CaGLK2) in LA RILs revealed that CaGLK2 expression affected both carotenoid and chlorophyll contents. The single nucleotide polymorphism (SNP) linkage maps were constructed using genotyping-by-sequencing (GBS) for the two populations, and QTL analysis was performed for green fruit color intensity and carotenoid content. The QTL (LA_BG-CST10) for capsanthin content in LA RILs located at 24.4 to 100.4 Mbp on chromosome 10 was overlapped with the QTL (FC15-Cap10) for capsanthin content in FC RILs. Three QTLs for capsanthin content, American spice trade association (ASTA) value, and immature green fruit color intensity were also overlapped from 178.2 to 204 Mbp on chromosome 10. At the location, 151.6 to 165 Mbp on chromosome 8, QTLs (FC15-tcar8, FC17-ASTA8.1, and FC17-ASTA8.2) for total carotenoid content and ASTA value were discovered, and this region contained 2-C-methyl-d-erythritol 4-phosphate cytidylyltransferase (MCT), which is involved in the MEP pathway. This result is the first report to show the correlation between carotenoid and chlorophyll biosynthesis in pepper. This research will expand our understanding of the mechanism of the chloroplast-to-chromoplast transition and the development of high pigment pepper varieties.

Genomic selection with fixed-effect markers improves the prediction accuracy for Capsaicinoid contents in Capsicum annuum

Capsaicinoids provide chili peppers (Capsicum spp.) with their characteristic pungency. Several structural and transcription factor genes are known to control capsaicinoid contents in pepper. However, many other genes also regulating capsaicinoid contents remain unknown, making it difficult to develop pepper cultivars with different levels of capsaicinoids. Genomic selection (GS) uses genome-wide random markers (including many in undiscovered genes) for a trait to improve selection efficiency. In this study, we predicted the capsaicinoid contents of pepper breeding lines using several GS models trained with genotypic and phenotypic data from a training population. We used a core collection of 351 Capsicum accessions and 96 breeding lines as training and testing populations, respectively. To obtain the optimal number of single nucleotide polymorphism (SNP) markers for GS, we tested various numbers of genome-wide SNP markers based on linkage disequilibrium. We obtained the highest mean prediction accuracy (0.550) for different models using 3294 SNP markers. Using this marker set, we conducted GWAS and selected 25 markers that were associated with capsaicinoid biosynthesis genes and quantitative trait loci for capsaicinoid contents. Finally, to develop more accurate prediction models, we obtained SNP markers from GWAS as fixed-effect markers for GS, where 3294 genome-wide SNPs were employed. When four to five fixed-effect markers from GWAS were used as fixed effects, the RKHS and RR-BLUP models showed accuracies of 0.696 and 0.689, respectively. Our results lay the foundation for developing pepper cultivars with various capsaicinoid levels using GS for capsaicinoid contents.

Linkage mapping combined with GWAS revealed the genetic structural relationship and candidate genes of maize flowering time-related traits

BACKGROUND

Flowering time is an important agronomic trait of crops and significantly affects plant adaptation and seed production. Flowering time varies greatly among maize (Zea mays) inbred lines, but the genetic basis of this variation is not well understood. Here, we report the comprehensive genetic architecture of six flowering time-related traits using a recombinant inbred line (RIL) population obtained from a cross between two maize genotypes, B73 and Abe2, and combined with genome-wide association studies to identify candidate genes that affect flowering time.

RESULTS

Our results indicate that these six traits showed extensive phenotypic variation and high heritability in the RIL population. The flowering time of this RIL population showed little correlation with the leaf number under different environmental conditions. A genetic linkage map was constructed by 10,114 polymorphic markers covering the whole maize genome, which was applied to QTL mapping for these traits, and identified a total of 82 QTLs that contain 13 flowering genes. Furthermore, a combined genome-wide association study and linkage mapping analysis revealed 17 new candidate genes associated with flowering time.

CONCLUSIONS

In the present study, by using genetic mapping and GWAS approaches with the RIL population, we revealed a list of genomic regions and candidate genes that were significantly associated with flowering time. This work provides an important resource for the breeding of flowering time traits in maize.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08718-y.

Combining QTL mapping and gene co-expression network analysis for prediction of candidate genes and molecular network related to yield in wheat

BACKGROUND

Wheat (Triticum aestivum L.) is an important cereal crop. Increasing grain yield for wheat is always a priority. Due to the complex genome of hexaploid wheat with 21 chromosomes, it is difficult to identify underlying genes by traditional genetic approach. The combination of genetics and omics analysis has displayed the powerful capability to identify candidate genes for major quantitative trait loci (QTLs), but such studies have rarely been carried out in wheat. In this study, candidate genes related to yield were predicted by a combined use of linkage mapping and weighted gene co-expression network analysis (WGCNA) in a recombinant inbred line population.

RESULTS

QTL mapping was performed for plant height (PH), spike length (SL) and seed traits. A total of 68 QTLs were identified for them, among which, 12 QTLs were stably identified across different environments. Using RNA sequencing, we scanned the 99,168 genes expression patterns of the whole spike for the recombinant inbred line population. By the combined use of QTL mapping and WGCNA, 29, 47, 20, 26, 54, 46 and 22 candidate genes were predicted for PH, SL, kernel length (KL), kernel width, thousand kernel weight, seed dormancy, and seed vigor, respectively. Candidate genes for different traits had distinct preferences. The known PH regulation genes Rht-B and Rht-D, and the known seed dormancy regulation genes TaMFT can be selected as candidate gene. Moreover, further experiment revealed that there was a SL regulatory QTL located in an interval of about 7 Mbp on chromosome 7A, named TaSL1, which also involved in the regulation of KL.

CONCLUSIONS

A combination of QTL mapping and WGCNA was applied to predicted wheat candidate genes for PH, SL and seed traits. This strategy will facilitate the identification of candidate genes for related QTLs in wheat. In addition, the QTL TaSL1 that had multi-effect regulation of KL and SL was identified, which can be used for wheat improvement. These results provided valuable molecular marker and gene information for fine mapping and cloning of the yield-related trait loci in the future.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Identification of Genetic Factors Controlling the Formation of Multiple Flowers Per Node in Pepper (Capsicum spp.)

Flower production provides the foundation for crop yield and increased profits. Capsicum annuum is a pepper species with a sympodial shoot structure with solitary flowers. By contrast, C. chinense produces multiple flowers per node. C. annuum accounts for 80% of pepper production worldwide. The identification of C. chinense genes that control multiple flowers and their transfer into C. annuum may open the way to increasing fruit yield. In this study, we dissected the genetic factors were dissected controlling the multiple-flower-per-node trait in Capsicum. 85 recombinant inbred lines (RILs) between the contrasting C. annuum 'TF68' and C. chinense 'Habanero' accessions were phenotyped and genotyped. Quantitative Trait Loci (QTL) analysis identified four novel QTLs on chromosomes 1, 2, 7, and 11 that accounted for 65% of the total phenotypic variation. Genome-wide association study was also performed on a panel of 276 genotyped and phenotyped C. annuum accessions, which revealed 28 regions significantly associated with the multiple-flower trait, of which three overlapped the identified QTLs. Five candidate genes involved in the development of the shoot and flower meristems were identified and these genes could cause multiple flowers per node in pepper. These results contribute to our understanding of multiple flower formation in Capsicum and will be useful to develop high-yielding cultivars.

Chile Pepper (Capsicum) Breeding and Improvement in the "Multi-Omics" Era

Chile pepper (Capsicum spp.) is a major culinary, medicinal, and economic crop in most areas of the world. For more than hundreds of years, chile peppers have "defined" the state of New Mexico, USA. The official state question, "Red or Green?" refers to the preference for either red or the green stage of chile pepper, respectively, reflects the value of these important commodities. The presence of major diseases, low yields, decreased acreages, and costs associated with manual labor limit production in all growing regions of the world. The New Mexico State University (NMSU) Chile Pepper Breeding Program continues to serve as a key player in the development of improved chile pepper varieties for growers and in discoveries that assist plant breeders worldwide. Among the traits of interest for genetic improvement include yield, disease resistance, flavor, and mechanical harvestability. While progress has been made, the use of conventional breeding approaches has yet to fully address producer and consumer demand for these traits in available cultivars. Recent developments in "multi-omics," that is, the simultaneous application of multiple omics approaches to study biological systems, have allowed the genetic dissection of important phenotypes. Given the current needs and production constraints, and the availability of multi-omics tools, it would be relevant to examine the application of these approaches in chile pepper breeding and improvement. In this review, we summarize the major developments in chile pepper breeding and present novel tools that can be implemented to facilitate genetic improvement. In the future, chile pepper improvement is anticipated to be more data and multi-omics driven as more advanced genetics, breeding, and phenotyping tools are developed.

Fruity, sticky, stinky, spicy, bitter, addictive, and deadly: evolutionary signatures of metabolic complexity in the Solanaceae

Covering: 2000-2022Plants collectively synthesize a huge repertoire of metabolites. General metabolites, also referred to as primary metabolites, are conserved across the plant kingdom and are required for processes essential to growth and development. These include amino acids, sugars, lipids, and organic acids. In contrast, specialized metabolites, historically termed secondary metabolites, are structurally diverse, exhibit lineage-specific distribution and provide selective advantage to host species to facilitate reproduction and environmental adaptation. Due to their potent bioactivities, plant specialized metabolites attract considerable attention for use as flavorings, fragrances, pharmaceuticals, and bio-pesticides. The Solanaceae (Nightshade family) consists of approximately 2700 species and includes crops of significant economic, cultural, and scientific importance: these include potato, tomato, pepper, eggplant, tobacco, and petunia. The Solanaceae has emerged as a model family for studying the biochemical evolution of plant specialized metabolism and multiple examples exist of lineage-specific metabolites that influence the senses and physiology of commensal and harmful organisms, including humans. These include, alcohols, phenylpropanoids, and carotenoids that contribute to fruit aroma and color in tomato (fruity), glandular trichome-derived terpenoids and acylsugars that contribute to plant defense (stinky & sticky, respectively), capsaicinoids in chilli-peppers that influence seed dispersal (spicy), and steroidal glycoalkaloids (bitter) from Solanum, nicotine (addictive) from tobacco, as well as tropane alkaloids (deadly) from Deadly Nightshade that deter herbivory. Advances in genomics and metabolomics, coupled with the adoption of comparative phylogenetic approaches, resulted in deeper knowledge of the biosynthesis and evolution of these metabolites. This review highlights recent progress in this area and outlines opportunities for - and challenges of-developing a more comprehensive understanding of Solanaceae metabolism.

Identification of Fruit Traits Related QTLs and a Candidate Gene, CaBRX, Controlling Locule Number in Pepper (Capsicum annuum L.)

Fruit traits are important in pepper (Capsicum annuum L.) and affect its quality and yield. These traits are controlled by quantitative trait loci (QTLs). In this study, we identified many major QTLs that control fruit length (Ftl), fruit diameter (Ftd), fruit shape (Fts), fruit weight (Ftw) and locule number (Lcn) in the F2 and F2:3 populations developed from the QTL mapping of GS6 (P1) and Qiemen (P2). A total of 111 simple sequence repeats and insertion/deletion markers were utilized to construct a linkage map with 12 linkage groups over a length of 1320.72 cM. An inclusive composite interval mapping analysis indicated that many QTLs were detected and included ftl2.1, ftd2.1, fts1.1, ftw2.1 and lcn1.1. As a novel QTL, lcn1.1 was located between HM1112 and EPMS709, and the genetic distance was 3.18 cM covering 60 predicted genes. Within the region, we identified Capana01g004285 as a candidate gene by functional annotation and expression analysis and found that it encodes the BREVIS RADIX (BRX) protein. Knockdown of CaBRX through the virus-induced gene silencing approach in GS6 reduced the number of locules and influenced the expressions of genes related to flower and locule development, suggesting that CaBRX plays an important function in the development of locules.

cis-prenyltransferase 3 and α/β-hydrolase are new determinants of dolichol accumulation in Arabidopsis

Dolichols (Dols), ubiquitous components of living organisms, are indispensable for cell survival. In plants, as well as other eukaryotes, Dols are crucial for post-translational protein glycosylation, aberration of which leads to fatal metabolic disorders in humans and male sterility in plants. Until now, the mechanisms underlying Dol accumulation remain elusive. In this study, we have analysed the natural variation of the accumulation of Dols and six other isoprenoids among more than 120 Arabidopsis thaliana accessions. Subsequently, by combining QTL and GWAS approaches, we have identified several candidate genes involved in the accumulation of Dols, polyprenols, plastoquinone and phytosterols. The role of two genes implicated in the accumulation of major Dols in Arabidopsis-the AT2G17570 gene encoding a long searched for cis-prenyltransferase (CPT3) and the AT1G52460 gene encoding an α/β-hydrolase-is experimentally confirmed. These data will help to generate Dol-enriched plants which might serve as a remedy for Dol-deficiency in humans.

Genetic Correlation Between Fe and Zn Biofortification and Yield Components in a Common Bean (Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) is the most important legume for direct human consumption worldwide. It is a rich and relatively inexpensive source of proteins and micronutrients, especially iron and zinc. Bean is a target for biofortification to develop new cultivars with high Fe/Zn levels that help to ameliorate malnutrition mainly in developing countries. A strong negative phenotypic correlation between Fe/Zn concentration and yield is usually reported, posing a significant challenge for breeders. The objective of this study was to investigate the genetic relationship between Fe/Zn. We used Quantitative Trait Loci (QTLs) mapping and Genome-Wide Association Studies (GWAS) analysis in three bi-parental populations that included biofortified parents, identifying genomic regions associated with yield and micromineral accumulation. Significant negative correlations were observed between agronomic traits (pod harvest index, PHI; pod number, PdN; seed number, SdN; 100 seed weight, 100SdW; and seed per pod, Sd/Pd) and micronutrient concentration traits (SdFe and SdZn), especially between pod harvest index (PHI) and SdFe and SdZn. PHI presented a higher correlation with SdN than PdN. Seventy-nine QTLs were identified for the three populations: 14 for SdFe, 12 for SdZn, 13 for PHI, 11 for SdN, 14 for PdN, 6 for 100SdW, and 9 for Sd/Pd. Twenty-three hotspot regions were identified in which several QTLs were co-located, of which 13 hotpots displayed QTL of opposite effect for yield components and Fe/Zn accumulation. In contrast, eight QTLs for SdFe and six QTLs for SdZn were observed that segregated independently of QTL of yield components. The selection of these QTLs will enable enhanced levels of Fe/Zn and will not affect the yield performance of new cultivars focused on biofortification.

Genome-Wide Association Study of Soybean Germplasm Derived From Canadian × Chinese Crosses to Mine for Novel Alleles to Improve Seed Yield and Seed Quality Traits

Genome-wide association study (GWAS) has emerged in the past decade as a viable tool for identifying beneficial alleles from a genomic diversity panel. In an ongoing effort to improve soybean [Glycine max (L.) Merr.], which is the third largest field crop in Canada, a GWAS was conducted to identify novel alleles underlying seed yield and seed quality and agronomic traits. The genomic panel consisted of 200 genotypes including lines derived from several generations of bi-parental crosses between modern Canadian × Chinese cultivars (CD-CH). The genomic diversity panel was field evaluated at two field locations in Ontario in 2019 and 2020. Genotyping-by-sequencing (GBS) was conducted and yielded almost 32 K high-quality SNPs. GWAS was conducted using Fixed and random model Circulating Probability Unification (FarmCPU) model on the following traits: seed yield, seed protein concentration, seed oil concentration, plant height, 100 seed weight, days to maturity, and lodging score that allowed to identify five QTL regions controlling seed yield and seed oil and protein content. A candidate gene search identified a putative gene for each of the three traits. The results of this GWAS study provide insight into potentially valuable genetic resources residing in Chinese modern cultivars that breeders may use to further improve soybean seed yield and seed quality traits.

Identification of an Allelic Variant of the CsOr Gene Controlling Fruit Endocarp Color in Cucumber (Cucumis sativus L.) Using Genotyping-By-Sequencing (GBS) and Whole-Genome Sequencing

The cucumber is a major vegetable crop around the world. Fruit flesh color is an important quality trait in cucumber and flesh color mainly depends on the relative content of β-carotene in the fruits. The β-carotene serves as a precursor of vitamin A, which has dietary benefits for human health. Cucumbers with orange flesh contain a higher amount of β-carotene than white fruit flesh. Therefore, development of orange-fleshed cucumber varieties is gaining attention for improved nutritional benefits. In this study, we performed genotyping-by-sequencing (GBS) based on genetic mapping and whole-genome sequencing to identify the orange endocarp color gene in the cucumber breeding line, CS-B. Genetic mapping, genetic sequencing, and genetic segregation analyses showed that a single recessive gene (CsaV3_6G040750) encodes a chaperone DnaJ protein (DnaJ) protein at the Cucumis sativus(CsOr) locus was responsible for the orange endocarp phenotype in the CS-B line. The Or gene harbored point mutations T13G and T17C in the first exon of the coding region, resulting in serine to alanine at position 13 and isoleucine to threonine at position 17, respectively. CS-B line displayed increased β-carotene content in the endocarp tissue, corresponding to elevated expression of CsOr gene at fruit developmental stages. Identifying novel missense mutations in the CsOr gene could provide new insights into the role of Or mechanism of action for orange fruit flesh in cucumber and serve as a valuable resource for developing β-carotene-rich cucumbers varieties with increased nutritional benefits.

Genome-wide association study dissects the genetic control of plant height and branch number in response to low-phosphorus stress in Brassica napus

BACKGROUND AND AIMS

Oilseed rape (Brassica napus) is one of the most important oil crops worldwide. Phosphorus (P) deficiency severely decreases the plant height and branch number of B. napus. However, the genetic bases controlling plant height and branch number in B. napus under P deficiency remain largely unknown. This study aims to mine candidate genes for plant height and branch number by genome-wide association study (GWAS) and determine low-P-tolerance haplotypes.

METHODS

An association panel of B. napus was grown in the field with a low P supply (P, 0 kg ha-1) and a sufficient P supply (P, 40 kg ha-1) across 2 years and plant height and branch number were investigated. More than five million single-nucleotide polymorphisms (SNPs) were used to conduct GWAS of plant height and branch number at two contrasting P supplies.

KEY RESULTS

A total of 2127 SNPs were strongly associated (P < 6·25 × 10-07) with plant height and branch number at two P supplies. There was significant correlation between phenotypic variation and the number of favourable alleles of associated loci on chromosomes A10 (chrA10_821671) and C08 (chrC08_27999846), which will contribute to breeding improvement by aggregating these SNPs. BnaA10g09290D and BnaC08g26640D were identified to be associated with chrA10_821671 and chrC08_27999846, respectively. Candidate gene association analysis and haplotype analysis showed that the inbred lines carrying ATT at BnaA10g09290Hap1 and AAT at BnaC08g26640Hap1 had greater plant height than lines carrying other haplotype alleles at low P supply.

CONCLUSION

Our results demonstrate the power of GWAS in identifying genes of interest in B. napus and provided insights into the genetic basis of plant height and branch number at low P supply in B. napus. Candidate genes and favourable haplotypes may facilitate marker-based breeding efforts aimed at improving P use efficiency in B. napus.

Genetic architecture of the metabolic pathway of salicylic acid biosynthesis in Populus

Salicylic acid (SA) is a vital hormone for adaptive responses to biotic and abiotic stresses, which facilitates growth-immunity trade-offs in plants. However, the genetic regulatory networks underlying the metabolic pathway of SA biosynthesis in perennial species remain unclear. Here, we integrated genome-wide association study (GWAS) with metabolite and expression profiling methodologies to dissect the genetic architecture of SA biosynthesis in Populus. First, we quantified nine intermediate metabolites of SA biosynthesis in 300 unrelated Populus tomentosa Carr. individuals. Then, we used a systematic genetic strategy to identify candidate genes for constructing the genetic regulatory network of SA biosynthesis. We focused on WRKY70, an efficient transcription factor, as the key causal gene in the regulatory network, and combined the novel genes coordinating the accumulation of SA. Finally, we identified eight GWAS signals and eight expression quantitative trait loci situated in a selective sweep, and showed the presence of large allele frequency differences among the three geographic populations, revealing that candidate genes subject to selection were involved in SA biosynthesis. This study provides an integrated strategy for dissecting the genetic architecture of the metabolic pathway of SA biosynthesis in Populus, thereby enhancing our understanding of genetic regulation of SA biosynthesis in trees, and accelerating marker-assisted breeding efforts toward high-resistance elite varieties of Populus.

Genotyping-by-Sequencing Derived Genetic Linkage Map and Quantitative Trait Loci for Sugar Content in Onion (Allium cepa L.)

Onion (2n = 2x = 16) has been a nutritional, medicinal and economically valuable vegetable crop all over the world since ancient times. To accelerate the molecular breeding in onion, genetic linkage maps are prerequisite. However, construction of genetic linkage maps of onion remains relatively rudimentary due to a large genome (about 16.3 Gbp) as well as biennial life cycle, cross-pollinated nature, and high inbreeding depression. In this study, we constructed single nucleotide polymorphism (SNP)-based genetic linkage map of onion in an F₂ segregating population derived from a cross between the doubled haploid line '16P118' and inbred line 'Sweet Green' through genotyping by sequencing (GBS). A total of 207.3 Gbp of raw sequences were generated using an Illumina HiSeq X system, and 24,341 SNPs were identified with the criteria based on three minimum depths, lower than 30% missing rate, and more than 5% minor allele frequency. As a result, an onion genetic linkage map consisting of 216 GBS-based SNPs were constructed comprising eight linkage groups spanning a genetic length of 827.0 cM. Furthermore, we identified the quantitative trait loci (QTLs) for the sucrose, glucose, fructose, and total sugar content across the onion genome. We identified a total of four QTLs associated with sucrose (qSC4.1), glucose (qGC5.1), fructose (qFC5.1), and total sugar content (qTSC5.1) explaining the phenotypic variation (R²%) ranging from 6.07-11.47%. This map and QTL information will contribute to develop the molecular markers to breed the cultivars with high sugar content in onion.

Capsicum chinense MYB Transcription Factor Genes: Identification, Expression Analysis, and Their Conservation and Diversification With Other Solanaceae Genomes

Myeloblastosis (MYB) genes are important transcriptional regulators of plant growth, development, and secondary metabolic biosynthesis pathways, such as capsaicinoid biosynthesis in Capsicum. Although MYB genes have been identified in Capsicum annuum, no comprehensive study has been conducted on other Capsicum species. We identified a total of 251 and 240 MYB encoding genes in Capsicum chinense MYBs (CcMYBs) and Capsicum baccatum MYBs (CbMYBs). The observation of twenty tandem and 41 segmental duplication events indicated expansion of the MYB gene family in the C. chinense genome. Five CcMYB genes, i.e., CcMYB101, CcMYB46, CcMYB6, CcPHR8, and CcRVE5, and two CaMYBs, i.e., CaMYB3 and CaHHO1, were found within the previously reported capsaicinoid biosynthesis quantitative trait loci. Based on phylogenetic analysis with tomato MYB proteins, the Capsicum MYBs were classified into 24 subgroups supported by conserved amino acid motifs and gene structures. Also, a total of 241 CcMYBs were homologous with 225 C. annuum, 213 C. baccatum, 125 potato, 79 tomato, and 23 Arabidopsis MYBs. Synteny analysis showed that all 251 CcMYBs were collinear with C. annuum, C. baccatum, tomato, potato, and Arabidopsis MYBs spanning over 717 conserved syntenic segments. Using transcriptome data from three fruit developmental stages, a total of 54 CcMYBs and 81 CaMYBs showed significant differential expression patterns. Furthermore, the expression of 24 CcMYBs from the transcriptome data was validated by quantitative real-time (qRT) PCR analysis. Eight out of the 24 CcMYBs validated by the qRT-PCR were highly expressed in fiery hot C. chinense than in the lowly pungent C. annuum. Furthermore, the co-expression analysis revealed several MYB genes clustered with genes from the capsaicinoid, anthocyanin, phenylpropanoid, carotenoid, and flavonoids biosynthesis pathways, and related to determining fruit shape and size. The homology modeling of 126 R2R3 CcMYBs showed high similarity with that of the Arabidopsis R2R3 MYB domain template, suggesting their potential functional similarity at the proteome level. Furthermore, we have identified simple sequence repeat (SSR) motifs in the CcMYB genes, which could be used in Capsicum breeding programs. The functional roles of the identified CcMYBs could be studied further so that they can be manipulated for Capsicum trait improvement.

Capsaicinoid biosynthesis in the pericarp of chili pepper fruits is associated with a placental septum-like transcriptome profile and tissue structure

CAP biosynthesis in the pericarp of chili pepper fruits occurs with an ambiguous boundary in the placental septum and pericarp. Capsaicinoid (CAP) is a pungent ingredient of chili pepper fruits. Generally, CAP biosynthesis is limited to the placental septum of fruits, but it has been reported that its biosynthesis occurs even in the pericarp of some extremely pungent varieties, resulting in a substantial increase in total content. To examine the mechanism of CAP biosynthesis in the pericarp, comparative transcriptome analysis of a variety that produces CAP in the pericarp (MY) and a variety that does not (HB) was carried out. RNA-seq revealed that 2264 genes were differentially expressed in the MY pericarp compared with the HB pericarp. PCA analysis and GO enrichment analysis indicated that the MY pericarp has a gene expression profile more like placental septum than the HB pericarp. The gene expression of CAP biosynthesis-related genes in the MY pericarp changed coordinately with the placental septum during fruit development. In most Capsicum accessions including HB, the distribution of slender epidermal cells producing CAP was limited to the placental septum, and the morphological boundary between the placental septum and pericarp was clear. In some extremely pungent varieties such as MY, slender epidermal cells ranged from the placental septum to the pericarp region, and the pericarp was morphologically similar to the placental septum, such as the absence of large sub-epidermal cells and abundant spaces in the parenchymal tissue. Our data suggest that CAP biosynthesis in the pericarp occurred with an ambiguous boundary in the placental septum and pericarp. These findings contribute to further enhancement of CAP production in chili pepper fruits.

Evidence for multiple colonisations and Wolbachia infections shaping the genetic structure of the widespread butterfly Polyommatus icarus in the British Isles

The paradigm of isolation in southern refugia during glacial periods followed by expansions during interglacials, producing limited genetic differentiation in northern areas, dominates European phylogeography. However, the existence of complex structured populations in formerly glaciated areas, and islands connected to mainland areas during glacial maxima, call for alternative explanations. We reconstructed the mtDNA phylogeography of the widespread Polyommatus Icarus butterfly with an emphasis on the formerly glaciated and connected British Isles. We found distinct geographical structuring of CO1 haplogroups, with an ancient lineage restricted to the marginal European areas, including Northern Scotland and Outer Hebrides. Population genomic analyses, using ddRADSeq genomic markers, also reveal substantial genetic structuring within Britain. However, there is negligble mito-nuclear concordance consistent with independent demographic histories of mitochondrial versus nuclear DNA. While mtDNA-Wolbachia associations in northern Britain could account for the geographic structuring of mtDNA across most of the British Isles, for nuclear DNA markers (derived from ddRADseq data) butterflies from France cluster between northern and southern British populations - an observation consistent with a scenario of multiple recolonisation. Taken together our results suggest that contemporary mtDNA structuring in the British Isles (and potentially elsewhere in Europe) largely results from Wolbachia infections, however, nuclear genomic structuring suggests a history of at least two distinct colonisations. This two-stage colonisation scenario has previously been put forth to explain genetic diversity and structuring in other British flora and fauna. Additionally, we also present preliminary evidence for potential Wolbachia-induced feminization in the Outer Hebrides.

A consensus and saturated genetic map provides insight into genome anchoring, synteny of Solanaceae and leaf- and fruit-related QTLs in wolfberry (Lycium Linn.)

BACKGROUND

Lycium Linn. (Solanaceae) is a genus of economically important plants producing fruits and leaves with high nutritional value and medicinal benefits. However, genetic analysis of this plant and molecular breeding for quality improvement are limited by the lack of sufficient molecular markers.

RESULTS

In this study, two parental strains, 'Ningqi No. 1' (Lycium barbarum L.) and 'Yunnan Gouqi' (Lycium yunnanense Kuang et A.M. Lu), and 200 F₁ hybrid individuals were resequenced for genetic analysis. In total, 8,507 well-selected SNPs were developed, and a high-density genetic map (NY map) was constructed with a total genetic distance of 2,122.24 cM. A consensus genetic map was established by integrating the NY map and a previously published genetic map (NC map) containing 15,240 SNPs, with a total genetic distance of 3,058.19 cM and an average map distance of 0.21 cM. The 12 pseudochromosomes of the Lycium reference genome were anchored using this consensus genetic map, with an anchoring rate of 64.3%. Moreover, weak collinearities between the consensus map and the pepper, potato, and tomato genomes were observed. Twenty-five stable QTLs were identified for leaf- and fruit-related phenotypes, including fruit weight, fruit longitude, leaf length, the fruit index, and the leaf index; these stable QTLs were mapped to four different linkage groups, with LOD scores ranging from 2.51 to 19.37 and amounts of phenotypic variance explained from 6.2% to 51.9%. Finally, 82 out of 188 predicted genes underlying stable QTLs for fruit-related traits were differentially expressed according to RNA-seq analysis.

CONCLUSIONS

A chromosome-level assembly can provide a foundation for further functional genomics research for wolfberry. The genomic regions of these stably expressed QTLs could be used as targets for further fine mapping and development of molecular markers for marker-assisted selection (MAS). The present study provided valuable information on saturated SNP markers and reliable QTLs for map-based cloning of functional genes related to yield and morphological traits in Lycium spp.

Fine Mapping and Candidate Gene Identification for the CapUp Locus Controlling Fruit Orientation in Pepper (Capsicum spp.)

The orientation of fruits is a distinguishing morphological feature of pepper (Capsicum spp.) varieties. The pendent (downward curved) growth of the fruit stalks, known as pedicels, is highly correlated with fruit weight and pedicel length. A previous genetic analysis revealed that the pendent fruit orientation is governed by a dominant gene, and incomplete inheritance is also observed in some Capsicum accessions. To identify and localize this gene, a single quantitative trait locus (QTL) analysis was performed on one F₂ and two recombinant inbred line (RIL) populations, and a genome-wide association study (GWAS) was performed using a core collection. Common QTL regions associated with fruit orientation were detected on chromosome 12. A total of 187,966 SNPs were identified in a genotyping-by-sequencing (GBS) for GWAS analysis of 196 Capsicum annuum, 25 Capsicum baccatum, 21 Capsicum chinense, and 14 Capsicum frutescens accessions, representing the germplasm collection of South Korea. The results of these analyses enabled us to narrow down the CapUp region of interest to 200-250 Mbp on chromosome 12. Seven candidate genes were found to be located between two markers that were completely cosegregated with the fruit orientation phenotype. The findings and markers developed in this study will be helpful for additional understanding of pepper fruit development and breeding for fruit orientation.

Recent Advances in Molecular Marker-Assisted Breeding for Quality Improvement of Traditional Chinese Medicine

BACKGROUND

The quality of Traditional Chinese Medicine (TCM), reflected by its bioactive compounds and associated contents, is directly linked to its clinical efficacy. Therefore, it is of great importance to improve the quality of TCM by increasing the bioactive compound content.

METHODS

Mapping the active component content-associated QTLs in TCM and further markerassisted breeding has enabled us to rapidly and effectively cultivate new varieties with high bioactive compound contents, which has opened the door for genetic breeding studies on medicinal plants.

RESULTS

In this paper, a strategy and technical molecular breeding method for TCM are discussed. The development of four methods and progress in functional marker development, as well as the applications of such markers in TCM, are reviewed.

CONCLUSION

The progress in, challenges of, and future of marker-assisted breeding for quality improvement of TCM are discussed, which provide valuable scientific references for future molecular breeding.

Prospects for Trifolium Improvement Through Germplasm Characterisation and Pre-breeding in New Zealand and Beyond

Trifolium is the most used pastoral legume genus in temperate grassland systems, and a common feature in meadows and open space areas in cities and parks. Breeding of Trifolium spp. for pastoral production has been going on for over a century. However, the breeding targets have changed over the decades in response to different environmental and production pressures. Relatively small gains have been made in Trifolium breeding progress. Trifolium breeding programmes aim to maintain a broad genetic base to maximise variation. New Zealand is a global hub in Trifolium breeding, utilising exotic germplasm imported by the Margot Forde Germplasm Centre. This article describes the history of Trifolium breeding in New Zealand as well as the role and past successes of utilising genebanks in forage breeding. The impact of germplasm characterisation and evaluation in breeding programmes is also discussed. The history and challenges of Trifolium breeding and its effect on genetic gain can be used to inform future pre-breeding decisions in this genus, as well as being a model for other forage legumes.

A mutation in Zeaxanthin epoxidase contributes to orange coloration and alters carotenoid contents in pepper fruit (Capsicum annuum)

Phytoene synthase (PSY1), capsanthin-capsorubin synthase (CCS), and pseudo-response regulator 2 (PRR2) are three major genes controlling fruit color in pepper (Capsicum spp.). However, the diversity of fruit color in pepper cannot be completely explained by these three genes. Here, we used an F₂ population derived from Capsicum annuum 'SNU-mini Orange' (SO) and C. annuum 'SNU-mini Yellow' (SY), both harboring functional PSY1 and mutated CCS, and observed that yellow color was dominant over orange color. We performed genotyping-by-sequencing and mapped the genetic locus to a 6.8-Mb region on chromosome 2, which we named CaOr. We discovered a splicing mutation in the zeaxanthin epoxidase (ZEP) gene within this region leading to a premature stop codon. HPLC analysis showed that SO contained higher amounts of zeaxanthin and total carotenoids in mature fruits than SY. A color complementation assay using Escherichia coli harboring carotenoid biosynthetic genes showed that the mutant ZEP protein had reduced enzymatic activity. Transmission electron microscopy of plastids revealed that the ZEP mutation affected plastid development with more rod-shaped inner membrane structures in chromoplasts of mature SO fruits. To validate the role of ZEP in fruit color formation, we performed virus-induced gene silencing of ZEP in the yellow-fruit cultivar C. annuum 'Micropep Yellow' (MY). The silencing of ZEP caused significant changes in the ratios of zeaxanthin to its downstream products and increased total carotenoid contents. Thus, we conclude that the ZEP genotype can determine orange or yellow mature fruit color in pepper.

Strategies to Modulate Specialized Metabolism in Mediterranean Crops: From Molecular Aspects to Field

Plant specialized metabolites (SMs) play an important role in the interaction with the environment and are part of the plant defense response. These natural products are volatile, semi-volatile and non-volatile compounds produced from common building blocks deriving from primary metabolic pathways and rapidly evolved to allow a better adaptation of plants to environmental cues. Specialized metabolites include terpenes, flavonoids, alkaloids, glucosinolates, tannins, resins, etc. that can be used as phytochemicals, food additives, flavoring agents and pharmaceutical compounds. This review will be focused on Mediterranean crop plants as a source of SMs, with a special attention on the strategies that can be used to modulate their production, including abiotic stresses, interaction with beneficial soil microorganisms and novel genetic approaches.

The pungent-variable sweet chili pepper 'Shishito' (Capsicum annuum) provides insights regarding the relationship between pungency, the number of seeds, and gene expression involving capsaicinoid biosynthesis

Pungent traits caused by capsaicinoids are characteristic of chili peppers (Capsicum spp.), and the pungent-variable sweet chili pepper 'Shishito' (Capsicum annuum) is unique in being known for the pungency in fruits with few seeds. In the present study, we tried to clarify the relationship between the number of seeds and pungency in 'Shishito'. First, we investigated the pungency of 'Shishito' by simple sensory evaluations and quantifications of capsaicinoids by high-performance liquid chromatography. As a result, few-seeded fruits had a larger fluctuation of capsaicinoid content than many-seeded ones. This indicates that the number of seeds, in particular a decrease of the seeds, has some sort of connection with the pungency of 'Shishito'. Then, we analyzed the relationship between pungency and gene expression involving capsaicinoid biosynthesis at the individual fruit level. We vertically separated the placental septum in which capsaicinoids are synthesized and performed quantitative reverse transcription-polymerase chain reaction (qRT-PCR) for 18 genes involved in capsaicinoid biosynthesis. We also used the placental septum for capsaicinoid quantification so that the gene expression levels and capsaicinoid contents in the same fruits were obtained, and their correlations were analyzed using 20 biological replicates. Among the 18 genes, expression levels of 11 genes (WRKY9, CaMYB31, AT3, BCAT, BCKDH, KAS I, KAS III, ACL, CaKR1, FAT, and pAMT) had a significant positive correlation with the capsaicinoid concentration, and they were considered to upregulate capsaicinoid biosynthesis. These results provide new insights regarding the environmental variation of the pungency traits in chili peppers and the relationship between pungency, the number of seeds, and gene expression involved in capsaicinoid biosynthesis.

An Integrated Approach of QTL Mapping and Genome-Wide Association Analysis Identifies Candidate Genes for Phytophthora Blight Resistance in Sesame (Sesamum indicum L.)

Phytophthora blight (PB) caused by Phytophthora nicotianae is a highly destructive disease in sesame (Sesamum indicum L.). In this study, we used linkage mapping and genome-wide association study (GWAS) to identify quantitative trait loci (QTL) and candidate genes associated with PB resistance. The QTL mapping in 90 RILs of the Goenbaek × Osan cross using genotyping-by-sequencing detected significant QTLs for PB resistance on chromosome 10, explaining 12.79%-13.34% of phenotypic variation. Association of this locus to PB resistance was also revealed through bulked segregant analysis in second RIL population (Goenbaek × Milsung cross) comprising 188 RILs. The GWAS of 87 sesame accessions evaluated against three P. nicotianae isolates identified 29 SNPs on chromosome 10 significantly associated with PB resistance. These SNPs were located within a 0.79 Mb region, which co-located with the QTL intervals identified in RIL populations, and hence scanned for identifying candidate genes. This region contained several defense-related candidate R genes, five of which were selected for quantitative expression analysis. One of these genes, SIN_1019016 was found to show significantly higher expression in the resistant parent compared to that in the susceptible parents and selected RILs. Paired-end sequencing of the gene SIN_1019016 in parental cultivars revealed two synonymous SNPs between Goenbaek and Osan in exon 2 of coding DNA sequence. These results suggested SIN_1019016 as one of the candidate gene conferring PB resistance in sesame. The findings from this study will be useful in the marker-assisted selection as well as the functional analysis of PB resistance candidate gene(s) in sesame.

Neighbor QTL: an interval mapping method for quantitative trait loci underlying plant neighborhood effects

Phenotypes of sessile organisms, such as plants, rely not only on their own genotypes but also on those of neighboring individuals. Previously, we incorporated such neighbor effects into a single-marker regression using the Ising model of ferromagnetism. However, little is known regarding how neighbor effects should be incorporated in quantitative trait locus (QTL) mapping. In this study, we propose a new method for interval QTL mapping of neighbor effects, designated "neighbor QTL," the algorithm of which includes: (1) obtaining conditional self-genotype probabilities with recombination fraction between flanking markers; (2) calculating conditional neighbor genotypic identity using the self-genotype probabilities; and (3) estimating additive and dominance deviations for neighbor effects. Our simulation using F2 and backcross lines showed that the power to detect neighbor effects increased as the effective range decreased. The neighbor QTL was applied to insect herbivory on Col × Kas recombinant inbred lines of Arabidopsis thaliana. Consistent with previous results, the pilot experiment detected a self-QTL effect on the herbivory at the GLABRA1 locus. Regarding neighbor QTL effects on herbivory, we observed a weak QTL on the top of chromosome 4, at which a weak self-bolting QTL was also identified. The neighbor QTL method is available as an R package (https://cran.r-project.org/package=rNeighborQTL), providing a novel tool to investigate neighbor effects in QTL studies.

Chromosome Level Assembly of Homozygous Inbred Line 'Wongyo 3115' Facilitates the Construction of a High-Density Linkage Map and Identification of QTLs Associated With Fruit Firmness in Octoploid Strawberry (Fragaria × ananassa)

Strawberry is an allo-octoploid crop with high genome heterozygosity and complexity, which hinders the sequencing and the assembly of the genome. However, in the present study, we have generated a chromosome level assembly of octoploid strawberry sourced from a highly homozygous inbred line 'Wongyo 3115', using long- and short-read sequencing technologies. The assembly of 'Wongyo 3115' produced 805.6 Mb of the genome with 323 contigs scaffolded into 208 scaffolds with an N50 of 27.3 Mb after further gap filling. The whole genome annotation resulted in 151,892 genes with a gene density of 188.52 (genes/Mb) and validation of a genome, using BUSCO analysis resulted in 94.10% complete BUSCOs. Firmness is one of the vital traits in strawberry, which facilitate the postharvest shelf-life qualities. The molecular and genetic mechanisms that contribute the firmness in strawberry remain unclear. We have constructed a high-density genetic map based on the 'Wongyo 3115' reference genome to identify loci associated with firmness in the present study. For the quantitative trait locus (QTL) identification, the 'BS F₂' populations developed from two inbred lines were genotyped, using an Axiom 35K strawberry chip, and marker positions were analyzed based on the 'Wongyo 3115' genome. Genetic maps were constructed with 1,049 bin markers, spanning the 3,861 cM. Using firmness data of 'BS F₂' obtained from 2 consecutive years, five QTLs were identified on chromosomes 3-3, 5-1, 6-1, and 6-4. Furthermore, we predicted the candidate genes associated with firmness in strawberries by utilizing transcriptome data and QTL information. Overall, we present the chromosome-level assembly and annotation of a homozygous octoploid strawberry inbred line and a linkage map constructed to identify QTLs associated with fruit firmness.

Population Structure and Genetic Diversity of Cucurbita moschata Based on Genome-Wide High-Quality SNPs

Pumpkins (Cucurbita moschata) are one of the most important economic crops in genus Cucurbita worldwide. They are a popular food resource and an important rootstock resource for various Cucurbitaceae. Especially, C. moschata is widely used as a rootstock for the commercial production of bloomless cucumbers in East Asia. Since the genetic diversity of the commercial rootstock varieties is narrow, there has been an increasing demand for the trait development of abiotic and biotic stress tolerance breeding. In this study, 2071 high-quality SNPs that were distributed evenly across 20 chromosomes of pumpkins were discovered through the genotyping-by-sequencing (GBS) analysis of 610 accessions of C. moschata germplasm with a global origin. Using these SNPs, various analyses of the genetic diversity and the population structure were performed. Three subgroups were clustered from the germplasm collection, which included East Asia, Africa, and America, and these areas were included the most in each subgroup. Among those groups, accessions from Africa and South Asia showed the highest genetic diversity, which was followed by the Mexico accessions. This result reflected that large gene pools that consist of various native landraces have been conserved in those of countries. Based on the genetic diversity, we finally constructed the C. moschata core collection, which included 67 representative accessions from the 610 germplasms. Five morphological traits that are important in commercial grafting and rootstock seed production, which include the cotyledon length, the cotyledon width, the hypocotyl length, the internode length, and the number of female flowers, were investigated for three years and used to confirm the validity of the core collection selection. The results are expected to provide valuable information about the genetic structure of the worldwide C. moschata germplasm and help to create new gene pools to develop genetically diverse rootstock breeding materials.

Genome-wide association studies and QTL mapping uncover the genetic architecture of ear tip-barrenness in maize

Ear tip-barrenness (ETB) phenotype threatens crop yield, because it reduces the kernel number per ear. The genetic basis of ETB in maize remains largely unknown. Herein, a genome-wide association study (GWAS) and quantitative trait loci (QTL) mapping were jointly applied to identify the significant genetic loci interrelated with ETB. Six significant SNPs were detected at a stringent P-value threshold (1.95 × 10^-6 ). Additionally, four environment-stable SNPs were co-detected across a single environment and best linear unbiased prediction (BLUP) model at a less stringent P-value threshold (1 × 10^-4 ). The above 10 SNPs were closely linked to 6 candidate genes, which mainly involved seed development, photosynthesis and ethylene response. Moreover, the ratio of superior allele at each significant SNP ranged from 0 to 83.33% in 30 investigated maize elite lines. QTL mapping identified 14 QTL with phenotypic variation explained (PVE) ranging from 3.64 to 7.09%, of which one QTL (qETB2-1) was repeatedly identified in two environments. Combined analysis of GWAS and QTL mapping showed that one SNP (PZE-102175229, chromosome 2: 217 66 Mb) was located in the QTL (qETB2-2, chromosome 2: 215 90-217 82 Mb). Eighteen gene models situated in the linkage disequilibrium (LD) region of the co-localized SNP were further used to evaluate their correlation with ETB by candidate gene association analysis. Two superior haplotypes and two superior alleles were detected among 74 lines for Zm00001d007195, Zm00001d007197 and Zm00001d007201. These results provide more information for clarifying the molecular mechanism of ETB and for speeding up the genetic improvement of maize varieties.