Comparative Transcriptomic Analysis Provides Insight into the Domestication and Improvement of Pear (P. pyrifolia) Fruit

The Transcriptional Landscape of Berry Skin in Red and White PIWI ("Pilzwiderstandsfähig") Grapevines Possessing QTLs for Partial Resistance to Downy and Powdery Mildews

PIWI, from the German word Pilzwiderstandsfähig, meaning "fungus-resistant", refers to grapevine cultivars bred for resistance to fungal pathogens such as Erysiphe necator (the causal agent of powdery mildew) and Plasmopara viticola (the causal agent of downy mildew), two major diseases in viticulture. These varieties are typically developed through traditional breeding, often crossbreeding European Vitis vinifera with American or Asian species that carry natural disease resistance. This study investigates the transcriptional profiles of exocarp tissues in mature berries from four PIWI grapevine varieties compared to their elite parental counterparts using RNA-seq analysis. We performed RNA-seq on four PIWI varieties (two red and two white) and their noble parents to identify differential gene expression patterns. Comprehensive analyses, including Differential Gene Expression (DEGs), Gene Set Enrichment Analysis (GSEA), Weighted Gene Co-expression Network Analysis (WGCNA), and tau analysis, revealed distinct gene clusters and individual genes characterizing the transcriptional landscape of PIWI varieties. Differentially expressed genes indicated significant changes in pathways related to organic acid metabolism and membrane transport, potentially contributing to enhanced resilience. WGCNA and k-means clustering highlighted co-expression modules linked to PIWI genotypes and their unique tolerance profiles. Tau analysis identified genes uniquely expressed in specific genotypes, with several already known for their defense roles. These findings offer insights into the molecular mechanisms underlying grapevine resistance and suggest promising avenues for breeding strategies to enhance disease resistance and overall grape quality in viticulture.

Increased DNA methylation contributes to the early ripening of pear fruits during domestication and improvement

BACKGROUND

DNA methylation is an essential epigenetic modification. However, its contribution to trait changes and diversity in the domestication of perennial fruit trees remains unknown.

RESULTS

Here, we investigate the variation in DNA methylation during pear domestication and improvement using whole-genome bisulfite sequencing in 41 pear accessions. Contrary to the significant decrease during rice domestication, we detect a global increase in DNA methylation during pear domestication and improvement. We find this specific increase in pear is significantly correlated with the downregulation of Demeter-like1 (DML1, encoding DNA demethylase) due to human selection. We identify a total of 5591 differentially methylated regions (DMRs). Methylation in the CG and CHG contexts undergoes co-evolution during pear domestication and improvement. DMRs have higher genetic diversity than selection sweep regions, especially in the introns. Approximately 97% of DMRs are not associated with any SNPs, and these DMRs are associated with starch and sucrose metabolism and phenylpropanoid biosynthesis. We also perform correlation analysis between DNA methylation and gene expression. We find genes close to the hypermethylated DMRs that are significantly associated with fruit ripening. We further verify the function of a hyper-DMR-associated gene, CAMTA2, and demonstrate that overexpression of CAMTA2 in tomato and pear callus inhibits fruit ripening.

CONCLUSIONS

Our study describes a specific pattern of DNA methylation in the domestication and improvement of a perennial pear tree and suggests that increased DNA methylation plays an essential role in the early ripening of pear fruits.

Salicylic Acid Spray Delays Sand Pear Fruit Senescence during Room Temperature Shelf Life by Regulating Antioxidant Capacity and Senescence-Related Genes

'Whangkeumbae' (Pyrus pyrifolia) is a variety of sand pear fruit well-known for its smooth surface and good taste. However, the fruit quality is adversely affected by postharvest ethylene production. Therefore, improving postharvest shelf life by regulating fruit senescence is critical to promoting the 'Whangkeumbae' fruit industry. Here, we investigated the effect of salicylic acid (SA) spray on fruit senescence in sand pears during room temperature shelf life. Exogenous SA reduced polyphenol oxidase (PPO) activity and malondialdehyde (MDA) content during room temperature shelf life. Additionally, SA effectively maintained the fruit skin coloration and increased the activity of antioxidant enzymes, such as superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), and ascorbate peroxidase (APX). SA treatment inhibited PpPPO1 expression and upregulated PpSOD1, PpAPX6, and PpGST2 expression. Furthermore, SA application downregulated the expression of PpACO2, PpEIN3a, PpNCED1, and PpAOC2, while upregulating PpNPR-1, PpTAR2, and PpCOMT1 during room temperature shelf life. SA treatment also influenced cell wall metabolism and modification genes by inhibiting PpPG1, PpPME2, and PpCEL3 and inducing PpPGIP1 expression. Additionally, SA treatment affected sugar and acid metabolism genes and increased the expression of PpSPS1, PpSUS1, PpSOT1, PpTMT4, PpSWEET15, and PpcyNAD-MDH, but suppressed the expression of PpcyNADP-ME. The Pearson correlation analysis indicated that PPO activity and MDA content were positively correlated with the expression of PpPPO1, PpACO2, PpEIN3a, PpNCED1, PpAOC2, PpPG1, PpPME2, PpCEL3, and PpcyNDA-MDH. Conversely, these factors were negatively associated with the activities of SOD, POD, CAT, and APX, as well as the expression levels of PpSOD1, PpPOD1, PpCAT1, PpAPX6, PpGST2, PpNPR-1, PpTAR2, PpCOMT1, PpPGIP1, PpSPS1, PpSUS1, PpSOT1, PpTMT4, PpSWEET15, and PpcyNAD-MDH. Our results reveal that exogenous SA could delay fruit senescence in sand pear fruit by regulating various biochemical and molecular mechanisms and can be used to effectively extend fruit shelf life during room temperature storage. However, further research is necessary to determine whether the fruits sprayed with SA are suitable for direct human consumption.

PbrbZIP15 promotes sugar accumulation in pear via activating the transcription of the glucose isomerase gene PbrXylA1

The composition and abundance of soluble sugars in mature pear (Pyrus) fruit are important for its acceptance by consumers. However, our understanding of the genes responsible for soluble sugar accumulation remains limited. In this study, a S1-group member of bZIP gene family, PbrbZIP15, was characterized from pear genome through the combined analyses of metabolite and transcriptome data followed by experimental validation. PbrbZIP15, located in nucleus, was found to function in fructose, sucrose, and total soluble sugar accumulation in pear fruit and calli. After analyzing the expression profiles of sugar-metabolism-related genes and the distribution of cis-acting elements in their promoters, the glucose isomerase 1 gene (PbrXylA1), whose corresponding protein catalyzed the isomerization of glucose and fructose in vitro, was identified as a downstream target gene of PbrbZIP15. PbrbZIP15 could directly bind to the G-box element in PbrXylA1 promoter and activate its transcription, as evidenced by chromatin immunoprecipitation-quantitative PCR, yeast one-hybrid, electrophoretic mobility shift assay, and dual-luciferase assay. PbrXylA1, featuring a leucine-rich signal peptide in its N-terminal, was localized to the endoplasmic reticulum. It was validated to play a significant role in fructose, sucrose, and total soluble sugar accumulation in pear fruit and calli, which was associated with the upregulated fructose/glucose ratio. Further studies revealed a positive correlation between the sucrose content and the expression levels of several sucrose-biosynthesis-related genes (PbrFRK3/8, PbrSPS1/3/4/8, and PbrSPP1) in PbrbZIP15-/PbrXylA1-transgenic fruit/calli. In conclusion, our results suggest that PbrbZIP15-induced soluble sugar accumulation during pear development is at least partly attributed to the activation of PbrXylA1 transcription.

Comparative physiological, metabolomic and transcriptomic analyses reveal the mechanisms of differences in pear fruit quality between distinct training systems

BACKGROUND

Canopy architecture is critical in determining the fruit-zone microclimate and, ultimately, in determining an orchard's success in terms of the quality and quantity of the fruit produced. However, few studies have addressed how the canopy environment leads to metabolomic and transcriptomic alterations in fruits. Designing strategies for improving the quality of pear nutritional components relies on uncovering the related regulatory mechanisms.

RESULTS

We performed an in-depth investigation of the impact of canopy architecture from physiological, metabolomic and transcriptomic perspectives by comparing pear fruits grown in a traditional freestanding system (SP) or a flat-type trellis system (DP). Physiological studies revealed relatively greater fruit sizes, soluble solid contents and titratable acidities in pear fruits from DP systems with open canopies. Nontargeted metabolite profiling was used to characterize fruits at the initial ripening stage. Significant differences in fruit metabolites, including carbohydrates, nucleic acids, alkaloids, glycerophospholipids, sterol lipids, and prenol lipids, were observed between the two groups. Transcriptomic analysis indicated that a series of organic substance catabolic processes (e.g., the glycerol-3-phosphate catabolic process, pectin catabolic process and glucan catabolic process) were overrepresented in fruits of the DP system. Moreover, integrative analysis of the metabolome and transcriptome at the pathway level showed that DP pear fruits may respond to the canopy microenvironment by upregulating phenylpropanoid biosynthesis pathway genes such as PpPOD. Transient assays revealed that the contents of malic acid and citric acid were lower in the pear flesh of PpPOD RNAi plants, which was associated with regulating the expression of organic acid metabolism-related genes.

CONCLUSIONS

Our results provide fundamental evidence that at the physiological and molecular levels, open-canopy architecture contributes to improving pear fruit quality and is correlated with increased levels of carbohydrates and lipid-like molecules. This study may lead to the development of rational culture practices for enhancing the nutritional traits of pear fruits.

The pear genomics database (PGDB): a comprehensive multi-omics research platform for Pyrus spp

BACKGROUND

Pears are among the most important temperate fruit trees in the world, with significant research efforts increasing over the last years. However, available omics data for pear cannot be easily and quickly retrieved to enable further studies using these biological data.

DESCRIPTION

Here, we present a publicly accessible multi-omics pear resource platform, the Pear Genomics Database (PGDB). We collected and collated data on genomic sequences, genome structure, functional annotation, transcription factor predictions, comparative genomics, and transcriptomics. We provide user-friendly functional modules to facilitate querying, browsing and usage of these data. The platform also includes basic and useful tools, including JBrowse, BLAST, phylogenetic tree building, and additional resources providing the possibility for bulk data download and quick usage guide services.

CONCLUSIONS

The Pear Genomics Database (PGDB, http://pyrusgdb.sdau.edu.cn ) is an online data analysis and query resource that integrates comprehensive multi-omics data for pear. This database is equipped with user-friendly interactive functional modules and data visualization tools, and constitutes a convenient platform for integrated research on pear.

A large-scale proteogenomic atlas of pear

Pear is an important fruit tree that is widely distributed around the world. The first pear genome map was reported from our laboratory approximately 10 years ago. To further study global protein expression patterns in pear, we generated pear proteome data based on 24 major tissues. The tissue-resolved profiles provided evidence of the expression of 17 953 proteins. We identified 4294 new coding events and improved the pear genome annotation via the proteogenomic strategy based on 18 090 peptide spectra with peptide spectrum matches >1. Among the eight randomly selected new short coding open reading frames that were expressed in the style, four promoted and one inhibited the growth of pear pollen tubes. Based on gene coexpression module analysis, we explored the key genes associated with important agronomic traits, such as stone cell formation in fruits. The network regulating the synthesis of lignin, a major component of stone cells, was reconstructed, and receptor-like kinases were implicated as core factors in this regulatory network. Moreover, we constructed the online database PearEXP (http://www.peardb.org.cn) to enable access to the pear proteogenomic resources. This study provides a paradigm for in-depth proteogenomic studies of woody plants.

Transcriptome-based variations effectively untangling the intraspecific relationships and selection signals in Xinyang Maojian tea population

As one of the world's top three popular non-alcoholic beverages, tea is economically and culturally valuable. Xinyang Maojian, this elegant green tea, is one of the top ten famous tea in China and has gained prominence for thousands of years. However, the cultivation history of Xinyang Maojian tea population and selection signals of differentiation from the other major variety Camellia sinensis var. assamica (CSA) remain unclear. We newly generated 94 Camellia sinensis (C. sinensis) transcriptomes including 59 samples in the Xinyang area and 35 samples collected from 13 other major tea planting provinces in China. Comparing the very low resolution of phylogeny inferred from 1785 low-copy nuclear genes with 94 C. sinensis samples, we successfully resolved the phylogeny of C. sinensis samples by 99,115 high-quality SNPs from the coding region. The sources of tea planted in the Xinyang area were extensive and complex. Specifically, Shihe District and Gushi County were the two earliest tea planting areas in Xinyang, reflecting a long history of tea planting. Furthermore, we identified numerous selection sweeps during the differentiation of CSA and CSS and these positive selection genes are involved in many aspects such as regulation of secondary metabolites synthesis, amino acid metabolism, photosynthesis, etc. Numerous specific selective sweeps of modern cultivars were annotated with functions in various different aspects, indicating the CSS and CSA populations possibly underwent independent specific domestication processes. Our study indicated that transcriptome-based SNP-calling is an efficient and cost-effective method in untangling intraspecific phylogenetic relationships. This study provides a significant understanding of the cultivation history of the famous Chinese tea Xinyang Maojian and unravels the genetic basis of physiological and ecological differences between the two major tea subspecies.

Rearrangement and domestication as drivers of Rosaceae mitogenome plasticity

Background

The mitochondrion is an important cellular component in plants and that functions in producing vital energy for the cell. However, the evolution and structure of mitochondrial genomes (mitogenomes) remain unclear in the Rosaceae family. In this study, we assembled 34 Rosaceae mitogenomes and characterized genome variation, rearrangement rate, and selection signal variation within these mitogenomes.

Results

Comparative analysis of six genera from the Amygdaloideae and five from the Rosoideae subfamilies of Rosaceae revealed that three protein-coding genes were absent from the mitogenomes of five Rosoideae genera. Positive correlations between genome size and repeat content were identified in 38 Rosaceae mitogenomes. Twenty repeats with high recombination frequency (> 50%) provided evidence for predominant substoichiometric conformation of the mitogenomes. Variations in rearrangement rates were identified between eleven genera, and within the Pyrus, Malus, Prunus, and Fragaria genera. Based on population data, phylogenetic inferences from Pyrus mitogenomes supported two distinct maternal lineages of Asian cultivated pears. A Pyrus-specific deletion (DEL-D) in selective sweeps was identified based on the assembled genomes and population data. After the DEL-D sequence fragments originally arose, they may have experienced a subsequent doubling event via homologous recombination and sequence transfer in the Amygdaloideae; afterwards, this variant sequence may have significantly expanded to cultivated groups, thereby improving adaptation during the domestication process.

Conclusions

This study characterizes the variations in gene content, genome size, rearrangement rate, and the impact of domestication in Rosaceae mitogenomes and provides insights into their structural variation patterns and phylogenetic relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01383-3.

Comparative transcriptomics of wild and commercial Citrus during early ripening reveals how domestication shaped fruit gene expression

BACKGROUND

Interspecific hybridizations and admixtures were key in Citrus domestication, but very little is known about their impact at the transcriptomic level. To determine the effects of genome introgressions on gene expression, the transcriptomes of the pulp and flavedo of three pure species (citron, pure mandarin and pummelo) and four derived domesticated genetic admixtures (sour orange, sweet orange, lemon and domesticated mandarin) have been analyzed at color break.

RESULTS

Many genes involved in relevant physiological processes for domestication, such sugar/acid metabolism and carotenoid/flavonoid synthesis, were differentially expressed among samples. In the low-sugar, highly acidic species lemon and citron, many genes involved in sugar metabolism, the TCA cycle and GABA shunt displayed a reduced expression, while the P-type ATPase CitPH5 and most subunits of the vacuolar ATPase were overexpressed. The red-colored species and admixtures were generally characterized by the overexpression in the flavedo of specific pivotal genes involved in the carotenoid biosynthesis, including phytoene synthase, ζ-carotene desaturase, β-lycopene cyclase and CCD4b, a carotenoid cleavage dioxygenase. The expression patterns of many genes involved in flavonoid modifications, especially the flavonoid and phenylpropanoid O-methyltransferases showed extreme diversity. However, the most noticeable differential expression was shown by a chalcone synthase gene, which catalyzes a key step in the biosynthesis of flavonoids. This chalcone synthase was exclusively expressed in mandarins and their admixed species, which only expressed the mandarin allele. In addition, comparisons between wild and domesticated mandarins revealed that the major differences between their transcriptomes concentrate in the admixed regions.

CONCLUSION

In this work we present a first study providing broad evidence that the genome introgressions that took place during citrus domestication largely shaped gene expression in their fruits.

Pear genetics: Recent advances, new prospects, and a roadmap for the future

Pear, belonging to the genus Pyrus, is one of the most economically important temperate fruit crops. Pyrus is an important genus of the Rosaceae family, subfamily Maloideae, and has at least 22 different species with over 5000 accessions maintained or identified worldwide. With the release of draft whole-genome sequences for Pyrus, opportunities for pursuing studies on the evolution, domestication, and molecular breeding of pear, as well as for conducting comparative genomics analyses within the Rosaceae family, have been greatly expanded. In this review, we highlight key advances in pear genetics, genomics, and breeding driven by the availability of whole-genome sequences, including whole-genome resequencing efforts, pear domestication, and evolution. We cover updates on new resources for undertaking gene identification and molecular breeding, as well as for pursuing functional validation of genes associated with desirable economic traits. We also explore future directions for "pear-omics".

Genome-Wide Identification of the 1-Aminocyclopropane-1-carboxylic Acid Synthase (ACS) Genes and Their Possible Role in Sand Pear (Pyrus pyrifolia) Fruit Ripening

Ethylene production is negatively associated with storage life in sand pear (Pyrus pyrifolia Nakai), particularly at the time of fruit harvest. 1-Aminocyclopropane-1-carboxylic acid synthase (ACS) is the rate-limiting enzyme in ethylene biosynthesis and is considered to be important for fruit storage life. However, the candidate ACS genes and their roles in sand pear remain unclear. The present study identified 13 ACS genes from the sand pear genome. Phylogenetic analysis categorized these ACS genes into four subgroups (type I, type II, type III and putative AAT), and indicated a close relationship between sand pear and Chinese white pear (P. bretschneideri). According to the RNA-seq data and qRT-PCR analysis, PpyACS1, PpyACS2, PpyACS3, PpyACS8, PpyACS9, PpyACS12 and PpyACS13 were differently expressed in climacteric and non-climacteric-type pear fruits, ‘Ninomiyahakuri’ and ‘Eli No.2’, respectively, during fruit ripening. In addition, the expressions of PpyACS2, PpyACS8, PpyACS12 and PpyACS13 were found to be associated with system 1 of ethylene production, while PpyACS1, PpyACS3, and PpyACS9 were found to be associated with system 2, indicating that these ACS genes have different roles in ethylene biosynthesis during fruit development. Overall, our study provides fundamental knowledge on the characteristics of the ACS gene family in sand pear, in addition to their possible roles in fruit ripening.

Transcriptome and Resequencing Analyses Provide Insight into Differences in Organic Acid Accumulation in Two Pear Varieties

Fruit acidity is one of the main determinants of fruit flavor and a target trait in fruit breeding. However, the genomic mechanisms governing acidity variation among different pear varieties remain poorly understood. In this study, two pear varieties with contrasting organic acid levels, ‘Dangshansuli’ (low-acidity) and ‘Amute’ (high-acidity), were selected, and a combination of transcriptome and population genomics analyses were applied to characterize their patterns of gene expression and genetic variation. Based on RNA-seq data analysis, differentially expressed genes (DEGs) involved in organic acid metabolism and accumulation were identified. Weighted correlation network analysis (WGCNA) revealed that nine candidate TCA (tricarboxylic acid)-related DEGs and three acid transporter-related DEGs were located in three key modules. The regulatory networks of the above candidate genes were also predicted. By integrating pear resequencing data, two domestication-related genes were found to be upregulated in ‘Amute’, and this trend was further validated for other pear varieties with high levels of organic acid, suggesting distinct selective sweeps during pear dissemination and domestication. Collectively, this study provides insight into organic acid differences related to expression divergence and domestication in two pear varieties, pinpointing several candidate genes for the genetic manipulation of acidity in pears.

The grapevine (Vitis vinifera L.) floral transcriptome in Pinot noir variety: identification of tissue-related gene networks and whorl-specific markers in pre- and post-anthesis phases

The comprehension of molecular processes underlying the development and progression of flowering in plants is a hot topic, not only because that often the products of interest for human and animal nutrition are linked to the development of fruits or seeds, but also because the processes of gametes formation occurring in sexual organs are at the basis of recombination and genetic variability which constitutes the matter on which evolution acts, whether understood as natural or human driven. In the present study, we used an NGS approach to produce a grapevine flower transcriptome snapshot in different whorls and tissues including calyx, calyptra, filament, anther, stigma, ovary, and embryo in both pre- and post-anthesis phases. Our investigation aimed at identifying hub genes that unequivocally distinguish the different tissues providing insights into the molecular mechanisms that are at the basis of floral whorls and tissue development. To this end we have used different analytical approaches, some now consolidated in transcriptomic studies on plants, such as pairwise comparison and weighted-gene coexpression network analysis, others used mainly in studies on animals or human's genomics, such as the tau (τ) analysis aimed at isolating highly and absolutely tissue-specific genes. The intersection of data obtained by these analyses allowed us to gradually narrow the field, providing evidence about the molecular mechanisms occurring in those whorls directly involved in reproductive processes, such as anther and stigma, and giving insights into the role of other whorls not directly related to reproduction, such as calyptra and calyx. We believe this work could represent an important genomic resource for functional analyses of grapevine floral organ growth and fruit development shading light on molecular networks underlying grapevine reproductive organ determination.

Identifying genomic regions targeted during eggplant domestication using transcriptome data

Identifying genes and traits that have diverged during domestication provides key information of importance for maintaining and even increasing yield and nutrients in existing crops. A “bottom-up” population genetics approach was used to identify signatures of selection across the eggplant genome, to better understand the process of domestication. RNA-seq data were obtained for 4 wild eggplants (Solanum insanum L.) and 16 domesticated eggplants (S. melongena L.) and mapped to the eggplant genome. Single-nucleotide polymorphism (SNPs) exhibiting signatures of selection in domesticates were identified as those exhibiting high F_ST between the 2 populations (evidence of significant divergence) and low π for the domesticated population (indicative of a selective sweep). Some of these regions appear to overlap with previously identified quantitative trait loci for domestication traits. Genes in regions of linkage disequilibrium surrounding these SNPs were searched against the Arabidopsis thaliana and tomato genomes to find orthologs. Subsequent gene ontology (GO) enrichment analysis identified over-representation of GO terms related to photosynthesis and response to the environment. This work reveals genomic changes involved in eggplant domestication and improvement, and how this compares to observed changes in the tomato genome, revealing shared chromosomal regions involved in the domestication of both species.

Comparative transcriptome analyses between cultivated and wild grapes reveal conservation of expressed genes but extensive rewiring of co-expression networks

The transcriptomes of wild and cultivated grapes consists of similar expressed genes but distinct wiring of co-expressed modules associated with environmental conditions. Grapevine is an important fruit crop worldwide, with high economic value and widespread distribution. Commercial production is based on Vitis vinifera, and, to a lesser extent, on hybrids with American grapes, such as V. labrusca. Wild grape relatives are important sources of resistance against biotic and abiotic factors; however, their global gene expression patterns remain poorly characterized. We associated genome-wide transcript profiling to phenotypic analyses to investigate the responses of cultivated and wild vines to vineyard conditions. The expressed genes in the Vitis reference transcriptome are largely shared by wild grapes, V. labrusca hybrids and vinifera cultivars. In contrast, significant differential regulation between wild and vinifera genotypes represents 80% of gene expression variation, regardless of the environment. In wild grapes, genes associated to regulatory processes are downregulated, whereas those involved in metabolic pathways are upregulated, in comparison to vinifera. Photosynthesis-related ontologies are overrepresented in the induced genes, in agreement with higher contents of chlorophyll in wild grapes. Co-regulated gene network analyses provide evidence of more complex transcriptome organization in vinifera. In wild grapes, genes involved in signaling pathways of stress-related hormones are overrepresented in modules associated with the environment. Consensus network analyses revealed high preservation within co-regulated gene modules between cultivated and wild grapes, but divergent relationships among the expression clusters. In conclusion, the distinct phenotypes of wild and cultivated grapes are underlain by differences in gene expression, but also by distinct higher-order organization of the transcriptome and contrasting association of co-expressed gene clusters with the environment.

Using RNA-seq to characterize pollen-stigma interactions for pollination studies

Insects are essential for the reproduction of pollinator-dependent crops and contribute to the pollination of 87% of wild plants and 75% of the world’s food crops. Understanding pollen flow dynamics between plants and pollinators is thus essential to manage and conserve wild plants and ensure yields are maximized in food crops. However, the determination of pollen transfer in the field is complex and laborious. We developed a field experiment in a pollinator-dependent crop and used high throughput RNA sequencing (RNA-seq) to quantify pollen flow by measuring changes in gene expression between pollination treatments across different apple (Malus domestica Borkh.) cultivars. We tested three potential molecular indicators of successful pollination and validated these results with field data by observing single and multiple visits by honey bees (Apis mellifera) to apple flowers and measured fruit set in a commercial apple orchard. The first indicator of successful outcrossing was revealed via differential gene expression in the cross-pollination treatments after 6 h. The second indicator of successful outcrossing was revealed by the expression of specific genes related to pollen tube formation and defense response at three different time intervals in the stigma and the style following cross-pollination (i.e. after 6, 24, and 48 h). Finally, genotyping variants specific to donor pollen could be detected in cross-pollination treatments, providing a third indicator of successful outcrossing. Field data indicated that one or five flower visits by honey bees were insufficient and at least 10 honey bee flower visits were required to achieve a 25% probability of fruit set under orchard conditions. By combining the genotyping data, the differential expression analysis, and the traditional fruit set field experiments, it was possible to evaluate the pollination effectiveness of honey bee visits under orchards conditions. This is the first time that pollen-stigma-style mRNA expression analysis has been conducted after a pollinator visit (honey bee) to a plant (in vivo apple flowers). This study provides evidence that mRNA sequencing can be used to address complex questions related to stigma–pollen interactions over time in pollination ecology.

Genome-wide association studies provide insights into the genetic determination of fruit traits of pear

Pear is a major fruit tree crop distributed worldwide, yet its breeding is a very time-consuming process. To facilitate molecular breeding and gene identification, here we have performed genome-wide association studies (GWAS) on eleven fruit traits. We identify 37 loci associated with eight fruit quality traits and five loci associated with three fruit phenological traits. Scans for selective sweeps indicate that traits including fruit stone cell content, organic acid and sugar contents might have been under continuous selection during breeding improvement. One candidate gene, PbrSTONE, identified in GWAS, has been functionally verified to be involved in the regulation of stone cell formation, one of the most important fruit quality traits in pear. Our study provides insights into the complex fruit related biology and identifies genes controlling important traits in pear through GWAS, which extends the genetic resources and basis for facilitating molecular breeding in perennial trees.

Contrasting genetic variation and positive selection followed the divergence of NBS-encoding genes in Asian and European pears

BACKGROUND

The NBS disease-related gene family coordinates the inherent immune system in plants in response to pathogen infections. Previous studies have identified NBS-encoding genes in Pyrus bretschneideri ('Dangshansuli', an Asian pear) and Pyrus communis ('Bartlett', a European pear) genomes, but the patterns of genetic variation and selection pressure on these genes during pear domestication have remained unsolved.

RESULTS

In this study, 338 and 412 NBS-encoding genes were identified from Asian and European pear genomes. This difference between the two pear species was the result of proximal duplications. About 15.79% orthologous gene pairs had Ka/Ks ratio more than one, indicating two pear species undergo strong positive selection after the divergence of Asian and European pear. We identified 21 and 15 NBS-encoding genes under fire blight and black spot disease-related QTL, respectively, suggesting their importance in disease resistance. Domestication caused decreased nucleotide diversity across NBS genes in Asian cultivars (cultivated 6.23E-03; wild 6.47E-03), but opposite trend (cultivated 6.48E-03; wild 5.91E-03) appeared in European pears. Many NBS-encoding coding regions showed Ka/Ks ratio of greater than 1, indicating the role of positive selection in shaping diversity of NBS-encoding genes in pear. Furthermore, we detected 295 and 122 significantly different SNPs between wild and domesticated accessions in Asian and European pear populations. Two NBS genes (Pbr025269.1 and Pbr019876.1) with significantly different SNPs showed >5x upregulation between wild and cultivated pear accessions, and > 2x upregulation in Pyrus calleryana after inoculation with Alternaria alternata. We propose that positively selected and significantly different SNPs of an NBS-encoding gene (Pbr025269.1) regulate gene expression differences in the wild and cultivated groups, which may affect resistance in pear against A. alternata.

CONCLUSION

Proximal duplication mainly led to the different number of NBS-encoding genes in P. bretschneideri and P. communis genomes. The patterns of genetic diversity and positive selection pressure differed between Asian and European pear populations, most likely due to their independent domestication events. This analysis helps us understand the evolution, diversity, and selection pressure in the NBS-encoding gene family in Asian and European populations, and provides opportunities to study mechanisms of disease resistance in pear.

Whole-genome resequencing reveals loci under selection during silkworm improvement

Breeding or genetic improvement refers to the process of artificial selection following domestication; as such, it has had a major influence on modern agriculture and animal production. Improvement generally focuses on traits that greatly affect the economic performance. Therefore, understanding the genetic basis underlying improvement will contribute to the identification of genes controlling economic traits and will facilitate future crop and animal breeding. However, genome-wide study of the molecular basis underlying improvement remains rare. The silkworm is a unique, entirely domesticated economically important invertebrate; genetic improvement has had a huge effect on the silkworm regarding silk-related traits. Herein, we performed whole-genomic sequencing on local and genetically improved silkworm lines to identify the genomic regions under strong selection in silkworm breeding/improvement. By genomic-wide selective sweeping analysis, we identified 24 genomic regions with strong selection signals, eight of which contained 13 candidate genes underlying silkworm breeding. Interestingly, six of these genes were annotated with functions related to neural signal response. Among the six genes, BGIBMGA004050 encodes silkworm CREB-regulated_transcription_coactivator_1 (BmCRTC1), which was reported to be involved in energy-sensing pathways. These results suggested that improvement may have affected the nervous system of the silkworm. This research will provide new insights into the genetic basis underlying the genetic improvement of silkworms and possibly of other species.

Cis- and Trans-Regulatory Variations in the Domestication of the Chili Pepper Fruit

The process of domestication requires the rapid transformation of the wild morphology into the cultivated forms that humans select for. This process often takes place through changes in the regulation of genes, yet, there is no definite pattern on the role of cis- and trans-acting regulatory variations in the domestication of the fruit among crops. Using allele-specific expression and network analyses, we characterized the regulatory patterns and the inheritance of gene expression in wild and cultivated accessions of chili pepper, a crop with remarkable fruit morphological variation. We propose that gene expression differences associated to the cultivated form are best explained by cis-regulatory hubs acting through trans-regulatory cascades. We show that in cultivated chili, the expression of genes associated with fruit morphology is partially recessive with respect to those in the wild relative, consistent with the hybrid fruit phenotype. Decreased expression of fruit maturation and growth genes in cultivated chili suggest that selection for loss-of-function took place in its domestication. Trans-regulatory changes underlie the majority of the genes showing regulatory divergence and had larger effect sizes on gene expression than cis-regulatory variants. Network analysis of selected cis-regulated genes, including ARP9 and MED25, indicated their interaction with many transcription factors involved in organ growth and fruit ripening. Differentially expressed genes linked to cis-regulatory variants and their interactions with downstream trans-acting genes have the potential to drive the morphological differences observed between wild and cultivated fruits and provide an attractive mechanism of morphological transformation during the domestication of the chili pepper.

Effects of Metaxenia on Stone Cell Formation in Pear (Pyrus bretschneideri) Based on Transcriptomic Analysis and Functional Characterization of the Lignin-Related Gene PbC4H2

The deposition of lignin in flesh parenchyma cells for pear stone cells, and excessive stone cells reduce the taste and quality of the fruit. The effect of metaxenia on the quality of fruit has been heavily studied, but the effect of metaxenia on stone cell formation has not been fully elucidated to date. This study used P. bretschneideri (Chinese white pear) cv. ‘Yali’ (high-stone cell content) and P. pyrifolia (Sand pear) cv. ‘Cuiguan’ (low-stone cell content) as pollination trees to pollinate P. bretschneideri cv. ‘Lianglizaosu’ separately to fill this gap in the literature. The results of quantitative determination, histochemical staining and electron microscopy indicated that the content of stone cells and lignin in YL fruit (‘Yali’ (pollen parent) × ‘Lianglizaosu’ (seed parent)) was significantly higher than that in CL fruit (‘Cuiguan’ (pollen parent) × ‘Lianglizaosu’ (seed parent)). The transcriptome sequencing results that were obtained from the three developmental stages of the two types of hybrid fruits indicated that a large number of differentially expressed genes (DEGs) related to auxin signal transduction (AUX/IAAs and ARFs), lignin biosynthesis, and lignin metabolism regulation (MYBs, LIMs, and KNOXs) between the CL and YL fruits at the early stage of fruit development. Therefore, metaxenia might change the signal transduction process of auxin in pear fruit, thereby regulating the expression of transcription factors (TFs) related to lignin metabolism, and ultimately affecting lignin deposition and stone cell development. In addition, we performed functional verification of a differentially expressed gene, PbC4H2 (cinnamate 4-hydroxylase). Heterologous expression of PbC4H2 in the c4h mutant not only restored its collapsed cell wall, but also significantly increased the lignin content in the inflorescence stem. The results of our research help to elucidate the metaxenia-mediated regulation of pear stone cell development and clarify the function of PbC4H2 in cell wall development and lignin synthesis, which establishes a foundation for subsequent molecular breeding.

Transcriptomic analysis of early fruit development in Chinese white pear (Pyrus bretschneideri Rehd.) and functional identification of PbCCR1 in lignin biosynthesis

BACKGROUND

The content of stone cells and lignin is one of the key factors affecting the quality of pear fruit. In a previous study, we determined the developmental regularity of stone cells and lignin in 'Dangshan Su' pear fruit 15-145 days after pollination (DAP). However, the development of fruit stone cells and lignin before 15 DAP has not been heavily researched.

RESULTS

In this study, we found that primordial stone cells began to appear at 7 DAP and that the fruit had formed a large number of stone cells at 15 DAP. Subsequently, transcriptome sequencing was performed on fruits at 0, 7, and 15 DAP and identified 3834 (0 vs. 7 DAP), 4049 (7 vs. 15 DAP) and 5763 (0 vs. 15 DAP) DEGs. During the 7-15 DAP period, a large number of key enzyme genes essential for lignin biosynthesis are gradually up-regulated, and their expression pattern is consistent with the accumulation of lignin in this period. Further analysis found that the biosynthesis of S-type lignin in 'Dangshan Su' pear does not depend on the catalytic activity of PbSAD but is primarily generated by the catalytic activity of caffeoyl-CoA through CCoAOMT, CCR, F5H, and CAD. We cloned PbCCR1, 2 and analysed their functions in Chinese white pear lignin biosynthesis. PbCCR1 and 2 have a degree of functional redundancy; both demonstrate the ability to participate in lignin biosynthesis. However, PbCCR1 may be the major gene for lignin biosynthesis, while PbCCR2 has little effect on lignin biosynthesis.

CONCLUSIONS

Our results revealed that 'Dangshan Su' pear began to form a large number of stone cells and produce lignin after 7 DAP and mainly accumulated materials from 0 to 7 DAP. PbCCR1 is mainly involved in the biosynthesis of lignin in 'Dangshan Su' pear and plays a positive role in lignin biosynthesis.