Human 5-lipoxygenase regulates transcription by association to euchromatin

Human 5-lipoxygenase (5-LO) is the key enzyme of leukotriene biosynthesis, mostly expressed in leukocytes and thus a crucial component of the innate immune system. In this study, we show that 5-LO, besides its canonical function as an arachidonic acid metabolizing enzyme, is a regulator of gene expression associated with euchromatin. By Crispr-Cas9-mediated 5-LO knockout (KO) in MonoMac6 (MM6) cells and subsequent RNA-Seq analysis, we identified 5-LO regulated genes which could be clustered to immune/defense response, cell adhesion, transcription and growth/developmental processes. Analysis of differentially expressed genes (DEG) identified cyclooxygenase-2 (COX2, PTGS2) and kynureninase (KYNU) as strongly regulated 5-LO target genes. 5-LO knockout affected MM6 cell adhesion and tryptophan metabolism via inhibition of the degradation of the immunoregulator kynurenine. By subsequent FAIRE-Seq and 5-LO ChIP-Seq analyses, we found an association of 5-LO with euchromatin, with prominent 5-LO binding to promoter regions in actively transcribed genes. By enrichment analysis of the ChIP-Seq results, we identified potential 5-LO interaction partners. Furthermore, 5-LO ChIP-Seq peaks resemble patterns of H3K27ac histone marks, suggesting that 5-LO recruitment mainly takes place at acetylated histones. In summary, we demonstrate a noncanonical function of 5-LO as transcriptional regulator in monocytic cells.

Selection following Gene Duplication Shapes Recent Genome Evolution in the Pea Aphid Acyrthosiphon pisum

Ecology of insects is as wide as their diversity, which reflects their high capacity of adaptation in most of the environments of our planet. Aphids, with over 4,000 species, have developed a series of adaptations including a high phenotypic plasticity and the ability to feed on the phloem sap of plants, which is enriched in sugars derived from photosynthesis. Recent analyses of aphid genomes have indicated a high level of shared ancestral gene duplications that might represent a basis for genetic innovation and broad adaptations. In addition, there are a large number of recent, species-specific gene duplications whose role in adaptation remains poorly understood. Here, we tested whether duplicates specific to the pea aphid Acyrthosiphon pisum are related to genomic innovation by combining comparative genomics, transcriptomics, and chromatin accessibility analyses. Consistent with large levels of neofunctionalization, we found that most of the recent pairs of gene duplicates evolved asymmetrically, showing divergent patterns of positive selection and gene expression. Genes under selection involved a plethora of biological functions, suggesting that neofunctionalization and tissue specificity, among other evolutionary mechanisms, have orchestrated the evolution of recent paralogs in the pea aphid and may have facilitated host-symbiont cooperation. Our comprehensive phylogenomics analysis allowed us to tackle the history of duplicated genes to pave the road toward understanding the role of gene duplication in ecological adaptation.

PeakMatcher facilitates updated Aedes aegypti embryonic cis-regulatory element map

BACKGROUND

The Aedes aegypti mosquito is a threat to human health across the globe. The A. aegypti genome was recently re-sequenced and re-assembled. Due to a combination of long-read PacBio and Hi-C sequencing, the AaegL5 assembly is chromosome complete and significantly improves the assembly in key areas such as the M/m sex-determining locus. Release of the updated genome assembly has precipitated the need to reprocess historical functional genomic data sets, including cis-regulatory element (CRE) maps that had previously been generated for A. aegypti.

RESULTS

We re-processed and re-analyzed the A. aegypti whole embryo FAIRE seq data to create an updated embryonic CRE map for the AaegL5 genome. We validated that the new CRE map recapitulates key features of the original AaegL3 CRE map. Further, we built on the improved assembly in the M/m locus to analyze overlaps of open chromatin regions with genes. To support the validation, we created a new method (PeakMatcher) for matching peaks from the same experimental data set across genome assemblies.

CONCLUSION

Use of PeakMatcher software, which is available publicly under an open-source license, facilitated the release of an updated and validated CRE map, which is available through the NIH GEO. These findings demonstrate that PeakMatcher software will be a useful resource for validation and transferring of previous annotations to updated genome assemblies.

Chromatin accessibility dynamics of Chlamydia-infected epithelial cells

Chlamydia are Gram-negative, obligate intracellular bacterial pathogens responsible for a broad spectrum of human and animal diseases. In humans, Chlamydia trachomatis is the most prevalent bacterial sexually transmitted infection worldwide and is the causative agent of trachoma (infectious blindness) in disadvantaged populations. Over the course of its developmental cycle, Chlamydia extensively remodels its intracellular niche and parasitises the host cell for nutrients, with substantial resulting changes to the host cell transcriptome and proteome. However, little information is available on the impact of chlamydial infection on the host cell epigenome and global gene regulation. Regions of open eukaryotic chromatin correspond to nucleosome-depleted regions, which in turn are associated with regulatory functions and transcription factor binding. We applied formaldehyde-assisted isolation of regulatory elements enrichment followed by sequencing (FAIRE-Seq) to generate temporal chromatin maps of C. trachomatis-infected human epithelial cells in vitro over the chlamydial developmental cycle. We detected both conserved and distinct temporal changes to genome-wide chromatin accessibility associated with C. trachomatis infection. The observed differentially accessible chromatin regions include temporally-enriched sets of transcription factors, which may help shape the host cell response to infection. These regions and motifs were linked to genomic features and genes associated with immune responses, re-direction of host cell nutrients, intracellular signalling, cell-cell adhesion, extracellular matrix, metabolism and apoptosis. This work provides another perspective to the complex response to chlamydial infection, and will inform further studies of transcriptional regulation and the epigenome in Chlamydia-infected human cells and tissues.

Pseudohyphal differentiation in Komagataella phaffii: investigating the FLO gene family

Many yeasts differentiate into multicellular phenotypes in adverse environmental conditions. Here, we investigate pseudohyphal growth in Komagataella phaffii and the involvement of the flocculin (FLO) gene family in its regulation. The K. phaffii FLO family consists of 13 members, and the conditions inducing pseudohyphal growth are different from Saccharomyces cerevisiae. So far, this phenotype was only observed when K. phaffii was cultivated at slow growth rates in glucose-limited chemostats, but not upon nitrogen starvation or the presence of fusel alcohols. Transcriptional analysis identified that FLO11, FLO400 and FLO5-1 are involved in the phenotype, all being controlled by the transcriptional regulator Flo8. The three genes exhibit a complex mechanism of expression and repression during transition from yeast to pseudohyphal form. Unlike in S. cerevisiae, deletion of FLO11 does not completely prevent the phenotype. In contrast, deletion of FLO400 or FLO5-1 prevents pseudohyphae formation, and hampers FLO11 expression. FAIRE-Seq data shows that the expression and repression of FLO400 and FLO5-1 are correlated to open or closed chromatin regions upstream of these genes, respectively. Our findings indicate that K. phaffii Flo400 and/or Flo5-1 act as upstream signals that lead to the induction of FLO11 upon glucose limitation in chemostats at slow growth and chromatin modulation is involved in the regulation of their expression.

Punctuated chromatin states regulate Plasmodium falciparum antigenic variation at the intron and 2 kb upstream regions

Background

Understanding the regulation mechanism of var gene expression is crucial for explaining antigenic variation in Plasmodium falciparum. Recent work observed that while all var genes produce transcripts, only a few var genes exhibit high expression levels. However, the global regulation of var expression and the relationship between epigenetic and genetic control remains to be established.

Result

We have systematically reanalyzed the existing genomic data including chromatin configurations and gene expressions; and for the first time used robust statistical methods to show that the intron and 2 kb upstream regions of each endogenous var gene always maintain high chromatin accessibility, with high potential to bind transcription factors (TFs). The levels of transcripts for different var gene family members are associated with this chromatin accessibility. Any given var gene thus shows punctuated chromatin states throughout the asexual life cycle. This is demonstrated by three independent transcript datasets. Chromatin accessibility in the var intron and 2 kb upstream regions are also positively correlated with their GC content, suggesting the level of var genes silencing might be encoded in their intron sequences. Interestingly, both var intron and 2 kb upstream regions exhibit higher chromatin accessibility when the genes have relatively lower transcription levels, suggesting a punctuated repressive function for these regulatory elements.

Conclusion

By integrating and analyzing epigenomic, genomic and transcriptomic data, our work reveals a novel distal element in var control. We found dynamic modulations of specific epigenetic marks around the var intron and distal upstream regions are involved in the general var gene expression patterns in malarial antigenic variation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3005-7) contains supplementary material, which is available to authorized users.

Identification of Lineage-Specific Cis-Regulatory Modules Associated with Variation in Transcription Factor Binding and Chromatin Activity Using Ornstein-Uhlenbeck Models

Scoring the impact of noncoding variation on the function of cis-regulatory regions, on their chromatin state, and on the qualitative and quantitative expression levels of target genes is a fundamental problem in evolutionary genomics. A particular challenge is how to model the divergence of quantitative traits and to identify relationships between the changes across the different levels of the genome, the chromatin activity landscape, and the transcriptome. Here, we examine the use of the Ornstein-Uhlenbeck (OU) model to infer selection at the level of predicted cis-regulatory modules (CRMs), and link these with changes in transcription factor binding and chromatin activity. Using publicly available cross-species ChIP-Seq and STARR-Seq data we show how OU can be applied genome-wide to identify candidate transcription factors for which binding site and CRM turnover is correlated with changes in regulatory activity. Next, we profile open chromatin in the developing eye across three Drosophila species. We identify the recognition motifs of the chromatin remodelers, Trithorax-like and Grainyhead as mostly correlating with species-specific changes in open chromatin. In conclusion, we show in this study that CRM scores can be used as quantitative traits and that motif discovery approaches can be extended towards more complex models of divergence.