Comparative and Expression Analysis of Ubiquitin Conjugating Domain-Containing Genes in Two Pyrus Species

The chromosome-level genome assembly of the dwarfing apple interstock Malus hybrid 'SH6'

Malus hybrid 'SH6' (M. honanensis × M. domestica)is a commonly used apple interstock in China, known for its excellent dwarfing characteristics and cold tolerance. In this study, a combined strategy utilizing PacBio HiFi, Hi-C and parental resequencing data were employed to assemble two haploid genomes for 'SH6'. After chromosome anchoring, the final hapH genome size was 596.63 Mb, with a contig N50 of 34.38 Mb. The hapR genome was 649.37 Mb, with a contig N50 of 36.84 Mb. Further analysis predicted that repeated sequences made up 59.69% and 62.52% of the entire genome, respectively. Gene annotations revealed 45,435 genes for hapH and 48,261 genes for hapR. Combined with genomic synteny we suggest that the hapR genome originates from its maternal parent M. domestica cv. Ralls Janet, while the hapH genome comes from its paternal parent, M. honanensis. The assembled genome significantly contributes to the discovery of genes associated with apple dwarfing and the molecular mechanisms governing them.

Evolution and function of ubiquitin-specific proteases (UBPs): Insight into seed development roles in tung tree (Vernicia fordii)

The tung oil produced by the tung tree (Vernicia fordii) provides resources for the manufacture of biodiesel. Ubiquitin-specific proteases (UBPs) are the largest group of deubiquitinases and play key roles in regulating development and stress responses. Here, 21 UBPs were identified in V. fordii, roughly one-half the number found in Manihot esculenta and Hevea brasiliensis. Most UBP duplications are produced from whole-genome duplication (WGD), and significant differences in gene retention existed among Euphorbiaceae. The great majority of UBP-containing blocks in V. fordii, V. montana, Ricinus communis, and Jatropha curcas exhibited extensive conservation with the duplicated regions of M. esculenta and H. brasiliensis. These blocks formed 14 orthologous groups, indicating they shared WGD with UBPs in M. esculenta and H. brasiliensis, but most of these UBPs copies were lost. The UBP orthologs contained significant functional divergence which explained the susceptibility of V. fordii to Fusarium wilt and the resistance of V. montana to Fusarium wilt. The expression patterns and experiments suggested that Vf03G1417 could affect the seed-related traits and positively regulate the seed oil accumulation. This study provided important insights into the evolution of UBPs in Euphorbiaceae and identified important candidate VfUBPs for marker-assisted breeding in V. fordii.

Whole-genome analysis of CGS, SAHH, SAMS gene families in five Rosaceae species and their expression analysis in Pyrus bretschneideri

Cystathionine γ-synthase (CGS), S-adenosyl-L-homocysteine hydrolase (SAHH), and S-adenosy-L-methionine synthetase (SAMS) play an important role in the regulation of plant growth, development, and secondary metabolism. In this study, a total of 6 CGS, 6 SAHH, and 28 SAMS genes were identified from five Rosaceae species (Pyrus bretschneideri, Prunus persica, Prunus mume, Fragaria vesca, and Malus domestica). The evolutionary relationship and microsynteny analysis in five Rosaceae species revealed that duplicated regions were conserved between three gene families (CGS, SAHH, SAMS). Moreover, the chromosomal locations, gene structures, conserved motifs, cis-elements, physicochemical properties, and Ka/Ks analysis were performed by using numerous bioinformatics tools. The expression of different organs showed that the CGS, SAHH and SAMS genes of pear have relatively high expression patterns in flowers and stems, except for PbCGS1. RNA-seq and qRT-PCR combined analysis showed that PbSAMS1 may be involved in the regulation of pear stone cell development. In summary, this study provides the basic information of CGS, SAHH and SAMS genes in five Rosaceae species, further revealing the expression patterns in the pear fruit, which provides the theoretical basis for the regulation of pear stone cells.

Characterization and regulation mechanism analysis of ubiquitin-conjugating family genes in strawberry reveals a potential role in fruit ripening

BACKGROUND

E2 ubiquitin-conjugating (UBC) enzymes are an integral component of the ubiquitin proteasome system that play an important role in plant development, growth, and external stress responses. Several UBC genes have been identified in various plants. However, no studies exploring the functions of UBC genes in regulating fruit of strawberry have been reported. In the present study, a systematic analysis of the entire UBC family members were conducted in the genome of strawberry (Fragaria ×ananassa) based on bioinformatics method, and the gene functioning in strawberry ripening was explored.

RESULTS

A total of 191 UBC genes were identified in the genome of cultivated strawberry. These genes were unevenly distributed across the 28 chromosomes from the 4 subgenomes of cultivated strawberry, ranging from 3 to 11 genes per chromosome. Moreover, the expansion of FaUBC genes in strawberry was mainly driven by WGD. All the FaUBC genes were clarified into 13 groups and most of them were included in the group VI. The gene structure analysis showed that the number of exons varied from 1 to 23, and the structure of genes had few differences within the same groups but a distinction in different groups. Identification of the cis-acting elements of the promoter revealed multiple regulatory elements that responded to plant growth and development, phytohormone responsive, and abiotic and biotic stress. Data from functional annotation indicated that FaUBC genes play a role in a variety of biological processes. The RNA-seq data showed that FaUBC genes displayed different expression pattern during the fruit ripening process and clarified into 6 clusters. In particular, cluster 3 exhibiting a sudden expression increase in the turning red stage were speculated to be involved in fruit ripening. Hence, two FaUBC genes (FaUBC76 and FaUBC78) were selected for gene function analysis by transient over-expression method. The results indicated that FaUBC76 has a positive effect on the fruit development and ripening in strawberry by up-regulating accumulation of anthocyanins. Moreover, expression of some maturity-related genes were also significantly increased, further supporting a role for FaUBC76 in the regulation of fruit ripening or softening. On the contrary, the overexpression of FaUBC78 significantly increased the firmness of strawberry fruit, indicating that FaUBC78 had a positive role in inhibiting the decrease of strawberry fruit firmness.

CONCLUSION

Our study not only provide comprehensive information on system evolution and function on UBC genes, but also give a new insight into explore the roles of FaUBC genes in the regulation of strawberry ripening.

Cloning and functional characterization of two cinnamate 4-hydroxylase genes from Pyrus bretschneideri

Cinnamate 4-hydroxylase (C4H) is a key enzyme in the phenylpropanoid pathway in plants and is involved in the biosynthesis of secondary metabolites such as lignin and flavonoids. However, the function of C4H in pear plants (Pyrus bretschneideri) has not yet been fully elucidated. By searching pear genome databases, we identified three C4H genes (PbC4H1, PbC4H2 and PbC4H3) encoding proteins that share higher identity with bonafide C4Hs from several species with typical cytochrome P450 domains, suggesting that all three PbC4Hs are also bonafide C4Hs that have close evolutionary relationships with C4Hs from other land plants. Quantitative real-time PCR (qRT-PCR) results indicated that the three PbC4Hs were specifically expressed in one or more tissues. The expression levels of PbC4H1 and PbC4H3 first increased and then decreased during pear fruit development. Treatment with exogenous hormones (ABA, MeJA, and SA) altered the expression of the three PbC4Hs to varying degrees. The expression levels of the PbC4Hs were first induced and then decreased under ABA treatment, while MeJA treatment significantly increased the expression levels of the PbC4Hs. Following treatment with SA, expression levels of PbC4H1 and PbC4H2 increased, while expression levels of PbC4H3 decreased. Enzymatic analysis of the recombinant proteins expressed in yeast indicated that PbC4H1 and PbC4H3 catalysed the conversion of trans-cinnamic acid to p-coumaric acid. Moreover, the expression of PbC4H1 and PbC4H3 in Arabidopsis resulted in an increase in both the lignin content and the thickness of cell walls for intervascular fibres and xylem cells. Taken together, the results of our study not only revealed the potential role of PbC4H1 and PbC4H3 in lignin biosynthesis but also established a foundation for future investigations of the regulation of lignin synthesis and stone cell development in pear fruit by molecular biological techniques.

Comprehensive Proteomic Analysis of Lysine Ubiquitination in Seedling Leaves of Nicotiana tabacum

Lysine ubiquitination, a widely studied posttranslational modification, plays vital roles in various biological processes in eukaryotic cells. Although several studies have examined the plant ubiquitylome, no such research has been performed in tobacco, a model plant for molecular biology. Here, we comprehensively analyzed lysine ubiquitination in tobacco (Nicotiana tabacum) using LC-MS/MS along with highly sensitive immune-affinity purification. In total, 964 lysine-ubiquitinated (K^ub) sites were identified in 572 proteins. Extensive bioinformatics studies revealed the distribution of these proteins in various cellular locations, including the cytoplasm, chloroplast, nucleus, and plasma membrane. Notably, 25% of the K^ub proteins were located in the chloroplast of which 21 were enzymatically involved in important pathways, that is, photosynthesis and carbon fixation. Western blot analysis indicated that TMV infection can cause changes in ubiquitination levels. This is the first comprehensive proteomic analysis of lysine ubiquitination in tobacco, illustrating the vital role of ubiquitination in various physiological and biochemical processes and representing a valuable addition to the existing landscape of lysine ubiquitination.

Functional divergence and adaptive selection of KNOX gene family in plants

KNOTTED-like homeodomain (KNOX) genes are transcriptional regulators that play an important role in morphogenesis. In the present study, a comparative analysis was performed to investigate the molecular evolution of the characteristics of the KNOX gene family in 10 different plant species. We identified 129 KNOX gene family members, which were categorized into two subfamilies based on multiple sequence alignment and phylogenetic tree reconstruction. Several segmental duplication pairs were found, indicating that different species share a common expansion model. Functional divergence analysis identified the 15 and 52 amino acid sites with significant changes in evolutionary rates and amino acid physicochemical properties as functional divergence sites. Additional selection analysis showed that 14 amino acid sites underwent positive selection during evolution, and two groups of co-evolutionary amino acid sites were identified by Coevolution Analysis using Protein Sequences software. These sites could play critical roles in the molecular evolution of the KNOX gene family in these species. In addition, the expression profiles of KNOX duplicated genes demonstrated functional divergence. Taken together, these results provide novel insights into the structural and functional evolution of the KNOX gene family.

A Chinese White Pear (Pyrus bretschneideri) BZR Gene PbBZR1 Act as a Transcriptional Repressor of Lignin Biosynthetic Genes in Fruits

BZR transcription factors play essential roles in plant growth and environmental stimuli, and they are also the positive regulators of Brassinosteroid (BR) signal transduction in diverse plants. In addition, BZR TFs, as crucial regulators of BR synthesis, may have multiple stress-resistance functions and their related regulatory mechanisms have been well illustrated in model plants. Here, we carried out a genome-wide identification of BZR members in Chinese pear (Pyrus bretschneideri) and identified 13 members. By comparative analysis in five Rosaceae genomes, BZR members in the pear genome may have undergone large-scale duplication events during evolution. Purifying selection played an important role in almost all of the orthologous and paralogous gene pairs. According to the expression analysis of the PbBZRs during fruit development, three PbBZRs were selected for detailed analysis. Transcriptional activation assays presented that PbBZR1 repressed the promoters of P. bretschneideri lignin biosynthetic genes, such as PbCES9, PbCOMT3, and PbHCT6. Our study traces the evolution of BZR gene family members in Rosaceae genomes and illustrates that the rates of gene loss and gain are far from equilibrium in different species. At the same time, our results suggest that PbBZR1 may be involved in the negative regulation of lignin biosynthesis.

A de novo genome assembly of the dwarfing pear rootstock Zhongai 1

'Zhongai 1' [(Pyrus ussuriensis × communis) × spp.] is an excellent pear dwarfing rootstock common in China. It is dwarf itself and has high dwarfing efficiency on most of main Pyrus cultivated species when used as inter-stock. Here we describe the draft genome sequences of 'Zhongai 1' which was assembled using PacBio long reads, Illumina short reads and Hi-C technology. We estimated the genome size is approximately 511.33 Mb by K-mer analysis and obtained a final genome of 510.59 Mb with a contig N50 size of 1.28 Mb. Next, 506.31 Mb (99.16%) of contigs were clustered into 17 chromosomes with a scaffold N50 size of 23.45 Mb. We further predicted 309.86 Mb (60.68%) of repetitive sequences and 43,120 protein-coding genes. The assembled genome will be a valuable resource and reference for future pear breeding, genetic improvement, and comparative genomics among related species. Moreover, it will help identify genes involved in dwarfism, early flowering, stress tolerance, and commercially desirable fruit characteristics.

Metacaspase gene family in Rosaceae genomes: Comparative genomic analysis and their expression during pear pollen tube and fruit development

Metacaspase (MC), which is discovered gene family with distant caspase homologs in plants, fungi, and protozoa, may be involved in programmed cell death (PCD) processes during plant development and respond abiotic and biotic stresses. To reveal the evolutionary relationship of MC gene family in Rosaceae genomes, we identified 8, 7, 8, 12, 12, and 23 MC genes in the genomes of Fragaria vesca, Prunus mume, Prunus persica, Pyrus communis, Pyrus bretschneideri and Malus domestica, respectively. Phylogenetic analysis suggested that the MC genes could be grouped into three clades: Type I*, Type I and Type II, which was supported by gene structure and conserved motif analysis. Microsynteny analysis revealed that MC genes present in the corresponding syntenic blocks of P. communis, P. bretschneideri and M. domestica, and further suggested that large-scale duplication events play an important role in the expansion of MC gene family members in these three genomes than other Rosaceae plants (F. vesca, P. mume and P. persica). RNA-seq data showed the specific expression patterns of PbMC genes in response to drought stress. The expression analysis of MC genes demonstrated that PbMC01 and PbMC03 were able to be detected in all four pear pollen tubes and seven fruit development stages. The current study highlighted the evolutionary relationship and duplication of the MC gene family in these six Rosaceae genomes and provided appropriate candidate genes for further studies in P. bretschneideri.

Evolutionary Analysis of GH3 Genes in Six Oryza Species/Subspecies and Their Expression under Salinity Stress in Oryza sativa ssp. japonica

Glycoside Hydrolase 3 (GH3), a member of the Auxin-responsive gene family, is involved in plant growth, the plant developmental process, and various stress responses. The GH3 gene family has been well-studied in Arabidopsis thaliana and Zea mays. However, the evolution of the GH3 gene family in Oryza species remains unknown and the function of the GH3 gene family in Oryza sativa is not well-documented. Here, a systematic analysis was performed in six Oryza species/subspecies, including four wild rice species and two cultivated rice subspecies. A total of 13, 13, 13, 13, 12, and 12 members were identified in O. sativa ssp. japonica, O. sativa ssp. indica, Oryza rufipogon, Oryza nivara, Oryza punctata, and Oryza glumaepatula, respectively. Gene duplication events, structural features, conserved motifs, a phylogenetic analysis, chromosome locations, and Ka/Ks ratios of this important family were found to be strictly conservative across these six Oryza species/subspecies, suggesting that the expansion of the GH3 gene family in Oryza species might be attributed to duplication events, and this expansion could occur in the common ancestor of Oryza species, even in common ancestor of rice tribe (Oryzeae) (23.07~31.01 Mya). The RNA-seq results of different tissues displayed that OsGH3 genes had significantly different expression profiles. Remarkably, the qRT-PCR result after NaCl treatment indicated that the majority of OsGH3 genes play important roles in salinity stress, especially OsGH3-2 and OsGH3-8. This study provides important insights into the evolution of the GH3 gene family in Oryza species and will assist with further investigation of OsGH3 genes’ functions under salinity stress.

Identification and Characterization of LRR-RLK Family Genes in Potato Reveal Their Involvement in Peptide Signaling of Cell Fate Decisions and Biotic/Abiotic Stress Responses

Leucine-rich repeat receptor-like kinases (LRR-RLKs) represent the largest subfamily of receptor-like kinases (RLKs) and play important roles in regulating growth, development, and stress responses in plants. In this study, 246 LRR-RLK genes were identified in the potato (Solanum tuberosum) genome, which were further classified into 14 subfamilies. Gene structure analysis revealed that genes within the same subgroup shared similar exon/intron structures. A signature small peptide recognition motif (RxR) was found to be largely conserved within members of subfamily IX, suggesting that these members may recognize peptide signals as ligands. 26 of the 246 StLRR-RLK genes were found to have arisen from tandem or segmental duplication events. Expression profiling revealed that StLRR-RLK genes were differentially expressed in various organs/tissues, and several genes were found to be responsive to different stress treatments. Furthermore, StLRR-RLK117 was found to be able to form homodimers and heterodimers with StLRR-RLK042 and StLRR-RLK052. Notably, the overlapping expression region of StLRR-RLK117 with Solanum tuberosumWUSCHEL (StWUS) suggested that the CLV3–CLV1/BAM–WUS feedback loop may be conserved in potato to maintain stem cell homeostasis within the shoot apical meristem.