Background: Understanding how traditional agriculture systems have been maintained would help design sustainable agriculture. In this study, we examined how farmers have used two types of local trees (Torreya grandis) for stable yield and maintaining genetic diversity in the “globally important agricultural heritage torreya tree system”. The two type of torreya trees are grafted torreya (GT) tree and non-grafted-torreya (NGT) tree. The GT tree has only female and was used to produced seed yields. The NGT tree has both male and female and was used to support GT tree by providing pollens and rootstocks. We first tested the ratio of GT tree to NGT tree, their age groups, ratio of female trees (including GT and NGT trees) to male, and the flowering period of GT and NGT trees. We then tested seed yields and genetic diversity of GT and NGT trees. We further tested gene flow among NGT trees, and the relationship of gene flow with exchange rates of pollens and seeds.Results: GT and NGT trees (male and female) were planted in a mosaic pattern with a ratio of 4:1 (GT:NGT). In this planting pattern, one NGT male trees provided pollen for 20 female trees of GT and NGT. The trees were classified into four age groups (I = 100–400 years old; II = 400–700 years old; III = 700–1000 years old; and IV = 1000–1300 years old) based on basal diameter. The entire flowering period was longer for NGT trees than for GT trees that ensured GT trees (which lack of males) being exposed to pollens. GT tree had high and stable seed yield that increased with age groups. High genetic diversity has been maintained in both rootstocks of the GT trees and NGT trees. There was a strong gene flow among NGT trees, which positive correlated with the exchange rates of pollens and seeds.Conclusions: Our results suggest that farmers obtain stable seed yields, and maintain high genetic diversity by ingeniously using the local GT tree as yield producer and NGT tree as supporter. These GT and NGT trees together ensure sustainable torreya production.

Camellia oleifera is a woody oil crop with the highest oil yield and the largest cultivation area in China, and C. oleifera seed oil is a high-quality edible oil recommended by the Food and Agriculture Organization of the United Nations. The objectives of this study were to investigate the variation in fruit yield traits and seed chemical compositions of wild C. oleifera in China and to identify the differences between wild C. oleifera and cultivated varieties. In this study, we collected wild C. oleifera samples from 13 sites covering the main distribution areas of wild C. oleifera to comprehensively evaluate 25 quantitative traits of wild C. oleifera fruit and seed chemical compositions and collected data of 10 quantitative traits from 434 cultivated varieties for a comparative analysis of the differences between wild and cultivars. The results showed that the coefficients of variation of the 25 quantitative traits of wild C. oleifera ranged from 2.605% to 156.641%, with an average of 38.569%. The phenotypic differentiation coefficients ranged from 25.003% to 99.911%, with an average of 77.894%. The Shannon–Wiener index (H’) ranged from 0.195 to 1.681. Based on the results of principal component analysis (PCA) and phenotypic differentiation coefficients, 10 traits differed significantly between wild C. oleifera and cultivated varieties, while the differentiation coefficients (V_ST) for fresh fruit weight, oleic acid, unsaturated fatty acids, stearic acid, and saturated fatty acids were more than 95%, of which fresh fruit weight and oleic acid content were potential domestication traits of C. oleifera. The results of this study can contribute to the efficient excavation and utilization of wild C. oleifera genetic resources for C. oleifera breeding.

Camellia obtusifolia is a wild relative of the economically important oil crop Camellia oleifera. Here, we sequenced and analyzed the complete chloroplast genome of C. obtusifolia to clarify its phylogenetic position. The cp genome was 156,540 bp in length and exhibited a typical quadripartite structure comprising a LSC (86,254 bp), a SSC (18,400 bp), and two IR (25,943 bp) regions. It contained 134 genes, including 88 protein-coding genes, 37 tRNA genes, and eight rRNA genes, with a GC content of 37.3%. Comparative analysis revealed a 452-bp deletion in the atpH–atpI intergenic region, which served as a distinguishing feature between C. obtusifolia and C. oleifera. Phylogenetic analysis showed that C. obtusifolia clustered into a clade with Camellia sasanqua, Camellia brevistyla, and Camellia meiocarpa, demonstrating the closest relationship with C. sasanqua (bootstrap = 98%), contradicting previous morphological classifications. This study provided essential genomic resources for future phylogenetic and conservation studies in genus Camellia.

Background Low temperature severely limits the growth and development of Camellia oleifera. The mitogen-activated protein kinase (MAPK) cascade plays a key role in the response to cold stress. Methods and results Our study aims to identify MAPK cascade genes in C. oleifera and reveal their roles in response to cold stress. In our study, we systematically identified and analyzed the MAPK cascade gene families of C. oleifera, including their physical and chemical properties, conserved motifs, and multiple sequence alignments. In addition, we characterized the interacting networks of MAPKKK-MAPKK-MAPK in C. oleifera. The molecular mechanism of cold stress resistance of MAPK cascade genes in wild C. oleifera was analyzed by differential gene expression and real-time quantitative reverse transcription-PCR (qRT-PCR). Conclusion In this study, 21 MAPKs, 4 MAPKKs and 55 MAPKKKs genes were identified in the leaf transcriptome of C. oleifera. According to the phylogenetic results, MAPKs were divided into 4 groups (A, B, C and D), MAPKKs were divided into 3 groups (A, B and D), and MAPKKKs were divided into 2 groups (MEKK and Raf). Motif analysis showed that the motifs in each subfamily were conserved, and most of the motifs in the same subfamily were basically the same. The protein interaction network based on Arabidopsis thaliana homologs revealed that MAPK, MAPKK, and MAPKKK genes were widely involved in C. oleifera growth and development and in responses to biotic and abiotic stresses. Gene expression analysis revealed that the CoMAPKKK5 / CoMAPKKK43 / CoMAPKKK49 - CoMAPKK4 - CoMAPK8 module may play a key role in the cold stress resistance of wild C. oleifera at a high-elevation site in Lu Mountain. This study can facilitate the mining and utilization of genetic resources of C. oleifera with low temperature tolerance.