PROBING THE ARCHITECTURE OF COMPLEX TRAITS: FUNCTIONAL GENOMICS METHODS AND APPLICATIONS

Genome wide association studies (GWA) have had tremendous success in identifying genetic variants associated with complex traits. However, the majority of associated loci fall outside of protein coding regions, suggesting a role in regulatory function. This has highlighted a critical need for understanding the regulatory architecture of the genome. Recent advances in high-throughput sequencing technology have enabled transcriptional profiling and mapping of epigenetic features across a broad range of cell types and conditions, both in human and model organisms. As a result, an increasingly higher-resolution genome-wide annotation of regulatory elements is now available. Additionally, expression quantitative trait loci (eQTL) studies mapping the genetic basis of gene expression have identified single nucleotide polymorphisms (SNPs) whose allelic variation correlates with gene expression levels. In conjunction with epigenetic annotations, these results have greatly improved interpretability of variants implicated in complex traits. However, a more comprehensive model of epigenetic regulation in disease can only be obtained by directly assaying disease-relevant tissue in affected individuals. Moreover, traditional eQTL methods often perform a prohibitive number of statistical tests, and are underpowered for detecting weaker associations between SNPs and distally-located genes. In the following chapters I present a novel statistical method that reduces eQTL testing burden and improves power to detect genetic variants associated with expression levels of distal genes. Applying this method to data sets in yeast, mouse, and human, I identified thousands of new eQTL and highlighted candidate master regulators, which were consistently enriched across species for metabolic function. Additionally, I present an analysis of the chromatin and transcriptional landscapes in colon tissue from 33 Crohn’s disease and non-IBD individuals. In ten samples, I found evidence of a molecular signature consistent with metaplasia, the prevalence of which was highly over-represented in CD patients. In an analysis of the remaining individuals, I identified thousands of regulatory regions implicated in disease, many of which co-localize with differentially expressed genes, and highlighted several candidate driver transcription factors. Together, these methods and applications provide a richer understanding of genetic and epigenetic variants implicated in complex traits and disease, and provide hypotheses for future follow up studies.Doctor of Philosoph

Development and assessment of a clinical calculator for estimating the likelihood of recurrence and survival among patients with locally advanced rectal cancer treated with chemotherapy, radiotherapy, and surgery

Importance: Predicting outcomes in patients receiving neoadjuvant therapy for rectal cancer is challenging because of tumor downstaging. Validated clinical calculators that can estimate recurrence-free survival (RFS) and overall survival (OS) among patients with rectal cancer who have received multimodal therapy are needed. Objective: To develop and validate clinical calculators providing estimates of rectal cancer recurrence and survival that are better for individualized decision-making than the American Joint Committee on Cancer (AJCC) staging system or the neoadjuvant rectal (NAR) score. Design, Setting, and Participants: This prognostic study developed risk models, graphically represented as nomograms, for patients with incomplete pathological response using Cox proportional hazards and multivariable regression analyses with restricted cubic splines. Because patients with complete pathological response to neoadjuvant therapy had uniformly favorable outcomes, their predictions were obtained separately. The study included 1400 patients with stage II or III rectal cancer who received treatment with chemotherapy, radiotherapy, and surgery at 2 comprehensive cancer centers (Memorial Sloan Kettering [MSK] Cancer Center and Siteman Cancer Center [SCC]) between January 1, 1998, and December 31, 2017. Patients from the MSK cohort received chemoradiation, surgery, and adjuvant chemotherapy from January 1, 1998, to December 31, 2014; these patients were randomly assigned to either a model training group or an internal validation group. Models were externally validated using data from the SCC cohort, who received either chemoradiation, surgery, and adjuvant chemotherapy (chemoradiotherapy group) or short-course radiotherapy, consolidation chemotherapy, and surgery (total neoadjuvant therapy with short-course radiotherapy group) from January 1, 2009, to December 31, 2017. Data were analyzed from March 1, 2020, to January 10, 2021. Exposures: Chemotherapy, radiotherapy, chemoradiotherapy, and surgery. Main Outcomes and Measures: Recurrence-free survival and OS were the outcome measures, and the discriminatory performance of the clinical calculators was measured with concordance index and calibration plots. The ability of the clinical calculators to predict RFS and OS was compared with that of the AJCC staging system and the NAR score. The models for RFS and OS among patients with incomplete pathological response included postoperative pathological tumor category, number of positive lymph nodes, tumor distance from anal verge, and large- and small-vessel venous and perineural invasion; age was included in the risk model for OS. The final clinical calculators provided RFS and OS estimates derived from Kaplan-Meier curves for patients with complete pathological response and from risk models for patients with incomplete pathological response. Results: Among 1400 total patients with locally advanced rectal cancer, the median age was 57.8 years (range, 18.0-91.9 years), and 863 patients (61.6%) were male, with tumors at a median distance of 6.7 cm (range, 0-15.0 cm) from the anal verge. The MSK cohort comprised 1069 patients; of those, 710 were assigned to the model training group and 359 were assigned to the internal validation group. The SCC cohort comprised 331 patients; of those, 200 were assigned to the chemoradiotherapy group and 131 were assigned to the total neoadjuvant therapy with short-course radiotherapy group. The concordance indices in the MSK validation data set were 0.70 (95% CI, 0.65-0.76) for RFS and 0.73 (95% CI, 0.65-0.80) for OS. In the external SCC data set, the concordance indices in the chemoradiotherapy group were 0.71 (95% CI, 0.62-0.81) for RFS and 0.72 (95% CI, 0.59-0.85) for OS; the concordance indices in the total neoadjuvant therapy with short-course radiotherapy group were 0.62 (95% CI, 0.49-0.75) for RFS and 0.67 (95% CI, 0.46-0.84) for OS. Calibration plots confirmed good agreement between predicted and observed events. These results compared favorably with predictions based on the AJCC staging system (concordance indices for MSK validation: RFS = 0.69 [95% CI, 0.64-0.74]; OS = 0.67 [95% CI, 0.58-0.75]) and the NAR score (concordance indices for MSK validation: RFS = 0.56 [95% CI, 0.50-0.63]; OS = 0.56 [95% CI, 0.46-0.66]). Furthermore, the clinical calculators provided more individualized outcome estimates compared with the categorical schemas (eg, estimated RFS for patients with AJCC stage IIIB disease ranged from 7% to 68%). Conclusions and Relevance: In this prognostic study, clinical calculators were developed and validated; these calculators provided more individualized estimates of the likelihood of RFS and OS than the AJCC staging system or the NAR score among patients with rectal cancer who received multimodal treatment. The calculators were easy to use and applicable to both short- and long-course radiotherapy regimens, and they may be used to inform surveillance strategies and facilitate future clinical trials and statistical power calculations

The early life microbiota protects neonatal mice from pathological small intestinal epithelial cell shedding

The early life gut microbiota plays a crucial role in regulating and maintaining the intestinal barrier, with disturbances in these communities linked to dysregulated renewal and replenishment of intestinal epithelial cells. Here we sought to determine pathological cell shedding outcomes throughout the postnatal developmental period, and which host and microbial factors mediate these responses. Surprisingly, neonatal mice (Day 14 and 21) were highly refractory to induction of cell shedding after intraperitoneal administration of liposaccharide (LPS), with Day 29 mice showing strong pathological responses, more similar to those observed in adult mice. These differential responses were not linked to defects in the cellular mechanisms and pathways known to regulate cell shedding responses. When we profiled microbiota and metabolites, we observed significant alterations. Neonatal mice had high relative abundances of Streptococcus, Escherichia, and Enterococcus and increased primary bile acids. In contrast, older mice were dominated by Candidatus Arthromitus, Alistipes, and Lachnoclostridium, and had increased concentrations of SCFAs and methyamines. Antibiotic treatment of neonates restored LPS-induced small intestinal cell shedding, whereas adult fecal microbiota transplant alone had no effect. Our findings further support the importance of the early life window for microbiota-epithelial interactions in the presence of inflammatory stimuli and highlights areas for further investigation

Conserved Mutations in the Pneumococcal Bacteriocin Transporter Gene, blpA, Result in a Complex Population Consisting of Producers and Cheaters

All fully sequenced strains of Streptococcus pneumoniae possess a version of the blp locus, which is responsible for bacteriocin production and immunity. Activation of the blp locus is stimulated by accumulation of the peptide pheromone, BlpC, following its secretion by the ABC transporter, BlpA. The blp locus is characterized by significant diversity in blpC type and in the region of the locus containing putative bacteriocin and immunity genes. In addition, the blpA gene can represent a single large open reading frame or be divided into several smaller fragments due to the presence of frameshift mutations. In this study, we use a collection of strains with blp-dependent inhibition and immunity to define the genetic changes that bring about phenotypic differences in bacteriocin production or immunity. We demonstrate that alterations in blpA, blpC, and bacteriocin/immunity content likely play an important role in competitive interactions between pneumococcal strains. Importantly, strains with a highly conserved frameshift mutation in blpA are unable to secrete bacteriocins or BlpC, but retain the ability to respond to exogenous peptide pheromone produced by cocolonizing strains, stimulating blp-mediated immunity. These “cheater” strains can only coexist with bacteriocin-producing strains that secrete their cognate BlpC and share the same immunity proteins. The variable outcome of these interactions helps to explain the heterogeneity of the blp pheromone, bacteriocin, and immunity protein content

Discovery and Biosynthesis of Gladiolin: A Burkholderia gladioli Antibiotic with Promising Activity against Mycobacterium tuberculosis.

An antimicrobial activity screen of Burkholderia gladioli BCC0238, a clinical isolate from a cystic fibrosis patient, led to the discovery of gladiolin, a novel macrolide antibiotic with potent activity against Mycobacterium tuberculosis H37Rv. Gladiolin is structurally related to etnangien, a highly unstable antibiotic from Sorangium cellulosum that is also active against Mycobacteria. Like etnangien, gladiolin was found to inhibit RNA polymerase, a validated drug target in M. tuberculosis. However, gladiolin lacks the highly labile hexaene moiety of etnangien and was thus found to possess significantly increased chemical stability. Moreover, gladiolin displayed low mammalian cytotoxicity and good activity against several M. tuberculosis clinical isolates, including four that are resistant to isoniazid and one that is resistant to both isoniazid and rifampicin. Overall, these data suggest that gladiolin may represent a useful starting point for the development of novel drugs to tackle multidrug-resistant tuberculosis. The B. gladioli BCC0238 genome was sequenced using Single Molecule Real Time (SMRT) technology. This resulted in four contiguous sequences: two large circular chromosomes and two smaller putative plasmids. Analysis of the chromosome sequences identified 49 putative specialized metabolite biosynthetic gene clusters. One such gene cluster, located on the smaller of the two chromosomes, encodes a trans-acyltransferase (trans-AT) polyketide synthase (PKS) multienzyme that was hypothesized to assemble gladiolin. Insertional inactivation of a gene in this cluster encoding one of the PKS subunits abrogated gladiolin production, confirming that the gene cluster is responsible for biosynthesis of the antibiotic. Comparison of the PKSs responsible for the assembly of gladiolin and etnangien showed that they possess a remarkably similar architecture, obfuscating the biosynthetic mechanisms responsible for most of the structural differences between the two metabolites

Molecular classification of Crohn's disease reveals two clinically relevant subtypes

The clinical presentation and course of Crohn’s disease (CD) is highly variable. We sought to better understand the cellular and molecular mechanisms that guide this heterogeneity, and characterize the cellular processes associated with disease phenotypes

Alterations to chromatin in intestinal macrophages link IL-10 deficiency to inappropriate inflammatory responses

Intestinal macrophages are uniquely programmed to tolerate exposure to bacteria without mounting potent inflammatory responses. The cytokine IL-10 maintains the macrophage anti-inflammatory response such that loss of IL-10 results in chronic intestinal inflammation. To investigate how IL-10-deficiency alters intestinal macrophage programming and bacterial tolerance, we studied changes in chromatin accessibility in response to bacteria in macrophages from two distinct niches, the intestine and bone-marrow, from both wild-type and IL-10-deficient mice. In both bone-marrow-derived and intestinal macrophages, we identified chromatin accessibility changes associated with bacterial exposure and IL-10-deficiency. Surprisingly, IL-10-deficient intestinal macrophages adopted chromatin and gene expression patterns characteristic of an inflammatory response, even in the absence of bacteria. Further, if IL-10 protein was added to cells that had previously been IL-10-deficient, it could not revert the chromatin landscape to a normal state. Our results demonstrate that IL-10 deficiency results in stable chromatin alterations in macrophages, even in the absence of bacteria. This supports a model where IL-10-deficiency leads to chromatin alterations that contribute to a loss of intestinal macrophage tolerance to bacteria, which is a primary initiating event in chronic intestinal inflammation

Conversion and Colonial History in Icíar Bollaín’s También la lluvia (2010)

This study concerns the representation of colonial Latin American history and the characterisation of Daniel/Hatuey in the 2010 film-about-a-film También la lluvia. A metacinematic work comprising historical study and political commentary, También la lluvia has received mixed critical reactions regarding its portrayal of the historical and social inequalities it analyses. This article examines the ambiguous nature of the work by analysing the motif of conversion. It argues that, by foregrounding the contemporary conversion story of Costa, the film sacrifices both nuanced historical attention to the colonial past it dramatises and sustained development of one of its apparently central characters: Daniel/Hatuey, who is repeatedly converted into narrative and symbolic figures of secondary prominence, despite their importance to the development and legibility of the work as a whole

DUSP5-mediated inhibition of smooth muscle cell proliferation suppresses pulmonary hypertension and right ventricular hypertrophy

Pulmonary hypertension (PH) is associated with structural remodeling of pulmonary arteries (PAs) because of excessive proliferation of fibroblasts, endothelial cells, and smooth muscle cells (SMCs). The peptide hormone angiotensin II (ANG II) contributes to pulmonary vascular remodeling, in part, through its ability to trigger extracellular signal-regulated kinase (ERK1/2) activation. Here, we demonstrate that the ERK1/2 phosphatase, dual-specificity phosphatase 5 (DUSP5), functions as a negative regulator of ANG II-mediated SMC proliferation and PH. In contrast to wild-type controls, Dusp5 null mice infused with ANG II developed PH and right ventricular (RV) hypertrophy. PH in Dusp5 null mice was associated with thickening of the medial layer of small PAs, suggesting an in vivo role for DUSP5 as a negative regulator of ANG II-dependent SMC proliferation. Consistent with this, overexpression of DUSP5 blocked ANG II-mediated proliferation of cultured human pulmonary artery SMCs (hPASMCs) derived from patients with idiopathic PH or from failed donor controls. Collectively, the data support a role for DUSP5 as a feedback inhibitor of ANG II-mediated ERK signaling and PASMC proliferation and suggest that disruption of this circuit leads to adverse cardiopulmonary remodeling. NEW & NOTEWORTHY Dual-specificity phosphatases (DUSPs) serve critical roles in the regulation of mitogen-activated protein kinases, but their functions in the cardiovascular system remain poorly defined. Here, we provide evidence that DUSP5, which resides in the nucleus and specifically dephosphorylates extracellular signal-regulated kinase (ERK1/2), blocks pulmonary vascular smooth muscle cell proliferation. In response to angiotensin II infusion, mice lacking DUSP5 develop pulmonary hypertension and right ventricular cardiac hypertrophy. These findings illustrate DUSP5-mediated suppression of ERK signaling in the lungs as a protective mechanism