Defining predictors of responsiveness to advanced therapies in Crohn's disease and ulcerative colitis: protocol for the IBD-RESPONSE and nested CD-metaRESPONSE prospective, multicentre, observational cohort study in precision medicine.

INTRODUCTION: Characterised by chronic inflammation of the gastrointestinal tract, inflammatory bowel disease (IBD) symptoms including diarrhoea, abdominal pain and fatigue can significantly impact patient's quality of life. Therapeutic developments in the last 20 years have revolutionised treatment. However, clinical trials and real-world data show primary non-response rates up to 40%. A significant challenge is an inability to predict which treatment will benefit individual patients.Current understanding of IBD pathogenesis implicates complex interactions between host genetics and the gut microbiome. Most cohorts studying the gut microbiota to date have been underpowered, examined single treatments and produced heterogeneous results. Lack of cross-treatment comparisons and well-powered independent replication cohorts hampers the ability to infer real-world utility of predictive signatures.IBD-RESPONSE will use multi-omic data to create a predictive tool for treatment response. Future patient benefit may include development of biomarker-based treatment stratification or manipulation of intestinal microbial targets. IBD-RESPONSE and downstream studies have the potential to improve quality of life, reduce patient risk and reduce expenditure on ineffective treatments. METHODS AND ANALYSIS: This prospective, multicentre, observational study will identify and validate a predictive model for response to advanced IBD therapies, incorporating gut microbiome, metabolome, single-cell transcriptome, human genome, dietary and clinical data. 1325 participants commencing advanced therapies will be recruited from ~40 UK sites. Data will be collected at baseline, week 14 and week 54. The primary outcome is week 14 clinical response. Secondary outcomes include clinical remission, loss of response in week 14 responders, corticosteroid-free response/remission, time to treatment escalation and change in patient-reported outcome measures. ETHICS AND DISSEMINATION: Ethical approval was obtained from the Wales Research Ethics Committee 5 (ref: 21/WA/0228). Recruitment is ongoing. Following study completion, results will be submitted for publication in peer-reviewed journals and presented at scientific meetings. Publications will be summarised at www.ibd-response.co.uk. TRIAL REGISTRATION NUMBER: ISRCTN96296121.

Baseline Expression of Immune Gene Modules in Blood is Associated With Primary Response to Anti-TNF Therapy in Crohn's Disease Patients.

BACKGROUND AND AIMS: Anti-tumour necrosis factor [anti-TNF] therapy is widely used for the treatment of inflammatory bowel disease, yet many patients are primary non-responders, failing to respond to induction therapy. We aimed to identify blood gene expression differences between primary responders and primary non-responders to anti-TNF monoclonal antibodies [infliximab and adalimumab], and to predict response status from blood gene expression and clinical data. METHODS: The Personalised Anti-TNF Therapy in Crohn's Disease [PANTS] study is a UK-wide prospective observational cohort study of anti-TNF therapy outcome in anti-TNF-naive Crohn's disease patients [ClinicalTrials.gov identifier: NCT03088449]. Blood gene expression in 324 unique patients was measured by RNA-sequencing at baseline [week 0], and at weeks 14, 30, and 54 after treatment initiation [total sample size = 814]. RESULTS: After adjusting for clinical covariates and estimated blood cell composition, baseline expression of major histocompatibility complex, antigen presentation, myeloid cell enriched receptor, and other innate immune gene modules was significantly higher in anti-TNF responders vs non-responders. Expression changes from baseline to week 14 were generally of consistent direction but greater magnitude [i.e. amplified] in responders, but interferon-related genes were upregulated uniquely in non-responders. Expression differences between responders and non-responders observed at week 14 were maintained at weeks 30 and 54. Prediction of response status from baseline clinical data, cell composition, and module expression was poor. CONCLUSIONS: Baseline gene module expression was associated with primary response to anti-TNF therapy in PANTS patients. However, these baseline expression differences did not predict response with sufficient sensitivity for clinical use.

Effects of psoriasis and psoralen exposure on the somatic mutation landscape of the skin.

Somatic mutations are hypothesized to play a role in many non-neoplastic diseases. We performed whole-exome sequencing of 1,182 microbiopsies dissected from lesional and nonlesional epidermis from 111 patients with psoriasis to search for evidence that somatic mutations in keratinocytes may influence the disease process. Lesional skin remained highly polyclonal, showing no evidence of large-scale spread of clones carrying potentially pathogenic mutations. The mutation rate of keratinocytes was similarly only modestly affected by the disease. We found evidence of positive selection in previously reported driver genes NOTCH1, NOTCH2, TP53, FAT1 and PPM1D and also identified mutations in four genes (GXYLT1, CHEK2, ZFP36L2 and EEF1A1) that we hypothesize are selected for in squamous epithelium irrespective of disease status. Finally, we describe a mutational signature of psoralens-a class of chemicals previously found in some sunscreens and which are used as part of PUVA (psoralens and ultraviolet-A) photochemotherapy treatment for psoriasis.

An Integrated Taxonomy for Monogenic Inflammatory Bowel Disease.

BACKGROUND & AIMS: Monogenic forms of inflammatory bowel disease (IBD) illustrate the essential roles of individual genes in pathways and networks safeguarding immune tolerance and gut homeostasis. METHODS: To build a taxonomy model, we assessed 165 disorders. Genes were prioritized based on penetrance of IBD and disease phenotypes were integrated with multi-omics datasets. Monogenic IBD genes were classified by (1) overlapping syndromic features, (2) response to hematopoietic stem cell transplantation, (3) bulk RNA-sequencing of 32 tissues, (4) single-cell RNA-sequencing of >50 cell subsets from the intestine of healthy individuals and patients with IBD (pediatric and adult), and (5) proteomes of 43 immune subsets. The model was validated by addition of newly identified monogenic IBD defects. As a proof-of-concept, we explore the intersection between immunometabolism and antimicrobial activity for a group of disorders (G6PC3/SLC37A4). RESULTS: Our quantitative integrated taxonomy defines the cellular landscape of monogenic IBD gene expression across 102 genes with high and moderate penetrance (81 in the model set and 21 genes in the validation set). We illustrate distinct cellular networks, highlight expression profiles across understudied cell types (e.g., CD8+ T cells, neutrophils, epithelial subsets, and endothelial cells) and define genotype-phenotype associations (perianal disease and defective antimicrobial activity). We illustrate processes and pathways shared across cellular compartments and phenotypic groups and highlight cellular immunometabolism with mammalian target of rapamycin activation as one of the converging pathways. There is an overlap of genes and enriched cell-specific expression between monogenic and polygenic IBD. CONCLUSION: Our taxonomy integrates genetic, clinical and multi-omic data; providing a basis for genomic diagnostics and testable hypotheses for disease functions and treatment responses.

Somatic mutations provide important and unique insights into the biology of complex diseases.

Somatic evolution of cells within the body is well known to lead to cancers. However, spread of somatic mutations within a tissue over time may also contribute to the pathogenesis of non-neoplastic diseases. Recent years have seen the publication of many studies aiming to characterize somatic evolution in healthy tissues. A logical next step is to extend such work to diseased conditions. As our understanding of the interplay between somatic mutations and non-neoplastic disease grows, opportunities for the joint study of germline and somatic variants will present themselves. Here, we present our thoughts on the utility of somatic mutations for understanding both the causes and consequences of common complex disease and the challenges that remain for the joint study of the soma and germline.

Common and Rare Variant Prediction and Penetrance of IBD in a Large, Multi-ethnic, Health System-based Biobank Cohort.

BACKGROUND AND AIMS: Polygenic risk scores (PRS) may soon be used to predict inflammatory bowel disease (IBD) risk in prevention efforts. We leveraged exome-sequence and single nucleotide polymorphism (SNP) array data from 29,358 individuals in the multiethnic, randomly ascertained health system-based BioMe biobank to define effects of common and rare IBD variants on disease prediction and pathophysiology. METHODS: PRS were calculated from European, African American, and Ashkenazi Jewish (AJ) reference case-control studies, and a meta-GWAS run using all three association datasets. PRS were then combined using regression to assess which combination of scores best predicted IBD status in European, AJ, Hispanic, and African American cohorts in BioMe. Additionally, rare variants were assessed in genes associated with very early-onset IBD (VEO-IBD), by estimating genetic penetrance in each BioMe population. RESULTS: Combining risk scores based on association data from distinct ancestral populations improved IBD prediction for every population in BioMe and significantly improved prediction among European ancestry UK Biobank individuals. Lower predictive power for non-Europeans was observed, reflecting in part substantially lower African IBD case-control reference sizes. We replicated associations for two VEO-IBD genes, ADAM17 and LRBA, with high dominant model penetrance in BioMe. Autosomal recessive LRBA risk alleles are associated with severe, early-onset autoimmunity; we show that heterozygous carriage of an African-predominant LRBA protein-altering allele is associated with significantly decreased LRBA and CTLA-4 expression with T-cell activation. CONCLUSIONS: Greater genetic diversity in African populations improves prediction across populations, and generalizes some VEO-IBD genes. Increasing African American IBD case-collections should be prioritized to reduce health disparities and enhance pathophysiological insight.

Somatic Evolution in Non-neoplastic IBD-Affected Colon.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease associated with increased risk of gastrointestinal cancers. We whole-genome sequenced 446 colonic crypts from 46 IBD patients and compared these to 412 crypts from 41 non-IBD controls from our previous publication on the mutation landscape of the normal colon. The average mutation rate of affected colonic epithelial cells is 2.4-fold that of healthy colon, and this increase is mostly driven by acceleration of mutational processes ubiquitously observed in normal colon. In contrast to the normal colon, where clonal expansions outside the confines of the crypt are rare, we observed widespread millimeter-scale clonal expansions. We discovered non-synonymous mutations in ARID1A, FBXW7, PIGR, ZC3H12A, and genes in the interleukin 17 and Toll-like receptor pathways, under positive selection in IBD. These results suggest distinct selection mechanisms in the colitis-affected colon and that somatic mutations potentially play a causal role in IBD pathogenesis.

Somatic mosaicism and common genetic variation contribute to the risk of very-early-onset inflammatory bowel disease.

Very-early-onset inflammatory bowel disease (VEO-IBD) is a heterogeneous phenotype associated with a spectrum of rare Mendelian disorders. Here, we perform whole-exome-sequencing and genome-wide genotyping in 145 patients (median age-at-diagnosis of 3.5 years), in whom no Mendelian disorders were clinically suspected. In five patients we detect a primary immunodeficiency or enteropathy, with clinical consequences (XIAP, CYBA, SH2D1A, PCSK1). We also present a case study of a VEO-IBD patient with a mosaic de novo, pathogenic allele in CYBB. The mutation is present in ~70% of phagocytes and sufficient to result in defective bacterial handling but not life-threatening infections. Finally, we show that VEO-IBD patients have, on average, higher IBD polygenic risk scores than population controls (99 patients and 18,780 controls; P < 4 × 10-10), and replicate this finding in an independent cohort of VEO-IBD cases and controls (117 patients and 2,603 controls; P < 5 × 10-10). This discovery indicates that a polygenic component operates in VEO-IBD pathogenesis.

HLA-DQA1*05 Carriage Associated With Development of Anti-Drug Antibodies to Infliximab and Adalimumab in Patients With Crohn's Disease.

BACKGROUND & AIMS: Anti-tumor necrosis factor (anti-TNF) therapies are the most widely used biologic drugs for treating immune-mediated diseases, but repeated administration can induce the formation of anti-drug antibodies. The ability to identify patients at increased risk for development of anti-drug antibodies would facilitate selection of therapy and use of preventative strategies. METHODS: We performed a genome-wide association study to identify variants associated with time to development of anti-drug antibodies in a discovery cohort of 1240 biologic-naïve patients with Crohn's disease starting infliximab or adalimumab therapy. Immunogenicity was defined as an anti-drug antibody titer ≥10 AU/mL using a drug-tolerant enzyme-linked immunosorbent assay. Significant association signals were confirmed in a replication cohort of 178 patients with inflammatory bowel disease. RESULTS: The HLA-DQA1*05 allele, carried by approximately 40% of Europeans, significantly increased the rate of immunogenicity (hazard ratio [HR], 1.90; 95% confidence interval [CI], 1.60-2.25; P = 5.88 × 10-13). The highest rates of immunogenicity, 92% at 1 year, were observed in patients treated with infliximab monotherapy who carried HLA-DQA1*05; conversely the lowest rates of immunogenicity, 10% at 1 year, were observed in patients treated with adalimumab combination therapy who did not carry HLA-DQA1*05. We confirmed this finding in the replication cohort (HR, 2.00; 95% CI, 1.35-2.98; P = 6.60 × 10-4). This association was consistent for patients treated with adalimumab (HR, 1.89; 95% CI, 1.32-2.70) or infliximab (HR, 1.92; 95% CI, 1.57-2.33), and for patients treated with anti-TNF therapy alone (HR, 1.75; 95% CI, 1.37-2.22) or in combination with an immunomodulator (HR, 2.01; 95% CI, 1.57-2.58). CONCLUSIONS: In an observational study, we found a genome-wide significant association between HLA-DQA1*05 and the development of antibodies against anti-TNF agents. A randomized controlled biomarker trial is required to determine whether pretreatment testing for HLA-DQA1*05 improves patient outcomes by helping physicians select anti-TNF and combination therapies. ClinicalTrials.gov ID: NCT03088449.

Association of Genetic Variants in NUDT15 With Thiopurine-Induced Myelosuppression in Patients With Inflammatory Bowel Disease.

Importance: Use of thiopurines may be limited by myelosuppression. TPMT pharmacogenetic testing identifies only 25% of at-risk patients of European ancestry. Among patients of East Asian ancestry, NUDT15 variants are associated with thiopurine-induced myelosuppression (TIM). Objective: To identify genetic variants associated with TIM among patients of European ancestry with inflammatory bowel disease (IBD). Design, Setting, and Participants: Case-control study of 491 patients affected by TIM and 679 thiopurine-tolerant unaffected patients who were recruited from 89 international sites between March 2012 and November 2015. Genome-wide association studies (GWAS) and exome-wide association studies (EWAS) were conducted in patients of European ancestry. The replication cohort comprised 73 patients affected by TIM and 840 thiopurine-tolerant unaffected patients. Exposures: Genetic variants associated with TIM. Main Outcomes and Measures: Thiopurine-induced myelosuppression, defined as a decline in absolute white blood cell count to 2.5 × 109/L or less or a decline in absolute neutrophil cell count to 1.0 × 109/L or less leading to a dose reduction or drug withdrawal. Results: Among 1077 patients (398 affected and 679 unaffected; median age at IBD diagnosis, 31.0 years [interquartile range, 21.2 to 44.1 years]; 540 [50%] women; 602 [56%] diagnosed as having Crohn disease), 919 (311 affected and 608 unaffected) were included in the GWAS analysis and 961 (328 affected and 633 unaffected) in the EWAS analysis. The GWAS analysis confirmed association of TPMT (chromosome 6, rs11969064) with TIM (30.5% [95/311] affected vs 16.4% [100/608] unaffected patients; odds ratio [OR], 2.3 [95% CI, 1.7 to 3.1], P = 5.2 × 10-9). The EWAS analysis demonstrated an association with an in-frame deletion in NUDT15 (chromosome 13, rs746071566) and TIM (5.8% [19/328] affected vs 0.2% [1/633] unaffected patients; OR, 38.2 [95% CI, 5.1 to 286.1], P = 1.3 × 10-8), which was replicated in a different cohort (2.7% [2/73] affected vs 0.2% [2/840] unaffected patients; OR, 11.8 [95% CI, 1.6 to 85.0], P = .03). Carriage of any of 3 coding NUDT15 variants was associated with an increased risk (OR, 27.3 [95% CI, 9.3 to 116.7], P = 1.1 × 10-7) of TIM, independent of TPMT genotype and thiopurine dose. Conclusions and Relevance: Among patients of European ancestry with IBD, variants in NUDT15 were associated with increased risk of TIM. These findings suggest that NUDT15 genotyping may be considered prior to initiation of thiopurine therapy; however, further study including additional validation in independent cohorts is required.

Combined Influence of B-Cell Receptor Rearrangement and Somatic Hypermutation on B-Cell Class-Switch Fate in Health and in Chronic Lymphocytic Leukemia.

A diverse B-cell receptor (BCR) repertoire is required to bind a wide range of antigens. BCRs are generated through genetic recombination and can be diversified through somatic hypermutation (SHM) or class-switch recombination (CSR). Patterns of repertoire diversity can vary substantially between different health conditions. We use isotype-resolved BCR sequencing to compare B-cell evolution and class-switch fate in healthy individuals and in patients with chronic lymphocytic leukemia (CLL). We show that the patterns of SHM and CSR in B-cells from healthy individuals are distinct from CLL. We identify distinct properties of clonal expansion that lead to the generation of antibodies of different classes in healthy, malignant, and non-malignant CLL BCR repertoires. We further demonstrate that BCR diversity is affected by relationships between antibody variable and constant regions leading to isotype-specific signatures of variable gene usage. This study provides powerful insights into the mechanisms underlying the evolution of the adaptive immune responses in health and their aberration during disease.

Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

As part of a broader collaborative network of exome sequencing studies, we developed a jointly called data set of 5,685 Ashkenazi Jewish exomes. We make publicly available a resource of site and allele frequencies, which should serve as a reference for medical genetics in the Ashkenazim (hosted in part at https://ibd.broadinstitute.org, also available in gnomAD at http://gnomad.broadinstitute.org). We estimate that 34% of protein-coding alleles present in the Ashkenazi Jewish population at frequencies greater than 0.2% are significantly more frequent (mean 15-fold) than their maximum frequency observed in other reference populations. Arising via a well-described founder effect approximately 30 generations ago, this catalog of enriched alleles can contribute to differences in genetic risk and overall prevalence of diseases between populations. As validation we document 148 AJ enriched protein-altering alleles that overlap with "pathogenic" ClinVar alleles (table available at https://github.com/macarthur-lab/clinvar/blob/master/output/clinvar.tsv), including those that account for 10-100 fold differences in prevalence between AJ and non-AJ populations of some rare diseases, especially recessive conditions, including Gaucher disease (GBA, p.Asn409Ser, 8-fold enrichment); Canavan disease (ASPA, p.Glu285Ala, 12-fold enrichment); and Tay-Sachs disease (HEXA, c.1421+1G>C, 27-fold enrichment; p.Tyr427IlefsTer5, 12-fold enrichment). We next sought to use this catalog, of well-established relevance to Mendelian disease, to explore Crohn's disease, a common disease with an estimated two to four-fold excess prevalence in AJ. We specifically attempt to evaluate whether strong acting rare alleles, particularly protein-truncating or otherwise large effect-size alleles, enriched by the same founder-effect, contribute excess genetic risk to Crohn's disease in AJ, and find that ten rare genetic risk factors in NOD2 and LRRK2 are enriched in AJ (p < 0.005), including several novel contributing alleles, show evidence of association to CD. Independently, we find that genomewide common variant risk defined by GWAS shows a strong difference between AJ and non-AJ European control population samples (0.97 s.d. higher, p<10-16). Taken together, the results suggest coordinated selection in AJ population for higher CD risk alleles in general. The results and approach illustrate the value of exome sequencing data in case-control studies along with reference data sets like ExAC (sites VCF available via FTP at ftp.broadinstitute.org/pub/ExAC_release/release0.3/) to pinpoint genetic variation that contributes to variable disease predisposition across populations.

Fine-mapping inflammatory bowel disease loci to single-variant resolution.

Inflammatory bowel diseases are chronic gastrointestinal inflammatory disorders that affect millions of people worldwide. Genome-wide association studies have identified 200 inflammatory bowel disease-associated loci, but few have been conclusively resolved to specific functional variants. Here we report fine-mapping of 94 inflammatory bowel disease loci using high-density genotyping in 67,852 individuals. We pinpoint 18 associations to a single causal variant with greater than 95% certainty, and an additional 27 associations to a single variant with greater than 50% certainty. These 45 variants are significantly enriched for protein-coding changes (n = 13), direct disruption of transcription-factor binding sites (n = 3), and tissue-specific epigenetic marks (n = 10), with the last category showing enrichment in specific immune cells among associations stronger in Crohn's disease and in gut mucosa among associations stronger in ulcerative colitis. The results of this study suggest that high-resolution fine-mapping in large samples can convert many discoveries from genome-wide association studies into statistically convincing causal variants, providing a powerful substrate for experimental elucidation of disease mechanisms.

Exploring the genetic architecture of inflammatory bowel disease by whole-genome sequencing identifies association at ADCY7.

To further resolve the genetic architecture of the inflammatory bowel diseases ulcerative colitis and Crohn's disease, we sequenced the whole genomes of 4,280 patients at low coverage and compared them to 3,652 previously sequenced population controls across 73.5 million variants. We then imputed from these sequences into new and existing genome-wide association study cohorts and tested for association at ∼12 million variants in a total of 16,432 cases and 18,843 controls. We discovered a 0.6% frequency missense variant in ADCY7 that doubles the risk of ulcerative colitis. Despite good statistical power, we did not identify any other new low-frequency risk variants and found that such variants explained little heritability. We detected a burden of very rare, damaging missense variants in known Crohn's disease risk genes, suggesting that more comprehensive sequencing studies will continue to improve understanding of the biology of complex diseases.

Genome-wide association study identifies distinct genetic contributions to prognosis and susceptibility in Crohn's disease.

For most immune-mediated diseases, the main determinant of patient well-being is not the diagnosis itself but instead the course that the disease takes over time (prognosis). Prognosis may vary substantially between patients for reasons that are poorly understood. Familial studies support a genetic contribution to prognosis, but little evidence has been found for a proposed association between prognosis and the burden of susceptibility variants. To better characterize how genetic variation influences disease prognosis, we performed a within-cases genome-wide association study in two cohorts of patients with Crohn's disease. We identified four genome-wide significant loci, none of which showed any association with disease susceptibility. Conversely, the aggregated effect of all 170 disease susceptibility loci was not associated with disease prognosis. Together, these data suggest that the genetic contribution to prognosis in Crohn's disease is largely independent of the contribution to disease susceptibility and point to a biology of prognosis that could provide new therapeutic opportunities.

Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations.

Ulcerative colitis and Crohn's disease are the two main forms of inflammatory bowel disease (IBD). Here we report the first trans-ancestry association study of IBD, with genome-wide or Immunochip genotype data from an extended cohort of 86,640 European individuals and Immunochip data from 9,846 individuals of East Asian, Indian or Iranian descent. We implicate 38 loci in IBD risk for the first time. For the majority of the IBD risk loci, the direction and magnitude of effect are consistent in European and non-European cohorts. Nevertheless, we observe genetic heterogeneity between divergent populations at several established risk loci driven by differences in allele frequency (NOD2) or effect size (TNFSF15 and ATG16L1) or a combination of these factors (IL23R and IRGM). Our results provide biological insights into the pathogenesis of IBD and demonstrate the usefulness of trans-ancestry association studies for mapping loci associated with complex diseases and understanding genetic architecture across diverse populations.

Generation of primary human intestinal T cell transcriptomes reveals differential expression at genetic risk loci for immune-mediated disease.

OBJECTIVE: Genome-wide association studies (GWAS) have identified genetic variants within multiple risk loci as predisposing to intestinal inflammatory diseases, including Crohn's disease, ulcerative colitis and coeliac disease. Most risk variants affect regulation of transcription, but a critical challenge is to identify which genes and which cell types these variants affect. We aimed to characterise whole transcriptomes for each common T lymphocyte subset resident within the gut mucosa, and use these to infer biological insights and highlight candidate genes of interest within GWAS risk loci. DESIGN: We isolated the four major intestinal T cell populations from pinch biopsies from healthy subjects and generated transcriptomes for each. We computationally integrated these transcriptomes with GWAS data from immune-related diseases. RESULTS: Robust, high quality transcriptomic data were generated from 1 ng of RNA from precisely sorted cell subsets. Gene expression patterns clearly differentiated intestinal T cells from counterparts in peripheral blood and revealed distinct signalling pathways for each intestinal T cell subset. Intestinal-specific T cell transcripts were enriched in GWAS risk loci for Crohn's disease, ulcerative colitis and coeliac disease, but also specific extraintestinal immune-mediated diseases, allowing prediction of novel candidate genes. CONCLUSIONS: This is the first report of transcriptomes for minimally manipulated intestinal T lymphocyte subsets in humans. We have demonstrated that careful processing of mucosal biopsies allows the generation of transcriptomes from as few as 1000 highly purified cells with minimal interindividual variation. Bioinformatic integration of transcriptomic data with recent GWAS data identified specific candidate genes and cell types for inflammatory pathologies.

Host genetic variants and gene expression patterns associated with Epstein-Barr virus copy number in lymphoblastoid cell lines.

Lymphoblastoid cell lines (LCLs) are commonly used in molecular genetics, supplying DNA for the HapMap and 1000 Genomes Projects, used to test chemotherapeutic agents, and informing the basis of a number of population genetics studies of gene expression. The process of transforming human B cells into LCLs requires the presence of Epstein-Barr virus (EBV), a double-stranded DNA virus which through B-cell immortalisation maintains an episomal virus genome in every cell of an LCL at variable copy numbers. Previous studies have reported that EBV alters host-gene expression and EBV copy number may be under host genetic control. We performed a genome-wide association study of EBV genome copy number in LCLs and found the phenotype to be highly heritable, although no individual SNPs achieved a significant association with EBV copy number. The expression of two host genes (CXCL16 and AGL) was positively correlated and expression of ADARB2 was negatively correlated with EBV copy number in a genotype-independent manner. This study shows an association between EBV copy number and the gene expression profile of LCLs, and suggests that EBV copy number should be considered as a covariate in future studies of host gene expression in LCLs.