Clinical Sequencing Uncovers Origins and Evolution of Lassa Virus.

The 2013-2015 West African epidemic of Ebola virus disease (EVD) reminds us of how little is known about biosafety level 4 viruses. Like Ebola virus, Lassa virus (LASV) can cause hemorrhagic fever with high case fatality rates. We generated a genomic catalog of almost 200 LASV sequences from clinical and rodent reservoir samples. We show that whereas the 2013-2015 EVD epidemic is fueled by human-to-human transmissions, LASV infections mainly result from reservoir-to-human infections. We elucidated the spread of LASV across West Africa and show that this migration was accompanied by changes in LASV genome abundance, fatality rates, codon adaptation, and translational efficiency. By investigating intrahost evolution, we found that mutations accumulate in epitopes of viral surface proteins, suggesting selection for immune escape. This catalog will serve as a foundation for the development of vaccines and diagnostics. VIDEO ABSTRACT.

Inflammatory triggers associated with exacerbations of COPD orchestrate plasticity of group 2 innate lymphoid cells in the lungs.

Innate lymphoid cells (ILCs) are critical mediators of mucosal immunity, and group 1 ILCs (ILC1 cells) and group 3 ILCs (ILC3 cells) have been shown to be functionally plastic. Here we found that group 2 ILCs (ILC2 cells) also exhibited phenotypic plasticity in response to infectious or noxious agents, characterized by substantially lower expression of the transcription factor GATA-3 and a concomitant switch to being ILC1 cells that produced interferon-γ (IFN-γ). Interleukin 12 (IL-12) and IL-18 regulated this conversion, and during viral infection, ILC2 cells clustered within inflamed areas and acquired an ILC1-like phenotype. Mechanistically, these ILC1 cells augmented virus-induced inflammation in a manner dependent on the transcription factor T-bet. Notably, IL-12 converted human ILC2 cells into ILC1 cells, and the frequency of ILC1 cells in patients with chronic obstructive pulmonary disease (COPD) correlated with disease severity and susceptibility to exacerbations. Thus, functional plasticity of ILC2 cells exacerbates anti-viral immunity, which may have adverse consequences in respiratory diseases such as COPD.

Cigarette smoke silences innate lymphoid cell function and facilitates an exacerbated type I interleukin-33-dependent response to infection.

Cigarette smoking is a major risk factor for chronic obstructive pulmonary disease and is presumed to be central to the altered responsiveness to recurrent infection in these patients. We examined the effects of smoke priming underlying the exacerbated response to viral infection in mice. Lack of interleukin-33 (IL-33) signaling conferred complete protection during exacerbation and prevented enhanced inflammation and exaggerated weight loss. Mechanistically, smoke was required to upregulate epithelial-derived IL-33 and simultaneously alter the distribution of the IL-33 receptor ST2. Specifically, smoke decreased ST2 expression on group 2 innate lymphoid cells (ILC2s) while elevating ST2 expression on macrophages and natural killer (NK) cells, thus altering IL-33 responsiveness within the lung. Consequently, upon infection and release, increased local IL-33 significantly amplified type I proinflammatory responses via synergistic modulation of macrophage and NK cell function. Therefore, in COPD, smoke alters the lung microenvironment to facilitate an alternative IL-33-dependent exaggerated proinflammatory response to infection, exacerbating disease.

The genomic substrate for adaptive radiation in African cichlid fish.

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Genetic and phenotypic intra-species variation in Candida albicans.

Candida albicans is a commensal fungus of the human gastrointestinal tract and a prevalent opportunistic pathogen. To examine diversity within this species, extensive genomic and phenotypic analyses were performed on 21 clinical C. albicans isolates. Genomic variation was evident in the form of polymorphisms, copy number variations, chromosomal inversions, subtelomeric hypervariation, loss of heterozygosity (LOH), and whole or partial chromosome aneuploidies. All 21 strains were diploid, although karyotypic changes were present in eight of the 21 isolates, with multiple strains being trisomic for Chromosome 4 or Chromosome 7. Aneuploid strains exhibited a general fitness defect relative to euploid strains when grown under replete conditions. All strains were also heterozygous, yet multiple, distinct LOH tracts were present in each isolate. Higher overall levels of genome heterozygosity correlated with faster growth rates, consistent with increased overall fitness. Genes with the highest rates of amino acid substitutions included many cell wall proteins, implicating fast evolving changes in cell adhesion and host interactions. One clinical isolate, P94015, presented several striking properties including a novel cellular phenotype, an inability to filament, drug resistance, and decreased virulence. Several of these properties were shown to be due to a homozygous nonsense mutation in the EFG1 gene. Furthermore, loss of EFG1 function resulted in increased fitness of P94015 in a commensal model of infection. Our analysis therefore reveals intra-species genetic and phenotypic differences in C. albicans and delineates a natural mutation that alters the balance between commensalism and pathogenicity.

In vivo multiplexed interrogation of amplified genes identifies GAB2 as an ovarian cancer oncogene.

High-grade serous ovarian cancers are characterized by widespread recurrent copy number alterations. Although some regions of copy number change harbor known oncogenes and tumor suppressor genes, the genes targeted by the majority of amplified or deleted regions in ovarian cancer remain undefined. Here we systematically tested amplified genes for their ability to promote tumor formation using an in vivo multiplexed transformation assay. We identified the GRB2-associated binding protein 2 (GAB2) as a recurrently amplified gene that potently transforms immortalized ovarian and fallopian tube secretory epithelial cells. Cancer cell lines overexpressing GAB2 require GAB2 for survival and show evidence of phosphatidylinositol 3-kinase (PI3K) pathway activation, which was required for GAB2-induced transformation. Cell lines overexpressing GAB2 were as sensitive to PI3K inhibition as cell lines harboring mutant PIK3CA. Together, these observations nominate GAB2 as an ovarian cancer oncogene, identify an alternative mechanism to activate PI3K signaling, and underscore the importance of PI3K signaling in this cancer.

Comparative analysis of RNA sequencing methods for degraded or low-input samples.

RNA-seq is an effective method for studying the transcriptome, but it can be difficult to apply to scarce or degraded RNA from fixed clinical samples, rare cell populations or cadavers. Recent studies have proposed several methods for RNA-seq of low-quality and/or low-quantity samples, but the relative merits of these methods have not been systematically analyzed. Here we compare five such methods using metrics relevant to transcriptome annotation, transcript discovery and gene expression. Using a single human RNA sample, we constructed and sequenced ten libraries with these methods and compared them against two control libraries. We found that the RNase H method performed best for chemically fragmented, low-quality RNA, and we confirmed this through analysis of actual degraded samples. RNase H can even effectively replace oligo(dT)-based methods for standard RNA-seq. SMART and NuGEN had distinct strengths for measuring low-quantity RNA. Our analysis allows biologists to select the most suitable methods and provides a benchmark for future method development.

The application of next-generation sequencing-based molecular diagnostics in endometrial stromal sarcoma.

AIMS: Endometrial stromal sarcomas (ESSs) are divided into low-grade and high-grade subtypes, with the latter showing more aggressive clinical behaviour. Although histology and immunophenotype can aid in the diagnosis of these tumours, genetic studies can provide additional diagnostic insights, as low-grade ESSs frequently harbour fusions involving JAZF1/SUZ12 and/or JAZF1/PHF1, whereas high-grade ESSs are defined by YWHAE-NUTM2A/B fusions. The aim of this study was to evaluate the utility of a next-generation sequencing (NGS)-based assay in identifying ESS fusions in archival formalin-fixed paraffin-embedded tumour samples. METHODS AND RESULTS: We applied an NGS-based fusion transcript detection assay (Archer FusionPlex Sarcoma Panel) that targets YWHAE and JAZF1 fusions in a series of low-grade ESSs (n = 11) and high-grade ESSs (n = 5) that were previously confirmed to harbour genetic rearrangements by fluorescence in-situ hybridization (FISH) and/or reverse transcription polymerase chain reaction (RT-PCR) analyses. The fusion assay identified junctional fusion transcript sequences that corresponded to the known FISH/RT-PCR results in all cases. Four low-grade ESSs harboured JAZF1-PHF1 fusions with different junctional sequences, and all were correctly identified because of the open-ended nature of the assay design, using anchored multiplex polymerase chain reaction. Seven non-ESS sarcomas were also included as negative controls, and no strong ESS fusion candidates were identified in these cases. CONCLUSIONS: Our findings demonstrate good sensitivity and specificity of an NGS-based gene fusion assay in the detection of ESS fusion transcripts.

Finished bacterial genomes from shotgun sequence data.

Exceptionally accurate genome reference sequences have proven to be of great value to microbial researchers. Thus, to date, about 1800 bacterial genome assemblies have been "finished" at great expense with the aid of manual laboratory and computational processes that typically iterate over a period of months or even years. By applying a new laboratory design and new assembly algorithm to 16 samples, we demonstrate that assemblies exceeding finished quality can be obtained from whole-genome shotgun data and automated computation. Cost and time requirements are thus dramatically reduced.

Paired-end sequencing of Fosmid libraries by Illumina.

Eliminating the bacterial cloning step has been a major factor in the vastly improved efficiency of massively parallel sequencing approaches. However, this also has made it a technical challenge to produce the modern equivalent of the Fosmid- or BAC-end sequences that were crucial for assembling and analyzing complex genomes during the Sanger-based sequencing era. To close this technology gap, we developed Fosill, a method for converting Fosmids to Illumina-compatible jumping libraries. We constructed Fosmid libraries in vectors with Illumina primer sequences and specific nicking sites flanking the cloning site. Our family of pFosill vectors allows multiplex Fosmid cloning of end-tagged genomic fragments without physical size selection and is compatible with standard and multiplex paired-end Illumina sequencing. To excise the bulk of each cloned insert, we introduced two nicks in the vector, translated them into the inserts, and cleaved them. Recircularization of the vector via coligation of insert termini followed by inverse PCR generates a jumping library for paired-end sequencing with 101-base reads. The yield of unique Fosmid-sized jumps is sufficiently high, and the background of short, incorrectly spaced and chimeric artifacts sufficiently low, to enable applications such as mapping of structural variation and scaffolding of de novo assemblies. We demonstrate the power of Fosill to map genome rearrangements in a cancer cell line and identified three fusion genes that were corroborated by RNA-seq data. Our Fosill-powered assembly of the mouse genome has an N50 scaffold length of 17.0 Mb, rivaling the connectivity (16.9 Mb) of the Sanger-sequencing based draft assembly.

The Naked Mole Rat Genome Resource: facilitating analyses of cancer and longevity-related adaptations.

MOTIVATION: The naked mole rat (Heterocephalus glaber) is an exceptionally long-lived and cancer-resistant rodent native to East Africa. Although its genome was previously sequenced, here we report a new assembly sequenced by us with substantially higher N50 values for scaffolds and contigs. RESULTS: We analyzed the annotation of this new improved assembly and identified candidate genomic adaptations which may have contributed to the evolution of the naked mole rat's extraordinary traits, including in regions of p53, and the hyaluronan receptors CD44 and HMMR (RHAMM). Furthermore, we developed a freely available web portal, the Naked Mole Rat Genome Resource (http://www.naked-mole-rat.org), featuring the data and results of our analysis, to assist researchers interested in the genome and genes of the naked mole rat, and also to facilitate further studies on this fascinating species.

Complete viral RNA genome sequencing of ultra-low copy samples by sequence-independent amplification.

RNA viruses are the causative agents for AIDS, influenza, SARS, and other serious health threats. Development of rapid and broadly applicable methods for complete viral genome sequencing is highly desirable to fully understand all aspects of these infectious agents as well as for surveillance of viral pandemic threats and emerging pathogens. However, traditional viral detection methods rely on prior sequence or antigen knowledge. In this study, we describe sequence-independent amplification for samples containing ultra-low amounts of viral RNA coupled with Illumina sequencing and de novo assembly optimized for viral genomes. With 5 million reads, we capture 96 to 100% of the viral protein coding region of HIV, respiratory syncytial and West Nile viral samples from as little as 100 copies of viral RNA. The methods presented here are scalable to large numbers of samples and capable of generating full or near full length viral genomes from clone and clinical samples with low amounts of viral RNA, without prior sequence information and in the presence of substantial host contamination.

Whole genome deep sequencing of HIV-1 reveals the impact of early minor variants upon immune recognition during acute infection.

Deep sequencing technologies have the potential to transform the study of highly variable viral pathogens by providing a rapid and cost-effective approach to sensitively characterize rapidly evolving viral quasispecies. Here, we report on a high-throughput whole HIV-1 genome deep sequencing platform that combines 454 pyrosequencing with novel assembly and variant detection algorithms. In one subject we combined these genetic data with detailed immunological analyses to comprehensively evaluate viral evolution and immune escape during the acute phase of HIV-1 infection. The majority of early, low frequency mutations represented viral adaptation to host CD8+ T cell responses, evidence of strong immune selection pressure occurring during the early decline from peak viremia. CD8+ T cell responses capable of recognizing these low frequency escape variants coincided with the selection and evolution of more effective secondary HLA-anchor escape mutations. Frequent, and in some cases rapid, reversion of transmitted mutations was also observed across the viral genome. When located within restricted CD8 epitopes these low frequency reverting mutations were sufficient to prime de novo responses to these epitopes, again illustrating the capacity of the immune response to recognize and respond to low frequency variants. More importantly, rapid viral escape from the most immunodominant CD8+ T cell responses coincided with plateauing of the initial viral load decline in this subject, suggestive of a potential link between maintenance of effective, dominant CD8 responses and the degree of early viremia reduction. We conclude that the early control of HIV-1 replication by immunodominant CD8+ T cell responses may be substantially influenced by rapid, low frequency viral adaptations not detected by conventional sequencing approaches, which warrants further investigation. These data support the critical need for vaccine-induced CD8+ T cell responses to target more highly constrained regions of the virus in order to ensure the maintenance of immunodominant CD8 responses and the sustained decline of early viremia.

High-quality draft assemblies of mammalian genomes from massively parallel sequence data.

Massively parallel DNA sequencing technologies are revolutionizing genomics by making it possible to generate billions of relatively short (~100-base) sequence reads at very low cost. Whereas such data can be readily used for a wide range of biomedical applications, it has proven difficult to use them to generate high-quality de novo genome assemblies of large, repeat-rich vertebrate genomes. To date, the genome assemblies generated from such data have fallen far short of those obtained with the older (but much more expensive) capillary-based sequencing approach. Here, we report the development of an algorithm for genome assembly, ALLPATHS-LG, and its application to massively parallel DNA sequence data from the human and mouse genomes, generated on the Illumina platform. The resulting draft genome assemblies have good accuracy, short-range contiguity, long-range connectivity, and coverage of the genome. In particular, the base accuracy is high (≥99.95%) and the scaffold sizes (N50 size = 11.5 Mb for human and 7.2 Mb for mouse) approach those obtained with capillary-based sequencing. The combination of improved sequencing technology and improved computational methods should now make it possible to increase dramatically the de novo sequencing of large genomes. The ALLPATHS-LG program is available at http://www.broadinstitute.org/science/programs/genome-biology/crd.

Thymic stromal lymphopoietin is induced by respiratory syncytial virus-infected airway epithelial cells and promotes a type 2 response to infection.

BACKGROUND: Respiratory viral infection, including respiratory syncytial virus (RSV) and rhinovirus, has been linked to respiratory disease in pediatric patients, including severe acute bronchiolitis and asthma exacerbation. OBJECTIVE: The study examined the role of the epithelial-derived cytokine thymic stromal lymphopoietin (TSLP) in the response to RSV infection. METHODS: Infection of human airway epithelial cells was used to examine TSLP induction after RSV infection. Air-liquid interface cultures from healthy children and children with asthma were also tested for TSLP production after infection. Finally, a mouse model was used to directly test the role of TSLP signaling in the response to RSV infection. RESULTS: Infection of airway epithelial cells with RSV led to the production of TSLP via activation of an innate signaling pathway that involved retinoic acid induced gene I, interferon promoter-stimulating factor 1, and nuclear factor-κB. Consistent with this observation, airway epithelial cells from asthmatic children a produced significantly greater levels of TSLP after RSV infection than cells from healthy children. In mouse models, RSV-induced TSLP expression was found to be critical for the development of immunopathology. CONCLUSION: These findings suggest that RSV can use an innate antiviral signaling pathway to drive a potentially nonproductive immune response and has important implications for the role of TSLP in viral immune responses in general.

IL-17-induced pulmonary pathogenesis during respiratory viral infection and exacerbation of allergic disease.

Severe respiratory syncytial virus (RSV) infections are characterized by airway epithelial cell damage, mucus hypersecretion, and Th2 cytokine production. Less is known about the role of IL-17. We observed increased IL-6 and IL-17 levels in tracheal aspirate samples from severely ill infants with RSV infection. In a mouse model of RSV infection, time-dependent increases in pulmonary IL-6, IL-23, and IL-17 expression were observed. Neutralization of IL-17 during infection and observations from IL-17(-/-) knockout mice resulted in significant inhibition of mucus production during RSV infection. RSV-infected animals treated with anti-IL-17 had reduced inflammation and decreased viral load, compared with control antibody-treated mice. Blocking IL-17 during infection resulted in significantly increased RSV-specific CD8 T cells. Factors associated with CD8 cytotoxic T lymphocytes, T-bet, IFN-γ, eomesodermin, and granzyme B were significantly up-regulated after IL-17 blockade. Additionally, in vitro analyses suggest that IL-17 directly inhibits T-bet, eomesodermin, and IFN-γ in CD8 T cells. The role of IL-17 was also investigated in RSV-induced exacerbation of allergic airway responses, in which neutralization of IL-17 led to a significant decrease in the exacerbated disease, including reduced mucus production and Th2 cytokines, with decreased viral proteins. Taken together, our data demonstrate that IL-17 plays a pathogenic role during RSV infections.

Notch ligand delta-like 4 regulates development and pathogenesis of allergic airway responses by modulating IL-2 production and Th2 immunity.

Activation of the canonical Notch pathways has been implicated in Th cell differentiation, but the role of specific Notch ligands in Th2-mediated allergic airway responses has not been completely elucidated. In this study, we show that delta-like ligand 4 (Dll4) was upregulated on dendritic cells in response to cockroach allergen. Blocking Dll4 in vivo during either the primary or secondary response enhanced allergen-induced pathogenic consequences including airway hyperresponsiveness and mucus production via increased Th2 cytokines. In vitro assays demonstrated that Dll4 regulates IL-2 in T cells from established Th2 responses as well as during primary stimulation. Notably, Dll4 blockade during the primary, but not the secondary, response increased IL-2 levels in lung and lymph node of allergic mice. The in vivo neutralization of Dll4 was associated with increased expansion and decreased apoptosis during the primary allergen sensitization. Moreover, Dll4-mediated Notch activation of T cells during primary stimulation in vitro increased apoptosis during the contraction/resting phase of the response, which could be rescued by exogenous IL-2. Consistent with the role for Dll4-mediated IL-2 regulation in overall T cell function, the frequency of IL-4-producing cells was also significantly altered by Dll4 both in vivo and in vitro. These data demonstrate a regulatory role of Dll4 both in initial Th2 differentiation and in Th2 cytokine production in established allergic responses.

Pulmonary IL-17E (IL-25) production and IL-17RB+ myeloid cell-derived Th2 cytokine production are dependent upon stem cell factor-induced responses during chronic allergic pulmonary disease.

In the present studies local neutralization of allergen-induced stem cell factor (SCF) leads to decreased production of Th2 cytokines, a reduction in inflammation, allergen-specific serum IgE/IgG1, and attenuation of severe asthma-like responses. The local blockade of pulmonary SCF also resulted in a significant reduction of IL-17E (IL-25). Sorted cell populations from the lung indicated that IL-25 was produced from c-kit(+) cells, whereas Th2 cytokine production was primarily from c-kit(-) cell populations. SCF stimulated c-kit(+) eosinophils produced IL-25, whereas bone marrow-derived mast cells did not. Using 4get mice that contain a IL-4-IRES-eGFP that when transcribed coexpress GFP and IL-4, our studies identified cells that comprised a CD11b(+), GR1(+), Ly6C(+/-), c-kit(-), CD4(-), CD11c(-), MHC class II(low) cell population as a source of IL-4 in the lung after chronic allergen challenge. In the bone marrow a similar cell was identified with approximately a third of the IL-4(+) cells also expressing c-kit(+). The pulmonary and bone marrow IL-4(+) cell populations were significantly reduced upon local pulmonary anti-SCF treatment. Subsequently, when IL-25R was examined during the chronic allergen responses the expression was found on the IL-4(+) myeloid cell population that expressed CD11b(+)GR1(+). Interestingly, the IL-25R(+) cells in the bone marrow were also all CD11b(+)GR1(+), similar to the lung cells, but they were also all c-kit(+), potentially suggesting a maturation of the bone marrow cell once it enters the lung and/or is stimulated by SCF. Overall, these studies suggest a complex relationship between SCF, bone marrow-derived IL-25-responsive myeloid cells, Th2 cytokines, and chronic allergic disease.