Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

Genetic risk converges on regulatory networks mediating early type 2 diabetes

Type 2 diabetes mellitus (T2D), a major cause of worldwide morbidity and mortality, is characterized by dysfunction of insulin-producing pancreatic islet β cells1,2. T2D genome-wide association studies (GWAS) have identified hundreds of signals in non-coding and β cell regulatory genomic regions, but deciphering their biological mechanisms remains challenging3-5. Here, to identify early disease-driving events, we performed traditional and multiplexed pancreatic tissue imaging, sorted-islet cell transcriptomics and islet functional analysis of early-stage T2D and control donors. By integrating diverse modalities, we show that early-stage T2D is characterized by β cell-intrinsic defects that can be proportioned into gene regulatory modules with enrichment in signals of genetic risk. After identifying the β cell hub gene and transcription factor RFX6 within one such module, we demonstrated multiple layers of genetic risk that converge on an RFX6-mediated network to reduce insulin secretion by β cells. RFX6 perturbation in primary human islet cells alters β cell chromatin architecture at regions enriched for T2D GWAS signals, and population-scale genetic analyses causally link genetically predicted reduced RFX6 expression with increased T2D risk. Understanding the molecular mechanisms of complex, systemic diseases necessitates integration of signals from multiple molecules, cells, organs and individuals, and thus we anticipate that this approach will be a useful template to identify and validate key regulatory networks and master hub genes for other diseases or traits using GWAS data.

An integrative single-cell multi-omics profiling of human pancreatic islets identifies T1D associated genes and regulatory signals

Genome wide association studies (GWAS) have identified over 100 signals associated with type 1 diabetes (T1D). However, translating any given T1D GWAS signal into mechanistic insights, including putative causal variants and the context (cell type and cell state) in which they function, has been limited. Here, we present a comprehensive multi-omic integrative analysis of single-cell/nucleus resolution profiles of gene expression and chromatin accessibility in healthy and autoantibody+ (AAB+) human islets, as well as islets under multiple T1D stimulatory conditions. We broadly nominate effector cell types for all T1D GWAS signals. We further nominated higher-resolution contexts, including effector cell types, regulatory elements, and genes for three independent T1D risk variants acting through islet cells within the pancreas at the DLK1/MEG3, RASGRP1, and TOX loci. Subsequently, we created isogenic gene knockouts DLK1-/-, RASGRP1-/-, and TOX-/-, and the corresponding regulatory region knockout, RASGRP1Δ, and DLK1Δ hESCs. Loss of RASGRP1 or DLK1, as well as knockout of the regulatory region of RASGRP1 or DLK1, increased β cell apoptosis. Additionally, pancreatic β cells derived from isogenic hESCs carrying the risk allele of rs3783355A/A exhibited increased β cell death. Finally, RNA-seq and ATAC-seq identified five genes upregulated in both RASGRP1-/- and DLK1-/- β-like cells, four of which are associated with T1D. Together, this work reports an integrative approach for combining single cell multi-omics, GWAS, and isogenic hESC-derived β-like cells to prioritize the T1D associated signals and their underlying context-specific cell types, genes, SNPs, and regulatory elements, to illuminate biological functions and molecular mechanisms.

FixItFelix: improving genomic analysis by fixing reference errors

The current version of the human reference genome, GRCh38, contains a number of errors including 1.2 Mbp of falsely duplicated and 8.04 Mbp of collapsed regions. These errors impact the variant calling of 33 protein-coding genes, including 12 with medical relevance. Here, we present FixItFelix, an efficient remapping approach, together with a modified version of the GRCh38 reference genome that improves the subsequent analysis across these genes within minutes for an existing alignment file while maintaining the same coordinates. We showcase these improvements over multi-ethnic control samples, demonstrating improvements for population variant calling as well as eQTL studies.

Dilated Optic Nerve Sheath in Mucopolysaccharidosis I: Common and Not Necessarily High Intracranial Pressure

Hydrocephalus is one of the earliest manifestations of mucopolysaccharidosis I-Hurler syndrome, and delayed treatment of hydrocephalus can lead to neurocognitive delay or even death. Optic nerve sheath diameter has been established as a noninvasive measurement to detect elevated intracranial pressure. This study aimed to establish correlations between optic nerve sheath diameter and opening pressure. Forty-nine MR images and opening pressures in patients with mucopolysaccharidosis I-Hurler syndrome were retrospectively reviewed from 2008 to 2020. The optic nerve sheath diameter was measured 3 mm posterior to the posterior margin of the globe (retrobulbar) and 10 mm anterior to the optic foramen (midpoint segment), and the average was taken between the 2 eyes. Opening pressure was measured with the patient in the lateral decubitus position with controlled end-tidal CO2 on the same day as the MR imaging. The average retrobulbar optic nerve sheath diameter was 5.33 mm, higher than the previously reported measurement in healthy controls, in patients with idiopathic intracranial hypertension, and there was a positive correlation between age and the optic nerve sheath diameter measured at the retrobulbar or midpoint segment (retrobulbar segment, R 2 = 0.27, P < .01; midpoint segment, R 2 = 0.20, P < .01). However, there was no correlation between retrobulbar or midpoint segment optic nerve sheath diameter and opening pressure (retrobulbar segment, R 2 = 0.02, P = .17; midpoint segment, R 2 = 0.03, P < .12). This study shows a higher average optic nerve sheath diameter in patients with mucopolysaccharidosis I-Hurler syndrome than in healthy controls regardless of the location of the measurement. However, the degree of optic nerve sheath dilation does not correlate with opening pressure, suggesting that increased optic nerve sheath diameter is an ocular manifestation of mucopolysaccharidosis I-Hurler syndrome itself rather than a marker of elevated intracranial pressure.

Comprehensive enhancer-target gene assignments improve gene set level interpretation of genome-wide regulatory data

BACKGROUND

Revealing the gene targets of distal regulatory elements is challenging yet critical for interpreting regulome data. Experiment-derived enhancer-gene links are restricted to a small set of enhancers and/or cell types, while the accuracy of genome-wide approaches remains elusive due to the lack of a systematic evaluation. We combined multiple spatial and in silico approaches for defining enhancer locations and linking them to their target genes aggregated across >500 cell types, generating 1860 human genome-wide distal enhancer-to-target gene definitions (EnTDefs). To evaluate performance, we used gene set enrichment (GSE) testing on 87 independent ENCODE ChIP-seq datasets of 34 transcription factors (TFs) and assessed concordance of results with known TF Gene Ontology annotations, and other benchmarks.

RESULTS

The top ranked 741 (40%) EnTDefs significantly outperform the common, naïve approach of linking distal regions to the nearest genes, and the top 10 EnTDefs perform well when applied to ChIP-seq data of other cell types. The GSE-based ranking of EnTDefs is highly concordant with ranking based on overlap with curated benchmarks of enhancer-gene interactions. Both our top general EnTDef and cell-type-specific EnTDefs significantly outperform seven independent computational and experiment-based enhancer-gene pair datasets. We show that using our top EnTDefs for GSE with either genome-wide DNA methylation or ATAC-seq data is able to better recapitulate the biological processes changed in gene expression data performed in parallel for the same experiment than our lower-ranked EnTDefs.

CONCLUSIONS

Our findings illustrate the power of our approach to provide genome-wide interpretation regardless of cell type.

Human and rat skeletal muscle single-nuclei multi-omic integrative analyses nominate causal cell types, regulatory elements, and SNPs for complex traits

Skeletal muscle accounts for the largest proportion of human body mass, on average, and is a key tissue in complex diseases and mobility. It is composed of several different cell and muscle fiber types. Here, we optimize single-nucleus ATAC-seq (snATAC-seq) to map skeletal muscle cell-specific chromatin accessibility landscapes in frozen human and rat samples, and single-nucleus RNA-seq (snRNA-seq) to map cell-specific transcriptomes in human. We additionally perform multi-omics profiling (gene expression and chromatin accessibility) on human and rat muscle samples. We capture type I and type II muscle fiber signatures, which are generally missed by existing single-cell RNA-seq methods. We perform cross-modality and cross-species integrative analyses on 33,862 nuclei and identify seven cell types ranging in abundance from 59.6% to 1.0% of all nuclei. We introduce a regression-based approach to infer cell types by comparing transcription start site-distal ATAC-seq peaks to reference enhancer maps and show consistency with RNA-based marker gene cell type assignments. We find heterogeneity in enrichment of genetic variants linked to complex phenotypes from the UK Biobank and diabetes genome-wide association studies in cell-specific ATAC-seq peaks, with the most striking enrichment patterns in muscle mesenchymal stem cells (∼3.5% of nuclei). Finally, we overlay these chromatin accessibility maps on GWAS data to nominate causal cell types, SNPs, transcription factor motifs, and target genes for type 2 diabetes signals. These chromatin accessibility profiles for human and rat skeletal muscle cell types are a useful resource for nominating causal GWAS SNPs and cell types.

A Transcription Start Site Map in Human Pancreatic Islets Reveals Functional Regulatory Signatures

Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), ∼20% of which are islet specific and occur mostly distal to known gene TSS. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validated the transcriptional enhancer activity for 2,279 of 3,378 (∼68%) tested islet CAGE elements (5% false discovery rate). TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting active enhancers at GWAS loci.

Genetic effects on liver chromatin accessibility identify disease regulatory variants

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

Chromatin information content landscapes inform transcription factor and DNA interactions

Interactions between transcription factors and chromatin are fundamental to genome organization and regulation and, ultimately, cell state. Here, we use information theory to measure signatures of organized chromatin resulting from transcription factor-chromatin interactions encoded in the patterns of the accessible genome, which we term chromatin information enrichment (CIE). We calculate CIE for hundreds of transcription factor motifs across human samples and identify two classes: low and high CIE. The 10-20% of common and tissue-specific high CIE transcription factor motifs, associate with higher protein-DNA residence time, including different binding site subclasses of the same transcription factor, increased nucleosome phasing, specific protein domains, and the genetic control of both chromatin accessibility and gene expression. These results show that variations in the information encoded in chromatin architecture reflect functional biological variation, with implications for cell state dynamics and memory.

Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

BAF60a deficiency uncouples chromatin accessibility and cold sensitivity from white fat browning

Brown and beige fat share a remarkably similar transcriptional program that supports fuel oxidation and thermogenesis. The chromatin-remodeling machinery that governs genome accessibility and renders adipocytes poised for thermogenic activation remains elusive. Here we show that BAF60a, a subunit of the SWI/SNF chromatin-remodeling complexes, serves an indispensable role in cold-induced thermogenesis in brown fat. BAF60a maintains chromatin accessibility at PPARγ and EBF2 binding sites for key thermogenic genes. Surprisingly, fat-specific BAF60a inactivation triggers more pronounced cold-induced browning of inguinal white adipose tissue that is linked to induction of MC2R, a receptor for the pituitary hormone ACTH. Elevated MC2R expression sensitizes adipocytes and BAF60a-deficient adipose tissue to thermogenic activation in response to ACTH stimulation. These observations reveal an unexpected dichotomous role of BAF60a-mediated chromatin remodeling in transcriptional control of brown and beige gene programs and illustrate a pituitary-adipose signaling axis in the control of thermogenesis.

Multiomic Profiling Identifies cis-Regulatory Networks Underlying Human Pancreatic β Cell Identity and Function

EndoC-βH1 is emerging as a critical human β cell model to study the genetic and environmental etiologies of β cell (dys)function and diabetes. Comprehensive knowledge of its molecular landscape is lacking, yet required, for effective use of this model. Here, we report chromosomal (spectral karyotyping), genetic (genotyping), epigenomic (ChIP-seq and ATAC-seq), chromatin interaction (Hi-C and Pol2 ChIA-PET), and transcriptomic (RNA-seq and miRNA-seq) maps of EndoC-βH1. Analyses of these maps define known (e.g., PDX1 and ISL1) and putative (e.g., PCSK1 and mir-375) β cell-specific transcriptional cis-regulatory networks and identify allelic effects on cis-regulatory element use. Importantly, comparison with maps generated in primary human islets and/or β cells indicates preservation of chromatin looping but also highlights chromosomal aberrations and fetal genomic signatures in EndoC-βH1. Together, these maps, and a web application we created for their exploration, provide important tools for the design of experiments to probe and manipulate the genetic programs governing β cell identity and (dys)function in diabetes.

EndoC-βH1 is becoming an important cellular model to study genes and pathways governing human β cell identity and function, but its (epi)genomic similarity to primary human islets is unknown. Lawlor et al. complete and compare extensive EndoC and primary human islet multiomic maps to identify shared and distinct genomic circuitry.

Genome-wide chromatin accessibility and transcriptome profiling show minimal epigenome changes and coordinated transcriptional dysregulation of hedgehog signaling in Danforth's short tail mice

Danforth's short tail (Sd) mice provide an excellent model for investigating the underlying etiology of human caudal birth defects, which affect 1 in 10 000 live births. Sd animals exhibit aberrant axial skeleton, urogenital and gastrointestinal development similar to human caudal malformation syndromes including urorectal septum malformation, caudal regression, vertebral-anal-cardiac-tracheo-esophageal fistula-renal-limb (VACTERL) association and persistent cloaca. Previous studies have shown that the Sd mutation results from an endogenous retroviral (ERV) insertion upstream of the Ptf1a gene resulting in its ectopic expression at E9.5. Though the genetic lesion has been determined, the resulting epigenomic and transcriptomic changes driving the phenotype have not been investigated. Here, we performed ATAC-seq experiments on isolated E9.5 tailbud tissue, which revealed minimal changes in chromatin accessibility in Sd/Sd mutant embryos. Interestingly, chromatin changes were localized to a small interval adjacent to the Sd ERV insertion overlapping a known Ptf1a enhancer region, which is conserved in mice and humans. Furthermore, mRNA-seq experiments revealed increased transcription of Ptf1a target genes and, importantly, downregulation of hedgehog pathway genes. Reduced sonic hedgehog (SHH) signaling was confirmed by in situ hybridization and immunofluorescence suggesting that the Sd phenotype results, in part, from downregulated SHH signaling. Taken together, these data demonstrate substantial transcriptome changes in the Sd mouse, and indicate that the effect of the ERV insertion on Ptf1a expression may be mediated by increased chromatin accessibility at a conserved Ptf1a enhancer. We propose that human caudal dysgenesis disorders may result from dysregulation of hedgehog signaling pathways.

Supporting Cancer Patients Who Are Also Parents: Establishment of a Cross-Sector Service for Families With Adolescent and Young Adult Children

Background and context: Traditionally adult hospitals focus on the patient and less on their family. Adolescent and young adult children of cancer patients (AYA offspring) have significant psychosocial burdens associated with their parent's cancer however they are often invisible within hospitals with no clear referral pathways to community-based support. AYA offspring are 3-6 times more likely than peers to have clinically elevated levels of distress which increases with age. Research shows that 1 of the greatest concerns for parents is how to communicate about cancer with their children, and 1 of the greatest needs for AYA offspring is information about their parent's cancer and talking with their parents about it. Aim: CanTeen, a national AYA cancer community support organization, sought to address these needs by embedding a Parent Support Worker within the social work teams of tertiary hospitals for patients who are parents of AYA children. This service aims to assist with specific parenting challenges that arise due to a cancer diagnosis and establish a referral pathway for AYA offspring to CanTeen for support. Strategy/Tactics: Cofunding with philanthropic organizations was sought and CanTeen executives engaged in advocacy work with senior hospital management, demonstrating the need for the service and a plan to embed it within the existing hospital social work team and services. Program/Policy process: A new Parent Support Worker role was established to provide social work care to parents of AYA children following a parent's cancer diagnosis. The service provides support with parenting issues that arise due to the diagnosis as well as staff education, secondary consultations, and referrals of AYA offspring to CanTeen. A service improvement approach has been established with the collection of monitoring data measuring volume of referrals/sessions, capacity building of other staff, information given to young people and referrals to CanTeen. A service evaluation seeking feedback from patients and relevant hospital staff is also underway. Outcomes: To date, philanthropic funding and hospital support was gained to establish a pilot program placing a Parent Support Worker in 3 hospitals. Early indications are that this novel service is integrating well into established hospital teams and processes, and adding considerable value with the provision of this focused family support. What was learned: This advocacy initiative is highlighting the benefits of a model of care in adult hospitals that centrally considers family and the strength of a well-planned cross-sector service initiative.

Improving Prediction of Risk of Hospital Admission in Chronic Obstructive Pulmonary Disease: Application of Machine Learning to Telemonitoring Data

BACKGROUND

Telemonitoring of symptoms and physiological signs has been suggested as a means of early detection of chronic obstructive pulmonary disease (COPD) exacerbations, with a view to instituting timely treatment. However, algorithms to identify exacerbations result in frequent false-positive results and increased workload. Machine learning, when applied to predictive modelling, can determine patterns of risk factors useful for improving prediction quality.

OBJECTIVE

Our objectives were to (1) establish whether machine learning techniques applied to telemonitoring datasets improve prediction of hospital admissions and decisions to start corticosteroids, and (2) determine whether the addition of weather data further improves such predictions.

METHODS

We used daily symptoms, physiological measures, and medication data, with baseline demography, COPD severity, quality of life, and hospital admissions from a pilot and large randomized controlled trial of telemonitoring in COPD. We linked weather data from the United Kingdom meteorological service. We used feature selection and extraction techniques for time series to construct up to 153 predictive patterns (features) from symptom, medication, and physiological measurements. We used the resulting variables to construct predictive models fitted to training sets of patients and compared them with common symptom-counting algorithms.

RESULTS

We had a mean 363 days of telemonitoring data from 135 patients. The two most practical traditional score-counting algorithms, restricted to cases with complete data, resulted in area under the receiver operating characteristic curve (AUC) estimates of 0.60 (95% CI 0.51-0.69) and 0.58 (95% CI 0.50-0.67) for predicting admissions based on a single day's readings. However, in a real-world scenario allowing for missing data, with greater numbers of patient daily data and hospitalizations (N=57,150, N+=55, respectively), the performance of all the traditional algorithms fell, including those based on 2 days' data. One of the most frequently used algorithms performed no better than chance. All considered machine learning models demonstrated significant improvements; the best machine learning algorithm based on 57,150 episodes resulted in an aggregated AUC of 0.74 (95% CI 0.67-0.80). Adding weather data measurements did not improve the predictive performance of the best model (AUC 0.74, 95% CI 0.69-0.79). To achieve an 80% true-positive rate (sensitivity), the traditional algorithms were associated with an 80% false-positive rate: our algorithm halved this rate to approximately 40% (specificity approximately 60%). The machine learning algorithm was moderately superior to the best symptom-counting algorithm (AUC 0.77, 95% CI 0.74-0.79 vs AUC 0.66, 95% CI 0.63-0.68) at predicting the need for corticosteroids.

CONCLUSIONS

Early detection and management of COPD remains an important goal given its huge personal and economic costs. Machine learning approaches, which can be tailored to an individual's baseline profile and can learn from experience of the individual patient, are superior to existing predictive algorithms and show promise in achieving this goal.

TRIAL REGISTRATION

International Standard Randomized Controlled Trial Number ISRCTN96634935; http://www.isrctn.com/ISRCTN96634935 (Archived by WebCite at http://www.webcitation.org/722YkuhAz).

Genome-wide association study of response to methotrexate in early rheumatoid arthritis patients

Methotrexate (MTX) monotherapy is a common first treatment for rheumatoid arthritis (RA), but many patients do not respond adequately. In order to identify genetic predictors of response, we have combined data from two consortia to carry out a genome-wide study of response to MTX in 1424 early RA patients of European ancestry. Clinical endpoints were change from baseline to 6 months after starting treatment in swollen 28-joint count, tender 28-joint count, C-reactive protein and the overall 3-component disease activity score (DAS28). No single nucleotide polymorphism (SNP) reached genome-wide statistical significance for any outcome measure. The strongest evidence for association was with rs168201 in NRG3 (p = 10^-7 for change in DAS28). Some support was also seen for association with ZMIZ1, previously highlighted in a study of response to MTX in juvenile idiopathic arthritis. Follow-up in two smaller cohorts of 429 and 177 RA patients did not support these findings, although these cohorts were more heterogeneous.

Population-based cohort study of variation in the use of emergency cholecystectomy for benign gallbladder diseases

Background

The aims of this prospective population-based cohort study were to identify the patient and hospital characteristics associated with emergency cholecystectomy, and the influences of these in determining variations between hospitals.

Methods

Data were collected for consecutive patients undergoing cholecystectomy in acute UK and Irish hospitals between 1 March and 1 May 2014. Potential explanatory variables influencing the performance of emergency cholecystectomy were analysed by means of multilevel, multivariable logistic regression modelling using a two-level hierarchical structure with patients (level 1) nested within hospitals (level 2).

Results

Data were collected on 4744 cholecystectomies from 165 hospitals. Increasing age, lower ASA fitness grade, biliary colic, the need for further imaging (magnetic retrograde cholangiopancreatography), endoscopic interventions (endoscopic retrograde cholangiopancreatography) and admission to a non-biliary centre significantly reduced the likelihood of an emergency cholecystectomy being performed. The multilevel model was used to calculate the probability of receiving an emergency cholecystectomy for a woman aged 40 years or over with an ASA grade of I or II and a BMI of at least 25·0 kg/m2, who presented with acute cholecystitis with an ultrasound scan showing a thick-walled gallbladder and a normal common bile duct. The mean predicted probability of receiving an emergency cholecystectomy was 0·52 (95 per cent c.i. 0·45 to 0·57). The predicted probabilities ranged from 0·02 to 0·95 across the 165 hospitals, demonstrating significant variation between hospitals.

Conclusion

Patients with similar characteristics presenting to different hospitals with acute gallbladder pathology do not receive comparable care.

Glycosylation of plasma IgG in colorectal cancer prognosis

In this study we demonstrate the potential value of Immunoglobulin G (IgG) glycosylation as a novel prognostic biomarker of colorectal cancer (CRC). We analysed plasma IgG glycans in 1229 CRC patients and correlated with survival outcomes. We assessed the predictive value of clinical algorithms and compared this to algorithms that also included glycan predictors. Decreased galactosylation, decreased sialylation (of fucosylated IgG glycan structures) and increased bisecting GlcNAc in IgG glycan structures were strongly associated with all-cause (q < 0.01) and CRC mortality (q = 0.04 for galactosylation and sialylation). Clinical algorithms showed good prediction of all-cause and CRC mortality (Harrell's C: 0.73, 0.77; AUC: 0.75, 0.79, IDI: 0.02, 0.04 respectively). The inclusion of IgG glycan data did not lead to any statistically significant improvements overall, but it improved the prediction over clinical models for stage 4 patients with the shortest follow-up time until death, with the median gain in the test AUC of 0.08. These glycan differences are consistent with significantly increased IgG pro-inflammatory activity being associated with poorer CRC prognosis, especially in late stage CRC. In the absence of validated biomarkers to improve upon prognostic information from existing clinicopathological factors, the potential of these novel IgG glycan biomarkers merits further investigation.

A new Australian online and phone mental health support service for young people living with cancer

OBJECTIVE

The purpose of CanTeen's E-Mental Health Service for Young People Living With Cancer (YPLWC) is to meet the unique psychosocial needs of young people (12-24 years) in Australia impacted by cancer (either as a patient or family member of someone with cancer).

CONCLUSIONS

This online platform will provide the primary site where all YPLWC can find information, connect with others going through a similar experience, express their feelings, utilise tools for support and access professional psychosocial support services that will meet their individual needs. The overall outcome of the service will be to ensure that the YPLWC visiting the site experience optimal psychological wellbeing. Ultimately, the service's value will be in improving the lives of young people who engage with it and the follow-on effect that this will have on their families and communities in the long-term.

Cardiac Ultrasound Findings in Infants with Severe (Hurler Phenotype) Untreated Mucopolysaccharidosis (MPS) Type I

BACKGROUND

Serious cardiac valve disease and left ventricular hypertrophy occur in most untreated older children with severe mucopolysaccharidosis type I. Although it is assumed that early intervention prevents these processes, evaluation of cardiac findings in these infants has not yet been reported.

METHODS

We reviewed echocardiograms of 13 untreated infants < 1 year of age with severe mucopolysaccharidosis type I who had undergone evaluation for hematopoietic cell transplantation. We recorded left ventricular chamber dimensions, septal and posterior wall thicknesses, ventricular function, and aortic sinus diameters. We evaluated mitral and aortic valves for increased thickness, regurgitation, and stenosis.

RESULTS

Average age (7M, 6F) was 221 (range 25-347) days. Left ventricular chamber dimension was ≥2 SD of normal in 3/13; wall thicknesses were ≥2 SD of normal in 2/13 infants. Systolic function was normal. Mitral valves were thickened in all infants; mitral regurgitation was present in 9/13, but significant in only three infants. Aortic valves were thickened in 10/13, but no infant had significant aortic regurgitation. Neither mitral nor aortic stenosis occurred. Aortic roots were dilated to ≥2 SD of normal in 5/13.

CONCLUSIONS

Characteristic cardiac features of severe mucopolysaccharidosis type I can be seen in infancy. Mitral and aortic valve thickening are nearly universally present, even in the youngest infants. In 20-30 % of infants, other abnormalities such as left ventricular dilation, increased wall thickness, and mild mitral/aortic regurgitation may occur. Aortic root dilation is a frequent finding. Early intervention with enzyme replacement therapy may minimize the incidence and severity of cardiac findings in these infants.

SUMMARY

Serious cardiac valve disease and left ventricular hypertrophy occur in most untreated older children with severe mucopolysaccharidosis type I. Although it is assumed that early intervention prevents these processes, evaluation of cardiac findings in these infants has not yet been reported. In our study of 13 infants with severe untreated MPS I < 1 year of age, mitral and aortic valve thickening was nearly universally present and aortic root dilation was frequent. Despite this, we found a lower incidence of left ventricular hypertrophy and both a lower incidence and milder expression of mitral and aortic valve dysfunction than previously reported in older children. These findings suggest that earlier intervention, including neonatal screening, may be of benefit to children with severe MPS I.

Research challenges in central nervous system manifestations of inborn errors of metabolism

The Research Challenges in CNS Manifestations of Inborn Errors of Metabolism workshop was designed to address challenges in translating potential therapies for these rare disorders, and to highlight novel therapeutic strategies and innovative approaches to CNS delivery, assessment of effects and directions for the future in the treatment of these diseases. Therapies for the brain in inborn errors represent some of the greatest challenges to translational research due to the special properties of the brain, and of inborn errors themselves. This review covers the proceedings of this workshop as submitted by participants. Scientific, ethical and regulatory issues are discussed, along with ways to measure outcomes and the conduct of clinical trials. Participants included regulatory and funding agencies, clinicians, scientists, industry and advocacy groups.

Engraftment and survival following hematopoietic stem cell transplantation for osteopetrosis using a reduced intensity conditioning regimen

Autosomal recessive osteopetrosis (OP) is a disease characterized by osteoclast dysfunction, leading to multisystem morbidity and death of most affected children. Hematopoietic stem cell transplantation (HSCT) is the treatment of choice for OP, but this patient population is particularly prone to post-transplant complications and death after myeloablative conditioning. To determine the potential of achieving improved overall outcomes in these patients by decreasing pre-transplant mortality, we investigated engraftment and survival following a reduced intensity regimen including busulfan, fludarabine and total lymphoid irradiation. We report outcomes in 11 patients. All six patients who received a bone marrow or peripheral stem cell graft engrafted with >75% donor chimerism. In contrast, all five recipients of unrelated cord blood as a stem cell source for a first graft failed to demonstrate donor hematopoietic chimerism. The day 100 and 6-month mortality was low at 9%. One year after HSCT, six of 11 patients (55%) were surviving. Our data suggest that this regimen results in low peri-transplant mortality without compromising engraftment when a marrow or peripheral stem cell graft is used. An umbilical cord blood graft, however, should be used with caution for patients with OP when this or a similar reduced intensity regimen is used.

Continued neurocognitive development and prevention of cardiopulmonary complications after successful BMT for I-cell disease: a long-term follow-up report

I-cell disease or mucolipidosis type II, a rare inherited storage disorder of lysosomal enzyme localization, is characterized by dysostosis multiplex, progressive severe psychomotor retardation and death by 5-8 years from congestive heart failure and recurrent pulmonary infections. A 19-month old girl with I-cell disease received a bone marrow transplant (BMT) from an HLA-identical carrier brother. At the age of 7 years, 5 years after BMT, she has no history of respiratory infections. Her cardiac function remains normal with a shortening fraction of 47%, and she continues to gain neurodevelopmental milestones, albeit at a very slow rate. Musculoskeletal deformities have worsened despite BMT. This is the first report describing neurodevelopmental gains and prevention of cardiopulmonary complications in I-cell disease after BMT.

The mutational spectrum of human malignant autosomal recessive osteopetrosis

Human malignant infantile osteopetrosis (arOP; MIM 259700) is a genetically heterogeneous autosomal recessive disorder of bone metabolism, which, if untreated, has a fatal outcome. Our group, as well as others, have recently identified mutations in the ATP6i (TCIRG1) gene, encoding the a3 subunit of the vacuolar proton pump, which mediates the acidification of the bone/osteoclast interface, are responsible for a subset of this condition. By sequencing the ATP6i gene in arOP patients from 44 unrelated families with a worldwide distribution we have now established that ATP6i mutations are responsible for approximately 50% of patients affected by this disease. The vast majority of these mutations (40 out of 42 alleles, including seven deletions, two insertions, 10 nonsense substitutions and 21 mutations in splice sites) are predicted to cause severe abnormalities in the protein product and are likely to represent null alleles. In addition, we have also analysed nine unrelated arOP patients from Costa Rica, where this disease is apparently much more frequent than elsewhere. All nine Costa Rican patients bore either or both of two missense mutations (G405R and R444L) in amino acid residues which are evolutionarily conserved from yeast to humans. The identification of ATP6i gene mutations in two families allowed us for the first time to perform prenatal diagnosis: both fetuses were predicted not to be affected and two healthy babies were born. This study contributes to the determination of genetic heterogeneity of arOP and allows further delineation of the other genetic defects causing this severe condition.

Transplantability and therapeutic effects of bone marrow-derived mesenchymal cells in children with osteogenesis imperfecta

In principle, transplantation of mesenchymal progenitor cells would attenuate or possibly correct genetic disorders of bone, cartilage and muscle, but clinical support for this concept is lacking. Here we describe the initial results of allogeneic bone marrow transplantation in three children with osteogenesis imperfecta, a genetic disorder in which osteoblasts produce defective type I collagen, leading to osteopenia, multiple fractures, severe bony deformities and considerably shortened stature. Three months after osteoblast engraftment (1.5-2.0% donor cells), representative specimens of trabecular bone showed histologic changes indicative of new dense bone formation. All patients had increases in total body bone mineral content ranging from 21 to 29 grams (median, 28), compared with predicted values of 0 to 4 grams (median, 0) for healthy children with similar changes in weight. These improvements were associated with increases in growth velocity and reduced frequencies of bone fracture. Thus, allogeneic bone marrow transplantation can lead to engraftment of functional mesenchymal progenitor cells, indicating the feasibility of this strategy in the treatment of osteogenesis imperfecta and perhaps other mesenchymal stem cell disorders as well.

Clinical importance of pre-morteum blood lymphocytes in cadaver donor tissue typing

We have refined our immunomagnetic bead (IM bead) procedures to isolate pure and viable lymphocyte subpopulation from pre-morteum (PM) blood for cadaver donor HLA typing, preliminary and final crossmatches (XMs). The results of 1220 XMs were compared using T/B lymphocytes isolated either from PM blood or spleen/lymphnode (SPLN) tissue. IM bead technique was used to isolate T/B cells from PM blood and nylon wool column (NWC) technique was used to isolate T/B cells from SPLN. When we compared the outcome of 800 T-cell crossmatches using T cells from PM blood or SPLN of 5 separate cadaver donors, NWC TXMs tended to be more falsenegative for high PRA (> 10%, total 500 XMs) as well as low PRA (< 10%, total 300 XMs) did not reach statistical significance. In contrast, NW BXM (420 B XM) were found to be far more false negative than IM bead BXM regardless of the PRA of the patients. In order to ensure that NWC BXMs were indeed false negative, 23 sera with known anti-DR antibodies were BXMed where antigen-specific B cells were isolated by both the techniques. Our results showed that IM bead BXM identified the DR specificities greater than 90% of the time, the titers of ab specificities were stronger (1:8). In comparison, NWB cell XMs were weak (titers 1:2), and the false negative rate for some ab was as high as 73%. Using IM bead and NWC techniques we compared our turnaround time (TAT) for cadaver donor typing, preliminary and final XMs.(ABSTRACT TRUNCATED AT 250 WORDS)

Primary hepatic lymphoma in an adolescent treated with hepatic lobectomy and chemotherapy

A case of B-cell non-Hodgkin's lymphoma, confined to the liver, in a 17-year-old boy is reported. The patient was treated with an extended left hepatectomy and combination chemotherapy: Cytoxan (cyclophosphamide), vincristine, prednisone, and methotrexate (COMP). The patient remains disease free at 1 year.