15 years of GDR: New data and functionality in the Genome Database for Rosaceae

Genome-wide identification, expression analysis, and potential roles under low-temperature stress of bHLH gene family in Prunus sibirica

The basic helix-loop-helix (bHLH) family is one of the most well-known transcription factor families in plants, and it regulates growth, development, and abiotic stress responses. However, systematic analyses of the bHLH gene family in Prunus sibirica have not been reported to date. In this study, 104 PsbHLHs were identified and classified into 23 subfamilies that were unevenly distributed on eight chromosomes. Nineteen pairs of segmental replication genes and ten pairs of tandem replication genes were identified, and all duplicated gene pairs were under purifying selection. PsbHLHs of the same subfamily usually share similar motif compositions and exon-intron structures. PsbHLHs contain multiple stress-responsive elements. PsbHLHs exhibit functional diversity by interacting and coordinating with other members. Twenty PsbHLHs showed varying degrees of expression. Eleven genes up-regulated and nine genes down-regulated in -4°C. The majority of PsbHLHs were highly expressed in the roots and pistils. Transient transfection experiments demonstrated that transgenic plants with overexpressed PsbHLH42 have better cold tolerance. In conclusion, the results of this study have significant implications for future research on the involvement of bHLH genes in the development and stress responses of Prunus sibirica.

Comprehensive Analysis of the Pectate Lyase Gene Family and the Role of FaPL1 in Strawberry Softening

Fruit softening is a crucial factor that controls shelf life and commercial value. Pectate lyase (PL) has a major role in strawberry fruit softening. However, the PL gene family in strawberry has not been comprehensively analyzed. In this study, 65 FaPL genes were identified in the octoploid strawberry genome. Subcellular localization prediction indicated that FaPLs are mostly localized to the extracellular and cytoplasmic spaces. Duplication event analysis suggested that FaPL gene family expansion is mainly driven by whole genome or segmental duplication. The FaPL family members were classified into six groups according to the phylogenetic analysis. Among them, FaPL1, 3, 5, 20, 25, 42, and 57 had gradually increased expressions during strawberry fruit development and ripening and higher expression levels in the fruits with less firmness than that in firmer fruit. This result suggested that these members are involved in strawberry softening. Furthermore, overexpression of FaPL1 significantly reduced the fruit firmness, ascorbic acid (AsA), and malondialdehyde (MDA) content but obviously increased the anthocyanins, soluble proteins, and titratable acidity (TA), while it had no apparent effects on flavonoids, phenolics, and soluble sugar content. These findings provide basic information on the FaPL gene family for further functional research and indicate that FaPL1 plays a vital role in strawberry fruit softening.

Genome-Wide Identification and Expression Profiling Analysis of SWEET Family Genes Involved in Fruit Development in Plum (Prunus salicina Lindl)

SWEETs (sugars will eventually be exported transporters) play a vital role in longer-distance sugar transportation, and thus control carbon flow and energy metabolism in plants. SWEET genes have been identified in various plant species, but their functions in fruit development remain uncharacterized. Here, we isolated 15 putative PsSWEETs from the Prunus salicina genome. For further analysis, comprehensive bioinformatics methods were applied to determine the gene structure, chromosome distribution, phylogeny, cis-acting regulatory elements, and expression profiles of PsSWEETs. qRT-PCR analysis suggested that these SWEETs might have diverse functions in the development of plum fruit. The relative expression levels of PsSWEET1 and PsSWEET9 were obviously higher in ripened fruit than the ones in other developmental stages, suggesting their possible roles in the transport and accumulation of sugars in plum fruit. Positive correlations were found between the expression level of PsSWEET3/10/13 and the content of sucrose, and the expression level of PsSWEET2 and the content of fructose, respectively, during the development of 'Furongli' fruit, suggesting their possible roles in the accumulation of sucrose and fructose. The current study investigated the initial genomic characterization and expression patterns of the SWEET gene family in plum, which could provide a foundation for the further understanding of the functional analysis of the SWEET gene family.

GWAS of adventitious root formation in roses identifies a putative phosphoinositide phosphatase (SAC9) for marker-assisted selection

Rose propagation by cuttings is limited by substantial genotypic differences in adventitious root formation. To identify possible genetic factors causing these differences and to develop a marker for marker-assisted selection for high rooting ability, we phenotyped 95 cut and 95 garden rose genotypes in a hydroponic rooting system over 6 weeks. Data on rooting percentage after 3 to 6 weeks, root number, and root fresh mass were highly variable among genotypes and used in association mappings performed on genotypic information from the WagRhSNP 68 K Axiom SNP array for roses. GWAS analyses revealed only one significantly associated SNP for rooting percentage after 3 weeks. Nevertheless, prominent genomic regions/peaks were observed and further analysed for rooting percentage after 6 weeks, root number and root fresh mass. Some of the SNPs in these peak regions were associated with large effects on adventitious root formation traits. Very prominent were ten SNPs, which were all located in a putative phosphoinositide phosphatase SAC9 on chromosome 2 and showed very high effects on rooting percentage after 6 weeks of more than 40% difference between nulliplex and quadruplex genotypes. SAC9 was reported to be involved in the regulation of endocytosis and in combination with other members of the SAC gene family to regulate the translocation of auxin-efflux PIN proteins via the dephosphorylation of phosphoinositides. For one SNP within SAC9, a KASP marker was successfully derived and used to select genotypes with a homozygous allele configuration. Phenotyping these homozygous genotypes for adventitious root formation verified the SNP allele dosage effect on rooting. Hence, the presented KASP derived from a SNP located in SAC9 can be used for marker-assisted selection in breeding programs for high rooting ability in the future.

Pedigree-based QTL analysis of flower size traits in two multi-parental diploid rose populations

Rose (Rosa spp.) is one of the most economically important ornamental species worldwide. Flower diameter, flower weight, and the number of petals and petaloids are key flower-size parameters and attractive targets for DNA-informed breeding. Pedigree-based analysis (PBA) using FlexQTL software was conducted using two sets of multi-parental diploid rose populations. Phenotypic data for flower diameter (Diam), flower weight (fresh (FWT)/dry (DWT)), number of petals (NP), and number of petaloids (PD) were collected over six environments (seasons) at two locations in Texas. The objectives of this study were to 1) identify new and/or validate previously reported QTL(s); 2) identify SNP haplotypes associated with QTL alleles (Q-/q-) of a trait and their sources; and 3) determine QTL genotypes for important rose breeding parents. Several new and previously reported QTLs for NP and Diam traits were identified. In addition, QTLs associated with flower weight and PD were identified for the first time. Two major QTLs with large effects were mapped for all traits. The first QTL was at the distal end of LG1 (60.44-60.95 Mbp) and was associated with Diam and DWT in the TX2WOB populations. The second QTL was consistently mapped in the middle region on LG3 (30.15-39.34 Mbp) and associated with NP, PD, and flower weight across two multi-parent populations (TX2WOB and TX2WSE). Haplotype results revealed a series of QTL alleles with differing effects at important loci for most traits. This work is distinct from previous studies by conducting co-factor analysis to account for the DOUBLE FLOWER locus while mapping QTL for NP. Sources of high-value (Q) alleles were identified, namely, 'Old Blush' and Rosa wichuraiana from J14-3 for Diam, while 'Violette' and PP-J14-3 were sources for other traits. In addition, the source of the low-value (q) alleles for Diam was 'Little Chief', and Rosa wichuraiana through J14-3 was the source for the remaining traits. Hence, our results can potentially inform parental/seedling selections as means to improve ornamental quality in roses and a step towards implementing DNA-informed techniques for use in rose breeding programs.

QTL discovery for resistance to black spot and cercospora leaf spot, and defoliation in two interconnected F1 bi-parental tetraploid garden rose populations

Garden roses are an economically important horticultural crop worldwide, and two major fungal pathogens, black spot (Diplocarpon rosae F.A. Wolf) and cercospora leaf spot of rose (Rosisphaerella rosicola Pass.), affect both the health and ornamental value of the plant. Most studies on black spot disease resistance have focused on diploid germplasm, and little work has been performed on cercospora leaf spot resistance. With the use of newly developed software tools for autopolyploid genetics, two interconnected tetraploid garden rose F1 populations (phenotyped over the course of 3 years) were used for quantitative trait locus (QTL) analysis of black spot and cercospora leaf spot resistance as well as plant defoliation. QTLs for black spot resistance were mapped to linkage groups (LGs) 1-6. QTLs for cercospora resistance and susceptibility were found in LGs 1, 4, and 5 and for defoliation in LGs 1, 3, and 5. The major locus on LG 5 for black spot resistance coincides with the previously discovered Rdr4 locus inherited from Rosa L. 'Radbrite' (Brite Eyes™), the common parent used in these mapping populations. This work is the first report of any QTL for cercospora resistance/susceptibility in tetraploid rose germplasm and the first report of defoliation QTL in roses. A major QTL for cercospora susceptibility coincides with the black spot resistance QTL on LG 5 (Rdr4). A major cercospora resistance QTL was found on LG 1. These populations provide a genetic resource that will further the knowledge base of rose genetics as more traits are studied. Studying more traits from these populations will allow for the stacking of various QTLs for desirable traits.

Metabolomic and transcriptomice analyses of flavonoid biosynthesis in apricot fruits

Introduction

Flavonoids, as secondary metabolites in plants, play important roles in many biological processes and responses to environmental factors.

Methods

Apricot fruits are rich in flavonoid compounds, and in this study, we performed a combined metabolomic and transcriptomic analysis of orange flesh (JN) and white flesh (ZS) apricot fruits.

Results and discussion

A total of 222 differentially accumulated flavonoids (DAFs) and 15855 differentially expressed genes (DEGs) involved in flavonoid biosynthesis were identified. The biosynthesis of flavonoids in apricot fruit may be regulated by 17 enzyme-encoding genes, namely PAL (2), 4CL (9), C4H (1), HCT (15), C3'H (4), CHS (2), CHI (3), F3H (1), F3'H (CYP75B1) (2), F3'5'H (4), DFR (4), LAR (1), FLS (3), ANS (9), ANR (2), UGT79B1 (6) and CYP81E (2). A structural gene-transcription factor (TF) correlation analysis yielded 3 TFs (2 bHLH, 1 MYB) highly correlated with 2 structural genes. In addition, we obtained 26 candidate genes involved in the biosynthesis of 8 differentially accumulated flavonoids metabolites in ZS by weighted gene coexpression network analysis. The candidate genes and transcription factors identified in this study will provide a highly valuable molecular basis for the in-depth study of flavonoid biosynthesis in apricot fruits.

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.

RoseAP: an analytical platform for gene function of Rosa rugosa

Rosa rugosa, a perennial shrub belonging to family Rosaceae, is a well-known ornamental plant. Its petals contain an abundance of essential oils and anthocyanins with enormous economic and health benefits when used as edible or cosmetic ingredients. The whole genome of R. rugosa was sequenced in 2021, which provided opportunities and challenges for gene regulation. However, many gene functions remain unknown. Therefore, an analytical platform named RoseAP (http://www.gzybioinformatics.cn/RoseAP/index.php) for the functional analysis of R. rugosa genes was constructed. It improved the gene annotation rate by integrating and analyzing genomic and transcriptomic datasets. First, 38,815 genes, covering 97.76% of the coding genes, were annotated functionally and structurally using a variety of algorithms and rules. Second, a total of 33 transcriptome samples were integrated, including 23 samples from our lab and 10 samples from the SRA database. A co-expression network containing approximately 29,657 positive or negative gene pairs, covering 74.7% of the coding genes, was constructed based on PCC and MR algorithms. Network analysis revealed that the DFR function was closely related to anthocyanin metabolism. It demonstrated the reliability of the network. Several SAUR genes of R. rugosa shared similar expression patterns. RoseAP was used to determine the sequence, structure, functional annotation, expression profile, regulatory network, and functional modules at the transcriptional and protein levels by inputting gene IDs. In addition, auxiliary analytical tools, including BLAST, gene set enrichment, orthologue conversion, gene sequence extraction, gene expression value extraction, and JBrowse, were utilized. Regular updates to RoseAP are expected to facilitate mining of gene function and promote genetic improvement in R. rugosa.

Genome of tetraploid sour cherry (Prunus cerasus L.) 'Montmorency' identifies three distinct ancestral Prunus genomes

Sour cherry (Prunus cerasus L.) is a valuable fruit crop in the Rosaceae family and a hybrid between progenitors closely related to extant Prunus fruticosa (ground cherry) and Prunus avium (sweet cherry). Here we report a chromosome-scale genome assembly for sour cherry cultivar Montmorency, the predominant cultivar grown in the USA. We also generated a draft assembly of P. fruticosa to use alongside a published P. avium sequence for syntelog-based subgenome assignments for 'Montmorency' and provide compelling evidence P. fruticosa is also an allotetraploid. Using hierarchal k-mer clustering and phylogenomics, we show 'Montmorency' is trigenomic, containing two distinct subgenomes inherited from a P. fruticosa-like ancestor (A and A') and two copies of the same subgenome inherited from a P. avium-like ancestor (BB). The genome composition of 'Montmorency' is AA'BB and little-to-no recombination has occurred between progenitor subgenomes (A/A' and B). In Prunus, two known classes of genes are important to breeding strategies: the self-incompatibility loci (S-alleles), which determine compatible crosses, successful fertilization, and fruit set, and the Dormancy Associated MADS-box genes (DAMs), which strongly affect dormancy transitions and flowering time. The S-alleles and DAMs in 'Montmorency' and P. fruticosa were manually annotated and support subgenome assignments. Lastly, the hybridization event 'Montmorency' is descended from was estimated to have occurred less than 1.61 million years ago, making sour cherry a relatively recent allotetraploid. The 'Montmorency' genome highlights the evolutionary complexity of the genus Prunus and will inform future breeding strategies for sour cherry, comparative genomics in the Rosaceae, and questions regarding neopolyploidy.

Evolution of the 14-3-3 gene family in monocotyledons and dicotyledons and validation of MdGRF13 function in transgenic Arabidopsis thaliana

KEY MESSAGE

The 14-3-3 family is more highly conserved among monocotyledons, and overexpression of MdGRF13 improved drought and salt tolerance in transgenic Arabidopsis thaliana. The 14-3-3 are highly conserved regulatory proteins found in eukaryotes and play an essential role in plant growth, development and stress response. However, the 14-3-3 gene family evolution in monocotyledons and dicotyledons and the biological functions of the MdGRF13 under abiotic stress remain unknown. In our study, 195 members of the 14-3-3 family were identified from 12 species and divided into ε group and the Non-ε group. Synteny analysis within the 14-3-3 family indicated that segmental duplication events contributed to the expansion of the family. Selective pressure analysis indicated that purifying selection was a vital force in the 14-3-3 genes evolution, and monocotyledons had a lower million years ago (Mya) mean values than dicotyledons. Meanwhile, the codon adaptation index (CAI) and frequency of optical codons (FOP) are higher and the effective number of codons (Nc) is lower in monocotyledons 14-3-3 genes compared to dicotyledons. Moreover, the yeast two-hybrid (Y2H) demonstrated that MdGRF13 interacts with MdRD22, MdLHP1a and MdMORF1. Significantly, the malondialdehyde (MDA) content and relative conductivity were decreased, while the superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities were increased in transgenic Arabidopsis compared to the wild type (WT) under drought and salt stress. These results suggest that overexpression of MdGRF13 significantly improved the tolerance to drought and salt stress in transgenic Arabidopsis. Thus, our results provide a theoretical basis for exploring the evolution and function of the 14-3-3 gene family in monocotyledons and dicotyledons.

Defying Gravity: WEEP promotes negative gravitropism in Prunus persica (peach) shoots and roots by establishing asymmetric auxin gradients

Trees with weeping shoot architectures are valued for their beauty and serve as tremendous resources for understanding how plants regulate posture control. The Prunus persica (peach) weeping phenotype, which has elliptical downward arching branches, is caused by a homozygous mutation in the WEEP gene. Until now, little was known about the function of WEEP protein despite its high conservation throughout Plantae. Here, we present the results of anatomical, biochemical, biomechanical, physiological, and molecular experiments that provide insight into WEEP function. Our data suggest that weeping peach does not have defects in branch structure. Rather, transcriptomes from the adaxial (upper) and abaxial (lower) sides of standard and weeping branch shoot tips revealed flipped expression patterns for genes associated with early auxin response, tissue patterning, cell elongation, and tension wood development. This suggests that WEEP promotes polar auxin transport toward the lower side during shoot gravitropic response, leading to cell elongation and tension wood development. In addition, weeping peach trees exhibited steeper root systems and faster root gravitropic response, just as barley and wheat with mutations in their WEEP homolog EGT2. This suggests that the role of WEEP in regulating lateral organ angles and orientations during gravitropism may be conserved. Additionally, size-exclusion chromatography indicated that WEEP proteins self-oligomerize, like other SAM-domain proteins. This oligomerization may be required for WEEP to function in formation of protein complexes during auxin transport. Collectively, our results from weeping peach provide new insight into polar auxin transport mechanisms associated with gravitropism and lateral shoot and root orientation.

Identification and characterization of cytochrome P450 CYP77A59 of loquat (Rhaphiolepis bibas) responsible for biosynthesis of phenylacetonitrile, a floral nitrile compound

MAIN CONCLUSION

Cytochrome P450 CYP77A59 is responsible for the biosynthesis of phenylacetonitrile in loquat flowers. Flowers of some plants emit volatile nitrile compounds, but the biosynthesis of these compounds is unclear. Loquat (Rhaphiolepis bibas) flowers emit characteristic N-containing volatiles, such as phenylacetonitrile (PAN), (E/Z)-phenylacetaldoxime (PAOx), and (2-nitroethyl)benzene (NEB). These volatiles likely play a defense role against pathogens and insects. PAN and NEB are commonly biosynthesized from L-phenylalanine via (E/Z)-PAOx. Two cytochrome P450s-CYP79D80 and "promiscuous fatty acid ω-hydroxylase" CYP94A90, which catalyze the formation of (E/Z)-PAOx from L-phenylalanine and NEB from (E/Z)-PAOx, respectively-are involved in NEB biosynthesis. However, the enzymes catalyzing the formation of PAN from (E/Z)-PAOx in loquat have not been identified. In this study, we aimed to identify candidate cytochrome P450s catalyzing PAN formation in loquat flowers. Yeast whole-cell biocatalyst assays showed that among nine candidate cytochrome P450s, CYP77A58 and CYP77A59 produced PAN from (E/Z)-PAOx. CYP77As catalyzed the dehydration of aldoximes, which is atypical of cytochrome P450; the reaction was NADPH-dependent, with an optimum temperature and pH of 40 °C and 8.0, respectively. CYP77As acted on (E/Z)-PAOx, (E/Z)-4-hydroxyphenylacetaldoxime, and (E/Z)-indole-3-acetaldoxime. Previously characterized CYP77As are known to hydroxylate fatty acids; loquat CYP77As did not act on tested fatty acids. We observed higher expression of CYP77A59 in flowers than in buds; expression of CYP77A58 was remarkably reduced in the flowers. Because the flowers, but not buds, emit PAN, CYP77A59 is likely responsible for the biosynthesis of PAN in loquat flowers. This study will help us understand the biosynthesis of floral nitrile compounds.

Insights into the Genes Involved in ABA Biosynthesis and Perception during Development and Ripening of the Chilean Strawberry Fruit

Hormones act as master ripening regulators. In non-climacteric fruit, ABA plays a key role in ripening. Recently, we confirmed in Fragaria chiloensis fruit that in response to ABA treatment the fruit induces ripening-associated changes such as softening and color development. In consequence of these phenotypic changes, transcriptional variations associated with cell wall disassembly and anthocyanins biosynthesis were reported. As ABA stimulates the ripening of F. chiloensis fruit, the molecular network involved in ABA metabolism was analyzed. Therefore, the expression level of genes involved in ABA biosynthesis and ABA perception was quantified during the development of the fruit. Four NCED/CCDs and six PYR/PYLs family members were identified in F. chiloensis. Bioinformatics analyses confirmed the existence of key domains related to functional properties. Through RT-qPCR analyses, the level of transcripts was quantified. FcNCED1 codifies a protein that displays crucial functional domains, and the level of transcripts increases as the fruit develops and ripens, in parallel with the increment in ABA. In addition, FcPYL4 codifies for a functional ABA receptor, and its expression follows an incremental pattern during ripening. The study concludes that FcNCED1 is involved in ABA biosynthesis; meanwhile, FcPYL4 participates in ABA perception during the ripening of F. chiloensis fruit.

Roles and evolution of four LEAFY homologs in floral patterning and leaf development in woodland strawberry

The plant-specific transcription factor LEAFY (LFY), generally maintained as a single-copy gene in most angiosperm species, plays critical roles in flower development. The woodland strawberry (Fragaria vesca) possesses four LFY homologs in the genome; however, their respective functions and evolution remain unknown. Here, we identified and validated that mutations in one of the four LFY homologs, FveLFYa, cause homeotic conversion of floral organs and reiterative outgrowth of ectopic flowers. In contrast to FveLFYa, FveLFYb/c/d appear dispensable under normal growth conditions, as fvelfyc mutants are indistinguishable from wild type and FveLFYb and FveLFYd are barely expressed. Transgenic analysis and yeast one-hybrid assay showed that FveLFYa and FveLFYb, but not FveLFYc and FveLFYd, are functionally conserved with AtLFY in Arabidopsis (Arabidopsis thaliana). Unexpectedly, LFY-binding site prediction and yeast one-hybrid assay revealed that the transcriptional links between LFY and the APETALA1 (AP1) promoter/the large AGAMOUS (AG) intron are missing in F. vesca, which is due to the loss of LFY-binding sites. The data indicate that mutations in cis-regulatory elements could contribute to LFY evolution. Moreover, we showed that FveLFYa is involved in leaf development, as approximately 30% of mature leaves have smaller or fewer leaflets in fvelfya. Phylogenetic analysis indicated that LFY homologs in Fragaria species may arise from recent duplication events in their common ancestor and are undergoing convergent gene loss. Together, these results provide insight into the role of LFY in flower and leaf development in strawberry and have important implications for the evolution of LFY.

Meta-Analysis of Rose Rosette Disease-Resistant Quantitative Trait Loci and a Search for Candidate Genes

Rose rosette disease (RRD), caused by the rose rosette emaravirus (RRV), is a major viral disease in roses (Rosa sp.) that threatens the rose industry. Recent studies have revealed quantitative trait loci (QTL) for reduced susceptibility to RRD in the linkage groups (LGs) 1, 5, 6, and 7 in tetraploid populations and the LGs 1, 3, 5, and 6 in diploid populations. In this study, we seek to better localize and understand the relationship between QTL identified in both diploid and tetraploid populations. We do so by remapping the populations found in these studies and performing a meta-analysis. This analysis reveals that the peaks and intervals for QTL using diploid and tetraploid populations co-localized on LG 1, suggesting that these are the same QTL. The same was seen on LG 3. Three meta-QTL were identified on LG 5, and two were discovered on LG 6. The meta-QTL on LG 1, MetaRRD1.1, had a confidence interval (CI) of 10.53 cM. On LG 3, MetaRRD3.1 had a CI of 5.94 cM. MetaRRD5.1 had a CI of 17.37 cM, MetaRRD5.2 had a CI of 4.33 cM, and MetaRRD5.3 had a CI of 21.95 cM. For LG 6, MetaRRD6.1 and MetaRRD6.2 had CIs of 9.81 and 8.81 cM, respectively. The analysis also led to the identification of potential disease resistance genes, with a primary interest in genes localized in meta-QTL intervals on LG 5 as this LG was found to explain the greatest proportion of phenotypic variance for RRD resistance. The results from this study may be used in the design of more robust marker-based selection tools to track and use a given QTL in a plant breeding context.

Promoter Variation of the Key Apple Fruit Texture Related Gene MdPG1 and the Upstream Regulation Analysis

MdPG1 encoding polygalacturonase in apple (Malus × domestica) is a key gene associated with fruit firmness and texture variations among apple cultivars. However, the causative variants of MdPG1 are still not known. In this study, we identified a SNP^A/C variant within an ERF-binding element located in the promoter region of MdPG1. The promoter containing the ERF-binding element with SNP^A, rather than the SNP^C, could be strongly bound and activated by MdCBF2, a member of the AP2/ERF transcription factor family, as determined by yeast-one-hybrid and dual-luciferase reporter assays. We also demonstrated that the presence of a novel long non-coding RNA, lncRNA_PG1, in the promoter of MdPG1 was a causative variant. lncRNA_PG1 was specifically expressed in fruit tissues postharvest. lncRNA_PG1 could reduce promoter activity when it was fused to the promoter of MdPG1 and a tobacco gene encoding Mg-chelatase H subunit (NtCHLH) in transgenic tobacco cells but could not reduce promoter activity when it was supplied in a separate gene construct, indicating a cis-regulatory effect. Our results provide new insights into genetic regulation of MdPG1 allele expression and are also useful for the development of elite apple cultivars.

Genome-wide identification and expression profiling reveal the regulatory role of U-box E3 ubiquitin ligase genes in strawberry fruit ripening and abiotic stresses resistance

The plant U-box (PUB) proteins are a type of E3 ubiquitin ligases well known for their functions in response to various stresses. They are also related to fruit development and ripening. However, PUB members possess such roles that remain unclear in strawberry. In this study, 155 PUB genes were identified in octoploid strawberry and classified into four groups. Their promoters possessed a variety of cis-acting elements, most of which are associated with abiotic stresses, followed by phytohormones response and development. Protein-protein interaction analysis suggested that FaU-box members could interact with each other as well as other proteins involved in hormone signaling and stress resistance. Transcriptome-based and RT-qPCR expression analysis revealed the potential involvement of FaU-box genes in resistance to stresses and fruit ripening. Of these, FaU-box98 and FaU-box136 were positively while FaU-box52 was negatively related to strawberry ripening. FaU-box98 comprehensively participated in resistance of ABA, cold, and salt, while FaU-box83 and FaU-box136 were broadly associated with drought and salt stresses. FaU-box18 and FaU-box52 were ABA-specific; FaU-box3 was specific to salt stress. In addition, the functional analysis of a randomly selected FaU-box (FaU-box127) showed that the transient overexpression of FaU-box127 promoted the ripening of strawberry fruit, along with significant changes in the expression levels of some ripening-related genes and the content of organic acid and soluble sugar. Overall, these findings provided comprehensive information about the FaU-box gene family and identified the potential FaU-box members participating in stress resistance and strawberry fruit ripening regulation.

Biodiversity Evaluation and Preservation of Italian Stone Fruit Germplasm (Peach and Apricot) in Southern Italy

The Prunus genus encompasses a group of economically important and closely related crops, sharing an essentially common genome and, thereby, a high level of conserved and transferable microsatellite (SSR) loci. In Southern Italy, many of the local and/or neglected varieties are abandoned and at risk of extinction due to the high degree of urbanization and agricultural intensification, despite their value as genetic resources for crop improvement. This research aimed to genetically and morphologically characterize the traditional apricot (P. armenica) and peach (P. persica) germplasms collected in old family orchards. Most of the official descriptor categories were scored, thus revealing a rather high level of phenotypic variation in both collections. Genetic data allowed the discovery of diversity masked by morphological traits. Genotyping in 15 and 18 SSRs, eight of which were transferable across both species, showed an average polymorphic informativeness (PIC) of 0.44 and 0.59 for apricot and peach, respectively, and a total of 70 and 144 alleles. A reliable identification of each genotype was achieved, and the presence of possible mislabeling and/or erroneous denominations was solved. These results are encouraging for the valorization of the still poorly explored Italian Prunus germplasm, with significant economic consequences for bioresource conservation and management.

Genome-Wide Identification and Expression Analysis of the TIR-NBS-LRR Gene Family and Its Response to Fungal Disease in Rose (Rosa chinensis)

Roses, which are one of the world's most important ornamental plants, are often damaged by pathogens, resulting in serious economic losses. As a subclass of the disease resistance gene family of plant nucleotide-binding oligomerization domain (NOD)-like receptors, TIR-NBS-LRR (TNL) genes play a vital role in identifying pathogen effectors and activating defense responses. However, a systematic analysis of the TNL gene family is rarely reported in roses. Herein, 96 intact TNL genes were identified in Rosa chinensis. Their phylogenies, physicochemical characteristics, gene structures, conserved domains and motifs, promoter cis-elements, microRNA binding sites, and intra- and interspecific collinearity relationships were analyzed. An expression analysis using transcriptome data revealed that RcTNL genes were dominantly expressed in leaves. Some RcTNL genes responded to gibberellin, jasmonic acid, salicylic acid, Botrytis cinerea, Podosphaera pannosa, and Marssonina rosae (M. rosae); the RcTNL23 gene responded significantly to three hormones and three pathogens, and exhibited an upregulated expression. Furthermore, the black spot pathogen was identified as M. rosae. After inoculating rose leaves, an expression pattern analysis of the RcTNL genes suggested that they act during different periods of pathogen infection. The present study lays the foundations for an in-depth investigation of the TNL gene function and the mining of disease resistance genes in roses.

Genome-wide identification of WOX family members in nine Rosaceae species and a functional analysis of MdWOX13-1 in drought resistance

WUSCHEL-related homeobox (WOX) transcription factors (TFs) are important in plant development processes and evolutionary novelties. In this study, a genome-wide comprehensive analysis of WOX genes from nine Rosaceae species was carried out, and their potential roles in Malus were subsequently investigated. 125 WOXs in 9 Rosaceae species were identified and classified into three clades, i.e., the ancient, intermediate, and WUS clades. Prunus. domestica contained the most intra-genomic collinearity among the nine Rosaceae species. Additionally, the cis-elements in WOX gene family members were compared and classified into three categories, including phytohormone-responsive, plant growth and development, and abiotic and biotic stresses. Overexpression (OE) of MdWOX13-1 also increased the callus weight and enhanced ROS scavenging against drought stress. Furthermore, via yeast-one hybrid assay and LUC analyses, MdWOX13-1 could directly bind to the MdMnSOD promoter. Therefore, our results will facilitate further study of the WOX genes' function in the Rosaceae family.

A chromosome-length genome assembly and annotation of blackberry (Rubus argutus, cv. "Hillquist")

Blackberries (Rubus spp.) are the fourth most economically important berry crop worldwide. Genome assemblies and annotations have been developed for Rubus species in subgenus Idaeobatus, including black raspberry (R. occidentalis), red raspberry (R. idaeus), and R. chingii, but very few genomic resources exist for blackberries and their relatives in subgenus Rubus. Here we present a chromosome-length assembly and annotation of the diploid blackberry germplasm accession "Hillquist" (R. argutus). "Hillquist" is the only known source of primocane-fruiting (annual-fruiting) in tetraploid fresh-market blackberry breeding programs and is represented in the pedigree of many important cultivars worldwide. The "Hillquist" assembly, generated using Pacific Biosciences long reads scaffolded with high-throughput chromosome conformation capture sequencing, consisted of 298 Mb, of which 270 Mb (90%) was placed on 7 chromosome-length scaffolds with an average length of 38.6 Mb. Approximately 52.8% of the genome was composed of repetitive elements. The genome sequence was highly collinear with a novel maternal haplotype-resolved linkage map of the tetraploid blackberry selection A-2551TN and genome assemblies of R. chingii and red raspberry. A total of 38,503 protein-coding genes were predicted, of which 72% were functionally annotated. Eighteen flowering gene homologs within a previously mapped locus aligning to an 11.2 Mb region on chromosome Ra02 were identified as potential candidate genes for primocane-fruiting. The utility of the "Hillquist" genome has been demonstrated here by the development of the first genotyping-by-sequencing-based linkage map of tetraploid blackberry and the identification of possible candidate genes for primocane-fruiting. This chromosome-length assembly will facilitate future studies in Rubus biology, genetics, and genomics and strengthen applied breeding programs.

Deducing genotypes for loci of interest from SNP array data via haplotype sharing, demonstrated for apple and cherry

Breeders, collection curators, and other germplasm users require genetic information, both genome-wide and locus-specific, to effectively manage their genetically diverse plant material. SNP arrays have become the preferred platform to provide genome-wide genetic profiles for elite germplasm and could also provide locus-specific genotypic information. However, genotypic information for loci of interest such as those within PCR-based DNA fingerprinting panels and trait-predictive DNA tests is not readily extracted from SNP array data, thus creating a disconnect between historic and new data sets. This study aimed to establish a method for deducing genotypes at loci of interest from their associated SNP haplotypes, demonstrated for two fruit crops and three locus types: quantitative trait loci Ma and Ma3 for acidity in apple, apple fingerprinting microsatellite marker GD12, and Mendelian trait locus Rf for sweet cherry fruit color. Using phased data from an apple 8K SNP array and sweet cherry 6K SNP array, unique haplotypes spanning each target locus were associated with alleles of important breeding parents. These haplotypes were compared via identity-by-descent (IBD) or identity-by-state (IBS) to haplotypes present in germplasm important to U.S. apple and cherry breeding programs to deduce target locus alleles in this germplasm. While IBD segments were confidently tracked through pedigrees, confidence in allele identity among IBS segments used a shared length threshold. At least one allele per locus was deduced for 64-93% of the 181 individuals. Successful validation compared deduced Rf and GD12 genotypes with reported and newly obtained genotypes. Our approach can efficiently merge and expand genotypic data sets, deducing missing data and identifying errors, and is appropriate for any crop with SNP array data and historic genotypic data sets, especially where linkage disequilibrium is high. Locus-specific genotypic information extracted from genome-wide SNP data is expected to enhance confidence in management of genetic resources.

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.

Identification of PLATZ genes in Malus and expression characteristics of MdPLATZs in response to drought and ABA stresses

Plant AT-rich sequences and zinc-binding proteins (PLATZ) play crucial roles in response to environmental stresses. Nevertheless, PLATZ gene family has not been systemically studied in Rosaceae species, such as in apple, pear, peach, or strawberry. In this study, a total of 134 PLATZ proteins were identified from nine Rosaceae genomes and were classified into seven phylogenetic groups. Subsequently, the chromosomal localization, duplication, and collinearity relationship for apple PLATZ genes were investigated, and segmental duplication is a major driving-force in the expansion of PLATZ in Malus. Expression profiles analysis showed that PLATZs had distinct expression patterns in different tissues, and multiple genes were significantly changed after drought and ABA treatments. Furthermore, the co-expression network combined with RNA-seq data showed that PLATZ might be involved in drought stress by regulating ABA signaling pathway. In summary, this study is the first in-depth and systematic identification of PLATZ gene family in Rosaceae species, especially for apple, and provided specific PLATZ gene resource for further functional research in response to abiotic stress.

Ppe.CR.1 DNA test for predicting chilling requirement in peach

Chilling requirement (CR) is an important agronomic trait controlling the floral bud break for proper flowering in peach. Even though it has been widely researched and several peach CR quantitative trait loci (QTLs) have been identified, no diagnostic DNA tests validated in the U.S. peach breeding germplasm are available for this trait. Breeders and growers need a simple DNA test to predict the CR of peach cultivars for their particular environment. Therefore, we developed a quick and reliable Kompetitive Allele Specific PCR (KASP) DNA test using haplotype information from 9K IPSC genotype data of the U.S. peach germplasm integrating four CR-associated SNP markers from the previously reported CR QTL region on linkage group 1. Four KASP assays (Ppe.CR.1-1 to -4) were developed and validated on 77 peach cultivars, and nine accessions from two F₂ populations, with 96 and 74% accuracy in determining expected CR genotype (compared to SNP array) and predicting phenotype, respectively. Furthermore, the Ppe.CR.1 showed 80% accuracy in predicting the precise CR phenotype in the Clemson University peach breeding material. Only one Ppe.CR.1 KASP assay (Ppe.CR.1-1) is needed to distinguish between haplotypes with CR lower and higher than 800 chilling hours, and two Ppe.CR.1 assays (Pp.CR.1-1 and -4), are capable of distinguishing low, moderate, and high CR alleles. Coupled with the crude DNA extraction, the Ppe.CR.1 DNA test provides a low-cost option for breeders and growers to predict CR in peach material with more than 70% accuracy.

Transcriptome analyses reveal the effects of mixed saline-alkali stress on indoleacetic acid and cytokinins in Malus hupehensis Rehd. leaves

Saline-alkali stress is a universal abiotic stress factor limiting fruit tree cultivation worldwide. Apple (Malus×domestica Borkh.) is one of the fruits with the largest yields worldwide. Tea crabapple (Malus hupehensis Rehd. var. pingyiensis Jiang) is a type of common apple rootstock in China. Because facultative apomixis occurs in this species, it is often used in molecular research. The present study investigated the molecular mechanism of the response of indoleacetic acid (IAA) and cytokinins [zeatin, trans-zeatin riboside (tZR), isopentenyladenine (iP), and isopentenyladenosine (iPA)] to mixed saline-alkali stress (MSAS) in tea crabapple leaves. The endogenous hormone content of tea crabapple leaves under MSAS was measured, and the expression of stress response-related genes was analyzed by RNA sequencing. The results showed that the concentration of IAA was initially higher and then lower than that in the control, whereas the concentration of zeatin, tZR, iP, and iPA was higher than that in the control. A total of 1262 differentially expressed genes were identified in the three comparison groups. Further analyses suggested that IAA and cytokinin biosynthetic genes were mostly upregulated in tea crabapple leaves, indicating that auxin and cytokinin signaling pathway regulation occurred in response to MSAS. These findings suggest that IAA and cytokinins play an important role in the response of tea crabapple to MSAS.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01275-4.

Functional characterization of MANNOSE-BINDING LECTIN 1, a G-type lectin gene family member, in response to fungal pathogens of strawberry

The mannose-binding lectin gene MANNOSE-BINDING LECTIN 1 (MBL1) is a member of the G-type lectin family and is involved in defense in strawberry (Fragaria × ananassa). Genome-wide identification of the G-type lectin family was carried out in woodland strawberry, F. vesca, and 133 G-lectin genes were found. Their expression profiles were retrieved from available databases and indicated that many are actively expressed during plant development or interaction with pathogens. We selected MBL1 for further investigation and generated stable transgenic FaMBL1-overexpressing plants of F. ×ananassa to examine the role of this gene in defense. Plants were selected and evaluated for their contents of disease-related phytohormones and their reaction to biotic stresses, and this revealed that jasmonic acid decreased in the overexpressing lines compared with the wild-type (WT). Petioles of the overexpressing lines inoculated with Colletotrichum fioriniae had lower disease incidence than the WT, and leaves of these lines challenged by Botrytis cinerea showed significantly smaller lesion diameters than the WT and higher expression of CLASS II CHITINASE 2-1. Our results indicate that FaMBL1 plays important roles in strawberry response to fungal diseases caused by C. fioriniae and B. cinerea.

Chromosome-scale genome sequence assemblies of the 'Autumn Bliss' and 'Malling Jewel' cultivars of the highly heterozygous red raspberry (Rubus idaeus L.) derived from long-read Oxford Nanopore sequence data

Red raspberry (Rubus idaeus L.) is an economically valuable soft-fruit species with a relatively small (~300 Mb) but highly heterozygous diploid (2n = 2x = 14) genome. Chromosome-scale genome sequences are a vital tool in unravelling the genetic complexity controlling traits of interest in crop plants such as red raspberry, as well as for functional genomics, evolutionary studies, and pan-genomics diversity studies. In this study, we developed genome sequences of a primocane fruiting variety ('Autumn Bliss') and a floricane variety ('Malling Jewel'). The use of long-read Oxford Nanopore Technologies sequencing data yielded long read lengths that permitted well resolved genome sequences for the two cultivars to be assembled. The de novo assemblies of 'Malling Jewel' and 'Autumn Bliss' contained 79 and 136 contigs respectively, and 263.0 Mb of the 'Autumn Bliss' and 265.5 Mb of the 'Malling Jewel' assembly could be anchored unambiguously to a previously published red raspberry genome sequence of the cultivar 'Anitra'. Single copy ortholog analysis (BUSCO) revealed high levels of completeness in both genomes sequenced, with 97.4% of sequences identified in 'Autumn Bliss' and 97.7% in 'Malling Jewel'. The density of repetitive sequence contained in the 'Autumn Bliss' and 'Malling Jewel' assemblies was significantly higher than in the previously published assembly and centromeric and telomeric regions were identified in both assemblies. A total of 42,823 protein coding regions were identified in the 'Autumn Bliss' assembly, whilst 43,027 were identified in the 'Malling Jewel' assembly. These chromosome-scale genome sequences represent an excellent genomics resource for red raspberry, particularly around the highly repetitive centromeric and telomeric regions of the genome that are less complete in the previously published 'Anitra' genome sequence.

Genome-Wide Identification and Characterization of Argonaute, Dicer-like and RNA-Dependent RNA Polymerase Gene Families and Their Expression Analyses in Fragaria spp

In the growth and development of plants, some non-coding small RNAs (sRNAs) not only mediate RNA interference at the post-transcriptional level, but also play an important regulatory role in chromatin modification at the transcriptional level. In these processes, the protein factors Argonaute (AGO), Dicer-like (DCL), and RNA-dependent RNA polymerase (RDR) play very important roles in the synthesis of sRNAs respectively. Though they have been identified in many plants, the information about these gene families in strawberry was poorly understood. In this study, using a genome-wide analysis and a phylogenetic approach, 13 AGO, six DCL, and nine RDR genes were identified in diploid strawberry Fragaria vesca. We also identified 33 AGO, 18 DCL, and 28 RDR genes in octoploid strawberry Fragaria × ananassa, studied the expression patterns of these genes in various tissues and developmental stages of strawberry, and researched the response of these genes to some hormones, finding that almost all genes respond to the five hormone stresses. This study is the first report of a genome-wide analysis of AGO, DCL, and RDR gene families in Fragaria spp., in which we provide basic genomic information and expression patterns for these genes. Additionally, this study provides a basis for further research on the functions of these genes and some evidence for the evolution between diploid and octoploid strawberries.

PIPs from Fragaria vesca: A structural analysis of native and mutated protein

Aquaporins are an ancient family of membrane channel proteins present in all eukaryotes and most prokaryotes, and apart from water, allow the transport of neutral solutes and organic compounds through the pore. These proteins are essential role differentially expressed during ripening in Fragaria vesca fruits. Fv PIP2-1a is intensively expressed in fruits, inclusive several other proteins member are differentially expressed in fruit but also in other plant tissues. Phylogenetic analysis shows that Fv PIP2-1a grouped with other Fragaria proteins and far apart from other F. vesca PIP proteins. A structural model for Fv PIP2-1a protein was built by comparative modelling methodology, which was validated and refined by molecular dynamics simulation. Fv PIP2-1a structure consists of 6 transmembrane regions and two NPA domains. The mobilization of water was analyzed by molecular docking simulations in wildtype and two mutants. Interestingly, the mutant FvPIP2-1a_H214G allowed the prediction of an increment in the flux of water molecules. On contrary, structural analysis predicted that H214E mutation blocked passage of water associated to constriction of the pore.

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.

Identification and comprehensive analysis of MIPSs in Rosaceae and their expression under abiotic stresses in rose (Rosa chinensis)

The Myo-Inositol-1-phosphate synthase (MIPS) gene family is involved in the myo-inositol synthesis and plays a significant role in signal transduction, membrane biogenesis, oligosaccharides synthesis, auxin storage and transport, programmed cell death and abiotic stress tolerance in plants. This study comprehensively identified the MIPS genes in Rosaceae plant species, and 51 MIPS genes were identified from 26 Rosaceae species. The phylogenetic analysis divided the MIPSs into two clades (clade I; subfamily Amygdaloideae specific, and clade II; subfamily Rosoideae specific). MIPS genes of all 26 Rosaceae species consist of similar gene structure, motif and domain composition, which shows their conserved nature. The cis-regulatory elements (CREs) analysis revealed that most Rosaceae MIPS genes play a role in growth, development, and stress responses. Furthermore, the qRT-PCR analysis also revealed the involvement of RcMIPS gene in plant development and response to abiotic stresses, including drought and heat. The results of the present study contribute to the understanding of the biological function of Rosaceae MIPS genes, and that could be used in further functional validations.

Spontaneous, Artificial, and Genome Editing-Mediated Mutations in Prunus

Mutation is a source of genetic diversity widely used in breeding programs for the acquisition of agronomically interesting characters in commercial varieties of the Prunus species, as well as in the rest of crop species. Mutation can occur in nature at a very low frequency or can be induced artificially. Spontaneous or bud sport mutations in somatic cells can be vegetatively propagated to get an individual with the mutant phenotype. Unlike animals, plants have unlimited growth and totipotent cells that let somatic mutations to be transmitted to the progeny. On the other hand, in vitro tissue culture makes it possible to induce mutation in plant material and perform large screenings for mutant’s selection and cleaning of chimeras. Finally, targeted mutagenesis has been boosted by the application of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas9 and Transcription activator-like effector nuclease (TALEN) editing technologies. Over the last few decades, environmental stressors such as global warming have been threatening the supply of global demand for food based on population growth in the near future. For this purpose, the release of new varieties adapted to such changes is a requisite, and selected or generated Prunus mutants by properly regulated mechanisms could be helpful to this task. In this work, we reviewed the most relevant mutations for breeding traits in Prunus species such as flowering time, self-compatibility, fruit quality, and disease tolerance, including new molecular perspectives in the present postgenomic era including CRISPR/Cas9 and TALEN editing technologies.

MdMYB52 regulates lignin biosynthesis upon the suberization process in apple

Our previous studies, comparing russeted vs. waxy apple skin, highlighted a MYeloBlastosys (Myb) transcription factor (MdMYB52), which displayed a correlation with genes associated to the suberization process. The present article aims to assess its role and function in the suberization process. Phylogenetic analyses and research against Arabidopsis thaliana MYBs database were first performed and the tissue specific expression of MdMYB52 was investigated using RT-qPCR. The function of MdMYB52 was further investigated using Agrobacterium-mediated transient overexpression in Nicotiana benthamiana leaves. An RNA-Seq analysis was performed to highlight differentially regulated genes in response MdMYB52. Transcriptomic data were supported by analytical chemistry and microscopy. A massive decreased expression of photosynthetic and primary metabolism pathways was observed with a concomitant increased expression of genes associated with phenylpropanoid and lignin biosynthesis, cell wall modification and senescence. Interestingly key genes involved in the synthesis of suberin phenolic components were observed. The analytical chemistry displayed a strong increase in the lignin content in the cell walls during MdMYB52 expression. More specifically, an enrichment in G-Unit lignin residues was observed, supporting transcriptomic data as well as previous work describing the suberin phenolic domain as a G-unit enriched lignin-like polymer. The time-course qPCR analysis revealed that the observed stress response, might be explain by this lignin biosynthesis and by a possible programmed senescence triggered by MdMYB52. The present work supports a crucial regulatory role for MdMYB52 in the biosynthesis of the suberin phenolic domain and possibly in the fate of suberized cells in russeted apple skins.

Integration of genome and transcriptome reveal molecular regulation mechanism of early flowering trait in Prunus genus (Prunus mume and Prunus persica)

Flowering time is crucial for the survival and reproduction. Prunus genus belongs to the Rosaceae family and includes several hundred species of flowering trees and shrubs with important ornamental and economic values. However, the molecular mechanism underlying early flowering in Prunus genus is unclear. Here, we utilized the genome and transcriptome of P. mume and P. persica to explore the transcriptional regulation mechanism of early flowering. Comparative genomics found that genes accounting for 92.4% of the total P. mume genome and 91.2% of the total P. persica genome belonged to orthogroups. A total of 19,169 orthogroups were found between P. mume and P. persica, including 20,431 corresponding orthologues and 20,080 collinearity gene pairs. A total of 305 differentially expressed genes (DEGs) associated with early flowering were found, among which FT, TLI65, and NAP57 were identified as hub genes in the early flowering regulation pathway. Moreover, we identified twenty-five transcription factors (TFs) from nine protein families, including MADS-box, AP2/ERF, and MYB. Our results provide insights into the underlying molecular model of flowering time regulation in Prunus genus and highlight the utility of multi-omics in deciphering the properties of the inter-genus plants.

Identification of Key Genes Related to Dormancy Control in Prunus Species by Meta-Analysis of RNAseq Data

Bud dormancy is a genotype-dependent mechanism observed in Prunus species in which bud growth is inhibited, and the accumulation of a specific amount of chilling (endodormancy) and heat (ecodormancy) is necessary to resume growth and reach flowering. We analyzed publicly available transcriptome data from fifteen cultivars of four Prunus species (almond, apricot, peach, and sweet cherry) sampled at endo- and ecodormancy points to identify conserved genes and pathways associated with dormancy control in the genus. A total of 13,018 genes were differentially expressed during dormancy transitions, of which 139 and 223 were of interest because their expression profiles correlated with endo- and ecodormancy, respectively, in at least one cultivar of each species. The endodormancy-related genes comprised transcripts mainly overexpressed during chilling accumulation and were associated with abiotic stresses, cell wall modifications, and hormone regulation. The ecodormancy-related genes, upregulated after chilling fulfillment, were primarily involved in the genetic control of carbohydrate regulation, hormone biosynthesis, and pollen development. Additionally, the integrated co-expression network of differentially expressed genes in the four species showed clusters of co-expressed genes correlated to dormancy stages and genes of breeding interest overlapping with quantitative trait loci for bloom time and chilling and heat requirements.

CarrotOmics: a genetics and comparative genomics database for carrot (Daucus carota)

CarrotOmics (https://carrotomics.org/) is a comprehensive database for carrot (Daucus carota L.) breeding and research. CarrotOmics was developed using resources available at the MainLab Bioinformatics core (https://www.bioinfo.wsu.edu/) and is implemented using Tripal with Drupal modules. The database delivers access to download or visualize the carrot reference genome with gene predictions, gene annotations and sequence assembly. Other genomic resources include information for 11 224 genetic markers from 73 linkage maps or genotyping-by-sequencing and descriptions of 371 mapped loci. There are records for 1601 Apiales species (or subspecies) and descriptions of 9408 accessions from 11 germplasm collections representing more than 600 of these species. Additionally, 204 Apiales species have phenotypic information, totaling 28 517 observations from 10 041 biological samples. Resources on CarrotOmics are freely available, search functions are provided to find data of interest and video tutorials are available to describe the search functions and genomic tools. CarrotOmics is a timely resource for the Apiaceae research community and for carrot geneticists developing improved cultivars with novel traits addressing challenges including an expanding acreage in tropical climates, an evolving consumer interested in sustainably grown vegetables and a dynamic environment due to climate change. Data from CarrotOmics can be applied in genomic-assisted selection and genetic research to improve basic research and carrot breeding efficiency.

Database URL

https://carrotomics.org/

Genome-Wide Analysis of the SPL Gene Family and Expression Analysis during Flowering Induction in Prunus × yedoensis 'Somei-yoshino'

SQUAMOSA Promoter-Binding Protein-Like (SPL) genes encode plant-specific transcription factors which bind to the SQUAMOSA promoter of the MADS-box genes to regulate its expression. It plays important regulatory roles in floral induction and development, fertility, light signals and hormonal transduction, and stress response in plants. In this study, 32 PySPL genes with complete SBP (squamosa promoter binding protein) conserved domain were identified from the genome of Prunus × yedoensis ‘Somei-yoshino’ and analyzed by bioinformatics. 32 PySPLs were distributed on 13 chromosomes, encoding 32 PySPL proteins with different physical and chemical properties. The phylogenetic tree constructed with Arabidopsis thaliana and Oryza sativa can be divided into 10 subtribes, indicating PySPLs of different clusters have different biological functions. The conserved motif prediction showed that the number and distribution of motifs on each PySPL is varied. The gene structure analysis revealed that PySPLs harbored exons ranging from 2 to 10. The predictive analysis of acting elements showed that the promoter of PySPLs contain a large number of light-responsive elements, as well as response elements related to hormone response, growth and development and stress response. The analysis of the PySPLs expressions in flower induction and flower organs based on qRT-PCR showed that PySPL06/22 may be the key genes of flower development, PySPL01/06 and PySPL22 may play a role in the development of sepal and pistil, respectively. The results provide a foundation for the study of SPL transcription factors of Prunus × yedoensis ‘Somei-yoshino’ and provide more reference information of the function of SPL gene in flowering.

Towards smart and sustainable development of modern berry cultivars in Europe

Fresh berries are a popular and important component of the human diet. The demand for high-quality berries and sustainable production methods is increasing globally, challenging breeders to develop modern berry cultivars that fulfill all desired characteristics. Since 1994, research projects have characterized genetic resources, developed modern tools for high-throughput screening, and published data in publicly available repositories. However, the key findings of different disciplines are rarely linked together, and only a limited range of traits and genotypes has been investigated. The Horizon2020 project BreedingValue will address these challenges by studying a broader panel of strawberry, raspberry and blueberry genotypes in detail, in order to recover the lost genetic diversity that has limited the aroma and flavor intensity of recent cultivars. We will combine metabolic analysis with sensory panel tests and surveys to identify the key components of taste, flavor and aroma in berries across Europe, leading to a high-resolution map of quality requirements for future berry cultivars. Traits linked to berry yields and the effect of environmental stress will be investigated using modern image analysis methods and modeling. We will also use genetic analysis to determine the genetic basis of complex traits for the development and optimization of modern breeding technologies, such as molecular marker arrays, genomic selection and genome-wide association studies. Finally, the results, raw data and metadata will be made publicly available on the open platform Germinate in order to meet FAIR data principles and provide the basis for sustainable research in the future.

QTL mapping and characterization of black spot disease resistance using two multi-parental diploid rose populations

Black spot disease (BSD) (Diplocarpon rosae) is the most common and damaging fungal disease in garden roses (Rosa sp.). Although qualitative resistance to BSD has been extensively investigated, the research on quantitative resistance lags behind. The goal of this research was to study the genetic basis of BSD resistance in two multi-parental populations (TX2WOB and TX2WSE) through a pedigree-based analysis approach (PBA). Both populations were genotyped and evaluated for BSD incidence over five years in three locations in Texas. A total of 28 QTLs, distributed over all linkage groups (LGs), were detected across both populations. Consistent minor effect QTLs included two on LG1 and LG3 (TX2WOB and TX2WSE), two on LG4 and LG5 (TX2WSE), and one QTL on LG7 (TX2WOB). In addition, one major QTL detected in both populations was consistently mapped on LG3. This QTL was localized to an interval ranging from 18.9 to 27.8 Mbp on the Rosa chinensis genome and explained 20 and 33% of the phenotypic variation. Furthermore, haplotype analysis showed that this QTL had three distinct functional alleles. The parent PP-J14-3 was the common source of the LG3 BSD resistance in both populations. Taken together, this research presents the characterization of new SNP-tagged genetic determinants of BSD resistance, the discovery of marker-trait associations to enable parental choice based on their BSD resistance QTL haplotypes, and substrates for the development of trait-predictive DNA tests for routine use in marker-assisted breeding for BSD resistance.

Inheritance of esters and other volatile compounds responsible for the fruity aroma in strawberry

Cultivated strawberry, Fragaria  ×  ananassa, has a complex aroma due to the presence of more than 350 volatile organic compounds (VOCs). However, a mixture of only 19 compounds, called Key Volatile Compounds (KVC), can impart the main strawberry aroma. The octoploid nature of the cultivated strawberry species (2n = 8x = 56) adds complexity to the heritance of the accumulation of the volatiles responsible for aroma. An F1 population cross between two breeding parental lines, FC50 and FD54, was phenotyped for aroma by SPME GCMS during six harvests. A total of 58 compounds were identified: 33 esters, nine terpenes, seven aldehydes, four lactones, two furans, one acid, one alkane and one alcohol, of which 16 were KVCs. A total of 179 QTLs were found, and 85 of these were detected in at least three harvests, of which 50 QTLs were considered major (LOD > 4.0) and detected in five or six analyzed harvests. Several clusters of ester QTLs associated with fruity aroma were discovered, such as QTLs for esters that share hexanoate group that were mapped in LG4A (Hexanoate_4A), those that share acetate and octyl groups in LG6A (Acetate_6A and Octyl_6A) or those with the same methyl group in LG7B (Methyl_7B). Different terpene QTLs associated with floral aroma appear grouped in a cluster in LG3C (Terpene_3C). Some of these clusters of QTLs were validated in a second F2 population, a cross of "Camarosa" and "Dover," that was also phenotyped for three years. Selected SNPs from floral and fruity aroma QTLs were tested in a third population, which will most likely be useful for marker-assisted breeding (MAB).

Genomic and Bioinformatic Resources for Perennial Fruit Species

In the post-genomic era, data management and development of bioinformatic tools are critical for the adequate exploitation of genomics data. In this review, we address the actual situation for the subset of crops represented by the perennial fruit species. The agronomical singularity of these species compared to plant and crop model species provides significant challenges on the implementation of good practices generally not addressed in other species. Studies are usually performed over several years in non-controlled environments, usage of rootstock is common, and breeders heavily rely on vegetative propagation. A reference genome is now available for all the major species as well as many members of the economically important genera for breeding purposes. Development of pangenome for these species is beginning to gain momentum which will require a substantial effort in term of bioinformatic tool development. The available tools for genome annotation and functional analysis will also be presented.

Signatures of polygenic adaptation align with genome-wide methylation patterns in wild strawberry plants

Epigenetic inheritance can drive adaptive evolution independently of DNA sequence variation. However, to what extent epigenetic variation represents an autonomous evolutionary force remains largely elusive. Through gene ontology and comparative analyses of genomic and epigenomic variation of wild strawberry plants raised in distinct drought settings, we characterised genome-wide covariation between single nucleotide polymorphisms (SNPs) and differentially methylated cytosines (DMCs). Covariation between SNPs and DMCs was independent of genomic proximity, but instead associated with fitness-related processes such as stress responses, genome regulation and reproduction. We expected this functional SNP-DMC covariation to be driven by adaptive evolution canalising SNP and DMC variation, but instead observed significantly lower covariation with DMCs for adaptive rather than for neutral SNPs. Drought-induced DMCs frequently co-varied with tens of SNPs, suggesting high genomic redundancy as a broad potential basis for polygenic adaptation of gene expression. Our findings suggest that stress-responsive DMCs initially co-vary with many SNPs under increased environmental stress, and that natural selection acting upon several of these SNPs subsequently reduces standing covariation with stress-responsive DMCs. Our study supports DNA methylation profiles that represent complex quantitative traits rather than autonomous evolutionary forces. We provide a conceptual framework for polygenic regulation and adaptation shaping genome-wide methylation patterns in plants.

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.

Chromosome-scale assembly with a phased sex-determining region resolves features of early Z and W chromosome differentiation in a wild octoploid strawberry

When sex chromosomes stop recombining, they start to accumulate differences. The sex-limited chromosome (Y or W) especially is expected to degenerate via the loss of nucleotide sequence and the accumulation of repetitive sequences. However, how early signs of degeneration can be detected in a new sex chromosome is still unclear. The sex-determining region of the octoploid strawberries is young, small, and dynamic. Using PacBio HiFi reads, we obtained a chromosome-scale assembly of a female (ZW) Fragaria chiloensis plant carrying the youngest and largest of the known sex-determining region on the W in strawberries. We fully characterized the previously incomplete sex-determining region, confirming its gene content, genomic location, and evolutionary history. Resolution of gaps in the previous characterization of the sex-determining region added 10 kb of sequence including a noncanonical long terminal repeat-retrotransposon; whereas the Z sequence revealed a Harbinger transposable element adjoining the sex-determining region insertion site. Limited genetic differentiation of the sex chromosomes coupled with structural variation may indicate an early stage of W degeneration. The sex chromosomes have a similar percentage of repeats but differ in their repeat distribution. Differences in the pattern of repeats (transposable element polymorphism) apparently precede sex chromosome differentiation, thus potentially contributing to recombination cessation as opposed to being a consequence of it.

Mechanism of fertilization-induced auxin synthesis in the endosperm for seed and fruit development

The dominance of flowering plants on earth is owed largely to the evolution of maternal tissues such as fruit and seedcoat that protect and disseminate the seeds. The mechanism of how fertilization triggers the development of these specialized maternal tissues is not well understood. A key event is the induction of auxin synthesis in the endosperm, and the mobile auxin subsequently stimulates seedcoat and fruit development. However, the regulatory mechanism of auxin synthesis in the endosperm remains unknown. Here, we show that a type I MADS box gene AGL62 is required for the activation of auxin synthesis in the endosperm in both Fragaria vesca, a diploid strawberry, and in Arabidopsis. Several strawberry FveATHB genes were identified as downstream targets of FveAGL62 and act to repress auxin biosynthesis. In this work, we identify a key mechanism for auxin induction to mediate fertilization success, a finding broadly relevant to flowering plants.

In flowering plants, fertilization triggers auxin synthesis in the endosperm to promote seed and fruit development. Here the authors show that an MADS-box transcription factor AGL62 is required to activate auxin synthesis in the endosperms of Fragaria vesca, a diploid strawberry, and Arabidopsis.

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.

Effect of Developmental Stages on Genes Involved in Middle and Downstream Pathway of Volatile Terpene Biosynthesis in Rose Petals

Terpenoids are economically and ecologically important compounds, and they are vital constituents in rose flower fragrance and rose essential oil. The terpene synthase genes (TPSs), trans-prenyltransferases genes (TPTs), NUDX1 are involved in middle and downstream pathway of volatile terpene biosynthesis in rose flowers. We identified 7 complete RcTPTs, 49 complete RcTPSs, and 9 RcNUDX1 genes in the genome of Rosachinensis. During the flower opening process of butterfly rose (Rosachinensis ‘Mutabilis’, MU), nine RcTPSs expressed in the petals of opening MU flowers exhibited two main expression trends, namely high and low, in old and fresh petals. Five short-chain petal-expressed RcTPTs showed expression patterns corresponding to RcTPSs. Analysis of differential volatile terpenes and differential expressed genes indicated that higher emission of geraniol from old MU petals might be related to the RcGPPS expression. Comprehensive analysis of volatile emission, sequence structure, micro-synteny and gene expression suggested that RcTPS18 may encode (E,E)-α-farnesene synthase. These findings may be useful for elucidating the molecular mechanism of terpenoid metabolism in rose and are vital for future studies on terpene regulation.