Ensembl 2020

Evolutionary screening of precision oncology biomarkers and its applications in prognostic model construction

Biomarker screening is critical for precision oncology. However, one of the main challenges in precision oncology is that the screened biomarkers often fail to achieve the expected clinical effects and are rarely approved by regulatory authorities. Considering the close association between cancer pathogenesis and the evolutionary events of organisms, we first explored the evolutionary feature underlying clinically approved biomarkers, and two evolutionary features of approved biomarkers (Ohnologs and specific evolutionary stages of genes) were identified. Subsequently, we utilized evolutionary features for screening potential prognostic biomarkers in four common cancers: head and neck squamous cell carcinoma, liver hepatocellular carcinoma, lung adenocarcinoma, and lung squamous cell carcinoma. Finally, we constructed an evolution-strengthened prognostic model (ESPM) for cancers. These models can predict cancer patients' survival time across different cancer cohorts effectively and perform better than conventional models. In summary, our study highlights the application potentials of evolutionary information in precision oncology biomarker screening.

The molecular evolution of cancer associated genes in mammals

Cancer is a disease that many multicellular organisms have faced for millions of years, and species have evolved various tumour suppression mechanisms to control oncogenesis. Although cancer occurs across the tree of life, cancer related mortality risks vary across mammalian orders, with Carnivorans particularly affected. Evolutionary theory predicts different selection pressures on genes associated with cancer progression and suppression, including oncogenes, tumour suppressor genes and immune genes. Therefore, we investigated the evolutionary history of cancer associated gene sequences across 384 mammalian taxa, to detect signatures of selection across categories of oncogenes (GRB2, FGL2 and CDC42), tumour suppressors (LITAF, Casp8 and BRCA2) and immune genes (IL2, CD274 and B2M). This approach allowed us to conduct a fine scale analysis of gene wide and site-specific signatures of selection across mammalian lineages under the lens of cancer susceptibility. Phylogenetic analyses revealed that for most species the evolution of cancer associated genes follows the species' evolution. The gene wide selection analyses revealed oncogenes being the most conserved, tumour suppressor and immune genes having similar amounts of episodic diversifying selection. Despite BRCA2's status as a key caretaker gene, episodic diversifying selection was detected across mammals. The site-specific selection analyses revealed that the two apoptosis associated domains of the Casp8 gene of bats (Chiroptera) are under opposing forces of selection (positive and negative respectively), highlighting the importance of site-specific selection analyses to understand the evolution of highly complex gene families. Our results highlighted the need to critically assess different types of selection pressure on cancer associated genes when investigating evolutionary adaptations to cancer across the tree of life. This study provides an extensive assessment of cancer associated genes in mammals with highly representative, and substantially large sample size for a comparative genomic analysis in the field and identifies various avenues for future research into the mechanisms of cancer resistance and susceptibility in mammals.

Resident and recruited macrophages differentially contribute to cardiac healing after myocardial ischemia

Cardiac macrophages are heterogenous in phenotype and functions, which has been associated with differences in their ontogeny. Despite extensive research, our understanding of the precise role of different subsets of macrophages in ischemia/reperfusion (I/R) injury remains incomplete. We here investigated macrophage lineages and ablated tissue macrophages in homeostasis and after I/R injury in a CSF1R-dependent manner. Genomic deletion of a fms-intronic regulatory element (FIRE) in the Csf1r locus resulted in specific absence of resident homeostatic and antigen-presenting macrophages, without affecting the recruitment of monocyte-derived macrophages to the infarcted heart. Specific absence of homeostatic, monocyte-independent macrophages altered the immune cell crosstalk in response to injury and induced proinflammatory neutrophil polarization, resulting in impaired cardiac remodeling without influencing infarct size. In contrast, continuous CSF1R inhibition led to depletion of both resident and recruited macrophage populations. This augmented adverse remodeling after I/R and led to an increased infarct size and deterioration of cardiac function. In summary, resident macrophages orchestrate inflammatory responses improving cardiac remodeling, while recruited macrophages determine infarct size after I/R injury. These findings attribute distinct beneficial effects to different macrophage populations in the context of myocardial infarction.

Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3' untranslated region (3'-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3'-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3'-UTR eQTLs in immune-related genes. Our approach identifies numerous 3'-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

Conserved regulatory switches for the transition from natal down to juvenile feather in birds

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We report that extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial-mesenchymal interactions for branching morphogenesis. α-SMA (ACTA2) compartmentalizes dermal papilla stem cells for feather renewal cycling. LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We show that this primary feather transition is largely conserved in chicken (precocial) and zebra finch (altricial) and discuss the possibility that this evolutionary adaptation process started in feathered dinosaurs.

Systematic mapping of TF-mediated cell fate changes by a pooled induction coupled with scRNA-seq and multi-omics approaches

Transcriptional regulation controls cellular functions through interactions between transcription factors (TFs) and their chromosomal targets. However, understanding the fate conversion potential of multiple TFs in an inducible manner remains limited. Here, we introduce iTF-seq as a method for identifying individual TFs that can alter cell fate toward specific lineages at a single-cell level. iTF-seq enables time course monitoring of transcriptome changes, and with biotinylated individual TFs, it provides a multi-omics approach to understanding the mechanisms behind TF-mediated cell fate changes. Our iTF-seq study in mouse embryonic stem cells identified multiple TFs that trigger rapid transcriptome changes indicative of differentiation within a day of induction. Moreover, cells expressing these potent TFs often show a slower cell cycle and increased cell death. Further analysis using bioChIP-seq revealed that GCM1 and OTX2 act as pioneer factors and activators by increasing gene accessibility and activating the expression of lineage specification genes during cell fate conversion. iTF-seq has utility in both mapping cell fate conversion and understanding cell fate conversion mechanisms.

Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the relative contributions of mRNA transcription and degradation to shaping those profiles remains challenging, especially within embryos with diverse cellular identities. Here, we combine single-cell RNA-Seq and metabolic labeling to capture temporal cellular transcriptomes of zebrafish embryos where newly-transcribed (zygotic) and pre-existing (maternal) mRNA can be distinguished. We introduce kinetic models to quantify mRNA transcription and degradation rates within individual cell types during their specification. These models reveal highly varied regulatory rates across thousands of genes, coordinated transcription and destruction rates for many transcripts, and link differences in degradation to specific sequence elements. They also identify cell-type-specific differences in degradation, namely selective retention of maternal transcripts within primordial germ cells and enveloping layer cells, two of the earliest specified cell types. Our study provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Bento: a toolkit for subcellular analysis of spatial transcriptomics data

The spatial organization of molecules in a cell is essential for their functions. While current methods focus on discerning tissue architecture, cell-cell interactions, and spatial expression patterns, they are limited to the multicellular scale. We present Bento, a Python toolkit that takes advantage of single-molecule information to enable spatial analysis at the subcellular scale. Bento ingests molecular coordinates and segmentation boundaries to perform three analyses: defining subcellular domains, annotating localization patterns, and quantifying gene-gene colocalization. We demonstrate MERFISH, seqFISH + , Molecular Cartography, and Xenium datasets. Bento is part of the open-source Scverse ecosystem, enabling integration with other single-cell analysis tools.

Cancer-associated fibroblasts rewire the estrogen receptor response in luminal breast cancer, enabling estrogen independence

Advanced breast cancers represent a major therapeutic challenge due to their refractoriness to treatment. Cancer-associated fibroblasts (CAFs) are the most abundant constituents of the tumor microenvironment and have been linked to most hallmarks of cancer. However, the influence of CAFs on therapeutic outcome remains largely unchartered. Here, we reveal that spatial coincidence of abundant CAF infiltration with malignant cells was associated with reduced estrogen receptor (ER)-α expression and activity in luminal breast tumors. Notably, CAFs mediated estrogen-independent tumor growth by selectively regulating ER-α signaling. Whereas most prototypical estrogen-responsive genes were suppressed, CAFs maintained gene expression related to therapeutic resistance, basal-like differentiation, and invasion. A functional drug screen in co-cultures identified effector pathways involved in the CAF-induced regulation of ER-α signaling. Among these, the Transforming Growth Factor-β and the Janus kinase signaling cascades were validated as actionable targets to counteract the CAF-induced modulation of ER-α activity. Finally, genes that were downregulated in cancer cells by CAFs were predictive of poor response to endocrine treatment. In conclusion, our work reveals that CAFs directly control the luminal breast cancer phenotype by selectively modulating ER-α expression and transcriptional function, and further proposes novel targets to disrupt the crosstalk between CAFs and tumor cells to reinstate treatment response to endocrine therapy in patients.

Investigating the impact of vitrification on bovine ovarian tissue morphology, follicle survival, and transcriptomic signature

PURPOSE

Ovarian tissue cryopreservation is vital for fertility preservation, yet its effect on ovarian tissue follicle survival and transcriptomic signature requires further investigation. This study delves into the effects of vitrification on tissue morphology, function, and transcriptomic changes, helping to find possibilities for vitrification protocol improvements.

METHODS

Ovarian cortex from 19 bovine animals were used to conduct pre- and post-vitrification culture followed by histological assessment, immunohistochemistry, and TUNEL assay. Follicles' functionality was assessed for viability and growth within the tissue and in isolated cultures. RNA-sequencing of ovarian tissue was used to explore the transcriptomic alterations caused by vitrification.

RESULTS

Follicle density, cell proliferation, and DNA damage in ovarian stroma were unaffected by vitrification. However, vitrified cultured tissue exhibited reduced follicle density of primordial/primary and antral follicles, while freshly cultured tissue manifested reduction of antral follicles. Increased stromal cell proliferation and DNA damage occurred in both groups post-culture. Isolated follicles from vitrified tissue exhibited similar viability to fresh follicles until day 4, after which the survival dropped. RNA-sequencing revealed minor effects of vitrification on transcriptomic signatures, while culture induced significant gene expression changes in both groups. The altered expression of WNT and hormonal regulation pathway genes post-vitrification suggests the molecular targets for vitrification protocol refinement.

CONCLUSION

Vitrification minimally affects tissue morphology, follicle density, and transcriptomic signature post-thawing. However, culture revealed notable changes in vitrified tissue samples, including reduced follicle density, decreased isolated follicle survival, and alteration in WNT signalling and ovarian hormonal regulation pathways, highlighted them as possible limitations of the current vitrification protocol.

Cytokine profile in peripheral blood mononuclear cells differs between embryo donor and potential recipient sows

Introduction

Pregnancy success relies on the establishment of a delicate immune balance that requires the early activation of a series of local and systemic immune mechanisms. The changes in the immunological profile that are normally occurring in the pregnant uterus does not take place in cyclic (non-pregnant) uterus, a fact that has been widely explored in pigs at the tissue local level. Such differences would be especially important in the context of embryo transfer (ET), where a growing body of literature indicates that immunological differences at the uterine level between donors and recipients may significantly impact embryonic mortality. However, whether components of peripheral immunity also play a role in this context remains unknown. Accordingly, our hypothesis is that the immune status of donor sows differs from potential recipients, not only at the tissue local level but also at the systemic level. These differences could contribute to the high embryonic mortality rates occurring in ET programs.

Methods

In this study differences in systemic immunity, based on cytokine gene expression profile in peripheral blood mononuclear cells (PBMCs), between embryo-bearing donor (DO group; N = 10) and potential recipient sows (RE group; N = 10) at Day 6 after the onset of the estrus were explored. Gene expression analysis was conducted for 6 proinflammatory (IL-1α, IL-1β, IL-2, GM-CSF, IFN-γ, and TNF-α) and 6 anti-inflammatory (IL-4, IL-6, IL-10, IL-13, TGF-β1, and LIF) cytokines.

Results and discussion

All cytokines were overexpressed in the DO group except for IL-4, suggesting that stimuli derived from the insemination and/or the resultant embryos modify the systemic immune profile in DO sows compared to RE (lacking these stimuli). Our results also suggest that certain cytokines (e.g., IL-1α and IL-1β) might have a predictive value for the pregnancy status.

Search for chromosomal instability aiding variants reveal naturally occurring kinetochore gene variants that perturb chromosome segregation

Chromosomal instability (CIN) is a hallmark of cancers, and CIN-promoting mutations are not fully understood. Here, we report 141 chromosomal instability aiding variant (CIVa) candidates by assessing the prevalence of loss-of-function (LoF) variants in 135 chromosome segregation genes from over 150,000 humans. Unexpectedly, we observe both heterozygous and homozygous CIVa in Astrin and SKA3, two evolutionarily conserved kinetochore and microtubule-associated proteins essential for chromosome segregation. To stratify harmful versus harmless variants, we combine live-cell microscopy and controlled protein expression. We find the naturally occurring Astrin p.Q1012∗ variant is harmful as it fails to localize normally and induces chromosome misalignment and missegregation, in a dominant negative manner. In contrast, the Astrin p.L7Qfs∗21 variant generates a shorter isoform that localizes and functions normally, and the SKA3 p.Q70Kfs∗7 variant allows wild-type SKA complex localisation and function, revealing distinct resilience mechanisms that render these variants harmless. Thus, we present a scalable framework to predict and stratify naturally occurring CIVa, and provide insight into resilience mechanisms that compensate for naturally occurring CIVa.

Gigaxonin Suppresses Epithelial-to-Mesenchymal Transition of Human Cancer Through Downregulation of Snail

Gigaxonin is an E3 ubiquitin ligase that plays a role in cytoskeletal stability. Its role in cancer is not yet clearly understood. Our previous studies of head and neck cancer had identified gigaxonin interacting with p16 for NFκB ubiquitination. To explore its role in cancer cell growth suppression, we analyzed normal and tumor DNA from cervical and head and neck cancers. There was a higher frequency of exon 8 SNP (c.1293 C>T, rs2608555) in the tumor (46% vs. 25% normal, P = 0.011) pointing to a relationship to cancer. Comparison of primary tumor with recurrence and metastasis did not reveal a statistical significance. Two cervical cancer cell lines, ME180 and HT3 harboring exon 8 SNP and showing T allele expression correlated with higher gigaxonin expression, reduced in vitro cell growth and enhanced cisplatin sensitivity in comparison with C allele expressing cancer cell lines. Loss of gigaxonin expression in ME180 cells through CRISPR-Cas9 or siRNA led to aggressive cancer cell growth including increased migration and Matrigel invasion. The in vitro cell growth phenotypes were reversed with re-expression of gigaxonin. Suppression of cell growth correlated with reduced Snail and increased e-cadherin expression. Mouse tail vein injection studies showed increased lung metastasis of cells with low gigaxonin expression and reduced metastasis with reexpression of gigaxonin. We have found an association between C allele expression and RNA instability and absence of multimeric protein formation. From our results, we conclude that gigaxonin expression is associated with suppression of epithelial-mesenchymal transition through inhibition of Snail.

SIGNIFICANCE

Our results suggest that GAN gene exon 8 SNP T allele expression correlates with higher gigaxonin expression and suppression of aggressive cancer cell growth. There is downregulation of Snail and upregulation of e-cadherin through NFκB ubiquitination. We hypothesize that exon 8 T allele and gigaxonin expression could serve as diagnostic markers of suppression of aggressive growth of head and neck cancer.

Evaluating the regulatory function of non-coding autism-associated single nucleotide polymorphisms on gene expression in human brain tissue

Common variants account for most of the estimated heritability associated with autism spectrum disorder (autism). Although several replicable single nucleotide polymorphisms (SNPs) for the condition have been detected using genome-wide association study (GWAS) methodologies, their pathophysiological relevance remains elusive. Examining this is complicated, however, as all detected loci are situated within non-coding regions of the genome. It is therefore likely that they possess roles of regulatory function as opposed to directly affecting gene coding sequences. To bridge the gap between SNP discovery and mechanistic insight, we applied a comprehensive bioinformatic pipeline to functionally annotate autism-associated polymorphisms and their non-coding linkage disequilibrium (i.e., non-randomly associated) partners. We identified 82 DNA variants of probable regulatory function that may contribute to autism pathogenesis. To validate these predictions, we measured the impact of 11 high-confidence candidates and their GWAS linkage disequilibrium partners on gene expression in human brain tissue from Autistic and non-Autistic donors. Although a small number of the surveyed variants exhibited measurable influence on gene expression as determined via quantitative polymerase chain reaction, these did not survive correction for multiple comparisons. Additionally, no significant genotype-by-diagnosis effects were observed for any of the SNP-gene associations. We contend that this may reflect an inability to effectively capture the modest, neurodevelopmental-specific impact of individual variants on biological dysregulation in available post-mortem tissue samples, as well as limitations in the existing autism GWAS data.

Quiescence enables unrestricted cell fate in naive embryonic stem cells

Quiescence in stem cells is traditionally considered as a state of inactive dormancy or with poised potential. Naive mouse embryonic stem cells (ESCs) can enter quiescence spontaneously or upon inhibition of MYC or fatty acid oxidation, mimicking embryonic diapause in vivo. The molecular underpinning and developmental potential of quiescent ESCs (qESCs) are relatively unexplored. Here we show that qESCs possess an expanded or unrestricted cell fate, capable of generating both embryonic and extraembryonic cell types (e.g., trophoblast stem cells). These cells have a divergent metabolic landscape comparing to the cycling ESCs, with a notable decrease of the one-carbon metabolite S-adenosylmethionine. The metabolic changes are accompanied by a global reduction of H3K27me3, an increase of chromatin accessibility, as well as the de-repression of endogenous retrovirus MERVL and trophoblast master regulators. Depletion of methionine adenosyltransferase Mat2a or deletion of Eed in the polycomb repressive complex 2 results in removal of the developmental constraints towards the extraembryonic lineages. Our findings suggest that quiescent ESCs are not dormant but rather undergo an active transition towards an unrestricted cell fate.

A human stomach cell type transcriptome atlas

BACKGROUND

The identification of cell type-specific genes and their modification under different conditions is central to our understanding of human health and disease. The stomach, a hollow organ in the upper gastrointestinal tract, provides an acidic environment that contributes to microbial defence and facilitates the activity of secreted digestive enzymes to process food and nutrients into chyme. In contrast to other sections of the gastrointestinal tract, detailed descriptions of cell type gene enrichment profiles in the stomach are absent from the major single-cell sequencing-based atlases.

RESULTS

Here, we use an integrative correlation analysis method to predict human stomach cell type transcriptome signatures using unfractionated stomach RNAseq data from 359 individuals. We profile parietal, chief, gastric mucous, gastric enteroendocrine, mitotic, endothelial, fibroblast, macrophage, neutrophil, T-cell, and plasma cells, identifying over 1600 cell type-enriched genes.

CONCLUSIONS

We uncover the cell type expression profile of several non-coding genes strongly associated with the progression of gastric cancer and, using a sex-based subset analysis, uncover a panel of male-only chief cell-enriched genes. This study provides a roadmap to further understand human stomach biology.

Transcriptomic analysis of intestinal organoids, derived from pigs divergent in feed efficiency, and their response to Escherichia coli

BACKGROUND

There is increasing interest in using intestinal organoids to study complex traits like feed efficiency (FE) and host-microbe interactions. The aim of this study was to investigate differences in the molecular phenotype of organoids derived from pigs divergent for FE as well as their responses to challenge with adherent and invasive Escherichia coli (E. coli).

RESULTS

Colon and ileum tissue from low and high FE pigs was used to generate 3D organoids and two dimensional (2D) monolayers of organoid cells for E. coli challenge. Genome-wide gene expression was used to investigate molecular differences between pigs that were phenotypically divergent for FE and to study the difference in gene expression after challenge with E. coli. We showed, (1) minor differences in gene expression of colon organoids from pigs with low and high FE phenotypes, (2) that an E. coli challenge results in a strong innate immune gene response in both colon and ileum organoids, (3) that the immune response seems to be less pronounced in the colon organoids of high FE pigs and (4) a slightly stronger immune response was observed in ileum than in colon organoids.

CONCLUSIONS

These findings demonstrate the potential for using organoids to gain insights into complex biological mechanisms such as FE.

Leukocyte Telomere Length and Cardiac Structure and Function: A Mendelian Randomization Study

BACKGROUND

Existing research demonstrates the association of shorter leukocyte telomere length with increased risk of age-related health outcomes including cardiovascular diseases. However, the direct causality of these relationships has not been definitively established. Cardiovascular aging at an organ level may be captured using image-derived phenotypes of cardiac anatomy and function.

METHODS AND RESULTS

In the current study, we use 2-sample Mendelian randomization to assess the causal link between leukocyte telomere length and 54 cardiac magnetic resonance imaging measures representing structure and function across the 4 cardiac chambers. Genetically predicted shorter leukocyte telomere length was causally linked to smaller ventricular cavity sizes including left ventricular end-systolic volume, left ventricular end-diastolic volume, lower left ventricular mass, and pulmonary artery. The association with left ventricular mass (β =0.217, Pfalse discovery rate=0.016) remained significant after multiple testing adjustment, whereas other associations were attenuated.

CONCLUSIONS

Our findings support a causal role for shorter leukocyte telomere length and faster cardiac aging, with the most prominent relationship with left ventricular mass.

Uncovering the architecture of selection in two Bos taurus cattle breeds

Directional selection alters the genome via hard sweeps, soft sweeps, and polygenic selection. However, mapping polygenic selection is difficult because it does not leave clear signatures on the genome like a selective sweep. In populations with temporally stratified genotypes, the Generation Proxy Selection Mapping (GPSM) method identifies variants associated with generation number (or appropriate proxy) and thus variants undergoing directional allele frequency changes. Here, we use GPSM on two large datasets of beef cattle to detect associations between an animal's generation and 11 million imputed SNPs. Using these datasets with high power and dense mapping resolution, GPSM detected a total of 294 unique loci actively under selection in two cattle breeds. We observed that GPSM has a high power to detect selection in the very recent past (<10 years), even when allele frequency changes are small. Variants identified by GPSM reside in genomic regions associated with known breed-specific selection objectives, such as fertility and maternal ability in Red Angus, and carcass merit and coat color in Simmental. Over 60% of the selected loci reside in or near (<50 kb) annotated genes. Using haplotype-based and composite selective sweep statistics, we identify hundreds of putative selective sweeps that likely occurred earlier in the evolution of these breeds; however, these sweeps have little overlap with recent polygenic selection. This makes GPSM a complementary approach to sweep detection methods when temporal genotype data are available. The selected loci that we identify across methods demonstrate the complex architecture of selection in domesticated cattle.

Protein interaction networks in the vasculature prioritize genes and pathways underlying coronary artery disease

Population-based association studies have identified many genetic risk loci for coronary artery disease (CAD), but it is often unclear how genes within these loci are linked to CAD. Here, we perform interaction proteomics for 11 CAD-risk genes to map their protein-protein interactions (PPIs) in human vascular cells and elucidate their roles in CAD. The resulting PPI networks contain interactions that are outside of known biology in the vasculature and are enriched for genes involved in immunity-related and arterial-wall-specific mechanisms. Several PPI networks derived from smooth muscle cells are significantly enriched for genetic variants associated with CAD and related vascular phenotypes. Furthermore, the networks identify 61 genes that are found in genetic loci associated with risk of CAD, prioritizing them as the causal candidates within these loci. These findings indicate that the PPI networks we have generated are a rich resource for guiding future research into the molecular pathogenesis of CAD.

Recent Developments and Applications of Biocatalytic and Chemoenzymatic Synthesis for the Generation of Diverse Classes of Drugs

Biocatalytic and chemoenzymatic biosynthesis are powerful methods of organic chemistry that use enzymes to execute selective reactions and allow the efficient production of organic compounds. The advantages of these approaches include high selectivity, mild reaction conditions, and the ability to work with complex substrates. The utilization of chemoenzymatic techniques for the synthesis of complicated compounds has lately increased dramatically in the area of organic chemistry. Biocatalytic technologies and modern synthetic methods are utilized synergistically in a multi-step approach to a target molecule under this paradigm. Chemoenzymatic techniques are promising for simplifying access to essential bioactive compounds because of the remarkable regio- and stereoselectivity of enzymatic transformations and the reaction diversity of modern organic chemistry. Enzyme kits may include ready-to-use, reproducible biocatalysts. Its use opens up new avenues for the synthesis of active therapeutic compounds and aids in drug development by synthesizing active components to construct scaffolds in a targeted and preparative manner. This study summarizes current breakthroughs as well as notable instances of biocatalytic and chemoenzymatic synthesis. To assist organic chemists in the use of enzymes for synthetic applications, it also provides some basic guidelines for selecting the most appropriate enzyme for a targeted reaction while keeping aspects like cofactor requirement, solvent tolerance, use of whole cell or isolated enzymes, and commercial availability in mind.

Olfactory ensheathing cells are hybrid glial cells that promote neural repair

Olfactory ensheathing cells (OECs) are unique glial cells found in both the central and peripheral nervous systems where they support the continuous axonal outgrowth of immature olfactory sensory neurons to their targets. Here we show that following severe spinal cord injury, olfactory bulb-derived OECs transplanted near the injury site modify the normally inhibitory glial scar and facilitate axon regeneration past the scar border and into the lesion center. To understand the mechanisms underlying the reparative properties of such transplanted OECs, we used single-cell RNA-sequencing to study their gene expression programs. Our analyses revealed five diverse subtypes of OECs, each expressing novel marker genes and pathways indicative of progenitor, axonal regeneration and repair, secreted molecules, or microglia-like functions. As expected, we found substantial overlap of OEC genes with those of Schwann cells, but also with astrocytes, oligodendrocytes and microglia. We confirmed established markers on cultured OECs, and then localized select top genes of OEC subtypes in rat olfactory bulb tissue. In addition, we present evidence that OECs secrete both Reelin and Connective tissue growth factor, extracellular matrix molecules which are important for neural repair and axonal outgrowth. Our results support that adult OECs are a unique hybrid glia, some with progenitor characteristics, and that their gene expression patterns indicate diverse functions related to wound healing, injury repair and axonal regeneration.

Germline gene fusions across species reveal the chromosomal instability regions and cancer susceptibility

The canine transmissible venereal tumor (CTVT) is a clonal cell-mediated cancer with a long evolutionary history and extensive karyotype rearrangements in its genome. However, little is known about its genetic similarity to human tumors. Here, using multi-omics data we identified 11 germline gene fusions (GGFs) in CTVT, which showed higher genetic susceptibility than others. Additionally, we illustrate a mechanism of a complex gene fusion of three gene segments (HSD17B4-DMXL1-TNFAIP8) that we refer to "greedy fusion". Our findings also provided evidence that expressions of GGFs are downregulated during the tumor regressive phase, which is associated with DNA methylation level. This study presents a comprehensive landscape of gene fusions (GFs) in CTVT, which offers a valuable genetic resource for exploring potential genetic mechanisms underlying the development of cancers in both dogs and humans.

TRAIL-dependent apoptosis of peritoneal mesothelial cells by NK cells promotes ovarian cancer invasion

A crucial requirement for metastasis formation in ovarian high-grade serous carcinoma (HGSC) is the disruption of the protective peritoneal mesothelium. Using co-culture systems of primary human cells, we discovered that tumor-associated NK cells induce TRAIL-dependent apoptosis in mesothelial cells via death receptors DR4 and DR5 upon encounter with activated T cells. Upregulation of TRAIL expression in NK cells concomitant with enhanced cytotoxicity toward mesothelial cells was driven predominantly by T-cell-derived TNFα, as shown by affinity proteomics-based analysis of the T cell secretome in conjunction with functional studies. Consistent with these findings, we detected apoptotic mesothelial cells in the peritoneal fluid of HGSC patients. In contrast to mesothelial cells, HGSC cells express negligible levels of both DR4 and DR5 and are TRAIL resistant, indicating cell-type-selective killing by NK cells. Our data point to a cooperative action of T and NK in breaching the mesothelial barrier in HGSC patients.

Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

Efficacy and safety of a SOD1-targeting artificial miRNA delivered by AAV9 in mice are impacted by miRNA scaffold selection

Toxic gain-of-function mutations in superoxide dismutase 1 (SOD1) contribute to approximately 2%-3% of all amyotrophic lateral sclerosis (ALS) cases. Artificial microRNAs (amiRs) delivered by adeno-associated virus (AAV) have been proposed as a potential treatment option to silence SOD1 expression and mitigate disease progression. Primary microRNA (pri-miRNA) scaffolds are used in amiRs to shuttle a hairpin RNA into the endogenous miRNA pathway, but it is unclear whether different primary miRNA (pri-miRNA) scaffolds impact the potency and safety profile of the expressed amiR in vivo. In our process to develop an AAV amiR targeting SOD1, we performed a preclinical characterization of two pri-miRNA scaffolds, miR155 and miR30a, sharing the same guide strand sequence. We report that, while the miR155-based vector, compared with the miR30a-based vector, leads to a higher level of the amiR and more robust suppression of SOD1 in vitro and in vivo, it also presents significantly greater risks for CNS-related toxicities in vivo. Despite miR30a-based vector showing relatively lower potency, it can significantly delay the development of ALS-like phenotypes in SOD1-G93A mice and increase survival in a dose-dependent manner. These data highlight the importance of scaffold selection in the pursuit of highly efficacious and safe amiRs for RNA interference gene therapy.

Dose rate dependent reduction in chromatin accessibility at transcriptional start sites long time after exposure to gamma radiation

Ionizing radiation (IR) impact cellular and molecular processes that require chromatin remodelling relevant for cellular integrity. However, the cellular implications of ionizing radiation (IR) delivered per time unit (dose rate) are still debated. This study investigates whether the dose rate is relevant for inflicting changes to the epigenome, represented by chromatin accessibility, or whether it is the total dose that is decisive. CBA/CaOlaHsd mice were whole-body exposed to either chronic low dose rate (2.5 mGy/h for 54 d) or the higher dose rates (10 mGy/h for 14 d and 100 mGy/h for 30 h) of gamma radiation (⁶⁰Co, total dose: 3 Gy). Chromatin accessibility was analysed in liver tissue samples using Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-Seq), both one day after and over three months post-radiation (>100 d). The results show that the dose rate contributes to radiation-induced epigenomic changes in the liver at both sampling timepoints. Interestingly, chronic low dose rate exposure to a high total dose (3 Gy) did not inflict long-term changes to the epigenome. In contrast to the acute high dose rate given to the same total dose, reduced accessibility at transcriptional start sites (TSS) was identified in genes relevant for the DNA damage response and transcriptional activity. Our findings link dose rate to essential biological mechanisms that could be relevant for understanding long-term changes after ionizing radiation exposure. However, future studies are needed to comprehend the biological consequence of these findings.

Predictive Factors and ACE-2 Gene Polymorphisms in Susceptibility to Long COVID-19 Syndrome

Long COVID-19 syndrome is present in 5-10% of patients infected with SARS-CoV-2, and there is still little information on the predisposing factors that lead to its development. The purpose of the study was to evaluate the predictive factors in early symptoms, clinical features and the role of Angiotensin-Converting Enzyme-2 (ACE-2) c.513-1451G>A (rs2106806) and c.15643279T>C (rs6629110) polymorphisms in the susceptibility to developing Long COVID-19 syndrome subsequent to COVID-19 infectionA total of 29 patients who suffered COVID-19 were recruited in a descriptive longitudinal study of two groups: Long COVID-19 (n = 16) and non-Long COVID-19 (n = 13). Early symptoms and clinical features during COVID-19 were classified by a medical service. ACE-2 polymorphisms were genotyped by using a Single Nucleotide Primer Extension (SNPE). Of the early symptoms, fatigue, myalgia and headache showed a high risk of increasing Long COVID-19 susceptibility. Clinical features such as emergency care, SARS-CoV-2 reinfection, previous diseases, respiratory disease and brain fog also had a high risk of increasing Long COVID-19 susceptibility. The A allele in the rs2106806 variant was associated with an odds ratio (OR) of 4.214 (95% CI 2.521-8.853; p < 0.001), and the T allele in the rs6629110 variant was associated with an OR of 3.754 (95% CI 1.785-6.105; p = 0.002) of increasing Long COVID-19 susceptibility. This study shows the risk of ACE-2 polymorphisms, different early symptoms and clinical features during SARS-CoV-2 infection in susceptibility to Long COVID-19.

Cellular reprogramming is driven by widespread rewiring of promoter-enhancer interactions

BACKGROUND

Long-range interactions between promoters and cis-regulatory elements, such as enhancers, play critical roles in gene regulation. However, the role of three-dimensional (3D) chromatin structure in orchestrating changes in transcriptional regulation during direct cell reprogramming is not fully understood.

RESULTS

Here, we performed integrated analyses of chromosomal architecture, epigenetics, and gene expression using Hi-C, promoter Capture Hi-C (PCHi-C), ChIP-seq, and RNA-seq during trans-differentiation of Pre-B cells into macrophages with a β-estradiol inducible C/EBPαER transgene. Within 1h of β-estradiol induction, C/EBPα translocated from the cytoplasm to the nucleus, binding to thousands of promoters and putative regulatory elements, resulting in the downregulation of Pre-B cell-specific genes and induction of macrophage-specific genes. Hi-C results were remarkably consistent throughout trans-differentiation, revealing only a small number of TAD boundary location changes, and A/B compartment switches despite significant changes in the expression of thousands of genes. PCHi-C revealed widespread changes in promoter-anchored loops with decreased interactions in parallel with decreased gene expression, and new and increased promoter-anchored interactions in parallel with increased expression of macrophage-specific genes.

CONCLUSIONS

Overall, our data demonstrate that C/EBPα-induced trans-differentiation involves few changes in genome architecture at the level of TADs and A/B compartments, in contrast with widespread reorganization of thousands of promoter-anchored loops in association with changes in gene expression and cell identity.

Cellular age explains variation in age-related cell-to-cell transcriptome variability

Organs and tissues age at different rates within a single individual. Such asynchrony in aging has been widely observed at multiple levels, from functional hallmarks, such as anatomical structures and physiological processes, to molecular endophenotypes, such as the transcriptome and metabolome. However, we lack a conceptual framework to understand why some components age faster than others. Just as demographic models explain why aging evolves, here we test the hypothesis that demographic differences among cell types, determined by cell-specific differences in turnover rate, can explain why the transcriptome shows signs of aging in some cell types but not others. Through analysis of mouse single-cell transcriptome data across diverse tissues and ages, we find that cellular age explains a large proportion of the variation in the age-related increase in transcriptome variance. We further show that long-lived cells are characterized by relatively high expression of genes associated with proteostasis and that the transcriptome of long-lived cells shows greater evolutionary constraint than short-lived cells. In contrast, in short-lived cell types, the transcriptome is enriched for genes associated with DNA repair. Based on these observations, we develop a novel heuristic model that explains how and why aging rates differ among cell types.

Fast alignment of mass spectra in large proteomics datasets, capturing dissimilarities arising from multiple complex modifications of peptides

BACKGROUND

In proteomics, the interpretation of mass spectra representing peptides carrying multiple complex modifications remains challenging, as it is difficult to strike a balance between reasonable execution time, a limited number of false positives, and a huge search space allowing any number of modifications without a priori. The scientific community needs new developments in this area to aid in the discovery of novel post-translational modifications that may play important roles in disease.

RESULTS

To make progress on this issue, we implemented SpecGlobX (SpecGlob eXTended to eXperimental spectra), a standalone Java application that quickly determines the best spectral alignments of a (possibly very large) list of Peptide-to-Spectrum Matches (PSMs) provided by any open modification search method, or generated by the user. As input, SpecGlobX reads a file containing spectra in MGF or mzML format and a semicolon-delimited spreadsheet describing the PSMs. SpecGlobX returns the best alignment for each PSM as output, splitting the mass difference between the spectrum and the peptide into one or more shifts while considering the possibility of non-aligned masses (a phenomenon resulting from many situations including neutral losses). SpecGlobX is fast, able to align one million PSMs in about 1.5 min on a standard desktop. Firstly, we remind the foundations of the algorithm and detail how we adapted SpecGlob (the method we previously developed following the same aim, but limited to the interpretation of perfect simulated spectra) to the interpretation of imperfect experimental spectra. Then, we highlight the interest of SpecGlobX as a complementary tool downstream to three open modification search methods on a large simulated spectra dataset. Finally, we ran SpecGlobX on a proteome-wide dataset downloaded from PRIDE to demonstrate that SpecGlobX functions just as well on simulated and experimental spectra. We then carefully analyzed a limited set of interpretations.

CONCLUSIONS

SpecGlobX is helpful as a decision support tool, providing keys to interpret peptides carrying complex modifications still poorly considered by current open modification search software. Better alignment of PSMs enhances confidence in the identification of spectra provided by open modification search methods and should improve the interpretation rate of spectra.

Polycomb protein binding and looping mediated by Polycomb Response Elements in the ON transcriptional state

Polycomb group proteins (PcG) mediate epigenetic silencing of important developmental genes and other targets. In Drosophila, canonical PcG-target genes contain Polycomb Response Elements (PREs) that recruit PcG protein complexes including PRC2 that trimethylates H3K27 forming large H3K27me3 domains. In the OFF transcriptional state, PREs loop with each other and this looping strengthens silencing. Here we address the question of what PcG proteins bind to PREs when canonical PcG target genes are expressed, and whether PREs loop when these genes are ON. Our data show that the answer to this question is PRE-specific but general conclusions can be made. First, within a PcG-target gene, some regulatory DNA can remain covered with H3K27me3 and PcG proteins remain bound to PREs in these regions. Second, when PREs are within H3K27ac domains, PcG-binding decreases, however, this depends on the protein and PRE. The DNA binding protein GAF, and the PcG protein Ph remain at PREs even when other PcG proteins are greatly depleted. In the ON state, PREs can still loop with each other, but also form loops with presumptive enhancers. These data support the model that, in addition to their role in PcG silencing, PREs can act as "promoter-tethering elements" mediating interactions between promoter proximal PREs and distant enhancers.

Proteogenomics Reveal the Overexpression of HLA-I in Cancer

An accurate quantification of HLA class I gene expression is important in understanding the interplay with the tumor microenvironment of antitumor cytotoxic T cell activities. Because HLA-I sequences are highly variable, standard RNAseq and mass spectrometry-based quantification workflows using common genome and protein sequence references do not provide HLA-I allele specific quantifications. Here, we used personalized HLA-I nucleotide and protein reference sequences based on the subjects' HLA-I genotypes and surveyed tumor and adjacent normal samples from patients across nine cancer types. Mass spectrometry using data dependent acquisition data was validated to be sufficient to estimate HLA-A protein expression at the allele level. We found that HLA-I proteins were present in significantly higher levels in tumors compared to adjacent normal tissues from 41 to 63% of head and neck squamous cell carcinoma, uterine corpus endometrial carcinoma, and clear cell renal cell carcinoma patients, and this was driven by increased levels of HLA-I gene transcripts. Most immune cell types are universally enriched in HLA-I high tumors, while endothelial and neuronal cells showed divergent relationships with HLA-I. Pathway analysis revealed that tumor senescence and autophagy activity influence the level of HLA-I proteins in glioblastoma. Genes correlated to HLA-I protein expression are mostly the ones directly involved in HLA-I function in immune response and cell death, while glycosylation genes are exclusively co-expressed with HLA-I at the protein level.

scDiffCom: a tool for differential analysis of cell-cell interactions provides a mouse atlas of aging changes in intercellular communication

Dysregulation of intercellular communication is a hallmark of aging. To better quantify and explore changes in intercellular communication, we present scDiffCom and scAgeCom. scDiffCom is an R package, relying on approximately 5,000 curated ligand-receptor interactions, that performs differential intercellular communication analysis between two conditions from single-cell transcriptomics data. Built upon scDiffCom, scAgeCom is an atlas of age-related cell-cell communication changes covering 23 mouse tissues from 58 single-cell RNA sequencing datasets from Tabula Muris Senis and the Calico murine aging cell atlas. It offers a comprehensive resource of tissue-specific and sex-specific aging dysregulations and highlights age-related intercellular communication changes widespread across the whole body, such as the upregulation of immune system processes and inflammation, the downregulation of developmental processes, angiogenesis and extracellular matrix organization and the deregulation of lipid metabolism. Our analysis emphasizes the relevance of the specific ligands, receptors and cell types regulating these processes. The atlas is available online ( https://scagecom.org ).

Circulating DNA reveals a specific and higher fragmentation of the Y chromosome

Chromosome stability is a key point in genome evolution, particularly that of the Y chromosome. The Y chromosome loss in blood and tumor cells is well established. Through processes that are common to other chromosomes too, the Y chromosome undergoes degradation and fragmentation in the blood stream before elimination. This process gives rise to circulating DNA (cirDNA) fragments, whose examination may provide potential insight into the role of DNA fragmentation in blood for the Y chromosome elimination. In this study, we employed shallow whole genome sequencing (sWGS) to comprehensively assess the total cirDNA and the individual chromosome fragment size profiles in the plasma of healthy male individuals. Here, we show that (i) the fragment size profiles of total circulating DNA (cirDNA) and DNA fragments originating from autosomes and the X chromosome in blood plasma are homogeneous, and have a remarkably low variability (mean CV = 7%) among healthy individuals, (ii) the Y chromosome has a distinct fragment size profile with the accumulation of the fragment < 145 bp and depletion of the dinucleosome-associated fragments (290-390 bp), and its fragment fraction in blood decreases with age. These results indicate a higher fragmentation of the Y chromosome compared to other chromosomes and this in turn might be due to its increased susceptibility to degradation. Our findings pave the way for an elucidation of the impact of chromosomal origin on DNA degradation and the Y chromosome biology.

Immune checkpoint inhibitor-induced colitis is mediated by polyfunctional lymphocytes and is dependent on an IL23/IFNγ axis

Immune checkpoint inhibitors (CPIs) are a relatively newly licenced cancer treatment, which make a once previously untreatable disease now amenable to a potential cure. Combination regimens of anti-CTLA4 and anti-PD-1 show enhanced efficacy but are prone to off-target immune-mediated tissue injury, particularly at the barrier surfaces. To probe the impact of immune checkpoints on intestinal homoeostasis, mice are challenged with anti-CTLA4 and anti-PD-1 immunotherapy and manipulation of the intestinal microbiota. The immune profile of the colon of these mice with CPI-colitis is analysed using bulk RNA sequencing, single-cell RNA sequencing and flow cytometry. CPI-colitis in mice is dependent on the composition of the intestinal microbiota and by the induction of lymphocytes expressing interferon-γ (IFNγ), cytotoxicity molecules and other pro-inflammatory cytokines/chemokines. This pre-clinical model of CPI-colitis could be attenuated following blockade of the IL23/IFNγ axis. Therapeutic targeting of IFNγ-producing lymphocytes or regulatory networks, may hold the key to reversing CPI-colitis.

Systematic perturbations of SETD2, NSD1, NSD2, NSD3 and ASH1L reveals their distinct contributions to H3K36 methylation

Background

Methylation of histone 3 lysine 36 (H3K36me) has emerged as an essential epigenetic component for the faithful regulation of gene expression. Despite its importance in development, disease, and cancer, how the molecular agents collectively shape the H3K36me landscape is unclear.

Results

We use a mouse mesenchymal stem cell model to perturb the H3K36me deposition machinery and infer the activities of the five most prominent players: SETD2, NSD1, NSD2, NSD3, and ASH1L. We find that H3K36me2 is the most abundant of the three methylation states and is predominantly deposited at intergenic regions by NSD1, and partly by NSD2. In contrast, H3K36me1/3 are most abundant within exons and are positively correlated with gene expression. We demonstrate that while SETD2 deposits most H3K36me3, it also deposits H3K36me2 within transcribed genes. Additionally, loss of SETD2 results in an increase of exonic H3K36me1, suggesting other H3K36 methyltransferases (K36MTs) prime gene bodies with lower methylation states ahead of transcription. Through a reductive approach, we uncover the distribution patterns of NSD3- and ASH1L-catalyzed H3K36me2. While NSD1/2 establish broad intergenic H3K36me2 domains, NSD3 deposits H3K36me2 peaks on active promoters and enhancers. Meanwhile, the activity of ASH1L is restricted to the regulatory elements of developmentally relevant genes, and our analyses implicate PBX2 as a potential recruitment factor.

Conclusions

Within genes, SETD2 deposits both H3K36me2/3, while the other K36MTs are capable of depositing H3K36me1/2 independently of SETD2 activity. For the deposition of H3K36me1/2, we find a hierarchy of K36MT activities where NSD1>NSD2>NSD3>ASH1L. While NSD1 and NSD2 are responsible for most genome-wide propagation of H3K36me2, the activities of NSD3 and ASH1L are confined to active regulatory elements.

Unraveling the relationships between alpha- and beta-adrenergic modulation and the risk of heart failure

Background

The effects of α and ß adrenergic receptor modulation on the risk of developing heart failure (HF) remains uncertain due to a lack of randomized controlled trials. This study aimed to estimate the effects of α and ß adrenergic receptors modulation on the risk of HF and to provide proof of principle for genetic target validation studies in HF.

Methods

Genetic variants within the cis regions encoding the adrenergic receptors α1A, α2B, ß1, and ß2 associated with blood pressure in a 757,601-participant genome-wide association study (GWAS) were selected as instruments to perform a drug target Mendelian randomization study. Effects of these variants on HF risk were derived from the HERMES GWAS (542,362 controls; 40,805 HF cases).

Results

Lower α1A or ß1 activity was associated with reduced HF risk: odds ratio (OR) 0.83 (95% CI 0.74-0.93, P = 0.001) and 0.95 (95% CI 0.93-0.97, P = 8 × 10-6). Conversely, lower α2B activity was associated with increased HF risk: OR 1.09 (95% CI 1.05-1.12, P = 3 × 10-7). No evidence of an effect of lower ß2 activity on HF risk was found: OR 0.99 (95% CI 0.92-1.07, P = 0.95). Complementary analyses showed that these effects were consistent with those on left ventricular dimensions and acted independently of any potential effect on coronary artery disease.

Conclusions

This study provides genetic evidence that α1A or ß1 receptor inhibition will likely decrease HF risk, while lower α2B activity may increase this risk. Genetic variant analysis can assist with drug development for HF prevention.

Immune-related Adverse Events after Immune Checkpoint Blockade-based Therapy Are Associated with Improved Survival in Advanced Sarcomas

The association between immune-related AEs (irAE) and outcome in patients with sarcoma is not known. We retrospectively reviewed a cohort of patients with advanced sarcoma treated with immune checkpoint blockade (ICB)-based therapy. Association of irAEs with survival was assessed using a Cox regression model that incorporated irAE occurrence as a time-dependent covariate. Tumor samples with available RNA sequencing data were stratified by presence of an irAE to identify patterns of differential gene expression. A total of 131 patients were included. Forty-two (32%) had at least one irAE of any grade and 16 (12%) had at least one grade ≥ 3 irAE. The most common irAEs were hypothyroidism (8.3%), arthralgias (5.3%), pneumonitis (4.6%), allergic reaction (3.8%), and elevated transaminases (3.8%). Median progression-free survival (PFS) and overall survival (OS) from the time of study entry were 11.4 [95% confidence interval (CI), 10.7-15.0) and 74.6 weeks (CI, 44.9-89.7), respectively. On Cox analysis adjusting for clinical covariates that were significant in the univariate setting, the HR for an irAE (HR, 0.662; CI, 0.421-1.041) approached, but did not reach statistical significance for PFS (P = 0.074). Patients had a significantly lower HR for OS (HR, 0.443; CI, 0.246-0.798; P = 0.007) compared with those without or before an irAE. Gene expression profiling on baseline tumor samples found that patients who had an irAE had higher numbers of tumor-infiltrating dendritic cells, CD8+ T cells, and regulatory T cells as well as upregulation of immune and inflammatory pathways.

SIGNIFICANCE

irAE after ICB therapy was associated with an improved OS; it also approached statistical significance for improved PFS. Patients who had an irAE were more likely to have an inflamed tumor microenvironment at baseline.

Niclosamide as a chemical probe for analyzing SARS-CoV-2 modulation of host cell lipid metabolism

Introduction

SARS-CoV-2 subverts host cell processes to facilitate rapid replication and dissemination, and this leads to pathological inflammation.

Methods

We used niclosamide (NIC), a poorly soluble anti-helminth drug identified initially for repurposed treatment of COVID-19, which activates the cells' autophagic and lipophagic processes as a chemical probe to determine if it can modulate the host cell's total lipid profile that would otherwise be either amplified or reduced during SARS-CoV-2 infection.

Results

Through parallel lipidomic and transcriptomic analyses we observed massive reorganization of lipid profiles of SARS-CoV-2 infected Vero E6 cells, especially with triglycerides, which were elevated early during virus replication, but decreased thereafter, as well as plasmalogens, which were elevated at later timepoints during virus replication, but were also elevated under normal cell growth. These findings suggested a complex interplay of lipid profile reorganization involving plasmalogen metabolism. We also observed that NIC treatment of both low and high viral loads does not affect virus entry. Instead, NIC treatment reduced the abundance of plasmalogens, diacylglycerides, and ceramides, which we found elevated during virus infection in the absence of NIC, resulting in a significant reduction in the production of infectious virions. Unexpectedly, at higher viral loads, NIC treatment also resulted in elevated triglyceride levels, and induced significant changes in phospholipid metabolism.

Discussion

We posit that future screens of approved or new partner drugs should prioritize compounds that effectively counter SARS-CoV-2 subversion of lipid metabolism, thereby reducing virus replication, egress, and the subsequent regulation of key lipid mediators of pathological inflammation.

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

Gene duplication is recognized as a critical process in genome evolution; however, many questions about this process remain unanswered. Although gene duplicability has been observed to differ by duplication mechanism and evolutionary rate, there is so far no broad characterization of its determinants. Many features correlate with this difference in duplicability; however, our ability to exploit these observations to advance our understanding of the role of duplication in evolution is hampered by limitations within existing work. In particular, the existence of methodological differences across studies impedes meaningful comparison. Here, we use consistent definitions of duplicability in the human lineage to explore these associations, allow resolution of the impact of confounding factors, and define the overall relevance of individual features. Using a classifier approach and controlling for the confounding effect of duplicate longevity, we find a subset of gene features important in differentiating genes duplicable by small-scale duplication from those duplicable by whole-genome duplication, revealing critical roles for gene dosage and expression costs in duplicability. We further delve into patterns of functional enrichment and find a lack of constraint on duplicate retention in any context for genes duplicable by small-scale duplication.

Assessing the Presence of Recent Adaptation in the Human Genome With Mixture Density Regression

How much genome differences between species reflect neutral or adaptive evolution is a central question in evolutionary genomics. In humans and other mammals, the presence of adaptive versus neutral genomic evolution has proven particularly difficult to quantify. The difficulty notably stems from the highly heterogeneous organization of mammalian genomes at multiple levels (functional sequence density, recombination, etc.) which complicates the interpretation and distinction of adaptive versus neutral evolution signals. In this study, we introduce mixture density regressions (MDRs) for the study of the determinants of recent adaptation in the human genome. MDRs provide a flexible regression model based on multiple Gaussian distributions. We use MDRs to model the association between recent selection signals and multiple genomic factors likely to affect the occurrence/detection of positive selection, if the latter was present in the first place to generate these associations. We find that an MDR model with two Gaussian distributions provides an excellent fit to the genome-wide distribution of a common sweep summary statistic (integrated haplotype score), with one of the two distributions likely enriched in positive selection. We further find several factors associated with signals of recent adaptation, including the recombination rate, the density of regulatory elements in immune cells, GC content, gene expression in immune cells, the density of mammal-wide conserved elements, and the distance to the nearest virus-interacting gene. These results support the presence of strong positive selection in recent human evolution and highlight MDRs as a powerful tool to make sense of signals of recent genomic adaptation.

Insights into the genetics of menopausal vasomotor symptoms: genome-wide analyses of routinely-collected primary care health records

BACKGROUND

Vasomotor symptoms (VMS) can often significantly impact women's quality of life at menopause. In vivo studies have shown that increased neurokinin B (NKB) / neurokinin 3 receptor (NK3R) signalling contributes to VMS, with previous genetic studies implicating the TACR3 gene locus that encodes NK3R. Large-scale genomic analyses offer the possibility of biological insights but few such studies have collected data on VMS, while proxy phenotypes such as hormone replacement therapy (HRT) use are likely to be affected by changes in clinical practice. We investigated the genetic basis of VMS by analysing routinely-collected health records.

METHODS

We performed a GWAS of VMS derived from linked primary-care records and cross-sectional self-reported HRT use in up to 153,152 women from UK Biobank, a population-based cohort. In a subset of this cohort (n = 39,356), we analysed exome-sequencing data to test the association with VMS of rare deleterious genetic variants. Finally, we used Mendelian randomisation analysis to investigate the reasons for HRT use over time.

RESULTS

Our GWAS of health-records derived VMS identified a genetic signal near TACR3 associated with a lower risk of VMS (OR=0.76 (95% CI 0.72,0.80) per A allele, P=3.7x10-27), which was consistent with previous studies, validating this approach. Conditional analyses demonstrated independence of genetic signals for puberty timing and VMS at the TACR3 locus, including a rare variant predicted to reduce functional NK3R levels that was associated with later menarche (P = 5 × 10-9) but showed no association with VMS (P = 0.6). Younger menopause age was causally-associated with greater HRT use before 2002 but not after.

CONCLUSIONS

We provide support for TACR3 in the genetic basis of VMS but unexpectedly find that rare genomic variants predicted to lower NK3R levels did not modify VMS, despite the proven efficacy of NK3R antagonists. Using genomics we demonstrate changes in genetic associations with HRT use over time, arising from a change in clinical practice since the early 2000s, which is likely to reflect a switch from preventing post-menopausal complications in women with earlier menopause to primarily treating VMS. Our study demonstrates that integrating routinely-collected primary care health records and genomic data offers great potential for exploring the genetic basis of symptoms.

From maternal glucocorticoid and thyroid hormones to epigenetic regulation of offspring gene expression: An experimental study in a wild bird species

Offspring phenotype at birth is determined by its genotype and the prenatal environment including exposure to maternal hormones. Variation in both maternal glucocorticoids and thyroid hormones can affect offspring phenotype, but the underlying molecular mechanisms, especially those contributing to long-lasting effects, remain unclear. Epigenetic changes (such as DNA methylation) have been postulated as mediators of long-lasting effects of early-life environment. In this study, we determined the effects of elevated prenatal glucocorticoid and thyroid hormones on handling stress response (breath rate) as well as DNA methylation and gene expression of glucocorticoid receptor (GR) and thyroid hormone receptor (THR) in great tits (Parus major). Eggs were injected before incubation onset with corticosterone (the main avian glucocorticoid) and/or thyroid hormones (thyroxine and triiodothyronine) to simulate variation in maternal hormone deposition. Breath rate during handling and gene expression of GR and THR were evaluated 14 days after hatching. Methylation status of GR and THR genes was analyzed from the longitudinal blood cells sampled 7 and 14 days after hatching, as well as the following autumn. Elevated prenatal corticosterone level significantly increased the breath rate during handling, indicating an enhanced metabolic stress response. Prenatal corticosterone manipulation had CpG-site-specific effects on DNA methylation at the GR putative promoter region, while it did not significantly affect GR gene expression. GR expression was negatively associated with earlier hatching date and chick size. THR methylation or expression did not exhibit any significant relationship with the hormonal treatments or the examined covariates, suggesting that TH signaling may be more robust due to its crucial role in development. This study provides some support to the hypothesis suggesting that maternal corticosterone may influence offspring metabolic stress response via epigenetic alterations, yet their possible adaptive role in optimizing offspring phenotype to the prevailing conditions, context-dependency, and the underlying molecular interplay needs further research.

A lamprey neural cell type atlas illuminates the origins of the vertebrate brain

The vertebrate brain emerged more than ~500 million years ago in common evolutionary ancestors. To systematically trace its cellular and molecular origins, we established a spatially resolved cell type atlas of the entire brain of the sea lamprey-a jawless species whose phylogenetic position affords the reconstruction of ancestral vertebrate traits-based on extensive single-cell RNA-seq and in situ sequencing data. Comparisons of this atlas to neural data from the mouse and other jawed vertebrates unveiled various shared features that enabled the reconstruction of cell types, tissue structures and gene expression programs of the ancestral vertebrate brain. However, our analyses also revealed key tissues and cell types that arose later in evolution. For example, the ancestral brain was probably devoid of cerebellar cell types and oligodendrocytes (myelinating cells); our data suggest that the latter emerged from astrocyte-like evolutionary precursors in the jawed vertebrate lineage. Altogether, our work illuminates the cellular and molecular architecture of the ancestral vertebrate brain and provides a foundation for exploring its diversification during evolution.

Regulation of the three-dimensional chromatin organization by transposable elements in pig spleen

Like other mammalian species, the pig genome is abundant with transposable elements (TEs). The importance of TEs for three-dimensional (3D) chromatin organization has been observed in species like human and mouse, yet current understanding about pig TEs is absent. Here, we investigated the contribution of TEs for the 3D chromatin organization in three pig tissues, focusing on spleen which is crucial for both adaptive and innate immunity. We identified dozens of TE families overrepresented with CTCF binding sites, including LTR22_SS, LTR15_SS and LTR16_SSc which are pig-specific families of endogenous retroviruses (ERVs). Interestingly, LTR22_SS elements harbor a CTCF motif and create hundreds of CTCF binding sites that are associated with adaptive immunity. We further applied Hi-C to profile the 3D chromatin structure in spleen and found that TE-derived CTCF binding sites correlate with chromatin insulation and frequently overlap TAD borders and loop anchors. Notably, one LTR22_SS-derived CTCF binding site demarcate a TAD boundary upstream of XCL1, which is a spleen-enriched chemokine gene important for lymphocyte trafficking and inflammation. Overall, this study represents a first step toward understanding the function of TEs on 3D chromatin organization regulation in pigs and expands our understanding about the functional importance of TEs in mammals.

Mobilisation of jerboa kidney gene networks during dehydration and opportunistic rehydration

Desert animals have evolved systems that enable them to thrive under dry conditions. Focusing on the kidney, we have investigated the transcriptomic adaptations that enable a desert rodent, the Lesser Egyptian Jerboa (Jaculus jaculus), to withstand water deprivation and opportunistic rehydration. Analysis of the whole kidney transcriptome showed many differentially expressed genes in the Jerboa kidney, 6.4% of genes following dehydration and an even greater number (36.2%) following rehydration compared to control. Genes correlated with the rehydration condition included many ribosomal protein coding genes suggesting a concerted effort to accelerate protein synthesis when water is made available. We identify an increase in TGF-beta signaling antagonists in dehydration (e.g., GREM2). We also describe expression of multiple aquaporin and solute carrier transporters mapped to specific nephron segments. The desert adapted renal transcriptome presented here is a valuable resource to expand our understanding of osmoregulation beyond that derived from model organisms.

Comprehensive single-cell genome analysis at nucleotide resolution using the PTA Analysis Toolbox

Detection of somatic mutations in single cells has been severely hampered by technical limitations of whole-genome amplification. Novel technologies including primary template-directed amplification (PTA) significantly improved the accuracy of single-cell whole-genome sequencing (WGS) but still generate hundreds of artifacts per amplification reaction. We developed a comprehensive bioinformatic workflow, called the PTA Analysis Toolbox (PTATO), to accurately detect single base substitutions, insertions-deletions (indels), and structural variants in PTA-based WGS data. PTATO includes a machine learning approach and filtering based on recurrence to distinguish PTA artifacts from true mutations with high sensitivity (up to 90%), outperforming existing bioinformatic approaches. Using PTATO, we demonstrate that hematopoietic stem cells of patients with Fanconi anemia, which cannot be analyzed using regular WGS, have normal somatic single base substitution burdens but increased numbers of deletions. Our results show that PTATO enables studying somatic mutagenesis in the genomes of single cells with unprecedented sensitivity and accuracy.

Genetic variability in sporadic amyotrophic lateral sclerosis

With the advent of gene therapies for amyotrophic lateral sclerosis (ALS), there is a surge in gene testing for this disease. Although there is ample experience with gene testing for C9orf72, SOD1, FUS and TARDBP in familial ALS, large studies exploring genetic variation in all ALS-associated genes in sporadic ALS (sALS) are still scarce. Gene testing in a diagnostic setting is challenging, given the complex genetic architecture of sALS, for which there are genetic variants with large and small effect sizes. Guidelines for the interpretation of genetic variants in gene panels and for counselling of patients are lacking. We aimed to provide a thorough characterization of genetic variability in ALS genes by applying the American College of Medical Genetics and Genomics (ACMG) criteria on whole genome sequencing data from a large cohort of 6013 sporadic ALS patients and 2411 matched controls from Project MinE. We studied genetic variation in 90 ALS-associated genes and applied customized ACMG-criteria to identify pathogenic and likely pathogenic variants. Variants of unknown significance were collected as well. In addition, we determined the length of repeat expansions in C9orf72, ATXN1, ATXN2 and NIPA1 using the ExpansionHunter tool. We found C9orf72 repeat expansions in 5.21% of sALS patients. In 50 ALS-associated genes, we did not identify any pathogenic or likely pathogenic variants. In 5.89%, a pathogenic or likely pathogenic variant was found, most commonly in SOD1, TARDBP, FUS, NEK1, OPTN or TBK1. Significantly more cases carried at least one pathogenic or likely pathogenic variant compared to controls (odds ratio 1.75; P-value 1.64 × 10-5). Isolated risk factors in ATXN1, ATXN2, NIPA1 and/or UNC13A were detected in 17.33% of cases. In 71.83%, we did not find any genetic clues. A combination of variants was found in 2.88%. This study provides an inventory of pathogenic and likely pathogenic genetic variation in a large cohort of sALS patients. Overall, we identified pathogenic and likely pathogenic variants in 11.13% of ALS patients in 38 known ALS genes. In line with the oligogenic hypothesis, we found significantly more combinations of variants in cases compared to controls. Many variants of unknown significance may contribute to ALS risk, but diagnostic algorithms to reliably identify and weigh them are lacking. This work can serve as a resource for counselling and for the assembly of gene panels for ALS. Further characterization of the genetic architecture of sALS is necessary given the growing interest in gene testing in ALS.

Enabling Broader Adoption of Biocatalysis in Organic Chemistry

Biocatalysis is becoming an increasingly impactful method in contemporary synthetic chemistry for target molecule synthesis. The selectivity imparted by enzymes has been leveraged to complete previously intractable chemical transformations and improve synthetic routes toward complex molecules. However, the implementation of biocatalysis in mainstream organic chemistry has been gradual to this point. This is partly due to a set of historical and technological barriers that have prevented chemists from using biocatalysis as a synthetic tool with utility that parallels alternative modes of catalysis. In this Perspective, we discuss these barriers and how they have hindered the adoption of enzyme catalysts into synthetic strategies. We also summarize tools and resources that already enable organic chemists to use biocatalysts. Furthermore, we discuss ways to further lower the barriers for the adoption of biocatalysis by the broader synthetic organic chemistry community through the dissemination of resources, demystifying biocatalytic reactions, and increasing collaboration across the field.