Progress in genome sequencing is proceeding at an exponential pace, and several new algal genomes are becoming available every year. One of the challenges facing the community is the association of protein sequences encoded in the genomes with biological function. While most genome assembly projects generate annotations for predicted protein sequences, they are usually limited and integrate functional terms from a limited number of databases. Another challenge is the use of annotations to interpret large lists of 'interesting' genes generated by genome-scale datasets. Previously, these gene lists had to be analyzed across several independent biological databases, often on a gene-by-gene basis. In contrast, several annotation databases, such as DAVID, integrate data from multiple functional databases and reveal underlying biological themes of large gene lists. While several such databases have been constructed for animals, none is currently available for the study of algae. Due to renewed interest in algae as potential sources of biofuels and the emergence of multiple algal genome sequences, a significant need has arisen for such a database to process the growing compendiums of algal genomic data.

Of the ∼ 1.3 million Alu elements in the human genome, only a tiny number are estimated to be active in transcription by RNA polymerase (Pol) III. Tracing the individual loci from which Alu transcripts originate is complicated by their highly repetitive nature. By exploiting RNA-Seq data sets and unique Alu DNA sequences, we devised a bioinformatic pipeline allowing us to identify Pol III-dependent transcripts of individual Alu elements. When applied to ENCODE transcriptomes of seven human cell lines, this search strategy identified ∼ 1300 Alu loci corresponding to detectable transcripts, with ∼ 120 of them expressed in at least three cell lines. In vitro transcription of selected Alus did not reflect their in vivo expression properties, and required the native 5'-flanking region in addition to internal promoter. We also identified a cluster of expressed AluYa5-derived transcription units, juxtaposed to snaR genes on chromosome 19, formed by a promoter-containing left monomer fused to an Alu-unrelated downstream moiety. Autonomous Pol III transcription was also revealed for Alus nested within Pol II-transcribed genes. The ability to investigate Alu transcriptomes at single-locus resolution will facilitate both the identification of novel biologically relevant Alu RNAs and the assessment of Alu expression alteration under pathological conditions.

Here we demonstrate that heat shock protein 90 (HSP90) interacts with calpain-1, but not with calpain-2, and forms a discrete complex in which the protease maintains its catalytic activity, although with a lower affinity for Ca2+. Equilibrium gel distribution experiments show that this complex is composed by an equal number of molecules of each protein partner. Moreover, in resting cells, cytosolic calpain-1 is completely associated with HSP90. Since calpain-1, in association with HSP90, retains its proteolytic activity, and the chaperone is displaced by calpastatin also in the absence of Ca2+, the catalytic cleft of the protease is not involved in this association. Thus, calpain-1 can form two distinct complexes depending on the availability of calpastatin in the cytosol. The occurrence of a complex between HSP90 and calpain-1, in which the protease is still activable, can prevent the complete inhibition of the protease even in the presence of high calpastatin levels. We also demonstrate that in basal cell conditions HSP90 and calpain-1, but not calpain-2, are inserted in the multi-protein N-Methyl-D-Aspartate receptor (NMDAR) complex. The amount of calpain-1 at the NMDAR cluster is not modified in conditions of increased [Ca2+]i, and this resident protease is involved in the processing of NMDAR components. Finally, the amount of calpain-1 associated with NMDAR cluster is independent from Ca2+-mediated translocation. Our findings show that HSP90 plays an important role in maintaining a given and proper amount of calpain-1 at the functional sites.

Scl/Tal1 confers hemogenic competence and prevents ectopic cardiomyogenesis in embryonic endothelium by unknown mechanisms. We discovered that Scl binds to hematopoietic and cardiac enhancers that become epigenetically primed in multipotent cardiovascular mesoderm, to regulate the divergence of hematopoietic and cardiac lineages. Scl does not act as a pioneer factor but rather exploits a pre-established epigenetic landscape. As the blood lineage emerges, Scl binding and active epigenetic modifications are sustained in hematopoietic enhancers, whereas cardiac enhancers are decommissioned by removal of active epigenetic marks. Our data suggest that, rather than recruiting corepressors to enhancers, Scl prevents ectopic cardiogenesis by occupying enhancers that cardiac factors, such as Gata4 and Hand1, use for gene activation. Although hematopoietic Gata factors bind with Scl to both activated and repressed genes, they are dispensable for cardiac repression, but necessary for activating genes that enable hematopoietic stem/progenitor cell development. These results suggest that a unique subset of enhancers in lineage-specific genes that are accessible for regulators of opposing fates during the time of the fate decision provide a platform where the divergence of mutually exclusive fates is orchestrated.

Microalgae hold promise for yielding a biofuel feedstock that is sustainable, carbon-neutral, distributed, and only minimally disruptive for the production of food and feed by traditional agriculture. Amongst oleaginous eukaryotic algae, the B race of Botryococcus braunii is unique in that it produces large amounts of liquid hydrocarbons of terpenoid origin. These are comparable to fossil crude oil, and are sequestered outside the cells in a communal extracellular polymeric matrix material. Biosynthetic engineering of terpenoid bio-crude production requires identification of genes and reconstruction of metabolic pathways responsible for production of both hydrocarbons and other metabolites of the alga that compete for photosynthetic carbon and energy.

Mechanisms by which regulatory T (Treg) cells fail to control inflammation in asthma remain poorly understood. We show that a severe asthma-associated polymorphism in the gene encoding the interleukin (IL)-4 receptor alpha chain (Il4ra(R576)) promotes conversion of induced Treg (iTreg) cells toward a T helper 17 (TH17) cell fate. This skewing is mediated by the recruitment by IL-4Rα(R576) of the growth-factor-receptor-bound protein 2 (GRB2) adaptor protein, which drives IL-17 expression by activating a pathway that involves extracellular-signal-regulated kinase, IL-6 and the transcription factor STAT3. Treg cell-specific deletion of genes that regulate TH17 cell differentiation, including Il6ra and RAR-related orphan receptor gamma (Rorc), but not of Il4 or Il13, prevented exacerbated airway inflammation in mice expressing Il4ra(R576) (hereafter referred to as Il4ra(R576) mice). Furthermore, treatment of Il4ra(R576) mice with a neutralizing IL-6-specific antibody prevented iTreg cell reprogramming into TH17-like cells and protected against severe airway inflammation. These findings identify a previously unknown mechanism for the development of mixed TH2-TH17 cell inflammation in genetically prone individuals and point to interventions that stabilize iTreg cells as potentially effective therapeutic strategies.

Upon fertilization Chlamydomonas reinhardtii zygotes undergo a program of differentiation into a diploid zygospore that is accompanied by transcription of hundreds of zygote-specific genes. We identified a distinct sequence motif we term a zygotic response element (ZYRE) that is highly enriched in promoter regions of C reinhardtii early zygotic genes. A luciferase reporter assay was used to show that native ZYRE motifs within the promoter of zygotic gene ZYS3 or intron of zygotic gene DMT4 are necessary for zygotic induction. A synthetic luciferase reporter with a minimal promoter was used to show that ZYRE motifs introduced upstream are sufficient to confer zygotic upregulation, and that ZYRE-controlled zygotic transcription is dependent on the homeodomain transcription factor GSP1. We predict that ZYRE motifs will correspond to binding sites for the homeodomain proteins GSP1-GSM1 that heterodimerize and activate zygotic gene expression in early zygotes.

Metastatic renal cell carcinoma (mRCC) is nearly incurable and accounts for most of the mortality associated with RCC. Von Hippel Lindau (VHL) is a tumour suppressor that is lost in the majority of clear cell RCC (ccRCC) cases. Its role in regulating hypoxia-inducible factors-1α (HIF-1α) and -2α (HIF-2α) is well-studied. Recent work has demonstrated that VHL knock down induces an epithelial-mesenchymal transition (EMT) phenotype. In this study we showed that a CRISPR/Cas9-mediated knock out of VHL in the RENCA model leads to morphologic and molecular changes indicative of EMT, which in turn drives increased metastasis to the lungs. RENCA cells deficient in HIF-1α failed to undergo EMT changes upon VHL knockout. RNA-seq revealed several HIF-1α-regulated genes that are upregulated in our VHL knockout cells and whose overexpression signifies an aggressive form of ccRCC in the cancer genome atlas (TCGA) database. Independent validation in a new clinical dataset confirms the upregulation of these genes in ccRCC samples compared to adjacent normal tissue. Our findings indicate that loss of VHL could be driving tumour cell dissemination through stabilization of HIF-1α in RCC. A better understanding of the mechanisms involved in this phenomenon can guide the search for more effective treatments to combat mRCC.

The genetic origins of chemotherapy resistance are well established; however, the role of epigenetics in drug resistance is less well understood. To investigate mechanisms of drug resistance, we performed systematic genetic, epigenetic, and transcriptomic analyses of an alkylating agent-sensitive murine lymphoma cell line and a series of resistant lines derived by drug dose escalation.

The HIV-1 Trans-Activator of Transcription (Tat) protein binds to multiple host cellular factors and greatly enhances the level of transcription of the HIV genome. While Tat's control of viral transcription is well-studied, much less is known about the interaction of Tat with the human genome. Here, we report the genome-wide binding map of Tat to the human genome in Jurkat T cells using chromatin immunoprecipitation combined with next-generation sequencing. Surprisingly, we found that ~53% of the Tat target regions are within DNA repeat elements, greater than half of which are Alu sequences. The remaining target regions are located in introns and distal intergenic regions; only ~7% of Tat-bound regions are near transcription start sites (TSS) at gene promoters. Interestingly, Tat binds to promoters of genes that, in Jurkat cells, are bound by the ETS1 transcription factor, the CBP histone acetyltransferase and/or are enriched for histone H3 lysine 4 tri-methylation (H3K4me3) and H3K27me3. Tat binding is associated with genes enriched with functions in T cell biology and immune response. Our data reveal that Tat's interaction with the host genome is more extensive than previously thought, with potentially important implications for the viral life cycle.

Adenovirus small e1a oncoprotein causes ~70% reduction in cellular levels of histone H3 lysine 18 acetylation (H3K18ac). It is unclear, however, where this dramatic reduction occurs genome-wide. ChIP-sequencing revealed that by 24 h after expression, e1a erases 95% of H3K18ac peaks in normal, contact-inhibited fibroblasts and replaces them with one-third as many at new genomic locations. The H3K18ac peaks at promoters and intergenic regions of genes with fibroblast-related functions are eliminated after infection, and new H3K18ac peaks are established at promoters of highly induced genes that regulate cell cycling and at new putative enhancers. Strikingly, the regions bound by the retinoblastoma family of proteins in contact-inhibited fibroblasts gain new peaks of H3K18ac in the e1a-expressing cells, including 55% of RB1-bound loci. In contrast, over half of H3K9ac peaks are similarly distributed before and after infection, independently of RB1. The strategic redistribution of H3K18ac by e1a highlights the importance of this modification for transcriptional activation and cellular transformation as well as functional differences between the RB-family member proteins.

Investigation of the human microbiome has revealed diverse and complex microbial communities at distinct anatomic sites. The microbiome of the human sebaceous follicle provides a tractable model in which to study its dominant bacterial inhabitant, Propionibacterium acnes, which is thought to contribute to the pathogenesis of the human disease acne. To explore the diversity of the bacteriophages that infect P. acnes, 11 P. acnes phages were isolated from the sebaceous follicles of donors with healthy skin or acne and their genomes were sequenced. Comparative genomic analysis of the P. acnes phage population, which spans a 30-year temporal period and a broad geographic range, reveals striking similarity in terms of genome length, percent GC content, nucleotide identity (>85%), and gene content. This was unexpected, given the far-ranging diversity observed in virtually all other phage populations. Although the P. acnes phages display a broad host range against clinical isolates of P. acnes, two bacterial isolates were resistant to many of these phages. Moreover, the patterns of phage resistance correlate closely with the presence of clustered regularly interspaced short palindromic repeat elements in the bacteria that target a specific subset of phages, conferring a system of prokaryotic innate immunity. The limited diversity of the P. acnes bacteriophages, which may relate to the unique evolutionary constraints imposed by the lipid-rich anaerobic environment in which their bacterial hosts reside, points to the potential utility of phage-based antimicrobial therapy for acne.

Many of the functional units in cells are multi-protein complexes such as RNA polymerase, the ribosome, and the proteasome. For such units to work together, one might expect a high level of regulation to enable co-appearance or repression of sets of complexes at the required time. However, this type of coordinated regulation between whole complexes is difficult to detect by existing methods for analyzing mRNA co-expression. We propose a new methodology that is able to detect such higher order relationships.

Phylogenetic profiles record the occurrence of homologs of genes across fully sequenced organisms. Proteins with similar profiles are typically components of protein complexes or metabolic pathways. Various existing methods measure similarity between two profiles and, hence, the likelihood that the two proteins co-evolve. Some methods ignore phylogenetic relationships between organisms while others account for such with metrics that explicitly model the likelihood of two proteins co-evolving on a tree. The latter methods more sensitively detect co-evolving proteins, but at a significant computational cost. Here we propose a novel heuristic to improve phylogenetic profile analysis that accounts for phylogenetic relationships between genomes in a computationally efficient fashion. We first order the genomes within profiles and then enumerate runs of consecutive matches and accurately compute the probability of observing these. We hypothesize that profiles with many runs are more likely to involve functionally related proteins than profiles in which all the matches are concentrated in one interval of the tree.

Interpretation of lists of genes or proteins with altered expression is a critical and time-consuming part of microarray and proteomics research, but relatively little attention has been paid to methods for extracting biological meaning from these output lists. One powerful approach is to examine the expression of predefined biological pathways and gene sets, such as metabolic and signaling pathways and macromolecular complexes. Although many methods for measuring pathway expression have been proposed, a systematic analysis of the performance of multiple methods over multiple independent data sets has not previously been reported.

During mammalian embryogenesis, de novo hematopoiesis occurs transiently in multiple anatomical sites including the yolk sac, dorsal aorta, and heart tube. A long-unanswered question is whether these local transient hematopoietic mechanisms are essential for embryonic growth. Here, we show that endocardial hematopoiesis is critical for cardiac valve remodeling as a source of tissue macrophages. Colony formation assay from explanted heart tubes and genetic lineage tracing with the endocardial specific Nfatc1-Cre mouse revealed that hemogenic endocardium is a de novo source of tissue macrophages in the endocardial cushion, the primordium of the cardiac valves. Surface marker characterization, gene expression profiling, and ex vivo phagocytosis assay revealed that the endocardially derived cardiac tissue macrophages play a phagocytic and antigen presenting role. Indeed, genetic ablation of endocardially derived macrophages caused severe valve malformation. Together, these data suggest that transient hemogenic activity in the endocardium is indispensable for the valvular tissue remodeling in the heart.

How artificial environmental cues are biologically integrated and transgenerationally inherited is still poorly understood. Here, we investigate the mechanisms of inheritance of reproductive outcomes elicited by the model environmental chemical Bisphenol A in C. elegans. We show that Bisphenol A (BPA) exposure causes the derepression of an epigenomically silenced transgene in the germline for 5 generations, regardless of ancestral response. Chromatin immunoprecipitation sequencing (ChIP-seq), histone modification quantitation, and immunofluorescence assays revealed that this effect is associated with a reduction of the repressive marks H3K9me3 and H3K27me3 in whole worms and in germline nuclei in the F3, as well as with reproductive dysfunctions, including germline apoptosis and embryonic lethality. Furthermore, targeting of the Jumonji demethylases JMJD-2 and JMJD-3/UTX-1 restores H3K9me3 and H3K27me3 levels, respectively, and it fully alleviates the BPA-induced transgenerational effects. Together, our results demonstrate the central role of repressive histone modifications in the inheritance of reproductive defects elicited by a common environmental chemical exposure.

Reversal reactions (RRs) in leprosy are characterized by a reduction in the number of bacilli in lesions associated with an increase in cell-mediated immunity against the intracellular bacterium Mycobacterium leprae, the causative pathogen of leprosy. To identify the mechanisms that contribute to cell-mediated immunity in leprosy, we measured changes in the whole blood-derived transcriptome of patients with leprosy before, during and after RR. We identified an 'RR signature' of 1017 genes that were upregulated at the time of the clinical diagnosis of RR. Using weighted gene correlated network analysis (WGCNA), we detected a module of 794 genes, bisque4, that was significantly correlated with RR, of which 434 genes were part of the RR signature. An enrichment for both IFN-γ and IFN-β downstream gene pathways was present in the RR signature as well as the RR upregulated genes in the bisque4 module, including those encoding proteins of the guanylate binding protein (GBP) family that contributes to antimicrobial responses against mycobacteria. Specifically, GBP1, GBP2, GBP3 and GBP5 mRNAs were upregulated in the RR peripheral blood transcriptome, with GBP1, GBP2 and GBP5 mRNAs also upregulated in the RR disease lesion transcriptome. These data indicate that RRs involve a systemic upregulation of IFN-γ downstream genes including GBP family members as part of the host antimicrobial response against mycobacteria.

The conserved hypoxia inducible factor 1 α (HIF1α) injury-response pro-survival pathway has recently been implicated in early beta cell dysfunction but slow beta cell loss in type 2 diabetes. We hypothesised that the unexplained prolonged prediabetes phase in type 1 diabetes may also be, in part, due to activation of the HIF1α signalling pathway.

We profiled soybean and Arabidopsis methylomes from the globular stage through dormancy and germination to understand the role of methylation in seed formation. CHH methylation increases significantly during development throughout the entire seed, targets primarily transposable elements (TEs), is maintained during endoreduplication, and drops precipitously within the germinating seedling. By contrast, no significant global changes in CG- and CHG-context methylation occur during the same developmental period. An Arabidopsis ddcc mutant lacking CHH and CHG methylation does not affect seed development, germination, or major patterns of gene expression, implying that CHH and CHG methylation does not play a significant role in seed development or in regulating seed gene activity. By contrast, over 100 TEs are transcriptionally de-repressed in ddcc seeds, suggesting that the increase in CHH-context methylation may be a failsafe mechanism to reinforce transposon silencing. Many genes encoding important classes of seed proteins, such as storage proteins, oil biosynthesis enzymes, and transcription factors, reside in genomic regions devoid of methylation at any stage of seed development. Many other genes in these classes have similar methylation patterns, whether the genes are active or repressed. Our results suggest that methylation does not play a significant role in regulating large numbers of genes important for programming seed development in both soybean and Arabidopsis. We conclude that understanding the mechanisms controlling seed development will require determining how cis-regulatory elements and their cognate transcription factors are organized in genetic regulatory networks.