Phylogenetic comparison of Epstein-Barr virus genomes

Technologies used for genome analysis and whole genome sequencing are useful for us to understand genomic characterization and divergence. The Epstein-Barr virus (EBV) is an oncogenic virus that causes diverse diseases such as Burkitt’s lymphoma (BL), nasopharyngeal carcinoma (NPC), Hodgkin’s lymphoma (HL), and gastric carcinoma (GC). EBV genomes found in these diseases can be classified either by phases of EBV latency (type-I, -II, and -III latency) or types of EBNA2 sequence difference (type-I EBV, type-II EBV or EBV-1, EBV-2). EBV from EBV-transformed lymphoblastoid cell line (LCL) establishes type-III latency, EBV from NPC establishes type-II latency, and EBV from GC establishes type-I latency. However, other important factors play key roles in classifying numerous EBV strains because EBV genomes are highly diverse and not phylogenetically related to types of EBV-associated diseases. Herein, we first reviewed previous studies to describe molecular characteristics of EBV genomes. Then, using comparative and phylogenetic analyses, we phylogenetically analyzed molecular variations of EBV genomes and proteins. The review of previous studies and our phylogenetic analysis showed that EBV genomes and proteins were highly diverse regardless of types of EBV-associated diseases. Other factors should be considered in determining EBV taxonomy. This review will be helpful to understand complicated phylogenetic relationships of EBV genomes.     

Epstein–Barr Virus Genetic Variation in Lymphoblastoid Cell Lines Derived from Kenyan Pediatric Population

Epstein-Barr virus (EBV) is associated with Burkitt’s lymphoma (BL), and in regions of sub-Saharan Africa where endemic BL is common, both the EBV Type 1 (EBV-1) and EBV Type 2 strains (EBV-2) are found. Little is known about genetic variation of EBV strains in areas of sub-Saharan Africa. In the present study, spontaneous lymphoblastoid cell lines (LCLs) were generated from samples obtained from Kenya. Polymerase chain reaction (PCR) amplification of the EBV genome was done using multiple primers and sequenced by next-generation sequencing (NGS). Phylogenetic analyses against the published EBV-1 and EBV-2 strains indicated that one sample, LCL10 was closely related to EBV-2, while the remaining 3 LCL samples were more closely related to EBV-1. Moreover, single nucleotide polymorphism (SNP) analyses showed clustering of LCL variants. We further show by analysis of EBNA-1, BLLF1, BPLF1, and BRRF2 that latent genes are less conserved than lytic genes in these LCLs from a single geographic region. In this study we have shown that NGS is highly useful for deciphering detailed inter and intra-variations in EBV genomes and that within a geographic region different EBV genetic variations can co-exist, the implications of which warrant further investigation. The findings will enhance our understanding of potential pathogenic variants critical to the development and maintenance of EBV-associated malignancies.     

Small‐ and large‐scale heterogeneity in genetic variation across the collard flycatcher genome: implications for estimating genetic diversity in nonmodel organisms

Population genetic studies in nonmodel organisms are often hampered by a lack of reference genomes that are essential for whole‐genome resequencing. In the light of this, genotyping methods have been developed to effectively eliminate the need for a reference genome, such as genotyping by sequencing or restriction site‐associated DNA sequencing (RAD‐seq). However, what remains relatively poorly studied is how accurately these methods capture both average and variation in genetic diversity across an organism's genome. In this issue of Molecular Ecology Resources, Dutoit et al. (2016) use whole‐genome resequencing data from the collard flycatcher to assess what factors drive heterogeneity in nucleotide diversity across the genome. Using these data, they then simulate how well different sequencing designs, including RAD sequencing, could capture most of the variation in genetic diversity. They conclude that for evolutionary and conservation‐related studies focused on the estimating genomic diversity, researchers should emphasize the number of loci analysed over the number of individuals sequenced.     

Draft Sequences of the Radish (Raphanus sativus L.) Genome

Radish (Raphanus sativus L., n = 9) is one of the major vegetables in Asia. Since the genomes of Brassica and related species including radish underwent genome rearrangement, it is quite difficult to perform functional analysis based on the reported genomic sequence of Brassica rapa. Therefore, we performed genome sequencing of radish. Short reads of genomic sequences of 191.1 Gb were obtained by next-generation sequencing (NGS) for a radish inbred line, and 76,592 scaffolds of ≥300 bp were constructed along with the bacterial artificial chromosome-end sequences. Finally, the whole draft genomic sequence of 402 Mb spanning 75.9% of the estimated genomic size and containing 61,572 predicted genes was obtained. Subsequently, 221 single nucleotide polymorphism markers and 768 PCR-RFLP markers were used together with the 746 markers produced in our previous study for the construction of a linkage map. The map was combined further with another radish linkage map constructed mainly with expressed sequence tag-simple sequence repeat markers into a high-density integrated map of 1,166 cM with 2,553 DNA markers. A total of 1,345 scaffolds were assigned to the linkage map, spanning 116.0 Mb. Bulked PCR products amplified by 2,880 primer pairs were sequenced by NGS, and SNPs in eight inbred lines were identified.     

Sequences of the Epstein-Barr Virus (EBV) large internal repeat form the center of a 16-kilobase-pair palindrome of EBV (P3HR-1) heterogeneous DNA.

We have previously characterized several genomic rearrangements of Epstein-Barr virus (EBV) DNA contained in one of the defective EBV genomes harbored by the P3HR-1 (HR-1) line (H. B. Jenson, M. S. Rabson, and G. Miller, J. Virol. 58:475-486, 1986). One recombinant clone of heterogeneous DNA (het DNA) from this defective genome is an EcoRI fragment of 16 kilobase pairs (kbp) which is a palindrome. DNA digestion fragments specific for the center of this palindrome were present in cells which contained het DNA but not in cells which lacked het DNA. Thus, the palindrome was not an artifact of DNA cloning. The organization of the center of this palindrome was studied by DNA sequencing. The comparable region of the parental HR-1 genome was also studied by DNA sequencing. The central 3,495 base pairs (bp) of the palindrome were composed of sequences derived exclusively from internal repeat 1 of EBV, represented by BamHI W fragment. At each end of the central 3,495 hp was a symmetrical recombination with sequences of BamHI-Z, located more than 50 kbp away on the standard EBV genome. The central 3,495 bp were composed of an unduplicated 341 bp flanked by two perfect palindromic repeats of 1,577 bp. The 341-bp unique region was a portion of a 387-bp region of standard HR-1 BamHI-W which was identical to the central 387 bp of the palindrome. This central 387-bp region contained numerous stretches of dyad symmetry capable of forming a large stem-and-loop structure. The palindromic rearrangement had created two novel open reading frames in het DNA derived from standard HR-1 BamHI-W sequences. These two het DNA open reading frames had different amino termini but identical carboxy termini derived from the large open reading frame in standard HR-1 BamHI-W (HR-1 BWRF1). The BamHI-W sequences found in het DNA did not include either the TATA box of standard HR-1 BamHI-W or the exons which are present in the potentially polycistronic latent mRNAs encoding EBV nuclear antigens. These marked alterations in genomic structure may relate to the unique biologic properties of virus stocks containing het DNA by creation of new polypeptides or by formation or deletion of regulatory or functional signals.     

Selective Whole Genome Amplification for Resequencing Target Microbial Species from Complex Natural Samples

Population genomic analyses have demonstrated power to address major questions in evolutionary and molecular microbiology. Collecting populations of genomes is hindered in many microbial species by the absence of a cost effective and practical method to collect ample quantities of sufficiently pure genomic DNA for next-generation sequencing. Here we present a simple method to amplify genomes of a target microbial species present in a complex, natural sample. The selective whole genome amplification (SWGA) technique amplifies target genomes using nucleotide sequence motifs that are common in the target microbe genome, but rare in the background genomes, to prime the highly processive phi29 polymerase. SWGA thus selectively amplifies the target genome from samples in which it originally represented a minor fraction of the total DNA. The post-SWGA samples are enriched in target genomic DNA, which are ideal for population resequencing. We demonstrate the efficacy of SWGA using both laboratory-prepared mixtures of cultured microbes as well as a natural host–microbe association. Targeted amplification of Borrelia burgdorferi mixed with Escherichia coli at genome ratios of 1:2000 resulted in >10⁵-fold amplification of the target genomes with <6.7-fold amplification of the background. SWGA-treated genomic extracts from Wolbachia pipientis-infected Drosophila melanogaster resulted in up to 70% of high-throughput resequencing reads mapping to the W. pipientis genome. By contrast, 2–9% of sequencing reads were derived from W. pipientis without prior amplification. The SWGA technique results in high sequencing coverage at a fraction of the sequencing effort, thus allowing population genomic studies at affordable costs.     

Whole genome amplification and sequencing of a Daphnia resting egg

Resting eggs banks are unique windows that allow us to directly observe shifts in population genetics, and phenotypes over time as natural populations evolve. Though a variety of planktonic organisms also produce resting stages, the keystone freshwater consumer, Daphnia, is a well‐known model for paleogenetics and resurrection ecology. Nevertheless, paleogenomic investigations are limited largely because resting eggs do not contain enough DNA for genomic sequencing. In fact, genomic studies even on extant populations include a laborious preparatory phase of batch culturing dozens of individuals to generate sufficient genomic DNA. Here, we furnish a protocol to generate whole genomes of single ephippial (resting) eggs and single daphniids. Whole genomes of single ephippial eggs and single adults were amplified using Qiagen REPLI‐g Single Cell kit reaction, followed by NEBNext Ultra DNA Library Prep Kit for library construction and Illumina sequencing. We compared the quality of the single‐egg and single‐individual amplified genomes to the standard batch genomic DNA extraction in the absence of genome amplification. At mean 20× depth, coverage was essentially identical for the amplified single individual relative to the unamplified batch extracted genome (>90% of the genome was covered and callable). Finally, while amplification resulted in the slight loss of heterozygosity for the amplified genomes, estimates were largely comparable and illustrate the utility and limitations of this approach in estimating population genetic parameters over long periods of time in natural populations of Daphnia and also other small species known to produce resting stages.     

Adventures in the enormous: a 1.8 million clone BAC library for the 21.7 Gb genome of loblolly pine

Loblolly pine (LP; Pinus taeda L.) is the most economically important tree in the U.S. and a cornerstone species in southeastern forests. However, genomics research on LP and other conifers has lagged behind studies on flowering plants due, in part, to the large size of conifer genomes. As a means to accelerate conifer genome research, we constructed a BAC library for the LP genotype 7-56. The LP BAC library consists of 1,824,768 individually-archived clones making it the largest single BAC library constructed to date, has a mean insert size of 96 kb, and affords 7.6X coverage of the 21.7 Gb LP genome. To demonstrate the efficacy of the library in gene isolation, we screened macroarrays with overgos designed from a pine EST anchored on LP chromosome 10. A positive BAC was sequenced and found to contain the expected full-length target gene, several gene-like regions, and both known and novel repeats. Macroarray analysis using the retrotransposon IFG-7 (the most abundant repeat in the sequenced BAC) as a probe indicates that IFG-7 is found in roughly 210,557 copies and constitutes about 5.8% or 1.26 Gb of LP nuclear DNA; this DNA quantity is eight times the Arabidopsis genome. In addition to its use in genome characterization and gene isolation as demonstrated herein, the BAC library should hasten whole genome sequencing of LP via next-generation sequencing strategies/technologies and facilitate improvement of trees through molecular breeding and genetic engineering. The library and associated products are distributed by the Clemson University Genomics Institute (www.genome.clemson.edu).     

An efficient and tunable parameter to improve variant calling for whole genome and exome sequencing data

Next generation sequencing (NGS) has traditionally been performed in various fields including agricultural to clinical and there are so many sequencing platforms available in order to obtain accurate and consistent results. However, these platforms showed amplification bias when facilitating variant calls in personal genomes. Here, we sequenced whole genomes and whole exomes from ten Korean individuals using Illumina and Ion Proton, respectively to find the vulnerability and accuracy of NGS platform in the GC rich/poor area. Overall, a total of 1013 Gb reads from Illumina and ~39.1 Gb reads from Ion Proton were analyzed using BWA-GATK variant calling pipeline. Furthermore, conjunction with the VQSR tool and detailed filtering strategies, we achieved high-quality variants. Finally, each of the ten variants from Illumina only, Ion Proton only, and intersection was selected for Sanger validation. The validation results revealed that Illumina platform showed higher accuracy than Ion Proton. The described filtering methods are advantageous for large population-based whole genome studies designed to identify common and rare variations associated with complex diseases.     

Sustained viral activity of epstein-Barr virus contributes to cellular immortalization of lymphoblastoid cell lines

EBV-transformed lymphoblastoid cell lines (LCLs) are used as a resource for human genetic, immunological, and pharmacogenomic studies. We investigated the biological activity of 20 LCL strains during continuous long-term subculture up to a passage number of 160. Out of 20 LCL strains, 17 proliferated up to a passage number of 160, at which point LCLs are generally considered as “immortalized”. The other three LCL strains lost the ability to proliferate at an average passage number of 41, during which these LCLs may have undergone cellular crisis. These non-immortal LCL strains exhibited no telomerase activity, decreased EBV gene expression, and a lower copy number of the EBV genome and mitochondrial DNA when compared with immortal LCLs. Thus, this study suggests that sustained EBV viral activity as well as telomerase activity may be required for complete LCL immortalization. These authors contributed equally to this work.     

Exploration of Genetic Variations through Single‐cell Whole‐genome Sequencing in the Model Ciliate Tetrahymena thermophila

Ciliates are unicellular eukaryotes with separate germline and somatic genomes and diverse life cycles, which make them a unique model to improve our understanding of population genetics through the detection of genetic variations. However, traditional sequencing methods cannot be directly applied to ciliates because the majority are uncultivated. Single‐cell whole‐genome sequencing (WGS) is a powerful tool for studying genetic variation in microbes, but no studies have been performed in ciliates. We compared the use of single‐cell WGS and bulk DNA WGS to detect genetic variation, specifically single nucleotide polymorphisms (SNPs), in the model ciliate Tetrahymena thermophila. Our analyses showed that (i) single‐cell WGS has excellent performance regarding mapping rate and genome coverage but lower sequencing uniformity compared with bulk DNA WGS due to amplification bias (which was reproducible); (ii) false‐positive SNP sites detected by single‐cell WGS tend to occur in genomic regions with particularly high sequencing depth and high rate of C:G to T:A base changes; (iii) SNPs detected in three or more cells should be reliable (an detection efficiency of 83.4–97.4% was obtained for combined data from three cells). This analytical method could be adapted to measure genetic variation in other ciliates and broaden research into ciliate population genetics.     

Fidelity of SNP Array Genotyping Using Epstein Barr Virus-Transformed B-Lymphocyte Cell Lines: Implications for Genome-Wide Association Studies

Background

As availability of primary cells can be limited for genetic studies of human disease, lymphoblastoid cell lines (LCL) are common sources of genomic DNA. LCL are created in a transformation process that entails in vitro infection of human B-lymphocytes with the Epstein-Barr Virus (EBV).

Methodology/Principal Findings

To test for genotypic errors potentially induced by the Epstein-Barr Virus transformation process, we compared single nucleotide polymorphism (SNP) genotype calls in peripheral blood mononuclear cells (PBMC) and LCL from the same individuals. The average mismatch rate across 19 comparisons was 0.12% for SNPs with a population call rate of at least 95%, and 0.03% at SNPs with a call rate of at least 99%. Mismatch rates were not correlated across genotype subarrays run on all sample pairs.

Conclusions/Significance

Genotypic discrepancies found in PBMC and LCL pairs were not significantly different than control pairs, and were not correlated across subarrays. These results suggest that mismatch rates are minimal with stringent quality control, and that most genotypic discrepancies are due to technical artifacts rather than the EBV transformation process. Thus, LCL likely constitute a reliable DNA source for host genotype analysis.     

Global mapping of transposon location

Transposable genetic elements are ubiquitous, yet their presence or absence at any given position within a genome can vary between individual cells, tissues, or strains. Transposable elements have profound impacts on host genomes by altering gene expression, assisting in genomic rearrangements, causing insertional mutations, and serving as sources of phenotypic variation. Characterizing a genome's full complement of transposons requires whole genome sequencing, precluding simple studies of the impact of transposition on interindividual variation. Here, we describe a global mapping approach for identifying transposon locations in any genome, using a combination of transposon-specific DNA extraction and microarray-based comparative hybridization analysis. We use this approach to map the repertoire of endogenous transposons in different laboratory strains of Saccharomyces cerevisiae and demonstrate that transposons are a source of extensive genomic variation. We also apply this method to mapping bacterial transposon insertion sites in a yeast genomic library. This unique whole genome view of transposon location will facilitate our exploration of transposon dynamics, as well as defining bases for individual differences and adaptive potential.     

Sequence analysis of genes and genomes

A major step towards understanding of the genetic basis of an organism is the complete sequence determination of all genes in its genome. The development of powerful techniques for DNA sequencing has enabled sequencing of large amounts of gene fragments and even complete genomes. Important new techniques for physical mapping, DNA sequencing and sequence analysis have been developed. To increase the throughput, automated procedures for sample preparation and new software for sequence analysis have been applied. This review describes the development of new sequencing methods and the optimisation of sequencing strategies for whole genome and cDNA analysis, as well as discusses issues regarding sequence analysis and annotation.     

Molecular characterization of transgenic shallots (Allium cepa L.) by adaptor ligation PCR (AL-PCR) and sequencing of genomic DNA flanking T-DNA borders

Genomic DNA blot hybridization is traditionally used to demonstrate that, via genetic transformation, foreign genes are integrated into host genomes. However, in large genome species, such as Allium cepa L., the use of genomic DNA blot hybridization is pushed towards its limits, because a considerable quantity of DNA is needed to obtain enough genome copies for a clear hybridization pattern. Furthermore, genomic DNA blot hybridization is a time-consuming method. Adaptor ligation PCR (AL-PCR) of genomic DNA flanking T-DNA borders does not have these drawbacks and seems to be an adequate alternative to genomic DNA blot hybridization. Using AL-PCR we proved that T-DNA was integrated into the A. cepa genome of three transgenic lines transformed with Agrobacterium tumefaciens EHA105 (pCAMBIA 1301). The AL-PCR patterns obtained were specific and reproducible for a given transgenic line. The results showed that T-DNA integration took place and gave insight in the number of T-DNA copies present. Comparison of AL-PCR and previously obtained genomic DNA blot hybridization results pointed towards complex T-DNA integration patterns in some of the transgenic plants. After cloning and sequencing the AL-PCR products, the junctions between plant genomic DNA and the T-DNA insert could be analysed in great detail. For example it was shown that upon T-DNA integration a 66bp genomic sequence was deleted, and no filler DNA was inserted. Primers located within the left and right flanking genomic DNA in transgenic shallot plants were used to recover the target site of T-DNA integration.     

Bread matters: a national initiative to profile the genetic diversity of Australian wheat

The large and complex genome of wheat makes genetic and genomic analysis in this important species both expensive and resource intensive. The application of next-generation sequencing technologies is particularly resource intensive, with at least 17?Gbp of sequence data required to obtain minimal (1×) coverage of the genome. A similar volume of data would represent almost 40× coverage of the rice genome. Progress can be made through the establishment of consortia to produce shared genomic resources. Australian wheat genome researchers, working with Bioplatforms Australia, have collaborated in a national initiative to establish a genetic diversity dataset representing Australian wheat germplasm based on whole genome next-generation sequencing data. Here, we describe the establishment and validation of this resource which can provide a model for broader international initiatives for the analysis of large and complex genomes.     

Human blymphocytes immortalization by epstein-barr virus in the presence of cyclosporin a

Summary A simple method was devised for the establishment of continuous lymphoblastoid, cell lines (LCL). Unfractionated mononuclear cells collected from healthy donors were infected in vitro by Epstein-Barr Virus (EBV) (strain B95-8) under specific, conditions: an immunosuppressive drug, Cyclosporin-A (CS-A, Sandoz), associated with the use of a feeder layer (MRC-5) led to 100% efficiency of LCL establishment. A bank of 400 LCL was set up for completion of genetic studies. Regression and kinetics of virus-induced transformation were monitored and related to donors' EBV immune status. Mean time of LCL establishment and probability of regression among seropositive donors were not linked to any given value titer of antibodies against Viral Capsid Antigen (VCA) or against Epstein-Barr Nuclear Antigen (EBNA). However, when the anti-VCA:anti-EBNA ratio was considered, this parameter, seemed to be linked to the kinetics of transformation but not to the probability of regression. Once LCL are established, large quantities of human cells can be produced. The complete cellular DNA is available so that any part of it can be scrutinized. Moreover, some of the phenotypic characteristics of these B cells be used for a wide range of investigations. The work was supported by Fondation de la Recherche Medicale (France) for Professor J. Dausset (Paris) and by the International Agency for Research on Cancer (Lyon).     

Capacity of robot handling for Epstein-Barr virus transformation

Objectives:  Epstein-Barr virus (EBV) transformation has been described as a routine method to establish human B lymphoblastoid cell lines. Each established lymphoblastoid cell line represents one unique genetic information carrier and can produce unlimited quantities of DNA materials available for downstream applications and research. Undoubtedly, it is of great value to human clinical and experimental genetic studies. However, the current process of EBV transformation requires much manpower in the routine renewal of medium, which is time-consuming. This situation can become a serious problem especially when establishing a human B lymphoblastoid cell bank. A modified and cost-effective protocol for EBV transformation should be considered.
Materials and methods:  In the present study, process in EBV transformation was modified to fit the requirements of robot handling.
Results:  1 mL of whole blood was demonstrated to be sufficient to perform EBV transformation. Additionally, EBV transformation can performed in 96-deep-well plates that are directly and widely used with automatic work platforms.
Conclusions:  Based on these facts, a process of EBV transformation can be modified to fit the requirements of robot handling.     

Isothermal strand-displacement amplification applications for high-throughput genomics

Amplification of source DNA is a nearly universal requirement for molecular biology applications. The primary methods currently available to researchers are limited to in vivo amplification in Escherichia coli hosts and the polymerase chain reaction. Rolling-circle DNA replication is a well-known method for synthesis of phage genomes and recently has been applied as rolling circle amplification (RCA) of specific target sequences as well as circular vectors used in cloning. Here, we demonstrate that RCA using random hexamer primers with 29 DNA polymerase can be used for strand-displacement amplification of different vector constructs containing a variety of insert sizes to produce consistently uniform template for end-sequencing reactions. We show this procedure to be especially effective in a high-throughput plasmid production sequencing process. In addition, we demonstrate that whole bacterial genomes can be effectively amplified from cells or small amounts of purified genomic DNA without apparent bias for use in downstream applications, including whole genome shotgun sequencing.