New insights into the phylogeny of the TMBIM superfamily across the tree of life: Comparative genomics and synteny networks reveal independent evolution of the BI and LFG families in plants

Ancient Duplication and Lineage-Specific Transposition Determine Evolutionary Trajectory of ERF Subfamily across Angiosperms

AP2/ERF transcription factor family plays an important role in plant development and stress responses. Previous studies have shed light on the evolutionary trajectory of the AP2 and DREB subfamilies. However, knowledge about the evolutionary history of the ERF subfamily in angiosperms still remains limited. In this study, we performed a comprehensive analysis of the ERF subfamily from 107 representative angiosperm species by combining phylogenomic and synteny network approaches. We observed that the expansion of the ERF subfamily was driven not only by whole-genome duplication (WGD) but also by tandem duplication (TD) and transposition duplication events. We also found multiple transposition events in Poaceae, Brassicaceae, Poales, Brassicales, and Commelinids. These events may have had notable impacts on copy number variation and subsequent functional divergence of the ERF subfamily. Moreover, we observed a number of ancient tandem duplications occurred in the ERF subfamily across angiosperms, e.g., in Subgroup IX, IXb originated from ancient tandem duplication events within IXa. These findings together provide novel insights into the evolution of this important transcription factor family.

The origin and evolution of RNase T2 family and gametophytic self-incompatibility system in plants

Ribonuclease (RNase) T2 genes are found widely in both eukaryotes and prokaryotes, and genes from this family have been revealed to have various functions in plants. In particular, S-RNase is known to be the female determinant in the S-RNase-based gametophytic self-incompatibility (GSI) system. However, the origin and evolution of the RNase T2 gene family and GSI system are not well understood. In this study, 785 RNase T2 genes were identified in 81 sequenced plant genomes representing broad-scale diversity and divided into three subgroups (Class I, II, and III) based on phylogenetic and synteny network analysis. Class I was found to be of ancient origin and to emerge in green algae, Class II was shown to originate with the appearance of angiosperms, while Class III was discovered to be eudicot-specific. Each of the three major classes could be further classified into several subclasses of which some subclasses were found to be lineage-specific. Furthermore, duplication, deletion, or inactivation of the S/S-like-locus was revealed to be linked to repeated loss and gain of self-incompatibility in different species from distantly related plant families with GSI. Finally, the origin and evolutionary history of S-locus in Rosaceae species was unraveled with independent loss and gain of S-RNase occurred in different subfamilies of Rosaceae. Our findings provide insights into the origin and evolution of the RNase T2 family and the GSI system in plants.

Applying Synteny Networks (SynNet) to Study Genomic Arrangements of Protein-Coding Genes in Plants

In comparative genomics, the study of synteny can be a powerful method for exploring genome rearrangements, inferring genomic ancestry, defining orthology relationships, determining gene and genome duplications, and inferring gene positional conservation patterns across taxa. In this chapter, we present a step-by-step protocol for microsynteny network (SynNet) analysis, as an alternative to traditional methods of synteny comparison, where nodes in the network represent protein-coding genes and edges represent the pairwise syntenic relationships. The SynNet pipeline consists of six main steps: (1) pairwise genome comparisons between all the genomes being analyzed, (2) detection of inter- and intrasynteny blocks, (3) generation of an entire synteny database (i.e., edgelist), (4) network clustering, (5) phylogenomic profiling of the gene family of interest, and (6) evolutionary inference. The SynNet approach facilitates the rapid analysis and visualization of synteny relationships (from specific genes, specific gene families up to all genes) across a large number of genomes.

The Evolutionary Conserved Transmembrane BAX Inhibitor Motif (TMBIM) Containing Protein Family Members 5 and 6 Are Essential for the Development and Survival of Drosophila melanogaster

The mammalian Transmembrane BAX Inhibitor Motif (TMBIM) protein family consists of six evolutionarily conserved hydrophobic proteins that affect programmed cell death and the regulation of intracellular calcium levels. The bacterial ortholog BsYetJ is a pH-dependent calcium channel. We here identified seven TMBIM family members in Drosophila melanogaster and describe their expression levels in diverse tissues and developmental stages. A phylogenetic analysis revealed that CG30379 represents the ortholog of human TMBIM4 although these two proteins are much less related than TMBIM5 (CG2076 and CG1287/Mics1) and TMBIM6 (CG7188/Bi-1) to their respective orthologs. For TMBIM1-3 the assignment is more dubious because the fly and the human proteins cluster together. We conducted a functional analysis based on expression levels and the availability of RNAi lines. This revealed that the ubiquitous knockdown of CG3798/Nmda1 and CG3814/Lfg had no effect on development while knockdown of CG2076/dTmbim5 resulted in death at the pupa stage and knockdown of CG7188/dTmbim6 in death at the embryonic stage. Ubiquitous knockdown of the second TMBIM5 paralog CG1287/Mics1 ensued in male sterility. Knockdown of dTmbim5 and 6 in muscle and neural tissue also greatly reduced lifespan through different mechanisms. Knockdown of the mitochondrial family member dTmbim5 resulted in reduced ATP production and a pro-apoptotic expression profile while knockdown of the ER protein dTmbim6 increased the ER calcium levels similar to findings in mammalian cells. Our data demonstrate that dTmbim5 and 6 are essential for fly development and survival but affect cell survival through different mechanisms.

Integrated Proteomic and Transcriptomic Analyses Reveal the Roles of Brucella Homolog of BAX Inhibitor 1 in Cell Division and Membrane Homeostasis of Brucella suis S2

BAX inhibitor 1 (BI-1) is an evolutionarily conserved transmembrane protein first identified in a screening process for human proteins that suppress BAX-induced apoptosis in yeast cells. Eukaryotic BI-1 is a cytoprotective protein that suppresses cell death induced by multiple stimuli in eukaryotes. Brucella, the causative agent of brucellosis that threatens public health and animal husbandry, contains a conserved gene that encodes BI-1-like protein. To explore the role of the Brucella homolog of BI-1, BrBI, in Brucella suis S2, we constructed the brbI deletion mutant strain and its complemented strain. brbI deletion altered the membrane properties of Brucella suis S2 and decreased its resistance to acidic pH, H₂O₂, polymyxin B, and lincomycin. Additionally, deleting brbI led to defective growth, cell division, and viability in Brucella suis S2. We then revealed the effect of brbI deletion on the physiological characteristics of Brucella suis S2 via integrated transcriptomic and proteomic analyses. The integrated analysis showed that brbI deletion significantly affected the expression of multiple genes at the mRNA and/or protein levels. Specifically, the affected divisome proteins, FtsB, FtsI, FtsL, and FtsQ, may be the molecular basis of the impaired cell division of the brbI mutant strain, and the extensively affected membrane proteins and transporter-associated proteins were consistent with the phenotype of the membrane properties’ alterations of the brbI mutant strain. In conclusion, our results revealed that BrBI is a bacterial cytoprotective protein involved in membrane homeostasis, cell division, and stress resistance in Brucella suis S2.

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

An aphid RNA transcript migrates systemically within plants and is a virulence factor

Aphids are sap-feeding insects that colonize a broad range of plant species and often cause feeding damage and transmit plant pathogens, including bacteria, viruses, and viroids. These insects feed from the plant vascular tissue, predominantly the phloem. However, it remains largely unknown how aphids, and other sap-feeding insects, establish intimate long-term interactions with plants. To identify aphid virulence factors, we took advantage of the ability of the green peach aphid Myzus persicae to colonize divergent plant species. We found that a M. persicae clone of near-identical females established stable colonies on nine plant species of five representative plant eudicot and monocot families that span the angiosperm phylogeny. Members of the novel aphid gene family Ya are differentially expressed in aphids on the nine plant species and are coregulated and organized as tandem repeats in aphid genomes. Aphids translocate Ya transcripts into plants, and some transcripts migrate to distal leaves within several plant species. RNAi-mediated knockdown of Ya genes reduces M. persicae fecundity, and M. persicae produces more progeny on transgenic plants that heterologously produce one of the systemically migrating Ya transcripts as a long noncoding (lnc) RNA. Taken together, our findings show that beyond a range of pathogens, M. persicae aphids translocate their own transcripts into plants, including a Ya lncRNA that migrates to distal locations within plants, promotes aphid fecundity, and is a member of a previously undescribed host-responsive aphid gene family that operate as virulence factors.

Phylogenomic synteny network analyses reveal ancestral transpositions of auxin response factor genes in plants

BACKGROUND

Auxin response factors (ARFs) have long been a research focus and represent a class of key regulators of plant growth and development. Integrated phylogenomic synteny network analyses were able to provide novel insights into the evolution of the ARF gene family.

RESULTS

Here, more than 3500 ARFs collected from plant genomes and transcriptomes covering major streptophyte lineages were used to reconstruct the broad-scale family phylogeny, where the early origin and diversification of ARF in charophytes was delineated. Based on the family phylogeny, we proposed a unified six-group classification system for angiosperm ARFs. Phylogenomic synteny network analyses revealed the deeply conserved genomic syntenies within each of the six ARF groups and the interlocking syntenic relationships connecting distinct groups. Recurrent duplication events, such as those that occurred in seed plants, angiosperms, core eudicots and grasses contributed to the expansion of ARF genes which facilitated functional diversification. Ancestral transposition activities in important plant families, including crucifers, legumes and grasses, were unveiled by synteny network analyses. Ancestral gene duplications along with transpositions have profound evolutionary significance which may have accelerated the functional diversification process of paralogues.

CONCLUSIONS

The broad-scale family phylogeny in combination with the state-of-art phylogenomic synteny network analyses not only allowed us to infer the evolutionary trajectory of ARF genes across distinct plant lineages, but also facilitated to generate a more robust classification regime for this transcription factor family. Our study provides insights into the evolution of ARFs which will enhance our current understanding of this important transcription factor family.

Deciphering GRINA/Lifeguard1: Nuclear Location, Ca2+ Homeostasis and Vesicle Transport

The Glutamate Receptor Ionotropic NMDA-Associated Protein 1 (GRINA) belongs to the Lifeguard family and is involved in calcium homeostasis, which governs key processes, such as cell survival or the release of neurotransmitters. GRINA is mainly associated with membranes of the endoplasmic reticulum, Golgi, endosome, and the cell surface, but its presence in the nucleus has not been explained yet. Here we dissect, with the help of different software tools, the potential roles of GRINA in the cell and how they may be altered in diseases, such as schizophrenia or celiac disease. We describe for the first time that the cytoplasmic N-terminal half of GRINA (which spans a Proline-rich domain) contains a potential DNA-binding sequence, in addition to cleavage target sites and probable PY-nuclear localization sequences, that may enable it to be released from the rest of the protein and enter the nucleus under suitable conditions, where it could participate in the transcription, alternative splicing, and mRNA export of a subset of genes likely involved in lipid and sterol synthesis, ribosome biogenesis, or cell cycle progression. To support these findings, we include additional evidence based on an exhaustive review of the literature and our preliminary data of the protein–protein interaction network of GRINA.

Ion and pH Sensitivity of a TMBIM Ca2+ Channel

The anti-apoptotic transmembrane Bax inhibitor motif (TMBIM) containing protein family regulates Ca²⁺ homeostasis, cell death, and the progression of diseases including cancers. The recent crystal structures of the TMBIM homolog BsYetJ reveal a conserved Asp171-Asp195 dyad that is proposed in regulating a pH-dependent Ca²⁺ translocation. Here we show that BsYetJ mediates Ca²⁺ fluxes in permeabilized mammalian cells, and its interaction with Ca²⁺ is sensitive to protons and other cations. We report crystal structures of BsYetJ in additional states, revealing the flexibility of the dyad in a closed state and a pore-opening mechanism. Functional studies show that the dyad is responsible for both Ca²⁺ affinity and pH dependence. Computational simulations suggest that protonation of Asp171 weakens its interaction with Arg60, leading to an open state. Our integrated analysis provides insights into the regulation of the BsYetJ Ca²⁺ channel that may inform understanding of human TMBIM proteins regarding their roles in cell death and diseases.