The Rho GTPase Family Genes in Bivalvia Genomes: Sequence, Evolution and Expression Analysis

The complex, dynamic SpliceOme of the small GTPase transcripts altered by technique, sex, genetics, tissue specificity, and RNA base editing

The small GTPase family is well-studied in cancer and cellular physiology. With 162 annotated human genes, the family has a broad expression throughout cells of the body. Members of the family have multiple exons that require splicing. Yet, the role of splicing within the family has been underexplored. We have studied the splicing dynamics of small GTPases throughout 41,671 samples by integrating Nanopore and Illumina sequencing techniques. Within this work, we have made several discoveries. 1). Using the GTEx long read data of 92 samples, each small GTPase gene averages two transcripts, with 83 genes (51%) expressing two or more isoforms. 2). Cross-tissue analysis of GTEx from 17,382 samples shows 41 genes (25%) expressing two or more protein-coding isoforms. These include protein-changing transcripts in genes such as RHOA, RAB37, RAB40C, RAB4B, RAB5C, RHOC, RAB1A, RAN, RHEB, RAC1, and KRAS. 3). The isolation and library technique of the RNAseq influences the abundance of non-sense-mediated decay and retained intron transcripts of small GTPases, which are observed more often in genes than appreciated. 4). Analysis of 16,243 samples of "Blood PAXgene" identified seven genes (3.7%; RHOA, RAB40C, RAB4B, RAB37, RAB5B, RAB5C, RHOC) with two or more transcripts expressed as the major isoform (75% of the total gene), suggesting a role of genetics in altering splicing. 5). Rare (ARL6, RAB23, ARL13B, HRAS, NRAS) and common variants (GEM, RHOC, MRAS, RAB5B, RERG, ARL16) can influence splicing and have an impact on phenotypes and diseases. 6). Multiple genes (RAB9A, RAP2C, ARL4A, RAB3A, RAB26, RAB3C, RASL10A, RAB40B, and HRAS) have sex differences in transcript expression. 7). Several exons are included or excluded for small GTPase genes (RASEF, KRAS, RAC1, RHEB, ARL4A, RHOA, RAB30, RHOBTB1, ARL16, RAP1A) in one or more forms of cancer. 8). Ten transcripts are altered in hypoxia (SAR1B, IFT27, ARL14, RAB11A, RAB10, RAB38, RAN, RIT1, RAB9A) with RHOA identified to have a transient 3'UTR RNA base editing at a conserved site found in all of its transcripts. Overall, we show a remarkable and dynamic role of splicing within the small GTPase family that requires future explorations.

Miro proteins and their role in mitochondrial transfer in cancer and beyond

Mitochondria are organelles essential for tumor cell proliferation and metastasis. Although their main cellular function, generation of energy in the form of ATP is dispensable for cancer cells, their capability to drive their adaptation to stress originating from tumor microenvironment makes them a plausible therapeutic target. Recent research has revealed that cancer cells with damaged oxidative phosphorylation import healthy (functional) mitochondria from surrounding stromal cells to drive pyrimidine synthesis and cell proliferation. Furthermore, it has been shown that energetically competent mitochondria are fundamental for tumor cell migration, invasion and metastasis. The spatial positioning and transport of mitochondria involves Miro proteins from a subfamily of small GTPases, localized in outer mitochondrial membrane. Miro proteins are involved in the structure of the MICOS complex, connecting outer and inner-mitochondrial membrane; in mitochondria-ER communication; Ca²⁺ metabolism; and in the recycling of damaged organelles via mitophagy. The most important role of Miro is regulation of mitochondrial movement and distribution within (and between) cells, acting as an adaptor linking organelles to cytoskeleton-associated motor proteins. In this review, we discuss the function of Miro proteins in various modes of intercellular mitochondrial transfer, emphasizing the structure and dynamics of tunneling nanotubes, the most common transfer modality. We summarize the evidence for and propose possible roles of Miro proteins in nanotube-mediated transfer as well as in cancer cell migration and metastasis, both processes being tightly connected to cytoskeleton-driven mitochondrial movement and positioning.

Molecular Mechanisms Underlying Vibrio Tolerance in Ruditapes philippinarum Revealed by Comparative Transcriptome Profiling

The clam Ruditapes philippinarum is an important species in the marine aquaculture industry in China. However, in recent years, the aquaculture of R. philippinarum has been negatively impacted by various bacterial pathogens. In this study, the transcriptome libraries of R. philippinarum showing different levels of resistance to challenge with Vibrio anguillarum were constructed and RNA-seq was performed using the Illumina sequencing platform. Host immune factors were identified that responded to V. anguillarum infection, including C-type lectin domain, glutathione S-transferase 9, lysozyme, methyltransferase FkbM domain, heat shock 70 kDa protein, Ras-like GTP-binding protein RHO, C1q, F-box and BTB/POZ domain protein zf-C2H2. Ten genes were selected and verified by RT-qPCR, and nine of the gene expression results were consistent with those of RNA-seq. The lectin gene in the phagosome pathway was expressed at a significantly higher level after V. anguillarum infection, which might indicate the role of lectin in the immune response to V. anguillarum. Comparing the results from R. philippinarum resistant and nonresistant to V. anguillarum increases our understanding of the resistant genes and key pathways related to Vibrio challenge in this species. The results obtained here provide a reference for future immunological research focusing on the response of R. philippinarum to V. anguillarum infection.

Genome-wide analysis of acute low salinity tolerance in the eastern oyster Crassostrea virginica and potential of genomic selection for trait improvement

As the global demand for seafood increases, research into the genetic basis of traits that can increase aquaculture production is critical. The eastern oyster (Crassostrea virginica) is an important aquaculture species along the Atlantic and Gulf Coasts of the United States, but increases in heavy rainfall events expose oysters to acute low salinity conditions, which negatively impact production. Low salinity survival is known to be a moderately heritable trait, but the genetic architecture underlying this trait is still poorly understood. In this study, we used ddRAD sequencing to generate genome-wide single-nucleotide polymorphism (SNP) data for four F₂ families to investigate the genomic regions associated with survival in extreme low salinity (<3). SNP data were also used to assess the feasibility of genomic selection (GS) for improving this trait. Quantitative trait locus (QTL) mapping and combined linkage disequilibrium analysis revealed significant QTL on eastern oyster chromosomes 1 and 7 underlying both survival and day to death in a 36-day experimental challenge. Significant QTL were located in genes related to DNA/RNA function and repair, ion binding and membrane transport, and general response to stress. GS was investigated using Bayesian linear regression models and prediction accuracies ranged from 0.48 to 0.57. Genomic prediction accuracies were largest using the BayesB prior and prediction accuracies did not substantially decrease when SNPs located within the QTL region on Chr1 were removed, suggesting that this trait is controlled by many genes of small effect. Our results suggest that GS will likely be a viable option for improvement of survival in extreme low salinity.

Miro proteins connect mitochondrial function and intercellular transport

Mitochondria are organelles present in most eukaryotic cells, where they play major and multifaceted roles. The classical notion of the main mitochondrial function as the powerhouse of the cell per se has been complemented by recent discoveries pointing to mitochondria as organelles affecting a number of other auxiliary processes. They go beyond the classical energy provision via acting as a relay point of many catabolic and anabolic processes, to signaling pathways critically affecting cell growth by their implication in de novo pyrimidine synthesis. These additional roles further underscore the importance of mitochondrial homeostasis in various tissues, where its deregulation promotes a number of pathologies. While it has long been known that mitochondria can move within a cell to sites where they are needed, recent research has uncovered that mitochondria can also move between cells. While this intriguing field of research is only emerging, it is clear that mobilization of mitochondria requires a complex apparatus that critically involves mitochondrial proteins of the Miro family, whose role goes beyond the mitochondrial transfer, as will be covered in this review.

Rho Family of Ras-Like GTPases in Early-Branching Animals

Non-bilaterian animals consist of four phyla; Porifera, Cnidaria, Ctenophora, and Placozoa. These early-diverging animals are crucial for understanding the evolution of the entire animal lineage. The Rho family of proteins make up a major branch of the Ras superfamily of small GTPases, which function as key molecular switches that play important roles in converting and amplifying external signals into cellular responses. This review represents a compilation of the current knowledge on Rho-family GTPases in non-bilaterian animals, the available experimental data about their biochemical characteristics and functions, as well as original bioinformatics analysis, in order to gain a general insight into the evolutionary history of Rho-family GTPases in simple animals.

Proteomic profile in the mussel Perna viridis after short-term exposure to the brown tide alga Aureococcus anophagefferens

Blooms of Aureococcus anophagefferens, referred to as brown tides are responsible for massive mortalities and recruitment failure of some bivalves. However, the molecular mechanisms underlying the toxicity remain elusive despite its biological significance, and the information currently available on the molecular effects is still insufficient. In this study, to evaluate the toxicity and associated mechanism of A. anophagefferens on bivalves, we analyzed the protein expression profiles in digestive glands of the A. anophagefferens-exposed Perna viridis by using iTRAQ. A total of 3138 proteins were identified in the digestive glands of A. anophagefferens-exposed P. viridis based on iTRAQ. Amongst, a repertoire of 236 proteins involved in cell, cell part, catalytic activity, metabolic process, biological regulation, immune system process, and response to stimulus were found to be differentially expressed. Functional analysis of the differentially expressed proteins demonstrated that innate immune system of P. viridis was activated, and some proteins associated with stress response and lipid metabolism were induced after exposure to A. anophagefferens. Additionally, MDA content, SOD activity and GSH-Px activity was increased significantly in the digestive gland of A. anophagefferens-exposed P. viridis. Taken together, our results indicated that the A. anophagefferens could induce oxidative stress, activate complement system and alter fat acid metabolism of P. viridis.

Functional and molecular responses of the blue mussel Mytilus edulis' hemocytes exposed to cadmium - An in vitro model and transcriptomic approach

The bivalve mollusk, Mytilus edulis, is used as a sentinel species in several monitoring programs due to its ability to bio-accumulate contaminants. Its immune system consists of hemocytes and humoral components, which constitute the main part of the hemolymph. The present study is aimed at understanding the effects of Cd on the differentially expressed genes involved in the phagocytosis of M. edulis' hemocytes. Our approach focuses on an in vitro model by exposing hemocytes to different concentrations of Cd ranging from 10^-9 M to 10^-3 M. Phagocytosis and cell viability as functional markers were measured using flow cytometry. The molecular mechanisms regulated by Cd were investigated using RNA-seq and DGE analysis. Results showed that viability and phagocytosis of hemocytes exposed to 10^-3 M of Cd were significantly decreased after 21 h of exposure. RNA sequencing data showed that 1112 transcripts (out of 352,976 contigs) were differentially regulated by the highest concentration of Cd. Among these identified transcripts, 1028 and 84 were up and down-regulated respectively. The induction of super oxide dismutase (SOD), glutathion-s-transferase (GST), cytochrome P450 2C8 (CYP2C8), multidrug resistance protein (MRP1) and heat shock protein 70 (HSP70) suggests that Cd can regulate key molecular mechanisms. In addition, several toll-like receptors (TLR) as well as genes involved in phagocytosis (actin and CDC42) and apoptosis (caspase 8 and XIAP/IAP) were induced by Cd. Thus, our model highlights the effect of Cd on the phagocytic function of M. edulis' hemocytes along with the regulation of gene expression involved in innate immunity, detoxification and apoptosis. Further investigations need to be pursued to unravel the effects of Cd on the molecular mechanisms identified in this study.