Glenn Oliveira

Zika virus evolution and spread in the Americas

Although the recent Zika virus (ZIKV) epidemic in the Americas and its link to birth defects have attracted a great deal of attention, much remains unknown about ZIKV disease epidemiology and ZIKV evolution, in part owing to a lack of genomic data. Here we address this gap in knowledge by using multiple sequencing approaches to generate 110 ZIKV genomes from clinical and mosquito samples from 10 countries and territories, greatly expanding the observed viral genetic diversity from this outbreak. We analysed the timing and patterns of introductions into distinct geographic regions; our phylogenetic evidence suggests rapid expansion of the outbreak in Brazil and multiple introductions of outbreak strains into Puerto Rico, Honduras, Colombia, other Caribbean islands, and the continental United States. We find that ZIKV circulated undetected in multiple regions for many months before the first locally transmitted cases were confirmed, highlighting the importance of surveillance of viral infections. We identify mutations with possible functional implications for ZIKV biology and pathogenesis, as well as those that might be relevant to the effectiveness of diagnostic tests.

Gain‐of‐function ADCY5 mutations in familial dyskinesia with facial myokymia

To identify the cause of childhood onset involuntary paroxysmal choreiform and dystonic movements in 2 unrelated sporadic cases and to investigate the functional effect of missense mutations in adenylyl cyclase 5 (ADCY5) in sporadic and inherited cases of autosomal dominant familial dyskinesia with facial myokymia (FDFM).

Genomic epidemiology reveals multiple introductions of Zika virus into the United States

Zika virus (ZIKV) is causing an unprecedented epidemic linked to severe congenital abnormalities. In July 2016, mosquito-borne ZIKV transmission was reported in the continental United States; since then, hundreds of locally acquired infections have been reported in Florida. To gain insights into the timing, source, and likely route(s) of ZIKV introduction, we tracked the virus from its first detection in Florida by sequencing ZIKV genomes from infected patients and Aedes aegypti mosquitoes. We show that at least 4 introductions, but potentially as many as 40, contributed to the outbreak in Florida and that local transmission is likely to have started in the spring of 2016—several months before its initial detection. By analysing surveillance and genetic data, we show that ZIKV moved among transmission zones in Miami. Our analyses show that most introductions were linked to the Caribbean, a finding corroborated by the high incidence rates and traffic volumes from the region into the Miami area. Our study provides an understanding of how ZIKV initiates transmission in new regions.

Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples

Genome sequencing has become a powerful tool for studying emerging infectious diseases; however, genome sequencing directly from clinical samples (i.e., without isolation and culture) remains challenging for viruses such as Zika, for which metagenomic sequencing methods may generate insufficient numbers of viral reads. Here we present a protocol for generating coding-sequence-complete genomes, comprising an online primer design tool, a novel multiplex PCR enrichment protocol, optimized library preparation methods for the portable MinION sequencer (Oxford Nanopore Technologies) and the Illumina range of instruments, and a bioinformatics pipeline for generating consensus sequences. The MinION protocol does not require an Internet connection for analysis, making it suitable for field applications with limited connectivity. Our method relies on multiplex PCR for targeted enrichment of viral genomes from samples containing as few as 50 genome copies per reaction. Viral consensus sequences can be achieved in 1–2 d by starting with clinical samples and following a simple laboratory workflow. This method has been successfully used by several groups studying Zika virus evolution and is facilitating an understanding of the spread of the virus in the Americas. The protocol can be used to sequence other viral genomes using the online Primal Scheme primer designer software. It is suitable for sequencing either RNA or DNA viruses in the field during outbreaks or as an inexpensive, convenient method for use in the lab.

Sensitive, multiplex and direct quantification of RNA sequences using a modified RASL assay

A sensitive and highly multiplex method to directly measure RNA sequence abundance without requiring reverse transcription would be of value for a number of biomedical applications, including high throughput small molecule screening, pathogen transcript detection and quantification of short/degraded RNAs. RNA Annealing, Selection and Ligation (RASL) assays, which are based on RNA template-dependent oligonucleotide probe ligation, have been developed to meet this need, but technical limitations have impeded their adoption. Whereas DNA ligase-based RASL assays suffer from extremely low and sequence-dependent ligation efficiencies that compromise assay robustness, Rnl2 can join a fully DNA donor probe to a 3′-diribonucleotide-terminated acceptor probe with high efficiency on an RNA template strand. Rnl2-based RASL exhibits sub-femtomolar transcript detection sensitivity, and permits the rational tuning of probe signals for optimal analysis by massively parallel DNA sequencing (RASL-seq). A streamlined Rnl2-based RASL-seq protocol was assessed in a small molecule screen using 77 probe sets designed to monitor complex human B cell phenotypes during antibody class switch recombination. Our data demonstrate the robustness, cost-efficiency and broad applicability of Rnl2-based RASL assays.

Rapid deep sequencing of patient-derived HIV with ion semiconductor technology

The development of next-generation sequencing technologies has facilitated the study of HIV drug resistance evolution. However, the high capacity and per-run cost of many sequencers is not ideal for viral sequencing unless many samples are analyzed simultaneously. Ion semiconductor sequencing has recently emerged as a flexible, lower-cost alternative with short runtime. This paper describes the use of Ion Torrent devices for deep sequencing of drug resistant HIV samples. High levels of sequencing coverage were obtained in HIV Gag and protease, allowing the detection of mutations at low frequencies.

Deep sequencing of protease inhibitor resistant HIV patient isolates reveals patterns of correlated mutations in Gag and protease.

While the role of drug resistance mutations in HIV protease has been studied comprehensively, mutations in its substrate, Gag, have not been extensively cataloged. Using deep sequencing, we analyzed a unique collection of longitudinal viral samples from 93 patients who have been treated with therapies containing protease inhibitors (PIs). Due to the high sequence coverage within each sample, the frequencies of mutations at individual positions were calculated with high precision. We used this information to characterize the variability in the Gag polyprotein and its effects on PI-therapy outcomes. To examine covariation of mutations between two different sites using deep sequencing data, we developed an approach to estimate the tight bounds on the two-site bivariate probabilities in each viral sample, and the mutual information between pairs of positions based on all the bounds. Utilizing the new methodology we found that mutations in the matrix and p6 proteins contribute to continued therapy failure and have a major role in the network of strongly correlated mutations in the Gag polyprotein, as well as between Gag and protease. Although covariation is not direct evidence of structural propensities, we found the strongest correlations between residues on capsid and matrix of the same Gag protein were often due to structural proximity. This suggests that some of the strongest inter-protein Gag correlations are the result of structural proximity. Moreover, the strong covariation between residues in matrix and capsid at the N-terminus with p1 and p6 at the C-terminus is consistent with residue-residue contacts between these proteins at some point in the viral life cycle.

ASXL1 and DNMT3A mutation in a cytogenetically normal B3 thymoma

The molecular drivers of thymoma are poorly understood. Outside of the identification of rarely occurring epidermal growth factor receptor and v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog mutations via candidate gene sequencing, mutations in common cancer genes have yet to be observed. Only a single thymoma genome sequence has been previously reported, with no mutations in known cancer genes identified. Thus, we attempted to identify somatic driver mutations in a cytogenetically normal thymoma. A stage IVB type B3 thymoma from a 47-year-old male of Asian descent with no history of myasthenia gravis or other autoimmune condition was genomically evaluated. Exome sequencing and low-pass whole-genome sequencing was performed to identify somatic point mutations, copy number changes and structural variants. Mutations in known tumor suppressors DNMT3A (p.G728D) and ASXL1 (p.E657fs), consistent with mutations of known consequence in acute myeloid leukemia, were identified. Contrary to a previous report, this finding suggests the genetic etiology of thymomas may not be fundamentally distinct from other tumor types. Rather, these findings suggest that further sequencing of cytogenetically normal thymoma samples should reveal the specific molecular drivers of thymoma.

Spectrum of mutations in monogenic diabetes genes identified from high-throughput DNA sequencing of 6888 individuals

BackgroundDiagnosis of monogenic as well as atypical forms of diabetes mellitus has important clinical implications for their specific diagnosis, prognosis, and targeted treatment. Single gene mutations that affect beta-cell function represent 1–2% of all cases of diabetes. However, phenotypic heterogeneity and lack of family history of diabetes can limit the diagnosis of monogenic forms of diabetes. Next-generation sequencing technologies provide an excellent opportunity to screen large numbers of individuals with a diagnosis of diabetes for mutations in disease-associated genes.MethodsWe utilized a targeted sequencing approach using the Illumina HiSeq to perform a case-control sequencing study of 22 monogenic diabetes genes in 4016 individuals with type 2 diabetes (including 1346 individuals diagnosed before the age of 40 years) and 2872 controls. We analyzed protein-coding variants identified from the sequence data and compared the frequencies of pathogenic variants (protein-truncating variants and missense variants) between the cases and controls.ResultsA total of 40 individuals with diabetes (1.8% of early onset sub-group and 0.6% of adult onset sub-group) were carriers of known pathogenic missense variants in the GCK, HNF1A, HNF4A, ABCC8, and INS genes. In addition, heterozygous protein truncating mutations were detected in the GCK, HNF1A, and HNF1B genes in seven individuals with diabetes. Rare missense mutations in the GCK gene were significantly over-represented in individuals with diabetes (0.5% carrier frequency) compared to controls (0.035%). One individual with early onset diabetes was homozygous for a rare pathogenic missense variant in the WFS1 gene but did not have the additional phenotypes associated with Wolfram syndrome.ConclusionTargeted sequencing of genes linked with monogenic diabetes can identify disease-relevant mutations in individuals diagnosed with type 2 diabetes not suspected of having monogenic forms of the disease. Our data suggests that GCK-MODY frequently masquerades as classical type 2 diabetes. The results confirm that MODY is under-diagnosed, particularly in individuals presenting with early onset diabetes and clinically labeled as type 2 diabetes; thus, sequencing of all monogenic diabetes genes should be routinely considered in such individuals. Genetic information can provide a specific diagnosis, inform disease prognosis and may help to better stratify treatment plans.

A feasibility study of colorectal cancer diagnosis via circulating tumor DNA derived CNV detection

Circulating tumor DNA (ctDNA) has shown great promise as a biomarker for early detection of cancer. However, due to the low abundance of ctDNA, especially at early stages, it is hard to detect at high accuracies while keeping sequencing costs low. Here we present a pilot stage study to detect large scale somatic copy numbers variations (CNVs), which contribute more molecules to ctDNA signal compared to point mutations, via cell free DNA sequencing. We show that it is possible to detect somatic CNVs in early stage colorectal cancer (CRC) patients and subsequently discriminate them from normal patients. With 25 normal and 24 CRC samples, we achieve 100% specificity (lower bound confidence interval: 86%) and ~79% sensitivity (95% confidence interval: 63% - 95%,), though the performance should be considered with caution given the limited sample size. We report a lack of concordance between the CNVs detected via cfDNA sequencing and CNVs identified in parent tissue samples. However, recent findings suggest that a lack of concordance is expected for CNVs in CRC because of their sub-clonal nature. Finally, the CNVs we detect very likely contribute to cancer progression as they lie in functionally important regions, and have been shown to be associated with CRC specifically. This study paves the path for a larger scale exploration of the potential of CNV detection for both diagnoses and prognoses of cancer.