Construction of a SNP-Based High-Density Genetic Map Using Genotyping by Sequencing (GBS) and QTL Analysis of Nut Traits in Chinese Chestnut (Castanea mollissima Blume)

Comprehensive curation and validation of genomic datasets for chestnut

The Chinese chestnut (Castanea mollissima) stands out as a plant with significant ecological and economic value, excellent nutritional quality and natural resistance to pests and diseases. Recent strides in high-throughput techniques have enabled the continuous accumulation of genomic data on chestnuts, presenting a promising future for genetic research and advancing traits in this species. To facilitate the accessibility and utility of this data, we have curated and analyzed a collection of genomic datasets for eight Castanea species, including functional annotations, 213 RNA-Seq samples, and 330 resequencing samples. These datasets are publicly available on Figshare and are also available through other platforms such as GEO and EVA, providing a valuable resource for researchers studying Castanea genetics, functional genomics, and evolutionary biology. Furthermore, the datasets are integrated into the Castanea Genome Database (CGD, http://castaneadb.net ), which serves as a complementary platform, offering advanced data mining and analysis tools, including BLAST, Batch Query, GO/KEGG Enrichment Analysis, and Synteny Viewer, to enhance the usability of the curated datasets.

Identification and characterization of CmPP2C31 playing a positive role in the abiotic stress resistance of Chinese chestnut via an integrated strategy

Chinese chestnut (Castanea mollissima Blume) is an important economic forest tree species and mainly cultivated in mountainous areas and wastelands, subjecting it to various abiotic stresses. The protein phosphatase 2C (PP2C) genes contributes largely to stress responses in plants. However, the characteristics and functions of PP2C genes in C. mollissima remain unknown. This study provides comprehensive analyses (including phylogenetic, synteny, RNA-seq, transgenic and yeast one-hybrid methods) revealing the characteristics of CmPP2C gene, which plays an important role in response to abiotic stress. Here, we identified 68 CmPP2Cs in the Chinese chestnut genome, and analyzed their characteristics and phylogenetic relationships. Furthermore, synteny analysis revealed that segmental and tandem duplication drove the expansion of the CmPP2C family to adapt to natural environmental pressures. RNA sequencing and co-expression analyses indicated that four hub CmPP2Cs in two key modules probably play important roles in the resistance to abiotic stress in chestnut. Among them, CmPP2C31 was significantly down-regulated under drought stress. Transgenic experiments via pollen magnetofection revealed that CmPP2C31 could positively and significantly regulate the drought resistance of Chinese chestnut seedlings. Subcellular localization showed that CmPP2C31 was a nuclear protein. Yeast one-hybrid assays suggested that EVM0007407 could regulate CmPP2C31 expression by binding to its promoter, thereby participating in abiotic stress resistance. These findings in our study provided detailed information on the CmPP2C family genes and laid a foundation for further elucidating the molecular mechanism of resistance to abiotic stress chestnut.

GWAS identifies two important genes involved in Chinese chestnut weight and leaf length regulation

There are many factors that affect the yield of Chinese chestnut (Castanea mollissima), with single nut weight (SNW) being one of the most important. Leaf length is also related to Chinese chestnut yield. However, the genetic architecture and gene function associated with Chinese chestnut nut yield have not been fully explored. In this study, we performed genotyping by sequencing 151 Chinese chestnut cultivars, followed by a genome-wide association study (GWAS) on six horticultural traits. First, we analyzed the phylogeny of the Chinese chestnut and found that the Chinese chestnut cultivars divided into two ecotypes, a northern and southern cultivar group. Differences between the cultivated populations were found in the pathways of plant growth and adaptation to the environment. In the selected regions, we also found interesting tandemly arrayed genes that may influence Chinese chestnut traits and environmental adaptability. To further investigate which horticultural traits were selected, we performed a GWAS using six horticultural traits from 151 cultivars. Forty-five loci that strongly associated with horticultural traits were identified, and six genes highly associated with these traits were screened. In addition, a candidate gene associated with SNW, APETALA2 (CmAP2), and another candidate gene associated with leaf length (LL), CRYPTOCHROME INTERACTING BASIC HELIX-LOOP-HELIX 1 (CmCIB1), were verified in Chinese chestnut and Arabidopsis (Arabidopsis thaliana). Our results showed that CmAP2 affected SNW by negatively regulating cell size. CmCIB1 regulated the elongation of new shoots and leaves by inducing cell elongation, potentially affecting photosynthesis. This study provided valuable information and insights for Chinese chestnut breeding research.

Marker-Assisted Selection in Breeding for Fruit Trait Improvement: A Review

Breeding fruit species is time-consuming and expensive. With few exceptions, trees are likely the worst species to work with in terms of genetics and breeding. Most are characterized by large trees, long juvenile periods, and intensive agricultural practice, and environmental variability plays an important role in the heritability evaluations of every single important trait. Although vegetative propagation allows for the production of a significant number of clonal replicates for the evaluation of environmental effects and genotype × environment interactions, the spaces required for plant cultivation and the intensity of work necessary for phenotypic surveys slow down the work of researchers. Fruit breeders are very often interested in fruit traits: size, weight, sugar and acid content, ripening time, fruit storability, and post-harvest practices, among other traits relevant to each individual species. The translation of trait loci and whole-genome sequences into diagnostic genetic markers that are effective and affordable for use by breeders, who must choose genetically superior parents and subsequently choose genetically superior individuals among their progeny, is one of the most difficult tasks still facing tree fruit geneticists. The availability of updated sequencing techniques and powerful software tools offered the opportunity to mine tens of fruit genomes to find out sequence variants potentially useful as molecular markers. This review is devoted to analysing what has been the role of molecular markers in assisting breeders in selection processes, with an emphasis on the fruit traits of the most important fruit crops for which examples of trustworthy molecular markers have been developed, such as the MDo.chr9.4 marker for red skin colour in apples, the CCD4-based marker CPRFC1, and LG3_13.146 marker for flesh colour in peaches, papayas, and cherries, respectively.

Genetic Diversity and Population Structure of Chinese Chestnut (Castanea mollissima Blume) Cultivars Revealed by GBS Resequencing

Chinese chestnut (Castanea mollissima Bl.) is one of the earliest domesticated and cultivated fruit trees, and it is widely distributed in China. Because of the high quality of its nuts and its high resistance to abiotic and biotic stresses, Chinese chestnut could be used to improve edible chestnut varieties worldwide. However, the unclear domestication history and highly complex genetic background of Chinese chestnut have prevented the efficiency of breeding efforts. To explore the genetic diversity and structure of Chinese chestnut populations and generate new insights that could aid chestnut breeding, heterozygosity statistics, molecular variance analysis, ADMIXTURE analysis, principal component analysis, and phylogenetic analysis were conducted to analyze single nucleotide polymorphism data from 185 Chinese chestnut landraces from five geographical regions in China via genotyping by sequencing. Results showed that the genetic diversity level of the five populations from different regions was relatively high, with an observed heterozygosity of 0.2796-0.3427. The genetic diversity level of the population in the mid-western regions was the highest, while the population north of the Yellow River was the lowest. Molecular variance analysis showed that the variation among different populations was only 2.07%, while the intra-group variation reached 97.93%. The Chinese chestnut samples could be divided into two groups: a northern and southern population, separated by the Yellow River; however, some samples from the southern population were genetically closer to samples from the northern population. We speculate that this might be related to the migration of humans during the Han dynasty due to the frequent wars that took place during this period, which might have led to the introduction of chestnut to southern regions. Some samples from Shandong Province and Beijing City were outliers that did not cluster with their respective groups, and this might be caused by the special geographical, political, and economic significance of these two regions. The findings of our study showed the complex genetic relationships among Chinese chestnut landraces and the high genetic diversity of these resources.

Genome-Wide SNP Markers Accelerate Perennial Forest Tree Breeding Rate for Disease Resistance through Marker-Assisted and Genome-Wide Selection

The ecological and economic importance of forest trees is evident and their survival is necessary to provide the raw materials needed for wood and paper industries, to preserve the diversity of associated animal and plant species, to protect water and soil, and to regulate climate. Forest trees are threatened by anthropogenic factors and biotic and abiotic stresses. Various diseases, including those caused by fungal pathogens, are one of the main threats to forest trees that lead to their dieback. Genomics and transcriptomics studies using next-generation sequencing (NGS) methods can help reveal the architecture of resistance to various diseases and exploit natural genetic diversity to select elite genotypes with high resistance to diseases. In the last two decades, QTL mapping studies led to the identification of QTLs related to disease resistance traits and gene families and transcription factors involved in them, including NB-LRR, WRKY, bZIP and MYB. On the other hand, due to the limitation of recombination events in traditional QTL mapping in families derived from bi-parental crosses, genome-wide association studies (GWAS) that are based on linkage disequilibrium (LD) in unstructured populations overcame these limitations and were able to narrow down QTLs to single genes through genotyping of many individuals using high-throughput markers. Association and QTL mapping studies, by identifying markers closely linked to the target trait, are the prerequisite for marker-assisted selection (MAS) and reduce the breeding period in perennial forest trees. The genomic selection (GS) method uses the information on all markers across the whole genome, regardless of their significance for development of a predictive model for the performance of individuals in relation to a specific trait. GS studies also increase gain per unit of time and dramatically increase the speed of breeding programs. This review article is focused on the progress achieved in the field of dissecting forest tree disease resistance architecture through GWAS and QTL mapping studies. Finally, the merit of methods such as GS in accelerating forest tree breeding programs is also discussed.

A High-Density Genetic Map Enables Genome Synteny and QTL Mapping of Vegetative Growth and Leaf Traits in Gardenia

The gardenia is a traditional medicinal horticultural plant in China, but its molecular genetic research has been largely hysteretic. Here, we constructed an F₁ population with 200 true hybrid individuals. Using the genotyping-by-sequencing method, a high-density sex-average genetic map was generated that contained 4,249 SNPs with a total length of 1956.28 cM and an average genetic distance of 0.46 cM. We developed 17 SNP-based Kompetitive Allele-Specific PCR markers and found that 15 SNPs were successfully genotyped, of which 13 single-nucleotide polymorphism genotypings of 96 F₁ individuals showed genotypes consistent with GBS-mined genotypes. A genomic collinearity analysis between gardenia and the Rubiaceae species Coffea arabica, Coffea canephora and Ophiorrhiza pumila showed the relativity strong conservation of LG11 with NC_039,919.1, HG974438.1 and Bliw01000011.1, respectively. Lastly, a quantitative trait loci analysis at three phenotyping time points (2019, 2020, and 2021) yielded 18 QTLs for growth-related traits and 31 QTLs for leaf-related traits, of which qBSBN7-1, qCD8 and qLNP2-1 could be repeatably detected. Five QTL regions (qCD8 and qSBD8, qBSBN7 and qSI7, qCD4-1 and qLLLS4, qLNP10 and qSLWS10-2, qSBD10 and qLLLS10) with potential pleiotropic effects were also observed. This study provides novel insight into molecular genetic research and could be helpful for further gene cloning and marker-assisted selection for early growth and development traits in the gardenia.

Chromosome-level genome assembly of Japanese chestnut (Castanea crenata Sieb. et Zucc.) reveals conserved chromosomal segments in woody rosids

Japanese chestnut (Castanea crenata Sieb. et Zucc.), unlike other Castanea species, is resistant to most diseases and wasps. However, genomic data of Japanese chestnut that could be used to determine its biotic stress resistance mechanisms have not been reported to date. In this study, we employed long-read sequencing and genetic mapping to generate genome sequences of Japanese chestnut at the chromosome level. Long reads (47.7 Gb; 71.6× genome coverage) were assembled into 781 contigs, with a total length of 721.2 Mb and a contig N50 length of 1.6 Mb. Genome sequences were anchored to the chestnut genetic map, comprising 14,973 single nucleotide polymorphisms (SNPs) and covering 1,807.8 cM map distance, to establish a chromosome-level genome assembly (683.8 Mb), with 69,980 potential protein-encoding genes and 425.5 Mb repetitive sequences. Furthermore, comparative genome structure analysis revealed that Japanese chestnut shares conserved chromosomal segments with woody plants, but not with herbaceous plants, of rosids. Overall, the genome sequence data of Japanese chestnut generated in this study is expected to enhance not only its genetics and genomics but also the evolutionary genomics of woody rosids.

Genome-wide identification of agronomically important genes in outcrossing crops using OutcrossSeq

Many important crops (e.g., tuber, root, and tree crops) are cross-pollinating. For these crops, no inbred lines are available for genetic study and breeding because they are self-incompatible, clonally propagated, or have a long generation time, making the identification of agronomically important genes difficult, particularly in crops with a complex autopolyploid genome. In this study, we developed a method, OutcrossSeq, for mapping agronomically important loci in outcrossing crops based on whole-genome low-coverage resequencing of a large genetic population, and designed three computation algorithms in OutcrossSeq for different types of outcrossing populations. We applied OutcrossSeq to a tuberous root crop (sweet potato, autopolyploid), a tree crop (walnut tree, highly heterozygous diploid), and hybrid crops (double-cross populations) to generate high-density genotype maps for the outcrossing populations, which enable precise identification of genomic loci underlying important agronomic traits. Candidate causative genes at these loci were detected based on functional clues. Taken together, our results indicate that OutcrossSeq is a robust and powerful method for identifying agronomically important genes in heterozygous species, including polyploids, in a cost-efficient way. The OutcrossSeq software and its instruction manual are available for downloading at www.xhhuanglab.cn/tool/OutcrossSeq.html.

Identification of transcriptome-wide, nut weight-associated SNPs in Castanea crenata

Nut weight is one of the most important traits that can affect a chestnut grower's returns. Due to the long juvenile phase of chestnut trees, the selection of desired characteristics at early developmental stages represents a major challenge for chestnut breeding. In this study, we identified single nucleotide polymorphisms (SNPs) in transcriptomic regions, which were significantly associated with nut weight in chestnuts (Castanea crenata), using a genome-wide association study (GWAS). RNA-sequencing (RNA-seq) data were generated from large and small nut-bearing trees, using an Illumina HiSeq. 2000 system, and 3,271,142 SNPs were identified. A total of 21 putative SNPs were significantly associated with chestnut weight (false discovery rate [FDR] < 10-5), based on further analyses. We also applied five machine learning (ML) algorithms, support vector machine (SVM), C5.0, k-nearest neighbour (k-NN), partial least squares (PLS), and random forest (RF), using the 21 SNPs to predict the nut weights of a second population. The average accuracy of the ML algorithms for the prediction of chestnut weights was greater than 68%. Taken together, we suggest that these SNPs have the potential to be used during marker-assisted selection to facilitate the breeding of large chestnut-bearing varieties.

Hybrid de novo genome assembly of Chinese chestnut (Castanea mollissima)

BACKGROUND

The Chinese chestnut (Castanea mollissima) is widely cultivated in China for nut production. This plant also plays an important ecological role in afforestation and ecosystem services. To facilitate and expand the use of C. mollissima for breeding and its genetic improvement, we report here the whole-genome sequence of C. mollissima.

FINDINGS

We produced a high-quality assembly of the C. mollissima genome using Pacific Biosciences single-molecule sequencing. The final draft genome is ∼785.53 Mb long, with a contig N50 size of 944 kb, and we further annotated 36,479 protein-coding genes in the genome. Phylogenetic analysis showed that C. mollissima diverged from Quercus robur, a member of the Fagaceae family, ∼13.62 million years ago.

CONCLUSIONS

The high-quality whole-genome assembly of C. mollissima will be a valuable resource for further genetic improvement and breeding for disease resistance and nut quality.

A novel high-density grapevine (Vitis vinifera L.) integrated linkage map using GBS in a half-diallel population

A half-diallel population involving five elite grapevine cultivars was generated and genotyped by GBS, and highly-informative segregation data was used to construct a high-density genetic map for Vitis vinifera L. Grapevine is one of the most relevant fruit crops in the world. Deeper genetic knowledge could assist modern grapevine breeding programs to develop new wine grape varieties able to face climate change effects. To assist in the rapid identification of markers for crop yield components, grape quality traits and adaptation potential, we generated a large Vitis vinifera L. population (N = 624) by crossing five red wine cultivars in a half-diallel scheme, which was subsequently sequenced by an efficient GBS procedure. A high number of fully informative genetic variants was detected using a novel mapping approach capable of reconstructing local haplotypes from adjacent biallelic SNPs, which were subsequently used to construct the densest consensus genetic map available for the cultivated grapevine to date. This 1378.3-cM map integrates 10 bi-parental consensus maps and orders 4437 markers in 3353 unique positions on 19 chromosomes. Markers are well distributed all along the grapevine reference genome, covering up to 98.8% of its genomic sequence. Additionally, a good agreement was observed between genetic and physical orders, adding confidence in the quality of this map. Collectively, our results pave the way for future genetic studies (such as fine QTL mapping) aimed to understand the complex relationship between genotypic and phenotypic variation in the cultivated grapevine. In addition, the method used (which efficiently delivers a high number of fully informative markers) could be of interest to other outbred organisms, notably perennial fruit crops.

Genetic diversity, linkage disequilibrium, and population structure analysis of the tea plant (Camellia sinensis) from an origin center, Guizhou plateau, using genome-wide SNPs developed by genotyping-by-sequencing

BACKGROUND

To efficiently protect and exploit germplasm resources for marker development and breeding purposes, we must accurately depict the features of the tea populations. This study focuses on the Camellia sinensis (C. sinensis) population and aims to (i) identify single nucleotide polymorphisms (SNPs) on the genome level, (ii) investigate the genetic diversity and population structure, and (iii) characterize the linkage disequilibrium (LD) pattern to facilitate next genome-wide association mapping and marker-assisted selection.

RESULTS

We collected 415 tea accessions from the Origin Center and analyzed the genetic diversity, population structure and LD pattern using the genotyping-by-sequencing (GBS) approach. A total of 79,016 high-quality SNPs were identified; the polymorphism information content (PIC) and genetic diversity (GD) based on these SNPs showed a higher level of genetic diversity in cultivated type than in wild type. The 415 accessions were clustered into three groups by STRUCTURE software and confirmed using principal component analyses (PCA)-wild type, cultivated type, and admixed wild type. However, unweighted pair group method with arithmetic mean (UPGMA) trees indicated the accessions should be grouped into more clusters. Further analyses identified four groups, the Pure Wild Type, Admixed Wild Type, ancient landraces and modern landraces using STRUCTURE, and the results were confirmed by PCA and UPGMA tree method. A higher level of genetic diversity was detected in ancient landraces and Admixed Wild Type than that in the Pure Wild Type and modern landraces. The highest differentiation was between the Pure Wild Type and modern landraces. A relatively fast LD decay with a short range (kb) was observed, and the LD decays of four inferred populations were different.

CONCLUSIONS

This study is, to our knowledge, the first population genetic analysis of tea germplasm from the Origin Center, Guizhou Plateau, using GBS. The LD pattern, population structure and genetic differentiation of the tea population revealed by our study will benefit further genetic studies, germplasm protection, and breeding.

Genetic Mapping of the Incompatibility Locus in Olive and Development of a Linked Sequence-Tagged Site Marker

The genetic control of self-incompatibility (SI) has been recently disclosed in olive. Inter-varietal crossing confirmed the presence of only two incompatibility groups (G1 and G2), suggesting a simple Mendelian inheritance of the trait. A double digest restriction associated DNA (ddRAD) sequencing of a biparental population segregating for incompatibility groups has been performed and high-density linkage maps were constructed in order to map the SI locus and identify gene candidates and linked markers. The progeny consisted of a full-sib family of 229 individuals derived from the cross 'Leccino' (G1) × 'Dolce Agogia' (G2) varieties, segregating 1:1 (G1:G2), in accordance with a diallelic self-incompatibility (DSI) model. A total of 16,743 single nucleotide polymorphisms was identified, 7,006 in the female parent 'Leccino' and 9,737 in the male parent 'Dolce Agogia.' Each parental map consisted of 23 linkage groups and showed an unusual large size (5,680 cM in 'Leccino' and 3,538 cM in 'Dolce Agogia'). Recombination was decreased across all linkage groups in pollen mother cells of 'Dolce Agogia,' the parent with higher heterozygosity, compared to megaspore mother cells of 'Leccino,' in a context of a species that showed exceptionally high recombination rates. A subset of 109 adult plants was assigned to either incompatibility group by a stigma test and the diallelic self-incompatibility (DSI) locus was mapped to an interval of 5.4 cM on linkage group 18. This region spanned a size of approximately 300 Kb in the olive genome assembly. We developed a sequence-tagged site marker in the DSI locus and identified five haplotypes in 57 cultivars with known incompatibility group assignment. A combination of two single-nucleotide polymorphisms (SNPs) was sufficient to predict G1 or G2 phenotypes in olive cultivars, enabling early marker-assisted selection of compatible genotypes and allowing for a rapid screening of inter-compatibility among cultivars in order to guarantee effective fertilization and increase olive production. The construction of high-density linkage maps has led to the development of the first functional marker in olive and provided positional candidate genes in the SI locus.