👤 Rebecca Smith

The challenge of developing effective pharmacodynamic biomarkers for preclinical and clinical testing of FGFR signaling inhibition is significant. Assays that rely on the measurement of phospho-protein epitopes can be limited by the availability of effective antibody detection reagents. Transcript profiling enables accurate quantification of many biomarkers and provides a broader representation of pathway modulation. To identify dynamic transcript biomarkers of FGFR signaling inhibition by AZD4547, a potent inhibitor of FGF receptors 1, 2, and 3, a gene expression profiling study was performed in FGFR2-amplified, drug-sensitive tumor cell lines. Consistent with known signaling pathways activated by FGFR, we identified transcript biomarkers downstream of the RAS-MAPK and PI3K/AKT pathways. Using different tumor cell lines in vitro and xenografts in vivo, we confirmed that some of these transcript biomarkers (DUSP6, ETV5, YPEL2) were modulated downstream of oncogenic FGFR1, 2, 3, whereas others showed selective modulation only by FGFR2 signaling (EGR1). These transcripts showed consistent time-dependent modulation, corresponding to the plasma exposure of AZD4547 and inhibition of phosphorylation of the downstream signaling molecules FRS2 or ERK. Combination of FGFR and AKT inhibition in an FGFR2-mutated endometrial cancer xenograft model enhanced modulation of transcript biomarkers from the PI3K/AKT pathway and tumor growth inhibition. These biomarkers were detected on the clinically validated nanoString platform. Taken together, these data identified novel dynamic transcript biomarkers of FGFR inhibition that were validated in a number of in vivo models, and which are more robustly modulated by FGFR inhibition than some conventional downstream signaling protein biomarkers. Mol Cancer Ther; 15(11); 2802-13. ©2016 AACR. Show less

Delta-5 desaturase (D5D) and delta-6 desaturase (D6D), encoded by fatty acid desaturase 1 (FADS1) and FADS2 genes, respectively, are enzymes in the synthetic pathways for ω3, ω6, and ω9 polyunsaturated fatty acids (PUFAs). Although PUFAs appear to be involved in mammalian metabolic pathways, the physiologic effect of isolated D5D deficiency on these pathways is unclear. After generating >4,650 knockouts (KOs) of independent mouse genes and analyzing them in our high-throughput phenotypic screen, we found that Fads1 KO mice were among the leanest of 3,651 chow-fed KO lines analyzed for body composition and were among the most glucose tolerant of 2,489 high-fat-diet-fed KO lines analyzed by oral glucose tolerance test. In confirmatory studies, chow- or high-fat-diet-fed Fads1 KO mice were leaner than wild-type (WT) littermates; when data from multiple cohorts of adult mice were combined, body fat was 38% and 31% lower in Fads1 male and female KO mice, respectively. Fads1 KO mice also had lower glucose and insulin excursions during oral glucose tolerance tests along with lower fasting glucose, insulin, triglyceride, and total cholesterol levels. In additional studies using a vascular injury model, Fads1 KO mice had significantly decreased femoral artery intima/media ratios consistent with a decreased inflammatory response in their arterial wall. Based on this result, we bred Fads1 KO and WT mice onto an ApoE KO background and fed them a Western diet for 14 weeks; in this atherogenic environment, aortic trees of Fads1 KO mice had 40% less atheromatous plaque compared to WT littermates. Importantly, PUFA levels measured in brain and liver phospholipid fractions of Fads1 KO mice were consistent with decreased D5D activity and normal D6D activity. The beneficial metabolic phenotype demonstrated in Fads1 KO mice suggests that selective D5D inhibitors may be useful in the treatment of human obesity, diabetes, and atherosclerotic cardiovascular disease. Show less

Lethal fetal akinesia deformation sequence (FADS) describes a clinically and genetically heterogeneous phenotype that includes fetal akinesia, intrauterine growth retardation, arthrogryposis and developmental anomalies. Affected babies die as a result of pulmonary hypoplasia. We aimed to identify the underlying genetic cause of this disorder in a family in which there were three affected individuals from two sibships. Autosomal-recessive inheritance was suggested by a family history of consanguinity and by recurrence of the phenotype between the two sibships. We performed exome sequencing of the affected individuals and their unaffected mother, followed by autozygosity mapping and variant filtering to identify the causative gene. Five autozygous regions were identified, spanning 31.7 Mb of genomic sequence and including 211 genes. Using standard variant filtering criteria, we excluded all variants as being the likely pathogenic cause, apart from a single novel nonsense mutation, c.188C>A p.(Ser63*) (NM₀₀₂₄₇₈.4), in MYOD1. This gene encodes an extensively studied transcription factor involved in muscle development, which has nonetheless not hitherto been associated with a hereditary human disease phenotype. We provide the first description of a human phenotype that appears to result from MYOD1 mutation. The presentation with FADS is consistent with a large body of data demonstrating that in the mouse, MyoD is a major controller of precursor cell commitment to the myogenic differentiation programme. Show less

White blood cells play diverse roles in innate and adaptive immunity. Genetic association analyses of phenotypic variation in circulating white blood cell (WBC) counts from large samples of otherwise healthy individuals can provide insights into genes and biologic pathways involved in production, differentiation, or clearance of particular WBC lineages (myeloid, lymphoid) and also potentially inform the genetic basis of autoimmune, allergic, and blood diseases. We performed an exome array-based meta-analysis of total WBC and subtype counts (neutrophils, monocytes, lymphocytes, basophils, and eosinophils) in a multi-ancestry discovery and replication sample of ∼157,622 individuals from 25 studies. We identified 16 common variants (8 of which were coding variants) associated with one or more WBC traits, the majority of which are pleiotropically associated with autoimmune diseases. Based on functional annotation, these loci included genes encoding surface markers of myeloid, lymphoid, or hematopoietic stem cell differentiation (CD69, CD33, CD87), transcription factors regulating lineage specification during hematopoiesis (ASXL1, IRF8, IKZF1, JMJD1C, ETS2-PSMG1), and molecules involved in neutrophil clearance/apoptosis (C10orf54, LTA), adhesion (TNXB), or centrosome and microtubule structure/function (KIF9, TUBD1). Together with recent reports of somatic ASXL1 mutations among individuals with idiopathic cytopenias or clonal hematopoiesis of undetermined significance, the identification of a common regulatory 3' UTR variant of ASXL1 suggests that both germline and somatic ASXL1 mutations contribute to lower blood counts in otherwise asymptomatic individuals. These association results shed light on genetic mechanisms that regulate circulating WBC counts and suggest a prominent shared genetic architecture with inflammatory and autoimmune diseases. Show less

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79 x 10(-13)), rs74506613 (JMJD1C, P = 1.17 x 10(-19)), rs4782371 (ZFPM1, P = 1.59 x 10(-9)) and rs2639990 (ZADH2, P = 1.72 x 10(-8)), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52 x 10(-18); rs7043199, VLDLR-AS1, P = 5.12 x 10(-14)) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39 x 10(-1467); rs1740073, C6orf223, P = 2.34 x 10(-17); rs6993770, ZFPM2, P = 2.44 x 10(-60); rs2375981, KCNV2, P = 1.48 x 10(-100)). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases. Show less

APOE ɛ4, the most significant genetic risk factor for Alzheimer disease (AD), may mask effects of other loci. We re-analyzed genome-wide association study (GWAS) data from the International Genomics of Alzheimer's Project (IGAP) Consortium in APOE ɛ4+ (10 352 cases and 9207 controls) and APOE ɛ4- (7184 cases and 26 968 controls) subgroups as well as in the total sample testing for interaction between a single-nucleotide polymorphism (SNP) and APOE ɛ4 status. Suggestive associations (P<1 × 10(-4)) in stage 1 were evaluated in an independent sample (stage 2) containing 4203 subjects (APOE ɛ4+: 1250 cases and 536 controls; APOE ɛ4-: 718 cases and 1699 controls). Among APOE ɛ4- subjects, novel genome-wide significant (GWS) association was observed with 17 SNPs (all between KANSL1 and LRRC37A on chromosome 17 near MAPT) in a meta-analysis of the stage 1 and stage 2 data sets (best SNP, rs2732703, P=5·8 × 10(-9)). Conditional analysis revealed that rs2732703 accounted for association signals in the entire 100-kilobase region that includes MAPT. Except for previously identified AD loci showing stronger association in APOE ɛ4+ subjects (CR1 and CLU) or APOE ɛ4- subjects (MS4A6A/MS4A4A/MS4A6E), no other SNPs were significantly associated with AD in a specific APOE genotype subgroup. In addition, the finding in the stage 1 sample that AD risk is significantly influenced by the interaction of APOE with rs1595014 in TMEM106B (P=1·6 × 10(-7)) is noteworthy, because TMEM106B variants have previously been associated with risk of frontotemporal dementia. Expression quantitative trait locus analysis revealed that rs113986870, one of the GWS SNPs near rs2732703, is significantly associated with four KANSL1 probes that target transcription of the first translated exon and an untranslated exon in hippocampus (P ⩽ 1.3 × 10(-8)), frontal cortex (P ⩽ 1.3 × 10(-9)) and temporal cortex (P⩽1.2 × 10(-11)). Rs113986870 is also strongly associated with a MAPT probe that targets transcription of alternatively spliced exon 3 in frontal cortex (P=9.2 × 10(-6)) and temporal cortex (P=2.6 × 10(-6)). Our APOE-stratified GWAS is the first to show GWS association for AD with SNPs in the chromosome 17q21.31 region. Replication of this finding in independent samples is needed to verify that SNPs in this region have significantly stronger effects on AD risk in persons lacking APOE ɛ4 compared with persons carrying this allele, and if this is found to hold, further examination of this region and studies aimed at deciphering the mechanism(s) are warranted. Show less

Though current functional genomic screening systems are useful for investigating human susceptibility to chemical toxicity, they have limitations. Well-established, high-throughput yeast mutant screens identify only evolutionarily conserved processes. RNA interference can be applied in human cells but is limited by incomplete gene knockout and off-target effects. Human haploid cell screening is advantageous as it requires knockdown of only a single copy of each gene. A human haploid cell mutant library (KBM7-Mu), derived from a chronic myeloid leukemia (CML) patient, was recently developed and has been used to identify genes that modulate sensitivity to infectious agents and pharmaceutical drugs. Here, we sought to improve the KBM7-Mu screening process to enable efficient screening of environmental chemicals. We developed a semi-solid medium based screening approach that cultures individual mutant colonies from chemically resistant cells, faster (by 2-3 weeks) and with less labor than the original liquid medium-based approach. As proof of principle, we identified genetic mutants that confer resistance to the carcinogen formaldehyde (FA, 12 genes, 18 hits) and the CML chemotherapeutic agent imatinib (6 genes, 13 hits). Validation experiments conducted on KBM7 mutants lacking each of the 18 genes confirmed resistance of 6 FA mutants (CTC1, FCRLA, GOT1, LPR5, M1AP, and MAP2K5) and 1 imatinib-resistant mutant (LYRM9). Despite the improvements to the method, it remains technically challenging to limit false positive findings. Nonetheless, our findings demonstrate the broad applicability of this optimized haploid approach to screen toxic chemicals to identify novel susceptibility genes and gain insight into potential mechanisms of toxicity. Show less

Mutations of cardiac sarcomere genes have been identified to cause HCM, but the molecular mechanisms that lead to cardiomyocyte hypertrophy and risk for sudden death are uncertain. The aim of this study was to examine HCM disease mechanisms at play during cardiac differentiation of human HCM specific pluripotent stem cells. We generated a human embryonic stem cell (hESC) line carrying a naturally occurring mutation of MYPBC3 (c.2905 +1 G >A) to study HCM pathogenesis during cardiac differentiation. HCM-specific hESC-derived cardiomyocytes (hESC-CMs) displayed hallmark aspects of HCM including sarcomere disarray, hypertrophy and impaired calcium impulse propagation. HCM hESC-CMs presented a transient haploinsufficiency of cMyBP-C during cardiomyocyte differentiation, but by day 30 post-differentiation cMyBP-C levels were similar to control hESC-CMs. Gene transfer of full-length MYBPC3 during differentiation prevented hypertrophy, sarcomere disarray and improved calcium impulse propagation in HCM hESC-CMs. These findings point to the critical role of MYBPC3 during sarcomere assembly in cardiac myocyte differentiation and suggest developmental influences of MYBPC3 truncating mutations on the mature hypertrophic phenotype. Show less

Cardiac myosin binding protein-C (cMyBP-C) is a structural and regulatory component of cardiac thick filaments. It is observed in electron micrographs as seven to nine transverse stripes in the central portion of each half of the A band. Its C-terminus binds tightly to the myosin rod and contributes to thick filament structure, while the N-terminus can bind both myosin S2 and actin, influencing their structure and function. Mutations in the MYBPC3 gene (encoding cMyBP-C) are commonly associated with hypertrophic cardiomyopathy (HCM). In cardiac cells there exists a population of myosin heads in the super-relaxed (SRX) state, which are bound to the thick filament core with a highly inhibited ATPase activity. This report examines the role cMyBP-C plays in regulating the population of the SRX state of cardiac myosin by using an assay that measures single ATP turnover of myosin. We report a significant decrease in the proportion of myosin heads in the SRX state in homozygous cMyBP-C knockout mice, however heterozygous cMyBP-C knockout mice do not significantly differ from the wild type. A smaller, non-significant decrease is observed when thoracic aortic constriction is used to induce cardiac hypertrophy in mutation negative mice. These results support the proposal that cMyBP-C stabilises the thick filament and that the loss of cMyBP-C results in an untethering of myosin heads. This results in an increased myosin ATP turnover, further consolidating the relationship between thick filament structure and the myosin ATPase. Show less

Several studies have investigated the effect of known adult body mass index (BMI) associated single nucleotide polymorphisms (SNPs) on BMI in childhood. There has been no genome-wide association study (GWAS) of BMI trajectories over childhood. We conducted a GWAS meta-analysis of BMI trajectories from 1 to 17 years of age in 9377 children (77,967 measurements) from the Avon Longitudinal Study of Parents and Children (ALSPAC) and the Western Australian Pregnancy Cohort (Raine) Study. Genome-wide significant loci were examined in a further 3918 individuals (48,530 measurements) from Northern Finland. Linear mixed effects models with smoothing splines were used in each cohort for longitudinal modelling of BMI. A novel SNP, downstream from the FAM120AOS gene on chromosome 9, was detected in the meta-analysis of ALSPAC and Raine. This association was driven by a difference in BMI at 8 years (T allele of rs944990 increased BMI; PSNP = 1.52 × 10(-8)), with a modest association with change in BMI over time (PWald(Change) = 0.006). Three known adult BMI-associated loci (FTO, MC4R and ADCY3) and one childhood obesity locus (OLFM4) reached genome-wide significance (PWald < 1.13 × 10(-8)) with BMI at 8 years and/or change over time. This GWAS of BMI trajectories over childhood identified a novel locus that warrants further investigation. We also observed genome-wide significance with previously established obesity loci, making the novel observation that these loci affected both the level and the rate of change in BMI. We have demonstrated that the use of repeated measures data can increase power to allow detection of genetic loci with smaller sample sizes. Show less

General cognitive function is substantially heritable across the human life course from adolescence to old age. We investigated the genetic contribution to variation in this important, health- and well-being-related trait in middle-aged and older adults. We conducted a meta-analysis of genome-wide association studies of 31 cohorts (N=53,949) in which the participants had undertaken multiple, diverse cognitive tests. A general cognitive function phenotype was tested for, and created in each cohort by principal component analysis. We report 13 genome-wide significant single-nucleotide polymorphism (SNP) associations in three genomic regions, 6q16.1, 14q12 and 19q13.32 (best SNP and closest gene, respectively: rs10457441, P=3.93 × 10(-9), MIR2113; rs17522122, P=2.55 × 10(-8), AKAP6; rs10119, P=5.67 × 10(-9), APOE/TOMM40). We report one gene-based significant association with the HMGN1 gene located on chromosome 21 (P=1 × 10(-6)). These genes have previously been associated with neuropsychiatric phenotypes. Meta-analysis results are consistent with a polygenic model of inheritance. To estimate SNP-based heritability, the genome-wide complex trait analysis procedure was applied to two large cohorts, the Atherosclerosis Risk in Communities Study (N=6617) and the Health and Retirement Study (N=5976). The proportion of phenotypic variation accounted for by all genotyped common SNPs was 29% (s.e.=5%) and 28% (s.e.=7%), respectively. Using polygenic prediction analysis, ~1.2% of the variance in general cognitive function was predicted in the Generation Scotland cohort (N=5487; P=1.5 × 10(-17)). In hypothesis-driven tests, there was significant association between general cognitive function and four genes previously associated with Alzheimer's disease: TOMM40, APOE, ABCG1 and MEF2C. Show less

Tissue concentrations of omega-3 fatty acids may reduce cardiovascular disease risk, and genetic variants are associated with circulating fatty acids concentrations. Whether dietary fatty acids interact with genetic variants to modify circulating omega-3 fatty acids is unclear. We evaluated interactions between genetic variants and fatty acid intakes for circulating alpha-linoleic acid, eicosapentaenoic acid, docosahexaenoic acid, and docosapentaenoic acid. We conducted meta-analyses (N = 11 668) evaluating interactions between dietary fatty acids and genetic variants (rs174538 and rs174548 in FADS1 (fatty acid desaturase 1), rs7435 in AGPAT3 (1-acyl-sn-glycerol-3-phosphate), rs4985167 in PDXDC1 (pyridoxal-dependent decarboxylase domain-containing 1), rs780094 in GCKR (glucokinase regulatory protein), and rs3734398 in ELOVL2 (fatty acid elongase 2)). Stratification by measurement compartment (plasma versus erthyrocyte) revealed compartment-specific interactions between FADS1 rs174538 and rs174548 and dietary alpha-linolenic acid and linoleic acid for docosahexaenoic acid and docosapentaenoic acid. Our findings reinforce earlier reports that genetically based differences in circulating fatty acids may be partially due to differences in the conversion of fatty acid precursors. Further, fatty acids measurement compartment may modify gene-diet relationships, and considering compartment may improve the detection of gene-fatty acids interactions for circulating fatty acid outcomes. Show less

Genome-wide association studies (GWAS) of BMI are mostly undertaken under the assumption that "kg/m(2) " is an index of weight fully adjusted for height, but in general this is not true. The aim here was to assess the contribution of common genetic variation to a adjusted version of that phenotype which appropriately accounts for covariation in height in children. A GWAS of height-adjusted BMI (BMI[x] = weight/height(x) ), calculated to be uncorrelated with height, in 5809 participants (mean age 9.9 years) from the Avon Longitudinal Study of Parents and Children (ALSPAC) was performed. GWAS based on BMI[x] yielded marked differences in genomewide results profile. SNPs in ADCY3 (adenylate cyclase 3) were associated at genome-wide significance level (rs11676272 (0.28 kg/m(3.1) change per allele G (0.19, 0.38), P = 6 × 10(-9) ). In contrast, they showed marginal evidence of association with conventional BMI [rs11676272 (0.25 kg/m(2) (0.15, 0.35), P = 6 × 10(-7) )]. Results were replicated in an independent sample, the Generation R study. Analysis of BMI[x] showed differences to that of conventional BMI. The association signal at ADCY3 appeared to be driven by a missense variant and it was strongly correlated with expression of this gene. Our work highlights the importance of well understood phenotype use (and the danger of convention) in characterising genetic contributions to complex traits. Show less

The analysis of rich catalogues of genetic variation from population-based sequencing provides an opportunity to screen for functional effects. Here we report a rare variant in APOC3 (rs138326449-A, minor allele frequency ~0.25% (UK)) associated with plasma triglyceride (TG) levels (-1.43 s.d. (s.e.=0.27 per minor allele (P-value=8.0 × 10(-8))) discovered in 3,202 individuals with low read-depth, whole-genome sequence. We replicate this in 12,831 participants from five additional samples of Northern and Southern European origin (-1.0 s.d. (s.e.=0.173), P-value=7.32 × 10(-9)). This is consistent with an effect between 0.5 and 1.5 mmol l(-1) dependent on population. We show that a single predicted splice donor variant is responsible for association signals and is independent of known common variants. Analyses suggest an independent relationship between rs138326449 and high-density lipoprotein (HDL) levels. This represents one of the first examples of a rare, large effect variant identified from whole-genome sequencing at a population scale. Show less

Common fragile sites (CFSs) are large regions with profound genomic instability that often span extremely large genes a number of which have been found to be important tumor suppressors. RNA sequencing previously revealed that there was a group of six large CFS genes which frequently had decreased expression in oropharyngeal squamous cell carcinomas (OPSCCs) and real-time reverse transcriptase polymerase chain reaction experiments validated that these six large CFS genes (PARK2, DLG2, NBEA, CTNNA3, DMD, and FHIT) had decreased expression in most of the tumor samples. In this study, we investigated whether the decreased expression of these genes has any clinical significance in OPSCCs. We analyzed the six CFS large genes in 45 OPSCC patients and found that 27 (60%) of the OPSCC tumors had decreased expression of these six genes. When we correlated the expression of these six genes to each patient's clinical records, for 11 patients who had tumor recurrence, 10 of them had decreased expression of almost all 6 genes. When we divided the patients into two groups, one group with decreased expression of the six genes and the other group with either slight changes or increased expression of the six genes, we found that there is significant difference in the incidence of tumor recurrence between these two groups by Kaplan-Meier plot analysis (P < .05). Our results demonstrated that those OPSCC tumors with decreased expression of this select group of six large CFS genes were much more likely to be associated with tumor recurrence and these genes are potential prognostic markers for predicting tumor recurrence in OPSCC. Show less

The common fragile sites (CFSs) are large regions of profound genomic instability found in all individuals. The frequent deletions and other alterations in these regions in multiple cancers has led to the discovery of a number of extremely large genes contained within these regions and several of the large CFS genes have already been demonstrated to function as tumor suppressors involved in the formation of many different cancers. To study the large CFS genes in oropharyngeal squamous cell carcinoma (OPSCC), we did RNA seq analysis from 11 head and neck cancer patients. This revealed that there are six large CFS genes which consistently had decreased expression in the tumor samples compared to their matched normal tissues. These six genes are PARK2, DLG2, NBEA, CTNNA3, DMD, and FHIT. PARK2 and FHIT are located within the two most frequently expressed CFSs and both have been demonstrated to function as tumor suppressors, while the other large genes are found to have frequent alterations in multiple cancers. Validation experiments using real time PCR indicated that over 60% of OPSCC tumors showed decreased expression for all six genes. Both HPV-positive and HPV-negative OPSCCs had similar proportions with loss of expression of these genes. Our results suggest that this selected group of large genes might serve as potential tumor suppressors involved in the development of OPSCCs. Further studies are needed to investigate whether the decreased expression observed is due to genomic instability within the CFS regions or the selection for alterations of specific large CFS genes. Show less

Forced vital capacity (FVC), a spirometric measure of pulmonary function, reflects lung volume and is used to diagnose and monitor lung diseases. We performed genome-wide association study meta-analysis of FVC in 52,253 individuals from 26 studies and followed up the top associations in 32,917 additional individuals of European ancestry. We found six new regions associated at genome-wide significance (P < 5 × 10(-8)) with FVC in or near EFEMP1, BMP6, MIR129-2-HSD17B12, PRDM11, WWOX and KCNJ2. Two loci previously associated with spirometric measures (GSTCD and PTCH1) were related to FVC. Newly implicated regions were followed up in samples from African-American, Korean, Chinese and Hispanic individuals. We detected transcripts for all six newly implicated genes in human lung tissue. The new loci may inform mechanisms involved in lung development and the pathogenesis of restrictive lung disease. Show less

Our understanding of how histone demethylation contributes to the regulation of basal gene expression in the brain is largely unknown in any injury model, and especially in the healthy adult brain. Although Jumonji genes are often regulated transcriptionally, cell-specific gene expression of Jumonji histone demethylases in the brain remains poorly understood. Thus, in the present study we profiled the mRNA levels of 26 Jumonji genes in microglia (CD11b+), neurons (NeuN+) and astrocytes (GFAP+) from the healthy adult rat brain. We optimized a method combining a mZBF (modified zinc-based fixative) and FCM (flow cytometry) to simultaneously sort cells from non-transgenic animals. We evaluated cell-surface, intracellular and nuclear proteins, including histones, as well as messenger- and micro-RNAs in different cell types simultaneously from a single-sorted sample. We found that 12 Jumonji genes were differentially expressed between adult microglia, neurons and astrocytes. While JMJD2D was neuron-restricted, PHF8 and JMJD1C were expressed in all three cell types although the expression was highest in neurons. JMJD3 and JMJD5 were expressed in all cell types, but were highly enriched in microglia; astrocytes had the lowest expression of UTX and JHDM1D. Levels of global H3K27 (H3 lysine 27) methylation varied among cell types and appeared to be lowest in microglia, indicating that differences in basal gene expression of specific Jumonji histone demethylases may contribute to cell-specific gene expression in the CNS (central nervous system). This multiparametric technique will be valuable for simultaneously assaying chromatin modifications and gene regulation in the adult CNS. Show less

Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation. Show less

Compared to DBA/2J (D2), C57BL/6J (B6) inbred mice exhibit strong morphine preference when tested using a two-bottle choice drinking paradigm. A morphine preference quantitative trait locus (QTL), Mop2, was originally mapped to proximal chromosome (Chr) 10 using a B6xD2 F2 intercross population, confirmed with reciprocal congenic strains and fine mapped with recombinant congenic strains. These efforts identified a ∼ 10-Million base pair (Mbp) interval, underlying Mop2, containing 35 genes. To further reduce the interval, mice from the D2.B6-Mop2-P1 congenic strain were backcrossed to parental D2 mice and two new recombinant strains of interest were generated: D2.B6-Mop2-P1.pD.dB and D2.B6-Mop2-P1.pD.dD. Results obtained from testing these strains in the two-bottle choice drinking paradigm suggest that the gene(s) responsible for the Mop2 QTL is one or more of 22 remaining within the newly defined interval (∼ 7.6 Mbp) which includes Oprm1 and several other genes related to opioid pharmacology. Real-time qRT-PCR analysis of Oprm1 and opioid-related genes Rgs17, Ppp1r14c, Vip, and Iyd revealed both between-strain and within-strain expression differences in comparisons of saline- and morphine-treated B6 and D2 mice. Analysis of Rgs17 protein levels also revealed both between-strain and within-strain differences in comparisons of saline- and morphine-treated B6 and D2 mice. Results suggest that the Mop2 QTL represents the combined influence of multiple genetic variants on morphine preference in these two strains. Relative contributions of each variant remain to be determined. Show less

The pubertal height growth spurt is a distinctive feature of childhood growth reflecting both the central onset of puberty and local growth factors. Although little is known about the underlying genetics, growth variability during puberty correlates with adult risks for hormone-dependent cancer and adverse cardiometabolic health. The only gene so far associated with pubertal height growth, LIN28B, pleiotropically influences childhood growth, puberty and cancer progression, pointing to shared underlying mechanisms. To discover genetic loci influencing pubertal height and growth and to place them in context of overall growth and maturation, we performed genome-wide association meta-analyses in 18 737 European samples utilizing longitudinally collected height measurements. We found significant associations (P < 1.67 × 10(-8)) at 10 loci, including LIN28B. Five loci associated with pubertal timing, all impacting multiple aspects of growth. In particular, a novel variant correlated with expression of MAPK3, and associated both with increased prepubertal growth and earlier menarche. Another variant near ADCY3-POMC associated with increased body mass index, reduced pubertal growth and earlier puberty. Whereas epidemiological correlations suggest that early puberty marks a pathway from rapid prepubertal growth to reduced final height and adult obesity, our study shows that individual loci associating with pubertal growth have variable longitudinal growth patterns that may differ from epidemiological observations. Overall, this study uncovers part of the complex genetic architecture linking pubertal height growth, the timing of puberty and childhood obesity and provides new information to pinpoint processes linking these traits. Show less

Apolipoproteins have an important role in lipid metabolism and transport. Polymorphisms in the APOA1/C3/A4/A5 gene cluster have been associated with lipid alterations and cardiovascular diseases. We investigated APOA1 XmnI, APOA5 S19W, and APOA5 -1131T>C polymorphisms in 377 individuals from a cohort of a longitudinal Brazilian elderly study. Allele frequencies, genotype distribution, and association with major morbidities as well as with lipids, creatinine, albumin, urea, glycated hemoglobin, and fasting glucose serum levels were investigated. Linkage disequilibrium and haplotype associations were also analyzed. This is the first time that haplotypes involving these polymorphisms were evaluated. Genotyping was performed by PCR-RFLP. Minor allele frequencies were 0.119, 0.071, and 0.158 for XmnI, S19W, and -1131T>C polymorphisms, respectively. We found a significant association of the -1131C allele with low LDL-C levels. We also observed that XmnI and S19W polymorphisms were in linkage disequilibrium. The C/G haplotype, which is composed of the wild-type allele of XmnI and the minor allele of S19W, was associated with high total cholesterol serum levels in this elderly population. We conclude that the -1131T>C polymorphism and the C/G haplotype, including XmnI and S19W polymorphisms, are associated with alterations in lipid levels and may be risk factors for cardiovascular disease in the Brazilian elderly. Show less

Fetal umbilical cord HDL concentration is lower in IUGR neonates as compared to gestational age matched controls (CTRL). The causes by now are unknown. A full apolipoprotein analysis of cord blood might help in understanding the changes in lipid metabolism seen in IUGR. To characterize cord blood apolipoprotein profile of IUGR neonates. Serum of venous umbilical cord blood (15 IUGR vs. 15 CTRL) was analyzed by Multiple Reaction Monitoring (MRM). 15 different known apolipoproteins were profiled. HDL and LDL were measured by colorimetric methods in fetal cord blood and their corresponding mothers. Fetal HDL (p<0.0001), ApoC1 (p<0.0001), and ApoE (p=0.0001) levels were lower in IUGR as compared to CTRL. Fetal HDL levels were positive correlated to ApoE, ApoC1, and ApoA2 (r=0.79, r=0.74, r=0.56). Fetal LDL levels were positive correlated to ApoB, ApoE, ApoA2, and ApoC3 (r=0.74, r=0.67, r=0.57, r=0.55). Maternal LDL concentrations correlated positive to fetal ApoC1, ApoC2, and LCAT-concentration (r=0.54, r=0.52, r=0.52). The results underlines the relevance of ApoE in fetal development. Moreover, we speculate that maternal lipid profile has an impact on fetal lipid metabolisms as evidenced by the association of maternal LDL levels and fetal ApoC1, ApoC2, and LCAT concentrations. This observation requires further confirmation and is worth to be analyzed since it provides a mechanistic link for therapeutic options. Show less

Circulating metabolites associated with insulin sensitivity may represent useful biomarkers, but their causal role in insulin sensitivity and diabetes is less certain. We previously identified novel metabolites correlated with insulin sensitivity measured by the hyperinsulinemic-euglycemic clamp. The top-ranking metabolites were in the glutathione and glycine biosynthesis pathways. We aimed to identify common genetic variants associated with metabolites in these pathways and test their role in insulin sensitivity and type 2 diabetes. With 1,004 nondiabetic individuals from the RISC study, we performed a genome-wide association study (GWAS) of 14 insulin sensitivity-related metabolites and one metabolite ratio. We replicated our results in the Botnia study (n = 342). We assessed the association of these variants with diabetes-related traits in GWAS meta-analyses (GENESIS [including RISC, EUGENE2, and Stanford], MAGIC, and DIAGRAM). We identified four associations with three metabolites-glycine (rs715 at CPS1), serine (rs478093 at PHGDH), and betaine (rs499368 at SLC6A12; rs17823642 at BHMT)-and one association signal with glycine-to-serine ratio (rs1107366 at ALDH1L1). There was no robust evidence for association between these variants and insulin resistance or diabetes. Genetic variants associated with genes in the glycine biosynthesis pathways do not provide consistent evidence for a role of glycine in diabetes-related traits. Show less

The WWP2 E3 ubiquitin ligase has previously been shown to regulate TGFβ/Smad signalling activity linked to epithelial-mesenchymal transition (EMT). Whilst inhibitory I-Smad7 was found to be the preferred substrate for full-length WWP2-FL and a WWP2-C isoform, WWP2-FL also formed a stable complex with an N-terminal WWP2 isoform (WWP2-N) in the absence of TGFβ, and rapidly stimulated activating Smad2/3 turnover. Here, using stable knockdown experiments we show that specific depletion of individual WWP2 isoforms impacts differentially on Smad protein levels, and in WWP2-N knockdown cells we unexpectedly find spontaneous expression of the EMT marker vimentin. Re-introduction of WWP2-N into WWP2-N knockout cells also repressed TGFβ-induced vimentin expression. In support of the unique role for WWP2-N in regulating TGFβ/Smad functional activity, we then show that a novel V717M-WWP2 mutant in the MZ7-mel melanoma cell line forms a stable complex with the WWP2-N isoform and promotes EMT by stabilizing Smad3 protein levels. Finally, we report the first analysis of WWP2 expression in cancer cDNA panel arrays using WWP2 isoform-specific probes and identify unique patterns of WWP2 isoform abundance associated with early/advanced disease stages. WWP2-N is significantly downregulated in stage IIIC melanoma and up-regulated in stage II/III prostate cancer, and we also find isolated examples of WWP2-FL and WWP2-C overexpression in early-stage breast cancer. Together, these data suggest that individual WWP2 isoforms, and particularly WWP2-N, could play central roles in tumourigenesis linked to aberrant TGFβ-dependent signalling function, and also have potential as both prognostic markers and molecular therapeutic targets. Show less

Genome wide association studies (GWAS), mostly in Europeans have identified several common variants as associated with key lipid traits. Replication of these genetic effects in South Asian populations is important since it would suggest wider relevance for these findings. Given the rising prevalence of metabolic disorders and heart disease in the Indian sub-continent, these studies could be of future clinical relevance. We studied seven common variants associated with a variety of lipid traits in previous GWASs. The study sample comprised of 3178 sib-pairs recruited as participants for the Indian Migration Study (IMS). Associations with various lipid parameters and quantitative traits were analyzed using the Fulker genetic association model. We replicated five of the 7 main effect associations with p-values ranging from 0.03 to 1.97x10(-7). We identified particularly strong association signals at rs662799 in APOA5 (beta=0.18 s.d, p=1.97 x 10(-7)), rs10503669 in LPL (beta =-0.18 s.d, p=1.0 x 10(-4)) and rs780094 in GCKR (beta=0.11 s.d, p=0.001) loci in relation to triglycerides. In addition, the GCKR variant was also associated with total cholesterol (beta=0.11 s.d, p=3.9x10(-4)). We also replicated the association of rs562338 in APOB (p=0.03) and rs4775041 in LIPC (p=0.007) with LDL-cholesterol and HDL-cholesterol respectively. We report associations of five loci with various lipid traits with the effect size consistent with the same reported in Europeans. These results indicate an overlap of genetic effects pertaining to lipid traits across the European and Indian populations. Show less

γ-Secretase inhibitors (GSIs) are being investigated for their potential to modify the progression of Alzheimer disease based on their ability to regulate amyloid-β (Aβ) accumulation. BMS-708163 (avagacestat) is an oral GSI designed for selective inhibition of Aβ synthesis currently in development for the treatment of mild to moderate and predementia AD. In addition to the desired effect on Aβ synthesis, GSIs affect Notch processing, which is thought to mediate some toxic adverse effects reported with this drug class. Avagacestat produced up to 190-fold greater selectivity for Aβ synthesis than Notch processing in preclinical studies and may therefore produce less toxic adverse events than other less selective compounds. Presented here are the results of the first in-human study for this new GSI compound. The goal of this study was to assess the tolerability profile, pharmacokinetic properties, and effects on pharmacodynamic markers (Aβ, trefoil factor family 3 protein, dual specificity phosphatase 6, and hairy and enhancer of split-1) of single, oral doses of avagacestat in healthy, young, male volunteers. This was a multicenter, randomized, double-blind, placebo-controlled, single-ascending dose study in 8 healthy young men (age, 18-45 years) per dosing panel. Each study participant was randomized to receive a single dose of placebo (n = 2) or avagacestat (n = 6 for each dose) as an oral solution in 1 of 9 sequential dose panels (0.3, 1.5, 5, 15, 50, 100, 200, 400, and 800 mg). For determination of avagacestat, blood samples were obtained before dosing and for up to 144 hours after dosing. For participants in the 800-mg avagacestat dose panel, additional samples were obtained at 216, 312, and 648 hours. For 40-amino acid isoform of Aβ (Aβ(1-40)) assessment, plasma samples were collected before avagacestat administration and up to 72 hours after dosing. Avagacestat concentrations peaked quickly after oral administration and then had a biphasic decrease in concentrations with a prolonged terminal phase. Exposures were proportional with doses up to 200 mg. Avagacestat was well tolerated at single oral doses up to 800 mg, with a biphasic effect on plasma Aβ(1-40). Adverse events were predominately mild to moderate in severity with no evidence of dose dependence up to 200 mg. Results from this single-ascending dose study suggest that avagacestat was tolerated at a single-dose range of 0.3 to 800 mg and suitable for further clinical development. Show less

The mitogen-activated protein kinase (MAPK) and phosphoinositide 3-kinase/AKT signaling pathways interact at multiple nodes in cancer, including at mTOR complexes, suggesting an increased likelihood of redundancy and innate resistance to any therapeutic effects of single pathway inhibition. In this study, we investigated the therapeutic effects of combining the MAPK extracellular signal-regulated kinase (MEK)1/2 inhibitor selumetinib (AZD6244) with the dual mTORC1 and mTORC2 inhibitor (AZD8055). Concurrent dosing in nude mouse xenograft models of human lung adenocarcinoma (non-small cell lung cancers) and colorectal carcinoma was well tolerated and produced increased antitumor efficacy relative to the respective monotherapies. Pharmacodynamic analysis documented reciprocal pathway inhibition associated with increased apoptosis and Bim expression in tumor tissue from the combination group, where key genes such as DUSP6 that are under MEK functional control were also modulated. Our work offers a strong rationale to combine selumetinib and AZD8055 in clinical trials as an attractive therapeutic strategy. Show less

Circulating androgen levels are often used as indicators of physiological or pathological conditions. More than half of the variance for circulating androgen levels is thought to be genetically influenced. A genome-wide association study (GWAS) has identified two loci, SHBG at 17p13 and FAM9B at Xp22, for serum testosterone (T) levels; however, these explain only a small fraction of inter-individual variability. To identify additional genetic determinants of androgen levels, a GWAS of baseline serum T and dihydrotestosterone (DHT) levels was conducted in 3225 men of European ancestry from the REduction by DUtasteride of Prostate Cancer Events (REDUCE) study. Cross-validation was used to confirm the observed associations between the drug (n = 1581) and placebo (n = 1644) groups of REDUCE. In addition to confirming the associations of two known loci with serum T levels (rs727428 in SHBG: P = 1.26 × 10(-12); rs5934505 in FAM9B: P = 1.61 × 10(-8)), we identified a new locus, JMJD1C at 10q21 that was associated with serum T levels at a genome-wide significance level (rs10822184: P = 1.12 × 10(-8)). We also observed that the SHBG locus was associated with serum DHT levels (rs727428: P = 1.47 × 10(-11)). Moreover, two additional variants in SHBG [rs72829446, in strong linkage equilibrium with the missense variant D356N (rs6259), and rs1799941] were also independently associated with circulating androgen levels in a statistical scale. These three loci (JMJD1C, SHBG and FAM9B) were estimated to account for ~5.3 and 4.1% of the variance of serum T and DHT levels. Our findings may provide new insights into the regulation of circulating androgens and potential targets for androgen-based therapy. Show less