GDF15 is a dynamic biomarker of the integrated stress response in the central nervous system.

AIM: Characterize Growth Differentiation Factor 15 (GDF15) as a secreted biomarker of the integrated stress response (ISR) within the central nervous system (CNS). METHODS: We determined GDF15 levels utilizing in vitro and in vivo neuronal systems wherein the ISR was activated. Primarily, we used the murine model of vanishing white matter disease (VWMD), a neurological disease driven by persistent ISR in the CNS, to establish a link between levels of GDF15 in the cerebrospinal fluid (CSF) and ISR gene expression signature in the CNS. GDF15 was also determined in the CSF of VWM patients. RESULTS: GDF15 expression was increased concomitant to ISR activation in stress-induced primary astrocytes as well as in retinal ganglion cells following optic nerve crush, while treatment with 2Bact, a specific eIF2B activator, suppressed both the ISR and GDF15. In the VWMD model, CSF GDF15 levels corresponded with the magnitude of the ISR and were reduced by 2BAct. In VWM patients, mean CSF GDF15 was elevated >20-fold as compared to healthy controls, whereas plasma GDF15 was undifferentiated. CONCLUSIONS: These data suggest that CSF GDF15 is a dynamic marker of ISR activation in the CNS and may serve as a pharmacodynamic biomarker for ISR-modulating therapies.

PGF2α signaling drives fibrotic remodeling and fibroblast population dynamics in mice.

Idiopathic pulmonary fibrosis (IPF) is a chronic parenchymal lung disease characterized by repetitive alveolar cell injury, myofibroblast proliferation, and excessive extracellular matrix deposition for which unmet need persists for effective therapeutics. The bioactive eicosanoid, prostaglandin F2α, and its cognate receptor FPr (Ptgfr) are implicated as a TGF-ß1-independent signaling hub for IPF. To assess this, we leveraged our published murine PF model (IER-SftpcI73T) expressing a disease-associated missense mutation in the surfactant protein C (Sftpc) gene. Tamoxifen-treated IER-SftpcI73T mice developed an early multiphasic alveolitis and transition to spontaneous fibrotic remodeling by 28 days. IER-SftpcI73T mice crossed to a Ptgfr-null (FPr-/-) line showed attenuated weight loss and gene dosage-dependent rescue of mortality compared with FPr+/+ cohorts. IER-SftpcI73T/FPr-/- mice also showed reductions in multiple fibrotic endpoints for which administration of nintedanib was not additive. Single-cell RNA-Seq, pseudotime analysis, and in vitro assays demonstrated Ptgfr expression predominantly within adventitial fibroblasts, which were reprogrammed to an "inflammatory/transitional" cell state in a PGF2α /FPr-dependent manner. Collectively, the findings provide evidence for a role for PGF2α signaling in IPF, mechanistically identify a susceptible fibroblast subpopulation, and establish a benchmark effect size for disruption of this pathway in mitigating fibrotic lung remodeling.

Temporal Patterns of Behavioral and Psychological Symptoms of Dementia and Caregiver Distress: Associated Daily and Individual Factors.

This study examined daily occurrences of behavioral and psychological symptoms of dementia (BPSD) and whether caregivers' perceived distress towards BPSD varies throughout four phases of the day (i.e., morning, daytime, evening, and night). Family caregivers residing with relatives who were using adult day services (ADS) participated in an 8-day daily diary study (caregiver N = 173; caregiver-day N = 1,359). BPSD occurred most frequently in the evenings. ADS use, sleep disturbances, and dementia severity were significantly associated with BPSD occurrence for some phases of the day. Caregivers' distress towards BPSD occurrences increased throughout the day (i.e., most stressful at night). However, caregivers showed lower reactivity to BPSD at night on days when their relatives used ADS. Evidence of temporal patterns of BPSD in community-dwelling older adults and caregiver distress demonstrated the importance of ADS use for BPSD reactivity and identified potential target windows and associated contextual factors for individualized interventions.

Disruption of Prostaglandin F2α Receptor Signaling Attenuates Fibrotic Remodeling and Alters Fibroblast Population Dynamics in A Preclinical Murine Model of Idiopathic Pulmonary Fibrosis.

Idiopathic Pulmonary Fibrosis (IPF) is a chronic parenchymal lung disease characterized by repetitive alveolar cell injury, myofibroblast proliferation, and excessive extracellular matrix deposition for which unmet need persists for effective therapeutics. The bioactive eicosanoid, prostaglandin F2α, and its cognate receptor FPr (Ptfgr) are implicated as a TGFß1 independent signaling hub for IPF. To assess this, we leveraged our published murine PF model (IER - SftpcI73T) expressing a disease-associated missense mutation in the surfactant protein C (Sftpc) gene. Tamoxifen treated IER-SftpcI73T mice develop an early multiphasic alveolitis and transition to spontaneous fibrotic remodeling by 28 days. IER-SftpcI73T mice crossed to a Ptgfr null (FPr-/-) line showed attenuated weight loss and gene dosage dependent rescue of mortality compared to FPr+/+ cohorts. IER-SftpcI73T/FPr-/- mice also showed reductions in multiple fibrotic endpoints for which administration of nintedanib was not additive. Single cell RNA sequencing, pseudotime analysis, and in vitro assays demonstrated Ptgfr expression predominantly within adventitial fibroblasts which were reprogrammed to an "inflammatory/transitional" cell state in a PGF2α/FPr dependent manner. Collectively, the findings provide evidence for a role for PGF2α signaling in IPF, mechanistically identify a susceptible fibroblast subpopulation, and establish a benchmark effect size for disruption of this pathway in mitigating fibrotic lung remodeling.

Protein kinase R and the integrated stress response drive immunopathology caused by mutations in the RNA deaminase ADAR1.

The RNA deaminase ADAR1 is an essential negative regulator of the RNA sensor MDA5, and loss of ADAR1 function triggers inappropriate activation of MDA5 by self-RNAs. Mutations in ADAR, the gene that encodes ADAR1, cause human immune diseases, including Aicardi-Goutières syndrome (AGS). However, the mechanisms of MDA5-dependent disease pathogenesis in vivo remain unknown. Here we generated mice with a single amino acid change in ADAR1 that models the most common human ADAR AGS mutation. These Adar mutant mice developed lethal disease that required MDA5, the RIG-I-like receptor LGP2, type I interferons, and the eIF2α kinase PKR. A small-molecule inhibitor of the integrated stress response (ISR) that acts downstream of eIF2α phosphorylation prevented immunopathology and rescued the mice from mortality. These findings place PKR and the ISR as central components of immunopathology in vivo and identify therapeutic targets for treatment of human diseases associated with the ADAR1-MDA5 axis.

Comprehensive characterization of RAS mutations in colon and rectal cancers in old and young patients.

Colorectal cancer (CRC) is increasingly appreciated as a heterogeneous disease, with factors such as microsatellite instability (MSI), cancer subsite within the colon versus rectum, and age of diagnosis associated with specific disease course and therapeutic response. Activating oncogenic mutations in KRAS and NRAS are common in CRC, driving tumor progression and influencing efficacy of both cytotoxic and targeted therapies. The RAS mutational spectrum differs substantially between tumors arising from distinct tissues. Structure-function analysis of relatively common somatic RAS mutations in G12, Q61, and other codons is characterized by differing potency and modes of action. Here we show the mutational profile of KRAS, NRAS, and the less common HRAS in 13,336 CRC tumors, comparing the frequency of specific mutations based on age of diagnosis, MSI status, and colon versus rectum subsite. We identify mutation hotspots, and unexpected differences in mutation spectrum, based on these clinical parameters.

Comprehensive Genomic Landscapes in Early and Later Onset Colorectal Cancer.

PURPOSE: The incidence rates of colorectal cancers are increasing in young adults. The objective of this study was to investigate genomic differences between tumor samples collected from younger and older patients with colorectal cancer. EXPERIMENTAL DESIGN: DNA was extracted from 18,218 clinical specimens, followed by hybridization capture of 3,769 exons from 403 cancer-related genes and 47 introns of 19 genes commonly rearranged in cancer. Genomic alterations (GA) were determined, and association with patient age and microsatellite stable/microsatellite instability high (MSS/MSI-H) status established. RESULTS: Overall genomic alteration rates in the younger (<40) and older (≥50) cohorts were similar in the majority of the genes analyzed. Gene alteration rates in the microsatellite stable (MSS) younger and older cohorts were largely similar, with several notable differences. In particular, TP53 (FDR < 0.01) and CTNNB1 (FDR = 0.01) alterations were more common in younger patients with colorectal cancer, and APC (FDR < 0.01), KRAS (FDR < 0.01), BRAF (FDR < 0.01), and FAM123B (FDR < 0.01) were more commonly altered in older patients with colorectal cancer. In the MSI-H cohort, the majority of genes showed similar rate of alterations in all age groups, but with significant differences seen in APC (FDR < 0.01), BRAF (FDR < 0.01), and KRAS (FDR < 0.01). CONCLUSIONS: Tumors from younger and older patients with colorectal cancer demonstrated similar overall rates of genomic alteration. However, differences were noted in several genes relevant to biology and response to therapy. Further study will need to be conducted to determine whether the differences in gene alteration rates can be leveraged to provide personalized therapies for young patients with early-onset sporadic colorectal cancer.

Prevalence of PDL1 Amplification and Preliminary Response to Immune Checkpoint Blockade in Solid Tumors.

Importance: Copy number alterations in programmed cell death ligand 1 (PDL1 or CD274), programmed cell death 1 ligand 2 (PDCD1LG2 or PDL2), and Janus kinase 2 (JAK2) genes (chromosome 9p24.1) characterize Hodgkin lymphoma, resulting in high response rates to programmed cell death 1 (PD-1)/programmed cell death ligand 1 (PD-L1) blockade. The prevalence and utility of PDL1 amplification as a response biomarker to PD-1/PD-L1 blockade are unknown in other tumors. Objectives: To examine the prevalence of PDL1 amplification and its utility as a response biomarker to PD-1/PD-L1 blockade in solid tumors. Design, Setting, and Participants: This retrospective study (October 1, 2012, to October 1, 2017) used a deidentified tumor database from a commercial company and annotated clinical records from a subset of patients treated at a university tertiary referral center. The study analyzed 118 187 tumors from the deidentified database, including a clinically annotated subgroup of 2039 malignant tumors. Interventions: Comprehensive genomic profiling was performed on all samples to determine PDL1 amplification, microsatellite instability, and tumor mutational burden (TMB). A subset of patients was treated with PD-1/PD-L1 blockade. Main Outcomes and Measures: The prevalence of PDL1 amplification was determined among 118 187 patient samples that underwent next-generation sequencing. Solid tumors treated with checkpoint blockade were evaluated for response and progression-free survival (PFS). Results: Of the 118 187 deidentified tumor samples, PDL1 amplifications were identified in 843 (0.7%), including more than 100 types of solid tumors. Most PDL1-amplified tumors (84.8%) had a low to intermediate TMB. PDL1 amplification did not always correlate with high-positive PD-L1 expression by immunohistochemical analysis. Six of 9 patients (66.7%) from 1 center with PDL1-amplified solid tumors had objective responses after checkpoint blockade administration. The median PFS among all treated patients was 15.2 months. Responders included 1 patient with glioblastoma (PFS, ≥5.2 months), 2 patients with head and neck squamous cell cancer (PFS, ≥9 and 15.2 months), 2 patients with metastatic basal cell cancer (PFS, 3.8 and ≥24.1 months), and 1 patient with urothelial cancer (PFS, ≥17.8 months). Conclusions and Relevance: The results of this study suggest that PDL1 amplification occurs in a small subset of malignant tumors. Additional large-scale, prospective studies of PDL1-amplified cancers are warranted to confirm the responses to checkpoint blockade described herein, even in the absence of microsatellite instability, high PD-L1 expression, and a high TMB.

A synthetic biology approach to probing nucleosome symmetry.

The repeating subunit of chromatin, the nucleosome, includes two copies of each of the four core histones, and several recent studies have reported that asymmetrically-modified nucleosomes occur at regulatory elements in vivo. To probe the mechanisms by which histone modifications are read out, we designed an obligate pair of H3 heterodimers, termed H3X and H3Y, which we extensively validated genetically and biochemically. Comparing the effects of asymmetric histone tail point mutants with those of symmetric double mutants revealed that a single methylated H3K36 per nucleosome was sufficient to silence cryptic transcription in vivo. We also demonstrate the utility of this system for analysis of histone modification crosstalk, using mass spectrometry to separately identify modifications on each H3 molecule within asymmetric nucleosomes. The ability to generate asymmetric nucleosomes in vivo and in vitro provides a powerful and generalizable tool to probe the mechanisms by which H3 tails are read out by effector proteins in the cell.

Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden.

BACKGROUND: High tumor mutational burden (TMB) is an emerging biomarker of sensitivity to immune checkpoint inhibitors and has been shown to be more significantly associated with response to PD-1 and PD-L1 blockade immunotherapy than PD-1 or PD-L1 expression, as measured by immunohistochemistry (IHC). The distribution of TMB and the subset of patients with high TMB has not been well characterized in the majority of cancer types. METHODS: In this study, we compare TMB measured by a targeted comprehensive genomic profiling (CGP) assay to TMB measured by exome sequencing and simulate the expected variance in TMB when sequencing less than the whole exome. We then describe the distribution of TMB across a diverse cohort of 100,000 cancer cases and test for association between somatic alterations and TMB in over 100 tumor types. RESULTS: We demonstrate that measurements of TMB from comprehensive genomic profiling are strongly reflective of measurements from whole exome sequencing and model that below 0.5 Mb the variance in measurement increases significantly. We find that a subset of patients exhibits high TMB across almost all types of cancer, including many rare tumor types, and characterize the relationship between high TMB and microsatellite instability status. We find that TMB increases significantly with age, showing a 2.4-fold difference between age 10 and age 90 years. Finally, we investigate the molecular basis of TMB and identify genes and mutations associated with TMB level. We identify a cluster of somatic mutations in the promoter of the gene PMS2, which occur in 10% of skin cancers and are highly associated with increased TMB. CONCLUSIONS: These results show that a CGP assay targeting ~1.1 Mb of coding genome can accurately assess TMB compared with sequencing the whole exome. Using this method, we find that many disease types have a substantial portion of patients with high TMB who might benefit from immunotherapy. Finally, we identify novel, recurrent promoter mutations in PMS2, which may be another example of regulatory mutations contributing to tumorigenesis.

Evolution and genetic architecture of chromatin accessibility and function in yeast.

Chromatin accessibility is an important functional genomics phenotype that influences transcription factor binding and gene expression. Genome-scale technologies allow chromatin accessibility to be mapped with high-resolution, facilitating detailed analyses into the genetic architecture and evolution of chromatin structure within and between species. We performed Formaldehyde-Assisted Isolation of Regulatory Elements sequencing (FAIRE-Seq) to map chromatin accessibility in two parental haploid yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus and their diploid hybrid. We show that although broad-scale characteristics of the chromatin landscape are well conserved between these species, accessibility is significantly different for 947 regions upstream of genes that are enriched for GO terms such as intracellular transport and protein localization exhibit. We also develop new statistical methods to investigate the genetic architecture of variation in chromatin accessibility between species, and find that cis effects are more common and of greater magnitude than trans effects. Interestingly, we find that cis and trans effects at individual genes are often negatively correlated, suggesting widespread compensatory evolution to stabilize levels of chromatin accessibility. Finally, we demonstrate that the relationship between chromatin accessibility and gene expression levels is complex, and a significant proportion of differences in chromatin accessibility might be functionally benign.

Integrative phenomics reveals insight into the structure of phenotypic diversity in budding yeast.

To better understand the quantitative characteristics and structure of phenotypic diversity, we measured over 14,000 transcript, protein, metabolite, and morphological traits in 22 genetically diverse strains of Saccharomyces cerevisiae. More than 50% of all measured traits varied significantly across strains [false discovery rate (FDR) = 5%]. The structure of phenotypic correlations is complex, with 85% of all traits significantly correlated with at least one other phenotype (median = 6, maximum = 328). We show how high-dimensional molecular phenomics data sets can be leveraged to accurately predict phenotypic variation between strains, often with greater precision than afforded by DNA sequence information alone. These results provide new insights into the spectrum and structure of phenotypic diversity and the characteristics influencing the ability to accurately predict phenotypes.

Population genomics and transcriptional consequences of regulatory motif variation in globally diverse Saccharomyces cerevisiae strains.

Noncoding genetic variation is known to significantly influence gene expression levels in a growing number of specific cases; however, the patterns of genome-wide noncoding variation present within populations, the evolutionary forces acting on noncoding variants, and the relative effects of regulatory polymorphisms on transcript abundance are not well characterized. Here, we address these questions by analyzing patterns of regulatory variation in motifs for 177 DNA binding proteins in 37 strains of Saccharomyces cerevisiae. Between S. cerevisiae strains, we found considerable polymorphism in regulatory motifs across strains (mean π = 0.005) as well as diversity in regulatory motifs (mean 0.91 motifs differences per regulatory region). Population genetics analyses reveal that motifs are under purifying selection, and there is considerable heterogeneity in the magnitude of selection across different motifs. Finally, we obtained RNA-Seq data in 22 strains and identified 49 polymorphic DNA sequence motifs in 30 distinct genes that are significantly associated with transcriptional differences between strains. In 22 of these genes, there was a single polymorphic motif associated with expression in the upstream region. Our results provide comprehensive insights into the evolutionary trajectory of regulatory variation in yeast and the characteristics of a compendium of regulatory alleles.

On the prospects of whole-genome association mapping in Saccharomyces cerevisiae.

Advances in sequencing technology have enabled whole-genome sequences to be obtained from multiple individuals within species, particularly in model organisms with compact genomes. For example, 36 genome sequences of Saccharomyces cerevisiae are now publicly available, and SNP data are available for even larger collections of strains. One potential use of these resources is mapping the genetic basis of phenotypic variation through genome-wide association (GWA) studies, with the benefit that associated variants can be studied experimentally with greater ease than in outbred populations such as humans. Here, we evaluate the prospects of GWA studies in S. cerevisiae strains through extensive simulations and a GWA study of mitochondrial copy number. We demonstrate that the complex and heterogeneous patterns of population structure present in yeast populations can lead to a high type I error rate in GWA studies of quantitative traits, and that methods typically used to control for population stratification do not provide adequate control of the type I error rate. Moreover, we show that while GWA studies of quantitative traits in S. cerevisiae may be difficult depending on the particular set of strains studied, association studies to map cis-acting quantitative trait loci (QTL) and Mendelian phenotypes are more feasible. We also discuss sampling strategies that could enable GWA studies in yeast and illustrate the utility of this approach in Saccharomyces paradoxus. Thus, our results provide important practical insights into the design and interpretation of GWA studies in yeast, and other model organisms that possess complex patterns of population structure.

Clusters of adaptive evolution in the human genome.

Considerable work has been devoted to identifying regions of the human genome that have been subjected to recent positive selection. Although detailed follow-up studies of putatively selected regions are critical for a deeper understanding of human evolutionary history, such studies have received comparably less attention. Recently, we have shown that ALMS1 has been the target of recent positive selection acting on standing variation in Eurasian populations. Here, we describe a careful follow-up analysis of genetic variation across the ALMS1 region, which unexpectedly revealed a cluster of substrates of positive selection. Specifically, through the analysis of SNP data from the HapMap and Human Genome Diversity Project-Centre d'Etude du Polymorphisme Humain samples as well sequence data from the region, we find compelling evidence for three independent and distinct signals of recent positive selection across this 3 Mb region surrounding ALMS1. Moreover, we analyzed the HapMap data to identify other putative clusters of independent selective events and conservatively discovered 19 additional clusters of adaptive evolution. This work has important implications for the interpretation of genome-scans for positive selection in humans and more broadly contributes to a better understanding of how recent positive selection has shaped genetic variation across the human genome.

Tracking footprints of artificial selection in the dog genome.

The size, shape, and behavior of the modern domesticated dog has been sculpted by artificial selection for at least 14,000 years. The genetic substrates of selective breeding, however, remain largely unknown. Here, we describe a genome-wide scan for selection in 275 dogs from 10 phenotypically diverse breeds that were genotyped for over 21,000 autosomal SNPs. We identified 155 genomic regions that possess strong signatures of recent selection and contain candidate genes for phenotypes that vary most conspicuously among breeds, including size, coat color and texture, behavior, skeletal morphology, and physiology. In addition, we demonstrate a significant association between HAS2 and skin wrinkling in the Shar-Pei, and provide evidence that regulatory evolution has played a prominent role in the phenotypic diversification of modern dog breeds. Our results provide a first-generation map of selection in the dog, illustrate how such maps can rapidly inform the genetic basis of canine phenotypic variation, and provide a framework for delineating the mechanistic basis of how artificial selection promotes rapid and pronounced phenotypic evolution.

Highly punctuated patterns of population structure on the X chromosome and implications for African evolutionary history.

It is well known that average levels of population structure are higher on the X chromosome compared to autosomes in humans. However, there have been surprisingly few analyses on the spatial distribution of population structure along the X chromosome. With publicly available data from the HapMap Project and Perlegen Sciences, we show a strikingly punctuated pattern of X chromosome population structure. Specifically, 87% of X-linked HapMap SNPs within the top 1% of F(ST) values cluster into five distinct loci. The largest of these regions spans 5.4 Mb and contains 66% of the most highly differentiated HapMap SNPs on the X chromosome. We demonstrate that the extreme clustering of highly differentiated SNPs on the X chromosome is not an artifact of ascertainment bias, nor is it specific to the populations genotyped in the HapMap Project. Rather, additional analyses and resequencing data suggest that these five regions have been substrates of recent and strong adaptive evolution. Finally, we discuss the implications that patterns of X-linked population structure have on the evolutionary history of African populations.

An ortholog of MIXTA-like2 controls epidermal cell shape in flowers of Thalictrum.

Here, we investigated the genetic underpinnings of pollination-related floral phenotypes in Thalictrum, a ranunculid with apetalous flowers. The variable presence of petaloid features in other floral organs correlates with distinct adaptations to insect vs. wind pollination. Conical cells are present in sepals or stamens of insect-pollinated species, and in stigmas. We characterized a Thalictrum ortholog of the Antirrhinum majus transcription factor MIXTA-like2, responsible for conical cells, from three species with distinct floral morphologies, representing two pollination syndromes. Genes were cloned by PCR and analysed phylogenetically. Expression analyses were conducted by quantitative PCR and in situ hybridization, followed by functional studies in transgenic tobacco. The cloned genes encode R2R3 MYB proteins closely related to Antirrhinum AmMYBML2 and Petunia hybrida PhMYB1. Spatial expression by in situ hybridization overlaps areas of conical cells. Overexpression in tobacco induces cell outgrowths in carpel epidermis and significantly increases the height of petal conical cells. We have described the first orthologs of AmMIXTA-like2 outside the core eudicots, likely ancestral to the MIXTA/MIXTA-like1 duplication. The conserved role in epidermal cell elongation results in conical cells, micromorphological markers for petaloidy. This adaptation to attract insect pollinators was apparently lost after the evolution of wind pollination in Thalictrum.