The Botrytis cinerea effector BcXYG1 suppresses immunity in Fragaria vesca by targeting FvBPL4 and FvACD11.

Botrytis cinerea is one of the most destructive pathogens in strawberry cultivation. Successful infection by B. cinerea requires releasing a large number of effectors that interfere with the plant's immune system. One of the effectors required by B. cinerea for optimal virulence is the secreted protein BcXYG1, which is thought to associate with proteins near the plasma membrane of the host plant to induce necrosis. However, the host proteins that associate with BcXYG1 at the plasma membrane are currently unknown. We found that BcXYG1 binds to FvBPL4 and FvACD11 at the plasma membrane. Both FvBPL4 and FvACD11 are negative regulators of plant immunity in strawberry. Our results demonstrate that degradation of FvBPL4 by BcXYG1 promotes disease resistance while stabilization of FvACD11 by BcXYG1 suppresses the immune response. These findings suggest that BcXYG1 suppresses plant immunity and promotes B. cinerea infection by regulating FvBPL4 and FvACD11 protein levels.

The high-quality genome of Cryptotaenia japonica and comparative genomics analysis reveals anthocyanin biosynthesis in Apiaceae.

Cryptotaenia japonica, a traditional medicinal and edible vegetable crops, is well-known for its attractive flavors and health care functions. As a member of the Apiaceae family, the evolutionary trajectory and biological properties of C. japonica are not clearly understood. Here, we first reported a high-quality genome of C. japonica with a total length of 427 Mb and N50 length 50.76 Mb, was anchored into 10 chromosomes, which confirmed by chromosome (cytogenetic) analysis. Comparative genomic analysis revealed C. japonica exhibited low genetic redundancy, contained a higher percentage of single-cope gene families. The homoeologous blocks, Ks, and collinearity were analyzed among Apiaceae species contributed to the evidence that C. japonica lacked recent species-specific WGD. Through comparative genomic and transcriptomic analyses of Apiaceae species, we revealed the genetic basis of the production of anthocyanins. Several structural genes encoding enzymes and transcription factor genes of the anthocyanin biosynthesis pathway in different species were also identified. The CjANSa, CjDFRb, and CjF3H gene might be the target of Cjaponica_2.2062 (bHLH) and Cjaponica_1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology.

Transcriptome and Metabolome Reveal Sugar and Organic Acid Accumulation in Rosa roxburghii Fruit.

Sugars and organic acids significantly impact fruit sensory quality, but their accumulation patterns and regulatory mechanisms during the development of Rosa roxburghii fruit are still unclear. We utilized transcriptomics and metabolomics to investigate genes related to sugar and organic acid metabolism in Rosa roxburghii. Metabolomics data revealed that sucrose, glucose and fructose were the primary sugars, whereas citric acid and malic acid were the primary organic acids in Rosa roxburghii fruit. We constructed the metabolic pathways of major sugars and organic acids in Rosa roxburghii and identified five key genes involved in sugar and organic acid synthesis. In addition, we identified a module containing 132 transcription factors that was significantly associated with sucrose, citric acid and malic acid. Based on quantitative polymerase chain reaction (qPCR), we identified 13 transcription factors involved in sugar and organic acid metabolism, including the transcription factor RrANL2 and the sucrose synthase gene RrSUS3. Further yeast one-hybrid and dual luciferase assays showed that RrANL2 could bind to the promoter of RrSUS3 to increase its expression. These results provide new insights into the metabolism of sugars and organic acids in Rosa roxburghii fruit.

Different evolutionary patterns of TIR1/AFBs and AUX/IAAs and their implications for the morphogenesis of land plants.

BACKGROUND: The plant hormone auxin is widely involved in plant growth, development, and morphogenesis, and the TIR1/AFB and AUX/IAA proteins are closely linked to rapid auxin response and signal transmission. However, their evolutionary history, historical patterns of expansion and contraction, and changes in interaction relationships are still unknown. RESULTS: Here, we analyzed the gene duplications, interactions, and expression patterns of TIR1/AFBs and AUX/IAAs to understand their underlying mechanisms of evolution. The ratios of TIR1/AFBs to AUX/IAAs range from 4:2 in Physcomitrium patens to 6:29 in Arabidopsis thaliana and 3:16 in Fragaria vesca. Whole-genome duplication (WGD) and tandem duplication have contributed to the expansion of the AUX/IAA gene family, but numerous TIR1/AFB gene duplicates were lost after WGD. We further analyzed the expression profiles of TIR1/AFBs and AUX/IAAs in different tissue parts of Physcomitrium patens, Selaginella moellendorffii, Arabidopsis thaliana and Fragaria vesca, and found that TIR1/AFBs and AUX/IAAs were highly expressed in all tissues in P. patens, S. moellendorffii. In A. thaliana and F. vesca, TIR1/AFBs maintained the same expression pattern as the ancient plants with high expression in all tissue parts, while AUX/IAAs appeared tissue-specific expression. In F. vesca, 11 AUX/IAAs interacted with TIR1/AFBs with different interaction strengths, and the functional specificity of AUX/IAAs was related to their ability to bind TIR1/AFBs, thus promoting the development of specific higher plant organs. Verification of the interactions among TIR1/AFBs and AUX/IAAs in Marchantia polymorpha and F. vesca also showed that the regulation of AUX/IAA members by TIR1/AFBs became more refined over the course of plant evolution. CONCLUSIONS: Our results indicate that specific interactions and specific gene expression patterns both contributed to the functional diversification of TIR1/AFBs and AUX/IAAs.

Development of a single-cell atlas for woodland strawberry (Fragaria vesca) leaves during early Botrytis cinerea infection using single cell RNA-seq.

Pathogen invasion leads to fast, local-to-systemic signal transduction that initiates plant defense responses. Despite tremendous progress in past decades, aspects of this process remain unknown, such as which cell types respond first and how signals are transferred among cell types. Here, we used single-cell RNA-seq of more than 50 000 single cells to document the gene expression landscape in leaves of woodland strawberry during infection by Botrytis cinerea and identify major cell types. We constructed a single-cell atlas and characterized the distinct gene expression patterns of hydathode, epidermal, and mesophyll cells during the incubation period of B. cinerea infection. Pseudotime trajectory analysis revealed signals of the transition from normal functioning to defense response in epidermal and mesophyll cells upon B. cinerea infection. Genes related to disease resistance showed different expression patterns among cell types: disease resistance-related genes and gene encoding transcription factors were highly expressed in individual cell types and interacted to trigger plant systemic immunity to B. cinerea. This is the first report to document the of single-cell transcriptional landscape of the plant pathogenic invasion process, it provides new insights into the wholistic dynamics of host-pathogen interactions and can guide the identification of genes and the formulation of strategies for resistant cultivar development.

A chromosome-level genome assembly of rugged rose (Rosa rugosa) provides insights into its evolution, ecology, and floral characteristics.

Rosa rugosa, commonly known as rugged rose, is a perennial ornamental shrub. It produces beautiful flowers with a mild fragrance and colorful seed pods. Unlike many other cultivated roses, R. rugosa adapts to a wide range of habitat types and harsh environmental conditions such as salinity, alkaline, shade, drought, high humidity, and frigid temperatures. Here, we produced and analyzed a high-quality genome sequence for R. rugosa to understand its ecology, floral characteristics and evolution. PacBio HiFi reads were initially used to construct the draft genome of R. rugosa, and then Hi-C sequencing was applied to assemble the contigs into 7 chromosomes. We obtained a 382.6 Mb genome encoding 39,704 protein-coding genes. The genome of R. rugosa appears to be conserved with no additional whole-genome duplication after the gamma whole-genome triplication (WGT), which occurred ~100 million years ago in the ancestor of core eudicots. Based on a comparative analysis of the high-quality genome assembly of R. rugosa and other high-quality Rosaceae genomes, we found a unique large inverted segment in the Chinese rose R. chinensis and a retroposition in strawberry caused by post-WGT events. We also found that floral development- and stress response signaling-related gene modules were retained after the WGT. Two MADS-box genes involved in floral development and the stress-related transcription factors DREB2A-INTERACTING PROTEIN 2 (DRIP2) and PEPTIDE TRANSPORTER 3 (PTR3) were found to be positively selected in evolution, which may have contributed to the unique ability of this plant to adapt to harsh environments. In summary, the high-quality genome sequence of R. rugosa provides a map for genetic studies and molecular breeding of this plant and enables comparative genomic studies of Rosa in the near future.

Genome-wide identification and expression analyses of Sm genes reveal their involvement in early somatic embryogenesis in Dimocarpus longan Lour.

The Sm proteins are a conserved protein family with Sm motifs. The family includes Sm and Sm-like proteins, which play important roles in pre-mRNA splicing. Most research on the Sm proteins have been conducted in herbaceous plants, and less in woody plants such as Dimocarpus longan (longan). And the embryo development status significantly affects the quality and yield of longan. In this study, we conducted a genome-wide analysis of longan Sm genes (DlSm) to clarify their roles during somatic embryogenesis (SE) and identified 29 Sm genes. Phylogenetic analysis deduced longan Sm proteins clustered into 17 phylogenetic groups with the homologous Sm proteins of Arabidopsis thaliana. We also analyzed the gene structures, motif compositions, and conserved domains of the longan Sm proteins. The promoter sequences of the DlSm genes contained many light, endosperm development, hormone, and temperature response elements, which suggested their possible functions. In the non-embryogenic callus(NEC) and during early SE in longan, the alternative splicing(AS) events of DlSm genes indicated that these genes may influence SE development by changing gene structures and sequences. The kinetin(KT) hormone, and blue and white light treatments affected the differentiation and growth of longan embryonic callus(EC) probably by affecting the AS events of DlSm genes. Expression profiles showed the possible functional divergence among Sm genes, and different hormones and light qualities affected their expression levels. The expression trends of the DlSm genes determined by RNA sequencing as fragments per kilobase of exon model per million mapped reads (FPKM) and by real-time quantitative PCR(qRT-PCR) during early SE in longan showed that the expression of the DlSm genes was affected by the growth and differentiation of longan SE, and decreased as the somatic embryo differentiation progressed. The results will contributed to understanding the longan Sm gene family and provide a basis for future functional validation studies.

Genome-wide identification of miRNAs and their targets during early somatic embryogenesis in Dimocarpus longan Lour.

miRNAs are endogenous regulatory factors that play pivotal roles in post-transcriptional regulation. However, their specific roles in early somatic embryogenesis (SE) remain unclear. Study of the SE system is fundamental for clarifying the molecular mechanisms in Dimocarpus longan. We identified 289 known miRNAs from 106 different miRNA families and 1087 novel miRNAs during early longan SE, including embryogenic callus (EC), incomplete pro-embryogenic culture (ICpEC), globular embryo (GE), and non-embryogenic callus (NEC). The abundances of known miRNAs were concentrated in GE. The differentially expression (DE) miRNAs showed five expression patterns during early SE. Largely miRNAs were expressed highly and specially in EC, ICpEC, and GE, respectively. Some miRNAs and putative target genes were enriched in lignin metabolism. Most potential targets were related to the pathways of plant hormone signal transduction, alternative splicing, tyrosine metabolism and sulfur metabolism in early longan SE. The regulatory relationships between dlo-miR166a-3p and DlHD-zip8, dlo-miR397a and DlLAC7, dlo-miR408-3p and DlLAC12 were confirmed by RNA ligase-mediated rapid amplification of cDNA ends. The expression patterns of eight DE miRNAs detected by qRT-PCR were consistent with RNA-seq. Finally, the miRNA regulatory network in early SE was constructed, which provided new insight into molecular mechanism of early SE in longan.

Genome-wide identification and characterization of long non-coding RNAs involved in the early somatic embryogenesis in Dimocarpus longan Lour.

BACKGROUND: Long non-coding RNAs (lncRNAs) are involved in variable cleavage, transcriptional interference, regulation of DNA methylation and protein modification. However, the regulation of lncRNAs in plant somatic embryos remains unclear. The longan (Dimocarpus longan) somatic embryogenesis (SE) system is a good system for research on longan embryo development. RESULTS: In this study, 7643 lncRNAs obtained during early SE in D. longan were identified by high-throughput sequencing, among which 6005 lncRNAs were expressed. Of the expressed lncRNAs, 4790 were found in all samples and 160 were specifically expressed in embryogenic callus (EC), 154 in incomplete embryogenic compact structures (ICpECs), and 376 in globular embryos (GEs). We annotated the 6005 expressed lncRNAs, and 1404 lncRNAs belonged to 506 noncoding RNA (ncRNA) families and 4682 lncRNAs were predicted to target protein-coding genes. The target genes included 5051 cis-regulated target genes (5712 pairs) and 1605 trans-regulated target genes (3618 pairs). KEGG analysis revealed that most of the differentially expressed target genes (mRNAs) of the lncRNAs were enriched in the "plant-pathogen interaction" and "plant hormone signaling" pathways during early longan SE. Real-time quantitative PCR confirmed that 20 selected lncRNAs showed significant differences in expression and that five lncRNAs were related to auxin response factors. Compared with the FPKM expression trends, 16 lncRNA expression trends were the same in qPCR. In lncRNA-miRNA-mRNA relationship prediction, 40 lncRNAs were predicted to function as eTMs for 15 miRNAs and 7 lncRNAs were identified as potential miRNA precursors. In addition, we verified the lncRNA-miRNA-mRNA regulatory relationships by transient expression of miRNAs (miR172a, miR159a.1 and miR398a). CONCLUSION: Analyses of lncRNAs during early longan SE showed that differentially expressed lncRNAs were involved in expression regulation at each SE stage, and may form a regulatory network with miRNAs and mRNAs. These findings provide new insights into lncRNAs and lay a foundation for future functional analysis of lncRNAs during early longan SE.