Bioinformatic Analysis of the BTB Gene Family in Paulownia fortunei and Functional Characterization in Response to Abiotic and Biotic Stresses.

To learn about the gene structure, phylogenetic evolution, and function under biotic and abiotic stresses of BTB (Bric-a-Brac/Tramtrack/Broad Complex) genes in Paulownia fortunei, a whole-genome sequence evaluation was carried out, and a total of 62 PfBTB genes were identified. The phylogenetic analysis showed that PfBTB proteins are divided into eight groups, and these proteins are highly conserved. PfBTB genes were unevenly distributed on 17 chromosomes. The colinearity analysis found that fragment replication and tandem replication are the main modes of gene amplification in the PfBTB family. The analysis of cis-acting elements suggests that PfBTB genes may be involved in a variety of biological processes. The transcriptomic analysis results showed that PfBTB3/12/14/16/19/36/44 responded to Paulownia witches' broom (PaWB), while PfBTB1/4/17/43 responded to drought stress, and the RT-qPCR results further support the reliability of transcriptome data. In addition, the association analysis between miRNA and transcriptome revealed a 91-pair targeting relationship between miRNAs and PfBTBs. In conclusion, the BTB genes in Paulownia are systematically identified in this research. This work provides useful knowledge to more fully appreciate the potential functions of these genes and their possible roles in the occurrence of PaWB and in response to stress.

N6-methyladenosine modification changes during the recovery processes for Paulownia witches' broom disease under the methyl methanesulfonate treatment.

Phytoplasmas induce diseases in more than 1000 plant species and cause substantial ecological damage and economic losses, but the specific pathogenesis of phytoplasma has not yet been clarified. N 6-methyladenosine (m6A) is the most common internal modification of the eukaryotic Messenger RNA (mRNA). As one of the species susceptible to phytoplasma infection, the pathogenesis and mechanism of Paulownia has been extensively studied by scholars, but the m6A transcriptome map of Paulownia fortunei (P. fortunei) has not been reported. Therefore, this study aimed to explore the effect of phytoplasma infection on m6A modification of P. fortunei and obtained the whole transcriptome m6A map in P. fortunei by m6A-seq. The m6A-seq results of Paulownia witches' broom (PaWB) disease and healthy samples indicate that PaWB infection increased the degree of m6A modification of P. fortunei. The correlation analysis between the RNA-seq and m6A-seq data detected that a total of 315 differentially methylated genes were predicted to be significantly differentially expressed at the transcriptome level. Moreover, the functions of PaWB-related genes were predicted by functional enrichment analysis, and two genes related to maintenance of the basic mechanism of stem cells in shoot apical meristem were discovered. One of the genes encodes the receptor protein kinase CLV2 (Paulownia_LG2G000076), and the other gene encodes the homeobox transcription factor STM (Paulownia_LG15G000976). In addition, genes F-box (Paulownia_LG17G000760) and MSH5 (Paulownia_LG8G001160) had exon skipping and mutually exclusive exon types of alternative splicing in PaWB-infected seedling treated with methyl methanesulfonate, and m6A modification was found in m6A-seq results. Moreover, Reverse Transcription-Polymerase Chain Reaction (RT-PCR) verified that the alternative splicing of these two genes was associated with m6A modification. This comprehensive map provides a solid foundation for revealing the potential function of the mRNA m6A modification in the process of PaWB. In future studies, we plan to verify genes directly related to PaWB and methylation-related enzymes in Paulownia to elucidate the pathogenic mechanism of PaWB caused by phytoplasma invasion.

Genomic insights into the fast growth of paulownias and the formation of Paulownia witches' broom.

Paulownias are among the fastest growing trees in the world, but they often suffer tremendous loss of wood production due to infection by Paulownia witches' broom (PaWB) phytoplasmas. In this study, we have sequenced and assembled a high-quality nuclear genome of Paulownia fortunei, a commonly cultivated paulownia species. The assembled genome of P. fortunei is 511.6 Mb in size, with 93.2% of its sequences anchored to 20 pseudo-chromosomes, and it contains 31 985 protein-coding genes. Phylogenomic analyses show that the family Paulowniaceae is sister to a clade composed of Phrymaceae and Orobanchaceae. Higher photosynthetic efficiency is achieved by integrating C3 photosynthesis and the crassulacean acid metabolism pathway, which may contribute to the extremely fast growth habit of paulownia trees. Comparative transcriptome analyses reveal modules related to cambial growth and development, photosynthesis, and defense responses. Additional genome sequencing of PaWB phytoplasma, combined with functional analyses, indicates that the effector PaWB-SAP54 interacts directly with Paulownia PfSPLa, which in turn causes the degradation of PfSPLa by the ubiquitin-mediated pathway and leads to the formation of witches' broom. Taken together, these results provide significant insights into the biology of paulownias and the regulatory mechanism for the formation of PaWB.