Garima Yagnik

A genome-wide association study identifies susceptibility loci for nonsyndromic sagittal craniosynostosis near BMP2 and within BBS9

Sagittal craniosynostosis is the most common form of craniosynostosis, affecting approximately one in 5,000 newborns. We conducted, to our knowledge, the first genome-wide association study for nonsyndromic sagittal craniosynostosis (sNSC) using 130 non-Hispanic case-parent trios of European ancestry (NHW). We found robust associations in a 120-kb region downstream of BMP2 flanked by rs1884302 (P = 1.13 × 10−14, odds ratio (OR) = 4.58) and rs6140226 (P = 3.40 × 10−11, OR = 0.24) and within a 167-kb region of BBS9 between rs10262453 (P = 1.61 × 10−10, OR = 0.19) and rs17724206 (P = 1.50 × 10−8, OR = 0.22). We replicated the associations to both loci (rs1884302, P = 4.39 × 10−31 and rs10262453, P = 3.50 × 10−14) in an independent NHW population of 172 unrelated probands with sNSC and 548 controls. Both BMP2 and BBS9 are genes with roles in skeletal development that warrant functional studies to further understand the etiology of sNSC.

Heterozygous mutations of FREM1 are associated with an increased risk of isolated metopic craniosynostosis in humans and mice

The premature fusion of the paired frontal bones results in metopic craniosynostosis (MC) and gives rise to the clinical phenotype of trigonocephaly. Deletions of chromosome 9p22.3 are well described as a cause of MC with variably penetrant midface hypoplasia. In order to identify the gene responsible for the trigonocephaly component of the 9p22.3 syndrome, a cohort of 109 patients were assessed by high-resolution arrays and MLPA for copy number variations (CNVs) involving 9p22. Five CNVs involving FREM1, all of which were de novo variants, were identified by array-based analyses. The remaining 104 patients with MC were then subjected to targeted FREM1 gene re-sequencing, which identified 3 further mutant alleles, one of which was de novo. Consistent with a pathogenic role, mouse Frem1 mRNA and protein expression was demonstrated in the metopic suture as well as in the pericranium and dura mater. Micro-computed tomography based analyses of the mouse posterior frontal (PF) suture, the human metopic suture equivalent, revealed advanced fusion in all mice homozygous for either of two different Frem1 mutant alleles, while heterozygotes exhibited variably penetrant PF suture anomalies. Gene dosage-related penetrance of midfacial hypoplasia was also evident in the Frem1 mutants. These data suggest that CNVs and mutations involving FREM1 can be identified in a significant percentage of people with MC with or without midface hypoplasia. Furthermore, we present Frem1 mutant mice as the first bona fide mouse model of human metopic craniosynostosis and a new model for midfacial hypoplasia.

ALX4 gain-of-function mutations in nonsyndromic craniosynostosis

Craniosynostosis is the early fusion of one or more sutures of the infant skull and is a common defect occurring in approximately 1 of every 2,500 live births. Nonsyndromic craniosynostosis (NSC) accounts for approximately 80% of all cases and is thought to have strong genetic determinants that are yet to be identified. ALX4 is a homeodomain transcription factor with known involvement in osteoblast regulation. By direct sequencing of the ALX4 coding region in sagittal or sagittal‐suture‐involved nonsyndromic craniosynostosis probands, we identified novel, nonsynonymous, familial variants in three of 203 individuals with NSC. Using dual‐luciferase assay we show that two of these variants (V7F and K211E) confer a significant gain‐of‐function effect on ALX4. Our results suggest that ALX4 variants may have an impact on the genetic etiology of NSC. Hum Mutat 33:1626–1629, 2012.

Genome-wide expression profiling of urinary bladder implicates desmosomal and cytoskeletal dysregulation in the bladder exstrophy-epispadias complex

The bladder exstrophy-epispadias complex (BEEC) represents a spectrum of urological abnormalities where part, or all, of the distal urinary tract fails to close during development, becoming exposed on the outer abdominal wall. While the etiology of BEEC remains unknown, strong evidence exists that genetic factors are implicated. To understand the pathways regulating embryonic bladder development and to identify high-probability BEEC candidate genes, we conducted a genome-wide expression profiling (GWEP) study using normal and exstrophic human urinary bladders and human and mouse embryologic bladder-precursor tissues. We identified 162 genes differentially expressed in both embryonic and postnatal human samples. Pathway analysis of these genes revealed 11 biological networks with top functions related to skeletal and muscular system development, cellular assembly and development, organ morphology, or connective tissue disorders. The two most down-regulated genes desmin (DES, fold-change, -74.7) and desmuslin (DMN, fold-change, -53.0) are involved in desmosome mediated cell-cell adhesion and cytoskeletal architecture. Intriguingly, the sixth most overexpressed gene was desmoplakin (DSP, fold-change, +48.8), the most abundant desmosomal protein. We found 30% of the candidate genes to be directly associated with desmosome structure/function or cytoskeletal assembly, pointing to desmosomal and/or cytoskeletal deregulation as an etiologic factor for BEEC. Further findings indicate that p63, PERP, SYNPO2 and the Wnt pathway may also contribute to BEEC etiology. This study provides the first expression profile of urogenital genes during bladder development and points to the high-probability candidate genes for BEEC.

Single-cell profiling of human gliomas reveals macrophage ontogeny as a basis for regional differences in macrophage activation in the tumor microenvironment

BackgroundTumor-associated macrophages (TAMs) are abundant in gliomas and immunosuppressive TAMs are a barrier to emerging immunotherapies. It is unknown to what extent macrophages derived from peripheral blood adopt the phenotype of brain-resident microglia in pre-treatment gliomas. The relative proportions of blood-derived macrophages and microglia have been poorly quantified in clinical samples due to a paucity of markers that distinguish these cell types in malignant tissue.ResultsWe perform single-cell RNA-sequencing of human gliomas and identify phenotypic differences in TAMs of distinct lineages. We isolate TAMs from patient biopsies and compare them with macrophages from non-malignant human tissue, glioma atlases, and murine glioma models. We present a novel signature that distinguishes TAMs by ontogeny in human gliomas. Blood-derived TAMs upregulate immunosuppressive cytokines and show an altered metabolism compared to microglial TAMs. They are also enriched in perivascular and necrotic regions. The gene signature of blood-derived TAMs, but not microglial TAMs, correlates with significantly inferior survival in low-grade glioma. Surprisingly, TAMs frequently co-express canonical pro-inflammatory (M1) and alternatively activated (M2) genes in individual cells.ConclusionsWe conclude that blood-derived TAMs significantly infiltrate pre-treatment gliomas, to a degree that varies by glioma subtype and tumor compartment. Blood-derived TAMs do not universally conform to the phenotype of microglia, but preferentially express immunosuppressive cytokines and show an altered metabolism. Our results argue against status quo therapeutic strategies that target TAMs indiscriminately and in favor of strategies that specifically target immunosuppressive blood-derived TAMs.

Macrophage migration inhibitory factor downregulation: a novel mechanism of resistance to anti-angiogenic therapy

Anti-angiogenic therapies for cancer such as VEGF neutralizing antibody bevacizumab have limited durability. While mechanisms of resistance remain undefined, it is likely that acquired resistance to anti-angiogenic therapy will involve alterations of the tumor microenvironment. We confirmed increased tumor-associated macrophages in bevacizumab-resistant glioblastoma patient specimens and two novel glioblastoma xenograft models of bevacizumab resistance. Microarray analysis suggested downregulated macrophage migration inhibitory factor (MIF) to be the most pertinent mediator of increased macrophages. Bevacizumab-resistant patient glioblastomas and both novel xenograft models of resistance had less MIF than bevacizumab-naive tumors, and harbored more M2/protumoral macrophages that specifically localized to the tumor edge. Xenografts expressing MIF-shRNA grew more rapidly with greater angiogenesis and had macrophages localizing to the tumor edge which were more prevalent and proliferative, and displayed M2 polarization, whereas bevacizumab-resistant xenografts transduced to upregulate MIF exhibited the opposite changes. Bone marrow-derived macrophage were polarized to an M2 phenotype in the presence of condition-media derived from bevacizumab-resistant xenograft-derived cells, while recombinant MIF drove M1 polarization. Media from macrophages exposed to bevacizumab-resistant tumor cell conditioned media increased glioma cell proliferation compared with media from macrophages exposed to bevacizumab-responsive tumor cell media, suggesting that macrophage polarization in bevacizumab-resistant xenografts is the source of their aggressive biology and results from a secreted factor. Two mechanisms of bevacizumab-induced MIF reduction were identified: (1) bevacizumab bound MIF and blocked MIF-induced M1 polarization of macrophages; and (2) VEGF increased glioma MIF production in a VEGFR2-dependent manner, suggesting that bevacizumab-induced VEGF depletion would downregulate MIF. Site-directed biopsies revealed enriched MIF and VEGF at the enhancing edge in bevacizumab-naive patients. This MIF enrichment was lost in bevacizumab-resistant glioblastomas, driving a tumor edge M1-to-M2 transition. Thus, bevacizumab resistance is driven by reduced MIF at the tumor edge causing proliferative expansion of M2 macrophages, which in turn promotes tumor growth.

Leucine-rich repeat, immunoglobulin-like and transmembrane domain 3 (LRIT3) is a modulator of FGFR1

Fibroblast growth factor receptors (FGFRs) play critical roles in craniofacial and skeletal development via multiple signaling pathways including MAPK, PI3K/AKT, and PLC‐γ. FGFR‐mediated signaling is modulated by several regulators. Proteins with leucine‐rich repeat (LRR) and/or immunoglobulin (IG) superfamily domains have been suggested to interact with FGFRs. In addition, fibronectin leucine‐rich repeat transmembrane protein 3 (FLRT3) has been shown to modulate the FGFR‐mediated signaling via the fibronectin type III (FNIII) domain. Therefore proteins with LRR, IG, and FNIII are candidate regulators of the FGFRs. Here we identify leucine‐rich repeat, immunoglobulin‐like and transmembrane domain 3 (LRIT3) as a regulator of the FGFRs.

GLUT3 upregulation promotes metabolic reprogramming associated with antiangiogenic therapy resistance

Clinical trials revealed limited response duration of glioblastomas to VEGF-neutralizing antibody bevacizumab. Thriving in the devascularized microenvironment occurring after antiangiogenic therapy requires tumor cell adaptation to decreased glucose, with 50% less glucose identified in bevacizumab-treated xenografts. Compared with bevacizumab-responsive xenograft cells, resistant cells exhibited increased glucose uptake, glycolysis, 13C NMR pyruvate to lactate conversion, and survival in low glucose. Glucose transporter 3 (GLUT3) was upregulated in bevacizumab-resistant versus sensitive xenografts and patient specimens in a HIF-1α-dependent manner. Resistant versus sensitive cell mitochondria in oxidative phosphorylation-selective conditions produced less ATP. Despite unchanged mitochondrial numbers, normoxic resistant cells had lower mitochondrial membrane potential than sensitive cells, confirming poorer mitochondrial health, but avoided the mitochondrial dysfunction of hypoxic sensitive cells. Thin-layer chromatography revealed increased triglycerides in bevacizumab-resistant versus sensitive xenografts, a change driven by mitochondrial stress. A glycogen synthase kinase-3β inhibitor suppressing GLUT3 transcription caused greater cell death in bevacizumab-resistant than -responsive cells. Overexpressing GLUT3 in tumor cells recapitulated bevacizumab-resistant cell features: survival and proliferation in low glucose, increased glycolysis, impaired oxidative phosphorylation, and rapid in vivo proliferation only slowed by bevacizumab to that of untreated bevacizumab-responsive tumors. Targeting GLUT3 or the increased glycolysis reliance in resistant tumors could unlock the potential of antiangiogenic treatments.

Candidate gene association study implicates p63 in the etiology of nonsyndromic bladder-exstrophy-epispadias complex

BACKGROUND Bladder-exstrophy-epispadias complex (BEEC) is a severe congenital anomaly that represents a spectrum of urological abnormalities where parts or all of the distal urinary tract fail to close during development. Multiple lines of evidence strongly suggested p63 as a plausible candidate gene. We conducted a candidate gene association study to further investigate the role of p63 in human BEEC. METHODS We conducted a family-based association study of p63 using 154 Caucasian patients with nonsyndromic BEEC and their unaffected parents. High throughput single nucleotide polymorphism (SNP) genotyping was carried out using Illuminas Golden Gate Assay for 109 selected tagging SNPs localized within p63 with a minor allele frequency > 0.01. Individual and haplotype SNP transmission disequilibrium tests were conducted using Plink and Haploview, respectively. We also examined parent-of-origin effects using paternal asymmetry tests implemented in FAMHAP (http://famhap.meb.uni-bonn.de/index.html). RESULTS Nominally significant associations were identified between BEEC and six SNPs (rs17447782, rs1913720, rs6790167, rs9865857, rs1543969, rs4687100), and four haplotype blocks including or near these significant SNPs. Analysis of parent-of-origin effects showed significant results for seven SNPs (rs4118375, rs12696596, rs6779677, rs13091309, rs7642420, rs1913721, and rs1399774). None of these results remained significant after multiple testing correction. CONCLUSION The altered transmission of p63 variants in nonsyndromic BEEC patients may be suggestive of its involvement in the disease etiology. Further and large multi-institutional collaborative studies are required to elucidate the role of p63 in nonsyndromic BEEC.

Cross-activating c-Met/β1 integrin complex drives metastasis and invasive resistance in cancer

Significance Invasion is a major cause of cancer mortality, as exemplified by metastatic spread of peripheral malignancies or local intracranial invasion of glioblastoma. While individual mediators of invasion are identified, functional or structural interactions between these mediators remain undefined. We identified a structural cross-activating c-Met/β1 integrin complex that promotes breast cancer metastases and invasive resistance of glioblastoma to the antiangiogenic therapy bevacizumab. We show that tumor cells adapt to their microenvironmental stressors by usurping c-Met and β1 integrin, with c-Met displacing α5 integrin from β1 integrin to form a c-Met/β1 complex with far greater fibronectin affinity than α5β1 integrin. These findings challenge conventional thinking about integrin–ligand interactions and define a molecular target for disrupting metastases or invasive oncologic resistance. The molecular underpinnings of invasion, a hallmark of cancer, have been defined in terms of individual mediators but crucial interactions between these mediators remain undefined. In xenograft models and patient specimens, we identified a c-Met/β1 integrin complex that formed during significant invasive oncologic processes: breast cancer metastases and glioblastoma invasive resistance to antiangiogenic VEGF neutralizing antibody, bevacizumab. Inducing c-Met/β1 complex formation through an engineered inducible heterodimerization system promoted features crucial to overcoming stressors during metastases or antiangiogenic therapy: migration in the primary site, survival under hypoxia, and extravasation out of circulation. c-Met/β1 complex formation was up-regulated by hypoxia, while VEGF binding VEGFR2 sequestered c-Met and β1 integrin, preventing their binding. Complex formation promoted ligand-independent receptor activation, with integrin-linked kinase phosphorylating c-Met and crystallography revealing the c-Met/β1 complex to maintain the high-affinity β1 integrin conformation. Site-directed mutagenesis verified the necessity for c-Met/β1 binding of amino acids predicted by crystallography to mediate their extracellular interaction. Far-Western blotting and sequential immunoprecipitation revealed that c-Met displaced α5 integrin from β1 integrin, creating a complex with much greater affinity for fibronectin (FN) than α5β1. Thus, tumor cells adapt to microenvironmental stressors induced by metastases or bevacizumab by coopting receptors, which normally promote both cell migration modes: chemotaxis, movement toward concentrations of environmental chemoattractants, and haptotaxis, movement controlled by the relative strengths of peripheral adhesions. Tumor cells then redirect these receptors away from their conventional binding partners, forming a powerful structural c-Met/β1 complex whose ligand-independent cross-activation and robust affinity for FN drive invasive oncologic processes.