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

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

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.

Activation of ß-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 the 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-wk-old C57BL/6NRj mice were treated with corticosterone via drinking water or placebo for 4 wk 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 being 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.NEW & NOTEWORTHY This preclinical study in mice shows that the ß3-adrenergic receptor can be a potential therapeutic approach to counteracting glucocorticoid (GC)-induced obesity and metabolic dysfunction. Both cold acclimation and ß3-adrenergic receptor stimulation in a mouse model of excess glucocorticoids were adequate in not only preventing obesity, adiposity, and adipose tissue dysfunction but also correcting hyperinsulinemia, hyperleptinemia, and dyslipidemia.

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 (PGE2), and cysteinyl leukotrienes, e.g., Leukotriene C4 (LTC4). 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 CD48hi 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.

Convergent and lineage-specific genomic differences in limb regulatory elements in limbless reptile lineages.

Loss of limbs evolved many times in squamate reptiles. Here we investigated the genomic basis of convergent limb loss in reptiles. We sequenced the genomes of a closely related pair of limbless-limbed gymnophthalmid lizards and performed a comparative genomic analysis including five snakes and the limbless glass lizard. Our analysis of these three independent limbless lineages revealed that signatures of shared sequence or transcription factor binding site divergence in individual limb regulatory elements are generally rare. Instead, shared divergence occurs more often at the level of signaling pathways, involving different regulatory elements associated with the same limb genes (such as Hand2 or Hox) and/or patterning mechanisms (such as Shh signaling). Interestingly, although snakes are known to have mutations in the Shh ZRS limb enhancer, this enhancer lacks relevant mutations in limbless lizards. Thus, different mechanisms could contribute to limb loss, and there are likely multiple evolutionary paths to limblessness in reptiles.

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.

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.

Large-scale genetic perturbations reveal regulatory networks and an abundance of gene-specific repressors.

To understand regulatory systems, it would be useful to uniformly determine how different components contribute to the expression of all other genes. We therefore monitored mRNA expression genome-wide, for individual deletions of one-quarter of yeast genes, focusing on (putative) regulators. The resulting genetic perturbation signatures reflect many different properties. These include the architecture of protein complexes and pathways, identification of expression changes compatible with viability, and the varying responsiveness to genetic perturbation. The data are assembled into a genetic perturbation network that shows different connectivities for different classes of regulators. Four feed-forward loop (FFL) types are overrepresented, including incoherent type 2 FFLs that likely represent feedback. Systematic transcription factor classification shows a surprisingly high abundance of gene-specific repressors, suggesting that yeast chromatin is not as generally restrictive to transcription as is often assumed. The data set is useful for studying individual genes and for discovering properties of an entire regulatory system.

A central role for TFIID in the pluripotent transcription circuitry.

Embryonic stem (ES) cells are pluripotent and characterized by open chromatin and high transcription levels, achieved through auto-regulatory and feed-forward transcription factor loops. ES-cell identity is maintained by a core of factors including Oct4 (also known as Pou5f1), Sox2, Klf4, c-Myc (OSKM) and Nanog, and forced expression of the OSKM factors can reprogram somatic cells into induced pluripotent stem cells (iPSCs) resembling ES cells. These gene-specific factors for RNA-polymerase-II-mediated transcription recruit transcriptional cofactors and chromatin regulators that control access to and activity of the basal transcription machinery on gene promoters. How the basal transcription machinery is involved in setting and maintaining the pluripotent state is unclear. Here we show that knockdown of the transcription factor IID (TFIID) complex affects the pluripotent circuitry in mouse ES cells and inhibits reprogramming of fibroblasts. TFIID subunits and the OSKM factors form a feed-forward loop to induce and maintain a stable transcription state. Notably, transient expression of TFIID subunits greatly enhanced reprogramming. These results show that TFIID is critical for transcription-factor-mediated reprogramming. We anticipate that, by creating plasticity in gene expression programs, transcription complexes such as TFIID assist reprogramming into different cellular states.

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.

The specificity and topology of chromatin interaction pathways in yeast.

Packaging of DNA into chromatin has a profound impact on gene expression. To understand how changes in chromatin influence transcription, we analyzed 165 mutants of chromatin machinery components in Saccharomyces cerevisiae. mRNA expression patterns change in 80% of mutants, always with specific effects, even for loss of widespread histone marks. The data are assembled into a network of chromatin interaction pathways. The network is function based, has a branched, interconnected topology, and lacks strict one-to-one relationships between complexes. Chromatin pathways are not separate entities for different gene sets, but share many components. The study evaluates which interactions are important for which genes and predicts additional interactions, for example between Paf1C and Set3C, as well as a role for Mediator in subtelomeric silencing. The results indicate the presence of gene-dependent effects that go beyond context-dependent binding of chromatin factors and provide a framework for understanding how specificity is achieved through regulating chromatin.

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.

Stratification of pediatric ALL by in vitro cellular responses to DNA double-strand breaks provides insight into the molecular mechanisms underlying clinical response.

The molecular basis of different outcomes in pediatric acute lymphoblastic leukemia (ALL) remains poorly understood. We addressed the clinical significance and mechanisms behind in vitro cellular responses to ionizing radiation (IR)-induced DNA double-strand breaks in 74 pediatric patients with ALL. We found an apoptosis-resistant response in 36% of patients characterized by failure to cleave caspase-3, -7, -9, and PARP1 by 24 hours after IR and an apoptosis-sensitive response with the cleavage of the same substrates in the remaining 64% of leukemias. Resistance to IR in vitro was associated with poor early blast clearance at day 7 or 15 and persistent minimal residual disease (MRD) at day 28 of induction treatment. Global gene expression profiling revealed abnormal up-regulation of multiple prosurvival pathways in response to IR in apoptosis-resistant leukemias and differential posttranscriptional activation of the PI3-Akt pathway was observed in representative resistant cases. Importantly, pharmacologic inhibition of selected prosurvival pathways sensitized apoptosis-resistant ALL cells to IR in vitro. We suggest that abnormal prosurvival responses to DNA damage provide one of the mechanisms of primary resistance in ALL, and that they should be considered as therapeutic targets in children with aggressive disease.

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.