MBEToolbox: a Matlab toolbox for sequence data analysis in molecular biology and evolution

BACKGROUND: MATLAB is a high-performance language for technical computing, integrating computation, visualization, and programming in an easy-to-use environment. It has been widely used in many areas, such as mathematics and computation, algorithm development, data acquisition, modeling, simulation, and scientific and engineering graphics. However, few functions are freely available in MATLAB to perform the sequence data analyses specifically required for molecular biology and evolution. RESULTS: We have developed a MATLAB toolbox, called MBEToolbox, aimed at filling this gap by offering efficient implementations of the most needed functions in molecular biology and evolution. It can be used to manipulate aligned sequences, calculate evolutionary distances, estimate synonymous and nonsynonymous substitution rates, and infer phylogenetic trees. Moreover, it provides an extensible, functional framework for users with more specialized requirements to explore and analyze aligned nucleotide or protein sequences from an evolutionary perspective. The full functions in the toolbox are accessible through the command-line for seasoned MATLAB users. A graphical user interface, that may be especially useful for non-specialist end users, is also provided. CONCLUSION: MBEToolbox is a useful tool that can aid in the exploration, interpretation and visualization of data in molecular biology and evolution. The software is publicly available at and

MBEToolbox 2.0: An enhanced version of a MATLAB toolbox for Molecular Biology and Evolution

MBEToolbox is an extensible MATLAB-based software package for analysis of DNA and protein sequences. MBEToolbox version 2.0 includes enhanced functions for phylogenetic analyses by the maximum likelihood method. For example, it is capable of estimating the synonymous and nonsynonymous substitution rates using a novel or several known codon substitution models. MBEToolbox 2.0 introduces new functions for estimating site-specific evolutionary rates by using a maximum likelihood method or an empirical Bayesian method. It also incorporates several different methods for recombination detection. Multi-platform versions of the software are freely available at http://www.bioinformatics.org/mbetoolbox/

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

Variants of SARS Coronavirus-2 and Their Potential Impact on the Future of the COVID-19 Pandemic

The emergence of SARS-CoV-2 variants of concern (VOCs), especially the sweeping spread of the delta variant, and differing public health management strategies, have rendered global eradication of SARS-CoV-2 unlikely. The currently available COVID-19 vaccines, including the inactivated whole virus vaccines, mRNA vaccines, and adenovirus-vectored vaccines, are effective in protecting people from severe disease and death from COVID-19, but they may not confer good mucosal immunity to prevent the establishment of infection and subsequent viral shedding and transmission. Mucosal vaccines delivered via intranasal route may provide a promising direction, which, if given as a third dose after a two-dose series of intramuscular vaccination, likely promotes mucosal immunity in addition to boosting the systemic cell-mediated immunity and antibody response. However, immunity induced by vaccination, and natural infection as well, is likely to wane followed by re-infection as in the case of human coronaviruses OC43, 229E, NL63, and HKU1. It is a challenge to prevent and control COVID-19 worldwide with the increasing number of VOCs associated with increased transmissibility and changing antigenicity. Nevertheless, we may seek to end the current pandemic situation through mass vaccination and gradual relaxation of non-pharmaceutical measures, which would limit the incidence of severe COVID-19. Repeated doses of booster vaccine will likely be required, similar to influenza virus, especially for the elderly and the immunocompromised patients who are most vulnerable to infection

General metabolism of Laribacter hongkongensis: a genome-wide analysis

<p>Abstract</p> <p>Background</p> <p><it>Laribacter hongkongensis </it>is associated with community-acquired gastroenteritis and traveler's diarrhea. In this study, we performed an in-depth annotation of the genes and pathways of the general metabolism of <it>L. hongkongensis </it>and correlated them with its phenotypic characteristics.</p> <p>Results</p> <p>The <it>L. hongkongensis </it>genome possesses the pentose phosphate and gluconeogenesis pathways and tricarboxylic acid and glyoxylate cycles, but incomplete Embden-Meyerhof-Parnas and Entner-Doudoroff pathways, in agreement with its asaccharolytic phenotype. It contains enzymes for biosynthesis and β-oxidation of saturated fatty acids, biosynthesis of all 20 universal amino acids and selenocysteine, the latter not observed in <it>Neisseria gonorrhoeae</it>, <it>Neisseria meningitidis </it>and <it>Chromobacterium violaceum</it>. The genome contains a variety of dehydrogenases, enabling it to utilize different substrates as electron donors. It encodes three terminal cytochrome oxidases for respiration using oxygen as the electron acceptor under aerobic and microaerophilic conditions and four reductases for respiration with alternative electron acceptors under anaerobic conditions. The presence of complete tetrathionate reductase operon may confer survival advantage in mammalian host in association with diarrhea. The genome contains CDSs for incorporating sulfur and nitrogen by sulfate assimilation, ammonia assimilation and nitrate reduction. The existence of both glutamate dehydrogenase and glutamine synthetase/glutamate synthase pathways suggests an importance of ammonia metabolism in the living environments that it may encounter.</p> <p>Conclusions</p> <p>The <it>L. hongkongensis </it>genome possesses a variety of genes and pathways for carbohydrate, amino acid and lipid metabolism, respiratory chain and sulfur and nitrogen metabolism. These allow the bacterium to utilize various substrates for energy production and survive in different environmental niches.</p

More and More Coronaviruses: Human Coronavirus HKU1

After human coronaviruses OC43, 229E and NL63, human coronavirus HKU1 (HCoV-HKU1) is the fourth human coronavirus discovered. HCoV-HKU1 is a group 2a coronavirus that is still not cultivable. The G + C contents of HCoV-HKU1 genomes are 32%, the lowest among all known coronaviruses with complete genome sequences available. Among all coronaviruses, HCoV-HKU1 shows the most extreme codon usage bias, attributed most importantly to severe cytosine deamination. All HCoV-HKU1 genomes contain unique tandem copies of a 30-base acidic tandem repeat of unknown function at the N-terminus of nsp3 inside the acidic domain upstream of papain-like protease 1. Three genotypes, A, B and C, of HCoV-HKU1 and homologous recombination among their genomes, are observed. The incidence of HCoV-HKU1 infections is the highest in winter. Similar to other human coronaviruses, HCoV-HKU1 infections have been reported globally, with a median (range) incidence of 0.9 (0 – 4.4) %. HCoV-HKU1 is associated with both upper and lower respiratory tract infections that are mostly self-limiting. The most common method for diagnosing HCoV-HKU1 infection is RT-PCR or real-time RT-PCR using RNA extracted from respiratory tract samples such as nasopharyngeal aspirates (NPA). Both the pol and nucleocapsid genes have been used as the targets for amplification. Monoclonal antibodies have been generated for direct antigen detection in NPA. For antibody detection, Escherichia coli BL21 and baculovirus-expressed recombinant nucleocapsid of HCoV-HKU1 have been used for IgG and IgM detection in sera of patients and normal individuals, using Western blot and enzyme-linked immunoassay

Possible Central Nervous System Infection by SARS Coronavirus

On day 22 of illness, generalized tonic-clonic convulsion developed in a 32-year-old woman with severe acute respiratory syndrome (SARS). Cerebrospinal fluid tested positive for SARS coronavirus (SARS-CoV) by reverse transcriptase–polymerase chain reaction. SARS-CoV may have caused an infection in the central nervous system in this patient

Transport genes and chemotaxis in Laribacter hongkongensis: a genome-wide analysis

<p>Abstract</p> <p>Background</p> <p><it>Laribacter hongkongensis </it>is a Gram-negative, sea gull-shaped rod associated with community-acquired gastroenteritis. The bacterium has been found in diverse freshwater environments including fish, frogs and drinking water reservoirs. Using the complete genome sequence data of <it>L. hongkongensis</it>, we performed a comprehensive analysis of putative transport-related genes and genes related to chemotaxis, motility and quorum sensing, which may help the bacterium adapt to the changing environments and combat harmful substances.</p> <p>Results</p> <p>A genome-wide analysis using Transport Classification Database TCDB, similarity and keyword searches revealed the presence of a large diversity of transporters (n = 457) and genes related to chemotaxis (n = 52) and flagellar biosynthesis (n = 40) in the <it>L. hongkongensis </it>genome. The transporters included those from all seven major transporter categories, which may allow the uptake of essential nutrients or ions, and extrusion of metabolic end products and hazardous substances. <it>L. hongkongensis </it>is unique among closely related members of <it>Neisseriaceae </it>family in possessing higher number of proteins related to transport of ammonium, urea and dicarboxylate, which may reflect the importance of nitrogen and dicarboxylate metabolism in this assacharolytic bacterium. Structural modeling of two C^₄-dicarboxylate transporters showed that they possessed similar structures to the determined structures of other DctP-TRAP transporters, with one having an unusual disulfide bond. Diverse mechanisms for iron transport, including hemin transporters for iron acquisition from host proteins, were also identified. In addition to the chemotaxis and flagella-related genes, the <it>L. hongkongensis </it>genome also contained two copies of <it>qseB/qseC </it>homologues of the AI-3 quorum sensing system.</p> <p>Conclusions</p> <p>The large number of diverse transporters and genes involved in chemotaxis, motility and quorum sensing suggested that the bacterium may utilize a complex system to adapt to different environments. Structural modeling will provide useful insights on the transporters in <it>L. hongkongensis</it>.</p