Coronavirus Genomics and Bioinformatics Analysis

The drastic increase in the number of coronaviruses discovered and coronavirus genomes being sequenced have given us an unprecedented opportunity to perform genomics and bioinformatics analysis on this family of viruses. Coronaviruses possess the largest genomes (26.4 to 31.7 kb) among all known RNA viruses, with G + C contents varying from 32% to 43%. Variable numbers of small ORFs are present between the various conserved genes (ORF1ab, spike, envelope, membrane and nucleocapsid) and downstream to nucleocapsid gene in different coronavirus lineages. Phylogenetically, three genera, Alphacoronavirus, Betacoronavirus and Gammacoronavirus, with Betacoronavirus consisting of subgroups A, B, C and D, exist. A fourth genus, Deltacoronavirus, which includes bulbul coronavirus HKU11, thrush coronavirus HKU12 and munia coronavirus HKU13, is emerging. Molecular clock analysis using various gene loci revealed that the time of most recent common ancestor of human/civet SARS related coronavirus to be 1999–2002, with estimated substitution rate of 4×10−4 to 2×10−2 substitutions per site per year. Recombination in coronaviruses was most notable between different strains of murine hepatitis virus (MHV), between different strains of infectious bronchitis virus, between MHV and bovine coronavirus, between feline coronavirus (FCoV) type I and canine coronavirus generating FCoV type II, and between the three genotypes of human coronavirus HKU1 (HCoV-HKU1). Codon usage bias in coronaviruses were observed, with HCoV-HKU1 showing the most extreme bias, and cytosine deamination and selection of CpG suppressed clones are the two major independent biological forces that shape such codon usage bias in coronaviruses

CoVDB: a comprehensive database for comparative analysis of coronavirus genes and genomes

The recent SARS epidemic has boosted interest in the discovery of novel human and animal coronaviruses. By July 2007, more than 3000 coronavirus sequence records, including 264 complete genomes, are available in GenBank. The number of coronavirus species with complete genomes available has increased from 9 in 2003 to 25 in 2007, of which six, including coronavirus HKU1, bat SARS coronavirus, group 1 bat coronavirus HKU2, groups 2c and 2d coronaviruses, were sequenced by our laboratory. To overcome the problems we encountered in the existing databases during comparative sequence analysis, we built a comprehensive database, CoVDB (http://covdb.microbiology.hku.hk), of annotated coronavirus genes and genomes. CoVDB provides a convenient platform for rapid and accurate batch sequence retrieval, the cornerstone and bottleneck for comparative gene or genome analysis. Sequences can be directly downloaded from the website in FASTA format. CoVDB also provides detailed annotation of all coronavirus sequences using a standardized nomenclature system, and overcomes the problems of duplicated and identical sequences in other databases. For complete genomes, a single representative sequence for each species is available for comparative analysis such as phylogenetic studies. With the annotated sequences in CoVDB, more specific blast search results can be generated for efficient downstream analysis

PacBio But Not Illumina Technology Can Achieve Fast, Accurate and Complete Closure of the High GC, Complex Burkholderia pseudomallei Two-Chromosome Genome

Although PacBio third-generation sequencers have improved the read lengths of genome sequencing which facilitates the assembly of complete genomes, no study has reported success in using PacBio data alone to completely sequence a two-chromosome bacterial genome from a single library in a single run. Previous studies using earlier versions of sequencing chemistries have at most been able to finish bacterial genomes containing only one chromosome with de novo assembly. In this study, we compared the robustness of PacBio RS II, using one SMRT cell and the latest P6-C4 chemistry, with Illumina HiSeq 1500 in sequencing the genome of Burkholderia pseudomallei, a bacterium which contains two large circular chromosomes, very high G+C content of 68–69%, highly repetitive regions and substantial genomic diversity, and represents one of the largest and most complex bacterial genomes sequenced, using a reference genome generated by hybrid assembly using PacBio and Illumina datasets with subsequent manual validation. Results showed that PacBio data with de novo assembly, but not Illumina, was able to completely sequence the B. pseudomallei genome without any gaps or mis-assemblies. The two large contigs of the PacBio assembly aligned unambiguously to the reference genome, sharing >99.9% nucleotide identities. Conversely, Illumina data assembled using three different assemblers resulted in fragmented assemblies (201–366 contigs), sharing only 92.2–100% and 92.0–100% nucleotide identities to chromosomes I and II reference sequences, respectively, with no indication that the B. pseudomallei genome consisted of two chromosomes with four copies of ribosomal operons. Among all assemblies, the PacBio assembly recovered the highest number of core and virulence proteins, and housekeeping genes based on whole-genome multilocus sequence typing (wgMLST). Most notably, assembly solely based on PacBio outperformed even hybrid assembly using both PacBio and Illumina datasets. Hybrid approach generated only 74 contigs, while the PacBio data alone with de novo assembly achieved complete closure of the two-chromosome B. pseudomallei genome without additional costly bench work and further sequencing. PacBio RS II using P6-C4 chemistry is highly robust and cost-effective and should be the platform of choice in sequencing bacterial genomes, particularly for those that are well-known to be difficult-to-sequence

Genetic Diversity of PCR-Positive, Culture-Negative and Culture-Positive Mycobacterium ulcerans Isolated from Buruli Ulcer Patients in Ghana.

Culture of Mycobacterium ulcerans from Buruli ulcer patients has very low sensitivity. Thus confirmation of M. ulcerans infection is primarily based on PCR directed against IS2404. In this study we compare the genotypes obtained by variable number of tandem repeat analysis of DNA from IS2404-PCR positive cultures with that obtained from IS2404 positive, culture-negative tissue. A significantly greater genetic heterogeneity was found among culture-negative samples compared with that found in cultured strains but a single genotype is over-represented in both sample sets. This study provides evidence that both the focal location of bacteria in a lesion as well as differences in the ability to culture a particular genotype may underlie the low sensitivity of culture. Though preliminary, data from this work also suggests that mycobacteria previously associated with fish disease (M. pseudoshottsii) may be pathogenic for humans

Quasi-Periodic Occultation by a Precessing Accretion Disk and Other Variabilities of SMC X-1

We have investigated the variability of the binary X-ray pulsar, SMC X-1, in data from several X-ray observatories. We confirm the ~60-day cyclic variation of the X-ray flux in the long-term monitoring data from the RXTE and CGRO observatories. X-ray light curves and spectra from the ROSAT, Ginga, and ASCA observatories show that the uneclipsed flux varies by as much as a factor of twenty between a high-flux state when 0.71 second pulses are present and a low-flux state when pulses are absent. In contrast, during eclipses when the X-rays consist of radiation scattered from circumsource matter, the fluxes and spectra in the high and low states are approximately the same. These observations prove that the low state of SMC X-1 is not caused by a reduction in the intrinsic luminosity of the source, or a spectral redistribution thereof, but rather by a quasi-periodic blockage of the line of sight, most likely by a precessing tilted accretion disk. In each of two observations in the midst of low states a brief increase in the X-ray flux and reappearance of 0.71 second pulses occurred near orbital phase 0.2. These brief increases result from an opening of the line of sight to the pulsar that may be caused by wobble in the precessing accretion disk. The records of spin up of the neutron star and decay of the binary orbit are extended during 1991-1996 by pulse-timing analysis of ROSAT, ASCA, and RXTE PCA data. The pulse profiles in various energy ranges from 0.1 to >21 keV are well represented as a combination of a pencil beam and a fan beam. Finally, there is a marked difference between the power spectra of random fluctuations in the high-state data from the RXTE PCA below and above 3.4 keV. Deviation from the fitted power law around 0.06 Hz may be QPO.Comment: Accepted to ApJ. 33 pages including 11 figure

Wild Type and Mutant 2009 Pandemic Influenza A (H1N1) Viruses Cause More Severe Disease and Higher Mortality in Pregnant BALB/c Mice

BACKGROUND: Pregnant women infected by the pandemic influenza A (H1N1) 2009 virus had more severe disease and higher mortality but its pathogenesis is still unclear. PRINCIPAL FINDINGS: We showed that higher mortality, more severe pneumonitis, higher pulmonary viral load, lower peripheral blood T lymphocytes and antibody responses, higher levels of proinflammatory cytokines and chemokines, and worse fetal development occurred in pregnant mice than non-pregnant controls infected by either wild type (clinical isolate) or mouse-adapted mutant virus with D222G substitution in hemagglutinin. These disease-associated changes and the lower respiratory tract involvement were worse in pregnant mice challenged by mutant virus. Though human placental origin JEG-3 cell line could be infected and proinflammatory cytokines or chemokines were elevated in amniotic fluid of some mice, no placental or fetal involvement by virus were detected by culture, real-time reverse transcription polymerase chain reaction or histopathological changes. Dual immunofluorescent staining of viral nucleoprotein and type II alveolar cell marker SP-C protein suggested that the majority of infected alveolar epithelial cells were type II pneumocytes. CONCLUSION: The adverse effect of this pandemic virus on maternal and fetal outcome is largely related to the severe pulmonary disease and the indirect effect of inflammatory cytokine spillover into the systemic circulation

Cryptic Polyketide Synthase Genes in Non-Pathogenic Clostridium SPP

Modular type I polyketide synthases (PKS) produce a vast array of bacterial metabolites with highly diverse biological functions. Notably, all known polyketides were isolated from aerobic bacteria, and yet no example has been reported for strict anaerobes. In this study we explored the diversity and distribution of PKS genes in the genus Clostridium. In addition to comparative genomic analyses combined with predictions of modular type I polyketide synthase (PKS) gene clusters in sequenced genomes of Clostridium spp., a representative selection of other species inhabiting a variety of ecological niches was investigated by PCR screening for PKS genes. Our data reveal that all studied pathogenic Clostridium spp. are devoid of putative PKS genes. In stark contrast, cryptic PKS genes are widespread in genomes of non-pathogenic Clostridium species. According to phylogenetic analyses, the Clostridium PKS genes have unusual and diverse origins. However, reverse transcription quantitative PCR demonstrates that these genes are silent under standard cultivation conditions, explaining why the related metabolites have been overlooked until now. This study presents clostridia as a putative source for novel bioactive polyketides

Embryonic Lethality in Mice Lacking the Nuclear Factor of Activated T Cells 5 Protein Due to Impaired Cardiac Development and Function

Nuclear factor of activated T cells 5 protein (NFAT5) is thought to be important for cellular adaptation to osmotic stress by regulating the transcription of genes responsible for the synthesis or transport of organic osmolytes. It is also thought to play a role in immune function, myogenesis and cancer invasion. To better understand the function of NFAT5, we developed NFAT5 gene knockout mice. Homozygous NFAT5 null (NFAT5−/−) mouse embryos failed to develop normally and died after 14.5 days of embryonic development (E14.5). The embryos showed peripheral edema, and abnormal heart development as indicated by thinner ventricular wall and reduced cell density at the compact and trabecular areas of myocardium. This is associated with reduced level of proliferating cell nuclear antigen and increased caspase-3 in these tissues. Cardiomyocytes from E14.5 NFAT5−/− embryos showed a significant reduction of beating rate and abnormal Ca2+ signaling profile as a consequence of reduced sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) and ryanodine receptor (RyR) expressions. Expression of NFAT5 target genes, such as HSP 70 and SMIT were reduced in NFAT5−/− cardiomyocytes. Our findings demonstrated an essential role of NFAT5 in cardiac development and Ca2+ signaling. Cardiac failure is most likely responsible for the peripheral edema and death of NFAT5−/− embryos at E14.5 days

Integrating Communities of Practice in Technology Development Projects

Technology development projects usually benefit when knowledge and expertise are drawn from a variety of sources, including potential users. Orchestrating the involvement of people from disparate groups is a crucial task for project managers. It requires finding a balance between differentiation, when teams work in isolation, and integration, when groups come together to exchange knowledge. This article argues that a “community of practice” perspective can help project managers to achieve this balance, by drawing attention to the assumptions, interests, skills, and formal and tacit knowledge of the different groups involved. Successful integration can be achieved by ensuring that the developing technology is comprehensible to all the groups concerned, and making sure that it satisfies their various interests