Whole-genome resequencing of Cucurbita pepo morphotypes to discover genomic variants associated with morphology and horticulturally valuable traits

The effect of environmental factors on the genetic differentiation of Cucurbita ficifolia populations based on whole-genome resequencing

BACKGROUND

Cucurbita ficifolia is one of the squash species most resistant to fungal pathogens, and has especially high resistance to melon Fusarium wilt. This species is therefore an important germplasm resource for the breeding of squash and melon cultivars.

RESULTS

Whole-genome resequencing of 223 individuals from 32 populations in Yunnan Province, the main cucurbit production area in China, was performed and 3,855,120 single-nucleotide polymorphisms (SNPs) and 1,361,000 InDels were obtained. SNP analysis suggested that levels of genetic diversity in C. ficifolia were high, but that different populations showed no significant genetic differentiation or geographical structure, and that individual C. ficifolia plants with fruit rinds of a similar color did not form independent clusters. A Mantel test conducted in combination with geographical distance and environmental factors suggested that genetic distance was not correlated with geographical distance, but had a significant correlation with environmental distance. Further associations between the genetic data and five environmental factors were analyzed using whole-genome association analysis. SNPs associated with each environmental factor were investigated and genes 250 kb upstream and downstream from associated SNPs were annotated. Overall, 15 marker-trait-associated SNPs (MTAs) and 293 genes under environmental selection were identified. The identified genes were involved in cell membrane lipid metabolism, macromolecular complexes, catalytic activity and other related aspects. Ecological niche modeling was used to simulate the distribution of C. ficifolia across time, from the present and into the future. We found that the area suitable for C. ficifolia changed with the changing climate in different periods.

CONCLUSIONS

Resequencing of the C. ficifolia accessions has allowed identification of genetic markers, such as SNPs and InDels. The SNPs identified in this study suggest that environmental factors mediated the formation of the population structure of C. ficifolia in China. These SNPs and Indels might also contribute to the variation in important pathways of genes for important agronomic traits such as yield, disease resistance and stress tolerance. Moreover, the genome resequencing data and the genetic markers identified from 223 accessions provide insight into the genetic variation of the C. ficifolia germplasm and will facilitate a broad range of genetic studies.

Characterization of the USDA Cucurbita pepo, C. moschata, and C. maxima germplasm collections

The Cucurbita genus is home to a number of economically and culturally important species. We present the analysis of genotype data generated through genotyping-by-sequencing of the USDA germplasm collections of Cucurbita pepo, C. moschata, and C. maxima. These collections include a mixture of wild, landrace, and cultivated specimens from all over the world. Roughly 1,500 - 32,000 high-quality single nucleotide polymorphisms (SNPs) were called in each of the collections, which ranged in size from 314 to 829 accessions. Genomic analyses were conducted to characterize the diversity in each of the species. Analysis revealed extensive structure corresponding to a combination of geographical origin and morphotype/market class. Genome-wide associate studies (GWAS) were conducted using both historical and contemporary data. Signals were observed for several traits, but the strongest was for the bush (Bu) gene in C. pepo. Analysis of genomic heritability, together with population structure and GWAS results, was used to demonstrate a close alignment of seed size in C. pepo, maturity in C. moschata, and plant habit in C. maxima with genetic subgroups. These data represent a large, valuable collection of sequenced Cucurbita that can be used to direct the maintenance of genetic diversity, for developing breeding resources, and to help prioritize whole-genome re-sequencing.

Phytophthora capsici Populations Are Structured by Host, Geography, and Fluopicolide Sensitivity

Phytophthora capsici epidemics are propelled by warm temperatures and wet conditions. With temperatures and inland flooding in many locations worldwide expected to rise as a result of global climate change, understanding of population structure can help to inform management of P. capsici in the field and prevent devastating epidemics. Thus, we investigated the effect of host crop, geographical origin, fungicide sensitivity, and mating type on shaping the population structure of P. capsici in the eastern United States. Our fungicide in vitro assays identified the emergence of insensitive isolates for fluopicolide and mefenoxam. A set of 12 microsatellite markers proved informative to assign 157 P. capsici isolates to five distinct genetic clusters. Implementation of Bayesian structure, population differentiation, genetic diversity statistics, and index of association analysis, allowed us to identify population structure by host with some correspondence with genetic clusters for cucumber and squash isolates. We found weak population structure by state for geographically close isolates. In this study, we discovered that North Carolina populations stratify by fluopicolide sensitivity with insensitive isolates experiencing nonrandom mating. Our findings highlight the need for careful monitoring of local field populations, improved selection of relevant isolates for breeding efforts, and hypervigilant surveillance of resistance to different fungicides.

Ethylene plays a dual role in sex determination and fruit shape in cucurbits

Shapes of vegetables and fruits are the result of adaptive evolution and human selection. Modules controlling organ shape have been identified. However, little is known about signals coordinating organ development and shape. Here, we describe the characterization of a melon mutation rf1, leading to round fruit. Histological analysis of rf1 flower and fruits revealed fruit shape is determined at flower stage 8, after sex determination and before flower fertilization. Using positional cloning, we identified the causal gene as the monoecy sex determination gene CmACS7, and survey of melon germplasms showed strong association between fruit shape and sexual types. We show that CmACS7-mediated ethylene production in carpel primordia enhances cell expansion and represses cell division, leading to elongated fruit. Cell size is known to rise as a result of endoreduplication. At stage 8 and anthesis, we found no variation in ploidy levels between female and hermaphrodite flowers, ruling out endoreduplication as a factor in fruit shape determination. To pinpoint the gene networks controlling elongated versus round fruit phenotype, we analyzed the transcriptomes of laser capture microdissected carpels of wild-type and rf1 mutant. These high-resolution spatiotemporal gene expression dynamics revealed the implication of two regulatory modules. The first module implicates E2F-DP transcription factors, controlling cell elongation versus cell division. The second module implicates OVATE- and TRM5-related proteins, controlling cell division patterns. Our finding highlights the dual role of ethylene in the inhibition of the stamina development and the elongation of ovary and fruit in cucurbits.

Multiple Stressors in Vegetable Production: Insights for Trait-Based Crop Improvement in Cucurbits

Vegetable production is a key determinant of contribution from the agricultural sector toward national Gross Domestic Product in a country like India, the second largest producer of fresh vegetables in the world. This calls for a careful scrutiny of the threats to vegetable farming in the event of climate extremes, environmental degradation and incidence of plant pests/diseases. Cucurbits are a vast group of vegetables grown almost throughout the world, which contribute to the daily diet on a global scale. Increasing food supply to cater to the ever-increasing world population, calls for intensive, off-season and year-round cultivation of cucurbits. Current situation predisposes these crops to a multitude of stressors, often simultaneously, under field conditions. This scenario warrants a systematic understanding of the different stress specific traits/mechanisms/pathways and their crosstalk that have been examined in cucurbits and identification of gaps and formulation of perspectives on prospective research directions. The careful dissection of plant responses under specific production environments will help in trait identification for genotype selection, germplasm screens to identify superior donors or for direct genetic manipulation by modern tools for crop improvement. Cucurbits exhibit a wide range of acclimatory responses to both biotic and abiotic stresses, among which a few like morphological characters like waxiness of cuticle; primary and secondary metabolic adjustments; membrane thermostability, osmoregulation and, protein and reactive oxygen species homeostasis and turnover contributing to cellular tolerance, appear to be common and involved in cross talk under combinatorial stress exposures. This is assumed to have profound influence in triggering system level acclimation responses that safeguard growth and metabolism. The possible strategies attempted such as grafting initiatives, molecular breeding, novel genetic manipulation avenues like gene editing and ameliorative stress mitigation approaches, have paved way to unravel the prospects for combined stress tolerance. The advent of next generation sequencing technologies and big data management of the omics output generated have added to the mettle of such emanated concepts and ideas. In this review, we attempt to compile the progress made in deciphering the biotic and abiotic stress responses of cucurbits and their associated traits, both individually and in combination.

Population Evolution, Genetic Diversity and Structure of the Medicinal Legume, Glycyrrhiza uralensis and the Effects of Geographical Distribution on Leaves Nutrient Elements and Photosynthesis

Glycyrrhiza uralensis is a valuable medicinal legume, which occurs widely in arid and semi-arid regions. G. uralensis demand has risen steeply due to its high medical and commercial value. Interpret genome-wide information can stimulate the G. uralensis development as far as its increased bioactive compound levels, and plant yield are concerned. In this study, leaf nutrient concentration and photosynthetic chlorophyll index of G. uralensis were evaluated to determine the G. uralensis growth physiology in three habitats. We observed that G. uralensis nutrient levels and photosynthesis differed significantly in three regions (p < 0.05). Whole-genome re-sequencing of the sixty G. uralensis populations samples from different habitats was performed using an Illumina HiSeq sequencing platform to elucidate the distribution patterns, population evolution, and genetic diversity of G. uralensis. 150.06 Gb high-quality clean data was obtained after strict filtering. The 895237686 reads were mapped against the reference genome, with an average 89.7% mapping rate and 87.02% average sample coverage rate. A total of 6985987 SNPs were identified, and 117970 high-quality SNPs were obtained after filtering, which were subjected to subsequent analysis. Principal component analysis (PCA) based on interindividual SNPs and phylogenetic analysis based on interindividual SNPs showed that the G. uralensis samples could be categorized into central, southern, and eastern populations, which reflected strong genetic differentiation due to long periods of geographic isolation. In this study, a total of 131 candidate regions were screened, and 145 candidate genes (such as Glyur001802s00036258, Glyur003702s00044485, Glyur001802s00036257, Glyur007364s00047495, Glyur000028s00003476, and Glyur000398s00034457) were identified by selective clearance analysis based on Fst and θπ values. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed significant enrichment of 110 GO terms including carbohydrate metabolic process, carbohydrate biosynthetic process, carbohydrate derivative biosynthetic process, and glucose catabolic process (p < 0.05). Alpha-linolenic acid metabolism, biosynthesis of unsaturated fatty acids, and fatty acid degradation pathways were significantly enriched (p < 0.05). This study provides information on the genetic diversity, genetic structure, and population adaptability of the medicinal legumes, G. uralensis. The data obtained in this study provide valuable information for plant development and future optimization of breeding programs for functional genes.

Large-scale gene gains and losses molded the NLR defense arsenal during the Cucurbita evolution

Genome-wide annotation reveals that the gene birth-death process of the Cucurbita R family is associated with a species-specific diversification of TNL and CNL protein classes. The Cucurbitaceae family includes nearly 1000 plant species known universally as cucurbits. Cucurbita genus includes many economically important worldwide crops vulnerable to more than 200 pathogens. Therefore, the identification of pathogen-recognition genes is of utmost importance for this genus. The major class of plant-resistance (R) genes encodes nucleotide-binding site and leucine-rich repeat (NLR) proteins, and is divided into three sub-classes namely, TIR-NB-LRR (TNL), CC-NB-LRR (CNL) and RPW8-NB-LRR (RNL). Although the characterization of the NLR gene family has been carried out in important Cucurbita species, this information is still linked to the availability of sequenced genomes. In this study, we analyzed 40 de novo transcriptomes and 5 genome assemblies, which were explored to investigate the Cucurbita expressed-NLR (eNLR) and NLR repertoires using an ad hoc gene annotation approach. Over 1850 NLR-encoding genes were identified, finely characterized and compared to 96 well-characterized plant R-genes. The maximum likelihood analyses revealed an unusual diversification of CNL/TNL genes and a strong RNL conservation. Indeed, several gene gain and loss events have shaped the Cucurbita NLR family. Finally, to provide a first validation step Cucurbita, eNLRs were explored by real-time PCR analysis. The NLR repertories of the 12 Cucurbita species presented in this paper will be useful to discover novel R-genes.

Genetic Resources and Vulnerabilities of Major Cucurbit Crops

The Cucurbitaceae family provides numerous important crops including watermelons (Citrullus lanatus), melons (Cucumis melo), cucumbers (Cucumis sativus), and pumpkins and squashes (Cucurbita spp.). Centers of domestication in Africa, Asia, and the Americas were followed by distribution throughout the world and the evolution of secondary centers of diversity. Each of these crops is challenged by multiple fungal, oomycete, bacterial, and viral diseases and insects that vector disease and cause feeding damage. Cultivated varieties are constrained by market demands, the necessity for climatic adaptations, domestication bottlenecks, and in most cases, limited capacity for interspecific hybridization, creating narrow genetic bases for crop improvement. This analysis of crop vulnerabilities examines the four major cucurbit crops, their uses, challenges, and genetic resources. ex situ germplasm banks, the primary strategy to preserve genetic diversity, have been extensively utilized by cucurbit breeders, especially for resistances to biotic and abiotic stresses. Recent genomic efforts have documented genetic diversity, population structure, and genetic relationships among accessions within collections. Collection size and accessibility are impacted by historical collections, current ability to collect, and ability to store and maintain collections. The biology of cucurbits, with insect-pollinated, outcrossing plants, and large, spreading vines, pose additional challenges for regeneration and maintenance. Our ability to address ongoing and future cucurbit crop vulnerabilities will require a combination of investment, agricultural, and conservation policies, and technological advances to facilitate collection, preservation, and access to critical Cucurbitaceae diversity.

Combined analysis of carotenoid metabolites and the transcriptome to reveal the molecular mechanism underlying fruit colouration in zucchini (Cucurbita pepo L.)

•

14 Carotenoids were detected and quantified in three Zucchini fruits by LC-MS/MS.

•

Plenty Lutein and less Chlorophyll were identified in yellow and orange fruits.

•

A total of 664 DEGs were identified by transcriptome analysis.

•

Hub DEGs related to Chlorophyll and Carotenoids were associated by WGCNA.

•

A model of gene regulatory network was proposed for Zucchini peel coloration.

To reveal the molecular mechanism underlying peel colouration, carotenoid metabolites and the transcriptome were jointly analysed in zucchini peels with three different colours: light green (Lg), yellow (Y), and orange (O). Our results showed that the carotenoid levels in O (157.075 μg/g) and Y (22.734 μg/g) were both significantly higher than in Lg (7.435 μg/g), while the chlorophyll content was highest in Lg (32.326 μg/g), followed by O (7.294 μg/g) and Y (4.617 μg/g). A total of 14 carotenoids were detected in zucchini peels, primarily lutein (103.167 μg/g in Lg, 509.667 μg/g in Y, and 1543.333 μg/g in O). In particular, significant accumulation of antheraxanthin, zeaxanthin, neoxanthin, and β-cryptoxanthin was first reported in orange zucchini in this study. Furthermore, two modules with hub genes related to carotenoid or chlorophyll content were identified through weighted gene coexpression network analysis. Additionally, the transcription level of some hub genes (PIF4, APRR2, bHLH128, ERF4, PSY1, LCYE2, and RCCR3) was highly correlated with pigment content in the peel, which may be responsible for carotenoid accumulation and chlorophyll degradation in the Y and O varieties. Taken together, the results obtained in this study help to provide a novel mechanism underlying peel colouration in zucchini.

The Distribution and Origins of Pyrus hopeiensis-"Wild Plant With Tiny Population" Using Whole Genome Resequencing

Pyrus hopeiensis is a valuable but endangered wild resource in the genus Pyrus. It has been listed as one of the 120 wild species with tiny population in China. The specie has been little studied. A preliminary study of propagation modes in P. hopeiensis was performed through seed propagation, hybridization, self-crossing trials, bud grafting, branch grafting, and investigations of natural growth. The results showed that the population size of P. hopeiensis was very small, the distribution range was limited, and the habitat was extremely degraded. In the wild population, natural hybridization and root tiller production were the major modes of propagation. Whole genome re-sequencing of the 23 wild and cultivated accessions from Pyrus species collected was performed using an Illumina HiSeq sequencing platform. The sequencing depth range was 26.56x-44.85x and the average sequencing depth was 32x. Phylogenetic tree and principal component analyses (PCA) based on SNPs showed that the wild Pyrus species, such as PWH06, PWH07, PWH09, PWH10, PWH13, and PWH17, were closely related to both P. hopeiensis HB-1 and P. hopeiensis HB-2. Using these results in combination with morphological characteristics, it speculated that P. hopeiensis populations may form a natural hybrid group with frequent gene exchanges between and within groups. A selective elimination analysis on the P. hopeiensis population were performed using Fst and π radio and a total of 381 overlapping genes including SAUR72, IAA20, HSFA2, and RKP genes were obtained. These genes were analyzed by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) function enrichment. And four KEGG pathways, including lysine degradation, sphingolipid metabolism, other glycan degradation, and betaine biosynthesis were significantly enriched in the P. hopeiensis population. Our study provided information on genetic variation, evolutionary relationships, and gene enrichment in P. hopeiensis population. These data will help reveal the evolutionary history and origin of P. hopeiensis and provide guidelines for subsequent research on the locations of functional genes.

Recent and Historical Gene Flow in Cultivars, Landraces, and a Wild Taxon of Cucurbita pepo in Mexico

Gene flow among crops and their wild relatives is an active study area in evolutionary biology and horticulture, because genetic exchange between them may impact their evolutionary trajectories and increase the genetic variation of the cultivated lineages. Mexico is a center of diversity for the genus Cucurbita that includes pumpkins, squash and gourds. Gene flow between domesticated and wild species has been reported as common in Cucurbita; but gene flow among populations of C. pepo ssp. pepo from Mexico and its wild relative has not been studied. We used 2,061 SNPs, derived from tunable genotyping by sequencing (tGBS) to estimate gene flow among 14 Mexican traditional landraces of C. pepo ssp. pepo, also including individuals from five improved cultivars of C. pepo ssp. pepo and C. pepo ssp. ovifera var. ovifera, and individuals of their wild relative C. pepo ssp. fraterna. We found moderate to high levels of genetic diversity, and low to moderate genetic differentiation. In the test of introgression between lineages, we found that all possible arrangements for ancestral and derived sites between the lineages showed similar frequencies; thus, incomplete lineage sorting, but also gene flow, might be taking place in C. pepo. Overall, our results suggest that gene flow between these subspecies and cultigens, incomplete lineage sorting and the retention of ancestral characters shaped the evolutionary trajectory of C. pepo in its area of origin and diversification. In addition, we found evidence of the use of Mexican landraces as genetic material for the improvement of commercial cultivars. The landraces of Mexico are an important source of genetic diversity for C. pepo, which has been preserved both by management practices of small farmers and by the natural gene flow that exists between the different crop fields of the region.

A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms

BACKGROUND

Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots.

RESULTS

In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity.

CONCLUSIONS

These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

The population genomics of adaptive loss of function

Discoveries of adaptive gene knockouts and widespread losses of complete genes have in recent years led to a major rethink of the early view that loss-of-function alleles are almost always deleterious. Today, surveys of population genomic diversity are revealing extensive loss-of-function and gene content variation, yet the adaptive significance of much of this variation remains unknown. Here we examine the evolutionary dynamics of adaptive loss of function through the lens of population genomics and consider the challenges and opportunities of studying adaptive loss-of-function alleles using population genetics models. We discuss how the theoretically expected existence of allelic heterogeneity, defined as multiple functionally analogous mutations at the same locus, has proven consistent with empirical evidence and why this impedes both the detection of selection and causal relationships with phenotypes. We then review technical progress towards new functionally explicit population genomic tools and genotype-phenotype methods to overcome these limitations. More broadly, we discuss how the challenges of studying adaptive loss of function highlight the value of classifying genomic variation in a way consistent with the functional concept of an allele from classical population genetics.

Mining Grapevine Downy Mildew Susceptibility Genes: A Resource for Genomics-Based Breeding and Tailored Gene Editing

Several pathogens continuously threaten viticulture worldwide. Until now, the investigation on resistance loci has been the main trend to understand the interaction between grapevine and the mildew causal agents. Dominantly inherited gene-based resistance has shown to be race-specific in some cases, to confer partial immunity, and to be potentially overcome within a few years since its introgression. Recently, on the footprint of research conducted in Arabidopsis, putative genes associated with downy mildew susceptibility have been discovered also in the grapevine genome. In this work, we deep-sequenced four putative susceptibility genes-namely VvDMR6.1, VvDMR6.2, VvDLO1, VvDLO2-in 190 genetically diverse grapevine genotypes to discover new sources of broad-spectrum and recessively inherited resistance. Identified Single Nucleotide Polymorphisms were screened in a bottleneck analysis from the genetic sequence to their impact on protein structure. Fifty-five genotypes showed at least one impacting mutation in one or more of the scouted genes. Haplotypes were inferred for each gene and two of them at the VvDMR6.2 gene were found significantly more represented in downy mildew resistant genotypes. The current results provide a resource for grapevine and plant genetics and could corroborate genomic-assisted breeding programs as well as tailored gene editing approaches for resistance to biotic stresses.

QTL Mapping for Leaf Area of Tea Plants (Camellia sinensis) Based on a High-Quality Genetic Map Constructed by Whole Genome Resequencing

High-quality genetic maps play important roles in QTL mapping and molecular marker-assisted breeding. Tea leaves are not only important vegetative organs but are also the organ for harvest with important economic value. However, the key genes and genetic mechanism of regulating leaf area have not been clarified. In this study, we performed whole-genome resequencing on "Jinxuan," "Yuncha 1" and their 96 F1 hybrid offspring. From the 1.84 Tb of original sequencing data, abundant genetic variation loci were identified, including 28,144,625 SNPs and 2,780,380 indels. By integrating the markers of a previously reported genetic map, a high-density genetic map consisting of 15 linkage groups including 8,956 high-quality SNPs was constructed. The total length of the genetic map is 1,490.81 cM, which shows good collinearity with the genome. A total of 25 representative markers (potential QTLs) related to leaf area were identified, and there were genes differentially expressed in large and small leaf samples near these markers. GWAS analysis further verified the reliability of QTL mapping. Thirty-one pairs of newly developed indel markers located near these potential QTLs showed high polymorphism and had good discrimination between large and small leaf tea plant samples. Our research will provide necessary support and new insights for tea plant genetic breeding, quantitative trait mapping and yield improvement.

Quantitative Trait Locus Analysis in Squash (Cucurbita moschata) Based on Simple Sequence Repeat Markers and Restriction Site-Associated DNA Sequencing Analysis

Squash (Cucurbita moschata) displays wide morphological and genetic variations; however, limited information is available regarding the genetic loci of squash that control its agronomic traits. To obtain basic genetic information for C. moschata, an F2 population was prepared derived from a cross between the Vietnamese cultivar ‘Bí Hồ Lô TN 6 (TN 6)’ and the Japanese cultivar ‘Shishigatani’, and flowering and fruit traits were examined. Overall, the traits showed a continuous distribution in the F2 population, suggesting that they were quantitative traits. A linkage map was constructed based on simple sequence repeat and restriction site-associated DNA (RAD) markers to detect quantitative trait loci (QTLs). Twelve QTLs for flowering and fruit traits, as well as one phenotypic trait locus, were successfully localized on the map. The present QTLs explained the phenotypic variations at a moderate to relatively high level (16.0%–47.3%). RAD markers linked to the QTLs were converted to codominant cleaved amplified polymorphic sequence (CAPS) and derived CAPS markers for the easy detection of alleles. The information reported here provides useful information for understanding the genetics of Cucurbita and other cucurbit species, and for the selection of individuals with ideal traits during the breeding of Cucurbita vegetables.

Genome-wide SNP discovery and core marker sets for assessment of genetic variations in cultivated pumpkin (Cucurbita spp.)

Three pumpkin species Cucurbita maxima, C. moschata, and C. pepo are commonly cultivated worldwide. To identify genome-wide SNPs in these cultivated pumpkin species, we collected 48 F1 cultivars consisting of 40 intraspecific hybrids (15 C. maxima, 18 C. moschata, and 7 C. pepo) and 8 interspecific hybrids (C. maxima x C. moschata). Genotyping by sequencing identified a total of 37,869 confident SNPs in this collection. These SNPs were filtered to generate a subset of 400 SNPs based on polymorphism and genome distribution. Of the 400 SNPs, 288 were used to genotype an additional 188 accessions (94 F1 cultivars, 50 breeding lines, and 44 landraces) with a SNP array-based platform. Reliable polymorphisms were observed in 224 SNPs (78.0%) and were used to assess genetic variations between and within the four predefined populations in 223 cultivated pumpkin accessions. Both principal component analysis and UPGMA clustering found four major clusters representing three pumpkin species and interspecific hybrids. This genetic differentiation was supported by pairwise Fst and Nei's genetic distance. The interspecific hybrids showed a higher level of genetic diversity relative to the other three populations. Of the 224 SNPs, five subsets of 192, 96, 48, 24, and 12 markers were evaluated for variety identification. The 192, 96, and 48 marker sets identified 204 (91.5%), 190 (85.2%), and 141 (63.2%) of the 223 accessions, respectively, while other subsets showed <25% of variety identification rates. These SNP markers provide a molecular tool with many applications for genetics and breeding in cultivated pumpkin.

Toward the 'golden' era: The status in uncovering the regulatory control of carotenoid accumulation in plants

Carotenoids are essential pigments to plants and important natural products to humans. Carotenoids as both primary and specialized metabolites fulfill multifaceted functions in plants. As such, carotenoid accumulation (a net process of biosynthesis, degradation and sequestration) is subjected to complicated regulation throughout plant life cycle in response to developmental and environmental signals. Investigation of transcriptional regulation of carotenoid metabolic genes remains the focus in understanding the regulatory control of carotenoid accumulation. While discovery of bona fide carotenoid metabolic regulators is still challenging, the recent progress of identification of various transcription factors and regulators helps us to construct hierarchical regulatory network of carotenoid accumulation. The elucidation of carotenoid regulatory mechanisms at protein level and in chromoplast provides some insights into post-translational regulation of carotenogenic enzymes and carotenoid sequestration in plastid sink. This review briefly describes the pathways and main flux-controlling steps for carotenoid accumulation in plants. It highlights our recent understanding of the regulatory mechanisms underlying carotenoid accumulation at both transcriptional and post-translational levels. It also discusses the opportunities to expand toolbox for further shedding light upon the intrinsic regulation of carotenoid accumulation in plants.

Revealing Distinctions in Genetic Diversity and Adaptive Evolution Between Two Varieties of Camellia sinensis by Whole-Genome Resequencing

Camellia sinensis var. sinensis (CSS) and C. sinensis var. assamica (CSA) are the two most economically important tea varieties. They have different characteristics and geographical distribution. Their genetic diversity and differentiation are unclear. Here, we identified 18,903,625 single nucleotide polymorphisms (SNPs) and 7,314,133 insertion-deletion mutations (indels) by whole-genome resequencing of 30 cultivated and three wild related species. Population structure and phylogenetic tree analyses divided the cultivated accessions into CSS and CSA containing 6,440,419 and 6,176,510 unique variations, respectively. The CSS subgroup possessed higher genetic diversity and was enriched for rare alleles. The CSA subgroup had more non-synonymous mutations and might have experienced a greater degree of balancing selection. The evolution rate (dN/dS) and KEGG enrichment indicated that genes involved in the synthesis and metabolism of flavor substances were positively selected in both CSS and CSA subpopulations. However, there are extensive genome differentiation regions (2959 bins and approximately 148 M in size) between the two subgroups. Compared with CSA (141 selected regions containing 124 genes), the CSS subgroup (830 selected regions containing 687 genes) displayed more selection regions potentially related to environmental adaptability. Fifty-three pairs of polymorphic indel markers were developed. Some markers were located in hormone-related genes with distinct alleles in the two cultivated subgroups. These identified variations and selected regions provide clues for the differentiation and adaptive evolution of tea varieties. The newly developed indel markers will be valuable in further genetic research on tea plants.