Pseudo-chromosome-length genome assembly of a double haploid "Bartlett" pear (Pyrus communis L.)

Gap-free genome assemblies of two Pyrus bretschneideri cultivars and GWAS analyses identify a CCCH zinc finger protein as a key regulator of stone cell formation in pear fruit

The Chinese white pear (Pyrus bretschneideri) is an economically significant fruit crop worldwide. Previous versions of the P. bretschneideri genome assembly contain numerous gaps and unanchored genetic regions. Here, we generated two high-quality, gap-free genome assemblies for 'Dangshansu' (DS; 503.92 Mb) and 'Lianglizaosu' (ZS; 509.01 Mb), each anchored to 17 chromosomes, achieving a benchmarking universal single-copy ortholog completeness score of nearly 99.0%. Our genome-wide association studies explored the associations between genetic variations and stone cell traits, revealing a significant association peak on DS chromosome 3 and identifying a novel non-tandem CCCH-type zinc finger gene, designated PbdsZF. Through genetic transformation, we verified the pivotal role of PbdsZF in regulation of both lignin biosynthesis and stone cell formation, as it transcriptionally activates multiple genes involved in these processes. By binding to the CT-rich motifs CT1 (CTTTTTTCT) and CT2 (CTCTTTTT), PbdsZF significantly influences the transcription of genes essential for lignin production, underscoring its regulatory importance in plant lignin metabolism. Our study illuminates the complex biology of fruit development and delineates the gene regulatory networks that influence stone cell and lignocellulose formation, thereby enriching genetic resources and laying the groundwork for the molecular breeding of perennial trees.

Long-read genome sequencing reveals the sequence characteristics of pear self-incompatibility locus

The S-RNase-based self-incompatibility locus (S-locus) in Petunia species contains 16-20 F-box genes, which collaboratively function in the recognition and subsequent degradation of non-self S-RNases, while distinguishing them from self S-RNase. However, the number of S-locus F-box genes (SFBBs) physically interacted with non-self S-RNases remains uncertain in Pyrus species. Utilizing Pacbio long-read sequencing, we successfully assembled the genome of pear cultivar 'Yali' (Pyrus bretschneideri), and identified 19 SFBBs from the Pyrus S17-locus spanning approximately 1.78 Mb. Additionally, we identified 17-21 SFBBs from other Pyrus and Malus S-loci spanning a range of 1.35 to 2.64 Mb. Based on the phylogenetic analysis, it was determined that Pyrus and Malus SFBBs could be classified into 22 groups, denoted as I to XXII. At amino acid level, SFBBs within a given group exhibited average identities ranged from 88.9% to 97.9%. Notably, all 19 SFBBs from the S17-locus co-segregated with S17-RNase, with 18 of them being specifically expressed in pollen. Consequently, these 18 pollen-specifically expressed SFBBs are considered potential candidates for the pollen-S determinant. Intriguingly, out of the 18 pollen-specifically expressed SFBBs, eight demonstrated interactions with at least one non-self S-RNase, while the remaining SFBBs failed to recognize any S-RNase. These findings provide compelling evidence supporting the existence of a collaborative non-self-recognition system governing self-incompatibility in pear species.

Comprehensive insights on association mapping in perennial fruit crops breeding - Its implications, current status and future perspectives

In order to provide food and nutritional security for the world's rapidly expanding population, fruit crop researchers have identified two critical priorities: increasing production and preserving fruit quality during the pre- and post-harvest periods. The genetic basis of these complex, commercially important fruit traits which are uniquely regulated by polygenes or multi-allelic genes that interact with one another and the environment can be analyzed with the aid of trait mapping tools. The most interesting trait mapping approach that offers the genetic level investigation for marker-trait associations (MTAs) for these complex fruit traits, without the development of mapping population, is association mapping. This approach was used during the genetic improvement program, emphasizing the obstacles (breeding strategies adopted, generation interval, and their genomic status) pertaining to perennial fruit crops. This method of studying population diversity and linkage disequilibrium in perennial fruit crops has been made possible by recent developments in genotyping, phenotyping, and statistical analysis. Thus, the purpose of this review is to provide an overview of different trait mapping techniques, with a focus on association mapping (method, essential components, viability, constraints, and future perspective) and its advantages, disadvantages, and possibilities for breeding perennial fruit crops.

The genome sequence of the pear, Pyrus communis L

We present a genome assembly from a specimen of Pyrus communis (the pear; Streptophyta; Magnoliopsida; Rosales; Rosaceae). The genome sequence has a total length of 487.30 megabases. Most of the assembly is scaffolded into 17 chromosomal pseudomolecules. The mitochondrial and plastid genome assemblies have lengths of 443.53 kilobases and 159.93 kilobases, respectively. Gene annotation of this assembly on Ensembl identified 37,713 protein-coding genes.

Genome-wide identification of the Dof gene family: How it plays a part in mediating cold stress response in Prunus mume

DNA binding with a finger transcription factor (Dof) takes part in several plant physiological activities such as seed germination, flowering time, cold and drought resistance. Although the function, molecular phylogeny and expression pattern of Dof genes in Prunus mume was not clear yet. Here, the gene structure, motif, chromosome location and phylogenetic relationship of the Dof gene family in Prunus species was explored. We identified 24 members of the Dof gene family from P. mume, which were divided into 3 different subgroups. All these PmDof genes can be mapped to the pseudochromosome. Only one pair of tandem duplication genes are located in Chr3, whereas 8 pairs of segmentally duplicated PmDof genes located in Chr1, Chr2, Chr4, Chr5, and Chr7. Motif and gene structure analysis showed that each group had a similar conservative motif and similar exon/intron composition. Cis-acting elements analysis indicate that PmDofs may be involved in regulating abiotic stress response. Gene expression patterns showed that most PmDofs genes were specifically expressed in different tissues and at different stages. We next found that PmDofs genes display an obvious expression preference or specificity in cold stress response according to qRT-PCR analysis. We further observe a great cold resistance in PmDof10/11/20 OE lines, they showed lower electrolyte leakage rate, MDA content and higher soluble sugar/protein, POD/SOD/proline content than WT after -5 °C 6h freezing treatment. This research offers fresh perspectives on the development of PmDofs, enhancing our comprehension of the structure and role of plant Dof gene families.

Telomere-to-telomere Genome Assembly of two representative Asian and European pear cultivars

As the third most important temperate fruit, Pear (Pyrus spp.) exhibits a remarkable genetic diversity and is classified into two mainly categories known as Asian pear and European pear. Although several pear genomes are available, most of the released versions are fragmented and not chromosome-level high-quality. In this study, we report two high-quality genomes for Pyrus bretschneideri Rhed. cv. 'Danshansuli' (DS) and Pyrus communis L. cv. 'Conference' (KFL), which represent the predominant Asian and European cultivars, respectively, with nearly telomere-to-telomere (T2T) gap-free level. The finally assembled genome sizes for DS and KFL were 510.98 Mb and 510.71 Mb, respectively, with Contig N50 of 29.47 Mb and 30.47 Mb, where each chromosome was represented by a single contig. The DS and KFL genomes yielded a total of 46,394 and 44,702 protein-coding genes, respectively. Among these genes, the functional annotation accounted for 96.47% and 96.46% in the DS and KFL genomes. The two novels nearly T2T genomic information offers an invaluable resource for comparative genomics, genetic diversity analysis, molecular breeding strategies, and functional exploration.

Haplotype-resolved T2T genome assemblies and pangenome graph of pear reveal diverse patterns of allele-specific expression and the genomic basis of fruit quality traits

Hybrid crops often exhibit increased yield and greater resilience, yet the genomic mechanism(s) underlying hybrid vigor or heterosis remain unclear, hindering our ability to predict the expression of phenotypic traits in hybrid breeding. Here, we generated haplotype-resolved T2T genome assemblies of two pear hybrid varieties, 'Yuluxiang' (YLX) and 'Hongxiangsu' (HXS), which share the same maternal parent but differ in their paternal parents. We then used these assemblies to explore the genome-scale landscape of allele-specific expression (ASE) and create a pangenome graph for pear. ASE was observed for close to 6000 genes in both hybrid cultivars. A subset of ASE genes related to aspects of fruit quality such as sugars, organic acids, and cuticular wax were identified, suggesting their important contributions to heterosis. Specifically, Ma1, a gene regulating fruit acidity, is absent in the paternal haplotypes of HXS and YLX. A pangenome graph was built based on our assemblies and seven published pear genomes. Resequencing data for 139 cultivated pear genotypes (including 97 genotypes sequenced here) were subsequently aligned to the pangenome graph, revealing numerous structural variant hotspots and selective sweeps during pear diversification. As predicted, the Ma1 allele was found to be absent in varieties with low organic acid content, and this association was functionally validated by Ma1 overexpression in pear fruit and calli. Overall, these results reveal the contributions of ASE to fruit-quality heterosis and provide a robust pangenome reference for high-resolution allele discovery and association mapping.

Advancing apple genetics research: Malus coronaria and Malus ioensis genomes and a gene family-based pangenome of native North American apples

Wild Malus species flourished in North America long before Europeans introduced domesticated apples. Malus coronaria and M. ioensis are native to the mid-western and eastern United States, while M. angustifolia and M. fusca grow in the southeast and west, respectively. They offer disease resistance, climate and soil adaptability, and horticultural traits for apple breeding. However, their utilization remains limited due to insufficient genomic resources and specific genetics. We report high-quality phased chromosome-scale assemblies of M. coronaria and M. ioensis, generated using long-read and conformation capture sequencing. Phylogenetic and synteny analysis indicated high relatedness between these 2 genomes and previously published genome of M. angustifolia, and lower relatedness with M. fusca. Gene family-based pangenome of North American Malus identified 60,211 orthogroups containing 340,087 genes. Genes involved in basic cellular and metabolic processes, growth, and development were core to the existence of these species, whereas genes involved in secondary metabolism, stress response, and interactions with other organisms were accessory and are likely associated with adaptation to specific environments. Structural variation hotspots were mostly overlapping with high gene density. This study offers novel native North American Malus genome resources that can be used to identify genes for apple breeding and understand their evolution and adaptation.

Genomic data provides insights into the evolutionary history and adaptive differentiation of two tetraploid strawberries

Over the decades, evolutionists and ecologists have shown intense interest in the role of polyploidization in plant evolution. Without clear knowledge of the diploid ancestor(s) of polyploids, we would not be able to answer fundamental ecological questions such as the evolution of niche differences between them or its underlying genetic basis. Here, we explored the evolutionary history of two Fragaria tetraploids, Fragaria corymbosa and Fragaria moupinensis. We de novo assembled five genomes including these two tetraploids and three diploid relatives. Based on multiple lines of evidence, we found no evidence of subgenomes in either of the two tetraploids, suggesting autopolyploid origins. We determined that Fragaria chinensis was the diploid ancestor of F. corymbosa while either an extinct species affinitive to F. chinensis or an unsampled population of F. chinensis could be the progenitor of F. moupinensis. Meanwhile, we found introgression signals between F. chinensis and Fragaria pentaphylla, leading to the genomic similarity between these two diploids. Compared to F. chinensis, gene families related to high ultraviolet (UV)-B and DNA repair were expanded, while those that responded towards abiotic and biotic stresses (such as salt stress, wounding, and various pathogens) were contracted in both tetraploids. Furthermore, the two tetraploids tended to down-regulate defense response genes but up-regulate UV-B response, DNA repairing, and cell division gene expression compared to F. chinensis. These findings may reflect adaptions toward high-altitude habitats. In summary, our work provides insights into the genome evolution of wild Fragaria tetraploids and opens up an avenue for future works to answer deeper evolutionary and ecological questions regarding the strawberry genus.

Four glycosyltransferase genes are responsible for synthesis and accumulation of different flavonol glycosides in apple tissues

Flavonols are widely synthesized throughout the plant kingdom, playing essential roles in plant physiology and providing unique health benefits for humans. Their glycosylation plays significant role in improving their stability and solubility, thus their accumulation and function. However, the genes encoding the enzymes catalyze this glycosylation remain largely unknown in apple. This study utilized a combination of methods to identify genes encoding such enzymes. Initially, candidate genes were selected based on their potential to encode UDP-dependent glycosyltransferases (UGTs) and their expression patterns in response to light induction. Subsequently, through testing the in vitro enzyme activity of the proteins produced in Escherichia coli cells, four candidates were confirmed to encode a flavonol 3-O-galactosyltransferase (UGT78T6), flavonol 3-O-glucosyltransferase (UGT78S1), flavonol 3-O-xylosyltransferase/arabinosyltransferase (UGT78T5), and flavonol 3-O-rhamnosyltransferase (UGT76AE22), respectively. Further validation of these genes' functions was conducted by modulating their expression levels in stably transformed apple plants. As anticipated, a positive correlation was observed between the expression levels of these genes and the content of specific flavonol glycosides corresponding to each gene. Moreover, overexpression of a flavonol synthase gene, MdFLS, resulted in increased flavonol glycoside content in apple roots and leaves. These findings provide valuable insights for breeding programs aimed at enriching apple flesh with flavonols and for identifying flavonol 3-O-glycosyltransferases of other plant species.

Haplotype-resolved genome assembly for tetraploid Chinese cherry (Prunus pseudocerasus) offers insights into fruit firmness

Chinese cherry (Prunus pseudocerasus) holds considerable importance as one of the primary stone fruit crops in China. However, artificially improving its traits and genetic analysis are challenging due to lack of high-quality genomic resources, which mainly result from difficulties associated with resolving its tetraploid and highly heterozygous genome. Herein, we assembled a chromosome-level, haplotype-resolved genome of the cultivar 'Zhuji Duanbing', comprising 993.69 Mb assembled into 32 pseudochromosomes using PacBio HiFi, Oxford Nanopore, and Hi-C. Intra-haplotype comparative analyses revealed extensive intra-genomic sequence and expression consistency. Phylogenetic and comparative genomic analyses demonstrated that P. pseudocerasus was a stable autotetraploid species, closely related to wild P. pusilliflora, with the two diverging ~18.34 million years ago. Similar to other Prunus species, P. pseudocerasus underwent a common whole-genome duplication event that occurred ~139.96 million years ago. Because of its low fruit firmness, P. pseudocerasus is unsuitable for long-distance transportation, thereby restricting its rapid development throughout China. At the ripe fruit stage, P. pseudocerasus cv. 'Zhuji Duanbing' was significantly less firm than P. avium cv. 'Heizhenzhu'. The difference in firmness is attributed to the degree of alteration in pectin, cellulose, and hemicellulose contents. In addition, comparative transcriptomic analyses identified GalAK-like and Stv1, two genes involved in pectin biosynthesis, which potentially caused the difference in firmness between 'Zhuji Duanbing' and 'Heizhenzhu'. Transient transformations of PpsGalAK-like and PpsStv1 increase protopectin content and thereby enhance fruit firmness. Our study lays a solid foundation for functional genomic studies and the enhancement of important horticultural traits in Chinese cherries.

Assessment and Partial Characterization of Candidate Genes in Dihydrochalcone and Arbutin Biosynthesis in an Apple-Pear Hybrid by De Novo Transcriptome Assembly

Apples (Malus × domestica Borkh.) and pears (Pyrus communis L.) are valuable crops closely related within the Rosaceae family with reported nutraceutical properties derived from secondary metabolites including phloridzin and arbutin, which are distinctive phenolic metabolites characterizing apples and pears, respectively. Here, we generated a de novo transcriptome assembly of an intergeneric hybrid between apple and pear, accumulating intermediate levels of phloridzin and arbutin. Combining RNA-seq, in silico functional annotation prediction, targeted gene expression analysis, and expression-metabolite correlations, we identified candidate genes for functional characterization, resulting in the identification of active arbutin synthases in the hybrid and parental genotypes. Despite exhibiting an active arbutin synthase in vitro, the natural lack of arbutin in apples is reasoned by the absence of the substrate and broad substrate specificity. Altogether, our study serves as the basis for future assessment of potential physiological roles of identified genes by genome editing of hybrids and pears.

Haplotype-resolved genome of Prunus zhengheensis provides insight into its evolution and low temperature adaptation in apricot

Prunus zhengheensis, an extremely rare population of apricots, originated in warm South-East China and is an excellent material for genetic breeding. However, most apricots and two related species (P. sibirica, P. mandshurica) are found in the cold northern regions in China and the mechanism of their distribution is still unclear. In addition, the classification status of P. zhengheensis is controversial. Thus, we generated a high-quality haplotype-resolved genome for P. zhengheensis, exploring key genetic variations in its adaptation and the causes of phylogenetic incongruence. We found extensive phylogenetic discordances between the nuclear and organelle phylogenies of P. zhengheensis, which could be explained by incomplete lineage sorting. A 242.22-Mb pan-genome of the Armeniaca section was developed with 13 chromosomal genomes. Importantly, we identified a 566-bp insertion in the promoter of the HSFA1d gene in apricot and showed that the activity of the HSFA1d promoter increased under low temperatures. In addition, HSFA1d overexpression in Arabidopsis thaliana indicated that HSFA1d positively regulated plant growth under chilling. Therefore, we hypothesized that the insertion in the promoter of HSFA1d in apricot improved its low-temperature adaptation, allowing it to thrive in relatively cold locations. The findings help explain the weather adaptability of Armeniaca plants.

Multi-omics provide insights into the regulation of DNA methylation in pear fruit metabolism

BACKGROUND

Extensive research has been conducted on fruit development in crops, but the metabolic regulatory networks underlying perennial fruit trees remain poorly understood. To address this knowledge gap, we conduct a comprehensive analysis of the metabolome, proteome, transcriptome, DNA methylome, and small RNAome profiles of pear fruit flesh at 11 developing stages, spanning from fruitlet to ripening. Here, we systematically investigate the metabolic landscape and regulatory network involved.

RESULTS

We generate an association database consisting of 439 metabolites and 14,399 genes to elucidate the gene regulatory network of pear flesh metabolism. Interestingly, we detect increased DNA methylation in the promoters of most genes within the database during pear flesh development. Application of a DNA methylation inhibitor to the developing fruit represses chlorophyll degradation in the pericarp and promotes xanthophyll, β-carotene, and abscisic acid (ABA) accumulation in the flesh. We find the gradual increase in ABA production during pear flesh development is correlated with the expression of several carotenoid pathway genes and multiple transcription factors. Of these transcription factors, the zinc finger protein PbZFP1 is identified as a positive mediator of ABA biosynthesis in pear flesh. Most ABA pathway genes and transcription factors are modified by DNA methylation in the promoters, although some are induced by the DNA methylation inhibitor. These results suggest that DNA methylation inhibits ABA accumulation, which may delay fruit ripening.

CONCLUSION

Our findings provide insights into epigenetic regulation of metabolic regulatory networks during pear flesh development, particularly with regard to DNA methylation.

Comparative transcriptome profiling and co-expression network analysis uncover the key genes associated with pear petal defense responses against Monilinia laxa infection

Pear brown rot and blossom blight caused by Monilinia laxa seriously affect pear production worldwide. Here, we compared the transcriptomic profiles of petals after inoculation with M. laxa using two pear cultivars with different levels of sensitivity to disease (Sissy, a relatively tolerant cultivar, and Kristalli, a highly susceptible cultivar). Physiological indexes were also monitored in the petals of both cultivars at 2 h and 48 h after infection (2 HAI and 48 HAI). RNA-seq data and weighted gene co-expression network analysis (WGCNA) allowed the identification of key genes and pathways involved in immune- and defense-related responses that were specific for each cultivar in a time-dependent manner. In particular, in the Kristalli cultivar, a significant transcriptome reprogramming occurred early at 2 HAI and was accompanied either by suppression of key differentially expressed genes (DEGs) involved in the modulation of any defense responses or by activation of DEGs acting as sensitivity factors promoting susceptibility. In contrast to the considerably high number of DEGs induced early in the Kristalli cultivar, upregulation of specific DEGs involved in pathogen perception and signal transduction, biosynthesis of secondary and primary metabolism, and other defense-related responses was delayed in the Sissy cultivar, occurring at 48 HAI. The WGCNA highlighted one module that was significantly and highly correlated to the relatively tolerant cultivar. Six hub genes were identified within this module, including three WRKY transcription factor-encoding genes: WRKY 65 (pycom05g27470), WRKY 71 (pycom10g22220), and WRKY28 (pycom17g13130), which may play a crucial role in enhancing the tolerance of pear petals to M. laxa. Our results will provide insights into the interplay of the molecular mechanisms underlying immune responses of petals at the pear-M. laxa pathosystem.

A chromosome-scale assembly for 'd'Anjou' pear

Cultivated pear consists of several Pyrus species with Pyrus communis (European pear) representing a large fraction of worldwide production. As a relatively recently domesticated crop and perennial tree, pear can benefit from genome-assisted breeding. Additionally, comparative genomics within Rosaceae promises greater understanding of evolution within this economically important family. Here, we generate a fully phased chromosome-scale genome assembly of P. communis 'd'Anjou.' Using PacBio HiFi and Dovetail Omni-C reads, the genome is resolved into the expected 17 chromosomes, with each haplotype totaling nearly 540 Megabases and a contig N50 of nearly 14 Mb. Both haplotypes are highly syntenic to each other and to the Malus domestica 'Honeycrisp' apple genome. Nearly 45,000 genes were annotated in each haplotype, over 90% of which have direct RNA-seq expression evidence. We detect signatures of the known whole-genome duplication shared between apple and pear, and we estimate 57% of d'Anjou genes are retained in duplicate derived from this event. This genome highlights the value of generating phased diploid assemblies for recovering the full allelic complement in highly heterozygous crop species.

A self-compatible pear mutant derived from γ-irradiated pollen carries an 11-Mb duplication in chromosome 17

Self-compatibility is a highly desirable trait for pear breeding programs. Our breeding program previously developed a novel self-compatible pollen-part Japanese pear mutant (Pyrus pyrifolia Nakai), '415-1', by using γ-irradiated pollen. '415-1' carries the S-genotype S4dS5S5, with "d" indicating a duplication of S 5 responsible for breakdown of self-incompatibility. Until now, the size and inheritance of the duplicated segment was undetermined, and a reliable detection method was lacking. Here, we examined genome duplications and their inheritance in 140 F1 seedlings resulting from a cross between '515-20' (S1S3) and '415-1'. Amplicon sequencing of S-RNase and SFBB18 clearly detected S-haplotype duplications in the seedlings. Intriguingly, 30 partially triploid seedlings including genotypes S1S4dS5, S3S4dS5, S1S5dS5, S3S5dS5, and S3S4dS4 were detected among the 140 seedlings. Depth-of-coverage analysis using ddRAD-seq showed that the duplications in those individuals were limited to chromosome 17. Further analysis through resequencing confirmed an 11-Mb chromosome duplication spanning the middle to the end of chromosome 17. The duplicated segment remained consistent in size across generations. The presence of an S3S4dS4 seedling provided evidence for recombination between the duplicated S5 segment and the original S4haplotype, suggesting that the duplicated segment can pair with other parts of chromosome 17. This research provides valuable insights for improving pear breeding programs using partially triploid individuals.

Introduction of a diverse genetic background of Pyrus into Malus through intergeneric hybridization

Wide hybridizations across species and genera have been employed to enhance agriculturally important traits in crops. Within the tribe Maleae of the Rosaceae family, different genera and species exhibit several traits useful for increasing diversity and gene pool through hybridization. This study aimed to develop and characterize intergeneric hybrid individuals between Malus and Pyrus. Through seed germination, shoot multiplication, and rooting in vitro, acclimatized seedlings showing vegetative growth on their own roots were obtained from crosses of Malus × domestica pollinated by Pyrus communis, P. bretschneideri, and the Pyrus interspecific hybrid (P. communis × P. pyrifolia). Comparative analysis of leaf morphology, flow cytometry, and molecular genotyping confirmed the hybrid status of the individuals. Genome-wide genotyping revealed that all the hybrid individuals inherited genomic fragments symmetrically from the Malus and Pyrus parents. To the best of our knowledge, this is the first report on the development of intergeneric hybrid seedlings between Malus × domestica and P. bretschneideri. Furthermore, the Pyrus interspecific hybrid individual served as a bridge plant for introducing the genetic background of P. pyrifolia into Malus × domestica. The results of this study provided a crucial foundation for breeding through intergeneric hybridization between Malus and Pyrus, facilitating the incorporation of valuable traits from diverse gene pools.

Near-gapless and haplotype-resolved apple genomes provide insights into the genetic basis of rootstock-induced dwarfing

Dwarfing rootstocks have transformed the production of cultivated apples; however, the genetic basis of rootstock-induced dwarfing remains largely unclear. We have assembled chromosome-level, near-gapless and haplotype-resolved genomes for the popular dwarfing rootstock 'M9', the semi-vigorous rootstock 'MM106' and 'Fuji', one of the most commonly grown apple cultivars. The apple orthologue of auxin response factor 3 (MdARF3) is in the Dw1 region of 'M9', the major locus for rootstock-induced dwarfing. Comparing 'M9' and 'MM106' genomes revealed a 9,723-bp allele-specific long terminal repeat retrotransposon/gypsy insertion, DwTE, located upstream of MdARF3. DwTE is cosegregated with the dwarfing trait in two segregating populations, suggesting its prospective utility in future dwarfing rootstock breeding. In addition, our pipeline discovered mobile mRNAs that may contribute to the development of dwarfed scion architecture. Our research provides valuable genomic resources and applicable methodology, which have the potential to accelerate breeding dwarfing rootstocks for apple and other perennial woody fruit trees.

Genome-wide identification of the Pyrus R2R3-MYB gene family and PhMYB62 regulation analysis in Pyrus hopeiensis flowers at low temperature

The R2R3-MYB gene family play an important role in plant growth, development and stress responses. In this study, a total of 122 PcoR2R3-MYB genes were identified and grouped into 26 clades in pear. And these PcoMYBs were unevenly distributed among 17 chromosomes. The sequence characteristics, conversed motifs, exon/intron structures, classification, duplication events and cis-acting elements were also investigated. The gene duplication events showed that segmental duplication may play key roles in expansion of the PcoMYB gene family. Pyrus hopeiensis, which is a valuable wild resource, has strong cold resistance. An integrative analyses of miRNA and mRNA showed that PhMYB62 was involved in regulating low-temperature stress in P. hopeiensis flower organs. Subcellular localization analysis showed that PhMYB62 protein was specifically localized to the nucleus. The result of DAP-seq showed that PhMYB62 responded to low-temperature stress in P. hopeiensis by regulating TFs, which were associated with plant stress resistance, and POD, GAUT12, AUX28 and CHS genes. Subsequently, yeast one-hybrid verified that PhMYB62 could bind and activate the promoter of POD gene. The current study would provide a comprehensive information for further functional research on the stress-responsive R2R3-MYB gene candidates in pear, and may help to identify the genes associated with cold resistance for the cultivation of cold-resistant pear varieties.

A haplotype resolved chromosome-scale assembly of North American wild apple Malus fusca and comparative genomics of the fire blight Mfu10 locus

The Pacific crabapple (Malus fusca) is a wild relative of the commercial apple (Malus × domestica). With a range extending from Alaska to Northern California, M. fusca is extremely hardy and disease resistant. The species represents an untapped genetic resource for the development of new apple cultivars with enhanced stress resistance. However, gene discovery and utilization of M. fusca have been hampered by the lack of genomic resources. Here, we present a high‐quality, haplotype‐resolved, chromosome‐scale genome assembly and annotation for M. fusca. The genome was assembled using high‐fidelity long‐reads and scaffolded using genetic maps and high‐throughput chromatin conformation capture sequencing, resulting in one of the most contiguous apple genomes to date. We annotated the genome using public transcriptomic data from the same species taken from diverse plant structures and developmental stages. Using this assembly, we explored haplotypic structural variation within the genome of M. fusca, identifying thousands of large variants. We further showed high sequence co‐linearity with other domesticated and wild Malus species. Finally, we resolve a known quantitative trait locus associated with resistance to fire blight (Erwinia amylovora). Insights gained from the assembly of a reference‐quality genome of this hardy wild apple relative will be invaluable as a tool to facilitate DNA‐informed introgression breeding.

Telomere-to-telomere pear (Pyrus pyrifolia) reference genome reveals segmental and whole genome duplication driving genome evolution

Previously released pear genomes contain a plethora of gaps and unanchored genetic regions. Here, we report a telomere-to-telomere (T2T) gap-free genome for the red-skinned pear, 'Yunhong No. 1' (YH1; Pyrus pyrifolia), which is mainly cultivated in Yunnan Province (southwest China), the pear's primary region of origin. The YH1 genome is 501.20 Mb long with a contig N50 length of 29.26 Mb. All 17 chromosomes were assembled to the T2T level with 34 characterized telomeres. The 17 centromeres were predicted and mainly consist of centromeric-specific monomers (CEN198) and long terminal repeat (LTR) Gypsy elements (≥74.73%). By filling all unclosed gaps, the integrity of YH1 is markedly improved over previous P. pyrifolia genomes ('Cuiguan' and 'Nijisseiki'). A total of 1531 segmental duplication (SD) driven duplicated genes were identified and enriched in stress response pathways. Intrachromosomal SDs drove the expansion of disease resistance genes, suggesting the potential of SDs in adaptive pear evolution. A large proportion of duplicated gene pairs exhibit dosage effects or sub-/neo-functionalization, which may affect agronomic traits like stone cell content, sugar content, and fruit skin russet. Furthermore, as core regulators of anthocyanin biosynthesis, we found that MYB10 and MYB114 underwent various gene duplication events. Multiple copies of MYB10 and MYB114 displayed obvious dosage effects, indicating role differentiation in the formation of red-skinned pear fruit. In summary, the T2T gap-free pear genome provides invaluable resources for genome evolution and functional genomics.

A telomere-to-telomere reference genome provides genetic insight into the pentacyclic triterpenoid biosynthesis in Chaenomeles speciosa

Chaenomeles speciosa (2n = 34), a medicinal and edible plant in the Rosaceae, is commonly used in traditional Chinese medicine. To date, the lack of genomic sequence and genetic studies has impeded efforts to improve its medicinal value. Herein, we report the use of an integrative approach involving PacBio HiFi (third-generation) sequencing and Hi-C scaffolding to assemble a high-quality telomere-to-telomere genome of C. speciosa. The genome comprised 650.4 Mb with a contig N50 of 35.5 Mb. Of these, 632.3 Mb were anchored to 17 pseudo-chromosomes, in which 12, 4, and 1 pseudo-chromosomes were represented by a single contig, two contigs, and four contigs, respectively. Eleven pseudo-chromosomes had telomere repeats at both ends, and four had telomere repeats at a single end. Repetitive sequences accounted for 49.5% of the genome, while a total of 45 515 protein-coding genes have been annotated. The genome size of C. speciosa was relatively similar to that of Malus domestica. Expanded or contracted gene families were identified and investigated for their association with different plant metabolisms or biological processes. In particular, functional annotation characterized gene families that were associated with the biosynthetic pathway of oleanolic and ursolic acids, two abundant pentacyclic triterpenoids in the fruits of C. speciosa. Taken together, this telomere-to-telomere and chromosome-level genome of C. speciosa not only provides a valuable resource to enhance understanding of the biosynthesis of medicinal compounds in tissues, but also promotes understanding of the evolution of the Rosaceae.

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.

Transcriptomic approach to uncover dynamic events in the development of mid-season sunburn in apple fruit

Apples grown in high heat, high light, and low humidity environments are at risk for sun injury disorders like sunburn and associated crop losses. Understanding the physiological and molecular mechanisms underlying sunburn will support improvement of mitigation strategies and breeding for more resilient varieties. Numerous studies have highlighted key biochemical processes involved in sun injury, such as the phenylpropanoid and reactive oxygen species (ROS) pathways, demonstrating both enzyme activities and expression of related genes in response to sunburn conditions. Most previous studies have focused on at-harvest activity of a small number of genes in response to heat stress. Thus, it remains unclear how stress events earlier in the season affect physiology and gene expression. Here, we applied heat stress to mid-season apples in the field and collected tissue along a time course-24, 48, and 72 h following a heat stimulus-to investigate dynamic gene expression changes using a transcriptomic lens. We found a relatively small number of differentially expressed genes (DEGs) and enriched functional terms in response to heat treatments. Only a few of these belonged to pathways previously described to be involved in sunburn, such as the AsA-GSH pathway, while most DEGs had not yet been implicated in sunburn or heat stress in pome fruit.

PearMODB: a multiomics database for pear (Pyrus) genomics, genetics and breeding study

Pear (Pyrus ssp.) belongs to Rosaceae and is an important fruit tree widely cultivated around the world. Currently, challenges to cope with the burgeoning sets of multiomics data are rapidly increasing. Here, we constructed the Pear Multiomics Database (PearMODB) by integrating genome, transcriptome, epigenome and population variation data, and aimed to provide a portal for accessing and analyzing pear multiomics data. A variety of online tools were built including gene search, BLAST, JBrowse, expression heatmap, synteny analysis and primer design. The information of DNA methylation sites and single-nucleotide polymorphisms can be retrieved through the custom JBrowse, providing an opportunity to explore the genetic polymorphisms linked to phenotype variation. Moreover, different gene families involving transcription factors, transcription regulators and disease resistance (nucleotide-binding site leucine-rich repeat) were identified and compiled for quick search. In particular, biosynthetic gene clusters (BGCs) were identified in pear genomes, and specialized webpages were set up to show detailed information of BGCs, laying a foundation for studying metabolic diversity among different pear varieties. Overall, PearMODB provides an important platform for pear genomics, genetics and breeding studies. Database URL http://pearomics.njau.edu.cn.

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.

Elucidation of the GAUT gene family in eight Rosaceae species and function analysis of PbrGAUT22 in pear pollen tube growth

The phylogenetic relationship and evolutionary history of the GAUT gene family were identified in 8 Rosaseae species. PbrGAUT22 was involved in controlling pollen tube growth by regulating the content of pectins. In plants, galacturonosyltransferases (GAUTs) were involved in homogalacturonan biosynthesis and functioned in maintaining pollen tube cell wall integrity. However, the feature and evolutionary history of the GAUT gene family in Rosaceae species and candidates in pear pollen tube growth remain unclear. Here, we identified 190 GAUT genes in 8 Rosaceae species, including Chinese white pear (Pyrus bretschneideri), European pear (Pyrus communis), apple (Malus × domestica), peach (Prunus persica), Japanese apricot (Prunus mume), sweet cherry (Prunus avium), woodland strawberry (Fragaria vesca) and black raspberry (Rubus occidentalis). Members in GAUT gene family were divided into 4 subfamilies according to the phylogenetic and structural analysis. Whole-genome duplication events and dispersed duplicates drove the expansion of the GAUT gene family. Among 23 pollen-expressed PbrGAUT genes in pear, PbrGAUT22 showed increased expression level during 1-6 h post-cultured pollen tubes. PbrGAUT22 was localized to the cytoplasm and plasma membrane. Knockdown of PbrGAUT22 expression in pollen tubes caused the decrease of pectin content and inhibited pear pollen tubes growth. Taken together, we investigated the identification and evolution of the GAUT gene family in Rosaceae species, and found that PbrGAUT22 played an essential role in the synthesis of pectin and the growth of pear pollen tubes.

PlantTribes2: Tools for comparative gene family analysis in plant genomics

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

The Role of Italy in the Use of Advanced Plant Genomic Techniques on Fruit Trees: State of the Art and Future Perspectives

Climate change is deeply impacting the food chain production, lowering quality and yield. In this context, the international scientific community has dedicated many efforts to enhancing resilience and sustainability in agriculture. Italy is among the main European producers of several fruit trees; therefore, national research centers and universities undertook several initiatives to maintain the specificity of the 'Made in Italy' label. Despite their importance, fruit crops are suffering from difficulties associated with the conventional breeding approaches, especially in terms of financial commitment, land resources availability, and long generation times. The 'new genomic techniques' (NGTs), renamed in Italy as 'technologies for assisted evolution' (TEAs), reduce the time required to obtain genetically improved cultivars while precisely targeting specific DNA sequences. This review aims to illustrate the role of the Italian scientific community in the use of NGTs, with a specific focus on Citrus, grapevine, apple, pear, chestnut, strawberry, peach, and kiwifruit. For each crop, the key genes and traits on which the scientific community is working, as well as the technological improvements and advancements on the regeneration of local varieties, are presented. Lastly, a focus is placed on the legal aspects in the European and in Italian contexts.

Comparison of structural variants detected by PacBio-CLR and ONT sequencing in pear

BACKGROUND

Structural variations (SVs) have recently become a topic of great interest in the area of genetic diversity and trait regulation. As genomic sequencing technologies have rapidly advanced, longer reads have been used to identify SVs at high resolution and with increased accuracy. It is important to choose a suitable sequencing platform and appropriate sequencing depth for SV detection in the pear genome.

RESULTS

In this study, two types of long reads from sequencing platforms, continuous long reads from Pacific Biosciences (PB-CLR) and long reads from Oxford Nanopore Technologies (ONT), were used to comprehensively analyze and compare SVs in the pear genome. The mapping rate of long reads was higher when the program Minimap2 rather than the other three mapping tools (NGMLR, LRA and Winnowmap2) was used. Three SV detection programs (Sniffles_v2, CuteSV, and Nanovar) were compared, and Nanovar had the highest sensitivity in detecting SVs at low sequencing depth (10-15×). A sequencing depth of 15× was suitable for SV detection in the pear genome using Nanovar. SVs detected by Sniffles_v2 and CuteSV with ONT reads had the high overlap with presence/absence variations (PAVs) in the pear cultivars 'Bartlett' and 'Dangshansuli', both of them with 38% of insertions and 55% of deletions overlapping with PAVs at sequencing depth of 30×. For the ONT sequencing data, over 37,526 SVs spanning ~ 28 Mb were identified by all three software packages for the 'Bartlett' and 'Dangshansuli' genomes. Those SVs were annotated and combined with transcriptome profiles derived from 'Bartlett' and 'Dangshansuli' fruit flesh at 60 days after cross-pollination. Several genes related to levels of sugars, acid, stone cells, and aromatic compounds were identified among the SVs. Transcription factors were then predicted among those genes, and results included bHLH, ERF, and MYB genes.

CONCLUSION

SV detection is of great significance in exploring phenotypic differences between pear varieties. Our study provides a framework for assessment of different SV software packages and sequencing platforms that can be applied in other plant genome studies. Based on these analyses, ONT sequencing data was determined to be more suitable than PB-CLR for SV detection in the pear genome. This analysis model will facilitate screening of genes related to agronomic traits in other crops.

Building a foundation for gene family analysis in Rosaceae genomes with a novel workflow: A case study in Pyrus architecture genes

The rapid development of sequencing technologies has led to a deeper understanding of plant genomes. However, direct experimental evidence connecting genes to important agronomic traits is still lacking in most non-model plants. For instance, the genetic mechanisms underlying plant architecture are poorly understood in pome fruit trees, creating a major hurdle in developing new cultivars with desirable architecture, such as dwarfing rootstocks in European pear (Pyrus communis). An efficient way to identify genetic factors for important traits in non-model organisms can be to transfer knowledge across genomes. However, major obstacles exist, including complex evolutionary histories and variable quality and content of publicly available plant genomes. As researchers aim to link genes to traits of interest, these challenges can impede the transfer of experimental evidence across plant species, namely in the curation of high-quality, high-confidence gene models in an evolutionary context. Here we present a workflow using a collection of bioinformatic tools for the curation of deeply conserved gene families of interest across plant genomes. To study gene families involved in tree architecture in European pear and other rosaceous species, we used our workflow, plus a draft genome assembly and high-quality annotation of a second P. communis cultivar, 'd'Anjou.' Our comparative gene family approach revealed significant issues with the most recent 'Bartlett' genome - primarily thousands of missing genes due to methodological bias. After correcting assembly errors on a global scale in the 'Bartlett' genome, we used our workflow for targeted improvement of our genes of interest in both P. communis genomes, thus laying the groundwork for future functional studies in pear tree architecture. Further, our global gene family classification of 15 genomes across 6 genera provides a valuable and previously unavailable resource for the Rosaceae research community. With it, orthologs and other gene family members can be easily identified across any of the classified genomes. Importantly, our workflow can be easily adopted for any other plant genomes and gene families of interest.

Molecular characterization of intergeneric hybrids between Malus and Pyrus

Apple (Malus) and pear (Pyrus) are economically important fruit crops well known for their unique textures, flavours, and nutritional qualities. Both genera are characterised by a distinct pattern of secondary metabolites, which directly affect not only resistance to certain diseases, but also have significant impacts on the flavour and nutritional value of the fruit. The identical chromosome numbers, similar genome size, and their recent divergence date, together with DNA markers have shown that apple and pear genomes are highly co-linear. This study utilized comparative genomic approaches, including simple sequence repeats, high resolution single nucleotide polymorphism melting analysis, and single nucleotide polymorphism chip analysis to identify genetic differences among hybrids of Malus and Pyrus, and F2 offspring. This research has demonstrated and validated that these three marker types, along with metabolomics analysis are very powerful tools to detect and confirm hybridity of progeny derived from crosses between apple and pear in both cross directions. Furthermore, this work analysed the genus-specific metabolite patterns and the resistance to fire blight (Erwinia amylovora) in progeny. The findings of this work will enhance and accelerate the breeding of novel tree fruit crops that benefit producers and consumers, by enabling marker assisted selection of desired traits introgressed between pear and apple.

The mechanism analysis of exogenous melatonin in limiting pear fruit aroma decrease under low temperature storage

Exogenous melatonin (MT) is widely used in fruit preservation, and can increase the storage time and delay the quality deterioration. Firstly, it was found that 150 μM MT was the optimal concentration to treat 'Xinli No.7' under storage at 4 °C for 60 days. MT could significantly improve oxidase activity and inhibit the reduction of physiological indexes, including pulp hardness, weight loss, titratable acid and soluble solid content. MT could also reduce ethylene release and limit the reduction of fruit aroma. The average content of fruit aroma substance increased by 43.53%. A relevant RNA-Seq database was built to further explore the regulation mechanism of MT. A total of 2,761 differentially expressed genes (DEGs) were identified. DEGs were enriched in 64 functional groups and 191 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. DEGs were mainly enriched in alpha-linolenic acid metabolism, fatty acid metabolism and plant hormone signal transduction pathway. The gene pycom09g05270 belonging to long chain acyl-CoA synthetase family and participating in fatty acid metabolism pathway was identified, and its expression level was consistent with fragments per kilobase per million mapped reads (FPKM) values, implying that pycom09g05270 might play a vital role in maintaining quality during the storage process.

A phased, chromosome-scale genome of 'Honeycrisp' apple (Malus domestica)

The apple cultivar 'Honeycrisp' has superior fruit quality traits, cold hardiness, and disease resistance, making it a popular breeding parent. However, it suffers from several physiological disorders, production, and postharvest issues. Despite several available apple genome sequences, understanding of the genetic mechanisms underlying cultivar-specific traits remains lacking. Here, we present a highly contiguous, fully phased, chromosome-level genome of 'Honeycrisp' apples, using PacBio HiFi, Omni-C, and Illumina sequencing platforms, with two assembled haplomes of 674 Mbp and 660 Mbp, and contig N50 values of 32.8 Mbp and 31.6 Mbp, respectively. Overall, 47,563 and 48,655 protein-coding genes were annotated from each haplome, capturing 96.8-97.4% complete BUSCOs in the eudicot database. Gene family analysis reveals most 'Honeycrisp' genes are assigned into orthogroups shared with other genomes, with 121 'Honeycrisp'-specific orthogroups. This resource is valuable for understanding the genetic basis of important traits in apples and related Rosaceae species to enhance breeding efforts.

Cultivated hawthorn (Crataegus pinnatifida var. major) genome sheds light on the evolution of Maleae (apple tribe)

Cultivated hawthorn (Crataegus pinnatifida var. major) is an important medicinal and edible plant with a long history of use for health protection in China. Herein, we provide a de novo chromosome-level genome sequence of the hawthorn cultivar "Qiu Jinxing." We assembled an 823.41 Mb genome encoding 40 571 genes and further anchored the 779.24 Mb sequence into 17 pseudo-chromosomes, which account for 94.64% of the assembled genome. Phylogenomic analyses revealed that cultivated hawthorn diverged from other species within the Maleae (apple tribe) at approximately 35.4 Mya. Notably, genes involved in the flavonoid and triterpenoid biosynthetic pathways have been significantly amplified in the hawthorn genome. In addition, our results indicated that the Maleae share a unique ancient tetraploidization event; however, no recent independent whole-genome duplication event was specifically detected in hawthorn. The amplification of non-specific long terminal repeat retrotransposons contributed the most to the expansion of the hawthorn genome. Furthermore, we identified two paleo-sub-genomes in extant species of Maleae and found that these two sub-genomes showed different rearrangement mechanisms. We also reconstructed the ancestral chromosomes of Rosaceae and discussed two possible paleo-polyploid origin patterns (autopolyploidization or allopolyploidization) of Maleae. Overall, our study provides an improved context for understanding the evolution of Maleae species, and this new high-quality reference genome provides a useful resource for the horticultural improvement of hawthorn.

Integrative Analysis of the DICER-like (DCL) Genes From Peach (Prunus persica): A Critical Role in Response to Drought Stress

DICER-likes (DCLs) proteins are the core component for non-coding RNA (ncRNA) biogenesis, playing essential roles in some biological processes. The DCL family has been characterized in model plants, such as Arabidopsis, rice, and poplar. However, the evolutionary aspect and the expression mechanism under drought stress were scarce and have never been reported and characterized in one of the most important worldwide cultivated fruit trees, peach (Prunus persica). Eight DCLs genes in the Prunus persica genome were detected, in addition to 51 DCLs in the other seven Rosaceae genomes. The phylogenetic analysis with Arabidopsis thaliana and RTL1 gene as outgroups suggested that DCL members are divided into four clades: DCL1, DCL2, DCL3, and DCL4 with several gene gain/loss events of DCL gene copies through the evolutionary tract of the Rosacea family. The number of homologous DCL copies within each clade, along with the chromosomal location indicated gene duplication event of the DCL2 gene occurred once for the subfamily Amygdaloideae and twice for Pyrus communis and Prunus dulics and trice for the P. persica on Chromosome number 7 genes. Another duplication event was found for the DCL3 gene that occurred once for all the eight Rosaceae species with no match in A. thaliana. The DCL genetic similarity and activity was evaluated using BLASTp and previously published RNA-seq data among different tissues and over different time points of peach trees exposed to drought conditions. Finally, the expression pattern of PrupeDCLs in response to drought stress was identified, and two of these members, Prupe.7G047900 and Prupe.6G363600, were found as main candidate genes for response to drought stress. Our data presented here provide useful information for a better understanding of the molecular evolution of DCL genes in Rosaceae genomes, and the function of DCLs in P. persica.

The chromosome-level genome provides insight into the molecular mechanism underlying the tortuous-branch phenotype of Prunus mume

Plant with naturally twisted branches is referred to as a tortuous-branch plant, which have extremely high ornamental value due to their zigzag shape and the natural twisting of their branches. Prunus mume is an important woody ornamental plant. However, the molecular mechanism underlying this unique trait in Prunus genus is unknown. Here, we present a chromosome-level genome assembly of the cultivated P. mume var. tortuosa created using Oxford Nanopore combined with Hi-C scaffolding, which resulted in a 237.8 Mb genome assembly being anchored onto eight pseudochromosomes. Molecular dating indicated that P. mume is the most recently differentiated species in Prunus. Genes associated with cell division, development and plant hormones play essential roles in the formation of tortuous branch trait. A putative regulatory pathway for the tortuous branch trait was constructed based on gene expression levels. Furthermore, after transferring candidate PmCYCD genes into Arabidopsis thaliana, we found that seedlings overexpressing these genes exhibited curled rosette leaves. Our results provide insights into the evolutionary history of recently differentiated species in Prunus genus, the molecular basis of stem morphology, and the molecular mechanism underlying the tortuous branch trait and highlight the utility of multi-omics in deciphering the properties of P. mume plant architecture.

Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.)

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

Construction of a high-density bin-map and identification of fruit quality-related quantitative trait loci and functional genes in pear

Pear (Pyrus spp.) is one of the most common fruit crops grown in temperate regions worldwide. Genetic enhancement of fruit quality is a fundamental goal of pear breeding programs. The genetic control of pear fruit quality traits is highly quantitative, and development of high-density genetic maps can facilitate fine-mapping of quantitative trait loci (QTLs) and gene identification. Bin-mapping is a powerful method of constructing high-resolution genetic maps from large-scale genotyping datasets. We performed whole-genome sequencing of pear cultivars 'Niitaka' and 'Hongxiangsu' and their 176 F ₁ progeny to identify genome-wide single-nucleotide polymorphism (SNP) markers for constructing a high-density bin-map of pear. This analysis yielded a total of 1.93 million SNPs and a genetic bin-map of 3190 markers spanning 1358.5 cM, with an average adjacent interval of 0.43 cM. This bin-map, along with other high-density genetic maps in pear, improved the reference genome assembly from 75.5 to 83.7% by re-anchoring the scaffolds. A quantitative genetic analysis identified 148 QTLs for 18 fruit-related traits; among them, QTLs for stone cell content, several key monosaccharides, and fruit pulp acids were identified for the first time in pear. A gene expression analysis of six pear cultivars identified 399 candidates in the identified QTL regions, which showed expression specific to fruit developmental stages in pear. Finally, we confirmed the function of PbrtMT1, a tonoplast monosaccharide transporter-related gene responsible for the enhancement of fructose accumulation in pear fruit on linkage group 16, in a transient transformation experiment. This study provides genomic and genetic resources as well as potential candidate genes for fruit quality improvement in pear.

Genome-wide clonal variability in European pear "Rocha" using high-throughput sequencing

Pears (Pyrus) are one of the most economically important fruits worldwide. The Pyrus genus is characterized by a high degree of genetic variability between species and interspecific hybrids, and several studies have been performed to assess this variability for both cultivated and wild accessions. These studies have mostly been limited by the resolving power of traditional molecular markers, although in the recent past the availability of reference genome sequences or SNP arrays for pear have enhanced the capability of high-resolution genomics studies. These tools can also be applied to better understand the intra-varietal (or clonal) variability in pear. Here we report the first high resolution genomics analysis of a pear clonal population using whole genome sequencing (WGS). Results showed unique signatures for the accumulation of mutations and transposable element insertions in each clone, which are likely related to their history of propagation and cultivation. The nucleotide diversity remained low in the clonal collection with the exception of few genomic windows, suggesting that balancing selection may be occurring. These windows included mainly genes related to plant fertility. Regions with higher mutational load were partially associated with transcription factors, probably reflecting the distinctive phenotypes in the collection. The annotation of variants also revealed the theoretical disruption of relevant genes in pear. Taken together, the results from this study show that pear clones accumulate mutations differently, and that those mutations can play a role on pear phenotypes, meaning that the study of pear clonal populations can be relevant in genetic studies, mainly when comparing with traditional association studies.

Genome-wide identification of the SWEET gene family mediating the cold stress response in Prunus mume

The Sugars Will Eventually be Exported Transporter (SWEET) gene family encodes a family of sugar transporters that play essential roles in plant growth, reproduction, and biotic and abiotic stresses. Prunus mume is a considerable ornamental wood plant with high edible and medicinal values; however, its lack of tolerance to low temperature has severely limited its geographical distribution. To investigate whether this gene family mediates the response of P. mume to cold stress, we identified that the P. mume gene family consists of 17 members and divided the family members into four groups. Sixteen of these genes were anchored on six chromosomes, and one gene was anchored on the scaffold with four pairs of segmental gene duplications and two pairs of tandem gene duplications. Cis-acting regulatory element analysis indicated that the PmSWEET genes are potentially involved in P. mume development, including potentially regulating roles in procedure, such as circadian control, abscisic acid-response and light-response, and responses to numerous stresses, such as low-temperature and drought. We performed low-temperature treatment in the cold-tolerant cultivar 'Songchun' and cold-sensitive cultivar 'Zaolve' and found that the expression of four of 17 PmSWEETs was either upregulated or downregulated with prolonged treatment times. This finding indicates that these family members may potentially play a role in cold stress responses in P. mume. Our study provides a basis for further investigation of the role of SWEET proteins in the development of P. mume and its responses to cold stress.

Chromosome-scale and haplotype-resolved genome assembly of a tetraploid potato cultivar

Potato is the most widely produced tuber crop worldwide. However, reconstructing the four haplotypes of its autotetraploid genome remained an unsolved challenge. Here, we report the 3.1 Gb haplotype-resolved (at 99.6% precision), chromosome-scale assembly of the potato cultivar ‘Otava’ based on high-quality long reads, single-cell sequencing of 717 pollen genomes and Hi-C data. Unexpectedly, ~50% of the genome was identical-by-descent due to recent inbreeding, which was contrasted by highly abundant structural rearrangements involving ~20% of the genome. Among 38,214 genes, only 54% were present in all four haplotypes with an average of 3.2 copies per gene. Taking the leaf transcriptome as an example, 11% of the genes were differently expressed in at least one haplotype, where 25% of them were likely regulated through allele-specific DNA methylation. Our work sheds light on the recent breeding history of potato, the functional organization of its tetraploid genome and has the potential to strengthen the future of genomics-assisted breeding.

Haplotype-resolved genome assembly of the tetraploid potato cultivar ‘Otava’ sheds light on functional organization of the tetraploid genome and provides the potential for genomics-assisted breeding.

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".

The Gillenia trifoliata genome reveals dynamics correlated with growth and reproduction in Rosaceae

The Rosaceae family has striking phenotypic diversity and high syntenic conservation. Gillenia trifoliata is sister species to the Maleae tribe of apple and ~1000 other species. Gillenia has many putative ancestral features, such as herb/sub-shrub habit, dry fruit-bearing and nine base chromosomes. This coalescence of ancestral characters in a phylogenetically important species, positions Gillenia as a 'rosetta stone' for translational science within Rosaceae. We present genomic and phenological resources to facilitate the use of Gillenia for this purpose. The Gillenia genome is the first fully annotated chromosome-level assembly with an ancestral genome complement (x = 9), and with it we developed an improved model of the Rosaceae ancestral genome. MADS and NAC gene family analyses revealed genome dynamics correlated with growth and reproduction and we demonstrate how Gillenia can be a negative control for studying fleshy fruit development in Rosaceae.

Advances in genomics and genome editing for breeding next generation of fruit and nut crops

Fruit tree crops are an essential part of the food production systems and are key to achieve food and nutrition security. Genetic improvement of fruit trees by conventional breeding has been slow due to the long juvenile phase. Advancements in genomics and molecular biology have paved the way for devising novel genetic improvement tools like genome editing, which can accelerate the breeding of these perennial crops to a great extent. In this article, advancements in genomics of fruit trees covering genome sequencing, transcriptome sequencing, genome editing technologies (GET), CRISPR-Cas system based genome editing, potential applications of CRISPR-Cas9 in fruit tree crops improvement, the factors influencing the CRISPR-Cas editing efficiency and the challenges for CRISPR-Cas9 applications in fruit tree crops improvement are reviewed. Besides, base editing, a recently emerging more precise editing system, and the future perspectives of genome editing in the improvement of fruit and nut crops are covered.

High-quality genome assembly of 'Cuiguan' pear (Pyrus pyrifolia) as a reference genome for identifying regulatory genes and epigenetic modifications responsible for bud dormancy

Dormancy-associated MADS-box (DAM) genes serve as crucial regulators of the endodormancy cycle in rosaceous plants. Although pear DAM genes have been identified previously, the lack of a high-quality reference genome and techniques to study gene function have prevented accurate genome-wide analysis and functional verification of such genes. Additionally, the contribution of other genes to the regulation of endodormancy release remains poorly understood. In this study, a high-quality genome assembly for 'Cuiguan' pear (Pyrus pyrifolia), which is a leading cultivar with a low chilling requirement cultivated in China, was constructed using PacBio and Hi-C technologies. Using this genome sequence, we revealed that pear DAM genes were tandemly clustered on Chr8 and Chr15 and were differentially expressed in the buds between 'Cuiguan' and the high-chilling-requirement cultivar 'Suli' during the dormancy cycle. Using a virus-induced gene silencing system, we determined the repressive effects of DAM genes on bud break. Several novel genes potentially involved in the regulation of endodormancy release were identified by RNA sequencing and H3K4me3 chromatin immunoprecipitation sequencing analyses of 'Suli' buds during artificial chilling using the new reference genome. Our findings enrich the knowledge of the regulatory mechanism underlying endodormancy release and chilling requirements and provide a foundation for the practical regulation of dormancy release in fruit trees as an adaptation to climate change.

Novel and emerging biotechnological crop protection approaches

Traditional breeding or genetically modified organisms (GMOs) have for a long time been the sole approaches to effectively cope with biotic and abiotic stresses and implement the quality traits of crops. However, emerging diseases as well as unpredictable climate changes affecting agriculture over the entire globe force scientists to find alternative solutions required to quickly overcome seasonal crises. In this review, we first focus on cisgenesis and genome editing as challenging biotechnological approaches for breeding crops more tolerant to biotic and abiotic stresses. In addition, we take into consideration a toolbox of new techniques based on applications of RNA interference and epigenome modifications, which can be adopted for improving plant resilience. Recent advances in these biotechnological applications are mainly reported for non-model plants and woody crops in particular. Indeed, the characterization of RNAi machinery in plants is fundamental to transform available information into biologically or biotechnologically applicable knowledge. Finally, here we discuss how these innovative and environmentally friendly techniques combined with traditional breeding can sustain a modern agriculture and be of potential contribution to climate change mitigation.

An exploration of assembly strategies and quality metrics on the accuracy of the rewarewa (Knightia excelsa) genome

We used long read sequencing data generated from Knightia excelsa, a nectar-producing Proteaceae tree endemic to Aotearoa (New Zealand), to explore how sequencing data type, volume and workflows can impact final assembly accuracy and chromosome reconstruction. Establishing a high-quality genome for this species has specific cultural importance to Māori and commercial importance to honey producers in Aotearoa. Assemblies were produced by five long read assemblers using data subsampled based on read lengths, two polishing strategies and two Hi-C mapping methods. Our results from subsampling the data by read length showed that each assembler tested performed differently depending on the coverage and the read length of the data. Subsampling highlighted that input data with longer read lengths but perhaps lower coverage constructed more contiguous, kmers and gene-complete assemblies than short read length input data with higher coverage. The final genome assembly was constructed into 14 pseudochromosomes using an initial flye long read assembly, a racon/medaka/pilon combined polishing strategy, salsa2 and allhic scaffolding, juicebox curation, and Macadamia linkage map validation. We highlighted the importance of developing assembly workflows based on the volume and read length of sequencing data and established a robust set of quality metrics for generating high-quality assemblies. Scaffolding analyses highlighted that problems found in the initial assemblies could not be resolved accurately by Hi-C data and that assembly scaffolding was more successful when the underlying contig assembly was of higher accuracy. These findings provide insight into how quality assessment tools can be implemented throughout genome assembly pipelines to inform the de novo reconstruction of a high-quality genome assembly for nonmodel organisms.

Relative performance of customized and universal probe sets in target enrichment: A case study in subtribe Malinae

PREMISE

Custom probe design for target enrichment in phylogenetics is tedious and often hinders broader phylogenetic synthesis. The universal angiosperm probe set Angiosperms353 may be the solution. Here, we test the relative performance of Angiosperms353 on the Rosaceae subtribe Malinae in comparison with custom probes that we specifically designed for this clade. We then address the impact of bioinformatically altering the performance of Angiosperms353 by replacing the original probe sequences with orthologs extracted from the Malus domestica genome.

METHODS

To evaluate the relative performance of these probe sets, we compared the enrichment efficiency, locus recovery, alignment length, proportion of parsimony-informative sites, proportion of potential paralogs, the topology and support of the resulting species trees, and the gene tree discordance.

RESULTS

Locus recovery was highest for our custom Malinae probe set, and replacing the original Angiosperms353 sequences with a Malus representative improved the locus recovery relative to Angiosperms353. The proportion of parsimony-informative sites was similar between all probe sets, while the gene tree discordance was lower in the case of the custom probes.

DISCUSSION

A custom probe set benefits from data completeness and can be tailored toward the specificities of the project of choice; however, Angiosperms353 was equally as phylogenetically informative as the custom probes. We therefore recommend using both a custom probe set and Angiosperms353 to facilitate large-scale systematic studies, where financially possible.