iBBiG: iterative binary bi-clustering of gene sets

Motivation: Meta-analysis of genomics data seeks to identify genes associated with a biological phenotype across multiple datasets; however, merging data from different platforms by their features (genes) is challenging. Meta-analysis using functionally or biologically characterized gene sets simplifies data integration is biologically intuitive and is seen as having great potential, but is an emerging field with few established statistical methods. Results: We transform gene expression profiles into binary gene set profiles by discretizing results of gene set enrichment analyses and apply a new iterative bi-clustering algorithm (iBBiG) to identify groups of gene sets that are coordinately associated with groups of phenotypes across multiple studies. iBBiG is optimized for meta-analysis of large numbers of diverse genomics data that may have unmatched samples. It does not require prior knowledge of the number or size of clusters. When applied to simulated data, it outperforms commonly used clustering methods, discovers overlapping clusters of diverse sizes and is robust in the presence of noise. We apply it to meta-analysis of breast cancer studies, where iBBiG extracted novel gene set—phenotype association that predicted tumor metastases within tumor subtypes

Functional classification analysis of somatically mutated genes in human breast and colorectal cancers

A recent study published by Sjoblom and colleagues performed comprehensive sequencing of 13,023 human genes and identified mutations in genes specific to breast and colorectal tumors, providing insight into organ-specific tumor biology. Here we present a systematic analysis of the functional classifications of Sjoblom’s “CAN” genes, a subset of these validated mutant genes that identify novel organ-specific biological themes and molecular pathways associated with diseasespecific etiology. This analysis links four somatically mutated genes associated with diverse oncological types to colorectal and breast cancers through established TGF-β1 regulated interactions, revealing mechanistic differences in these cancers and providing potential diagnostic and therapeutic targets

MicroRNA paraffin-based studies in osteosarcoma reveal reproducible independent prognostic profiles at 14q32

Background: Although microRNAs (miRNAs) are implicated in osteosarcoma biology and chemoresponse, miRNA prognostic models are still needed, particularly because prognosis is imperfectly correlated with chemoresponse. Formalin-fixed, paraffin-embedded tissue is a necessary resource for biomarker studies in this malignancy with limited frozen tissue availability. Methods: We performed miRNA and mRNA microarray formalin-fixed, paraffin-embedded assays in 65 osteosarcoma biopsy and 26 paired post-chemotherapy resection specimens and used the only publicly available miRNA dataset, generated independently by another group, to externally validate our strongest findings (n = 29). We used supervised principal components analysis and logistic regression for survival and chemoresponse, and miRNA activity and target gene set analysis to study miRNA regulatory activity. Results: Several miRNA-based models with as few as five miRNAs were prognostic independently of pathologically assessed chemoresponse (median recurrence-free survival: 59 months versus not-yet-reached; adjusted hazards ratio = 2.90; P = 0.036). The independent dataset supported the reproducibility of recurrence and survival findings. The prognostic value of the profile was independent of confounding by known prognostic variables, including chemoresponse, tumor location and metastasis at diagnosis. Model performance improved when chemoresponse was added as a covariate (median recurrence-free survival: 59 months versus not-yet-reached; hazard ratio = 3.91; P = 0.002). Most prognostic miRNAs were located at 14q32 - a locus already linked to osteosarcoma - and their gene targets display deregulation patterns associated with outcome. We also identified miRNA profiles predictive of chemoresponse (75% to 80% accuracy), which did not overlap with prognostic profiles. Conclusions: Formalin-fixed, paraffin-embedded tissue-derived miRNA patterns are a powerful prognostic tool for risk-stratified osteosarcoma management strategies. Combined miRNA and mRNA analysis supports a possible role of the 14q32 locus in osteosarcoma progression and outcome. Our study creates a paradigm for formalin-fixed, paraffin-embedded-based miRNA biomarker studies in cancer

Genome-to-genome analysis highlights the effect of the human innate and adaptive immune systems on the hepatitis C virus

Outcomes of hepatitis C virus (HCV) infection and treatment depend on viral and host genetic factors. Here we use human genome-wide genotyping arrays and new whole-genome HCV viral sequencing technologies to perform a systematic genome-to-genome study of 542 individuals who were chronically infected with HCV, predominantly genotype 3. We show that both alleles of genes encoding human leukocyte antigen molecules and genes encoding components of the interferon lambda innate immune system drive viral polymorphism. Additionally, we show that IFNL4 genotypes determine HCV viral load through a mechanism dependent on a specific amino acid residue in the HCV NS5A protein. These findings highlight the interplay between the innate immune system and the viral genome in HCV control

Viral genome wide association study identifies novel hepatitis C virus polymorphisms associated with sofosbuvir treatment failure

Persistent hepatitis C virus (HCV) infection is a major cause of chronic liver disease, worldwide. With the development of direct-acting antivirals, treatment of chronically infected patients has become highly effective, although a subset of patients responds less well to therapy. Sofosbuvir is a common component of current de novo or salvage combination therapies, that targets the HCV NS5B polymerase. We use pre-treatment whole-genome sequences of HCV from 507 patients infected with HCV subtype 3a and treated with sofosbuvir containing regimens to detect viral polymorphisms associated with response to treatment. We find three common polymorphisms in non-targeted HCV NS2 and NS3 proteins are associated with reduced treatment response. These polymorphisms are enriched in post-treatment HCV sequences of patients unresponsive to treatment. They are also associated with lower reductions in viral load in the first week of therapy. Using in vitro short-term dose-response assays, these polymorphisms do not cause any reduction in sofosbuvir potency, suggesting an indirect mechanism of action in decreasing sofosbuvir efficacy. The identification of polymorphisms in NS2 and NS3 proteins associated with poor treatment outcomes emphasises the value of systematic genome-wide analyses of viruses in uncovering clinically relevant polymorphisms that impact treatment

Therapeutic Implications of GIPC1 Silencing in Cancer

GIPC1 is a cytoplasmic scaffold protein that interacts with numerous receptor signaling complexes, and emerging evidence suggests that it plays a role in tumorigenesis. GIPC1 is highly expressed in a number of human malignancies, including breast, ovarian, gastric, and pancreatic cancers. Suppression of GIPC1 in human pancreatic cancer cells inhibits in vivo tumor growth in immunodeficient mice. To better understand GIPC1 function, we suppressed its expression in human breast and colorectal cancer cell lines and human mammary epithelial cells (HMECs) and assayed both gene expression and cellular phenotype. Suppression of GIPC1 promotes apoptosis in MCF-7, MDA-MD231, SKBR-3, SW480, and SW620 cells and impairs anchorage-independent colony formation of HMECs. These observations indicate GIPC1 plays an essential role in oncogenic transformation, and its expression is necessary for the survival of human breast and colorectal cancer cells. Additionally, a GIPC1 knock-down gene signature was used to interrogate publically available breast and ovarian cancer microarray datasets. This GIPC1 signature statistically correlates with a number of breast and ovarian cancer phenotypes and clinical outcomes, including patient survival. Taken together, these data indicate that GIPC1 inhibition may represent a new target for therapeutic development for the treatment of human cancers

Avoiding dynastic, assortative mating, and population stratification biases in Mendelian randomization through within-family analyses.

Estimates from Mendelian randomization studies of unrelated individuals can be biased due to uncontrolled confounding from familial effects. Here we describe methods for within-family Mendelian randomization analyses and use simulation studies to show that family-based analyses can reduce such biases. We illustrate empirically how familial effects can affect estimates using data from 61,008 siblings from the Nord-Trøndelag Health Study and UK Biobank and replicated our findings using 222,368 siblings from 23andMe. Both Mendelian randomization estimates using unrelated individuals and within family methods reproduced established effects of lower BMI reducing risk of diabetes and high blood pressure. However, while Mendelian randomization estimates from samples of unrelated individuals suggested that taller height and lower BMI increase educational attainment, these effects were strongly attenuated in within-family Mendelian randomization analyses. Our findings indicate the necessity of controlling for population structure and familial effects in Mendelian randomization studies

Oral abstracts 1: SpondyloarthropathiesO1. Detecting axial spondyloarthritis amongst primary care back pain referrals

Background: Inflammatory back pain (IBP) is an early feature of ankylosing spondylitis (AS) and its detection offers the prospect of early diagnosis of AS. However, since back pain is very common but only a very small minority of back pain sufferers have ASpA or AS, screening of back pain sufferers for AS is problematic. In early disease radiographs are often normal so that fulfilment of diagnostic criteria for AS is impossible though a diagnosis of axial SpA can be made if MRI evidence of sacroiliitis is present. This pilot study was designed to indicate whether a cost-effective pick up rate for ASpA/early AS could be achieved by identifying adults with IBP stratified on the basis of age. Methods: Patients aged between 18 and 45 years who were referred to a hospital physiotherapy service with back pain of more than 3 months duration were assessed for IBP. All were asked to complete a questionnaire based on the Berlin IBP criteria. Those who fulfilled IBP criteria were also asked to complete a second short questionnaire enquiring about SpA comorbidities, to have a blood test for HLA-B27 and CRP level and to undergo an MRI scan of the sacroiliac joints. This was a limited scan, using STIR, diffusion-weighted, T1 and T2 sequences of the sacroiliac joints to minimize time in the scanner and cost. The study was funded by a research grant from Abbott Laboratories Ltd. Results: 50 sequential patients agreed to participate in the study and completed the IBP questionnaire. Of these 27 (54%) fulfilled criteria for IBP. Of these, 2 patients reported a history of an SpA comorbidity - 1 psoriasis; 1 ulcerative colitis - and 3 reported a family history of an SpA comorbidity - 2 psoriasis; 1 Crohn's disease. 4 were HLA-B27 positive, though results were not available for 7. Two patients had marginally raised CRP levels (6, 10 -NR ≤ 5). 19 agreed to undergo MRI scanning of the sacroiliac joints and lumbar spine; 4 scans were abnormal, showing evidence of bilateral sacroiliitis on STIR sequences. In all cases the changes met ASAS criteria but were limited. Of these 4 patients 3 were HLA-B27 positive but none gave a personal or family history of an SpA-associated comorbidity and all had normal CRP levels. Conclusions: This was a pilot study yielding only limited conclusions. However, it is clear that: Screening of patients referred for physiotherapy for IBP is straightforward, inexpensive and quick. It appears that IBP is more prevalent in young adults than overall population data suggest so that targeting this population may be efficient. IBP questionnaires could be administered routinely during a physiotherapy assessment. HLA-B27 testing in this group of patients with IBP is a suitable screening tool. The sacroiliac joint changes identified were mild and their prognostic significance is not yet clear so that the value of early screening needs further evaluation. Disclosure statement: C.H. received research funding for this study from Abbott. A.K. received research funding for this study, and speaker and consultancy fees, from Abbott. All other authors have declared no conflicts of interes

Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity.

Global dispersal and increasing frequency of the SARS-CoV-2 spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large dataset, well represented by both spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant