Functional Divergence in Orthologous Transcription Factors: Insights from AtCBF2/3/1 and OsDREB1C.

Despite traditional beliefs of orthologous genes maintaining similar functions across species, growing evidence points to their potential for functional divergence. C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) are critical in cold acclimation, with their overexpression enhancing stress tolerance but often constraining plant growth. In contrast, a recent study unveiled a distinctive role of rice OsDREB1C in elevating nitrogen use efficiency (NUE), photosynthesis, and grain yield, implying functional divergence within the CBF/DREB1 orthologs across species. Here, we delve into divergent molecular mechanisms of OsDREB1C and AtCBF2/3/1 by exploring their evolutionary trajectories across rice and Arabidopsis genomes, regulatomes, and transcriptomes. Evolutionary scrutiny shows discrete clades for OsDREB1C and AtCBF2/3/1, with the Poaceae-specific DREB1C clade mediated by a transposon event. Genome-wide binding profiles highlight OsDREB1C's preference for GCCGAC compared to AtCBF2/3/1's preference for A/GCCGAC, a distinction determined by R12 in the OsDREB1C AP2/ERF domain. Cross-species multiomic analyses reveal shared gene orthogroups (OGs) and underscore numerous specific OGs uniquely bound and regulated by OsDREB1C, implicated in NUE, photosynthesis, and early flowering, or by AtCBF2/3/1, engaged in hormone and stress responses. This divergence arises from gene gains/losses (∼16.7% to 25.6%) and expression reprogramming (∼62.3% to 66.2%) of OsDREB1C- and AtCBF2/3/1-regulated OGs during the extensive evolution following the rice-Arabidopsis split. Our findings illustrate the regulatory evolution of OsDREB1C and AtCBF2/3/1 at a genomic scale, providing insights on the functional divergence of orthologous transcription factors following gene duplications across species.

Ectopic expression of BOTRYTIS SUSCEPTIBLE1 reveals its function as a positive regulator of wound-induced cell death and plant susceptibility to Botrytis.

Programmed cell death (PCD) is integral to plant life and required for stress responses, immunity, and development. Our understanding of the regulation of PCD is incomplete, especially concerning regulators involved in multiple divergent processes. The botrytis-susceptible (bos1) mutant of Arabidopsis is highly susceptible to fungal infection by Botrytis cinerea (Botrytis). BOS1 (also known as MYB108) regulates cell death propagation during plant responses to wounding. The bos1-1 allele contains a T-DNA insertion in the 5'-untranslated region upstream of the start codon. This insertion results in elevated expression of BOS1/MYB108. We used clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated nuclease 9 (Cas9) system (CRISPR/Cas9) to create new bos1 alleles with disrupted exons, and found that these lines lacked the typical bos1-1 wounding and Botrytis phenotypes. They did exhibit reduced fertility, as was previously observed in other bos1 alleles. Resequencing of the bos1-1 genome confirmed the presence of a mannopine synthase (MAS) promoter at the T-DNA left border. Expression of the BOS1 gene under control of the MAS promoter in wild-type plants conferred the characteristic phenotypes of bos1-1: Botrytis sensitivity and response to wounding. Multiple overexpression lines demonstrated that BOS1 was involved in regulation of cell death propagation in a dosage-dependent manner. Our data indicate that bos1-1 is a gain-of-function mutant and that BOS1 function in regulation of fertility and Botrytis response can both be understood as misregulated cell death.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization.

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

Role of lncRNAs in cis- and trans-regulatory responses to salt in Populus trichocarpa.

Long non-coding RNAs (lncRNAs) are emerging as versatile regulators in diverse biological processes. However, little is known about their cis- and trans-regulatory contributions in gene expression under salt stress. Using 27 RNA-seq data sets from Populus trichocarpa leaves, stems and roots, we identified 2988 high-confidence lncRNAs, including 1183 salt-induced differentially expressed lncRNAs. Among them, 301 lncRNAs have potential for positively affecting their neighboring genes, predominantly in a cis-regulatory manner rather than by co-transcription. Additionally, a co-expression network identified six striking salt-associated modules with a total of 5639 genes, including 426 lncRNAs, and in these lncRNA sequences, the DNA/RNA binding motifs are enriched. This suggests that lncRNAs might contribute to distant gene expression of the salt-associated modules in a trans-regulatory manner. Moreover, we found 30 lncRNAs that have potential to simultaneously cis- and trans-regulate salt-responsive homologous genes, and Ptlinc-NAC72, significantly induced under long-term salt stress, was selected for validating its regulation of the expression and functional roles of the homologs PtNAC72.A and PtNAC72.B (PtNAC72.A/B). The transient transformation of Ptlinc-NAC72 and a dual-luciferase assay of Ptlinc-NAC72 and PtNAC72.A/B promoters confirmed that Ptlinc-NAC72 can directly upregulate PtNAC72.A/B expression, and a presence/absence assay was further conducted to show that the regulation is probably mediated by Ptlinc-NAC72 recognizing the tandem elements (GAAAAA) in the PtNAC72.A/B 5' untranslated region (5'-UTR). Finally, the overexpression of Ptlinc-NAC72 produces a hypersensitive phenotype under salt stress. Altogether, our results shed light on the cis- and trans-regulation of gene expression by lncRNAs in Populus and provides an example of long-term salt-induced Ptlinc-NAC72 that could be used to mitigate growth costs by conferring plant resilience to salt stress.

Integrative multi-omics analysis of three early diverged rosid species reveals an ancient hierarchical cold-responsive regulatory network.

Elucidating regulators, including transcription factors (TFs) and RNA-binding proteins (RBPs), underlying gene transcriptional and post-transcriptional co-regulatory network is key to understand plant cold responses. Previous studies were mainly conducted on single species, and whether the regulators are conserved across different species remains elusive. Here, we selected three species that diverged at the early evolution of rosids (~99-113 million years ago), performed cold-responsive phylotranscriptome experiments, and integrated chromatin immunoprecipitation- and DNA affinity purification-sequencing (ChIP/DAP-seq) analysis to explore cold-responsive regulators and their regulatory networks. First, we detected over 10,000 cold-induced differentially expressed genes (DEGs) and alternative splicing genes (DASGs) in each species. Among the DEGs, a set of TFs and RBPs were conserved in rosid cold response. Compared to TFs, RBPs displayed a delayed cold-responsive pattern, implying a hierarchical regulation of DEGs and DASGs. By integrating DEGs and DASGs, we identified 259 overlapping DE-DASG orthogroups (closely-related homologs) that were cold-regulated at both transcriptional and post-transcriptional levels in all three studied species. Notably, pathway analysis on each of the DEGs, DASGs, and DE-DASGs in the three species showed a common enrichment connected to the circadian rhythm. Evidently, 226 cold-responsive genes were directly targeted by at least two circadian rhythm components (CCA1, LHY, RV4, RVE7, and RVE8). Finally, we revealed an ancient hierarchy of cold-responsive regulatory networks at transcriptional and post-transcriptional levels launched by circadian components in rosids. Altogether, this study sheds light on conserved regulators underlying cold-responsive regulatory networks across rosid species, despite a long evolutionary history after their divergence.

RNA-Binding Protein MAC5A Is Required for Gibberellin-Regulated Stamen Development.

The development of floral organs is coordinated by an elaborate network of homeotic genes, and gibberellin (GA) signaling is involved in floral organ development; however, the underlying molecular mechanisms remain elusive. In the present study, we found that MOS4-ASSOCIATED COMPLEX 5A (MAC5A), which is a protein containing an RNA-binding motif, was involved in the development of sepals, petals, and stamens; either the loss or gain of MAC5A function resulted in stamen malformation and a reduced seed set. The exogenous application of GA considerably exacerbated the defects in mac5a null mutants, including fewer stamens and male sterility. MAC5A was predominantly expressed in pollen grains and stamens, and overexpression of MAC5A affected the expression of homeotic genes such as APETALA1 (AP1), AP2, and AGAMOUS (AG). MAC5A may interact with RABBIT EARS (RBE), a repressor of AG expression in Arabidopsis flowers. The petal defect in rbe null mutants was at least partly rescued in mac5a rbe double mutants. These findings suggest that MAC5A is a novel factor that is required for the normal development of stamens and depends on the GA signaling pathway.

Innovations and stepwise evolution of CBFs/DREB1s and their regulatory networks in angiosperms.

The C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) have been identified as major regulators of cold acclimation in many angiosperm plants. However, their origin and evolutionary process associated to cold responsiveness are still lacking. By integrating multi-omics data of genomes, transcriptomes, and CBFs/DREB1s genome-wide binding profiles, we unveil the origin and evolution of CBFs/DREB1s and their regulatory network. Gene collinearity and phylogeny analyses show that CBF/DREB1 is an innovation evolved from tandem duplication-derived DREB III gene. A subsequent event of ε-whole genome duplication led to two CBF/DREB1 archetypes (Clades I and II) in ancient angiosperms. In contrast to cold-insensitivity of Clade I and their parent DREB III genes, Clade II evolved a further innovation in cold-sensitive response and was stepwise expanded in eudicots and monocots by independent duplications. In geological time, the duplication events were mainly enriched around the Cretaceous-Paleogene (K-Pg) boundary and/or in the Late Cenozoic Ice Age, when the global average temperature significantly decreased. Consequently, the duplicated CBF/DREB1 genes contributed to the rewiring of CBFs/DREB1s-regulatory network for cold tolerance. Altogether, our results highlight an origin and convergent evolution of CBFs/DREB1s and their regulatory network probably for angiosperms adaptation to global cooling.

The SR Splicing Factors: Providing Perspectives on Their Evolution, Expression, Alternative Splicing, and Function in Populus trichocarpa.

Serine/arginine-rich (SR) proteins are important splicing factors in plant development and abiotic/hormone-related stresses. However, evidence that SR proteins contribute to the process in woody plants has been lacking. Using phylogenetics, gene synteny, transgenic experiments, and RNA-seq analysis, we identified 24 PtSR genes and explored their evolution, expression, and function in Popolus trichocarpa. The PtSR genes were divided into six subfamilies, generated by at least two events of genome triplication and duplication. Notably, they were constitutively expressed in roots, stems, and leaves, demonstrating their fundamental role in P. trichocarpa. Additionally, most PtSR genes (~83%) responded to at least one stress (cold, drought, salt, SA, MeJA, or ABA), and, especially, cold stress induced a dramatic perturbation in the expression and/or alternative splicing (AS) of 18 PtSR genes (~75%). Evidentially, the overexpression of PtSCL30 in Arabidopsis decreased freezing tolerance, which probably resulted from AS changes of the genes (e.g., ICE2 and COR15A) critical for cold tolerance. Moreover, the transgenic plants were salt-hypersensitive at the germination stage. These indicate that PtSCL30 may act as a negative regulator under cold and salt stress. Altogether, this study sheds light on the evolution, expression, and AS of PtSR genes, and the functional mechanisms of PtSCL30 in woody plants.

TPST is involved in fructose regulation of primary root growth in Arabidopsis thaliana.

KEY MESSAGE: TPST is involved in fructose signaling to regulate the root development and expression of genes in biological processes including auxin biosynthesis and accumulation in Arabidopsis. Sulfonation of proteins by tyrosine protein sulfotransferases (TPST) has been implicated in many important biological processes in eukaryotic organisms. Arabidopsis possesses a single TPST gene and its role in auxin homeostasis and root development has been reported. Here we show that the Arabidopsis tpst mutants are hypersensitive to fructose. In contrast to sucrose and glucose, fructose represses primary root growth of various ecotypes of Arabidopsis at low concentrations. RNA-seq analysis identified 636 differentially expressed genes (DEGs) in Col-0 seedlings in response to fructose verses glucose. GO and KEGG analyses of the DEGs revealed that fructose down-regulates genes involved in photosynthesis, glucosinolate biosynthesis and IAA biosynthesis, but up-regulates genes involved in the degradation of branched amino acids, sucrose starvation response, and dark response. The fructose responsive DEGs in the tpst mutant largely overlapped with that in Col-0, and most DEGs in tpst displayed larger changes than in Col-0. Interestingly, the fructose up-regulated DEGs includes genes encoding two AtTPST substrate proteins, Phytosulfokine 2 (PSK2) and Root Meristem Growth Factor 7 (RGF7). Synthesized peptides of PSK-α and RGF7 could restore the fructose hypersensitivity of tpst mutant plants. Furthermore, auxin distribution and accumulation at the root tip were affected by fructose and the tpst mutation. Our findings suggest that fructose serves as a signal to regulate the expression of genes involved in various biological processes including auxin biosynthesis and accumulation, and that modulation of auxin accumulation and distribution in roots by fructose might be partly mediated by the TPST substrate genes PSK-α and RGF7.

CASH: a constructing comprehensive splice site method for detecting alternative splicing events.

RNA-sequencing (RNA-seq) can generate millions of reads to provide clues for analyzing novel or abnormal alternative splicing (AS) events in cells. However, current methods for exploring AS events are still far from being satisfactory. Here, we present Comprehensive AS Hunting (CASH), which constructs comprehensive splice sites including known and novel AS sites in cells, and identifies differentially AS events between cells. We illuminated the versatility of CASH on RNA-seq data from a wide range of species and also on simulated in silico data, validated the advantages of CASH over other AS predictors and exhibited novel differentially AS events. Moreover, we used CASH to identify SRSF10-regulated AS events and investigated the evolution of SRSF10-regulated splicing. The results showed that SRSF10-regulated splicing events are highly evolvable from chickens, mice to humans. However, SRSF10-regulated splicing model was observed to be immutable, in which SRSF10 binding to cassette exon promotes exon inclusion while binding to downstream exon induces exon skipping. Altogether, CASH can significantly improve the detection of AS events and facilitate the study of AS regulation and function in cells; the SRSF10 data first demonstrate a flexibility of SRSF10 with their regulated splicing events but an immutability of SRSF10-regulated splicing model to produce opposite AS outcomes in vertebrates.

Transcriptome profiling of Puccinellia tenuiflora during seed germination under a long-term saline-alkali stress.

BACKGROUND: Puccinellia tenuiflora is the most saline-alkali tolerant plant in the Songnen Plain, one of the three largest soda saline-alkali lands worldwide. Here, we investigated the physicochemical properties of saline-alkali soils from the Songnen Plain and sequenced the transcriptomes of germinated P. tenuiflora seedlings under long-term treatment (from seed soaking) with saline-alkali soil extracts. RESULTS: We found that the soils from Songnen Plain were reasonably rich in salts and alkali; moreover, the soils were severely deficient in nitrogen [N], phosphorus [P], potassium [K] and several other mineral elements. This finding demonstrated that P. tenuiflora can survive from not only high saline-alkali stress but also a lack of essential mineral elements. To explore the saline-alkali tolerance mechanism, transcriptional analyses of P. tenuiflora plants treated with water extracts from the saline-alkali soils was performed. Interestingly, unigenes involved in the uptake of N, P, K and the micronutrients were found to be significantly upregulated, which indicated the existence of an efficient nutrition-uptake system in P. tenuiflora. Compared with P. tenuiflora, the rice Oryza sativa was hypersensitive to saline-alkali stress. The results obtained using a noninvasive microtest techniques confirmed that the uptake of NO3- and NH4+ and the regulatory flux of Na+ and H+ were significantly higher in the roots of P. tenuiflora than in those of O. sativa. In the corresponding physiological experiments, the application of additional nutrition elements significantly eliminated the sensitive symptoms of rice to saline-alkali soil extracts. CONCLUSIONS: Our results imply that the survival of P. tenuiflora in saline-alkali soils is due to a combination of at least two regulatory mechanisms and the high nutrient uptake capacity of P. tenuiflora plays a pivotal role in its adaptation to those stress. Taken together, our results highlight the role of nutrition uptake in saline-alkali stress tolerance in plants.

U2-related proteins CHERP and SR140 contribute to colorectal tumorigenesis via alternative splicing regulation.

Dysregulation of calcium homeostasis endoplasmic reticulum protein (CHERP) has been implicated in several cancers, but it remains elusive how CHERP contributes to cancer cell proliferation and cancer development. Here, we observed that CHERP and its binding partner SR140 are significantly upregulated in human clinical colorectal cancer tissues (CRC). CHERP and SR140 could form a protein complex to stabilize each other. Knockdown of CHERP or SR140 triggers double-stranded DNA breaks and cell death. Furthermore, UPF3A, the RNA surveillance factor, was identified as a splicing target of CHERP and SR140, which bind specifically to the regulated exon4 and modulate UPF3A splicing. UPF3A knockdown recapitulates CHERP/SR140 depletion both in vitro and in mice. Importantly, overexpression of UPF3A significantly rescues proliferation defect of CHERP/SR140-depleted cells. These results confirmed that the effect of CHERP/SR140 in promoting tumorigenesis was partially mediated by UPF3A. Extending these results, upregulation of CHERP/SR140 observed in CRC remarkably parallels increased inclusion of UPF3A exon4. Together, our study clarifies how CHERP/SR140 exert an oncogenic role in CRC development partially through regulating expression of UPF3A variants.

Cell death regulation but not abscisic acid signaling is required for enhanced immunity to Botrytis in Arabidopsis cuticle-permeable mutants.

Prevailing evidence indicates that abscisic acid (ABA) negatively influences immunity to the fungal pathogen Botrytis cinerea in most but not all cases. ABA is required for cuticle biosynthesis, and cuticle permeability enhances immunity to Botrytis via unknown mechanisms. This complex web of responses obscures the role of ABA in Botrytis immunity. Here, we addressed the relationships between ABA sensitivity, cuticle permeability, and Botrytis immunity in the Arabidopsis thaliana ABA-hypersensitive mutants protein phosphatase2c quadruple mutant (pp2c-q) and enhanced response to aba1 (era1-2). Neither pp2c-q nor era1-2 exhibited phenotypes predicted by the known roles of ABA; conversely, era1-2 had a permeable cuticle and was Botrytis resistant. We employed RNA-seq analysis in cuticle-permeable mutants of differing ABA sensitivities and identified a core set of constitutively activated genes involved in Botrytis immunity and susceptibility to biotrophs, independent of ABA signaling. Furthermore, botrytis susceptible1 (bos1), a mutant with deregulated cell death and enhanced ABA sensitivity, suppressed the Botrytis immunity of cuticle permeable mutants, and this effect was linearly correlated with the extent of spread of wound-induced cell death in bos1. Overall, our data demonstrate that Botrytis immunity conferred by cuticle permeability can be genetically uncoupled from PP2C-regulated ABA sensitivity, but requires negative regulation of a parallel ABA-dependent cell-death pathway.

Sulforaphane, a Natural Isothiocyanate Compound, Improves Cardiac Function and Remodeling by Inhibiting Oxidative Stress and Inflammation in a Rabbit Model of Chronic Heart Failure.

BACKGROUND The aim of this study was to investigate the effects of sulforaphane (SFN), a natural isothiocyanate compound, in a rabbit ascending aortic cerclage model of chronic heart failure (CHF). MATERIAL AND METHODS Thirty New Zealand White rabbits were divided into the sham operation group (n=10), the CHF group (n=10), and the CHF + SFN group (n=10) treated with subcutaneous SFN (0.5 mg/kg) for five days per week for 12 weeks. After 12 weeks, echocardiography and biometric analysis were performed, followed by the examination of the rabbit hearts. Enzyme-linked immunosorbent assay (ELISA) and Western blot were used to detect levels of inflammatory cytokines, superoxide dismutase (SOD), and malondialdehyde (MDA). RESULTS In the CHF group, compared with the sham operation group, there was an increase in the heart weight to body weight ratio (HW/BW), the left ventricular weight to body weight ratio (LVW/BW), the left ventricular end diastolic diameter (LVEDD), the left ventricular end systolic diameter (LVESD), plasma brain natriuretic peptide (BNP) and atrial natriuretic peptide (ANP) levels, the cardiac collagen volume fraction (CVF), apoptotic index, expression levels of collagen I, collagen III, tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6) and malondialdehyde (MDA) in the myocardial tissue, and a decrease in the left ventricular shortening fraction (LVFS) and left ventricular ejection fraction (LVEF), and cardiac superoxide dismutase (SOD) activity. These changes were corrected in the SFN-treated group. CONCLUSIONS In a rabbit model of CHF, treatment with SFN improved cardiac function and remodeling by inhibiting oxidative stress and inflammation.

The inhibition of protein translation mediated by AtGCN1 is essential for cold tolerance in Arabidopsis thaliana.

In yeast, the interaction of General Control Non-derepressible 1 (GCN1) with GCN2 enables GCN2 to phosphorylate eIF2α (the alpha subunit of eukaryotic translation initiation factor 2) under a variety of stresses. Here, we cloned AtGCN1, an Arabidopsis homologue of GCN1. We show that AtGCN1 directly interacts with GCN2 and is essential for the phosphorylation of eIF2α under salicylic acid (SA), ultraviolet (UV), cold stress and amino acid deprivation conditions. Two mutant alleles, atgcn1-1 and atgcn1-2, which are defective in the phosphorylation of eIF2α, showed increased sensitivity to cold stress, compared with the wild type. Ribosome-bound RNA profiles showed that the translational state of mRNA was higher in atgcn1-1 than in the wild type. Our result also showed that cold treatment reduced the tendency of the tor mutant seedlings to produce purple hypocotyls. In addition, the kinase activity of TOR was transiently inhibited when plants were exposed to cold stress, suggesting that the inhibition of TOR is another pathway important for plants to respond to cold stress. In conclusion, our results indicate that the AtGCN1-mediated phosphorylation of eIF2α, which is required for inhibiting the initiation of protein translation, is essential for cold tolerance in Arabidopsis.

Transcriptome analysis of alternative splicing events regulated by SRSF10 reveals position-dependent splicing modulation.

Splicing factor SRSF10 is known to function as a sequence-specific splicing activator. Here, we used RNA-seq coupled with bioinformatics analysis to identify the extensive splicing network regulated by SRSF10 in chicken cells. We found that SRSF10 promoted both exon inclusion and exclusion. Motif analysis revealed that SRSF10 binding to cassette exons was associated with exon inclusion, whereas the binding of SRSF10 within downstream constitutive exons was associated with exon exclusion. This positional effect was further demonstrated by the mutagenesis of potential SRSF10 binding motifs in two minigene constructs. Functionally, many of SRSF10-verified alternative exons are linked to pathways of stress and apoptosis. Consistent with this observation, cells depleted of SRSF10 expression were far more susceptible to endoplasmic reticulum stress-induced apoptosis than control cells. Importantly, reconstituted SRSF10 in knockout cells recovered wild-type splicing patterns and considerably rescued the stress-related defects. Together, our results provide mechanistic insight into SRSF10-regulated alternative splicing events in vivo and demonstrate that SRSF10 plays a crucial role in cell survival under stress conditions.

A new computational strategy for predicting essential genes.

BACKGROUND: Determination of the minimum gene set for cellular life is one of the central goals in biology. Genome-wide essential gene identification has progressed rapidly in certain bacterial species; however, it remains difficult to achieve in most eukaryotic species. Several computational models have recently been developed to integrate gene features and used as alternatives to transfer gene essentiality annotations between organisms. RESULTS: We first collected features that were widely used by previous predictive models and assessed the relationships between gene features and gene essentiality using a stepwise regression model. We found two issues that could significantly reduce model accuracy: (i) the effect of multicollinearity among gene features and (ii) the diverse and even contrasting correlations between gene features and gene essentiality existing within and among different species. To address these issues, we developed a novel model called feature-based weighted Naïve Bayes model (FWM), which is based on Naïve Bayes classifiers, logistic regression, and genetic algorithm. The proposed model assesses features and filters out the effects of multicollinearity and diversity. The performance of FWM was compared with other popular models, such as support vector machine, Naïve Bayes model, and logistic regression model, by applying FWM to reciprocally predict essential genes among and within 21 species. Our results showed that FWM significantly improves the accuracy and robustness of essential gene prediction. CONCLUSIONS: FWM can remarkably improve the accuracy of essential gene prediction and may be used as an alternative method for other classification work. This method can contribute substantially to the knowledge of the minimum gene sets required for living organisms and the discovery of new drug targets.

Differences in alternative splicing and their potential underlying factors between animals and plants.

BACKGROUND: Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed. AIM OF REVIEW: This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants. KEY SCIENTIFIC CONCEPTS OF REVIEW: We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

A phylotranscriptomic dataset of angiosperm species under cold stress.

Angiosperms are one of the most diverse and abundant plant groups that are widely distributed on Earth, from tropical to temperate and polar zones. The wide distribution of angiosperms may be attributed to the evolution of sophisticated mechanisms of environmental adaptability, including cold tolerance. Since the development of high-throughput sequencing, transcriptome has been widely utilized to gain insights into the molecular mechanisms of plants in response to cold stress. However, previous studies generally focused on single or two species, and comparative transcriptome analyses for multispecies responding to cold stress were limited. In this study, we selected 11 representative angiosperm species, performed phylotranscriptome experiments at four time points before and after cold stress, and presented a profile of cold-induced transcriptome changes in angiosperms. Our multispecies cold-responsive RNA-seq datasets provide valuable references for exploring conserved and evolutionary mechanisms of angiosperms in adaptation to cold stress.

Orthogroup and phylotranscriptomic analyses identify transcription factors involved in the plant cold response: A case study of Arabidopsis BBX29.

C-repeat binding factors (CBFs) are well-known transcription factors (TFs) that regulate plant cold acclimation. RNA sequencing (RNA-seq) data from diverse plant species provide opportunities to identify other TFs involved in the cold response. However, this task is challenging because gene gain and loss has led to an intertwined community of co-orthologs and in-paralogs between and within species. Using orthogroup (closely related homologs) analysis, we identified 10,549 orthogroups in five representative eudicots. A phylotranscriptomic analysis of cold-treated seedlings from eudicots identified 35 high-confidence conserved cold-responsive transcription factor orthogroups (CoCoFos). These 35 CoCoFos included the well-known cold-responsive regulators CBFs, HSFC1, ZAT6/10, and CZF1 among others. We used Arabidopsis BBX29 for experimental validation. Expression and genetic analyses showed that cold-induction of BBX29 is CBF- and abscisic acid-independent, and BBX29 is a negative regulator of cold tolerance. Integrative RNA-seq and Cleavage Under Targets and Tagmentation followed by sequencing analyses revealed that BBX29 represses a set of cold-induced TFs (ZAT12, PRR9, RVE1, MYB96, etc.). Altogether, our analysis yielded a library of eudicot CoCoFos and demonstrated that BBX29 is a negative regulator of cold tolerance in Arabidopsis.