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

Plate-Based Assays for the Characterization of Mitochondrial and Cellular Phenotypes

The mitochondria are essential to eukaryotic life, acting as key drivers of energy generation while also being involved in the regulation of many cellular processes including apoptosis, cell proliferation, calcium homeostasis, and metabolism. Mitochondrial diseases which disrupt these processes lead to a diverse range of pathologies and lack consistency in symptom presentation. In disease, mitochondrial activity and energy homeostasis can be adapted to cellular requirements, and studies using Dictyostelium and human lymphoblastoid cell lines have shown that such changes can be facilitated by the key cellular and energy regulators, TORC1 and AMPK. Fluorescence-based assays are increasingly utilized to measure mitochondrial and cell signalling function in mitochondrial disease research. Here, we describe a streamlined method for the simultaneous measurement of mitochondrial mass, membrane potential, and reactive oxygen species production using MitoTracker Green™ FM, MitoTracker Red™ CMXRos, and DCFH-DA probes. This protocol has been adapted for both Dictyostelium and human lymphoblastoid cell lines. We also describe a method for assessing TORC1 and AMPK activity simultaneously in lymphoblastoid cells. These techniques allow for the characterization of mitochondrial defects in a rapid and easy to implement manner.

OPUSeq simplifies detection of low-frequency DNA variants and uncovers fragmentase-associated artifacts

Detection of low-frequency DNA variants (below 1%) is becoming increasingly important in biomedical research and clinical practice, but is challenging to do with standard sequencing approaches due to high error rates. The use of double-stranded unique molecular identifiers (dsUMIs) allows correction of errors by comparing reads arising from the same original DNA duplex. However, the implementation of such approaches is still challenging. Here, we present a novel method, one-pot dsUMI sequencing (OPUSeq), which allows incorporation of dsUMIs in the same reaction as the library PCR. This obviates the need for adapter pre-synthesis or additional enzymatic steps. OPUSeq can be incorporated into standard DNA library preparation approaches and coupled with hybridization target capture. We demonstrate successful error correction and detection of variants down to allele frequency of 0.01%. Using OPUSeq, we also show that the use of enzymatic fragmentation can lead to the appearance of spurious double-stranded variants, interfering with detection of variant fractions below 0.1%.

The three-dimensional structure of Epstein-Barr virus genome varies by latency type and is regulated by PARP1 enzymatic activity

Epstein-Barr virus (EBV) persists in human B-cells by maintaining its chromatinized episomes within the nucleus. We have previously shown that cellular factor Poly [ADP-ribose] polymerase 1 (PARP1) binds the EBV genome, stabilizes CTCF binding at specific loci, and that PARP1 enzymatic activity correlates with maintaining a transcriptionally active latency program. To better understand PARP1's role in regulating EBV latency, here we functionally characterize the effect of PARP enzymatic inhibition on episomal structure through in situ HiC mapping, generating a complete 3D structure of the EBV genome. We also map intragenomic contact changes after PARP inhibition to global binding of chromatin looping factors CTCF and cohesin across the EBV genome. We find that PARP inhibition leads to fewer total unique intragenomic interactions within the EBV episome, yet new chromatin loops distinct from the untreated episome are also formed. This study also illustrates that PARP inhibition alters gene expression at the regions where chromatin looping is most effected. We observe that PARP1 inhibition does not alter cohesin binding sites but does increase its frequency of binding at those sites. Taken together, these findings demonstrate that PARP has an essential role in regulating global EBV chromatin structure and latent gene expression.

Admixture with indigenous people helps local adaptation: admixture-enabled selection in Polynesians

Background

Homo sapiens have experienced admixture many times in the last few thousand years. To examine how admixture affects local adaptation, we investigated genomes of modern Polynesians, who are shaped through admixture between Austronesian-speaking people from Southeast Asia (Asian-related ancestors) and indigenous people in Near Oceania (Papuan-related ancestors).

Methods

In this study local ancestry was estimated across the genome in Polynesians (23 Tongan subjects) to find the candidate regions of admixture-enabled selection contributed by Papuan-related ancestors.

Results

The mean proportion of Papuan-related ancestry across the Polynesian genome was estimated as 24.6% (SD = 8.63%), and two genomic regions, the extended major histocompatibility complex (xMHC) region on chromosome 6 and the ATP-binding cassette transporter sub-family C member 11 (ABCC11) gene on chromosome 16, showed proportions of Papuan-related ancestry more than 5 SD greater than the mean (> 67.8%). The coalescent simulation under the assumption of selective neutrality suggested that such signals of Papuan-related ancestry enrichment were caused by positive selection after admixture (false discovery rate = 0.045). The ABCC11 harbors a nonsynonymous SNP, rs17822931, which affects apocrine secretory cell function. The approximate Bayesian computation indicated that, in Polynesian ancestors, a strong positive selection (s = 0.0217) acted on the ancestral allele of rs17822931 derived from Papuan-related ancestors.

Conclusions

Our results suggest that admixture with Papuan-related ancestors contributed to the rapid local adaptation of Polynesian ancestors. Considering frequent admixture events in human evolution history, the acceleration of local adaptation through admixture should be a common event in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01900-y.

Lymphoblastoid Cell Lines as Models to Study Mitochondrial Function in Neurological Disorders

Neurological disorders, including neurodegenerative diseases, are collectively a major cause of death and disability worldwide. Whilst the underlying disease mechanisms remain elusive, altered mitochondrial function has been clearly implicated and is a key area of study in these disorders. Studying mitochondrial function in these disorders is difficult due to the inaccessibility of brain tissue, which is the key tissue affected in these diseases. To overcome this issue, numerous cell models have been used, each providing unique benefits and limitations. Here, we focussed on the use of lymphoblastoid cell lines (LCLs) to study mitochondrial function in neurological disorders. LCLs have long been used as tools for genomic analyses, but here we described their use in functional studies specifically in regard to mitochondrial function. These models have enabled characterisation of the underlying mitochondrial defect, identification of altered signalling pathways and proteins, differences in mitochondrial function between subsets of particular disorders and identification of biomarkers of the disease. The examples provided here suggest that these cells will be useful for development of diagnostic tests (which in most cases do not exist), identification of drug targets and testing of pharmacological agents, and are a worthwhile model for studying mitochondrial function in neurological disorders.

Antibody-Mediated Immune Subset Depletion Modulates the Immune Response in a Rabbit (Oryctolagus cuniculus) Model of Epstein-Barr Virus Infection

Epstein-Barr Virus (EBV) is a γ-herpesvirus which infects over 90% of the adult human population. Most notably, this virus causes infectious mononucleosis but it is also associated with cancers such as Hodgkin and Burkitt lymphoma. EBV is a species-specific virus and has been studied in many animal models, including nonhuman primates, guinea pigs, humanized mice, and tree shrews. However, none of these animal models are considered the "gold standard" for EBV research. Recently, rabbits have emerged as a viable alternative model, as they are susceptible to EBV infection. In addition, the EBV infection progresses after immune suppression with cyclosporine A (CsA), modeling the reactivation of EBV after latency. We sought to refine this model for acute or active EBV infections by performing antibody-mediated depletion of certain immune subsets in rabbits. Fourteen 16 to 20-wk old, NZW rabbits were intravenously inoculated with EBV and concurrently treated with either anti-CD4 T-cell antibody, anti-pan-T-cell antibody (anti CD45), CSA, or, as a control, anti-HPV antibody. Rabbits that received the depleting antibodies were treated with CsA 3 times at a dose of 15 mg/kg SC once per day for 4 d starting at the time of EBV inoculation then the dose was increased to 20 mg/kg SC twice weekly for 2 wk. Weights, temperatures, and clinical signs were monitored, and rabbits were anesthetized once weekly for blood collection. When compared with the control group, anti-CD4-treated rabbits had fewer clinical signs and displayed higher levels of viral DNA via qPCR in splenocytes; however, flow cytometry results showed only a partial depletion of CD4 T-cells. Treatment with anti-pan-T-cell antibody did not result in noticeable T-cell depletion. These data suggest the EBV-infected rabbit is a promising model for testing antiviral medications and prophylactic vaccines for EBV.

The interaction of Multiple Sclerosis risk loci with Epstein-Barr virus phenotypes implicates the virus in pathogenesis

Translating the findings of genome wide association studies (GWAS) to new therapies requires identification of the relevant immunological contexts to interrogate for genetic effects. In one of the largest GWAS, more than 200 risk loci have been identified for Multiple Sclerosis (MS) susceptibility. Infection with Epstein-Barr virus (EBV) appears to be necessary for the development of Multiple Sclerosis (MS). Many MS risk loci are associated with altered gene expression in EBV infected B cells (LCLs). We have interrogated this immunological context to identify interaction between MS risk loci and EBV DNA copy number, intrinsic growth rate and EBV encoded miRNA expression. The EBV DNA copy number was associated with significantly more risk alleles for MS than for other diseases or traits. EBV miRNAs BART4-3p and BART3-5p were highly associated with EBV DNA copy number and MS risk loci. The poliovirus receptor (PVR) risk SNP was associated with EBV DNA copy number, PVR and miRNA expression. Targeting EBV miRNAs BART4-3p and BART3-5p, and the gene PVR, may provide therapeutic benefit in MS. This study also indicates how immunological context and risk loci interactions can be exploited to validate and develop novel therapeutic approaches.

Evidence from genome wide association studies implicates reduced control of Epstein-Barr virus infection in multiple sclerosis susceptibility

BACKGROUND

Genome wide association studies have identified > 200 susceptibility loci accounting for much of the heritability of multiple sclerosis (MS). Epstein-Barr virus (EBV), a memory B cell tropic virus, has been identified as necessary but not sufficient for development of MS. The molecular and immunological basis for this has not been established. Infected B cell proliferation is driven by signalling through the EBV produced cell surface protein LMP1, a homologue of the MS risk gene CD40.

METHODS

We have investigated transcriptomes of B cells and EBV-infected B cells at Latency III (LCLs) and identified MS risk genes with altered expression on infection and with expression levels associated with the MS risk genotype (LCLeQTLs). The association of LCLeQTL genomic burden with EBV phenotypes in vitro and in vivo was examined. The risk genotype effect on LCL proliferation with CD40 stimulation was assessed.

RESULTS

These LCLeQTL MS risk SNP:gene pairs (47 identified) were over-represented in genes dysregulated between B and LCLs (p < 1.53 × 10^-4), and as target loci of the EBV transcription factor EBNA2 (p < 3.17 × 10^-16). Overall genetic burden of LCLeQTLs was associated with some EBV phenotypes but not others. Stimulation of the CD40 pathway by CD40L reduced LCL proliferation (p < 0.001), dependent on CD40 and TRAF3 MS risk genotypes. Both CD40 and TRAF3 risk SNPs are in binding sites for the EBV transcription factor EBNA2, with expression of each correlated with EBNA2 expression dependent on genotype.

CONCLUSIONS

These data indicate targeting EBV may be of therapeutic benefit in MS.

Extensive epigenetic and transcriptomic variability between genetically identical human B-lymphoblastoid cells with implications in pharmacogenomics research

Genotyped human B-lymphoblastoid cell lines (LCLs) are widely used models in mapping quantitative trait loci for chromatin features, gene expression, and drug response. The extent of genotype-independent functional genomic variability of the LCL model, although largely overlooked, may inform association study design. In this study, we use flow cytometry, chromatin immunoprecipitation sequencing and mRNA sequencing to study surface marker patterns, quantify genome-wide chromatin changes (H3K27ac) and transcriptome variability, respectively, among five isogenic LCLs derived from the same individual. Most of the studied LCLs were non-monoclonal and had mature B cell phenotypes. Strikingly, nearly one-fourth of active gene regulatory regions showed significantly variable H3K27ac levels, especially enhancers, among which several were classified as clustered enhancers. Large, contiguous genomic regions showed signs of coordinated activity change. Regulatory differences were mirrored by mRNA expression changes, preferentially affecting hundreds of genes involved in specialized cellular processes including immune and drug response pathways. Differential expression of DPYD, an enzyme involved in 5-fluorouracil (5-FU) catabolism, was associated with variable LCL growth inhibition mediated by 5-FU. The extent of genotype-independent functional genomic variability might highlight the need to revisit study design strategies for LCLs in pharmacogenomics.

Extended methods for gene-environment-wide interaction scans in studies of admixed individuals with varying degrees of relationships

The etiology of many complex diseases involves both environmental exposures and inherited genetic predisposition as well as interactions between them. Gene-environment-wide interaction studies (GEWIS) provide a means to identify the interactions between genetic variation and environmental exposures that underlie disease risk. However, current GEWIS methods lack the capability to adjust for the potentially complex correlations in studies with varying degrees of relationships (both known and unknown) among individuals in admixed populations. We developed novel generalized estimating equation (GEE) based methods-GEE-adaptive and GEE-joint-to account for phenotypic correlations due to kinship while accounting for covariates, including, measures of genome-wide ancestry. In simulation studies of admixed individuals, both methods controlled family-wise error rates, an advantage over the case-only approach. They demonstrated higher power than traditional case-control methods across a wide range of underlying alternative hypotheses, especially where both marginal and interaction effects were present. We applied the proposed method to conduct a GEWIS of a known sarcoidosis risk factor (insecticide exposure) and risk of sarcoidosis in African Americans and identified two novel loci with suggestive evidence of G × E interaction.

Viral microRNA effects on persistent infection of human lymphoid cells by polyomavirus SV40

BACKGROUND

Polyomaviruses, including simian virus 40 (SV40), display evidence of lymphotropic properties. This study analyzed the nature of SV40-human lymphocyte interactions in established cell lines and in primary lymphocytes. The effects of viral microRNA and the structure of the viral regulatory region on SV40 persistence were examined.

RESULTS

SV40 DNA was maintained in infected B cell and myeloid cell lines during cell growth for at least 28 days. Limiting dilution analysis showed that low amounts of SV40 DNA (~2 copies per cell) were retained over time. Infected B cells remained viable and able to proliferate. Genome copies of the SV40 microRNA-null mutant persisted at higher levels than the DNA of wild-type viruses. Complex viral regulatory regions produced modestly higher DNA levels than simple regulatory regions. Viral large T-antigen protein was detected at low frequency and at low levels in infected B cells. Following infection of primary lymphocytes, SV40 DNA was detected in CD19+ B cells and CD14+ monocytes, but not in CD3+ T cells. Rescue attempts using either lysates of SV40-infected B lymphocytes, coculture of live cells, or infectious center assays all showed that replication-competent SV40 could be recovered on rare occasions. SV40 infections altered the expression of several B cell surface markers, with more pronounced changes following infections with the microRNA-null mutant.

CONCLUSION

These findings indicate that SV40 can establish persistent infections in human B lymphocytes. The cells retain low copy numbers of viral DNA; the infections are nonproductive and noncytolytic but can occasionally produce infectious virus. SV40 microRNA negatively regulates the degree of viral effects on B cells.

SIGNIFICANCE

Lymphocytes may serve as viral reservoirs and may function to disseminate polyomaviruses to different tissues in a host. To our knowledge, this report is the first extensive analysis of viral microRNA effects on SV40 infection of human lymphocytes.

Genetic factors affecting EBV copy number in lymphoblastoid cell lines derived from the 1000 Genome Project samples

Epstein-Barr virus (EBV), human herpes virus 4, has been classically associated with infectious mononucleosis, multiple sclerosis and several types of cancers. Many of these diseases show marked geographical differences in prevalence, which points to underlying genetic and/or environmental factors. Those factors may include a different susceptibility to EBV infection and viral copy number among human populations. Since EBV is commonly used to transform B-cells into lymphoblastoid cell lines (LCLs) we hypothesize that differences in EBV copy number among individual LCLs may reflect differential susceptibility to EBV infection. To test this hypothesis, we retrieved whole-genome sequenced EBV-mapping reads from 1,753 LCL samples derived from 19 populations worldwide that were sequenced within the context of the 1000 Genomes Project. An in silico methodology was developed to estimate the number of EBV copy number in LCLs and validated these estimations by real-time PCR. After experimentally confirming that EBV relative copy number remains stable over cell passages, we performed a genome wide association analysis (GWAS) to try detecting genetic variants of the host that may be associated with EBV copy number. Our GWAS has yielded several genomic regions suggestively associated with the number of EBV genomes per cell in LCLs, unraveling promising candidate genes such as CAND1, a known inhibitor of EBV replication. While this GWAS does not unequivocally establish the degree to which genetic makeup of individuals determine viral levels within their derived LCLs, for which a larger sample size will be needed, it potentially highlighted human genes affecting EBV-related processes, which constitute interesting candidates to follow up in the context of EBV related pathologies.

The copy number of Epstein-Barr virus latent genome correlates with the oncogenicity by the activation level of LMP1 and NF-κB

A tumor model that Epstein-Barr virus (EBV) latent infection facilitated the tumorigenicity was previously established using the Maxi-EBV system. In the present approach, EBV-lost cell clones demonstrated significantly decreased tumorigenesis. On the other hand, the LMP1 gene in Maxi-EBV genome was replaced by that of nasopharyngeal carcinoma origin. The resultant cell line, 293–1/NL showed much lower malignancy than the original 293-EBV. The result was opposite to our expectation. The change of 293 sublineage cells for EBV harboring also got similar result. To seek the underlying reason, the copy number of EBV genome in all the cell lines was detected. The result indicated that 293-EBV contained about 4.5-fold higher EBV copies than 293–1/NL did. Parallel EBV genomes led to relatively stable copies in different 293 sublineages, suggesting the viral genome structure is a factor for the sustainability of EBV's copy number. Moreover, the LMP1 transcription in high copy-containing cells showed abnormally high level. Furthermore, the main LMP1-driven pathway, transcription factor NF-κB, was highly activated in high-copy cells. Here we first manifest by experimental model that the copy number of EBV latent genome correlates with the viral pathogenesis, which depends on the activation level of LMP1 and NF-κB. Overall, both the presence and amount of EBV genome are crucial for the viral oncogenicity.

Coding and noncoding expression patterns associated with rare obesity-related disorders: Prader-Willi and Alström syndromes

Obesity is accompanied by hyperphagia in several classical genetic obesity-related syndromes that are rare, including Prader-Willi syndrome (PWS) and Alström syndrome (ALMS). We compared coding and noncoding gene expression in adult males with PWS, ALMS, and nonsyndromic obesity relative to nonobese males using readily available lymphoblastoid cells to identify disease-specific molecular patterns and disturbed mechanisms in obesity. We found 231 genes upregulated in ALMS compared with nonobese males, but no genes were found to be upregulated in obese or PWS males and 124 genes were downregulated in ALMS. The metallothionein gene (MT1X) was significantly downregulated in ALMS, in common with obese males. Only the complex SNRPN locus was disturbed (downregulated) in PWS along with several downregulated small nucleolar RNAs (snoRNAs) in the 15q11-q13 region (SNORD116, SNORD109B, SNORD109A, SNORD107). Eleven upregulated and ten downregulated snoRNAs targeting multiple genes impacting rRNA processing, developmental pathways, and associated diseases were found in ALMS. Fifty-two miRNAs associated with multiple, overlapping gene expression disturbances were upregulated in ALMS, and four were shared with obese males but not PWS males. For example, seven passenger strand microRNAs (miRNAs) (miR-93*, miR-373*, miR-29b-2*, miR-30c-1*, miR27a*, miR27b*, and miR-149*) were disturbed in association with six separate downregulated target genes (CD68, FAM102A, MXI1, MYO1D, TP53INP1, and ZRANB1). Cell cycle (eg, PPP3CA), transcription (eg, POLE2), and development may be impacted by upregulated genes in ALMS, while downregulated genes were found to be involved with metabolic processes (eg, FABP3), immune responses (eg, IL32), and cell signaling (eg, IL1B). The high number of gene and noncoding RNA disturbances in ALMS contrast with observations in PWS and males with nonsyndromic obesity and may reflect the progressing multiorgan pathology of the ALMS disease process.