Extracellular Vesicles Secreted by Pre-Hatching Bovine Embryos Produced In Vitro and In Vivo Alter the Expression of IFNtau-Stimulated Genes in Bovine Endometrial Cells.

The embryo-maternal interaction occurs during the early stages of embryo development and is essential for the implantation and full-term development of the embryo. In bovines, the secretion of interferon Tau (IFNT) during elongation is the main signal for pregnancy recognition, but its expression starts around the blastocyst stage. Embryos release extracellular vesicles (EVs) as an alternative mechanism of embryo-maternal communication. The aim of the study was to determine whether EVs secreted by bovine embryos during blastulation (D5-D7) could induce transcriptomic modifications, activating IFNT signaling in endometrial cells. Additionally, it aims to assess whether the EVs secreted by embryos produced in vivo (EVs-IVV) or in vitro (EVs-IVP) have different effects on the transcriptomic profiles of the endometrial cells. In vitro- and in vivo-produced bovine morulae were selected and individually cultured for 48 h to collect embryonic EVs (E-EVs) secreted during blastulation. E-EVs stained with PKH67 were added to in vitro-cultured bovine endometrial cells to assess EV internalization. The effect of EVs on the transcriptomic profile of endometrial cells was determined by RNA sequencing. EVs from both types of embryos induced several classical and non-classical IFNT-stimulated genes (ISGs) and other pathways related to endometrial function in epithelial endometrial cells. Higher numbers of differentially expressed genes (3552) were induced by EVs released by IVP embryos compared to EVs from IVV (1838). Gene ontology analysis showed that EVs-IVP/IVV induced the upregulation of the extracellular exosome pathway, the cellular response to stimulus, and the protein modification processes. This work provides evidence regarding the effect of embryo origin (in vivo or in vitro) on the early embryo-maternal interaction mediated by extracellular vesicles.

Bioinformatic analysis predicts that ethanol exposure during early development causes alternative splicing alterations of genes involved in RNA post-transcriptional regulation.

Prenatal ethanol exposure is associated with neurodevelopmental defects and long-lasting cognitive deficits, which are grouped as fetal alcohol spectrum disorders (FASD). The molecular mechanisms underlying FASD are incompletely characterized. Alternative splicing, including the insertion of microexons (exons of less than 30 nucleotides in length), is highly prevalent in the nervous system. However, whether ethanol exposure can have acute or chronic deleterious effects in this process is poorly understood. In this work, we used the bioinformatic tools VAST-TOOLS, rMATS, MAJIQ, and MicroExonator to predict alternative splicing events affected by ethanol from available RNA sequencing data. Experimental protocols of ethanol exposure included human cortical tissue development, human embryoid body differentiation, and mouse development. We found common genes with predicted differential alternative splicing using distinct bioinformatic tools in different experimental designs. Notably, Gene Ontology and KEGG analysis revealed that the alternative splicing of genes related to RNA processing and protein synthesis was commonly affected in the different ethanol exposure schemes. In addition, the inclusion of microexons was also affected by ethanol. This bioinformatic analysis provides a reliable list of candidate genes whose splicing is affected by ethanol during nervous system development. Furthermore, our results suggest that ethanol particularly modifies the alternative splicing of genes related to post-transcriptional regulation, which probably affects neuronal proteome complexity and brain function.

Bovine embryos release extracellular vesicles with differential miRNA signature during the compaction and blastulation stages.

Pre-implantation embryos release extracellular vesicles (EVs) to extracellular environment. In this work it is hypothesized that the EVs miRNA cargo will vary during pre-implantation development due to the constant changes in gene expression that take place through this period. The concentration, size and miRNA cargo of EVs secreted by competent bovine embryos during the period from compaction to blastulation (Day 3-7) were analyzed. For this analysis tow developmental windows were defined: W2 from 8-cells (D3) to morula (D5) and W3 from morula (D5) to blastocyst (D7). For W2, in vitro produced embryos were individually cultured in EVs-depleted medium from D3 to D5; culture media were collected and assigned to Group W2. Morulae were kept in culture up to blastocyst stage to determine the developmental competence. For W3, D5 morulae were collected and cultured individually in EVs-depleted medium up to blastocyst stage; culture media were assigned to Group W3, and blastocysts were kept in culture up to day 11 to define their competence. The mean size of EVs was similar between groups, however, EVs concentration was lower in W2. A total of 140 miRNAs were identified. From them, 79 were differentially expressed between the groups, 28 upregulated and 51 downregulated. miRNAs differentially detected between both developmental windows participate in the regulation of signaling pathways which crucial for embryonic development. It was concluded that the secretion of EVs is regulated by the developmental progress of the embryo during the pre-implantation period.

Lymph Leakage Promotes Immunosuppression by Enhancing Anti-Inflammatory Macrophage Polarization.

Lymphatic vasculature is a network of capillaries and vessels capable of draining extracellular fluid back to blood circulation and to facilitate immune cell migration. Although the role of the lymphatic vasculature as coordinator of fluid homeostasis has been extensively studied, the consequences of abnormal lymphatic vasculature function and impaired lymph drainage have been mostly unexplored. Here, by using the Prox1+/- mice with defective lymphatic vasculature and lymphatic leakage, we provide evidence showing that lymph leakage induces an immunosuppressive environment by promoting anti-inflammatory M2 macrophage polarization in different inflammatory conditions. In fact, by using a mouse model of tail lymphedema where lymphatic vessels are thermal ablated leading to lymph accumulation, an increasing number of anti-inflammatory M2 macrophages are found in the lymphedematous tissue. Moreover, RNA-seq analysis from different human tumors shows that reduced lymphatic signature, a hallmark of lymphatic dysfunction, is associated with increased M2 and reduced M1 macrophage signatures, impacting the survival of the patients. In summary, we show that lymphatic vascular leakage promotes an immunosuppressive environment by enhancing anti-inflammatory macrophage differentiation, with relevance in clinical conditions such as inflammatory bowel diseases or cancer.

Hidden Markov Modeling with HMMTeacher.

Is it possible to learn and create a first Hidden Markov Model (HMM) without programming skills or understanding the algorithms in detail? In this concise tutorial, we present the HMM through the 2 general questions it was initially developed to answer and describe its elements. The HMM elements include variables, hidden and observed parameters, the vector of initial probabilities, and the transition and emission probability matrices. Then, we suggest a set of ordered steps, for modeling the variables and illustrate them with a simple exercise of modeling and predicting transmembrane segments in a protein sequence. Finally, we show how to interpret the results of the algorithms for this particular problem. To guide the process of information input and explicit solution of the basic HMM algorithms that answer the HMM questions posed, we developed an educational webserver called HMMTeacher. Additional solved HMM modeling exercises can be found in the user's manual and answers to frequently asked questions. HMMTeacher is available at https://hmmteacher.mobilomics.org, mirrored at https://hmmteacher1.mobilomics.org. A repository with the code of the tool and the webpage is available at https://gitlab.com/kmilo.f/hmmteacher.

A brain-specific pgc1α fusion transcript affects gene expression and behavioural outcomes in mice.

PGC1α is a transcriptional coactivator in peripheral tissues, but its function in the brain remains poorly understood. Various brain-specific Pgc1α isoforms have been reported in mice and humans, including two fusion transcripts (FTs) with non-coding repetitive sequences, but their function is unknown. The FTs initiate at a simple sequence repeat locus ∼570 Kb upstream from the reference promoter; one also includes a portion of a short interspersed nuclear element (SINE). Using publicly available genomics data, here we show that the SINE FT is the predominant form of Pgc1α in neurons. Furthermore, mutation of the SINE in mice leads to altered behavioural phenotypes and significant up-regulation of genes in the female, but not male, cerebellum. Surprisingly, these genes are largely involved in neurotransmission, having poor association with the classical mitochondrial or antioxidant programs. These data expand our knowledge on the role of Pgc1α in neuronal physiology and suggest that different isoforms may have distinct functions. They also highlight the need for further studies before modulating levels of Pgc1α in the brain for therapeutic purposes.

Locus-specific analysis of Transposable Elements during the progression of ALS in the SOD1G93A mouse model.

Transposable Elements (TEs) are ubiquitous genetic elements with the ability to move within a genome. TEs contribute to a large fraction of the repetitive elements of a genome, and because of their nature, they are not routinely analyzed in RNA-Seq gene expression studies. Amyotrophic Lateral Sclerosis (ALS) is a lethal neurodegenerative disease, and a well-accepted model for its study is the mouse harboring the human SOD1G93A mutant. In this model, landmark stages of the disease can be recapitulated at specific time points, making possible to understand changes in gene expression across time. While there are several works reporting TE activity in ALS models, they have not explored their activity through the disease progression. Moreover, they have done it at the expense of losing their locus of expression. Depending on their genomic location, TEs can regulate genes in cis and in trans, making locus-specific analysis of TEs of importance in order to understand their role in modulating gene expression. Particularly, the locus-specific role of TEs in ALS has not been fully elucidated. In this work, we analyzed publicly available RNA-Seq datasets of the SOD1G93A mouse model, to understand the locus-specific role of TEs. We show that TEs become up-regulated at the early stages of the disease, and via statistical associations, we speculate that they can regulate several genes, which in turn might be contributing to the genetic dysfunction observed in ALS.

GSER (a Genome Size Estimator using R): a pipeline for quality assessment of sequenced genome libraries through genome size estimation.

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.

How to Grow a Tree: Plant Voltage-Dependent Cation Channels in the Spotlight of Evolution.

Phylogenetic analysis can be a powerful tool for generating hypotheses regarding the evolution of physiological processes. Here, we provide an updated view of the evolution of the main cation channels in plant electrical signalling: the Shaker family of voltage-gated potassium channels and the two-pore cation (K+) channel (TPC1) family. Strikingly, the TPC1 family followed the same conservative evolutionary path as one particular subfamily of Shaker channels (Kout) and remained highly invariant after terrestrialisation, suggesting that electrical signalling was, and remains, key to survival on land. We note that phylogenetic analyses can have pitfalls, which may lead to erroneous conclusions. To avoid these in the future, we suggest guidelines for analyses of ion channel evolution in plants.

MicroRNAs from Extracellular Vesicles Secreted by Bovine Embryos as Early Biomarkers of Developmental Competence.

During early development, embryos secrete extracellular vesicles (EVs) that participate in embryo-maternal communication. Among other molecules, EVs carry microRNAs (miRNAs) that interfere with gene expression in target cells; miRNAs participate in embryo-maternal communication. Embryo selection based on secreted miRNAs may have an impact on bovine breeding programs. This research aimed to evaluate the size, concentration, and miRNA content of EVs secreted by bovine embryos with different developmental potential, during the compaction period (days 3.5-5). Individual culture media from in vitro-produced embryos were collected at day 5, while embryos were further cultured and classified at day 7, as G1 (conditioned-culture media by embryos arrested in the 8-16-cells stage) and G2 (conditioned-culture media by embryos that reached blastocyst stages at day 7). Collected nanoparticles from embryo conditioned culture media were cataloged as EVs by their morphology and the presence of classical molecular markers. Size and concentration of EVs from G1 were higher than EVs secreted by G2. We identified 95 miRNAs; bta-miR-103, bta-miR-502a, bta-miR-100, and bta-miR-1 were upregulated in G1, whereas bta-miR-92a, bta-miR-140, bta-miR-2285a, and bta-miR-222 were downregulated. The most significant upregulated pathways were fatty acid biosynthesis and metabolism, lysine degradation, gap junction, and signaling pathways regulating pluripotency of stem cells. The characteristics of EVs secreted by bovine embryos during the compaction period vary according to embryo competence. Embryos that reach the blastocyst stage secrete fewer and smaller vesicles. Furthermore, the loading of specific miRNAs into the EVs depends on embryo developmental competence.

Distinctive Cellular Transcriptomic Signature and MicroRNA Cargo of Extracellular Vesicles of Horse Adipose and Endometrial Mesenchymal Stem Cells from the Same Donors.

Equine endometrial and adipose mesenchymal stem cells (eMSCs and aMSCs, respectively) were isolated from the same donors of thoroughbred mares. The cells displayed characteristic features of MSCs, including trilineage mesodermal and also neurogenic differentiation. We evaluated the influence of cellular origin on their transcriptome profile. Cellular RNA was isolated and sequenced and extracellular vesicles (EVs) were obtained from conditioned medium of cells cultured in medium depleted of EVs, and their microRNA (miRNA) cargo analyzed by sequencing. Differential expression of mRNAs and EV-miRNA was analyzed, as well as pathways and processes most represented in each cell origin. mRNA reads from all expressed genes clustered according to the cellular origin. A total of 125 up- and 51 downregulated genes were identified and 31 differentially expressed miRNAs. Based on mRNA sequencing, endometrial MSCs strongly upregulated genes involved in the Hippo, transforming growth factor beta, and pluripotency signaling pathways. Alongside with this, pathways involved in extracellular matrix reorganization were the most represented in the miRNA cargo of EVs secreted by eMSCs. The niche from which MSCs originated defined the transcriptomic signature of the cells, including the secretion of lineage-specific loaded EV to ensure proper communication and homeostasis. Identification and testing their biological functions can provide new tools for the therapeutic use of horse MSC.

Mitochondrial acetyl-CoA reversibly regulates locus-specific histone acetylation and gene expression.

The impact of mitochondrial dysfunction in epigenetics is emerging, but our understanding of this relationship and its effect on gene expression remains incomplete. We previously showed that acute mitochondrial DNA (mtDNA) loss leads to histone hypoacetylation. It remains to be defined if these changes are maintained when mitochondrial dysfunction is chronic and if they alter gene expression. To fill these gaps of knowledge, we here studied a progressive and a chronic model of mtDNA depletion using biochemical, pharmacological, genomics, and genetic assays. We show that histones are primarily hypoacetylated in both models. We link these effects to decreased histone acetyltransferase activity unrelated to changes in ATP citrate lyase, acetyl coenzyme A synthetase 2, or pyruvate dehydrogenase activities, which can be reversibly modulated by altering the mitochondrial pool of acetyl-coenzyme A. Also, we determined that the accompanying changes in histone acetylation regulate locus-specific gene expression and physiological outcomes, including the production of prostaglandins. These results may be relevant to the pathophysiology of mtDNA depletion syndromes and to understanding the effects of environmental agents that lead to physical or functional mtDNA loss.

TEcandidates: prediction of genomic origin of expressed transposable elements using RNA-seq data.

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 receptor-like pseudokinase MRH1 interacts with the voltage-gated potassium channel AKT2.

The potassium channel AKT2 plays important roles in phloem loading and unloading. It can operate as inward-rectifying channel that allows H+-ATPase-energized K+ uptake. Moreover, through reversible post-translational modifications it can also function as an open, K+-selective channel, which taps a 'potassium battery', providing additional energy for transmembrane transport processes. Knowledge about proteins involved in the regulation of the operational mode of AKT2 is very limited. Here, we employed a large-scale yeast two-hybrid screen in combination with fluorescence tagging and null-allele mutant phenotype analysis and identified the plasma membrane localized receptor-like kinase MRH1/MDIS2 (AT4G18640) as interaction partner of AKT2. The phenotype of the mrh1-1 knockout plant mirrors that of akt2 knockout plants in energy limiting conditions. Electrophysiological analyses showed that MRH1/MDIS2 failed to exert any functional regulation on AKT2. Using structural protein modeling approaches, we instead gathered evidence that the putative kinase domain of MRH1/MDIS2 lacks essential sites that are indispensable for a functional kinase suggesting that MRH1/MDIS2 is a pseudokinase. We propose that MRH1/MDIS2 and AKT2 are likely parts of a bigger protein complex. MRH1 might help to recruit other, so far unknown partners, which post-translationally regulate AKT2. Additionally, MRH1 might be involved in the recognition of chemical signals.

A Novel Analytical Strategy to Identify Fusion Transcripts between Repetitive Elements and Protein Coding-Exons Using RNA-Seq.

Repetitive elements (REs) comprise 40-60% of the mammalian genome and have been shown to epigenetically influence the expression of genes through the formation of fusion transcript (FTs). We previously showed that an intracisternal A particle forms an FT with the agouti gene in mice, causing obesity/type 2 diabetes. To determine the frequency of FTs genome-wide, we developed a TopHat-Fusion-based analytical pipeline to identify FTs with high specificity. We applied it to an RNA-seq dataset from the nucleus accumbens (NAc) of mice repeatedly exposed to cocaine. Cocaine was previously shown to increase the expression of certain REs in this brain region. Using this pipeline that can be applied to single- or paired-end reads, we identified 438 genes expressing 813 different FTs in the NAc. Although all types of studied repeats were present in FTs, simple sequence repeats were underrepresented. Most importantly, reverse-transcription and quantitative PCR validated the expression of selected FTs in an independent cohort of animals, which also revealed that some FTs are the prominent isoforms expressed in the NAc by some genes. In other RNA-seq datasets, developmental expression as well as tissue specificity of some FTs differed from their corresponding non-fusion counterparts. Finally, in silico analysis predicted changes in the structure of proteins encoded by some FTs, potentially resulting in gain or loss of function. Collectively, these results indicate the robustness of our pipeline in detecting these new isoforms of genes, which we believe provides a valuable tool to aid in better understanding the broad role of REs in mammalian cellular biology.

Towards the bridging of molecular genetics data across Xenopus species.

BACKGROUND: The clawed African frog Xenopus laevis has been one of the main vertebrate models for studies in developmental biology. However, for genetic studies, Xenopus tropicalis has been the experimental model of choice because it shorter life cycle and due to a more tractable genome that does not result from genome duplication as in the case of X. laevis. Today, although still organized in a large number of scaffolds, nearly 85% of X. tropicalis and 89% of X. laevis genomes have been sequenced. There is expectation for a comparative physical map that can be used as a Rosetta Stone between X. laevis genetic studies and X. tropicalis genomic research. RESULTS: In this work, we have mapped using coarse-grained alignment the 18 chromosomes of X. laevis, release 9.1, on the 10 reference scaffolds representing the haploid genome of X. tropicalis, release 9.0. After validating the mapping with theoretical data, and estimating reference averages of genome sequence identity, 37 to 44% between the two species, we have carried out a synteny analysis for 2,112 orthologous genes. We found that 99.6% of genes are in the same organization. CONCLUSIONS: Taken together, our results make possible to establish the correspondence between 62 and 65.5% of both genomes, percentage of identity, synteny and automatic annotation of transcripts of both species, providing a new and more comprehensive tool for comparative analysis of these two species, by allowing to bridge molecular genetics data among them.

Novel missense mutation of the FAM83H gene causes retention of amelogenin and a mild clinical phenotype of hypocalcified enamel.

OBJECTIVE: Amelogenesis imperfecta (AI) is a group of clinically and genetically heterogeneous inherited conditions, causing alterations in the structure of enamel and chemical composition of enamel matrix during development. The objective of this study was to compare the clinical, radiographic, histological and immunohistochemical phenotypes of subjects affected with hypocalcified AI from three Chilean families and identify causal mutations in the FAM83H gene. DESIGN: The diagnosis was made using clinical, radiographic, histological and genealogical data from the patients, who were evaluated according to the classification criteria by Witkop. PCR and Sanger sequencing of the complete coding sequence and surrounding intron regions of the FAM83H gene were conducted. The structural study of the affected teeth was performed with light microscopy, scanning electron microscopy and immunohistochemistry. RESULTS: The probands of the three families were diagnosed with hypocalcified AI, but in only one of them the missense variant p.Gly557Cys was identified. This variant was not present in the SNP database or in 100 healthy controls and segregated with the disease in the affected family. Using light microscopy, a normal prismatic structure was observed in all three cases. However, the ultrastructure was found to be affected in two of the cases, showing persistence of organic matter including amelogenins. CONCLUSIONS: These results suggest that FAM83H missense mutation reported in one of the families analyzed in this study might cause a phenotype of hypocalcified enamel more attenuated with retention of amelogenin.

K2p channels in plants and animals.

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.

Easy Identification of Residues Involved on Structural Differences Between Nonphosphorylated and Phosphorylated CDK2Cyclin A Complexes Using Two-Dimensional Networks.

The structures of proteins in Protein Data Bank (PDB) contain a lot of information that can be revealed through the use of tools to facilitate their organization and analysis. The increase in available structural data of nonphosphorylated and phosphorylated CDK2cyclin A (npCDK2cycA and pCDK2cycA) complexes has enabled a more realistic description of the fine structural details of the interface residues of these proteins. This work reports the application of two-dimensional network representations (TDNRs) to the structures deposited in PDB to distinguish the differences in the surface between both complexes due to phosphorylation. As a result, a detailed map of the hydrogen bonds (HBs) and hydrophobic interactions between the T-loop residues of CDK2 and the residues of cycA that are different among nonphosphorylated and phosphorylated complexes were described. In addition, we found some interesting subtle differences in the CDK2cycA interface between nonphosphorylated and phosphorylated complexes due to residues that are not located at the T-loop of CDK2. We noted that some HB interactions in CDK2cycA complex are reinforced when the CDK2 is phosphorylated.