Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history

Vitamin C and suicidal ideation: A cross-sectional and Mendelian randomization study

BACKGROUND

Studies have shown that a history of suicide attempts is associated with low levels of antioxidant vitamins. However, the specific relationship between vitamin C and suicidal ideation remains poorly understood. The purpose of this study is to investigate the effect of vitamin C on the development of suicidal ideation and to explore further the causal relationship using Mendelian randomization (MR).

METHODS

Using the 2005-2018 National Health and Nutrition Examination Survey (NHANES) dataset, a comprehensive, cross-sectional, stratified survey of 28,623 participants aged 20 years and older was conducted. Suicidal ideation was assessed using item 9 of the Patient Health Questionnaire-9 (PHQ-9) item. The relationship between vitamin C and suicidal ideation was revealed by weighted multiple linear regression, smooth curve fitting, hierarchical analysis, and interaction test. To determine causality, MR was performed on vitamin C and suicidal ideation.

RESULTS

In cross-sectional studies, all models showed a negative association between vitamin C and suicidal ideation, with a 37 % reduction in the risk of suicidal ideation for each unit increase in vitamin C in the highest quartile in the fully adjusted model. The results of the subgroup analysis showed that the association between vitamin C and suicidal ideation was stable across gender, age, education, smoking status, and marital status. In the MR study, our bivariate model showed no significant causal relationship between vitamin C and the development of suicidal ideation.

CONCLUSIONS

Dietary vitamin C intake was negatively associated with the development of suicidal ideation. However, genetic evidence does not support causation of these associations.

An ultra-dense linkage map identified quantitative trait loci corresponding to fruit quality- and size-related traits in red goji berry

Goji berries are a small-fruited shrub with industrial importance whose fruit considered beneficial in both fresh and dried forms. Current germplasms of goji berries include small fruits with a short shelf life, less sweet and bitter taste, and a lack of appropriate genetic information. This study aimed to employ whole genome resequencing to generate an ultra-dense bin linkage map and to elucidate the genetic basis of goji fruit quality and size using quantitative trait loci (QTL) mapping analysis in a cross-pollinated hybrid population. To achieve this goal, human sensory tests were carried out to determine the bitter taste (BT) and sweet taste (ST), and to quantify the soluble solid content (SSC), fruit firmness (FF), and fruit size-related traits of fresh goji fruits over three or four years. The results revealed that the goji bin linkage map based on resequencing spanned a total length of 966.42 cM and an average bin interval of 0.03 cM. Subsequent variant calling and ordering resulted in 3,058 bins containing 35,331 polymorphic markers across 12 chromosomes. A total of 99 QTLs, with individual loci in different environments explaining a phenotypic variance of 1.21-16.95% were identified for the studied traits. Ten major effects, including colocalized QTLs corresponding to different traits, were identified on chromosomes 1, 3, 5, 6, 7, and 8, with a maximum Logarithm of Odds (LOD) of 29.25 and 16.95% of explained phenotypic variance (PVE). In addition, four stable loci, one for FF, one for fruit weight (FW), and two for fruit shape index (FSI), were mainly mapped on chromosomes 5, 6, and 7, elucidating 2.10-16.95% PVE. These findings offer valuable insights into the genetic architecture of goji fruit traits along with identified specific loci and markers to further improve and develop sweeter, less bitter and larger fruited goji berry cultivars with extended shelf life.

Potential drug targets for osteoporosis identified: A Mendelian randomization study

Background

Osteoporosis is a prevalent global health condition, primarily affecting the aging population, and several therapies for osteoporosis have been widely used. However, available drugs for osteoporosis are far from satisfactory because they cannot alleviate disease progression. This study aimed to explore potential drug targets for osteoporosis through Mendelian randomization analysis.

Methods

Using cis-expression quantitative trait loci (cis-eQTL) data of druggable genes and two genome-wide association studies (GWAS) datasets related to osteoporosis (UK Biobank and FinnGen cohorts), we employed mendelian randomization (MR) analysis to identify the druggable genes with causal relationships with osteoporosis. Subsequently, a series of follow-up analyses were conducted, such as colocalization analysis, cell-type specificity analysis, and correlation analysis with risk factors. The association between potential drug targets and osteoporosis was validated by qRT-PCR.

Results

Six druggable genes with causal relationships with osteoporosis were identified and successfully replicated, including ACPP, DNASE1L3, IL32, PPOX, ST6GAL1, and TGM3. Cell-type specificity analysis revealed that PPOX and ST6GAL1 were expressed in all cell types in the bone samples, while IL32, ACPP, DNASE1L3, and TGM3 were expressed in specific cell types. The GWAS data showed there were seven risk factors for osteoporosis, including vitamin D deficiency, COPD, physical activity, BMI, MMP-9, ALP and PTH. Furthermore, ACPP was associated with vitamin D deficiency and COPD; DNASE1L3 was linked to physical activity; IL32 correlated with BMI and MMP-9; and ST6GAL1 was related to ALP, physical activity, and MMP-9. Among these risk factors, only MMP-9 had a high genetic correlation with osteoporosis. The results of qRT-PCR demonstrated that IL32 was upregulated while ST6GAL1 was downregulated in peripheral blood of osteoporosis patients.

Conclusion

Our findings suggested that those six druggable genes offer potential drug targets for osteoporosis and require further clinical investigation, especially IL32 and ST6GAL1.

Multi-omics analysis unravels chemical roadmap and genetic basis for peach fruit aroma improvement

Selection of fruits with enhanced health benefits and superior flavor is an important aspect of peach breeding. Understanding the genetic interplay between appearance and flavor chemicals remains a major challenge. We identify the most important volatiles contributing to consumer preferences for peach, thus establishing priorities for improving flavor quality. We quantify volatiles of a peach population consisting of 184 accessions and demonstrate major reductions in the important flavor volatiles linalool and Z-3-hexenyl acetate in red-fleshed accessions. We identify 474 functional gene regulatory networks (GRNs), among which GRN05 plays a crucial role in controlling both red flesh and volatile content through the NAM/ATAF1/2/CUC (NAC) transcription factor PpBL. Overexpressing PpBL results in reduced expression of PpNAC1, a positive regulator for Z-3-hexenyl acetate and linalool synthesis. Additionally, we identify haplotypes for three tandem PpAATs that are significantly correlated with reduced gene expression and ester content. We develop genetic resources for improvement of fruit quality.

Genomic variation, environmental adaptation, and feralization in ramie, an ancient fiber crop

Feralization is an important evolutionary process, but the mechanisms behind it remain poorly understood. Here, we use the ancient fiber crop ramie (Boehmeria nivea (L.) Gaudich.) as a model to investigate genomic changes associated with both domestication and feralization. We first produced a chromosome-scale de novo genome assembly of feral ramie and investigated structural variations between feral and domesticated ramie genomes. Next, we gathered 915 accessions from 23 countries, comprising cultivars, major landraces, feral populations, and the wild progenitor. Based on whole-genome resequencing of these accessions, we constructed the most comprehensive ramie genomic variation map to date. Phylogenetic, demographic, and admixture signal detection analyses indicated that feral ramie is of exoferal or exo-endo origin, i.e., descended from hybridization between domesticated ramie and the wild progenitor or ancient landraces. Feral ramie has higher genetic diversity than wild or domesticated ramie, and genomic regions affected by natural selection during feralization differ from those under selection during domestication. Ecological analyses showed that feral and domesticated ramie have similar ecological niches that differ substantially from the niche of the wild progenitor, and three environmental variables are associated with habitat-specific adaptation in feral ramie. These findings advance our understanding of feralization, providing a scientific basis for the excavation of new crop germplasm resources and offering novel insights into the evolution of feralization in nature.

Genome-wide association study of the loci and candidate genes associated with agronomic traits in Amomum villosum Lour

Amomum villosum Lour. (A. villosum) is a valuable herbaceous plant that produces the famous traditional Chinese medicine Amori Fructus. Identifying molecular markers associated with the growth of A. villosum can facilitate molecular marker-assisted breeding of the plant. This study employed 75 A. villosum accessions as the test material and utilized 71 pairs of polymorphic simple sequence repeat (SSR) molecular markers to genotype the population. The study analyzed the association between SSR markers and phenotypic traits through the linkage imbalance and population structure analysis. Candidate genes associated with the molecular markers were also identified. The results showed that the phenotypic diversity index range of the 12 agronomic traits was 4.081-4.312 and conformed to a normal distribution. Moreover, 293 allelic variations were detected in the 75 accessions, with an average of 5.32 amplified alleles per loci, ranging from 3 to 8. The maximum number of amplified alleles for AVL12 was 8. The population structure and cluster analysis indicated that the accessions could be divided into two subgroups. Using the mixed linear model (MLM) model of population structure (Q)+kinship matrix (K) for association analysis, three SSR molecular markers significantly associated with the agronomic traits were detected. Fluorescence quantification was used to analyze the expression levels of six candidate genes, and it was found that three of the genes were differentially expressed in phenotypically different accessions. This study is the first to use SSR markers for genome-wide association study (GWAS) mapping and identification of the associated agronomic traits in A. villosum. The results of this study provide a basis for identifying genetic markers for growth traits for marker-assisted breeding in A. villosum.

Origins and Genetic Characteristics of Egyptian Peach

Peach (Prunus persica), a significant economic fruit tree in the Rosaceae family, is extensively cultivated in temperate and subtropical regions due to its abundant genetic diversity, robust adaptability, and high nutritional value. Originating from China over 4000 years ago, peaches were introduced to Persia through the Silk Road during the Han Dynasty and gradually spread to India, Greece, Rome, Egypt, Europe, and America. Currently grown in more than 80 countries worldwide, the expansion of peach cultivation in Egypt is mainly due to the development and utilization of peach varieties with low chilling requirements. These varieties exhibit unique phenotypic characteristics such as early maturity, reduced need for winter cold temperatures, low water requirements, and high economic value. In this study, a systematic analysis was conducted on the genetic characteristics and kinship relationships of peaches with low chilling requirements in Egypt. We conducted a comprehensive evolutionary and Identity-by-Descent (IBD) analysis on over 300 peach core germplasm resources, including Egyptian cultivars with low chilling requirements, to investigate their origin and genetic characteristics. The evolutionary analysis revealed that 'Bitter almond' is closely related to China's wild relative species Prunus tangutica Batal, while 'Early grand' shares one branch with Chinese ornamental peach cultivars, and 'Nemaguard' clusters with some ancient local varieties from China. The IBD analysis also indicated similar genetic backgrounds, suggesting a plausible origin from China. Similarly, the analysis suggested that 'Swelling' may have originated from the Czech Republic while 'Met ghamr' has connections to South Africa. 'Desert red', 'Early swelling', and 'Florida prince' are likely derived from Brazil. These findings provide valuable insights into the genetic characteristics of Egyptian peach cultivars. They offer a significant foundation for investigating the origin and spread of cultivated peaches worldwide and serve as a valuable genetic resource for breeding low chilling requirement cultivars, which is of considerable significance for the advancement of peach cultivation in Egypt.

PPGV: a comprehensive database of peach population genome variation

Peach tree is one of the most important fruit trees in the world, and it has been cultivated for more than 7,500 years. In recent years, the genome and population resequencing of peach trees have been published continuously, which has effectively promoted the research of peach tree genetics and breeding. In order to promote the further mining and utilization of these data, we integrated and constructed a comprehensive peach genome and variation database (PPGV, http://peachtree.work/home ). The PPGV contains 10 sets of published peach tree genome data, as well as genomic variation information for 1,378 peach tree samples (the resequencing data of 1,378 samples were aligned with the high-quality genomes of Lovell, CN14 and Chinesecling, respectively, for mutation detection). A variety of useful and flexible tools, such as BLAST, Gene ID Convert, KEGG/GO Enrichment, Primer Design and Gene function, were also specially designed for searching data and assisting in breeding.

Genomic diversity, population structure, and genome-wide association reveal genetic differentiation and trait improvements in mango

Mango (Mangifera indica L.) has been widely cultivated as a culturally and economically significant fruit tree for roughly 4000 years. Despite its rich history, little is known about the crop's domestication, genomic variation, and the genetic loci underlying agronomic traits. This study employs the whole-genome re-sequencing of 224 mango accessions sourced from 22 countries, with an average sequencing depth of 16.37×, to explore their genomic variation and diversity. Through phylogenomic analysis, M. himalis J.Y. Liang, a species grown in China, was reclassified into the cultivated mango group known as M. indica. Moreover, our investigation of mango population structure and differentiation revealed that Chinese accessions could be divided into two distinct gene pools, indicating the presence of independent genetic diversity ecotypes. By coupling genome-wide association studies with analyses of genotype variation patterns and expression patterns, we identified several candidate loci and dominant genotypes associated with mango flowering capability, fruit weight, and volatile compound production. In conclusion, our study offers valuable insights into the genetic differentiation of mango populations, paving the way for future agronomic improvements through genomic-assisted breeding.

Comparative population genomics reveals convergent and divergent selection in the apricot-peach-plum-mei complex

The economically significant genus Prunus includes fruit and nut crops that have been domesticated for shared and specific agronomic traits; however, the genomic signals of convergent and divergent selection have not been elucidated. In this study, we aimed to detect genomic signatures of convergent and divergent selection by conducting comparative population genomic analyses of the apricot-peach-plum-mei (APPM) complex, utilizing a haplotype-resolved telomere-to-telomere (T2T) genome assembly and population resequencing data. The haplotype-resolved T2T reference genome for the plum cultivar was assembled through HiFi and Hi-C reads, resulting in two haplotypes 251.25 and 251.29 Mb in size, respectively. Comparative genomics reveals a chromosomal translocation of ~1.17 Mb in the apricot genomes compared with peach, plum, and mei. Notably, the translocation involves the D locus, significantly impacting titratable acidity (TA), pH, and sugar content. Population genetic analysis detected substantial gene flow between plum and apricot, with introgression regions enriched in post-embryonic development and pollen germination processes. Comparative population genetic analyses revealed convergent selection for stress tolerance, flower development, and fruit ripening, along with divergent selection shaping specific crop, such as somatic embryogenesis in plum, pollen germination in mei, and hormone regulation in peach. Notably, selective sweeps on chromosome 7 coincide with a chromosomal collinearity from the comparative genomics, impacting key fruit-softening genes such as PG, regulated by ERF and RMA1H1. Overall, this study provides insights into the genetic diversity, evolutionary history, and domestication of the APPM complex, offering valuable implications for genetic studies and breeding programs of Prunus crops.

Metabolite-based genome-wide association studies enable the dissection of the genetic bases of flavonoids, betaine and spermidine in wolfberry (Lycium)

Wolfberry is a plant with medicinal and food values. However, its bioactive ingredients and the corresponding genetic bases have not been determined. Here, we de novo generated a chromosome-level genome assembly for wolfberry, yielding a genome sequence of ~1.77 Gb with contig N50 of 50.55 Mb and 39 224 predicted gene models. A variation map, using 307 re-sequenced accessions, was called based on this genome assembly. Furthermore, the fruit metabolome of these accessions was profiled using 563 annotated metabolites, which separated Lycium barbarum L. and non-L. barbarum L. The flavonoids, coumarins, alkaloids and nicotinic acid contents were higher in the former than in the latter. A metabolite-based genome-wide association study mapped 156 164 significant single nucleotide polymorphisms corresponding to 340 metabolites. This included 19 219 unique lead single nucleotide polymorphisms in 1517 significant association loci, of which three metabolites, flavonoids, betaine and spermidine, were highlighted. Two candidate genes, LbUGT (evm.TU.chr07.2692) and LbCHS (evm.TU.chr07.2738), with non-synonymous mutations, were associated with the flavonoids content. LbCHS is a structural gene that interacts with a nearby MYB transcription factor (evm.TU.chr07.2726) both in L. barbarum and L. ruthenicum. Thus, these three genes might be involved in the biosynthesis/metabolism of flavonoids. LbSSADH (evm.TU.chr09.627) was identified as possibly participating in betaine biosynthesis/metabolism. Four lycibarbarspermidines (E-G and O) were identified, and only the lycibarbarspermidines O content was higher in L. barbarum varieties than in non-L. barbarum varieties. The evm.TU.chr07.2680 gene associated with lycibarbarspermidines O was annotated as an acetyl-CoA-benzylalcohol acetyltransferase, suggesting that it is a candidate gene for spermidine biosynthesis. These results provide novel insights into the specific metabolite profile of non-L. barbarum L. and the genetic bases of flavonoids, betaine and spermidine biosynthesis/metabolism.

Natural selection and genetic diversity maintenance in a parasitic wasp during continuous biological control application

Aphidius gifuensis is a parasitoid wasp and primary endoparasitoid enemy of the peach potato aphid, Myzus persicae. Artificially reared, captive wasps of this species have been extensively and effectively used to control populations of aphids and limit crop loss. However, the consequences of large-scale releasing of captive A. gifuensis, such as genetic erosion and reduced fitness in wild populations of this species, remains unclear. Here, we sequence the genomes of 542 A. gifuensis individuals collected across China, including 265 wild and 277 human-intervened samples. Population genetic analyses on wild individuals recovered Yunnan populations as the ancestral group with the most complex genetic structure. We also find genetic signature of environmental adaptation during the dispersal of wild populations from Yunnan to other regions. While comparative genomic analyses of captive wasps revealed a decrease in genetic diversity during long-term rearing, population genomic analyses revealed signatures of natural selection by several biotic (host plants) or abiotic (climate) factors, which support maintenance of the gene pool of wild populations in spite of the introduction of captive wasps. Therefore, the impact of large-scale release is reduced. Our study suggests that A. gifuensis is a good system for exploring the genetic and evolutionary effects of mass rearing and release on species commonly used as biocontrol agents.

Combined ATAC-seq, RNA-seq, and GWAS analysis reveals glycogen metabolism regulatory network in Jinjiang oyster ( Crassostrea ariakensis)

Glycogen serves as the principal energy reserve for metabolic processes in aquatic shellfish and substantially contributes to the flavor and quality of oysters. The Jinjiang oyster ( Crassostrea ariakensis) is an economically and ecologically important species in China. In the present study, RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) were performed to investigate gene expression and chromatin accessibility variations in oysters with different glycogen contents. Analysis identified 9 483 differentially expressed genes (DEGs) and 7 215 genes with significantly differential chromatin accessibility (DCAGs) were obtained, with an overlap of 2 600 genes between them. Notably, a significant proportion of these genes were enriched in pathways related to glycogen metabolism, including "Glycogen metabolic process" and "Starch and sucrose metabolism". In addition, genome-wide association study (GWAS) identified 526 single nucleotide polymorphism (SNP) loci associated with glycogen content. These loci corresponded to 241 genes, 63 of which were categorized as both DEGs and DCAGs. This study enriches basic research data and provides insights into the molecular mechanisms underlying the regulation of glycogen metabolism in C. ariakensis.

Genome variation and LTR-RT analyses of an ancient peach landrace reveal mechanism of blood-flesh fruit color formation and fruit maturity date advancement

Peach (Prunus persica) landrace has typical regional characteristics, strong environmental adaptability, and contains many valuable genes that provide the foundation for breeding excellent varieties. Therefore, it is necessary to assemble the genomes of specific landraces to facilitate the localization and utilization of these genes. Here, we de novo assembled a high-quality genome from an ancient blood-fleshed Chinese landrace Tianjin ShuiMi (TJSM) that originated from the China North Plain. The assembled genome size was 243.5 Mb with a contig N50 of 23.7 Mb and a scaffold N50 of 28.6 Mb. Compared with the reported peach genomes, our assembled TJSM genome had the largest number of specific structural variants (SVs) and long terminal repeat-retrotransposons (LTR-RTs). Among the LTR-RTs with the potential to regulate their host genes, we identified a 6688 bp LTR-RT (named it blood TE) in the promoter of NAC transcription factor-encoding PpBL, a gene regulating peach blood-flesh formation. The blood TE was not only co-separated with the blood-flesh phenotype but also associated with fruit maturity date advancement and different intensities of blood-flesh color formation. Our findings provide new insights into the mechanism underlying the development of the blood-flesh color and determination of fruit maturity date and highlight the potential of the TJSM genome to mine more variations related to agronomic traits in peach fruit.

Soluble sugar and organic acid composition and flavor evaluation of Chinese cherry fruits

Chinese cherry is an economically important fruit crop native to China. Flavor quality is greatly influenced by compositions of soluble sugars and organic acids. To better understand the flavor quality of Chinese cherry, we determined sugar and acid components in thirty-eight landrace and cultivar collections, and two wild resources using the HPLC method. Glucose and fructose were the main components, accounting for 85.91% of soluble sugars. Malic acid was the predominant organic acid, with an average proportion of 65.73% of total acids. Correlation and PCA analysis revealed seven key indicators for evaluating fruit flavor. Compared with wild Chinese cherry, the cultivated collections exhibited higher levels of soluble sugars, especially fructose, and lower levels of organic acid, particularly malic acid in fruits. Finally, we have established grading criteria for seven flavor indicators in Chinese cherry. Our study provides valuable references for identifying flavor compounds and improving flavor quality of Chinese cherry.

Genome-wide association study on resistance of cultivated soybean to Fusarium oxysporum root rot in Northeast China

BACKGROUND

Fusarium oxysporum is a prevalent fungal pathogen that diminishes soybean yield through seedling disease and root rot. Preventing Fusarium oxysporum root rot (FORR) damage entails on the identification of resistance genes and developing resistant cultivars. Therefore, conducting fine mapping and marker development for FORR resistance genes is of great significance for fostering the cultivation of resistant varieties. In this study, 350 soybean germplasm accessions, mainly from Northeast China, underwent genotyping using the SoySNP50K Illumina BeadChip, which includes 52,041 single nucleotide polymorphisms (SNPs). Their resistance to FORR was assessed in a greenhouse. Genome-wide association studies utilizing the general linear model, mixed linear model, compressed mixed linear model, and settlement of MLM under progressively exclusive relationship models were conducted to identify marker-trait associations while effectively controlling for population structure.

RESULTS

The results demonstrated that these models effectively managed population structure. Eight SNP loci significantly associated with FORR resistance in soybean were detected, primarily located on Chromosome 6. Notably, there was a strong linkage disequilibrium between the large-effect SNPs ss715595462 and ss715595463, contributing substantially to phenotypic variation. Within the genetic interval encompassing these loci, 28 genes were present, with one gene Glyma.06G088400 encoding a protein kinase family protein containing a leucine-rich repeat domain identified as a potential candidate gene in the reference genome of Williams82. Additionally, quantitative real-time reverse transcription polymerase chain reaction analysis evaluated the gene expression levels between highly resistant and susceptible accessions, focusing on primary root tissues collected at different time points after F. oxysporum inoculation. Among the examined genes, only this gene emerged as the strongest candidate associated with FORR resistance.

CONCLUSIONS

The identification of this candidate gene Glyma.06G088400 improves our understanding of soybean resistance to FORR and the markers strongly linked to resistance can be beneficial for molecular marker-assisted selection in breeding resistant soybean accessions against F. oxysporum.

Genome-wide association study of powdery mildew resistance in cultivated soybean from Northeast China

Powdery mildew (PMD), caused by the pathogen Microsphaera diffusa, leads to substantial yield decreases in susceptible soybean under favorable environmental conditions. Effective prevention of soybean PMD damage can be achieved by identifying resistance genes and developing resistant cultivars. In this study, we genotyped 331 soybean germplasm accessions, primarily from Northeast China, using the SoySNP50K BeadChip, and evaluated their resistance to PMD in a greenhouse setting. To identify marker-trait associations while effectively controlling for population structure, we conducted genome-wide association studies utilizing factored spectrally transformed linear mixed models, mixed linear models, efficient mixed-model association eXpedited, and compressed mixed linear models. The results revealed seven single nucleotide polymorphism (SNP) loci strongly associated with PMD resistance in soybean. Among these, one SNP was localized on chromosome (Chr) 14, and six SNPs with low linkage disequilibrium were localized near or in the region of previously mapped genes on Chr 16. In the reference genome of Williams82, we discovered 96 genes within the candidate region, including 17 resistance (R)-like genes, which were identified as potential candidate genes for PMD resistance. In addition, we performed quantitative real-time reverse transcription polymerase chain reaction analysis to evaluate the gene expression levels in highly resistant and susceptible genotypes, focusing on leaf tissues collected at different times after M. diffusa inoculation. Among the examined genes, three R-like genes, including Glyma.16G210800, Glyma.16G212300, and Glyma.16G213900, were identified as strong candidates associated with PMD resistance. This discovery can significantly enhance our understanding of soybean resistance to PMD. Furthermore, the significant SNPs strongly associated with resistance can serve as valuable markers for genetic improvement in breeding M. diffusa-resistant soybean cultivars.

Genome-wide presence/absence variation discovery and its application in Peach (Prunus persica)

Gene presence/absence variation (PAV) is an important contributor to the studies of genetic diversity, gene identification, and molecular marker development in plants. In the present study, 100 peach (Prunus persica) accessions were used for genome resequencing to identify PAVs. Alignmentwith a reference genome yielded a total of 2.52 Mb non-reference sequences and 923 novel genes were identified. The dispensable PAVs were enriched in resistance, perhaps reflecting their roles in plant adaptation to various environments. Furthermore, selection sweeps associated with peach domestication and improvement were identified based on PAV data. Only 4.3% and 13.4% of domestication and improvement sweeps, respectively, were identified simultaneously using single nucleotide polymorphism (SNP) data, suggesting flexible identification between the different methods. To further verify the applicability of PAV identification, a genome-wide association study was conducted using 21 agronomic traits. Some of the identified loci were consistent with those reported in previous studies, while some were mapped for the first time; the latter included petiole length, petiole gland shape, and petiole gland number. Through tissue-specific expression analysis and gene transformation experiments, a novel gene, evm.model.Contig322_A94.1, was identified and found to be involved in chilling requirements. We speculated that this novel gene might regulate the trait by participating in the ABA signaling pathway. The PAVs identified in P. persica provide valuable resources for mapping the entire gene set and identifying optional markers for molecular selection in future studies.

Identifying underlying patterns in Alzheimer's disease trajectory: a deep learning approach and Mendelian randomization analysis

Background

Alzheimer's disease (AD) is a heterogeneously progressive neurodegeneration disorder with varied rates of deterioration, either between subjects or within different stages of a certain subject. Estimating the course of AD at early stages has treatment implications. We aimed to analyze disease progression to identify distinct patterns in AD trajectory.

Methods

We proposed a deep learning model to identify underlying patterns in the trajectory from cognitively normal (CN) to a state of mild cognitive impairment (MCI) to AD dementia, by jointly predicting time-to-conversion and clustering out distinct subgroups characterized by comprehensive features as well as varied progression rates. We designed and validated our model on the ADNI dataset (1370 participants). Prediction of time-to-conversion in AD trajectory was used to validate the expression of the identified patterns. Causality between patterns and time-to-conversion was further inferred using Mendelian randomization (MR) analysis. External validation was performed on the AIBL dataset (233 participants).

Findings

The proposed model clustered out patterns characterized by significantly different biomarkers and varied progression rates. The discovered patterns also showed a strong prediction ability, as indicated by hazard ratio (CN→MCI, HR = 3.51, p < 0.001; MCI→AD, HR = 8.11, p < 0.001), C-Index (CN→MCI, 0.618; MCI→AD, 0.718), and AUC (CN→MCI, 3 years 0.802, 5 years 0.876; MCI→AD, 3 years 0.914, 5 years 0.957). In the external validation cohort, our model demonstrated competitive performance on conversion time prediction (CN→MCI, C-Index = 0.693; MCI→AD, C-Index = 0.752). Moreover, suggestive associations between CN→MCI/MCI→AD patterns with four/three SNPs were mediated and MR analysis indicated a causal link between MCI→AD patterns and time-to-conversion in the first three years.

Interpretation

Our proposed model identifies biologically and clinically meaningful patterns from real-world data and provides promising performance on time-to-conversion prediction in AD trajectory, which could promote the understanding of disease progression, facilitate clinical trial design, and provide potential for decision-making.

Funding

The National Key Research and Development Program of China, the Key R&D Program of Zhejiang, and the National Nature Science Foundation of China.

MADS-box protein PpDAM6 regulates chilling requirement-mediated dormancy and bud break in peach

Bud dormancy is crucial for winter survival and is characterized by the inability of the bud meristem to respond to growth-promotive signals before the chilling requirement (CR) is met. However, our understanding of the genetic mechanism regulating CR and bud dormancy remains limited. This study identified PpDAM6 (DORMANCY-ASSOCIATED MADS-box) as a key gene for CR using a genome-wide association study analysis based on structural variations in 345 peach (Prunus persica (L.) Batsch) accessions. The function of PpDAM6 in CR regulation was demonstrated by transiently silencing the gene in peach buds and stably overexpressing the gene in transgenic apple (Malus × domestica) plants. The results showed an evolutionarily conserved function of PpDAM6 in regulating bud dormancy release, followed by vegetative growth and flowering, in peach and apple. The 30-bp deletion in the PpDAM6 promoter was substantially associated with reducing PpDAM6 expression in low-CR accessions. A PCR marker based on the 30-bp indel was developed to distinguish peach plants with non-low and low CR. Modification of the H3K27me3 marker at the PpDAM6 locus showed no apparent change across the dormancy process in low- and non-low- CR cultivars. Additionally, H3K27me3 modification occurred earlier in low-CR cultivars on a genome-wide scale. PpDAM6 could mediate cell-cell communication by inducing the expression of the downstream genes PpNCED1 (9-cis-epoxycarotenoid dioxygenase 1), encoding a key enzyme for ABA biosynthesis, and CALS (CALLOSE SYNTHASE), encoding callose synthase. We shed light on a gene regulatory network formed by PpDAM6-containing complexes that mediate CR underlying dormancy and bud break in peach. A better understanding of the genetic basis for natural variations of CR can help breeders develop cultivars with different CR for growing in different geographical regions.

Genome-wide association study of maize resistance to Pythium aristosporum stalk rot

Stalk rot, a severe and widespread soil-borne disease in maize, globally reduces yield and quality. Recent documentation reveals that Pythium aristosporum has emerged as one of the dominant causal agents of maize stalk rot. However, a previous study of maize stalk rot disease resistance mechanisms and breeding had mainly focused on other pathogens, neglecting P. aristosporum. To mitigate crop loss, resistance breeding is the most economical and effective strategy against this disease. This study involved characterizing resistance in 295 inbred lines using the drilling inoculation method and genotyping them via sequencing. By combining with population structure, disease resistance phenotype, and genome-wide association study (GWAS), we identified 39 significant single-nucleotide polymorphisms (SNPs) associated with P. aristosporum stalk rot resistance by utilizing six statistical methods. Bioinformatics analysis of these SNPs revealed 69 potential resistance genes, among which Zm00001d051313 was finally evaluated for its roles in host defense response to P. aristosporum infection. Through virus-induced gene silencing (VIGS) verification and physiological index determination, we found that transient silencing of Zm00001d051313 promoted P. aristosporum infection, indicating a positive regulatory role of this gene in maize's antifungal defense mechanism. Therefore, these findings will help advance our current understanding of the underlying mechanisms of maize defense to Pythium stalk rot.

Plant biochemical genetics in the multiomics era

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel's milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

Plant genome resequencing and population genomics: Current status and future prospects

Advances in DNA sequencing technology have sparked a genomics revolution, driving breakthroughs in plant genetics and crop breeding. Recently, the focus has shifted from cataloging genetic diversity in plants to exploring their functional significance and delivering beneficial alleles for crop improvement. This transformation has been facilitated by the increasing adoption of whole-genome resequencing. In this review, we summarize the current progress of population-based genome resequencing studies and how these studies affect crop breeding. A total of 187 land plants from 163 countries have been resequenced, comprising 54 413 accessions. As part of resequencing efforts 367 traits have been surveyed and 86 genome-wide association studies have been conducted. Economically important crops, particularly cereals, vegetables, and legumes, have dominated the resequencing efforts, leaving a gap in 49 orders, including Lycopodiales, Liliales, Acorales, Austrobaileyales, and Commelinales. The resequenced germplasm is distributed across diverse geographic locations, providing a global perspective on plant genomics. We highlight genes that have been selected during domestication, or associated with agronomic traits, and form a repository of candidate genes for future research and application. Despite the opportunities for cross-species comparative genomics, many population genomic datasets are not accessible, impeding secondary analyses. We call for a more open and collaborative approach to population genomics that promotes data sharing and encourages contribution-based credit policy. The number of plant genome resequencing studies will continue to rise with the decreasing DNA sequencing costs, coupled with advances in analysis and computational technologies. This expansion, in terms of both scale and quality, holds promise for deeper insights into plant trait genetics and breeding design.

Genetic diversity and genome-wide association study of 13 agronomic traits in 977 Beta vulgaris L. germplasms

BACKGROUND

Sugar beet (Beta vulgaris L.) is an economically essential sugar crop worldwide. Its agronomic traits are highly diverse and phenotypically plastic, influencing taproot yield and quality. The National Beet Medium-term Gene Bank in China maintains more than 1700 beet germplasms with diverse countries of origin. However, it lacks detailed genetic background associated with morphological variability and diversity.

RESULTS

Here, a comprehensive genome-wide association study (GWAS) of 13 agronomic traits was conducted in a panel of 977 sugar beet accessions. Almost all phenotypic traits exhibited wide genetic diversity and high coefficient of variation (CV). A total of 170,750 high-quality single-nucleotide polymorphisms (SNPs) were obtained using the genotyping-by-sequencing (GBS). Neighbour-joining phylogenetic analysis, principal component analysis, population structure and kinship showed no obvious relationships among these genotypes based on subgroups or regional sources. GWAS was carried out using a mixed linear model, and 159 significant associations were detected for these traits. Within the 25 kb linkage disequilibrium decay of the associated markers, NRT1/PTR FAMILY 6.3 (BVRB_5g097760); nudix hydrolase 15 (BVRB_8g182070) and TRANSPORT INHIBITOR RESPONSE 1 (BVRB_8g181550); transcription factor MYB77 (BVRB_2g023500); and ethylene-responsive transcription factor ERF014 (BVRB_1g000090) were predicted to be strongly associated with the taproot traits of root groove depth (RGD); root shape (RS); crown size (CS); and flesh colour (FC), respectively. For the aboveground traits, UDP-glycosyltransferase 79B6 (BVRB_9g223780) and NAC domain-containing protein 7 (BVRB_5g097990); F-box protein At1g10780 (BVRB_6g140760); phosphate transporter PHO1 (BVRB_3g048660); F-box protein CPR1 (BVRB_8g181140); and transcription factor MYB77 (BVRB_2g023500) and alcohol acyltransferase 9 (BVRB_2g023460) might be associated with the hypocotyl colour (HC); plant type (PT); petiole length (PL); cotyledon size (C); and fascicled leaf type (FLT) of sugar beet, respectively. AP-2 complex subunit mu (BVRB_5g106130), trihelix transcription factor ASIL2 (BVRB_2g041790) and late embryogenesis abundant protein 18 (BVRB_5g106150) might be involved in pollen quantity (PQ) variation. The candidate genes extensively participated in hormone response, nitrogen and phosphorus transportation, secondary metabolism, fertilization and embryo maturation.

CONCLUSIONS

The genetic basis of agronomical traits is complicated in heterozygous diploid sugar beet. The putative valuable genes found in this study will help further elucidate the molecular mechanism of each phenotypic trait for beet breeding.

Multiple-statistical genome-wide association analysis and genomic prediction of fruit aroma and agronomic traits in peaches

'Chinese Cling' is an important founder in peach breeding history due to the pleasant flavor. Genome-wide association studies (GWAS) combined with genomic selection are promising tools in fruit tree breeding, as there is a considerable time lapse between crossing and release of a cultivar. In this study, 242 peaches from Shanghai germplasm were genotyped with 145 456 single-nucleotide polymorphisms (SNPs). The six agronomic traits of fruit flesh color, fruit shape, fruit hairiness, flower type, pollen sterility, and soluble solids content, along with 14 key volatile odor compounds (VOCs), were recorded for multiple-statistical GWAS. Except the reported candidate genes, six novel genes were identified as associated with these traits. Thirty-nine significant SNPs were associated with eight VOCs. The putative candidate genes were confirmed for VOCs by RNA-seq, including three genes in the biosynthesis pathway found to be associated with linalool, soluble solids content, and cis-3-hexenyl acetate. Multiple-trait genomic prediction enhanced the predictive ability for γ-decalactone to 0.7415 compared with the single-trait model value of 0.1017. One PTS1-SSR marker was designed to predict the linalool content, and the favorable genotype 187/187 was confirmed, mainly existing in the 'Shanghai Shuimi' landrace. Overall, our findings will be helpful in determining peach accessions with the ideal phenotype and show the potential of multiple-trait genomic prediction to improve accuracy for highly correlated genetic traits. The diagnostic marker will be valuable for the breeder to bridge the gap between quantitative trait loci and marker-assisted selection for developing strong-aroma cultivars.

Genome-Wide Association Analysis of Muscle pH in Texel Sheep × Altay Sheep F2 Resource Population

pH was one of the important meat quality traits, which was an important factor affecting the storage/shelf life and quality of meat in meat production. In order to find a way to extend the storage/shelf life, the pH values (pH_45min, pH_24h, pH_48h and pH_72h) of the longissimus dorsi muscles in F₂ individuals of 462 Texel sheep × Altay sheep were determined, genotyping was performed using Illumina Ovine SNP 600 K BeadChip and whole genome resequencing technology, a genome-wide association analysis (GWAS) was used to screen the candidate genes and molecular markers for pH values related to the quality traits of mutton, and the effects of population stratification were detected by Q-Q plots. The results showed that the pH population stratification analysis did not find significant systemic bias, and there was no obvious population stratification effect. The results of the association analysis showed that 28 SNPs significantly associated with pH reached the level of genomic significance. The candidate gene associated with pH_45min was identified as the CCDC92 gene by gene annotation and a search of the literature. Candidate genes related to pH_24h were KDM4C, TGFB2 and GOT2 genes. The candidate genes related to pH_48h were MMP12 and MMP13 genes. The candidate genes related to pH_72h were HILPDA and FAT1 genes. Further bioinformatics analyses showed 24 gene ontology terms and five signaling pathways that were significantly enriched (p ≤ 0.05). Many terms and pathways were related to cellular components, processes of protein modification, the activity of protein dimerization and hydrolase activity. These identified SNPs and genes could provide useful information about meat and the storage/shelf life of meat, thereby extending the storage/shelf life and quality of meat.

Genotypic variation of tocopherol content in a representative genetic population and genome-wide association study on tocopherol in rapeseed (Brassica napus)

Tocopherols (Tocs) are a kind of lipid-soluble substance required for the normal physiological function of mammals, particularly their antioxidant capacity. Rapeseed (Brassica napus) oil is an important source of exogenous Tocs. However, the genotypic differences in the total Toc contents, the Toc composition in the seeds, and the molecular markers associated with the seed Toc remain largely unknown. Here, we selected 290 rapeseed accessions based on the resequencing of 991 genomes in a worldwide collection of rapeseed germplasm. The contents of the four Toc isoforms, namely, α-, β-, γ-, and δ-Tocs, were also measured. Results show that the total Toc content and the γ-/α-Toc ratio varied greatly across the accessions from 85.34 to 387.00 mg/mg and 0.65 to 5.03, respectively. Furthermore, we conducted genome-wide association studies on the Tocs, which identified 28 and 73 single nucleotide polymorphisms significantly associated with the variation of total Toc content and γ-/α-Toc ratio, respectively. Bna.C02.VTE4, a putative orthologue of Arabidopsis VITAMIN E DEFICIENT 4, was tightly associated with the γ-/α-Toc ratio. This study recommends specific genetic materials with particularly high total Toc and/or low γ-/α-Toc ratio and the molecular markers and haplotypes associated with these quality traits for rapeseed breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01394-0.

Comparative Anatomical and Transcriptomics Reveal the Larger Cell Size as a Major Contributor to Larger Fruit Size in Apricot

Fruit size is one of the essential quality traits and influences the economic value of apricots. To explore the underlying mechanisms of the formation of differences in fruit size in apricots, we performed a comparative analysis of anatomical and transcriptomics dynamics during fruit growth and development in two apricot cultivars with contrasting fruit sizes (large-fruit Prunus armeniaca 'Sungold' and small-fruit P. sibirica 'F43'). Our analysis identified that the difference in fruit size was mainly caused by the difference in cell size between the two apricot cultivars. Compared with 'F43', the transcriptional programs exhibited significant differences in 'Sungold', mainly in the cell expansion period. After analysis, key differentially expressed genes (DEGs) most likely to influence cell size were screened out, including genes involved in auxin signal transduction and cell wall loosening mechanisms. Furthermore, weighted gene co-expression network analysis (WGCNA) revealed that PRE6/bHLH was identified as a hub gene, which interacted with 1 TIR1, 3 AUX/IAAs, 4 SAURs, 3 EXPs, and 1 CEL. Hence, a total of 13 key candidate genes were identified as positive regulators of fruit size in apricots. The results provide new insights into the molecular basis of fruit size control and lay a foundation for future breeding and cultivation of larger fruits in apricot.

Genome assembly of wild loquat (Eriobotrya japonica) and resequencing provide new insights into the genomic evolution and fruit domestication in loquat

Wild loquats (Eriobotrya japonica Lindl.) provide remarkable genetic resources for studying domestication and breeding improved varieties. Herein, we generate the first high-quality chromosome-level genome assembly of wild loquat, with 45 791 predicted protein-coding genes. Analysis of comparative genomics indicated that loquat shares a common ancestor with apple and pear, and a recent whole-genome duplication event occurred in loquat prior to its divergence. Genome resequencing showed that the loquat germplasms can be distinctly classified into wild and cultivated groups, and the commercial cultivars have experienced allelic admixture. Compared with cultivated loquats, the wild loquat genome showed very few selected genomic regions and had higher levels of genetic diversity. However, whole-genome scans of selective sweeps were mainly related to fruit quality, size, and flesh color during the domestication process. Large-scale transcriptome and metabolome analyses were further performed to identify differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) in wild and cultivated loquats at various fruit development stages. Unlike those in wild loquat, the key DEGs and DAMs involved in carbohydrate metabolism, plant hormone signal transduction, flavonoid biosynthesis, and carotenoid biosynthesis were significantly regulated in cultivated loquats during fruit development. These high-quality reference genome, resequencing, and large-scale transcriptome/metabolome data provide valuable resources for elucidating fruit domestication and molecular breeding in loquat.

Genome-wide association study of 23 flowering phenology traits and 4 floral agronomic traits in tree peony (Paeonia section Moutan DC.) reveals five genes known to regulate flowering time

Tree peony is a unique traditional flower in China, with large, fragrant, and colorful flowers. However, a relatively short and concentrated flowering period limits the applications and production of tree peony. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of flowering phenology traits and ornamental phenotypes in tree peony. A diverse panel of 451 tree peony accessions was phenotyped for 23 flowering phenology traits and 4 floral agronomic traits over 3 years. Genotyping by sequencing (GBS) was used to obtain a large number of genome-wide single-nucleotide polymorphisms (SNPs) (107 050) for the panel genotypes, and 1047 candidate genes were identified by association mapping. Eighty-two related genes were observed during at least 2 years for flowering, and seven SNPs repeatedly identified for multiple flowering phenology traits over multiple years were highly significantly associated with five genes known to regulate flowering time. We validated the temporal expression profiles of these candidate genes and highlighted their possible roles in the regulation of flower bud differentiation and flowering time in tree peony. This study shows that GWAS based on GBS can be used to identify the genetic determinants of complex traits in tree peony. The results expand our understanding of flowering time control in perennial woody plants. Identification of markers closely related to these flowering phenology traits can be used in tree peony breeding programs for important agronomic traits.

Whole-genome resequencing reveals genetic differences and the genetic basis of parapodium number in Russian and Chinese Apostichopus japonicus

BACKGROUND

Apostichopus japonicus is an economically important species in the global aquaculture industry. Russian A. japonicus, mainly harvested in the Vladivostok region, exhibits significant phenotypic differentiation, including in many economically important traits, compared with Chinese A. japonicus owing to differences in their habitat. However, both the genetic basis for the phenotypic divergence and the population genetic structure of Russian and Chinese A. japonicus are unknown.

RESULT

In this study, 210 individuals from seven Russian and Chinese A. japonicus populations were sampled for whole-genome resequencing. The genetic structure analysis differentiated the Russian and Chinese A. japonicus into two groups. Population genetic analyses indicated that the Russian population showed a high degree of allelic linkage and had undergone stronger positive selection compared with the Chinese populations. Gene ontology terms enriched among candidate genes with group selection analysis were mainly involved in immunity, such as inflammatory response, antimicrobial peptides, humoral immunity, and apoptosis. Genome-wide association analysis yielded eight single-nucleotide polymorphism loci significantly associated with parapodium number, and these loci are located in regions with a high degree of genomic differentiation between the Chinese and Russia populations. These SNPs were associated with five genes. Gene expression validation revealed that three of these genes were significantly differentially expressed in individuals differing in parapodium number. AJAP08772 and AJAP08773 may directly affect parapodium production by promoting endothelial cell proliferation and metabolism, whereas AJAP07248 indirectly affects parapodium production by participating in immune responses.

CONCLUSIONS

This study, we performed population genetic structure and GWAS analysis on Chinese and Russian A. japonicus, and found three candidate genes related to the number of parapodium. The results provide an in-depth understanding of the differences in the genetic structure of A. japonicus populations in China and Russia, and provide important information for subsequent genetic analysis and breeding of this species.

Genome-wide association analysis revealed genetic variation and candidate genes associated with the yield traits of upland cotton under drought conditions

Drought is one of the major abiotic stresses seriously affecting cotton yield. At present, the main cotton-producing areas in China are primarily arid and semiarid regions. Therefore, the identification of molecular markers and genes associated with cotton yield traits under drought conditions is of great importance for stabilize cotton yield under such conditions. In this study, resequencing data were used to conduct a genome-wide association study (GWAS) on 8 traits of 150 cotton germplasms. Under drought stress, 18 SNPs were significantly correlated with yield traits (single-boll weight (SBW) and seed (SC)), and 8 SNPs were identified as significantly correlated with effective fruit shoot number (EFBN) traits (a trait that is positively correlated with yield). Finally, a total of 15 candidate genes were screened. The combined results of the GWAS and transcriptome data analysis showed that four genes were highly expressed after drought stress, and these genes had significantly increased expression at 10, 15 and 25 DPA of fiber development. qRT-PCR was performed on two samples with drought tolerance extremes (drought-resistant Xinluzao 45 and drought-sensitive Xinluzao 26), revealing that three of the genes had the same differential expression pattern. This study provides a theoretical basis for the genetic analysis of cotton yield traits under drought stress, and provides gene resources for improved breeding of cotton yield traits under drought stress.

Leaf physiology variations are modulated by natural variations that underlie stomatal morphology in Populus

Stomata are essential for photosynthesis and abiotic stress tolerance. Here, we used multiomics approaches to dissect the genetic architecture and adaptive mechanisms that underlie stomatal morphology in Populus tomentosa juvenile natural population (303 accessions). We detected 46 candidate genes and 15 epistatic gene-pairs, associated with 5 stomatal morphologies and 18 leaf development and photosynthesis traits, through genome-wide association studies. Expression quantitative trait locus mapping revealed that stomata-associated gene loci were significantly associated with the expression of leaf-related genes; selective sweep analysis uncovered significant differentiation in the allele frequencies of genes that underlie stomatal variations. An allelic regulatory network operating under drought stress and adequate precipitation conditions, with three key regulators (DUF538, TRA2 and AbFH2) and eight interacting genes, was identified that might regulate leaf physiology via modulation of stomatal shape and density. Validation of candidate gene variations in drought-tolerant and F₁ hybrid populations of P. tomentosa showed that the DUF538, TRA2 and AbFH2 loci cause functional stabilisation of spatiotemporal regulatory, whose favourable alleles can be faithfully transmitted to offspring. This study provides insights concerning leaf physiology and stress tolerance via the regulation of stomatal determination in perennial plants.

Domestication and selection footprints in Persian walnuts (Juglans regia)

Walnut (Juglans) species are economically important hardwood trees cultivated worldwide for both edible nuts and high-quality wood. Broad-scale assessments of species diversity, evolutionary history, and domestication are needed to improve walnut breeding. In this study, we sequenced 309 walnut accessions from around the world, including 55 Juglans relatives, 98 wild Persian walnuts (J. regia), 70 J. regia landraces, and 86 J. regia cultivars. The phylogenetic tree indicated that J. regia samples (section Dioscaryon) were monophyletic within Juglans. The core areas of genetic diversity of J. regia germplasm were southwestern China and southern Asia near the Qinghai-Tibet Plateau and the Himalayas, and the uplift of the Himalayas was speculated to be the main factor leading to the current population dynamics of Persian walnut. The pattern of genomic variation in terms of nucleotide diversity, linkage disequilibrium, single nucleotide polymorphisms, and insertions/deletions revealed the domestication and selection footprints in Persian walnut. Selective sweep analysis, GWAS, and expression analysis further identified two transcription factors, JrbHLH and JrMYB6, that influence the thickness of the nut diaphragm as loci under selection during domestication. Our results elucidate the domestication and selection footprints in Persian walnuts and provide a valuable resource for the genomics-assisted breeding of this important crop.

Genome-wide association studies reveal genetic basis of ionomic variation in cassava

As one of the most important food crops, cassava (Manihot esculenta) is the main dietary source of micronutrients for about 1 billion people. However, the ionomic variation in cassava and the underlying genetic mechanisms remain unclear so far. Herein, genome-wide association studies were performed to reveal the specific single nucleotide polymorphisms (SNPs) that affect the ionomic variation in cassava. We identified 164 SNPs with P-values lower than the threshold located in 88 loci associated with divergent ionomic variations. Among them, 13 SNPs are related to both calcium (Ca) and magnesium (Mg), and many loci for different ionomic traits seem to be clustered on specific chromosome regions. Moreover, we identified the peak SNPs in the promoter regions of Sc10g003170 (encoding methionyl-tRNA synthetase [MetRS]) and Sc18g015190 (encoding the transcriptional regulatory protein AlgP) for nitrogen (N) and phosphorus (P) accumulation, respectively. Notably, these two SNPs (chr10_32807962 and chr18_31343738) were directly correlated with the transcript levels of Sc10g003170 (MetRS) and Sc18g015190 (AlgP), which positively modulated N accumulation and P concentration in cassava, respectively. Taken together, this study provides important insight into the genetic basis of cassava natural ionomic variation, which will promote genetic breeding to improve nutrient use and accumulation of elements in cassava.

Exploring the Genetic Association between Obesity and Serum Lipid Levels Using Bivariate Methods

It is crucial to understand the genetic mechanisms and biological pathways underlying the relationship between obesity and serum lipid levels. Structural equation models (SEMs) were constructed to calculate heritability for body mass index (BMI), total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and the genetic connections between BMI and the four classes of lipids using 1197 pairs of twins from the Chinese National Twin Registry (CNTR). Bivariate genomewide association studies (GWAS) were performed to identify genetic variants associated with BMI and lipids using the records of 457 individuals, and the results were further validated in 289 individuals. The genetic background affecting BMI may differ by gender, and the heritability of males and females was 71% (95% CI [.66, .75]) and 39% (95% CI [.15, .71]) respectively. BMI was positively correlated with TC, TG and LDL-C in phenotypic and genetic correlation, while negatively correlated with HDL-C. There were gender differences in the correlation between BMI and lipids. Bivariate GWAS analysis and validation stage found 7 genes (LOC105378740, LINC02506, CSMD1, MELK, FAM81A, ERAL1 and MIR144) that were possibly related to BMI and lipid levels. The significant biological pathways were the regulation of cholesterol reverse transport and the regulation of high-density lipoprotein particle clearance (p < .001). BMI and blood lipid levels were affected by genetic factors, and they were genetically correlated. There might be gender differences in their genetic correlation. Bivariate GWAS analysis found MIR144 gene and its related biological pathways may influence obesity and lipid levels.

Potential of traditional Chilean blood-fleshed peach to support livelihood opportunities in local agriculture

The blood-flesh peach or vineyard peach is an older heritage cultivar with juicy red-flesh and tart-sweet flavor. They are popular in France, where more than 200 years ago wine growers used to plant them on the vineyards as biological markers to detect the presence of powdery mildew. It is present in countries such as China, Italy, New Zealand, Australia and USA however, it remains a very rare variety worldwide. In Chile, the blood-flesh peach has a centenary presence in rural orchards where is called “Durazno Betarraga.” Reproduced by seeds, it has pass through generations of family farmers and has been adapted to local environmental conditions. This red-flesh peach is a local variety considered part of their traditional diets, however, cultural changes in food consumption, short postharvest life and water scarcity due to climate change are threatening its conservation. One of the objectives of the International Year of Fruits and Vegetables, as defined by the FAO, is to integrate small holders and family farmers into value chains for sustainable production and consumption of fruits and vegetables recognizing the contributions of farmer's landraces to their food security, nutrition, livelihoods and income. To promote this objective, we present the work we have been carry out for several years with a farming community. We have conducted ethnographic research to provide a qualitative description of the agricultural value of the blood peach in a limited territory of the Maule Region defined as the study area. For the quantitative section of our research we analyzed the antioxidant capacity (ORAC) and total polyphenol content and compared them with those of other fruits. To gather information on the presence of the blood-fleshed peach in other regions of Chile, we used a citizen science approach through social networks. We propose that this local variety is an innovative raw material to develop healthy fruit-based food, thus encouraging its conservation and consumption with a positive social and economic impact for the community and the local food system.

Genome-Wide Analysis of SAUR Gene Family Identifies a Candidate Associated with Fruit Size in Loquat (Eriobotrya japonica Lindl.)

Fruit size is an important fruit quality trait that influences the production and commodity values of loquats (Eriobotrya japonica Lindl.). The Small Auxin Upregulated RNA (SAUR) gene family has proven to play a vital role in the fruit development of many plant species. However, it has not been comprehensively studied in a genome-wide manner in loquats, and its role in regulating fruit size remains unknown. In this study, we identified 95 EjSAUR genes in the loquat genome. Tandem duplication and segmental duplication contributed to the expansion of this gene family in loquats. Phylogenetic analysis grouped the SAURs from Arabidopsis, rice, and loquat into nine clusters. By analyzing the transcriptome profiles in different tissues and at different fruit developmental stages and comparing two sister lines with contrasting fruit sizes, as well as by functional predictions, a candidate gene (EjSAUR22) highly expressed in expanding fruits was selected for further functional investigation. A combination of Indoleacetic acid (IAA) treatment and virus-induced gene silencing revealed that EjSAUR22 was not only responsive to auxin, but also played a role in regulating cell size and fruit expansion. The findings from our study provide a solid foundation for understanding the molecular mechanisms controlling fruit size in loquats, and also provide potential targets for manipulation of fruit size to accelerate loquat breeding.

Integrative analysis of the metabolome and transcriptome reveals the molecular mechanism of chlorogenic acid synthesis in peach fruit

As the most abundant phenolic acid in peach fruit, chlorogenic acid (CGA) is an important entry point for the development of natural dietary supplements and functional foods. However, the metabolic and regulation mechanisms underlying its accumulation in peach fruits remain unclear. In this study, we evaluated the composition and content of CGAs in mature fruits of 205 peach cultivars. In peach fruits, three forms of CGA (52.57%), neochlorogenic acid (NCGA, 47.13%), and cryptochlorogenic acid (CCGA, 0.30%) were identified. During the growth and development of peach fruits, the content of CGAs generally showed a trend of rising first and then decreasing. Notably, the contents of quinic acid, shikimic acid, p-coumaroyl quinic acid, and caffeoyl shikimic acid all showed similar dynamic patterns to that of CGA, which might provide the precursor material basis for the accumulation of CGA in the later stage. Moreover, CGA, lignin, and anthocyanins might have a certain correlation and these compounds work together to maintain a dynamic balance. By the comparative transcriptome analysis, 8 structural genes (Pp4CL, PpCYP98A, and PpHCT) and 15 regulatory genes (PpMYB, PpWRKY, PpERF, PpbHLH, and PpWD40) were initially screened as candidate genes of CGA biosynthesis. Our findings preliminarily analyzed the metabolic and molecular regulation mechanisms of CGA biosynthesis in peach fruit, which provided a theoretical basis for developing high-CGA content peaches in future breeding programs.

Combined nature and human selections reshaped peach fruit metabolome

Background

Plant metabolites reshaped by nature and human beings are crucial for both their lives and human health. However, which metabolites respond most strongly to selection pressure at different evolutionary stages and what roles they undertake on perennial fruit crops such as peach remain unclear.

Results

Here, we report 18,052 significant locus-trait associations, 12,691 expression-metabolite correlations, and 294,676 expression quantitative trait loci (eQTLs) for peach. Our results indicate that amino acids accumulated in landraces may be involved in the environmental adaptation of peaches by responding to low temperature and drought. Moreover, the contents of flavonoids, the major nutrients in fruits, have kept decreasing accompanied by the reduced bitter flavor during both domestication and improvement stages. However, citric acid, under the selection of breeders’ and consumers’ preference for flavor, shows significantly different levels between eastern and western varieties. This correlates with differences in activity against cancer cells in vitro in fruit from these two regions. Based on the identified key genes regulating flavonoid and acid contents, we propose that more precise and targeted breeding technologies should be designed to improve peach varieties with rich functional contents because of the linkage of genes related to bitterness and acid taste, antioxidant and potential anti-cancer activity that are all located at the top of chromosome 5.

Conclusions

This study provides powerful data for future improvement of peach flavor, nutrition, and resistance in future and expands our understanding of the effects of natural and artificial selection on metabolites.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02719-6.

Genetic diversity and population dynamic of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow in Central China

Phylogeographic research concerning Central China has been rarely conducted. Population genetic and phylogeography of Ziziphus jujuba var. spinosa (also called sour jujube) were investigated to improve our understanding of plant phylogeographic patterns in Central China. Single-copy nuclear gene markers and complete chloroplast genome data were applied to 328 individuals collected from 21 natural populations of sour jujube in China. Nucleotide variation of sour jujube was relatively high (π = 0.00720, θ w = 0.00925), which resulted from the mating system and complex population dynamics. Analysis of molecular variation analysis revealed that most of the total variation was attributed to variation within populations, and a high level of genetic differentiation among populations was detected (F st = 0.197). Relatively low long-distance dispersal capability and vitality of pollen contributed to high genetic differentiation among populations. Differences in the environmental conditions and long distance among populations further restricted gene flow. Structure clustering analysis uncovered intraspecific divergence between central and marginal populations. Migrate analysis found a high level of gene flow between these two intraspecific groups. Bayesian skyline plot detected population expansion of these two intraspecific groups. Network and phylogeny analysis of chloroplast haplotypes also found intraspecific divergence, and the divergence time was estimated to occur at about 55.86 Ma. Haplotype native to the Loess Plateau was more ancient, and multiple glacial refugia of sour jujube were found to locate at the Loess Plateau, areas adjacent to the Qinling Mountains and Tianmu Mountains. Species distribution model analysis found a typical contraction-expansion model corresponding to the Quaternary climatic oscillations. In the future, the distribution of sour jujube may shift to high-latitude areas. This study provides new insights for phylogeographic research of temperate plant species distributed in Central China and sets a solid foundation for the application of the scientific management strategy of Z. jujuba var. spinosa.

A de novo genome assembly of cultivated Prunus persica cv. ‘Sovetskiy’

Prunus persica is one of the main stone fruit crops in Crimea and southern Russia. The P. persica genome has recently been sequenced and annotated in good quality. However, for a deeper assessment of the peach genome, it is necessary to include in the research other cultivars that are in the collection of the Nikitsky Botanical Garden. The cultivars of the Nikitsky Botanical Garden are unique and differ from Western European and American ones, as they are derived from cultivars and forms originating from Central Asian, North Caucasian, Transcaucasian and Eastern European countries. In this paper, we present the assembly of the P. persica cv. ’Sovetskiy’ genome obtained using Oxford Nanopore long reads and Illumina short reads by hybrid assembly methods. The assembled genome of P. persica cv. ’Sovetskiy’ is 206.26 MB in 226 scaffolds, with N50 24 Mb, including 8 chromosomes. It contains 27140 coding genes, 26973 (99.38%) of which are annotated in at least one functional database. More than 36.05% of the genome regions were identified as repeating elements.

Chromosome-level genome assemblies of four wild peach species provide insights into genome evolution and genetic basis of stress resistance

Background

Peach (Prunus persica) is an economically important stone fruit crop in Rosaceae and widely cultivated in temperate and subtropical regions, emerging as an excellent material to study the interaction between plant and environment. During its genus, there are four wild species of peach, all living in harsh environments. For example, one of the wild species, P. mira, originates from the Qinghai-Tibet Plateau (QTP) and exhibits strong cold/ultraviolet ray environmental adaptations. Although remarkable progresses in the gene discovery of fruit quality-related traits in peach using previous assembled genome were obtained, genomic basis of the response of these wild species to different geographical environments remains unclear.

Results

To uncover key genes regulating adaptability in different species and analyze the role of genetic variations in resistance formation, we performed de novo genome assembling of four wild relatives of peach (P. persica), P. mira, P. davidiana, P. kansuensis, and P. ferganensis and resequenced 175 peach varieties. The phylogenetic tree showed that the divergence time of P. mira and other wild relatives of peach was 11.5 million years ago, which was consistent with the drastic crustal movement of QTP. Abundant genetic variations were identified in four wild species when compared to P. persica, and the results showed that plant-pathogen interaction pathways were enriched in genes containing small insertions and deletions and copy number variations in all four wild relatives of peach. Then, the data were used to identify new genes and variations regulating resistance. For example, presence/absence variations which result from a hybridization event that occurred between P. mira and P. dulcis enhanced the resistance of their putative hybrid, P. davidiana. Using bulked segregant analysis, we located the nematode resistance locus of P. kansuensis in chromosome 2. Within the mapping region, a deletion in the promoter of one NBS-LRR gene was found to involve the resistance by regulating gene expression. Furthermore, combined with RNA-seq and selective sweeps analysis, we proposed that a deletion in the promoter of one CBF gene was essential for high-altitude adaptation of P. mira through increasing its resistance to low temperature.

Conclusions

In general, the reference genomes assembled in the study facilitate our understanding of resistance mechanism of perennial fruit crops, and provide valuable resources for future breeding and improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01342-y.

Identification of QTLs Associated With Agronomic Traits in Tobacco via a Biparental Population and an Eight-Way MAGIC Population

Agronomic traits such as plant height (PH), leaf number (LN), leaf length (LL), and leaf width (LW), which are closely related to yield and quality, are important in tobacco (Nicotiana tabacum L.). To identify quantitative trait loci (QTLs) associated with agronomic traits in tobacco, 209 recombinant inbred lines (RILs) and 537 multiparent advanced generation intercross (MAGIC) lines were developed. The biparental RIL and MAGIC lines were genotyped using a 430 K single-nucleotide polymorphism (SNP) chip assay, and their agronomic traits were repeatedly evaluated under different conditions. A total of 43 QTLs associated with agronomic traits were identified through a combination of linkage mapping (LM) and association mapping (AM) methods. Among these 43 QTLs, three major QTLs, namely qPH13-3, qPH17-1, and qLW20-1, were repeatedly identified by the use of various genetically diverse populations across different environments. The candidate genes for these major QTLs were subsequently predicted. Validation and utilization of the major QTL qLW20-1 for the improvement of LW in tobacco were investigated. These results could be applied to molecular marker-assisted selection (MAS) for breeding important agronomic traits in tobacco.

Impacts of Chronic Habitat Fragmentation on Genetic Diversity of Natural Populations of Prunus persica in China

Wild peach is an important resource for improving existing peach varieties. However, the extant populations of wild peach show fragmented distribution due to human disturbance and geographic isolation. In this study, we used natural populations (or wild populations) of Prunus persica (Rosaceae) to assess the genetic effects of habitat fragmentation. A total of 368 individuals sampled from 16 natural populations were analyzed using 23 polymorphic simple sequence repeat (SSR) markers. Prunus persica maintained low within-population genetic variation and high level of genetic differentiation. Two genetic clusters were revealed based on three different methods (UPGMA, PCoA, and STRUCTURE). All populations showed a significant heterozygosity deficiency and most extant populations experienced recent reduction in population size. A significant isolation by distance (IBD) was observed with Mantel's test. Compared to historical gene flow, contemporary gene flow was restricted among the studied populations, suggesting a decrease in gene flow due to habitat fragmentation. Habitat fragmentation has impacted population genetic variation and genetic structure of P. persica. For breeding and conservation purpose, collecting as many individuals as possible from multiple populations to maximize genetic diversity was recommended during the process of germplasm collection. In addition, populations from central China had higher genetic diversity, suggesting these populations should be given priority for conservation and germplasm collection.

Identification of Genetic Loci for Sugarcane Leaf Angle at Different Developmental Stages by Genome-Wide Association Study

Sugarcane (Saccharum spp.) is an efficient crop mainly used for sugar and bioethanol production. High yield and high sucrose of sugarcane are always the fundamental demands in sugarcane growth worldwide. Leaf angle and size of sugarcane can be attributed to planting density, which was associated with yield. In this study, we performed genome-wide association studies (GWAS) with a panel of 216 sugarcane core parents and their derived lines (natural population) to determine the genetic basis of leaf angle and key candidate genes with +2, +3, and +4 leaf at the seedling, elongation, and mature stages. A total of 288 significantly associated loci of sugarcane leaf angle at different developmental stages (eight phenotypes) were identified by GWAS with 4,027,298 high-quality SNP markers. Among them, one key locus and 11 loci were identified in all three stages and two stages, respectively. An InDel marker (SNP Ss6A_102766953) linked to narrow leaf angle was obtained. Overall, 4,089 genes were located in the confidence interval of significant loci, among which 3,892 genes were functionally annotated. Finally, 13 core parents and their derivatives tagged with SNPs were selected for marker-assisted selection (MAS). These candidate genes are mainly related to MYB transcription factors, auxin response factors, serine/threonine protein kinases, etc. They are directly or indirectly associated with leaf angle in sugarcane. This research provided a large number of novel genetic resources for the improvement of leaf angles and simultaneously to high yield and high bioethanol production.

De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size

Peach (Prunus persica) is one of the most important fruit crops globally, but its cultivation can be hindered by large tree size. 'Zhongyoutao 14' (CN14) is a temperature-sensitive semi-dwarf (TSSD) cultivar which might be useful as breeding stock. The genome of CN14 was sequenced and assembled de novo using single-molecule real-time sequencing and chromosome conformation capture assembly. A high-quality genome was assembled and annotated, with 228.82 Mb mapped to eight chromosomes. Eighty-six re-sequenced F₁ individuals and 334 previously re-sequenced accessions were used to identify candidate genes controlling TSSD and flower type and size. An aquaporin tonoplast intrinsic protein (PpTIP2) was a strong candidate gene for control of TSSD. Sequence variations in the upstream regulatory region of PpTIP2 correlated with different transcriptional activity at different temperatures. PpB3-1, a candidate gene for flower type (SH) and flower size, contributed to petal development and promoted petal enlargement. The locus of another 12 agronomic traits was identified through genome-wide association study. Most of these loci exhibited consistent and precise association signals, except for flesh texture and flesh adhesion. A 6015-bp insertion in exon 3 and a 26-bp insertion upstream of PpMYB25 were associated with fruit hairless. Along with a 70.5-Kb gap at the F-M locus in CN14, another two new alleles were identified in peach accessions. Our findings will not only promote genomic research and agronomic breeding in peach but also provide a foundation for the peach pan-genome.

Integrating Genome-wide association and whole transcriptome analysis to reveal genetic control of leaf traits in Gossypium arboreum L

Leaves are important organs for crop photosynthesis and transpiration, and their morphological characteristics can directly reflect the growth state of plants. Accurate measurement of leaf traits and mining molecular markers are of great significance to the study of cotton growth. Here, we performed a Genome-wide association study on 7 leaf traits in 213 Asian cotton accessions. 32 significant SNPs and 44 genes were identified. A field experiment showed significant difference in leaf hair and leaf area between DPL971 and its natural mutant DPL972. We also compared the leaf transcriptome difference between DPL971 and DPL972, and found a batch of differentially expressed genes and non-coding RNAs (including lncRNAs, microRNAs, and circRNAs). After integrating the GWAS and transcriptome results, we finally selected two coding genes (Ga03G2383 and Ga05G3412) and two microRNAs (hbr-miR156, unconservative_Chr03_contig343_2364) as the candidate for leaf traits. Those findings will provide important genomic resources for cotton leaf improvement breeding.

De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size

Peach (Prunus persica) is one of the most important fruit crops globally, but its cultivation can be hindered by large tree size. 'Zhongyoutao 14' (CN14) is a temperature-sensitive semi-dwarf (TSSD) cultivar which might be useful as breeding stock. The genome of CN14 was sequenced and assembled de novo using single-molecule real-time sequencing and chromosome conformation capture assembly. A high-quality genome was assembled and annotated, with 228.82 Mb mapped to eight chromosomes. Eighty-six re-sequenced F₁ individuals and 334 previously re-sequenced accessions were used to identify candidate genes controlling TSSD and flower type and size. An aquaporin tonoplast intrinsic protein (PpTIP2) was a strong candidate gene for control of TSSD. Sequence variations in the upstream regulatory region of PpTIP2 correlated with different transcriptional activity at different temperatures. PpB3-1, a candidate gene for flower type (SH) and flower size, contributed to petal development and promoted petal enlargement. The locus of another 12 agronomic traits was identified through genome-wide association study. Most of these loci exhibited consistent and precise association signals, except for flesh texture and flesh adhesion. A 6015-bp insertion in exon 3 and a 26-bp insertion upstream of PpMYB25 were associated with fruit hairless. Along with a 70.5-Kb gap at the F-M locus in CN14, another two new alleles were identified in peach accessions. Our findings will not only promote genomic research and agronomic breeding in peach but also provide a foundation for the peach pan-genome.

Multi-omics approaches identify a key gene, PpTST1, for organic acid accumulation in peach

Organic acid content in fruit is an important determinant of peach organoleptic quality, which undergoes considerable variations during development and maturation. However, its molecular mechanism remains largely unclear. In this study, an integrative approach of genome-wide association studies and comparative transcriptome analysis were applied to identify candidate genes involved in organic acid accumulation in peach. A key gene PpTST1, encoding tonoplast sugar transporter, was identified and the genotype of PpTST1 with a single-base transversion (G1584T) in the third exon which leads to a single amino acid substitution (Q528H) was associated with low level of organic acid content in peach. Overexpression of PpTST1His resulted in reduced organic acid content along with increased sugar content both in peach and tomato fruits, suggesting its dual function in sugar accumulation and organic acid content reduction. Two V-type proton ATPases interact with PpTST1 in yeast two-hybridization assay. In addition, the G1584T transversion appeared and gradually accumulated during domestication and improvement, which indicated that PpTST1 was under selection. The identification and characterization of PpTST1 would facilitate the improvement of peach fruit quality.