Cohesin SMC1β promotes closed chromatin and controls TERRA expression at spermatocyte telomeres

Previous data showed that meiotic cohesin SMC1β protects spermatocyte telomeres from damage. The underlying reason, however, remained unknown as the expressions of telomerase and shelterin components were normal in Smc1β -/- spermatocytes. Here. we report that SMC1β restricts expression of the long noncoding RNA TERRA (telomeric repeat containing RNA) in spermatocytes. In somatic cell lines increased TERRA was reported to cause telomere damage through altering telomere chromatin structure. In Smc1β -/- spermatocytes, we observed strongly increased levels of TERRA which accumulate on damaged chromosomal ends, where enhanced R-loop formation was found. This suggested a more open chromatin configuration near telomeres in Smc1β -/- spermatocytes, which was confirmed by ATAC-seq. Telomere-distal regions were not affected by the absence of SMC1β but RNA-seq revealed increased transcriptional activity in telomere-proximal regions. Thus, SMC1β promotes closed chromatin specifically near telomeres and limits TERRA expression in spermatocytes.

Corrigendum: The aryl hydrocarbon receptor regulates lipid mediator production in alveolar macrophages

[This corrects the article DOI: 10.3389/fimmu.2023.1157373.].

Activation of beta-adrenergic receptor signaling prevents glucocorticoid-induced obesity and adipose tissue dysfunction in male mice

Elevated serum concentrations of glucocorticoids (GCs) result in excessive lipid accumulation in white adipose tissue (WAT) as well as dysfunction of thermogenic brown adipose tissue (BAT) - ultimately leading to the development of obesity and metabolic disease. Here, we hypothesized that activation of the sympathetic nervous system either via cold exposure or use of a selective β3-adrenergic receptor (β3-AR) agonist alleviates the adverse metabolic effects of chronic GC exposure in rodents. To this end, male 10-week-old C57BL/6NRj mice were treated with corticosterone via drinking water or placebo for 4 weeks while being maintained at 29°C (thermoneutrality), 22°C (room temperature) or 13°C (cold temperature); in a follow-up study mice received a selective β3-AR agonist or placebo with and without corticosterone while maintained at room temperature. Body weight and food intake were monitored throughout the study. Histological and molecular analyses were performed on white and brown adipose depots. Cold exposure not only preserved the thermogenic function of brown adipose tissue, but also reversed GC-induced lipid accumulation in white adipose tissue and corrected GC-driven obesity, hyperinsulinemia and hyperglycemia. The metabolic benefits of cold exposure were associated with enhanced sympathetic activity in adipose tissue, thus potentially linking an increase in sympathetic signaling to the observed metabolic benefits. In line with this concept, chronic administration of a selective β3-AR agonist reproduced the beneficial metabolic effects of cold adaption during exposure to exogenous GCs. This preclinical study demonstrates the potential of β3-AR as a therapeutic target in the management and prevention of GC-induced metabolic disease.

The aryl hydrocarbon receptor regulates lipid mediator production in alveolar macrophages

Allergic inflammation of the airways such as allergic asthma is a major health problem with growing incidence world-wide. One cardinal feature in severe type 2-dominated airway inflammation is the release of lipid mediators of the eicosanoid family that can either promote or dampen allergic inflammation. Macrophages are key producers of prostaglandins and leukotrienes which play diverse roles in allergic airway inflammation and thus require tight control. Using RNA- and ATAC-sequencing, liquid chromatography coupled to mass spectrometry (LC-MS/MS), enzyme immunoassays (EIA), gene expression analysis and in vivo models, we show that the aryl hydrocarbon receptor (AhR) contributes to this control via transcriptional regulation of lipid mediator synthesis enzymes in bone marrow-derived as well as in primary alveolar macrophages. In the absence or inhibition of AhR activity, multiple genes of both the prostaglandin and the leukotriene pathway were downregulated, resulting in lower synthesis of prostanoids, such as prostaglandin E2 (PGE₂), and cysteinyl leukotrienes, e.g., Leukotriene C4 (LTC₄). These AhR-dependent genes include PTGS1 encoding for the enzyme cyclooxygenase 1 (COX1) and ALOX5 encoding for the arachidonate 5-lipoxygenase (5-LO) both of which major upstream regulators of the prostanoid and leukotriene pathway, respectively. This regulation is independent of the activation stimulus and partially also detectable in unstimulated macrophages suggesting an important role of basal AhR activity for eicosanoid production in steady state macrophages. Lastly, we demonstrate that AhR deficiency in hematopoietic but not epithelial cells aggravates house dust mite induced allergic airway inflammation. These results suggest an essential role for AhR-dependent eicosanoid regulation in macrophages during homeostasis and inflammation.

Fate mapping of hematopoietic stem cells reveals two pathways of native thrombopoiesis

Hematopoietic stem cells (HSCs) produce highly diverse cell lineages. Here, we chart native lineage pathways emanating from HSCs and define their physiological regulation by computationally integrating experimental approaches for fate mapping, mitotic tracking, and single-cell RNA sequencing. We find that lineages begin to split when cells leave the tip HSC population, marked by high Sca-1 and CD201 expression. Downstream, HSCs either retain high Sca-1 expression and the ability to generate lymphocytes, or irreversibly reduce Sca-1 level and enter into erythro-myelopoiesis or thrombopoiesis. Thrombopoiesis is the sum of two pathways that make comparable contributions in steady state, a long route via multipotent progenitors and CD48^hi megakaryocyte progenitors (MkPs), and a short route from HSCs to developmentally distinct CD48^−/lo MkPs. Enhanced thrombopoietin signaling differentially accelerates the short pathway, enabling a rapid response to increasing demand. In sum, we provide a blueprint for mapping physiological differentiation fluxes from HSCs and decipher two functionally distinct pathways of native thrombopoiesis.

Hematopoietic stem cells produce diverse cell lineages. Here, the authors apply single-cell RNA-seq, computational integration of non-perturbative approaches for fate-mapping, and mitotic tracking to chart lineage decisions in native hematopoiesis and identify megakaryocyte progenitors that directly link HSCs to megakaryocytes.

The genome of the tegu lizard Salvator merianae: combining Illumina, PacBio, and optical mapping data to generate a highly contiguous assembly

Background

Reptiles are a species-rich group with great phenotypic and life history diversity but are highly underrepresented among the vertebrate species with sequenced genomes.

Results

Here, we report a high-quality genome assembly of the tegu lizard, Salvator merianae, the first lacertoid with a sequenced genome. We combined 74X Illumina short-read, 29.8X Pacific Biosciences long-read, and optical mapping data to generate a high-quality assembly with a scaffold N50 value of 55.4 Mb. The contig N50 value of this assembly is 521 Kb, making it the most contiguous reptile assembly so far. We show that the tegu assembly has the highest completeness of coding genes and conserved non-exonic elements (CNEs) compared to other reptiles. Furthermore, the tegu assembly has the highest number of evolutionarily conserved CNE pairs, corroborating a high assembly contiguity in intergenic regions. As in other reptiles, long interspersed nuclear elements comprise the most abundant transposon class. We used transcriptomic data, homology- and de novo gene predictions to annotate 22,413 coding genes, of which 16,995 (76%) likely have human orthologs as inferred by CESAR-derived gene mappings. Finally, we generated a multiple genome alignment comprising 10 squamates and 7 other amniote species and identified conserved regions that are under evolutionary constraint. CNEs cover 38 Mb (1.8%) of the tegu genome, with 3.3 Mb in these elements being squamate specific. In contrast to placental mammal-specific CNEs, very few of these squamate-specific CNEs (<20 Kb) overlap transposons, highlighting a difference in how lineage-specific CNEs originated in these two clades.

Conclusions

The tegu lizard genome together with the multiple genome alignment and comprehensive conserved element datasets provide a valuable resource for comparative genomic studies of reptiles and other amniotes.

Phenotype loss is associated with widespread divergence of the gene regulatory landscape in evolution

Detecting the genomic changes underlying phenotypic changes between species is a main goal of evolutionary biology and genomics. Evolutionary theory predicts that changes in cis-regulatory elements are important for morphological changes. We combined genome sequencing, functional genomics and genome-wide comparative analyses to investigate regulatory elements in lineages that lost morphological traits. We first show that limb loss in snakes is associated with widespread divergence of limb regulatory elements. We next show that eye degeneration in subterranean mammals is associated with widespread divergence of eye regulatory elements. In both cases, sequence divergence results in an extensive loss of transcription factor binding sites. Importantly, diverged regulatory elements are associated with genes required for normal limb patterning or normal eye development and function, suggesting that regulatory divergence contributed to the loss of these phenotypes. Together, our results show that genome-wide decay of the phenotype-specific cis-regulatory landscape is a hallmark of lost morphological traits.

Cis-regulatory elements are important factors for morphological changes. Here, the authors show widespread divergence of limb and eye regulatory elements in limb loss in snakes and eye degeneration in subterranean mammals respectively.

Iterative error correction of long sequencing reads maximizes accuracy and improves contig assembly

Next-generation sequencers such as Illumina can now produce reads up to 300 bp with high throughput, which is attractive for genome assembly. A first step in genome assembly is to computationally correct sequencing errors. However, correcting all errors in these longer reads is challenging. Here, we show that reads with remaining errors after correction often overlap repeats, where short erroneous k-mers occur in other copies of the repeat. We developed an iterative error correction pipeline that runs the previously published String Graph Assembler (SGA) in multiple rounds of k-mer-based correction with an increasing k-mer size, followed by a final round of overlap-based correction. By combining the advantages of small and large k-mers, this approach corrects more errors in repeats and minimizes the total amount of erroneous reads. We show that higher read accuracy increases contig lengths two to three times. We provide SGA-Iteratively Correcting Errors (https://github.com/hillerlab/IterativeErrorCorrection/) that implements iterative error correction by using modules from SGA.

A high-resolution gene expression atlas of epistasis between gene-specific transcription factors exposes potential mechanisms for genetic interactions

BACKGROUND

Genetic interactions, or non-additive effects between genes, play a crucial role in many cellular processes and disease. Which mechanisms underlie these genetic interactions has hardly been characterized. Understanding the molecular basis of genetic interactions is crucial in deciphering pathway organization and understanding the relationship between genotype, phenotype and disease.

RESULTS

To investigate the nature of genetic interactions between gene-specific transcription factors (GSTFs) in Saccharomyces cerevisiae, we systematically analyzed 72 GSTF pairs by gene expression profiling double and single deletion mutants. These pairs were selected through previously published growth-based genetic interactions as well as through similarity in DNA binding properties. The result is a high-resolution atlas of gene expression-based genetic interactions that provides systems-level insight into GSTF epistasis. The atlas confirms known genetic interactions and exposes new ones. Importantly, the data can be used to investigate mechanisms that underlie individual genetic interactions. Two molecular mechanisms are proposed, "buffering by induced dependency" and "alleviation by derepression".

CONCLUSIONS

These mechanisms indicate how negative genetic interactions can occur between seemingly unrelated parallel pathways and how positive genetic interactions can indirectly expose parallel rather than same-pathway relationships. The focus on GSTFs is important for understanding the transcription regulatory network of yeast as it uncovers details behind many redundancy relationships, some of which are completely new. In addition, the study provides general insight into the complex nature of epistasis and proposes mechanistic models for genetic interactions, the majority of which do not fall into easily recognizable within- or between-pathway relationships.

Yeast glucose pathways converge on the transcriptional regulation of trehalose biosynthesis

Background

Cellular glucose availability is crucial for the functioning of most biological processes. Our understanding of the glucose regulatory system has been greatly advanced by studying the model organism Saccharomyces cerevisiae, but many aspects of this system remain elusive. To understand the organisation of the glucose regulatory system, we analysed 91 deletion mutants of the different glucose signalling and metabolic pathways in Saccharomyces cerevisiae using DNA microarrays.

Results

In general, the mutations do not induce pathway-specific transcriptional responses. Instead, one main transcriptional response is discerned, which varies in direction to mimic either a high or a low glucose response. Detailed analysis uncovers established and new relationships within and between individual pathways and their members. In contrast to signalling components, metabolic components of the glucose regulatory system are transcriptionally more frequently affected. A new network approach is applied that exposes the hierarchical organisation of the glucose regulatory system.

Conclusions

The tight interconnection between the different pathways of the glucose regulatory system is reflected by the main transcriptional response observed. Tps2 and Tsl1, two enzymes involved in the biosynthesis of the storage carbohydrate trehalose, are predicted to be the most downstream transcriptional components. Epistasis analysis of tps2Δ double mutants supports this prediction. Although based on transcriptional changes only, these results suggest that all changes in perceived glucose levels ultimately lead to a shift in trehalose biosynthesis.

Functional overlap and regulatory links shape genetic interactions between signaling pathways

To understand relationships between phosphorylation-based signaling pathways, we analyzed 150 deletion mutants of protein kinases and phosphatases in S. cerevisiae using DNA microarrays. Downstream changes in gene expression were treated as a phenotypic readout. Double mutants with synthetic genetic interactions were included to investigate genetic buffering relationships such as redundancy. Three types of genetic buffering relationships are identified: mixed epistasis, complete redundancy, and quantitative redundancy. In mixed epistasis, the most common buffering relationship, different gene sets respond in different epistatic ways. Mixed epistasis arises from pairs of regulators that have only partial overlap in function and that are coupled by additional regulatory links such as repression of one by the other. Such regulatory modules confer the ability to control different combinations of processes depending on condition or context. These properties likely contribute to the evolutionary maintenance of paralogs and indicate a way in which signaling pathways connect for multiprocess control.

Functional modules integrating essential cellular functions are predictive of the response of leukaemia cells to DNA damage

MOTIVATION

Childhood B-precursor lymphoblastic leukaemia (ALL) is the most common paediatric malignancy. Despite the fact that 80% of ALL patients respond to anti-cancer drugs, the patho-physiology of this disease is still not fully understood. mRNA expression-profiling studies that have been performed have not yet provided novel insights into the mechanisms behind cellular response to DNA damage. More powerful data analysis techniques may be required for identifying novel functional pathways involved in the cellular responses to DNA damage.

RESULTS

In order to explore the possibility that unforeseen biological processes may be involved in the response to DNA damage, we have developed and applied a novel procedure for the identification of functional modules in ALL cells. We have discovered that the overall activity of functional modules integrating protein degradation and mRNA processing is predictive of response to DNA damage.

AVAILABILITY

Supplementary material including R code, additional results, experimental datasets, as well as a detailed description of the methodology are available at http://www.bip.bham.ac.uk/vivo/fumo.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Models and computational strategies linking physiological response to molecular networks from large-scale data

An important area of research in systems biology involves the analysis and integration of genome-wide functional datasets. In this context, a major goal is the identification of a putative molecular network controlling physiological response from experimental data. With very fragmentary mechanistic information, this is a challenging task. A number of methods have been developed, each one with the potential to address an aspect of the problem. Here, we review some of the most widely used methodologies and report new results in support of the usefulness of modularization and other modelling techniques in identifying components of the molecular networks that are predictive of physiological response. We also discuss how system identification in biology could be approached, using a combination of methodologies that aim to reconstruct the relationship between molecular pathways and physiology at different levels of the organizational complexity of the molecular network.

Fragrep: an efficient search tool for fragmented patterns in genomic sequences

Many classes of non-coding RNAs (ncRNAs; including Y RNAs, vault RNAs, RNase P RNAs, and MRP RNAs, as well as a novel class recently discovered in Dictyostelium discoideum) can be characterized by a pattern of short but well-conserved sequence elements that are separated by poorly conserved regions of sometimes highly variable lengths. Local alignment algorithms such as BLAST are therefore ill-suited for the discovery of new homologs of such ncRNAs in genomic sequences. The Fragrep tool instead implements an efficient algorithm for detecting the pattern fragments that occur in a given order. For each pattern fragment, the mismatch tolerance and bounds on the length of the intervening sequences can be specified separately. Furthermore, matches can be ranked by a statistically well-motivated scoring scheme.