Jean Chang

A high-resolution map of human evolutionary constraint using 29 mammals

The comparison of related genomes has emerged as a powerful lens for genome interpretation. Here we report the sequencing and comparative analysis of 29 eutherian genomes. We confirm that at least 5.5% of the human genome has undergone purifying selection, and locate constrained elements covering ∼4.2% of the genome. We use evolutionary signatures and comparisons with experimental data sets to suggest candidate functions for ∼60% of constrained bases. These elements reveal a small number of new coding exons, candidate stop codon readthrough events and over 10,000 regions of overlapping synonymous constraint within protein-coding exons. We find 220 candidate RNA structural families, and nearly a million elements overlapping potential promoter, enhancer and insulator regions. We report specific amino acid residues that have undergone positive selection, 280,000 non-coding elements exapted from mobile elements and more than 1,000 primate- and human-accelerated elements. Overlap with disease-associated variants indicates that our findings will be relevant for studies of human biology, health and disease.

Comparison of Francisella tularensis genomes reveals evolutionary events associated with the emergence of human pathogenic strains

BackgroundFrancisella tularensis subspecies tularensis and holarctica are pathogenic to humans, whereas the two other subspecies, novicida and mediasiatica, rarely cause disease. To uncover the factors that allow subspecies tularensis and holarctica to be pathogenic to humans, we compared their genome sequences with the genome sequence of Francisella tularensis subspecies novicida U112, which is nonpathogenic to humans.ResultsComparison of the genomes of human pathogenic Francisella strains with the genome of U112 identifies genes specific to the human pathogenic strains and reveals pseudogenes that previously were unidentified. In addition, this analysis provides a coarse chronology of the evolutionary events that took place during the emergence of the human pathogenic strains. Genomic rearrangements at the level of insertion sequences (IS elements), point mutations, and small indels took place in the human pathogenic strains during and after differentiation from the nonpathogenic strain, resulting in gene inactivation.ConclusionThe chronology of events suggests a substantial role for genetic drift in the formation of pseudogenes in Francisella genomes. Mutations that occurred early in the evolution, however, might have been fixed in the population either because of evolutionary bottlenecks or because they were pathoadaptive (beneficial in the context of infection). Because the structure of Francisella genomes is similar to that of the genomes of other emerging or highly pathogenic bacteria, this evolutionary scenario may be shared by pathogens from other species.

Methylobacterium Genome Sequences: A Reference Blueprint to Investigate Microbial Metabolism of C1 Compounds from Natural and Industrial Sources

Background Methylotrophy describes the ability of organisms to grow on reduced organic compounds without carbon-carbon bonds. The genomes of two pink-pigmented facultative methylotrophic bacteria of the Alpha-proteobacterial genus Methylobacterium, the reference species Methylobacterium extorquens strain AM1 and the dichloromethane-degrading strain DM4, were compared. Methodology/Principal Findings The 6.88 Mb genome of strain AM1 comprises a 5.51 Mb chromosome, a 1.26 Mb megaplasmid and three plasmids, while the 6.12 Mb genome of strain DM4 features a 5.94 Mb chromosome and two plasmids. The chromosomes are highly syntenic and share a large majority of genes, while plasmids are mostly strain-specific, with the exception of a 130 kb region of the strain AM1 megaplasmid which is syntenic to a chromosomal region of strain DM4. Both genomes contain large sets of insertion elements, many of them strain-specific, suggesting an important potential for genomic plasticity. Most of the genomic determinants associated with methylotrophy are nearly identical, with two exceptions that illustrate the metabolic and genomic versatility of Methylobacterium. A 126 kb dichloromethane utilization (dcm) gene cluster is essential for the ability of strain DM4 to use DCM as the sole carbon and energy source for growth and is unique to strain DM4. The methylamine utilization (mau) gene cluster is only found in strain AM1, indicating that strain DM4 employs an alternative system for growth with methylamine. The dcm and mau clusters represent two of the chromosomal genomic islands (AM1: 28; DM4: 17) that were defined. The mau cluster is flanked by mobile elements, but the dcm cluster disrupts a gene annotated as chelatase and for which we propose the name “island integration determinant” (iid). Conclusion/Significance These two genome sequences provide a platform for intra- and interspecies genomic comparisons in the genus Methylobacterium, and for investigations of the adaptive mechanisms which allow bacterial lineages to acquire methylotrophic lifestyles.

Middle East respiratory syndrome coronavirus (MERS-CoV) causes transient lower respiratory tract infection in rhesus macaques

Significance The Middle East respiratory syndrome coronavirus (MERS-CoV) is the latest emerged coronavirus causing severe respiratory disease with a high case fatality rate in humans. To better understand the disease caused by MERS-CoV, we developed a rhesus macaque model. Infection of rhesus macaques with MERS-CoV resulted in the rapid development of a transient pneumonia, with MERS-CoV replication largely restricted to the lower respiratory tract. This affinity of MERS-CoV for the lungs partly explains the severity of the disease observed in humans. The MERS-CoV rhesus macaque model will be instrumental in developing and testing vaccine and treatment options for an emerging viral pathogen with pandemic potential. In 2012, a novel betacoronavirus, designated Middle East respiratory syndrome coronavirus or MERS-CoV and associated with severe respiratory disease in humans, emerged in the Arabian Peninsula. To date, 108 human cases have been reported, including cases of human-to-human transmission. The availability of an animal disease model is essential for understanding pathogenesis and developing effective countermeasures. Upon a combination of intratracheal, ocular, oral, and intranasal inoculation with 7 × 106 50% tissue culture infectious dose of the MERS-CoV isolate HCoV-EMC/2012, rhesus macaques developed a transient lower respiratory tract infection. Clinical signs, virus shedding, virus replication in respiratory tissues, gene expression, and cytokine and chemokine profiles peaked early in infection and decreased over time. MERS-CoV caused a multifocal, mild to marked interstitial pneumonia, with virus replication occurring mainly in alveolar pneumocytes. This tropism of MERS-CoV for the lower respiratory tract may explain the severity of the disease observed in humans and the, up to now, limited human-to-human transmission.

Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer.

BackgroundPhylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction.ResultsBayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallaces Line, a biogeographic boundary between Southeast Asia and Australia.ConclusionWe describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

Evolution of Burkholderia pseudomallei in Recurrent Melioidosis

Burkholderia pseudomallei, the etiologic agent of human melioidosis, is capable of causing severe acute infection with overwhelming septicemia leading to death. A high rate of recurrent disease occurs in adult patients, most often due to recrudescence of the initial infecting strain. Pathogen persistence and evolution during such relapsing infections are not well understood. Bacterial cells present in the primary inoculum and in late infections may differ greatly, as has been observed in chronic disease, or they may be genetically similar. To test these alternative models, we conducted whole-genome comparisons of clonal primary and relapse B. pseudomallei isolates recovered six months to six years apart from four adult Thai patients. We found differences within each of the four pairs, and some, including a 330 Kb deletion, affected substantial portions of the genome. Many of the changes were associated with increased antibiotic resistance. We also found evidence of positive selection for deleterious mutations in a TetR family transcriptional regulator from a set of 107 additional B. pseudomallei strains. As part of the study, we sequenced to base-pair accuracy the genome of B. pseudomallei strain 1026b, the model used for genetic studies of B. pseudomallei pathogenesis and antibiotic resistance. Our findings provide new insights into pathogen evolution during long-term infections and have important implications for the development of intervention strategies to combat recurrent melioidosis.

Potential Source of Francisella tularensis Live Vaccine Strain Attenuation Determined by Genome Comparison

ABSTRACT Francisella tularensis is a bacterial pathogen that causes the zoonotic disease tularemia and is important to biodefense. Currently, the only vaccine known to confer protection against tularemia is a specific live vaccine strain (designated LVS) derived from a virulent isolate of Francisella tularensis subsp. holarctica. The origin and source of attenuation of this strain are not known. To assist with the design of a defined live vaccine strain, we sought to determine the genetic basis of the attenuation of LVS. This analysis relied primarily on the comparison between the genome of LVS and Francisella tularensis holarctica strain FSC200, which differ by only 0.08% of their nucleotide sequences. Under the assumption that the attenuation was due to a loss of function(s), only coding regions were examined in this comparison. To complement this analysis, the coding regions of two slightly more distantly related Francisella tularensis strains were also compared against the LVS coding regions. Thirty-five genes show unique sequence variations predicted to alter the protein sequence in LVS compared to the other Francisella tularensis strains. Due to these polymorphisms, the functions of 15 of these genes are very likely lost or impaired. Seven of these genes were demonstrated to be under stronger selective constraints, suggesting that they are the most probable to be the source of LVS attenuation and useful for a newly defined vaccine.

Pathogenic Influenza Viruses and Coronaviruses Utilize Similar and Contrasting Approaches To Control Interferon-Stimulated Gene Responses

ABSTRACT The broad range and diversity of interferon-stimulated genes (ISGs) function to induce an antiviral state within the host, impeding viral pathogenesis. While successful respiratory viruses overcome individual ISG effectors, analysis of the global ISG response and subsequent viral antagonism has yet to be examined. Employing models of the human airway, transcriptomics and proteomics datasets were used to compare ISG response patterns following highly pathogenic H5N1 avian influenza (HPAI) A virus, 2009 pandemic H1N1, severe acute respiratory syndrome coronavirus (SARS-CoV), and Middle East respiratory syndrome CoV (MERS-CoV) infection. The results illustrated distinct approaches utilized by each virus to antagonize the global ISG response. In addition, the data revealed that highly virulent HPAI virus and MERS-CoV induce repressive histone modifications, which downregulate expression of ISG subsets. Notably, influenza A virus NS1 appears to play a central role in this histone-mediated downregulation in highly pathogenic influenza strains. Together, the work demonstrates the existence of unique and common viral strategies for controlling the global ISG response and provides a novel avenue for viral antagonism via altered histone modifications. IMPORTANCE This work combines systems biology and experimental validation to identify and confirm strategies used by viruses to control the immune response. Using a novel screening approach, specific comparison between highly pathogenic influenza viruses and coronaviruses revealed similarities and differences in strategies to control the interferon and innate immune response. These findings were subsequently confirmed and explored, revealing both a common pathway of antagonism via type I interferon (IFN) delay as well as a novel avenue for control by altered histone modification. Together, the data highlight how comparative systems biology analysis can be combined with experimental validation to derive novel insights into viral pathogenesis. This work combines systems biology and experimental validation to identify and confirm strategies used by viruses to control the immune response. Using a novel screening approach, specific comparison between highly pathogenic influenza viruses and coronaviruses revealed similarities and differences in strategies to control the interferon and innate immune response. These findings were subsequently confirmed and explored, revealing both a common pathway of antagonism via type I interferon (IFN) delay as well as a novel avenue for control by altered histone modification. Together, the data highlight how comparative systems biology analysis can be combined with experimental validation to derive novel insights into viral pathogenesis.

Chloroplast genome sequencing analysis of Heterosigma akashiwo CCMP452 (West Atlantic) and NIES293 (West Pacific) strains

BackgroundHeterokont algae form a monophyletic group within the stramenopile branch of the tree of life. These organisms display wide morphological diversity, ranging from minute unicells to massive, bladed forms. Surprisingly, chloroplast genome sequences are available only for diatoms, representing two (Coscinodiscophyceae and Bacillariophyceae) of approximately 18 classes of algae that comprise this taxonomic cluster.A universal challenge to chloroplast genome sequencing studies is the retrieval of highly purified DNA in quantities sufficient for analytical processing. To circumvent this problem, we have developed a simplified method for sequencing chloroplast genomes, using fosmids selected from a total cellular DNA library. The technique has been used to sequence chloroplast DNA of two Heterosigma akashiwo strains. This raphidophyte has served as a model system for studies of stramenopile chloroplast biogenesis and evolution.ResultsH. akashiwo strain CCMP452 (West Atlantic) chloroplast DNA is 160,149 bp in size with a 21,822-bp inverted repeat, whereas NIES293 (West Pacific) chloroplast DNA is 159,370 bp in size and has an inverted repeat of 21,665 bp. The fosmid cloning technique reveals that both strains contain an isomeric chloroplast DNA population resulting from an inversion of their single copy domains. Both strains contain multiple small inverted and tandem repeats, non-randomly distributed within the genomes. Although both CCMP452 and NIES293 chloroplast DNAs contains 197 genes, multiple nucleotide polymorphisms are present in both coding and intergenic regions. Several protein-coding genes contain large, in-frame inserts relative to orthologous genes in other plastids. These inserts are maintained in mRNA products. Two genes of interest in H. akashiwo, not previously reported in any chloroplast genome, include tyr C, a tyrosine recombinase, which we hypothesize may be a result of a lateral gene transfer event, and an unidentified 456 amino acid protein, which we hypothesize serves as a G-protein-coupled receptor. The H. akashiwo chloroplast genomes share little synteny with other algal chloroplast genomes sequenced to date.ConclusionThe fosmid cloning technique eliminates chloroplast isolation, does not require chloroplast DNA purification, and reduces sequencing processing time. Application of this method has provided new insights into chloroplast genome architecture, gene content and evolution within the stramenopile cluster.

Mechanisms of Severe Acute Respiratory Syndrome Coronavirus-Induced Acute Lung Injury

ABSTRACT Systems biology offers considerable promise in uncovering novel pathways by which viruses and other microbial pathogens interact with host signaling and expression networks to mediate disease severity. In this study, we have developed an unbiased modeling approach to identify new pathways and network connections mediating acute lung injury, using severe acute respiratory syndrome coronavirus (SARS-CoV) as a model pathogen. We utilized a time course of matched virologic, pathological, and transcriptomic data within a novel methodological framework that can detect pathway enrichment among key highly connected network genes. This unbiased approach produced a high-priority list of 4 genes in one pathway out of over 3,500 genes that were differentially expressed following SARS-CoV infection. With these data, we predicted that the urokinase and other wound repair pathways would regulate lethal versus sublethal disease following SARS-CoV infection in mice. We validated the importance of the urokinase pathway for SARS-CoV disease severity using genetically defined knockout mice, proteomic correlates of pathway activation, and pathological disease severity. The results of these studies demonstrate that a fine balance exists between host coagulation and fibrinolysin pathways regulating pathological disease outcomes, including diffuse alveolar damage and acute lung injury, following infection with highly pathogenic respiratory viruses, such as SARS-CoV. IMPORTANCE Severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in 2002 and 2003, and infected patients developed an atypical pneumonia, acute lung injury (ALI), and acute respiratory distress syndrome (ARDS) leading to pulmonary fibrosis and death. We identified sets of differentially expressed genes that contribute to ALI and ARDS using lethal and sublethal SARS-CoV infection models. Mathematical prioritization of our gene sets identified the urokinase and extracellular matrix remodeling pathways as the most enriched pathways. By infecting Serpine1-knockout mice, we showed that the urokinase pathway had a significant effect on both lung pathology and overall SARS-CoV pathogenesis. These results demonstrate the effective use of unbiased modeling techniques for identification of high-priority host targets that regulate disease outcomes. Similar transcriptional signatures were noted in 1918 and 2009 H1N1 influenza virus-infected mice, suggesting a common, potentially treatable mechanism in development of virus-induced ALI. Severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in 2002 and 2003, and infected patients developed an atypical pneumonia, acute lung injury (ALI), and acute respiratory distress syndrome (ARDS) leading to pulmonary fibrosis and death. We identified sets of differentially expressed genes that contribute to ALI and ARDS using lethal and sublethal SARS-CoV infection models. Mathematical prioritization of our gene sets identified the urokinase and extracellular matrix remodeling pathways as the most enriched pathways. By infecting Serpine1-knockout mice, we showed that the urokinase pathway had a significant effect on both lung pathology and overall SARS-CoV pathogenesis. These results demonstrate the effective use of unbiased modeling techniques for identification of high-priority host targets that regulate disease outcomes. Similar transcriptional signatures were noted in 1918 and 2009 H1N1 influenza virus-infected mice, suggesting a common, potentially treatable mechanism in development of virus-induced ALI.