Search

The murine submandibular gland (SMG) is a model organ to study development, because it follows a branching morphogenesis pattern that is similar to that of lung, kidney, and other systems. It has been speculated that through its study, insights into regeneration and cancer could be obtained. Previously, using bulk RNA-Seq data, we reported that Transposable Elements (TEs) become activated during the SMG development. However, an outstanding question was as to whether their activity influenced different cell populations. Here, taking advantage of a single cell RNA-Seq atlas of the developing SMG, I studied TE expression to find out whether their activity can be recapitulated across its development, and if so, how they influenced cell types and cell fate specification. In this work, I found a total of 339 TEs that are markers of different cell populations, and then, through the modeling of the SMG development using Trajectory Inference methods, I found 2 TEs that could be potentially influencing differentiation processes. In sum, this short report reveals that TEs may be involved in the normal development of the SMG, and it highlights the importance of considering them in scRNA-Seq studies.

Amyotrophic lateral sclerosis (ALS) is a disease that progressively annihilates spinal cord motor neurons, causing severe motor decline and death. The disease is divided into familial and sporadic ALS. Mutations in the TAR DNA binding protein 43 (TDP-43) have been involved in the pathological emergence and progression of ALS, although the molecular mechanisms eliciting the disease are unknown. Transposable elements (TEs) and DNA sequences capable of transposing within the genome become dysregulated and transcribed in the presence of TDP-43 mutations. We performed RNA-Seq in human motor neurons (iMNs) derived from induced pluripotent stem cells (iPSCs) from TDP-43 wild-type—iMNs-TDP-43WT—and mutant—iMNs-TDP-43M337V—genotypes at 7 and 14 DIV, and, with state-of-the-art bioinformatic tools, analyzed whether TDP-43M337V alters both gene expression and TE activity. Our results show that TDP-43M337V induced global changes in the gene expression and TEs levels at all in vitro stages studied. Interestingly, many genetic pathways overlapped with that of the TEs activity, suggesting that TEs control the expression of several genes. TEs correlated with genes that played key roles in the extracellular matrix and RNA processing: all the regulatory pathways affected in ALS. Thus, the loss of TE regulation is present in TDP-43 mutations and is a critical determinant of the disease in human motor neurons. Overall, our results support the evidence that indicates TEs are critical regulatory sequences contributing to ALS neurodegeneration.

High cholesterol levels have been associated with cardiovascular diseases, and lowering them has been a key focus in the treatment of such diseases. Statins are drugs used with that aim, and can be divided in the lipophilic Simvastatin and the hydrophilic Rosuvastatin. Regardless of the statin type, a high proportion (∼70%) of patients stop using statins due to suffering from side effects on skeletal muscle, such as myalgia, and muscle cramps. Thus, there has been a considerable effort in understanding how statins contribute to these side effects. A catalogue of genes and molecular pathways that change upon statin treatment has been recently published, allowing further understanding how the side effects occur. However, Transposable Elements (TEs) were not studied. TEs can move within a genome, and they are highly repetitive, representing about half of the human genome. Currently, most TEs in the human genome are inactive, but it has been shown that TEs can still transcribe, and that either via their transposition or their transcriptional activity, can influence gene expression. Here, using novel computational tools to accurately estimate TE expression, we studied their activity and predicted their potential impact on gene expression. We developed a catalogue of TEs expressed upon statin treatment, and the putative genes whose expression might be influenced by TEs. Overall, we speculate that based on our findings, TEs might be a key target in order to understand statin-mediated side effects.

Spatial transcriptomics (ST) is transforming the way we can study gene expression and its regulation through position-specific resolution within tissues. However, as in bulk RNA-Seq, transposable elements (TEs) are not being studied due to their highly repetitive nature. In recent years, TEs have been recognized as important regulators of gene expression, and thus, TE expression analysis in a spatially resolved manner could further help to understand their role in gene regulation within tissues. We present SpatialTE, a tool to analyze TE expression from ST datasets and show its application in somatic and diseased tissues. The results indicate that TEs have spatially regulated expression patterns and that their expression profiles are spatially altered in ALS disease, indicating that TEs might perform differential regulatory functions within tissue organs. We have made SpatialTE publicly available as open-source software under an MIT license.

The submandibular gland (SG) is a relatively simple organ formed by three cell types: acinar, myoepithelial, and an intricate network of duct-forming epithelial cells, that together fulfills several physiological functions from assisting food digestion to acting as an immune barrier against pathogens. Successful SG organogenesis is the product of highly controlled and orchestrated genetic and transcriptional programs. Mounting evidence links Transposable Elements (TEs), originally thought to be selfish genetic elements, to different aspects of gene regulation in mammalian development and disease. To our knowledge, the role of TEs during murine SG organogenesis has not been studied. Using novel bioinformatic tools and publicly available RNA-Seq datasets, our results indicate that a significant number of genic and intergenic TEs are differentially expressed during the SG development. Furthermore, changes in expression of specific TEs correlated with that of genes involved in cellular division and differentiation, critical aspects for SG maturation. Altogether, we propose that TEs modulate gene networks that operate during SG development. Supplementary Information The online version contains supplementary material available at 10.1186/s13100-021-00251-1.

The first step in any genome research after obtaining the read data is to perform a due quality control of the sequenced reads. In a de novo genome assembly project, the second step is to estimate two important features, the genome size and ‘best k -mer’, to start the assembly tests with different de novo assembly software and its parameters. However, the quality control of the sequenced genome libraries as a whole, instead of focusing on the reads only, is frequently overlooked and realized to be important only when the assembly tests did not render the expected results. We have developed GSER, a Genome Size Estimator using R, a pipeline to evaluate the relationship between k -mers and genome size, as a means for quality assessment of the sequenced genome libraries. GSER generates a set of charts that allow the analyst to evaluate the library datasets before starting the assembly. The script which runs the pipeline can be downloaded from http://www.mobilomics.org/GSER/downloads or http://github.com/mobilomics/GSER .

HKT channels are a plant protein family involved in sodium (Na+) and potassium (K+) uptake and Na+-K+ homeostasis. Some HKTs underlie salt tolerance responses in plants, while others provide a mechanism to cope with short-term K+ shortage by allowing increased Na+ uptake under K+ starvation conditions. HKT channels present a functionally versatile family divided into two classes, mainly based on a sequence polymorphism found in the sequences underlying the selectivity filter of the first pore loop. Physiologically, most class I members function as sodium uniporters, and class II members as Na+/K+ symporters. Nevertheless, even within these two classes, there is a high functional diversity that, to date, cannot be explained at the molecular level. The high complexity is also reflected at the regulatory level. HKT expression is modulated at the level of transcription, translation, and functionality of the protein. Here, we summarize and discuss the structure and conservation of the HKT channel family from algae to angiosperms. We also outline the latest findings on gene expression and the regulation of HKT channels.

The early development stages of the frog Xenopus laevis are known as the regenerative (R) stage, because during that time, it is able to regenerate. After its metamorphosis into a frog, X. laevis goes into its non-regenerative (NR) stage. It has been reported that there are differences in the gene repertoire expressed at each stage after spinal cord injury, and also in its timing. One molecular agent that might influence gene expression and its timing are small non-coding RNAs (snc-RNAs). To gain insights into the role of snc-RNAs, in this work we aimed to investigate their differential expression between the R and the NR stages and their origins in genome regions. We have performed small RNA sequencing and differential expression analysis between the R and the NR stages. The majority of sequenced snc-RNAs have high quality, and their lengths suggest that the majority belong to the small interfering or, to the micro RNA class. In average, 98% of the ~50 million sequences completely aligned to X. laevis genome, and ~45% of them represent novel snc-RNAs. The snc-RNAs were further classified in 2.238 families, using as criteria their genomic location: ~49% of novel snc-RNA are in TEs regions. ~53% of families have their origins in genes (~33% in intronic regions) and ~47% in intergenic regions. 374 novel snc-RNAs are differentially expressed. 289 come from TEs, 250 of which are in intergenic regions and 39 are in intronic regions. These results, taken together with the differentially expressed genes, will help understanding the spinal cord regeneration process in frogs.

Motivation In recent years, Transposable Elements (TEs) have been related to gene regulation. However, estimating the origin of expression of TEs through RNA-seq is complicated by multi-mapping reads coming from their repetitive sequences. Current approaches that address multi-mapping reads are focused in expression quantification and not in finding the origin of expression. Addressing the genomic origin of expressed TEs could further aid in understanding the role that TEs might have in the cell. Results We have developed a new pipeline called TEcandidates, based on de novo transcriptome assembly to assess the instances of TEs being expressed, along with their location, to include in downstream DE analysis. TEcandidates takes as input the RNA-seq data, the genome sequence and the TE annotation file and returns a list of coordinates of candidate TEs being expressed, the TEs that have been removed and the genome sequence with removed TEs as masked. This masked genome is suited to include TEs in downstream expression analysis, as the ambiguity of reads coming from TEs is significantly reduced in the mapping step of the analysis. Availability and implementation The script which runs the pipeline can be downloaded at http://www.mobilomics.org/tecandidates/downloads or http://github.com/TEcandidates/TEcandidates. Supplementary information Supplementary data are available at Bioinformatics online.

The electrostatic potential plays a key role in many biological processes like determining the affinity of a ligand to a given protein target, and they are responsible for the catalytic activity of many enzymes. Understanding the effect that amino acid mutations will have on the electrostatic potential of a protein, will allow a thorough understanding of which residues are the most important in a protein. MutantElec, is a friendly web application for in silico generation of site-directed mutagenesis of proteins and the comparison of electrostatic potential between the wild type protein and the mutant(s), based on the three-dimensional structure of the protein. The effect of the mutation is evaluated using different approach to the traditional surface map. MutantElec provides a graphical display of the results that allows the visualization of changes occurring at close distance from the mutation and thus uncovers the local and global impact of a specific change

In arbuscular mycorrhizal (AM) symbiosis, fungi and plants exchange nutrients (sugars and phosphate, for instance) for reciprocal benefit. Until now it is not clear how this nutrient exchange system works. Here, we used computational cell biology to simulate the dynamics of a network of proton pumps and proton-coupled transporters that are upregulated during AM formation. We show that this minimal network is sufficient to describe accurately and realistically the nutrient trade system. By applying basic principles of microeconomics, we link the biophysics of transmembrane nutrient transport with the ecology of organismic interactions and straightforwardly explain macroscopic scenarios of the relations between plant and AM fungus. This computational cell biology study allows drawing far reaching hypotheses about the mechanism and the regulation of nutrient exchange and proposes that the ‘cooperation’ between plant and fungus can be in fact the result of a competition between both for the same resources in the tiny periarbuscular space. The minimal model presented here may serve as benchmark to evaluate in future the performance of more complex models of AM nutrient exchange. As a first step towards this goal, we included SWEET sugar transporters in the model and show that their co-occurrence with proton-coupled sugar transporters results in a futile carbon cycle at the plant plasma membrane proposing that two different pathways for the same substrate should not be active at the same time.

Two-pore domain potassium (K2P) channels are membrane proteins widely identified in mammals, plants, and other organisms. A functional channel is a dimer with each subunit comprising two pore-forming loops and four transmembrane domains. The genome of the model plant Arabidopsis thaliana harbors five genes coding for K2P channels. Homologs of Arabidopsis K2P channels have been found in all higher plants sequenced so far. As with the K2P channels in mammals, plant K2P channels are targets of external and internal stimuli, which fine-tune the electrical properties of the membrane for specialized transport and/or signaling tasks. Plant K2P channels are modulated by signaling molecules such as intracellular H+ and calcium and physical factors like temperature and pressure. In this review, we ask the following: What are the similarities and differences between K2P channels in plants and animals in terms of their physiology? What is the nature of the last common ancestor (LCA) of these two groups of proteins? To answer these questions, we present physiological, structural, and phylogenetic evidence that discards the hypothesis proposing that the duplication and fusion that gave rise to the K2P channels occurred in a prokaryote LCA. Conversely, we argue that the K2P LCA was most likely a eukaryote organism. Consideration of plant and animal K2P channels in the same study is novel and likely to stimulate further exchange of ideas between students of these fields.

Two-pore domain potassium (K2P) channels are membrane proteins widely identified in mammals, plants, and other organisms. A functional channel is a dimer with each subunit comprising two pore-forming loops and four transmembrane domains. The genome of the model plant Arabidopsis thaliana harbors five genes coding for K2P channels. Homologs of Arabidopsis K2P channels have been found in all higher plants sequenced so far. As with the K2P channels in mammals, plant K2P channels are targets of external and internal stimuli, which fine-tune the electrical properties of the membrane for specialized transport and/or signaling tasks. Plant K2P channels are modulated by signaling molecules such as intracellular H(+) and calcium and physical factors like temperature and pressure. In this review, we ask the following: What are the similarities and differences between K2P channels in plants and animals in terms of their physiology? What is the nature of the last common ancestor (LCA) of these two groups of proteins? To answer these questions, we present physiological, structural, and phylogenetic evidence that discards the hypothesis proposing that the duplication and fusion that gave rise to the K2P channels occurred in a prokaryote LCA. Conversely, we argue that the K2P LCA was most likely a eukaryote organism. Consideration of plant and animal K2P channels in the same study is novel and likely to stimulate further exchange of ideas between students of these fields.