ITGB2 (Integrin ß2) Immunomodulatory Gene Variants in Premature Infants With Necrotizing Enterocolitis.

ABSTRACT: Aberrant toll-like receptor (TLR) activation is central to necrotizing enterocolitis (NEC) pathogenesis. ß2 integrins regulate TLR signaling, and integrin ß2 (ITGB2) deficiency causes TLR hyperresponsiveness. To test the hypothesis that ITGB2 genetic variants modulate NEC susceptibility, we sequenced the exonic ITGB2 locus to compare the prevalence of deleterious variants among 221 preterm infants with and without NEC. ITGB2 variants were not associated with NEC in our entire cohort (NEC [9/56] versus controls [16/165], P = 0.19) or in extremely low birthweight infants (ELBW, controls [7.9%] versus NEC [18.2%]; P = 0.11) but were increased compared to the populace (4.5%, gnomad.broadinstitute.org). Combined annotation-dependent depletion -predicted deleterious ITGB2 variants increased proportionately with increasing NEC severity in ELBW infants (controls [6.7%] versus medical NEC [16.7%] versus surgical NEC [19%] (P = 0.03). Although ITGB2 variants were not associated with NEC in our preterm cohort, subgroup analysis showed a trend towards enrichment with NEC severity in ELBW infants.

Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

Preterm birth is a leading cause of morbidity and mortality in infants. Genetic and environmental factors play a role in the susceptibility to preterm birth, but despite many investigations, the genetic basis for preterm birth remain largely unknown. Our objective was to identify rare, possibly damaging, nucleotide variants in mothers from families with recurrent spontaneous preterm births (SPTB). DNA samples from 17 Finnish mothers who delivered at least one infant preterm were subjected to whole exome sequencing. All mothers were of northern Finnish origin and were from seven multiplex families. Additional replication samples of European origin consisted of 93 Danish sister pairs (and two sister triads), all with a history of a preterm delivery. Rare exonic variants (frequency <1%) were analyzed to identify genes and pathways likely to affect SPTB susceptibility. We identified rare, possibly damaging, variants in genes that were common to multiple affected individuals. The glucocorticoid receptor signaling pathway was the most significant (p<1.7e-8) with genes containing these variants in a subgroup of ten Finnish mothers, each having had 2-4 SPTBs. This pathway was replicated among the Danish sister pairs. A gene in this pathway, heat shock protein family A (Hsp70) member 1 like (HSPA1L), contains two likely damaging missense alleles that were found in four different Finnish families. One of the variants (rs34620296) had a higher frequency in cases compared to controls (0.0025 vs. 0.0010, p = 0.002) in a large preterm birth genome-wide association study (GWAS) consisting of mothers of general European ancestry. Sister pairs in replication samples also shared rare, likely damaging HSPA1L variants. Furthermore, in silico analysis predicted an additional phosphorylation site generated by rs34620296 that could potentially affect chaperone activity or HSPA1L protein stability. Finally, in vitro functional experiment showed a link between HSPA1L activity and decidualization. In conclusion, rare, likely damaging, variants in HSPA1L were observed in multiple families with recurrent SPTB.

Correction: Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

[This corrects the article DOI: 10.1371/journal.pgen.1007394.].

Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis.

Monozygotic or 'identical' twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4(+) lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among approximately 3.6 million single nucleotide polymorphisms (SNPs) or approximately 0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of approximately 19,000 genes in CD4(+) T cells. Only 2 to 176 differences in the methylation of approximately 2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to approximately 800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

Alström Syndrome: Mutation Spectrum of ALMS1.

Alström Syndrome (ALMS), a recessive, monogenic ciliopathy caused by mutations in ALMS1, is typically characterized by multisystem involvement including early cone-rod retinal dystrophy and blindness, hearing loss, childhood obesity, type 2 diabetes mellitus, cardiomyopathy, fibrosis, and multiple organ failure. The precise function of ALMS1 remains elusive, but roles in endosomal and ciliary transport and cell cycle regulation have been shown. The aim of our study was to further define the spectrum of ALMS1 mutations in patients with clinical features of ALMS. Mutational analysis in a world-wide cohort of 204 families identified 109 novel mutations, extending the number of known ALMS1 mutations to 239 and highlighting the allelic heterogeneity of this disorder. This study represents the most comprehensive mutation analysis in patients with ALMS, identifying the largest number of novel mutations in a single study worldwide. Here, we also provide an overview of all ALMS1 mutations identified to date.

Renal systems biology of patients with systemic inflammatory response syndrome.

A systems biology approach was used to comprehensively examine the impact of renal disease and hemodialysis (HD) on patient response during critical illness. To achieve this, we examined the metabolome, proteome, and transcriptome of 150 patients with critical illness, stratified by renal function. Quantification of plasma metabolites indicated greater change as renal function declined, with the greatest derangements in patients receiving chronic HD. Specifically, 6 uremic retention molecules, 17 other protein catabolites, 7 modified nucleosides, and 7 pentose phosphate sugars increased as renal function declined, consistent with decreased excretion or increased catabolism of amino acids and ribonucleotides. Similarly, the proteome showed increased levels of low-molecular-weight proteins and acute-phase reactants. The transcriptome revealed a broad-based decrease in mRNA levels among patients on HD. Systems integration revealed an unrecognized association between plasma RNASE1 and several RNA catabolites and modified nucleosides. Further, allantoin, N1-methyl-4-pyridone-3-carboxamide, and N-acetylaspartate were inversely correlated with the majority of significantly downregulated genes. Thus, renal function broadly affected the plasma metabolome, proteome, and peripheral blood transcriptome during critical illness; changes were not effectively mitigated by hemodialysis. These studies allude to several novel mechanisms whereby renal dysfunction contributes to critical illness.

A highly annotated whole-genome sequence of a Korean individual.

Recent advances in sequencing technologies have initiated an era of personal genome sequences. To date, human genome sequences have been reported for individuals with ancestry in three distinct geographical regions: a Yoruba African, two individuals of northwest European origin, and a person from China. Here we provide a highly annotated, whole-genome sequence for a Korean individual, known as AK1. The genome of AK1 was determined by an exacting, combined approach that included whole-genome shotgun sequencing (27.8x coverage), targeted bacterial artificial chromosome sequencing, and high-resolution comparative genomic hybridization using custom microarrays featuring more than 24 million probes. Alignment to the NCBI reference, a composite of several ethnic clades, disclosed nearly 3.45 million single nucleotide polymorphisms (SNPs), including 10,162 non-synonymous SNPs, and 170,202 deletion or insertion polymorphisms (indels). SNP and indel densities were strongly correlated genome-wide. Applying very conservative criteria yielded highly reliable copy number variants for clinical considerations. Potential medical phenotypes were annotated for non-synonymous SNPs, coding domain indels, and structural variants. The integration of several human whole-genome sequences derived from several ethnic groups will assist in understanding genetic ancestry, migration patterns and population bottlenecks.

Molecular diagnosis of infantile onset inflammatory bowel disease by exome sequencing.

Pediatric-onset inflammatory bowel disease (IBD) is known to be associated with severe disease, poor response to therapy, and increased morbidity and mortality. We conducted exome sequencing of two brothers from a non-consanguineous relationship who presented before the age of one with severe infantile-onset IBD, failure to thrive, skin rash, and perirectal abscesses refractory to medical management. We examined the variants discovered in all known IBD-associated and primary immunodeficiency genes in both siblings. The siblings were identified to harbor compound heterozygous mutations in IL10RA (c.784C>T, p.Arg262Cys; c.349C>T, p.Arg117Cys). Upon molecular diagnosis, the proband underwent successful hematopoietic stem cell transplantation and demonstrated marked clinical improvement of all IBD-associated clinical symptoms. Exome sequencing can be an effective tool to aid in the molecular diagnosis of pediatric-onset IBD. We provide additional evidence of the safety and benefit of HSCT for patients with IBD due to mutations in the IL10RA gene.

Diagnosis of mitochondrial disorders by concomitant next-generation sequencing of the exome and mitochondrial genome.

Mitochondrial diseases are notoriously difficult to diagnose due to extreme locus and allelic heterogeneity, with both nuclear and mitochondrial genomes potentially liable. Using exome sequencing we demonstrate the ability to rapidly and cost effectively evaluate both the nuclear and mitochondrial genomes to obtain a molecular diagnosis for four patients with three distinct mitochondrial disorders. One patient was found to have Leigh syndrome due to a mutation in MT-ATP6, two affected siblings were discovered to be compound heterozygous for mutations in the NDUFV1 gene, which causes mitochondrial complex I deficiency, and one patient was found to have coenzyme Q10 deficiency due to compound heterozygous mutations in COQ2. In all cases conventional diagnostic testing failed to identify a molecular diagnosis. We suggest that additional studies should be conducted to evaluate exome sequencing as a primary diagnostic test for mitochondrial diseases, including those due to mtDNA mutations.

Modeling.

Modeling is coming to the fore as it is now widely accepted and indeed expected during drug discovery and development. Modeling integrates knowledge, increases understanding and provides the ability to predict an outcome either before it occurs or when it is not possible to measure. This makes modeling an attractive option for inhaled drugs as it is not possible to routinely measure what is occurring to the drug (pharmacokinetics) and what effect the drug is having (pharmacodynamics) at local microscopic sites of such a diverse and complex organ as the lung. Many pieces of information (data and knowledge) exist like the pieces of a jigsaw puzzle and modeling brings the pieces together in a scientific and mechanistically coherent manner to increase understanding of both the efficacy and safety of inhaled drugs.

Genome sequencing and mapping reveal loss of heterozygosity as a mechanism for rapid adaptation in the vegetable pathogen Phytophthora capsici.

The oomycete vegetable pathogen Phytophthora capsici has shown remarkable adaptation to fungicides and new hosts. Like other members of this destructive genus, P. capsici has an explosive epidemiology, rapidly producing massive numbers of asexual spores on infected hosts. In addition, P. capsici can remain dormant for years as sexually recombined oospores, making it difficult to produce crops at infested sites, and allowing outcrossing populations to maintain significant genetic variation. Genome sequencing, development of a high-density genetic map, and integrative genomic or genetic characterization of P. capsici field isolates and intercross progeny revealed significant mitotic loss of heterozygosity (LOH) in diverse isolates. LOH was detected in clonally propagated field isolates and sexual progeny, cumulatively affecting >30% of the genome. LOH altered genotypes for more than 11,000 single-nucleotide variant sites and showed a strong association with changes in mating type and pathogenicity. Overall, it appears that LOH may provide a rapid mechanism for fixing alleles and may be an important component of adaptability for P. capsici.

Physiologically Based Pharmacokinetic Modelling of Inhaled Nemiralisib: Mechanistic Components for Pulmonary Absorption, Systemic Distribution, and Oral Absorption.

BACKGROUND AND OBJECTIVES: Physiologically based pharmacokinetic (PBPK) modelling has evolved to accommodate different routes of drug administration and enables prediction of drug concentrations in tissues as well as plasma. The inhalation route of administration has proven successful in treating respiratory diseases but can also be used for rapid systemic delivery, holding great promise for treatment of diseases requiring systemic exposure. The objective of this work was to develop a PBPK model that predicts plasma and tissue concentrations following inhalation administration of the PI3Kδ inhibitor nemiralisib. METHODS: A PBPK model was built in GastroPlus® that includes a complete mechanistic description of pulmonary absorption, systemic distribution and oral absorption following inhalation administration of nemiralisib. The availability of clinical data obtained after intravenous, oral and inhalation administration enabled validation of the model with observed data and accurate assessment of pulmonary drug absorption. The PBPK model described in this study incorporates novel use of key parameters such as lung systemic absorption rate constants derived from human physiological lung blood flows, and implementation of the specific permeability-surface area product per millilitre of tissue cell volume (SpecPStc) to predict tissue distribution. RESULTS: The inhaled PBPK model was verified using plasma and bronchoalveolar lavage fluid concentration data obtained in human subjects. Prediction of tissue concentrations using the permeability-limited systemic disposition tissue model was further validated using tissue concentration data obtained in the rat following intravenous infusion administration to steady state. CONCLUSIONS: Fully mechanistic inhaled PBPK models such as the model described herein could be applied for cross molecule assessments with respect to lung retention and systemic exposure, both in terms of pharmacology and toxicology, and may facilitate clinical indication selection.

CYP2C8, CYP2C9, and CYP2C19 Characterization Using Next-Generation Sequencing and Haplotype Analysis: A GeT-RM Collaborative Project.

Pharmacogenetic tests typically target selected sequence variants to identify haplotypes that are often defined by star (∗) allele nomenclature. Due to their design, these targeted genotyping assays are unable to detect novel variants that may change the function of the gene product and thereby affect phenotype prediction and patient care. In the current study, 137 DNA samples that were previously characterized by the Genetic Testing Reference Material (GeT-RM) program using a variety of targeted genotyping methods were recharacterized using targeted and whole genome sequencing analysis. Sequence data were analyzed using three genotype calling tools to identify star allele diplotypes for CYP2C8, CYP2C9, and CYP2C19. The genotype calls from next-generation sequencing (NGS) correlated well to those previously reported, except when novel alleles were present in a sample. Six novel alleles and 38 novel suballeles were identified in the three genes due to identification of variants not covered by targeted genotyping assays. In addition, several ambiguous genotype calls from a previous study were resolved using the NGS and/or long-read NGS data. Diplotype calls were mostly consistent between the calling algorithms, although several discrepancies were noted. This study highlights the utility of NGS for pharmacogenetic testing and demonstrates that there are many novel alleles that are yet to be discovered, even in highly characterized genes such as CYP2C9 and CYP2C19.

Transcriptome sequencing of malignant pleural mesothelioma tumors.

Cancers arise by the gradual accumulation of mutations in multiple genes. We now use shotgun pyrosequencing to characterize RNA mutations and expression levels unique to malignant pleural mesotheliomas (MPMs) and not present in control tissues. On average, 266 Mb of cDNA were sequenced from each of four MPMs, from a control pulmonary adenocarcinoma (ADCA), and from normal lung tissue. Previously observed differences in MPM RNA expression levels were confirmed. Point mutations were identified by using criteria that require the presence of the mutation in at least four reads and in both cDNA strands and the absence of the mutation from sequence databases, normal adjacent tissues, and other controls. In the four MPMs, 15 nonsynonymous mutations were discovered: 7 were point mutations, 3 were deletions, 4 were exclusively expressed as a consequence of imputed epigenetic silencing, and 1 was putatively expressed as a consequence of RNA editing. Notably, each MPM had a different mutation profile, and no mutated gene was previously implicated in MPM. Of the seven point mutations, three were observed in at least one tumor from 49 other MPM patients. The mutations were in genes that could be causally related to cancer and included XRCC6, PDZK1IP1, ACTR1A, and AVEN.

PharmVar GeneFocus: CYP2B6.

The Pharmacogene Variation Consortium (PharmVar) catalogs star (*) allele nomenclature for the polymorphic human CYP2B6 gene. Genetic variation within the CYP2B6 gene locus impacts the metabolism or bioactivation of clinically important drugs. Of particular importance are efficacy and safety concerns regarding: efavirenz, which is used for the treatment of HIV type-1 infection; methadone, a mainstay in the treatment of opioid use disorder and as an analgesic; ketamine, used as an antidepressant and analgesic; and bupropion, which is prescribed to treat depression and for smoking cessation. This GeneFocus provides a comprehensive overview and summary of CYP2B6 and describes how haplotype information catalogued by PharmVar is utilized by the Pharmacogenomics Knowledgebase (PharmGKB) and the Clinical Pharmacogenetics Implementation Consortium (CPIC).

Long-Distance Phasing of a Tentative "Enhancer" Single-Nucleotide Polymorphism With CYP2D6 Star Allele Definitions.

BACKGROUND: The CYP2D6 gene locus has been extensively studied over decades, yet a portion of variability in CYP2D6 activity cannot be explained by known sequence variations within the gene, copy number variation, or structural rearrangements. It was proposed that rs5758550, located 116 kb downstream of the CYP2D6 gene locus, increases gene expression and thus contributes to variability in CYP2D6 activity. This finding has, however, not been validated. The purpose of the study was to address a major technological barrier, i.e., experimentally linking rs5758550, also referred to as the "enhancer" single-nucleotide polymorphism (SNP), to CYP2D6 haplotypes >100 kb away. To overcome this challenge is essential to ultimately determine the contribution of the "enhancer" SNP to interindividual variability in CYP2D6 activity. METHODS: A large ethnically mixed population sample (n=3,162) was computationally phased to determine linkage between the "enhancer" SNP and CYP2D6 haplotypes (or star alleles). To experimentally validate predicted linkages, DropPhase2D6, a digital droplet PCR (ddPCR)-based method was developed. 10X Genomics Linked-Reads were utilized as a proof of concept. RESULTS: Phasing predicted that the "enhancer" SNP can occur on numerous CYP2D6 haplotypes including CYP2D6*1, *2, *5, and *41 and suggested that linkage is incomplete, i.e., a portion of these alleles do not have the "enhancer" SNP. Phasing also revealed differences among the European and African ancestry data sets regarding the proportion of alleles with and without the "enhancer" SNP. DropPhase2D6 was utilized to confirm or refute the predicted "enhancer" SNP location for individual samples, e.g., of n=3 samples genotyped as *1/*41, rs5758550 was on the *41 allele of two samples and on the *1 allele of one sample. Our findings highlight that the location of the "enhancer" SNP must not be assigned by "default." Furthermore, linkage between the "enhancer" SNP and CYP2D6 star allele haplotypes was confirmed with 10X Genomics technology. CONCLUSIONS: Since the "enhancer" SNP can be present on a portion of normal, decreased, or no function alleles, the phase of the "enhancer" SNP must be considered when investigating the impact of the "enhancer" SNP on CYP2D6 activity.

Physiologically Based Pharmacokinetic Modelling for First-In-Human Predictions: An Updated Model Building Strategy Illustrated with Challenging Industry Case Studies.

Physiologically based pharmacokinetic modelling is well established in the pharmaceutical industry and is accepted by regulatory agencies for the prediction of drug-drug interactions. However, physiologically based pharmacokinetic modelling is valuable to address a much wider range of pharmaceutical applications, and new regulatory impact is expected as its full power is leveraged. As one example, physiologically based pharmacokinetic modelling is already routinely used during drug discovery for in-vitro to in-vivo translation and pharmacokinetic modelling in preclinical species, and this leads to the application of verified models for first-in-human pharmacokinetic predictions. A consistent cross-industry strategy in this application area would increase confidence in the approach and facilitate further learning. With this in mind, this article aims to enhance a previously published first-in-human physiologically based pharmacokinetic model-building strategy. Based on the experience of scientists from multiple companies participating in the GastroPlus™ User Group Steering Committee, new Absorption, Distribution, Metabolism and Excretion knowledge is integrated and decision trees proposed for each essential component of a first-in-human prediction. We have reviewed many relevant scientific publications to identify new findings and highlight gaps that need to be addressed. Finally, four industry case studies for more challenging compounds illustrate and highlight key components of the strategy.

The Pharmacogene Variation (PharmVar) Consortium: Incorporation of the Human Cytochrome P450 (CYP) Allele Nomenclature Database.

The Human Cytochrome P450 (CYP) Allele Nomenclature Database, a critical resource for the pharmacogenetics and genomics communities, has transitioned to the Pharmacogene Variation (PharmVar) Consortium. In this report we provide a summary of the current database, provide an overview of the PharmVar consortium, and highlight the PharmVar database which will serve as the new home for pharmacogene nomenclature.