Multiple Data Analyses and Statistical Approaches for Analyzing Data from Metagenomic Studies and Clinical Trials

Metagenomics, also known as environmental genomics, is the study of the genomic content of a sample of organisms (microbes) obtained from a common habitat. Metagenomics and other “omics” disciplines have captured the attention of researchers for several decades. The effect of microbes in our body is a relevant concern for health studies. There are plenty of studies using metagenomics which examine microorganisms that inhabit niches in the human body, sometimes causing disease, and are often correlated with multiple treatment conditions. No matter from which environment it comes, the analyses are often aimed at determining either the presence or absence of specific species of interest in a given metagenome or comparing the biological diversity and the functional activity of a wider range of microorganisms within their communities. The importance increases for comparison within different environments such as multiple patients with different conditions, multiple drugs, and multiple time points of same treatment or same patient. Thus, no matter how many hypotheses we have, we need a good understanding of genomics, bioinformatics, and statistics to work together to analyze and interpret these datasets in a meaningful way. This chapter provides an overview of different data analyses and statistical approaches (with example scenarios) to analyze metagenomics samples from different medical projects or clinical trials

De Novo Truncating Mutations in WASF1 Cause Intellectual Disability with Seizures.

Next-generation sequencing has been invaluable in the elucidation of the genetic etiology of many subtypes of intellectual disability in recent years. Here, using exome sequencing and whole-genome sequencing, we identified three de novo truncating mutations in WAS protein family member 1 (WASF1) in five unrelated individuals with moderate to profound intellectual disability with autistic features and seizures. WASF1, also known as WAVE1, is part of the WAVE complex and acts as a mediator between Rac-GTPase and actin to induce actin polymerization. The three mutations connected by Matchmaker Exchange were c.1516C>T (p.Arg506Ter), which occurs in three unrelated individuals, c.1558C>T (p.Gln520Ter), and c.1482delinsGCCAGG (p.Ile494MetfsTer23). All three variants are predicted to partially or fully disrupt the C-terminal actin-binding WCA domain. Functional studies using fibroblast cells from two affected individuals with the c.1516C>T mutation showed a truncated WASF1 and a defect in actin remodeling. This study provides evidence that de novo heterozygous mutations in WASF1 cause a rare form of intellectual disability

Long-Term Outcomes with Subcutaneous C1-Inhibitor Replacement Therapy for Prevention of Hereditary Angioedema Attacks

Background For the prevention of attacks of hereditary angioedema (HAE), the efficacy and safety of subcutaneous human C1-esterase inhibitor (C1-INH[SC]; HAEGARDA, CSL Behring) was established in the 16-week Clinical Study for Optimal Management of Preventing Angioedema with Low-Volume Subcutaneous C1-Inhibitor Replacement Therapy (COMPACT). Objective To assess the long-term safety, occurrence of angioedema attacks, and use of rescue medication with C1-INH(SC). Methods Open-label, randomized, parallel-arm extension of COMPACT across 11 countries. Patients with frequent angioedema attacks, either study treatment-naive or who had completed COMPACT, were randomly assigned (1:1) to 40 IU/kg or 60 IU/kg C1-INH(SC) twice per week, with conditional uptitration to optimize prophylaxis (ClinicalTrials.gov registration no. NCT02316353). Results A total of 126 patients with a monthly attack rate of 4.3 in 3 months before entry in COMPACT were enrolled and treated for a mean of 1.5 years; 44 patients (34.9%) had more than 2 years of exposure. Mean steady-state C1-INH functional activity increased to 66.6% with 60 IU/kg. Incidence of adverse events was low and similar in both dose groups (11.3 and 8.5 events per patient-year for 40 IU/kg and 60 IU/kg, respectively). For 40 IU/kg and 60 IU/kg, median annualized attack rates were 1.3 and 1.0, respectively, and median rescue medication use was 0.2 and 0.0 times per year, respectively. Of 23 patients receiving 60 IU/kg for more than 2 years, 19 (83%) were attack-free during months 25 to 30 of treatment. Conclusions In patients with frequent HAE attacks, long-term replacement therapy with C1-INH(SC) is safe and exhibits a substantial and sustained prophylactic effect, with the vast majority of patients becoming free from debilitating disease symptoms

Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

Telomere length is a risk factor in disease and the dynamics of telomere length are crucial to our understanding of cell replication and vitality. The proliferation of whole genome sequencing represents an unprecedented opportunity to glean new insights into telomere biology on a previously unimaginable scale. To this end, a number of approaches for estimating telomere length from whole-genome sequencing data have been proposed. Here we present Telomerecat, a novel approach to the estimation of telomere length. Previous methods have been dependent on the number of telomeres present in a cell being known, which may be problematic when analysing aneuploid cancer data and non-human samples. Telomerecat is designed to be agnostic to the number of telomeres present, making it suited for the purpose of estimating telomere length in cancer studies. Telomerecat also accounts for interstitial telomeric reads and presents a novel approach to dealing with sequencing errors. We show that Telomerecat performs well at telomere length estimation when compared to leading experimental and computational methods. Furthermore, we show that it detects expected patterns in longitudinal data, repeated measurements, and cross-species comparisons. We also apply the method to a cancer cell data, uncovering an interesting relationship with the underlying telomerase genotype

Loss of function NFKB1 variants are the most common monogenic cause of CVID in Europeans.

BACKGROUND: The genetic etiology of primary immunodeficiency disease (PID) carries prognostic information. OBJECTIVE: We conducted a whole-genome sequencing study assessing a large proportion of the NIHR-BioResource - Rare Disease cohort. METHODS: In the predominantly European study population of principally sporadic unrelated PID cases (n=846), a novel Bayesian method identified NFKB1 as one most strongly associated with PID, and the association was explained by 16 novel heterozygous truncating, missense and gene deletion variants. This accounted for 4% of common variable immunodeficiency (CVID) cases (n=390) in the cohort. Amino-acid substitutions predicted to be pathogenic were assessed by analysis of structural protein data. Immunophenotyping, immunoblotting and ex vivo stimulation of lymphocytes determined the functional effects of these variants. Detailed clinical and pedigree information was collected for genotype-phenotype co-segregation analyses. RESULTS: Both sporadic and familial cases demonstrated evidence of the non-infective complications of CVID, including massive lymphadenopathy (24%), unexplained splenomegaly (48%) and autoimmune disease (48%), features prior studies correlate with worse clinical prognosis. Although partial penetrance of clinical symptoms was noted in certain pedigrees, all carriers have a deficiency in B lymphocyte differentiation. Detailed assessment of B lymphocyte numbers, phenotype and function identifies the presence of a raised CD21lowB cell population: combined with identification of the disease-causing variant, this distinguishes between healthy individuals, asymptomatic carriers and clinically affected cases. CONCLUSION: We show that heterozygous loss-of-function variants in NFKB1 are the most common known monogenic cause of CVID that results in a temporally progressive defect in the formation of immunoglobulin-producing B cells.This study was supported by The National Institute for Health Research England (grant number RG65966), and by the Center of Immunodeficiencies Amsterdam (CIDA). JET is supported by an MRC Clinician Scientist Fellowship (MR/L006197/1). AJT is supported by both the Wellcome Trust (104807/Z/14/Z) and by the National Institute for Health Research Biomedical Research Centre at Great Ormond Street Hospital for Children NHS Foundation Trust and University College London. EO receives personal fees from CSL Behring and MSD

Biallelic Mutation of ARHGEF18, Involved in the Determination of Epithelial Apicobasal Polarity, Causes Adult-Onset Retinal Degeneration

Mutations in more than 250 genes are implicated in inherited retinal dystrophy; the encoded proteins are involved in a broad spectrum of pathways. The presence of unsolved families after highly parallel sequencing strategies suggests that further genes remain to be identified. Whole-exome and -genome sequencing studies employed here in large cohorts of affected individuals revealed biallelic mutations in ARHGEF18 in three such individuals. ARHGEF18 encodes ARHGEF18, a guanine nucleotide exchange factor that activates RHOA, a small GTPase protein that is a key component of tight junctions and adherens junctions. This biological pathway is known to be important for retinal development and function, as mutation of CRB1, encoding another component, causes retinal dystrophy. The retinal structure in individuals with ARHGEF18 mutations resembled that seen in subjects with CRB1 mutations. Five mutations were found on six alleles in the three individuals: c.808A>G (p.Thr270Ala), c.1617+5G>A (p.Asp540Glyfs∗63), c.1996C>T (p.Arg666∗), c.2632G>T (p.Glu878∗), and c.2738_2761del (p.Arg913_Glu920del). Functional tests suggest that each disease genotype might retain some ARHGEF18 activity, such that the phenotype described here is not the consequence of nullizygosity. In particular, the p.Thr270Ala missense variant affects a highly conserved residue in the DBL homology domain, which is required for the interaction and activation of RHOA. Previously, knock-out of Arhgef18 in the medaka fish has been shown to cause larval lethality which is preceded by retinal defects that resemble those seen in zebrafish Crumbs complex knock-outs. The findings described here emphasize the peculiar sensitivity of the retina to perturbations of this pathway, which is highlighted as a target for potential therapeutic strategies

Publisher Correction: Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper