Genomic and phenotypic analyses suggest moderate fitness differences among Zika virus lineages.

RNA viruses have short generation times and high mutation rates, allowing them to undergo rapid molecular evolution during epidemics. However, the extent of RNA virus phenotypic evolution within epidemics and the resulting effects on fitness and virulence remain mostly unknown. Here, we screened the 2015-2016 Zika epidemic in the Americas for lineage-specific fitness differences. We engineered a library of recombinant viruses representing twelve major Zika virus lineages and used them to measure replicative fitness within disease-relevant human primary cells and live mosquitoes. We found that two of these lineages conferred significant in vitro replicative fitness changes among human primary cells, but we did not find fitness changes in Aedes aegypti mosquitoes. Additionally, we found evidence for elevated levels of positive selection among five amino acid sites that define major Zika virus lineages. While our work suggests that Zika virus may have acquired several phenotypic changes during a short time scale, these changes were relatively moderate and do not appear to have enhanced transmission during the epidemic.

Cellular mRNA triggers structural transformation of Ebola virus matrix protein VP40 to its essential regulatory form.

The Ebola virus matrix protein VP40 forms distinct structures linked to distinct functions in the virus life cycle. Dimeric VP40 is a structural protein associated with virus assembly, while octameric, ring-shaped VP40 is associated with transcriptional control. In this study, we show that suitable nucleic acid is sufficient to trigger a dynamic transformation of VP40 dimer into the octameric ring. Deep sequencing reveals a binding preference of the VP40 ring for the 3' untranslated region of cellular mRNA and a guanine- and adenine-rich binding motif. Complementary analyses of the nucleic-acid-induced VP40 ring by native mass spectrometry, electron microscopy, and X-ray crystal structures at 1.8 and 1.4 Å resolution reveal the stoichiometry of RNA binding, as well as an interface involving a key guanine nucleotide. The host factor-induced structural transformation of protein structure in response to specific RNA triggers in the Ebola virus life cycle presents unique opportunities for therapeutic inhibition.

Two Sides of a Coin: a Zika Virus Mutation Selected in Pregnant Rhesus Macaques Promotes Fetal Infection in Mice but at a Cost of Reduced Fitness in Nonpregnant Macaques and Diminished Transmissibility by Vectors.

Although fetal death is now understood to be a severe outcome of congenital Zika syndrome, the role of viral genetics is still unclear. We sequenced Zika virus (ZIKV) from a rhesus macaque fetus that died after inoculation and identified a single intrahost substitution, M1404I, in the ZIKV polyprotein, located in nonstructural protein 2B (NS2B). Targeted sequencing flanking position 1404 in 9 additional macaque mothers and their fetuses identified M1404I at a subconsensus frequency in the majority (5 of 9, 56%) of animals and some of their fetuses. Despite its repeated presence in pregnant macaques, M1404I has occurred rarely in humans since 2015. Since the primary ZIKV transmission cycle is human-mosquito-human, mutations in one host must be retained in the alternate host to be perpetuated. We hypothesized that ZIKV I1404 increases viral fitness in nonpregnant macaques and pregnant mice but is less efficiently transmitted by vectors, explaining its low frequency in humans during outbreaks. By examining competitive fitness relative to that of ZIKV M1404, we observed that ZIKV I1404 produced lower viremias in nonpregnant macaques and was a weaker competitor in tissues. In pregnant wild-type mice, ZIKV I1404 increased the magnitude and rate of placental infection and conferred fetal infection, in contrast to ZIKV M1404, which was not detected in fetuses. Although infection and dissemination rates were not different, Aedes aegypti mosquitoes transmitted ZIKV I1404 more poorly than ZIKV M1404. Our data highlight the complexity of arbovirus mutation-fitness dynamics and suggest that intrahost ZIKV mutations capable of augmenting fitness in pregnant vertebrates may not necessarily spread efficiently via mosquitoes during epidemics.IMPORTANCE Although Zika virus infection of pregnant women can result in congenital Zika syndrome, the factors that cause the syndrome in some but not all infected mothers are still unclear. We identified a mutation that was present in some ZIKV genomes in experimentally inoculated pregnant rhesus macaques and their fetuses. Although we did not find an association between the presence of the mutation and fetal death, we performed additional studies with ZIKV with the mutation in nonpregnant macaques, pregnant mice, and mosquitoes. We observed that the mutation increased the ability of the virus to infect mouse fetuses but decreased its capacity to produce high levels of virus in the blood of nonpregnant macaques and to be transmitted by mosquitoes. This study shows that mutations in mosquito-borne viruses like ZIKV that increase fitness in pregnant vertebrates may not spread in outbreaks when they compromise transmission via mosquitoes and fitness in nonpregnant hosts.

Travel Surveillance and Genomics Uncover a Hidden Zika Outbreak during the Waning Epidemic.

The Zika epidemic in the Americas has challenged surveillance and control. As the epidemic appears to be waning, it is unclear whether transmission is still ongoing, which is exacerbated by discrepancies in reporting. To uncover locations with lingering outbreaks, we investigated travel-associated Zika cases to identify transmission not captured by reporting. We uncovered an unreported outbreak in Cuba during 2017, a year after peak transmission in neighboring islands. By sequencing Zika virus, we show that the establishment of the virus was delayed by a year and that the ensuing outbreak was sparked by long-lived lineages of Zika virus from other Caribbean islands. Our data suggest that, although mosquito control in Cuba may initially have been effective at mitigating Zika virus transmission, such measures need to be maintained to be effective. Our study highlights how Zika virus may still be "silently" spreading and provides a framework for understanding outbreak dynamics. VIDEO ABSTRACT.

The Autoimmune Protocol Diet Modifies Intestinal RNA Expression in Inflammatory Bowel Disease.

Researchers from the Scripps Clinic in La Jolla, CA recently looked at gene expression to better understand the role that diet plays in inflammatory bowel disease. Their findings suggest that diet may help modify inflammatory pathways in people with ulcerative colitis.

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.

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

BACKGROUND: Diagnosis 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. METHODS: We 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. RESULTS: A 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. CONCLUSION: Targeted 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 Whole Blood Molecular Signature for Acute Myocardial Infarction.

Chest pain is a leading reason patients seek medical evaluation. While assays to detect myocyte death are used to diagnose a heart attack (acute myocardial infarction, AMI), there is no biomarker to indicate an impending cardiac event. Transcriptional patterns present in circulating endothelial cells (CEC) may provide a window into the plaque rupture process and identify a proximal biomarker for AMI. Thus, we aimed to identify a transcriptomic signature of AMI present in whole blood, but derived from CECs. Candidate genes indicative of AMI were nominated from microarray of enriched CEC samples, and then verified for detectability and predictive potential via qPCR in whole blood. This signature was validated in an independent cohort. Our findings suggest that a whole blood CEC-derived molecular signature identifies patients with AMI and sets the framework to potentially identify the earlier stages of an impending cardiac event when used in concert with clinical history and other diagnostics where conventional biomarkers indicative of myonecrosis remain undetected.

Efficacy of the Autoimmune Protocol Diet for Inflammatory Bowel Disease.

INTRODUCTION: Data suggest dietary modification can improve clinical responses in inflammatory bowel disease (IBD). The goal of this study was to determine the efficacy of an autoimmune protocol diet in patients with Crohn's disease and ulcerative colitis. METHODS: We enrolled adults with active IBD (Harvey-Bradshaw index ≥ 5 or partial Mayo score ≥3 and erosions on endoscopy and/or elevated fecal calprotectin). For the autoimmune protocol, patients underwent 6-week elimination followed by 5-week maintenance phase. Clinical indices, laboratories, and biomarkers were assessed at baseline and weeks 6 and 11. Endoscopy was performed at study completion. RESULTS: The final cohort included 15 patients with IBD, with mean disease duration 19 years (SD 14.6) and active biological use in 7 (47%) patients. Nutrient repletion was initiated for deficiencies in vitamin D (n = 3) and iron (n = 6). From week 0 to weeks 6 and 11, mean partial Mayo score significantly improved from 5.8 (SD 1.2) to 1.2 (SD 2.0) and 1.0 (SD 2.0) for ulcerative colitis, and mean Harvey-Bradshaw index significantly improved from 7 (SD 1.5) to 3.6 (SD 2.1) and 3.4 (SD 2.6) for Crohn's disease. C-reactive protein did not significantly change during study. Mean fecal calprotectin improved from 471 (SD 562) to 112 (SD 104) at week 11 (P = 0.12). Among those with follow-up endoscopy at week 11 (n = 7), improvements were noted in simple endoscopic score for Crohn's disease (n = 1), Rutgeerts score (n = 1), and Mayo endoscopy subscore (n = 4). DISCUSSION: Dietary elimination can improve symptoms and endoscopic inflammation in patients with IBD. Randomized controlled trials are warranted.

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.

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.

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.

Mutation of WIF1: a potential novel cause of a Nail-Patella-like disorder.

PURPOSE: Nail-Patella syndrome is a dominantly inherited genetic disorder characterized by abnormalities of the nails, knees, elbows, and pelvis. Nail abnormalities are the most constant feature of Nail-Patella syndrome. Pathogenic mutations in a single gene, LMX1B, a mesenchymal determinant of dorsal-ventral patterning, explain approximately 95% of Nail-Patella syndrome cases. However, 5% of cases remain unexplained. METHODS: Here, we present exome sequencing and analysis of four generations of a family with a dominantly inherited Nail-Patella-like disorder (nail dysplasia with some features of Nail-Patella syndrome) who tested negative for LMX1B mutation. RESULTS: We identify a loss-of-function mutation in WIF1 (NM_007191 p.W15*), which is involved in mesoderm segmentation, as the suspected cause of the Nail-Patella-like disorder observed in this family. CONCLUSIONS: Mutation of WIF1 is a potential novel cause of a Nail-Patella-like disorder. Testing of additional patients negative for LMX1B mutation is needed to confirm this finding and further clarify the phenotype.Genet Med advance online publication 06 April 2017.

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.

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. R NA A: nnealing, S: election and L: igation (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.

Gain-of-function ADCY5 mutations in familial dyskinesia with facial myokymia.

OBJECTIVE: 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). METHODS: Whole exome sequencing was performed on 2 parent-child trios. The effect of mutations in ADCY5 was studied by measurement of cyclic adenosine monophosphate (cAMP) accumulation under stimulatory and inhibitory conditions. RESULTS: The same de novo mutation (c.1252C>T, p.R418W) in ADCY5 was found in both studied cases. An inherited missense mutation (c.2176G>A, p.A726T) in ADCY5 was previously reported in a family with FDFM. The significant phenotypic overlap with FDFM was recognized in both cases only after discovery of the molecular link. The inherited mutation in the FDFM family and the recurrent de novo mutation affect residues in different protein domains, the first cytoplasmic domain and the first membrane-spanning domain, respectively. Functional studies revealed a statistically significant increase in ß-receptor agonist-stimulated intracellular cAMP consistent with an increase in adenylyl cyclase activity for both mutants relative to wild-type protein, indicative of a gain-of-function effect. INTERPRETATION: FDFM is likely caused by gain-of-function mutations in different domains of ADCY5-the first definitive link between adenylyl cyclase mutation and human disease. We have illustrated the power of hypothesis-free exome sequencing in establishing diagnoses in rare disorders with complex and variable phenotype. Mutations in ADCY5 should be considered in patients with undiagnosed complex movement disorders even in the absence of a family history.

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.