Global Ring Study to Investigate the Comparability of Total Assay Performance of Commercial Claudin 18 Antibodies for Evaluation in Gastric Cancer

Claudin 18.2 (CLDN18.2), the dominant isoform of CLDN18 in gastric tissues, is a highly specific tight junction protein of the gastric mucosa with variably retained expressions in gastric and gastroesophageal junction cancers. Additionally, CLDN18.2-targeted treatment with zolbetuximab, in combination with chemotherapy, has recently been assessed in 2 phase-III studies of patients with HER2-negative, locally advanced, unresectable, or metastatic gastric or gastroesophageal junction adenocarcinoma. These trials used the investigational VENTANA CLDN18 (43-14A) RxDx immunohistochemistry (IHC) assay on the Ventana BenchMark platform to identify patients eligible for CLDN18.2-targeted treatment. We report the findings of a global ring study evaluating the analytical comparability of concordance of the results of 3 CLDN18 antibodies (Ventana, LSBio, and Novus) stained on 3 IHC-staining platforms (Ventana, Dako, and Leica). A tissue microarray (TMA), comprising 15 gastric cancer cases, was stained by 27 laboratories across 11 countries. Each laboratory stained the TMAs using at least 2 of the 3 evaluated CLDN18 antibodies. Stained TMAs were assessed and scored using an agreed IHC-scoring algorithm, and the results were collated for statistical analysis. The data confirmed a high level of concordance for the VENTANA CLDN18 (43-14A; Ventana platform only) and LSBio antibodies on both the Dako and Leica platforms, with accuracy, precision, sensitivity, and specificity rates all reaching a minimum acceptable ≥85% threshold and good-to-excellent levels of concordance as measured by Cohen's kappa coefficient. The Novus antibody showed the highest level of variability against the reference central laboratory results for the same antibody/platform combinations. It also failed to meet the threshold for accuracy and sensitivity when used on either the Dako or Leica platform. These results demonstrated the reliability of IHC testing for CLDN18 expression in gastric tumor samples when using commercially available platforms with an appropriate methodology and primary antibody selection.

Transcription factor combinations that define human astrocyte identity encode significant variation of maturity and function

Increasing evidence indicates that cellular identity can be reduced to the distinct gene regulatory networks controlled by transcription factors (TFs). However, redundancy exists in these states as different combinations of TFs can induce broadly similar cell types. We previously demonstrated that by overcoming gene silencing, it is possible to deterministically reprogram human pluripotent stem cells directly into cell types of various lineages. In the present study we leverage the consistency and precision of our approach to explore four different TF combinations encoding astrocyte identity, based on previously published reports. Analysis of the resulting induced astrocytes (iAs) demonstrated that all four cassettes generate cells with the typical morphology of in vitro astrocytes, which expressed astrocyte-specific markers. The transcriptional profiles of all four iAs clustered tightly together and displayed similarities with mature human astrocytes, although maturity levels differed between cells. Importantly, we found that the TF cassettes induced iAs with distinct differences with regards to their cytokine response and calcium signaling. In vivo transplantation of selected iAs into immunocompromised rat brains demonstrated long term stability and integration. In conclusion, all four TF combinations were able to induce stable astrocyte-like cells that were morphologically similar but showed subtle differences with respect to their transcriptome. These subtle differences translated into distinct differences with regards to cell function, that could be related to maturation state and/or regional identity of the resulting cells. This insight opens an opportunity to precision-engineer cells to meet functional requirements, for example, in the context of therapeutic cell transplantation.

High-throughput sequencing of IgH gene in minor salivary glands from Sjögren's syndrome patients reveals dynamic B cell recirculation between ectopic lymphoid structures

OBJECTIVES

B cells play a central role in Sjögren's syndrome (SS) whereby autoreactive B-cells populate ectopic germinal centres (GC) in SS salivary glands (SG) and undergo somatic hypermutation (SHM) and class-switch recombination of the immunoglobulin genes. However, the capacity of specific B cell clones to seed ectopic GC in different SG and undergo clonal diversification is unclear. To unravel the dynamics of B cell recirculation among minor SG biopsies, we investigated the immunoglobulin heavy chain (IgH) gene usage and the pattern of SHM using a high-throughput sequencing approach.

METHODS

We generated ~166,000 reads longer than 350bp and detected 1631 clonotypes across eight samples from four different SS patients, all characterised by the presence of functional ectopic GC as demonstrated by the expression of activation-induced cytidine deaminase.

RESULTS

A large number of shared clonotypes were observed among paired mSG biopsies from each patient but not across different patients. Lineage tree analysis revealed significant clonal expansion within the mSG with the identification of shared dominant B cell clones suggestive of extensive recirculation across different SG. Several shared clonotypes with high proliferating capacity displayed IgH-VH gene usage common in autoreactive B cells, including VH1-69, which is typical of rheumatoid factor+ B cells representing potential lymphoma precursors.

CONCLUSIONS

The complex dynamic recirculation of B cells that we observed within ectopic GC responses linked with their ability to independently proliferate, undergo ongoing SHM and Ig class-switching within individual glands may explain the difficulty in achieving consistent eradication of ectopic GCs following B cell depleting agents reported in different studies.

Development and validation of resource-driven risk prediction models for incident chronic kidney disease in type 2 diabetes

Prediction models for population-based screening need, for global usage, to be resource-driven, involving predictors that are affordably resourced. Here, we report the development and validation of three resource-driven risk models to identify people with type 2 diabetes (T2DM) at risk of stage 3 CKD defined by a decline in estimated glomerular filtration rate (eGFR) to below 60 mL/min/1.73m2. The observational study cohort used for model development consisted of data from a primary care dataset of 20,510 multi-ethnic individuals with T2DM from London, UK (2007-2018). Discrimination and calibration of the resulting prediction models developed using cox regression were assessed using the c-statistic and calibration slope, respectively. Models were internally validated using tenfold cross-validation and externally validated on 13,346 primary care individuals from Wales, UK. The simplest model was simplified into a risk score to enable implementation in community-based medicine. The derived full model included demographic, laboratory parameters, medication-use, cardiovascular disease history (CVD) and sight threatening retinopathy status (STDR). Two less resource-intense models were developed by excluding CVD and STDR in the second model and HbA1c and HDL in the third model. All three 5-year risk models had good internal discrimination and calibration (optimism adjusted C-statistics were each 0.85 and calibration slopes 0.999-1.002). In Wales, models achieved excellent discrimination(c-statistics ranged 0.82-0.83). Calibration slopes at 5-years suggested models over-predicted risks, however were successfully updated to accommodate reduced incidence of stage 3 CKD in Wales, which improved their alignment with the observed rates in Wales (E/O ratios near to 1). The risk score demonstrated similar model performance compared to direct evaluation of the cox model. These resource-driven risk prediction models may enable universal screening for Stage 3 CKD to enable targeted early optimisation of risk factors for CKD.

Evaluation of the IgG antibody response to SARS CoV-2 infection and performance of a lateral flow immunoassay: cross-sectional and longitudinal analysis over 11 months

OBJECTIVE

To evaluate the dynamics and longevity of the humoral immune response to SARS-CoV-2 infection and assess the performance of professional use of the UK-RTC AbC-19 Rapid Test lateral flow immunoassay (LFIA) for the target condition of SARS-CoV-2 spike protein IgG antibodies.

DESIGN

Nationwide serological study.

SETTING

Northern Ireland, UK, May 2020-February 2021.

PARTICIPANTS

Plasma samples were collected from a diverse cohort of individuals from the general public (n=279), Northern Ireland healthcare workers (n=195), pre-pandemic blood donations and research studies (n=223) and through a convalescent plasma programme (n=183). Plasma donors (n=101) were followed with sequential samples over 11 months post-symptom onset.

MAIN OUTCOME MEASURES

SARS-CoV-2 antibody levels in plasma samples using Roche Elecsys Anti-SARS-CoV-2 IgG/IgA/IgM, Abbott SARS-CoV-2 IgG and EuroImmun IgG SARS-CoV-2 ELISA immunoassays over time. UK-RTC AbC-19 LFIA sensitivity and specificity, estimated using a three-reference standard system to establish a characterised panel of 330 positive and 488 negative SARS-CoV-2 IgG samples.

RESULTS

We detected persistence of SARS-CoV-2 IgG antibodies for up to 10 months post-infection, across a minimum of two laboratory immunoassays. On the known positive cohort, the UK-RTC AbC-19 LFIA showed a sensitivity of 97.58% (95.28% to 98.95%) and on known negatives, showed specificity of 99.59% (98.53 % to 99.95%).

CONCLUSIONS

Through comprehensive analysis of a cohort of pre-pandemic and pandemic individuals, we show detectable levels of IgG antibodies, lasting over 46 weeks when assessed by EuroImmun ELISA, providing insight to antibody levels at later time points post-infection. We show good laboratory validation performance metrics for the AbC-19 rapid test for SARS-CoV-2 spike protein IgG antibody detection in a laboratory-based setting.

ImmunoCluster provides a computational framework for the nonspecialist to profile high-dimensional cytometry data

High-dimensional cytometry is an innovative tool for immune monitoring in health and disease, and it has provided novel insight into the underlying biology as well as biomarkers for a variety of diseases. However, the analysis of large multiparametric datasets usually requires specialist computational knowledge. Here, we describe ImmunoCluster (https://github.com/kordastilab/ImmunoCluster), an R package for immune profiling cellular heterogeneity in high-dimensional liquid and imaging mass cytometry, and flow cytometry data, designed to facilitate computational analysis by a nonspecialist. The analysis framework implemented within ImmunoCluster is readily scalable to millions of cells and provides a variety of visualization and analytical approaches, as well as a rich array of plotting tools that can be tailored to users' needs. The protocol consists of three core computational stages: (1) data import and quality control; (2) dimensionality reduction and unsupervised clustering; and (3) annotation and differential testing, all contained within an R-based open-source framework.

Systems level profiling of chemotherapy-induced stress resolution in cancer cells reveals druggable trade-offs

Cancer cells can survive chemotherapy-induced stress, but how they recover from it is not known. Using a temporal multiomics approach, we delineate the global mechanisms of proteotoxic stress resolution in multiple myeloma cells recovering from proteasome inhibition. Our observations define layered and protracted programs for stress resolution that encompass extensive changes across the transcriptome, proteome, and metabolome. Cellular recovery from proteasome inhibition involved protracted and dynamic changes of glucose and lipid metabolism and suppression of mitochondrial function. We demonstrate that recovering cells are more vulnerable to specific insults than acutely stressed cells and identify the general control nonderepressable 2 (GCN2)-driven cellular response to amino acid scarcity as a key recovery-associated vulnerability. Using a transcriptome analysis pipeline, we further show that GCN2 is also a stress-independent bona fide target in transcriptional signature-defined subsets of solid cancers that share molecular characteristics. Thus, identifying cellular trade-offs tied to the resolution of chemotherapy-induced stress in tumor cells may reveal new therapeutic targets and routes for cancer therapy optimization.

Assessment of the conjunctival microcirculation for patients presenting with acute myocardial infarction compared to healthy controls

Microcirculatory dysfunction occurs early in cardiovascular disease (CVD) development. Acute myocardial infarction (MI) is a late consequence of CVD. The conjunctival microcirculation is readily-accessible for quantitative assessment and has not previously been studied in MI patients. We compared the conjunctival microcirculation of acute MI patients and age/sex-matched healthy controls to determine if there were differences in microcirculatory parameters. We acquired images using an iPhone 6s and slit-lamp biomicroscope. Parameters measured included diameter, axial velocity, wall shear rate and blood volume flow. Results are for all vessels as they were not sub-classified into arterioles or venules. The conjunctival microcirculation was assessed in 56 controls and 59 inpatients with a presenting diagnosis of MI. Mean vessel diameter for the controls was 21.41 ± 7.57 μm compared to 22.32 ± 7.66 μm for the MI patients (p < 0.001). Axial velocity for the controls was 0.53 ± 0.15 mm/s compared to 0.49 ± 0.17 mm/s for the MI patients (p < 0.001). Wall shear rate was higher for controls than MI patients (162 ± 93 s-1 vs 145 ± 88 s-1, p < 0.001). Blood volume flow did not differ significantly for the controls and MI patients (153 ± 124 pl/s vs 154 ± 125 pl/s, p = 0.84). This pilot iPhone and slit-lamp assessment of the conjunctival microcirculation found lower axial velocity and wall shear rate in patients with acute MI. Further study is required to correlate these findings further and assess long-term outcomes in this patient group with a severe CVD phenotype.

Circulating miRNAs as Biomarkers for Mitochondrial Neuro-Gastrointestinal Encephalomyopathy

Mitochondrial neurogastrointestinal encephalomyopathy (MNGIE) is an ultra-rare disease for which there are currently no validated outcome measures for assessing therapeutic intervention efficacy. The aim of this study was to identify a plasma and/or serum microRNA (miRNA) biomarker panel for MNGIE. Sixty-five patients and 65 age and sex matched healthy controls were recruited and assigned to one of four study phases: (i) discovery for sample size determination; (ii) candidate screening; (iii) candidate validation; and (iv) verifying the performance of the validated miRNA panel in four patients treated with erythrocyte-encapsulated thymidine phosphorylase (EE-TP), an enzyme replacement under development for MNGIE. Quantitative PCR (qPCR) was used to profile miRNAs in serum and/or plasma samples collected for the discovery, validation and performance phases, and next generation sequencing (NGS) analysis was applied to serum samples assigned to the candidate screening phase. Forty-one differentially expressed candidate miRNAs were identified in the sera of patients (p < 0.05, log₂ fold change > 1). The validation cohort revealed that of those, 27 miRNAs were upregulated in plasma and three miRNAs were upregulated in sera (p < 0.05). Through binary logistic regression analyses, five plasma miRNAs (miR-192-5p, miR-193a-5p, miR-194-5p, miR-215-5p and miR-34a-5p) and three serum miRNAs (miR-192-5p, miR-194-5p and miR-34a-5p) were shown to robustly distinguish MNGIE from healthy controls. Reduced longitudinal miRNA expression of miR-34a-5p was observed in all four patients treated with EE-TP and coincided with biochemical and clinical improvements. We recommend the inclusion of the plasma exploratory miRNA biomarker panel in future clinical trials of investigational therapies for MNGIE; it may have prognostic value for assessing clinical status.

Calcium calmodulin kinase II activity is required for cartilage homeostasis in osteoarthritis

WNT ligands can activate several signalling cascades of pivotal importance during development and regenerative processes. Their de-regulation has been associated with the onset of different diseases. Here we investigated the role of the WNT/Calcium Calmodulin Kinase II (CaMKII) pathway in osteoarthritis. We identified Heme Oxygenase I (HMOX1) and Sox-9 as specific markers of the WNT/CaMKII signalling in articular chondrocytes through a microarray analysis. We showed that the expression of the activated form of CaMKII, phospho-CaMKII, was increased in human and murine osteoarthritis and the expression of HMOX1 was accordingly reduced, demonstrating the activation of the pathway during disease progression. To elucidate its function, we administered the CaMKII inhibitor KN93 to mice in which osteoarthritis was induced by resection of the anterior horn of the medial meniscus and of the medial collateral ligament in the knee joint. Pharmacological blockade of CaMKII exacerbated cartilage damage and bone remodelling. Finally, we showed that CaMKII inhibition in articular chondrocytes upregulated the expression of matrix remodelling enzymes alone and in combination with Interleukin 1. These results suggest an important homeostatic role of the WNT/CaMKII signalling in osteoarthritis which could be exploited in the future for therapeutic purposes.

Identification of genetic factors controlling phosphorus utilization efficiency in wheat by genome-wide association study with principal component analysis

Despite the economic importance of P utilization efficiency, information on genetic factors underlying this trait remains elusive. To address that, we performed a genome-wide association study in a spring wheat diversity panel ranging from landraces to elite varieties. We evaluated the phenotype variation for P utilization efficiency in controlled conditions and genotype variation using wheat 90 K SNP array. Phenotype variables were transformed into a smaller set of uncorrelated principal components that captured the most important variation data. We identified two significant loci associated with both P utilization efficiency and the 1st principal component on chromosomes 3A and 4A: qPE1-3A and qPE2-4A. Annotation of genes at these loci revealed 53 wheat genes, among which 6 were identified in significantly enriched pathways. The expression pattern of these 6 genes indicated that TraesCS4A02G481800, involved in pyruvate metabolism and TCA cycle, had a significantly higher expression in the P efficient variety under limited P conditions. Further characterization of these loci and candidate genes can help stimulate P utilization efficiency in wheat.

Molecular Portraits of Early Rheumatoid Arthritis Identify Clinical and Treatment Response Phenotypes

There is a current imperative to unravel the hierarchy of molecular pathways that drive the transition of early to established disease in rheumatoid arthritis (RA). Herein, we report a comprehensive RNA sequencing analysis of the molecular pathways that drive early RA progression in the disease tissue (synovium), comparing matched peripheral blood RNA-seq in a large cohort of early treatment-naive patients, namely, the Pathobiology of Early Arthritis Cohort (PEAC). We developed a data exploration website (https://peac.hpc.qmul.ac.uk/) to dissect gene signatures across synovial and blood compartments, integrated with deep phenotypic profiling. We identified transcriptional subgroups in synovium linked to three distinct pathotypes: fibroblastic pauci-immune pathotype, macrophage-rich diffuse-myeloid pathotype, and a lympho-myeloid pathotype characterized by infiltration of lymphocytes and myeloid cells. This is suggestive of divergent pathogenic pathways or activation disease states. Pro-myeloid inflammatory synovial gene signatures correlated with clinical response to initial drug therapy, whereas plasma cell genes identified a poor prognosis subgroup with progressive structural damage.

•

Deep phenotyping and RNA-seq of early rheumatoid arthritis individuals pre-treatment

•

Synovial plasma cell gene expression predicts future progressive joint damage on X-ray

•

Blood interferon gene signature associates with synovial B and plasma cell infiltration

•

Interactive website enables RNA-seq and clinical data to be fully explored

Mutation-Independent Allele-Specific Editing by CRISPR-Cas9, a Novel Approach to Treat Autosomal Dominant Disease

CRISPR-Cas9 provides a tool to treat autosomal dominant disease by non-homologous end joining (NHEJ) gene disruption of the mutant allele. In order to discriminate between wild-type and mutant alleles, Streptococcus pyogenes Cas9 (SpCas9) must be able to detect a single nucleotide change. Allele-specific editing can be achieved by using either a guide-specific approach, in which the missense mutation is found within the guide sequence, or a protospacer-adjacent motif (PAM)-specific approach, in which the missense mutation generates a novel PAM. While both approaches have been shown to offer allele specificity in certain contexts, in cases where numerous missense mutations are associated with a particular disease, such as TGFBI (transforming growth factor β-induced) corneal dystrophies, it is neither possible nor realistic to target each mutation individually. In this study, we demonstrate allele-specific CRISPR gene editing independent of the disease-causing mutation that is capable of achieving complete allele discrimination, and we propose it as a targeting approach for autosomal dominant disease. Our approach utilizes natural variants in the target region that contain a PAM on one allele that lies in cis with the causative mutation, removing the constraints of a mutation-dependent approach. Our innovative patient-specific guide design approach takes into account the patient's individual genetic make-up, allowing on- and off-target activity to be assessed in a personalized manner.

Comprehensive transcriptomic analysis of papillary thyroid cancer: potential biomarkers associated with tumor progression

PURPOSE

Identification of stage-specific prognostic/predictive biomarkers in papillary thyroid carcinoma (PTC) could lead to its more efficient clinical management. The main objective of this study was to characterize the stage-specific deregulation in genes and miRNA expression in PTC to identify potential prognostic biomarkers.

METHODS

495 RNASeq and 499 miRNASeq PTC samples (stage I-IV) as well as, respectively, 56 and 57 normal samples were retrieved from The Cancer Genome Atlas (TCGA). Differential expression analysis was performed using DESeq 2 to identify deregulation of genes and miRNAs between sequential stages. To identify the minority of patients who progress to higher stages, we performed clustering analysis on stage I RNASeq data. An independent PTC RNASeq data set (BioProject accession PRJEB11591) was also used for the validation of the results.

RESULTS

LTF and PLA2R1 were identified as two promising biomarkers down-regulated in a subgroup of stage I (both in TCGA and in the validation data set) and in the majority of stage IV of PTC (in TCGA data set). hsa-miR-205, hsa-miR-509-2, hsa-miR-514-1 and hsa-miR-514-2 were also detected as up-regulated miRNAs in both PTC patients with stage I and stage III. Hierarchical clustering of stage I samples showed substantial heterogeneity in the expression pattern of PTC indicating the necessity of categorizing stage I patients based on the expressional alterations of specific biomarkers.

CONCLUSION

Stage I PTC patients showed large amount of expressional heterogeneity. Therefore, risk stratification based on the expressional alterations of candidate biomarkers could be an important step toward personalized management of these patients.

The T cell differentiation landscape is shaped by tumour mutations in lung cancer

Tumour mutational burden (TMB) predicts immunotherapy outcome in non-small cell lung cancer (NSCLC), consistent with immune recognition of tumour neoantigens. However, persistent antigen exposure is detrimental for T cell function. How TMB affects CD4 and CD8 T cell differentiation in untreated tumours, and whether this affects patient outcomes is unknown. Here we paired high-dimensional flow cytometry, exome, single-cell and bulk RNA sequencing from patients with resected, untreated NSCLC to examine these relationships. TMB was associated with compartment-wide T cell differentiation skewing, characterized by loss of TCF7-expressing progenitor-like CD4 T cells, and an increased abundance of dysfunctional CD8 and CD4 T cell subsets, with significant phenotypic and transcriptional similarity to neoantigen-reactive CD8 T cells. A gene signature of redistribution from progenitor-like to dysfunctional states associated with poor survival in lung and other cancer cohorts. Single-cell characterization of these populations informs potential strategies for therapeutic manipulation in NSCLC.

Fish oil supplementation during pregnancy is protective against asthma/wheeze in offspring

The trial was approved by the Institutional Review Boards of the participating institutions and at Brigham and Women’s Hospital and is registered with NCT00920621 ). Written informed consent was obtained from all women at recruitment.

Lactate Buildup at the Site of Chronic Inflammation Promotes Disease by Inducing CD4+ T Cell Metabolic Rewiring

Accumulation of lactate in the tissue microenvironment is a feature of both inflammatory disease and cancer. Here, we assess the response of immune cells to lactate in the context of chronic inflammation. We report that lactate accumulation in the inflamed tissue contributes to the upregulation of the lactate transporter SLC5A12 by human CD4+ T cells. SLC5A12-mediated lactate uptake into CD4+ T cells induces a reshaping of their effector phenotype, resulting in increased IL17 production via nuclear PKM2/STAT3 and enhanced fatty acid synthesis. It also leads to CD4+ T cell retention in the inflamed tissue as a consequence of reduced glycolysis and enhanced fatty acid synthesis. Furthermore, antibody-mediated blockade of SLC5A12 ameliorates the disease severity in a murine model of arthritis. Finally, we propose that lactate/SLC5A12-induced metabolic reprogramming is a distinctive feature of lymphoid synovitis in rheumatoid arthritis patients and a potential therapeutic target in chronic inflammatory disorders.

The role of the 17q21 genotype in the prevention of early childhood asthma and recurrent wheeze by vitamin D

Evidence suggests vitamin D has preventive potential in asthma; however, not all children benefit from this intervention. This study aimed to investigate whether variation in the functional 17q21 single nucleotide polymorphism rs12936231 affects the preventive potential of vitamin D against asthma.A combined secondary analysis of two randomised controlled trials of prenatal vitamin D supplementation for the prevention of asthma in offspring (Vitamin D Antenatal Asthma Reduction Trial (VDAART) and Copenhagen Prospective Studies on Asthma in Childhood 2010 (COPSAC₂₀₁₀)) was performed, stratifying by genotype and integrating metabolite data to explore underlying mechanisms.The protective effect of vitamin D on asthma/wheeze was evident among children with the low-risk rs12936231 GG genotype (hazard ratio (HR) 0.49, 95% CI 0.26-0.94, p=0.032) but not the high-risk CC genotype (HR 1.08, 95% CI 0.69-1.69, p=0.751). In VDAART, in the GG genotype vitamin D supplementation was associated with increased plasma levels of sphingolipids, including sphingosine-1-phosphate (β 0.022, 95% CI 0.001-0.044, p=0.038), but this was not evident with the CC genotype, known to be associated with increased expression of ORMDL3 in bronchial epithelial cells. Sphingolipid levels were associated with decreased risk of asthma/wheeze, and there was evidence of interactions between sphingolipid levels, vitamin D and genotype (p-interaction_{vitaminD*genotype*sphingosine-1-phosphate}=0.035). In a cellular model, there was a significant difference in the induction of sphingosine-1-phosphate by vitamin D between a control human bronchial epithelial cell line and a cell line overexpressing ORMDL3 (p=0.002).Results suggest prenatal vitamin D supplementation may reduce the risk of early childhood asthma/wheeze via alterations of sphingolipid metabolism dependent on the 17q21 genotype.

Gene editing in the context of an increasingly complex genome

The reporting of the first draft of the human genome in 2000 brought with it much hope for the future in what was felt as a paradigm shift toward improved health outcomes. Indeed, we have now mapped the majority of variation across human populations with landmark projects such as 1000 Genomes; in cancer, we have catalogued mutations across the primary carcinomas; whilst, for other diseases, we have identified the genetic variants with strongest association. Despite this, we are still awaiting the genetic revolution in healthcare to materialise and translate itself into the health benefits for which we had hoped. A major problem we face relates to our underestimation of the complexity of the genome, and that of biological mechanisms, generally. Fixation on DNA sequence alone and a 'rigid' mode of thinking about the genome has meant that the folding and structure of the DNA molecule -and how these relate to regulation- have been underappreciated. Projects like ENCODE have additionally taught us that regulation at the level of RNA is just as important as that at the spatiotemporal level of chromatin.In this review, we chart the course of the major advances in the biomedical sciences in the era pre- and post the release of the first draft sequence of the human genome, taking a focus on technology and how its development has influenced these. We additionally focus on gene editing via CRISPR/Cas9 as a key technique, in particular its use in the context of complex biological mechanisms. Our aim is to shift the mode of thinking about the genome to that which encompasses a greater appreciation of the folding of the DNA molecule, DNA- RNA/protein interactions, and how these regulate expression and elaborate disease mechanisms.Through the composition of our work, we recognise that technological improvement is conducive to a greater understanding of biological processes and life within the cell. We believe we now have the technology at our disposal that permits a better understanding of disease mechanisms, achievable through integrative data analyses. Finally, only with greater understanding of disease mechanisms can techniques such as gene editing be faithfully conducted.

An Integrative Transcriptomic and Metabolomic Study of Lung Function in Children With Asthma

BACKGROUND

Single omic analyses have provided some insight into the basis of lung function in children with asthma, but the underlying biologic pathways are still poorly understood.

METHODS

Weighted gene coexpression network analysis (WGCNA) was used to identify modules of coregulated gene transcripts and metabolites in blood among 325 children with asthma from the Genetic Epidemiology of Asthma in Costa Rica study. The biology of modules associated with lung function as measured by FEV₁, the FEV₁/FVC ratio, bronchodilator response, and airway responsiveness to methacholine was explored. Significantly correlated gene-metabolite module pairs were then identified, and their constituent features were analyzed for biologic pathway enrichments.

RESULTS

WGCNA clustered 25,060 gene probes and 8,185 metabolite features into eight gene modules and eight metabolite modules, where four and six, respectively, were associated with lung function (P ≤ .05). The gene modules were enriched for immune, mitotic, and metabolic processes and asthma-associated microRNA targets. The metabolite modules were enriched for lipid and amino acid metabolism. Integration of correlated gene-metabolite modules expanded the single omic findings, linking the FEV₁/FVC ratio with ORMDL3 and dysregulated lipid metabolism. This finding was replicated in an independent population.

CONCLUSIONS

The results of this hypothesis-generating study suggest a mechanistic basis for multiple asthma genes, including ORMDL3, and a role for lipid metabolism. They demonstrate that integrating multiple omic technologies may provide a more informative picture of asthmatic lung function biology than single omic analyses.

Racial differences in endometrial cancer molecular portraits in The Cancer Genome Atlas

Endometrial cancer (EC) is now the most prevalent gynaecological malignancy in the Western world. Black or African American women (BoAA) have double the mortality of Caucasian women, and their tumours tend to be of higher grade. Despite these disparities, little is known regarding the mutational landscape of EC between races. Hence, we wished to investigate the molecular features of ECs within The Cancer Genome Atlas (TCGA) dataset by racial groupings. In total 374 Caucasian, 109 BoAA and 20 Asian patients were included in the analysis. Asian women were diagnosed at younger age, 54.2 years versus 64.5 years for Caucasian and 64.9 years for BoAA women (OR 3.432; p=0.011); BoAA women were more likely to have serous type tumors (OR 2.061; p=0.008). No difference in overall survival was evident. The most frequently mutated gene in Caucasian and Asian tumours was PTEN (63% and 85%), unlike BoAA cases where it was TP53 (49%). Mutation and somatic copy number alteration (SCNA) analysis revealed an enrichment of TP53 mutations in BoAAs; whereas POLE and RPL22 mutations were more frequent in Caucasians. Major recurrent SCNA racial differences were observed at chromosomes 3p, 8, 10, and 16, which clustered BoAA tumors into 4 distinct groups and Caucasian tumors into 5 groups. There was a significantly higher frequency of somatic mutations in DNA mismatch repair genes in Asian tumours, in particular PMS2 (p=0.0036). In conclusion, inherent racial disparities appear to be present in the molecular profile of EC, which could have potential implications on clinical management.

Vitamin D prenatal programming of childhood metabolomics profiles at age 3 y

Background: Vitamin D deficiency is implicated in a range of common complex diseases that may be prevented by gestational vitamin D repletion. Understanding the metabolic mechanisms related to in utero vitamin D exposure may therefore shed light on complex disease susceptibility.Objective: The goal was to analyze the programming role of in utero vitamin D exposure on children's metabolomics profiles.Design: First, unsupervised clustering was done with plasma metabolomics profiles from a case-control subset of 245 children aged 3 y with and without asthma from the Vitamin D Antenatal Asthma Reduction Trial (VDAART), in which pregnant women were randomly assigned to vitamin D supplementation or placebo. Thereafter, we analyzed the influence of maternal pre- and postsupplement vitamin D concentrations on cluster membership. Finally, we used the metabolites driving the clustering of children to identify the dominant metabolic pathways that were influential in each cluster.Results: We identified 3 clusters of children characterized by 1) high concentrations of fatty acids and amines and low maternal postsupplement vitamin D (mean ± SD; 27.5 ± 11.0 ng/mL), 2) high concentrations of amines, moderate concentrations of fatty acids, and normal maternal postsupplement vitamin D (34.0 ± 14.1 ng/mL), and 3) low concentrations of fatty acids, amines, and normal maternal postsupplement vitamin D (35.2 ± 15.9 ng/mL). Adjusting for sample storage time, maternal age and education, and both child asthma and vitamin D concentration at age 3 y did not modify the association between maternal postsupplement vitamin D and cluster membership (P = 0.0014). Maternal presupplement vitamin D did not influence cluster membership, whereas the combination of pre- and postsupplement concentrations did (P = 0.03).Conclusions: Young children can be clustered into distinct biologically meaningful groups by their metabolomics profiles. The clusters differed in concentrations of inflammatory mediators, and cluster membership was influenced by in utero vitamin D exposure, suggesting a prenatal programming role of vitamin D on the child's metabolome. This trial was registered at clinicaltrials.gov as NCT00920621.

Applications of Metabolomics in the Study and Management of Preeclampsia; A Review of the Literature

Introduction

Preeclampsia represents a major public health burden worldwide, but predictive and diagnostic biomarkers are lacking. Metabolomics is emerging as a valuable approach to generating novel biomarkers whilst increasing the mechanistic understanding of this complex condition.

Objectives

To summarize the published literature on the use of metabolomics as a tool to study preeclampsia.

Methods

PubMed and Web of Science were searched for articles that performed metabolomic profiling of human biosamples using either Mass-spectrometry or Nuclear Magnetic Resonance based approaches and which included preeclampsia as a primary endpoint.

Results

Twenty-eight studies investigating the metabolome of preeclampsia in a variety of biospecimens were identified. Individual metabolite and metabolite profiles were reported to have discriminatory ability to distinguish preeclamptic from normal pregnancies, both prior to and post diagnosis. Lipids and carnitines were among the most commonly reported metabolites. Further work and validation studies are required to demonstrate the utility of such metabolites as preeclampsia biomarkers.

Conclusion

Metabolomic-based biomarkers of preeclampsia have yet to be integrated into routine clinical practice. However, metabolomic profiling is becoming increasingly popular in the study of preeclampsia and is likely to be a valuable tool to better understand the pathophysiology of this disorder and to better classify its subtypes, particularly when integrated with other omic data.

Lineage-Specific Genes Are Prominent DNA Damage Hotspots during Leukemic Transformation of B Cell Precursors

In human leukemia, lineage-specific genes represent predominant targets of deletion, with lymphoid-specific genes frequently affected in lymphoid leukemia and myeloid-specific genes in myeloid leukemia. To investigate the basis of lineage-specific alterations, we analyzed global DNA damage in primary B cell precursors expressing leukemia-inducing oncogenes by ChIP-seq. We identified more than 1,000 sensitive regions, of which B lineage-specific genes constitute the most prominent targets. Identified hotspots at B lineage genes relate to DNA-DSBs, affect genes that harbor genomic lesions in human leukemia, and associate with ectopic deletion in successfully transformed cells. Furthermore, we show that most identified regions overlap with gene bodies of highly expressed genes and that induction of a myeloid lineage phenotype in transformed B cell precursors promotes de novo DNA damage at myeloid loci. Hence, we demonstrate that lineage-specific transcription predisposes lineage-specific genes in transformed B cell precursors to DNA damage, which is likely to promote the frequent alteration of lineage-specific genes in human leukemia.

KSR1 regulates BRCA1 degradation and inhibits breast cancer growth

Kinase suppressor of Ras-1 (KSR1) facilitates signal transduction in Ras-dependent cancers, including pancreatic and lung carcinomas but its role in breast cancer has not been well studied. Here, we demonstrate for the first time it functions as a tumor suppressor in breast cancer in contrast to data in other tumors. Breast cancer patients (n>1000) with high KSR1 showed better disease-free and overall survival, results also supported by Oncomine analyses, microarray data (n=2878) and genomic data from paired tumor and cell-free DNA samples revealing loss of heterozygosity. KSR1 expression is associated with high breast cancer 1, early onset (BRCA1), high BRCA1-associated ring domain 1 (BARD1) and checkpoint kinase 1 (Chk1) levels. Phospho-profiling of major components of the canonical Ras-RAF-mitogen-activated protein kinases pathway showed no significant changes after KSR1 overexpression or silencing. Moreover, KSR1 stably transfected cells formed fewer and smaller size colonies compared to the parental ones, while in vivo mouse model also demonstrated that the growth of xenograft tumors overexpressing KSR1 was inhibited. The tumor suppressive action of KSR1 is BRCA1 dependent shown by 3D-matrigel and soft agar assays. KSR1 stabilizes BRCA1 protein levels by reducing BRCA1 ubiquitination through increasing BARD1 abundance. These data link these proteins in a continuum with clinical relevance and position KSR1 in the major oncoprotein pathways in breast tumorigenesis.

Whole genome sequence analysis suggests intratumoral heterogeneity in dissemination of breast cancer to lymph nodes

BACKGROUND

Intratumoral heterogeneity may help drive resistance to targeted therapies in cancer. In breast cancer, the presence of nodal metastases is a key indicator of poorer overall survival. The aim of this study was to identify somatic genetic alterations in early dissemination of breast cancer by whole genome next generation sequencing (NGS) of a primary breast tumor, a matched locally-involved axillary lymph node and healthy normal DNA from blood.

METHODS

Whole genome NGS was performed on 12 µg (range 11.1-13.3 µg) of DNA isolated from fresh-frozen primary breast tumor, axillary lymph node and peripheral blood following the DNA nanoball sequencing protocol. Single nucleotide variants, insertions, deletions, and substitutions were identified through a bioinformatic pipeline and compared to CIN25, a key set of genes associated with tumor metastasis.

RESULTS

Whole genome sequencing revealed overlapping variants between the tumor and node, but also variants that were unique to each. Novel mutations unique to the node included those found in two CIN25 targets, TGIF2 and CCNB2, which are related to transcription cyclin activity and chromosomal stability, respectively, and a unique frameshift in PDS5B, which is required for accurate sister chromatid segregation during cell division. We also identified dominant clonal variants that progressed from tumor to node, including SNVs in TP53 and ARAP3, which mediates rearrangements to the cytoskeleton and cell shape, and an insertion in TOP2A, the expression of which is significantly associated with tumor proliferation and can segregate breast cancers by outcome.

CONCLUSION

This case study provides preliminary evidence that primary tumor and early nodal metastasis have largely overlapping somatic genetic alterations. There were very few mutations unique to the involved node. However, significant conclusions regarding early dissemination needs analysis of a larger number of patient samples.

Influence of plasma processing on recovery and analysis of circulating nucleic acids

Circulating nucleic acids (CNAs) are under investigation as a liquid biopsy in cancer. However there is wide variation in blood processing and methods for isolation of circulating free DNA (cfDNA) and microRNAs (miRNAs). Here we compare the extraction efficiency and reproducibility of 4 commercially available kits for cfDNA and 3 for miRNA using spike-in of reference templates. We also compare the effects of increasing time between venepuncture and centrifugation and differential centrifugation force on recovery of CNAs. cfDNA was quantified by TaqMan qPCR and targeted deep sequencing. miRNA profiles were assessed with TaqMan low-density arrays and assays. The QIAamp(®) DNA Blood Mini and Circulating nucleic acid kits gave the highest recovery of cfDNA and efficient recovery (>90%) of a 564bp spike-in. Moreover, targeted sequencing revealed overlapping cfDNA profiles and variant depth, including detection of HER2 gene amplification, using the Ion AmpliSeq™Cancer Hotspot Panel v2. Highest yields of miRNA and the synthetic Arabidopsis thaliana miR-159a spike-in were obtained using the miRNeasy Serum/Plasma kit, with saturation above 200 µl of plasma. miRNA profiles showed significant variation with increasing time before centrifugation (p<0.001) and increasing centrifugation force, with depletion of platelet associated miRNAs, whereas cfDNA was unaffected. However, sample replicates showed excellent reproducibility on TaqMan low density arrays (ρ = 0.96, p<0.0001). We also successfully generated miRNA profiles for plasma samples stored > 12 years, highlighting the potential for analysis of stored sample biobanks. In the era of the liquid biopsy, standardisation of methods is required to minimise variation, particularly for miRNA.

LMTK3 is implicated in endocrine resistance via multiple signaling pathways

Resistance to endocrine therapy in breast cancer is common. With the aim of discovering new molecular targets for breast cancer therapy, we have recently identified LMTK3 as a regulator of the estrogen receptor-alpha (ERα) and wished to understand its role in endocrine resistance. We find that inhibition of LMTK3 in a xenograft tamoxifen (Tam)-resistant (BT474) breast cancer mouse model results in re-sensitization to Tam as demonstrated by a reduction in tumor volume. A whole genome microarray analysis, using a BT474 cell line, reveals genes significantly modulated (positively or negatively) after LMTK3 silencing, including some that are known to be implicated in Tam resistance, notably c-MYC, HSPB8 and SIAH2. We show that LMTK3 is able to increase the levels of HSPB8 at a transcriptional and translational level thereby protecting MCF7 cells from Tam-induced cell death, by reducing autophagy. Finally, high LMTK3 levels at baseline in tumors are predictive for endocrine resistance; therapy does not lead to alteration in levels, whereas in patient's plasma samples, acquired LMTK3 gene amplification (copy number variation) was associated with relapse while receiving Tam. In aggregate, these data support a role for LMTK3 in both innate (intrinsic) and acquired (adaptive) endocrine resistance in breast cancer.

Hide and seek: tell-tale signs of breast cancer lurking in the blood

Breast cancer treatment is improving due to the introduction of new drugs, guided by molecular testing of the primary tumour for mutations/oncogenic drivers (e.g. HER2 gene amplification). However, tumour tissue is not always available for molecular analysis, intra-tumoural heterogeneity is common and the "cancer genome" is known to evolve with time, particularly following treatment as resistance develops. After resection, those patients with only residual micrometastases are likely to be cured but those with radiologically detectable overt disease are not. Thus, the discovery of blood test(s) that could (1) alert clinicians to early primary or recurrent disease and (2) monitor response to treatment could impact significantly on mortality. Towards this, we and others have focused on molecular profiling of circulating nucleic acids isolated from plasma, both cell-free DNA (cfDNA) and microRNAs, and the relationship of these to circulating tumour cells (CTCs). This review considers the utility of each as circulating biomarkers in breast cancer with particular emphasis on the bioinformatic tools available to support molecular profiling.

Genomic analysis of circulating cell-free DNA infers breast cancer dormancy

Biomarkers in breast cancer to monitor minimal residual disease have remained elusive. We hypothesized that genomic analysis of circulating free DNA (cfDNA) isolated from plasma may form the basis for a means of detecting and monitoring breast cancer. We profiled 251 genomes using Affymetrix SNP 6.0 arrays to determine copy number variations (CNVs) and loss of heterozygosity (LOH), comparing 138 cfDNA samples with matched primary tumor and normal leukocyte DNA in 65 breast cancer patients and eight healthy female controls. Concordance of SNP genotype calls in paired cfDNA and leukocyte DNA samples distinguished between breast cancer patients and healthy female controls (P < 0.0001) and between preoperative patients and patients on follow-up who had surgery and treatment (P = 0.0016). Principal component analyses of cfDNA SNP/copy number results also separated presurgical breast cancer patients from the healthy controls, suggesting specific CNVs in cfDNA have clinical significance. We identified focal high-level DNA amplification in paired tumor and cfDNA clustered in a number of chromosome arms, some of which harbor genes with oncogenic potential, including USP17L2 (DUB3), BRF1, MTA1, and JAG2. Remarkably, in 50 patients on follow-up, specific CNVs were detected in cfDNA, mirroring the primary tumor, up to 12 yr after diagnosis despite no other evidence of disease. These data demonstrate the potential of SNP/CNV analysis of cfDNA to distinguish between patients with breast cancer and healthy controls during routine follow-up. The genomic profiles of cfDNA infer dormancy/minimal residual disease in the majority of patients on follow-up.

The presence of disseminated tumour cells in the bone marrow is inversely related to circulating free DNA in plasma in breast cancer dormancy

BACKGROUND

The aim of this study was to gain insight into breast cancer dormancy by examining different measures of minimal residual disease (MRD) over time in relation to known prognostic factors.

METHODS

Sixty-four primary breast cancer patients on follow-up (a median of 8.3 years post surgery) who were disease free had sequential bone marrow aspirates and blood samples taken for the measurement of disseminated tumour cells (DTCs), circulating tumour cells (CTCs) by CellSearch and qPCR measurement of overlapping (96-bp and 291-bp) amplicons in circulating free DNA (cfDNA).

RESULTS

The presence of CTCs was correlated with the presence of DTCs measured by immunocytochemistry (P=0.01) but both were infrequently detected. Increasing cfDNA concentration correlated with ER, HER2 and triple-negative tumours and high tumour grade, and the 291-bp amplicon was inversely correlated with DTCs measured by CK19 qRT-PCR (P=0.047).

CONCLUSION

Our results show that breast cancer patients have evidence of MRD for many years after diagnosis despite there being no overt evidence of disease. The inverse relationship between bone marrow CK19 mRNA and the 291-bp amplicon in cfDNA suggests that an inverse relationship between a measure of cell viability in the bone marrow (DTCs) and cell death in the plasma occurs during the dormancy phase of breast cancer.

Detection of HER2 amplification in circulating free DNA in patients with breast cancer

BACKGROUND

Human epidermal growth factor receptor 2 (HER2) is amplified and overexpressed in 20-25% of breast cancers. This study investigated circulating free DNA (cfDNA) for detection of HER2 gene amplification in patients with breast cancer.

METHODS

Circulating free DNA was extracted from plasma of unselected patients with primary breast cancer (22 before surgery and 68 following treatment), 30 metastatic patients and 98 female controls using the QIAamp Blood DNA Mini Kit (Qiagen). The ratio of HER2 to an unamplified reference gene (contactin-associated protein 1 (CNTNAP1)) was measured in cfDNA samples by quantitative PCR (qPCR) using SK-BR-3 cell line DNA as a positive control.

RESULTS

We validated the qPCR assay with DNA extracted from 23 HER2 3+ and 40 HER2-negative tumour tissue samples; the results agreed for 60 of 63 (95.2%) tumours. Amplification was detected in cfDNA for 8 of 68 patients following primary breast cancer treatment and 5 of 30 metastatic patients, but was undetected in 22 patients with primary breast cancer and 98 healthy female controls. Of the patients with amplification in cfDNA, 10 had HER2 3+ tumour status by immunohistochemistry.

CONCLUSIONS

The results demonstrate for the first time the existence of amplified HER2 in cfDNA in the follow-up of breast cancer patients who are otherwise disease free. This approach could potentially provide a marker in patients with HER2-positive breast cancer.