A genome-wide approach for detecting novel insertion-deletion variants of mid-range size.

We present SWAN, a statistical framework for robust detection of genomic structural variants in next-generation sequencing data and an analysis of mid-range size insertion and deletions (<10 Kb) for whole genome analysis and DNA mixtures. To identify these mid-range size events, SWAN collectively uses information from read-pair, read-depth and one end mapped reads through statistical likelihoods based on Poisson field models. SWAN also uses soft-clip/split read remapping to supplement the likelihood analysis and determine variant boundaries. The accuracy of SWAN is demonstrated by in silico spike-ins and by identification of known variants in the NA12878 genome. We used SWAN to identify a series of novel set of mid-range insertion/deletion detection that were confirmed by targeted deep re-sequencing. An R package implementation of SWAN is open source and freely available.

Haplotyping germline and cancer genomes with high-throughput linked-read sequencing.

Haplotyping of human chromosomes is a prerequisite for cataloguing the full repertoire of genetic variation. We present a microfluidics-based, linked-read sequencing technology that can phase and haplotype germline and cancer genomes using nanograms of input DNA. This high-throughput platform prepares barcoded libraries for short-read sequencing and computationally reconstructs long-range haplotype and structural variant information. We generate haplotype blocks in a nuclear trio that are concordant with expected inheritance patterns and phase a set of structural variants. We also resolve the structure of the EML4-ALK gene fusion in the NCI-H2228 cancer cell line using phased exome sequencing. Finally, we assign genetic aberrations to specific megabase-scale haplotypes generated from whole-genome sequencing of a primary colorectal adenocarcinoma. This approach resolves haplotype information using up to 100 times less genomic DNA than some methods and enables the accurate detection of structural variants.

Sensing of latent EBV infection through exosomal transfer of 5'pppRNA.

Complex interactions between DNA herpesviruses and host factors determine the establishment of a life-long asymptomatic latent infection. The lymphotropic Epstein-Barr virus (EBV) seems to avoid recognition by innate sensors despite massive transcription of immunostimulatory small RNAs (EBV-EBERs). Here we demonstrate that in latently infected B cells, EBER1 transcripts interact with the lupus antigen (La) ribonucleoprotein, avoiding cytoplasmic RNA sensors. However, in coculture experiments we observed that latent-infected cells trigger antiviral immunity in dendritic cells (DCs) through selective release and transfer of RNA via exosomes. In ex vivo tonsillar cultures, we observed that EBER1-loaded exosomes are preferentially captured and internalized by human plasmacytoid DCs (pDCs) that express the TIM1 phosphatidylserine receptor, a known viral- and exosomal target. Using an EBER-deficient EBV strain, enzymatic removal of 5'ppp, in vitro transcripts, and coculture experiments, we established that 5'pppEBER1 transfer via exosomes drives antiviral immunity in nonpermissive DCs. Lupus erythematosus patients suffer from elevated EBV load and activated antiviral immunity, in particular in skin lesions that are infiltrated with pDCs. We detected high levels of EBER1 RNA in such skin lesions, as well as EBV-microRNAs, but no intact EBV-DNA, linking non-cell-autonomous EBER1 presence with skin inflammation in predisposed individuals. Collectively, our studies indicate that virus-modified exosomes have a physiological role in the host-pathogen stand-off and may promote inflammatory disease.

Emergence of Hemagglutinin Mutations During the Course of Influenza Infection.

Influenza remains a significant cause of disease mortality. The ongoing threat of influenza infection is partly attributable to the emergence of new mutations in the influenza genome. Among the influenza viral gene products, the hemagglutinin (HA) glycoprotein plays a critical role in influenza pathogenesis, is the target for vaccines and accumulates new mutations that may alter the efficacy of immunization. To study the emergence of HA mutations during the course of infection, we employed a deep-targeted sequencing method. We used samples from 17 patients with active H1N1 or H3N2 influenza infections. These patients were not treated with antivirals. In addition, we had samples from five patients who were analyzed longitudinally. Thus, we determined the quantitative changes in the fractional representation of HA mutations during the course of infection. Across individuals in the study, a series of novel HA mutations directly altered the HA coding sequence were identified. Serial viral sampling revealed HA mutations that either were stable, expanded or were reduced in representation during the course of the infection. Overall, we demonstrated the emergence of unique mutations specific to an infected individual and temporal genetic variation during infection.

Metastatic tumor evolution and organoid modeling implicate TGFBR2 as a cancer driver in diffuse gastric cancer.

BACKGROUND: Gastric cancer is the second-leading cause of global cancer deaths, with metastatic disease representing the primary cause of mortality. To identify candidate drivers involved in oncogenesis and tumor evolution, we conduct an extensive genome sequencing analysis of metastatic progression in a diffuse gastric cancer. This involves a comparison between a primary tumor from a hereditary diffuse gastric cancer syndrome proband and its recurrence as an ovarian metastasis. RESULTS: Both the primary tumor and ovarian metastasis have common biallelic loss-of-function of both the CDH1 and TP53 tumor suppressors, indicating a common genetic origin. While the primary tumor exhibits amplification of the Fibroblast growth factor receptor 2 (FGFR2) gene, the metastasis notably lacks FGFR2 amplification but rather possesses unique biallelic alterations of Transforming growth factor-beta receptor 2 (TGFBR2), indicating the divergent in vivo evolution of a TGFBR2-mutant metastatic clonal population in this patient. As TGFBR2 mutations have not previously been functionally validated in gastric cancer, we modeled the metastatic potential of TGFBR2 loss in a murine three-dimensional primary gastric organoid culture. The Tgfbr2 shRNA knockdown within Cdh1-/-; Tp53-/- organoids generates invasion in vitro and robust metastatic tumorigenicity in vivo, confirming Tgfbr2 metastasis suppressor activity. CONCLUSIONS: We document the metastatic differentiation and genetic heterogeneity of diffuse gastric cancer and reveal the potential metastatic role of TGFBR2 loss-of-function. In support of this study, we apply a murine primary organoid culture method capable of recapitulating in vivo metastatic gastric cancer. Overall, we describe an integrated approach to identify and functionally validate putative cancer drivers involved in metastasis.

A programmable method for massively parallel targeted sequencing.

We have developed a targeted resequencing approach referred to as Oligonucleotide-Selective Sequencing. In this study, we report a series of significant improvements and novel applications of this method whereby the surface of a sequencing flow cell is modified in situ to capture specific genomic regions of interest from a sample and then sequenced. These improvements include a fully automated targeted sequencing platform through the use of a standard Illumina cBot fluidics station. Targeting optimization increased the yield of total on-target sequencing data 2-fold compared to the previous iteration, while simultaneously increasing the percentage of reads that could be mapped to the human genome. The described assays cover up to 1421 genes with a total coverage of 5.5 Megabases (Mb). We demonstrate a 10-fold abundance uniformity of greater than 90% in 1 log distance from the median and a targeting rate of up to 95%. We also sequenced continuous genomic loci up to 1.5 Mb while simultaneously genotyping SNPs and genes. Variants with low minor allele fraction were sensitively detected at levels of 5%. Finally, we determined the exact breakpoint sequence of cancer rearrangements. Overall, this approach has high performance for selective sequencing of genome targets, configuration flexibility and variant calling accuracy.

Hierarchy in somatic mutations arising during genomic evolution and progression of follicular lymphoma.

Follicular lymphoma (FL) is currently incurable using conventional chemotherapy or immunotherapy regimes, compelling new strategies. Advances in high-throughput sequencing technologies that can reveal oncogenic pathways have stimulated interest in tailoring therapies toward actionable somatic mutations. However, for mutation-directed therapies to be most effective, the mutations must be uniformly present in evolved tumor cells as well as in the self-renewing tumor-cell precursors. Here, we show striking intratumoral clonal diversity within FL tumors in the representation of mutations in the majority of genes as revealed by whole exome sequencing of subpopulations. This diversity captures a clonal hierarchy, resolved using immunoglobulin somatic mutations and IGH-BCL2 translocations as a frame of reference and by comparing diagnosis and relapse tumor pairs, allowing us to distinguish early versus late genetic eventsduring lymphomagenesis. We provide evidence that IGH-BCL2 translocations and CREBBP mutations are early events, whereas MLL2 and TNFRSF14 mutations probably represent late events during disease evolution. These observations provide insight into which of the genetic lesions represent suitable candidates for targeted therapies.

LMP1 association with CD63 in endosomes and secretion via exosomes limits constitutive NF-κB activation.

The ubiquitous Epstein Barr virus (EBV) exploits human B-cell development to establish a persistent infection in ∼90% of the world population. Constitutive activation of NF-κB by the viral oncogene latent membrane protein 1 (LMP1) has an important role in persistence, but is a risk factor for EBV-associated lymphomas. Here, we demonstrate that endogenous LMP1 escapes degradation upon accumulation within intraluminal vesicles of multivesicular endosomes and secretion via exosomes. LMP1 associates and traffics with the intracellular tetraspanin CD63 into vesicles that lack MHC II and sustain low cholesterol levels, even in 'cholesterol-trapping' conditions. The lipid-raft anchoring sequence FWLY, nor ubiquitylation of the N-terminus, controls LMP1 sorting into exosomes. Rather, C-terminal modifications that retain LMP1 in Golgi compartments preclude assembly within CD63-enriched domains and/or exosomal discharge leading to NF-κB overstimulation. Interference through shRNAs further proved the antagonizing role of CD63 in LMP1-mediated signalling. Thus, LMP1 exploits CD63-enriched microdomains to restrain downstream NF-κB activation by promoting trafficking in the endosomal-exosomal pathway. CD63 is thus a critical mediator of LMP1 function in- and outside-infected (tumour) cells.

Exosomes: Fit to deliver small RNA.

Exosomes are specialized membranous nano-sized vesicles derived from endocytic compartments that are released by many cell types. Microvesicles are distinctive from exosomes in that they are produced by shedding of the plasmamembrane and usually larger in size (>1 µm). Exosome biogenesis involves the tightly controlled process of inward budding from the limiting membrane of multivesicular bodies (MVBs). This results in numerous intraluminal vesicles in the lumen of MVBs that contain distinct protein repertoires. It has been suggested that microvesicles shed by certain tumor cells hold functional messenger RNA (mRNA) that may promote tumor progression. We discovered that purified exosomes contain functional microRNAs (miRNAs) and small RNA, but detected little mRNA. Although a clear and decisive distinction between microvesicles and exosomes cannot be made and different subsets of exosomes exist, we speculate that exosomes are specialized in carrying small RNA including the class 22-25 nucleotide regulatory miRNAs. To demonstrate this we developed a co-culture system and found that exosomes are continuously secreted and transferred from Epstein Barr virus (EBV)-infected cells to uninfected neighboring cells. Throughout exosome transfer, the exogenous EBV-encoded miRNAs were delivered to subcellular sites of miRNA-mediated gene repression. Additionally, we found evidence that mature miRNAs are transferred between circulating cells in humans, since we detected EBV-miRNAs in non-infected cells in the peripheral blood of patients that include monocytes and T cells. In this addendum we discuss these findings in the context of recently published papers that advanced our current knowledge of exosome physiology, (mi)RNA function and intercellular RNA transfer. Based on this information we propose that an intercellular (miRNA-based) mode of signal transmission may be well suited in controlling space-confined processes such as the initiation of immune responses in the secondary (peripheral) lymphoid tissues or in a tumor microenvironment. Deciphering the molecular mechanism(s) that control small RNA loading into exosomes and transfer to recipient cells in vitro will provide new evidence for the physiological relevance of vesicle-mediated intercellular communication in vivo.

Functional delivery of viral miRNAs via exosomes.

Noncoding regulatory microRNAs (miRNAs) of cellular and viral origin control gene expression by repressing the translation of mRNAs into protein. Interestingly, miRNAs are secreted actively through small vesicles called "exosomes" that protect them from degradation by RNases, suggesting that these miRNAs may function outside the cell in which they were produced. Here we demonstrate that miRNAs secreted by EBV-infected cells are transferred to and act in uninfected recipient cells. Using a quantitative RT-PCR approach, we demonstrate that mature EBV-encoded miRNAs are secreted by EBV-infected B cells through exosomes. These EBV-miRNAs are functional because internalization of exosomes by MoDC results in a dose-dependent, miRNA-mediated repression of confirmed EBV target genes, including CXCL11/ITAC, an immunoregulatory gene down-regulated in primary EBV-associated lymphomas. We demonstrate that throughout coculture of EBV-infected B cells EBV-miRNAs accumulate in noninfected neighboring MoDC and show that this accumulation is mediated by transfer of exosomes. Thus, the exogenous EBV-miRNAs transferred through exosomes are delivered to subcellular sites of gene repression in recipient cells. Finally, we show in peripheral blood mononuclear cells from patients with increased EBV load that, although EBV DNA is restricted to the circulating B-cell population, EBV BART miRNAs are present in both B-cell and non-B-cell fractions, suggestive of miRNA transfer. Taken together our findings are consistent with miRNA-mediated gene silencing as a potential mechanism of intercellular communication between cells of the immune system that may be exploited by the persistent human gamma-herpesvirus EBV.