Accelerated evolution in the human lineage led to gain and loss of transcriptional enhancers in the RBFOX1 locus.

A long-standing goal of evolutionary biology is to decode how changes in gene regulatory networks contribute to human-specific traits. Human accelerated regions (HARs) are prime candidates for driving gene regulatory modifications in human development. The RBFOX1 locus is densely populated with HARs, providing a set of potential regulatory elements that could have changed its expression in the human lineage. Here, we examined the role of RBFOX1-HARs using transgenic zebrafish reporter assays and identified 15 transcriptional enhancers that are active in the developing nervous system, 9 of which displayed differential activity between the human and chimpanzee sequences. The engineered loss of two selected RBFOX1-HARs in knockout mouse models modified Rbfox1 expression at specific developmental stages and tissues in the brain, influencing the expression and splicing of a high number of Rbfox1 target genes. Our results provided insight into the spatial and temporal changes in gene expression driven by RBFOX1-HARs.

Improving genome-wide mapping of nucleosomes in Trypanosome cruzi.

In Trypanosoma cruzi DNA is packaged into chromatin by octamers of histone proteins that form nucleosomes. Transcription of protein coding genes in trypanosomes is constitutive producing polycistronic units and gene expression is primarily regulated post-transcriptionally. However, chromatin organization influences DNA dependent processes. Hence, determining nucleosome position is of uppermost importance to understand the peculiarities found in trypanosomes. To map nucleosomes genome-wide in several organisms, digestion of chromatin with micrococcal nuclease followed by deep sequencing has been applied. Nonetheless, the special requirements for cell manipulation and the uniqueness of the chromatin organization in trypanosomes entails a customized analytical approach. In this work, we adjusted this broadly used method to the hybrid reference strain, CL Brener. Particularly, we implemented an exhaustive and thorough computational workflow to overcome the difficulties imposed by this complex genome. We tested the performance of two aligners, Bowtie2 and HISAT2, and discuss their advantages and caveats. Specifically, we highlight the relevance of using the whole genome as a reference instead of the commonly used Esmeraldo-like haplotype to avoid spurious alignments. Additionally, we show that using the whole genome refines the average nucleosome representation, but also the quality of mapping for every region represented. Moreover, we show that the average nucleosome organization around trans-splicing acceptor site described before, is not just an average since the same chromatin pattern is detected for most of the represented regions. In addition, we extended the study to a non-hybrid strain applying the experimental and analytical approach to Sylvio-X10 strain. Furthermore, we provide a source code for the construction of 2D plots and heatmaps which are easy to adapt to any T. cruzi strain.

Expression of Inhibitory Receptors TIGIT, TIM-3, and LAG-3 on CD4+ T Cells from Patients with Different Clinical Forms of Chronic Chagas Disease.

T cells are central to the adaptive immune response against Trypanosoma cruzi infection. In chronic Chagas disease (CCD), circulating parasite-specific memory T cells show reduced functionality and increased expression of inhibitory receptors as a result of persistent antigenic stimulation. This phenotype has been linked to progression of cardiac pathology, whereas the presence of polyfunctional T cells shows association with therapeutic success. In this study, we demonstrate that T. cruzi-specific human CD4+ T cells can be identified by their expression of OX40 and CD25 upon in vitro stimulation. We characterized the expression of the inhibitory receptors T cell immunoreceptor with Ig and ITIM domains (TIGIT), T cell Ig and mucin-domain containing-3 (TIM-3), and lymphocyte activation gene 3 (LAG-3) in CD4+ T cells from CCD patients with and without cardiac alterations. Our results show that, independently of their clinical stage, CCD patients present an increased frequency of CD4+ T cells expressing TIGIT in comparison with non-T. cruzi-infected donors. Exposure to parasite Ags increases the expression of TIM-3 in CD4+ T cells from CCD patients, especially in those with cardiac compromise. Upregulation of LAG-3 was also detected in CCD individuals without cardiac manifestations, predominantly within the subpopulation of cells that did not become activated upon stimulation. Further differences were found between groups in the coexpression of these receptors. Blockade of each individual receptor did not affect activation or the production of IFN-γ and IL-10 by CD4+ T cells in response to parasite Ags. Our results suggest a role for TIGIT, TIM-3, and LAG-3 in the modulation of inflammatory phenomena thought to ultimately lead to tissue damage and cardiac pathology.

A novel strategy to uncover specific GO terms/phosphorylation pathways in phosphoproteomic data in Arabidopsis thaliana.

BACKGROUND: Proteins are the workforce of the cell and their phosphorylation status tailors specific responses efficiently. One of the main challenges of phosphoproteomic approaches is to deconvolute biological processes that specifically respond to an experimental query from a list of phosphoproteins. Comparison of the frequency distribution of GO (Gene Ontology) terms in a given phosphoproteome set with that observed in the genome reference set (GenRS) is the most widely used tool to infer biological significance. Yet, this comparison assumes that GO term distribution between the phosphoproteome and the genome are identical. However, this hypothesis has not been tested due to the lack of a comprehensive phosphoproteome database. RESULTS: In this study, we test this hypothesis by constructing three phosphoproteome databases in Arabidopsis thaliana: one based in experimental data (ExpRS), another based in in silico phosphorylation protein prediction (PredRS) and a third that is the union of both (UnRS). Our results show that the three phosphoproteome reference sets show default enrichment of several GO terms compared to GenRS, indicating that GO term distribution in the phosphoproteomes does not match that of the genome. Moreover, these differences overshadow the identification of GO terms that are specifically enriched in a particular condition. To overcome this limitation, we present an additional comparison of the sample of interest with UnRS to uncover GO terms specifically enriched in a particular phosphoproteome experiment. Using this strategy, we found that mRNA splicing and cytoplasmic microtubule compounds are important processes specifically enriched in the phosphoproteome of dark-grown Arabidopsis seedlings. CONCLUSIONS: This study provides a novel strategy to uncover GO specific terms in phosphoproteome data of Arabidopsis that could be applied to any other organism. We also highlight the importance of specific phosphorylation pathways that take place during dark-grown Arabidopsis development.

Creating 2D Occupancy Plots Using plot2DO.

Chromatin organization and epigenetic marks play a critical role in stem cell pluripotency and differentiation. Chromatin digestion by micrococcal nuclease (MNase) followed by high-throughput sequencing (MNase-seq) is the most widely used genome-wide method for studying nucleosome organization, that is, the first level of DNA packaging into chromatin. Combined with chromatin immunoprecipitation (ChIP), MNase-ChIP-seq represents a high-resolution method for investigating both chromatin organization and the distribution of epigenetic marks and histone variants. The plot2DO package presented here is a flexible tool for evaluating the quality of MNase-seq and MNase-ChIP-seq data, and for visualizing the distribution of nucleosomes near the functional regions of the genome. The plot2DO package is open-source software, and it is freely available from https://github.com/rchereji/plot2DO under the MIT license.

Convergent evolution of two mammalian neuronal enhancers by sequential exaptation of unrelated retroposons.

The proopiomelanocortin gene (POMC) is expressed in a group of neurons present in the arcuate nucleus of the hypothalamus. Neuron-specific POMC expression in mammals is conveyed by two distal enhancers, named nPE1 and nPE2. Previous transgenic mouse studies showed that nPE1 and nPE2 independently drive reporter gene expression to POMC neurons. Here, we investigated the evolutionary mechanisms that shaped not one but two neuron-specific POMC enhancers and tested whether nPE1 and nPE2 drive identical or complementary spatiotemporal expression patterns. Sequence comparison among representative genomes of most vertebrate classes and mammalian orders showed that nPE1 is a placental novelty. Using in silico paleogenomics we found that nPE1 originated from the exaptation of a mammalian-apparent LTR retrotransposon sometime between the metatherian/eutherian split (147 Mya) and the placental mammal radiation (≈ 90 Mya). Thus, the evolutionary origin of nPE1 differs, in kind and time, from that previously demonstrated for nPE2, which was exapted from a CORE-short interspersed nucleotide element (SINE) retroposon before the origin of prototherians, 166 Mya. Transgenic mice expressing the fluorescent markers tomato and EGFP driven by nPE1 or nPE2, respectively, demonstrated coexpression of both reporter genes along the entire arcuate nucleus. The onset of reporter gene expression guided by nPE1 and nPE2 was also identical and coincidental with the onset of Pomc expression in the presumptive mouse diencephalon. Thus, the independent exaptation of two unrelated retroposons into functional analogs regulating neuronal POMC expression constitutes an authentic example of convergent molecular evolution of cell-specific enhancers.