Comparative Genomic Analysis Reveals a Possible Novel Non-Tuberculous Mycobacterium Species with High Pathogenic Potential

Mycobacteria have been reported to cause a wide range of human diseases. We present the first whole-genome study of a Non-Tuberculous Mycobacterium, Mycobacterium sp. UM_CSW (referred to hereafter as UM_CSW), isolated from a patient diagnosed with bronchiectasis. Our data suggest that this clinical isolate is likely a novel mycobacterial species, supported by clear evidence from molecular phylogenetic, comparative genomic, ANI and AAI analyses. UM_CSW is closely related to the Mycobacterium avium complex. While it has characteristic features of an environmental bacterium, it also shows a high pathogenic potential with the presence of a wide variety of putative genes related to bacterial virulence and shares very similar pathogenomic profiles with the known pathogenic mycobacterial species. Thus, we conclude that this possible novel Mycobacterium species should be tightly monitored for its possible causative role in human infections.     

Comparative Genomic and Phylogenetic Approaches to Characterize the Role of Genetic Recombination in Mycobacterial Evolution

The genus Mycobacterium encompasses over one hundred named species of environmental and pathogenic organisms, including the causative agents of devastating human diseases such as tuberculosis and leprosy. The success of these human pathogens is due in part to their ability to rapidly adapt to their changing environment and host. Recombination is the fastest way for bacterial genomes to acquire genetic material, but conflicting results about the extent of recombination in the genus Mycobacterium have been reported. We examined a data set comprising 18 distinct strains from 13 named species for evidence of recombination. Genomic regions common to all strains (accounting for 10% to 22% of the full genomes of all examined species) were aligned and concatenated in the chromosomal order of one mycobacterial reference species. The concatenated sequence was screened for evidence of recombination using a variety of statistical methods, with each proposed event evaluated by comparing maximum-likelihood phylogenies of the recombinant section with the non-recombinant portion of the dataset. Incongruent phylogenies were identified by comparing the site-wise log-likelihoods of each tree using multiple tests. We also used a phylogenomic approach to identify genes that may have been acquired through horizontal transfer from non-mycobacterial sources. The most frequent associated lineages (and potential gene transfer partners) in the Mycobacterium lineage-restricted gene trees are other members of suborder Corynebacterinae, but more-distant partners were identified as well. In two examined cases of potentially frequent and habitat-directed transfer (M. abscessus to Segniliparus and M. smegmatis to Streptomyces), observed sequence distances were small and consistent with a hypothesis of transfer, while in a third case (M. vanbaalenii to Streptomyces) distances were larger. The analyses described here indicate that whereas evidence of recombination in core regions within the genus is relatively sparse, the acquisition of genes from non-mycobacterial lineages is a significant feature of mycobacterial evolution.     

Genomic,molecular evolution,and expression analysis of NOX genes in soybean (Glycine max)

Reactive oxygen species (ROS) are versatile signaling molecules in sensing stresses and play critical roles in signaling and development. Plasma membrane NADPH oxidases (NOXs) are key producers of ROS, and play important roles in the regulation of plant-pathogen interactions. Here, we performed a comprehensive analysis of the NOX gene family in the soybean genome (Glycine max) and 17 NOX (GmNOX) genes were identified. Structural analysis revealed that the GmNOX proteins in soybean were as conserved as those in other plants. 8 duplicated gene pairs were formed by a Glycine-specific whole-genome duplication (WGD) event approximately 13 million years ago (Mya). The Ka/Ks ratios of GmNOX genes ranged from 0.04 to 0.28, suggesting that the GmNOX family had undergone purifying selection in soybean. Gene expression patterns showed different expression of these duplicate genes, suggesting that the GmNOXs were retained by substantial subfunctionalization during the soybean evolutionary processes. Subsequently, the expression of GmNOXs in response to drought and phytohormones were characterized via qPCR. Importantly, four GmNOXs showed strong expression in nodules, pointing to their probable involvement in nodulation. Thus, our results shed light on the evolutionary history of this family in soybean and contribute to the functional characterization of GmNOX genes in soybean.     

Characterization of a Novel Plasmid,pMAH135, from Mycobacterium avium Subsp. hominissuis

Mycobacterium avium complex (MAC) causes mainly two types of disease. The first is disseminated disease in immunocompromised hosts, such as individuals infected by human immunodeficiency virus (HIV). The second is pulmonary disease in individuals without systemic immunosuppression, and the incidence of this type is increasing worldwide. M. avium subsp. hominissuis, a component of MAC, causes infection in pigs as well as in humans. Many aspects of the different modes of M. avium infection and its host specificity remain unclear. Here, we report the characteristics and complete sequence of a novel plasmid, designated pMAH135, derived from M. avium strain TH135 in an HIV-negative patient with pulmonary MAC disease. The pMAH135 plasmid consists of 194,711 nucleotides with an average G + C content of 66.5% and encodes 164 coding sequences (CDSs). This plasmid was unique in terms of its homology to other mycobacterial plasmids. Interestingly, it contains CDSs with sequence homology to mycobactin biosynthesis proteins and type VII secretion system-related proteins, which are involved in the pathogenicity of mycobacteria. It also contains putative conserved domains of the multidrug efflux transporter. Screening of isolates from humans and pigs for genes located on pMAH135 revealed that the detection rate of these genes was higher in clinical isolates from pulmonary MAC disease patients than in those from HIV-positive patients, whereas the genes were almost entirely absent in isolates from pigs. Moreover, variable number tandem repeats typing analysis showed that isolates carrying pMAH135 genes are grouped in a specific cluster. Collectively, the pMAH135 plasmid contains genes associated with M. avium’s pathogenicity and resistance to antimicrobial agents. The results of this study suggest that pMAH135 influence not only the pathological manifestations of MAC disease, but also the host specificity of MAC infection.     

Predators modify the evolutionary response of prey to temperature change

As climate regimes shift in many ecosystems worldwide, evolution may be a critical process allowing persistence in rapidly changing environments. Organisms regularly interact with other species, yet whether climate-mediated evolution can occur in the context of species interactions is not well understood. We tested whether a species interaction could modify evolutionary responses to temperature. We demonstrate that predation pressure by Dipteran larvae (Chaoborus americanus) modified the evolutionary response of a freshwater crustacean (Daphnia pulex) to its thermal environment over approximately seven generations in laboratory conditions. Daphnia kept at 21°C evolved higher population growth rates than those kept at 18°C, but only in those populations that were also reared with predators. Furthermore, predator-mediated selection resulted in the evolution of elevated Daphnia thermal plasticity. This laboratory natural selection experiment demonstrates that biotic interactions can modify evolutionary adaptation to temperature. Understanding the interplay between multiple selective forces can improve predictions of ecological and evolutionary responses of organisms to rapid environmental change.     

Molecular Evolution of Immune Genes in the Malaria Mosquito Anopheles gambiae

Background

As pathogens that circumvent the host immune response are favoured by selection, so are host alleles that reduce parasite load. Such evolutionary processes leave their signature on the genes involved. Deciphering modes of selection operating on immune genes might reveal the nature of host-pathogen interactions and factors that govern susceptibility in host populations. Such understanding would have important public health implications.

Methodology/Findings

We analyzed polymorphisms in four mosquito immune genes (SP14D1, GNBP, defensin, and gambicin) to decipher selection effects, presumably mediated by pathogens. Using samples of Anopheles arabiensis, An. quadriannulatus and four An. gambiae populations, as well as published sequences from other Culicidae, we contrasted patterns of polymorphisms between different functional units of the same gene within and between populations. Our results revealed selection signatures operating on different time scales. At the most recent time scale, within-population diversity revealed purifying selection. Between populations and between species variation revealed reduced differentiation (GNBP and gambicin) at coding vs. noncoding- regions, consistent with balancing selection. McDonald-Kreitman tests between An. quadriannulatus and both sibling species revealed higher fixation rate of synonymous than nonsynonymous substitutions (GNBP) in accordance with frequency dependent balancing selection. At the longest time scale (>100 my), PAML analysis using distant Culicid taxa revealed positive selection at one codon in gambicin. Patterns of genetic variation were independent of exposure to human pathogens.

Significance and Conclusions

Purifying selection is the most common form of selection operating on immune genes as it was detected on a contemporary time scale on all genes. Selection for “hypervariability” was not detected, but negative balancing selection, detected at a recent evolutionary time scale between sibling species may be rather common. Detection of positive selection at the deepest evolutionary time scale suggests that it occurs infrequently, possibly in association with speciation events. Our results provided no evidence to support the hypothesis that selection was mediated by pathogens that are transmitted to humans.     

Whole Genome Analyses of Marine Fish Pathogenic Isolate,Mycobacterium sp. 012931

Mycobacterium is a genus within the order Actinomycetales that comprises of a large number of well-characterized species, several of which includes pathogens known to cause serious disease in human and animal. Here, we report the whole genome sequence of Mycobacterium sp. strain 012931 isolated from the marine fish, yellowtail (Seriola quinqueradiata). Mycobacterium sp. 012931 is a fish pathogen causing serious damage to aquaculture farms in Japan. DNA dot plot analysis showed that Mycobacterium sp. 012931 was more closely related to Mycobacterium marinum when compared across several Mycobacterium species. However, little conservation of the gene order was observed between Mycobacterium sp. 012931 and M. marinum genome. The annotated 5,464 genes of Mycobacterium sp. 012931 was classified into 26 subsystems. The insertion/deletion gene analysis shows Mycobacterium sp. 012931 had 643 unique genes that were not found in the M. marinum strains. In the virulence, disease, and defense subsystem, both insertion and deletion genes of Mycobacterium sp. 012931 were associated with the PPE gene cluster of Mycobacteria. Of seven plcB genes in Mycobacterium sp. 012931, plcB_2 and plcB_3 showed low identities with those of M. marinum strains. Therefore, Mycobacterium sp. 012931 has differences on genetic and virulence from M. marinum and may induce different interaction mechanisms between host and pathogen.     

Comprehensive profiling of functional attributes,virulence potential and evolutionary dynamics in mycobacterial secretomes

Mycobacterium is an interesting genus which not only includes intimidating pathogens, associated with severe devastations globally, but also comprises of non-pathogenic eco-friendly members that detoxify environmental pollutants. Secretory proteins of the mycobacterial communities are essential components which are firmly believed to facilitate proper cross-talk and apt communication with host cellular surroundings and environmental niche. Secretory elements also play vital roles in mycobacterial pathogenesis. In the present endeavor, an extensive profiling of mycobacterial secretomes, considering both pathogenic and non-pathogenic members, has been executed. Thorough analysis on amino acid composition and functional behavior of the mycobacterial secretory proteins has also been performed. In-depth scrutiny of biosynthetic cost of the secretory proteins with respect to the non-secretory ones indicated that the genus Mycobacterium strictly follows the policy of cost-minimization among the sets of imperative secretory proteins. Comprehensive assessment of potential virulence among the key secretory components signified that the pathogenic mycobacterial members possess a larger share of potentially virulent secretory elements in comparison to their non-pathogenic counterparts. Present analysis also revealed contrasted evolutionary features of the secretomes wherein secretory proteins were found to evolve faster than non-secretory proteins in mycobacterial pathogens but not in the concerned non-pathogens. Outcomes of present investigation promise to provide novel insights into the enigma of mycobacterial pathogenesis, bioremediation and adaptation in diverse niche and aid further scientific investigations associated with concerned research area.     

The characterization of complete mitochondrial genome and phylogenetic relationship within Rattus genus (Rodentia: Muridae)

The genus Rattus is one of the main pest genus of rodent. Most species of the genus carry all kinds of pathogenic bacteria to human being. They are traditionally considered to be a least understood group. The complete mitochondrial genome of the White-Footed Indochinese Rat, Rattus nitidus was determined in this study. The characterization of mitochondrial genomes of Rattus genus was also analyzed based on comprehensive comparison. The result of evolutionary patterns of protein-coding genes (PCGs) suggested purifying selection was the predominant evolutionary forces in the mitochondrial genomes of Rattus genus. The NADH dehydrogenase 4 gene (ND4) showed a highly elevated Ka/Ks ratio compared to the other protein-coding genes, which indicated ND4 was most likely under relaxed selection pressure. Phylogenetic analysis provided a well-supported outline of Rattus genus, and revealed two groups in the genus. R. nitidus had a sister relationship with R. norvegicus.     

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein‐coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein‐coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non‐conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non‐conifers. It is more obvious that in non‐conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6‐ to 3.1‐fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2‐ to 3.6‐fold of non‐conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein‐coding genes due to gene function, plant habitat, or newly acquired IRs.     

Phylogenetic and molecular evolution of the ADAM (A Disintegrin And Metalloprotease) gene family from Xenopus tropicalis, to Mus musculus, Rattus norvegicus, and Homo sapiens

ADAM (a disintegrin and metalloprotease) genes have been identified in various tissues and species, and recently associated with several important human diseases such as tumor and asthma. Although various biological processes have been known for the ADAM family in different species including fertilization, neurogenesis, infection and inflammation, little is known about its detailed phylogenetic and molecular evolutionary history. In this study, the ADAMs of Xenopus (Silurana) tropicalis, Mus musculus, Rattus norvegicus, and Homo sapiens were collected and analyzed by using the Bayesian analysis and gene synteny analysis to establish a comprehensive phylogenetic relationship and evolutionary drive of this gene family. It was found that there were more ADAMs in the two rodents than in the amphibian, suggesting an expansion of the ADAM gene family during the early evolution of mammals. All ADAMs from this expansion were retained in both the rodents, but other duplication events occurred subsequently in the two rodents, respectively, leading to the classification of rodent ADAMs as classes I, II and III. Moreover, these duplicated ADAM genes in the rodents were found to be driven by positive selection, which might be the major force to retain them in the genome. Importantly, it was also found that orthologs of ADAM3 and 5 have been lost in humans. These results not only provide valuable information of the evolution of ADAM genes, but may also help in understanding the role of ADAM genes in the pathobiology of relevant diseases.     

The evolutionary history of the genes involved in the biosynthesis of the antioxidant ergothioneine

Ergothioneine (EGT) is a histidine betaine derivative that exhibits antioxidant action in humans. EGT is primarily synthesized by fungal species and a number of bacterial species. A five-gene cluster (egtA, egtB, egtC, egtD & egtE) responsible for EGT production in Mycobacteria smegmatis has recently been identified. The first fungal biosynthetic EGT gene (NcEgt-1) has also been identified in Neurospora crassa. NcEgt-1 contains domains similar to those found in M. smegmatis egtB and egtD. EGT is biomembrane impermeable. Here we inferred the evolutionary history of the EGT cluster in prokaryotes as well as examining the phyletic distribution of Egt-1 in the fungal kingdom. A genomic survey of 2509 prokaryotes showed that the five-gene EGT cluster is only found in the Actinobacteria. Our survey identified more than 400 diverse prokaryotes that contain genetically linked orthologs of egtB and egtD. Phylogenetic analyses of Egt proteins show a complex evolutionary history and multiple incidences of horizontal gene transfer. Our analysis also identified two independent incidences of a fusion event of egtB and egtD in bacterial species. A genomic survey of over 100 fungal genomes shows that Egt-1 is found in all fungal phyla, except species that belong to the Saccharomycotina subphylum. This analysis provides a comprehensive analysis of the distribution of the key genes involved in the synthesis of EGT in prokaryotes and fungi. Our phylogenetic inferences illuminate the complex evolutionary history of the genes involved in EGT synthesis in prokaryotes. The potential to synthesize EGT is a fungal trait except for species belonging to the Saccharomycotina subphylum.     

Distribution of Environmental Mycobacteria in Karonga District, Northern Malawi

The genus Mycobacterium includes many species that are commonly found in the environment (in soil and water or associated with plants and animals), as well as species that are responsible for two major human diseases, tuberculosis (Mycobacterium tuberculosis) and leprosy (Mycobacterium leprae). The distribution of environmental mycobacteria was investigated in the context of a long-term study of leprosy, tuberculosis, Mycobacterium bovis BCG vaccination, and the responses of individuals to various mycobacterial antigens in Karonga District, northern Malawi, where epidemiological studies had indicated previously that people may be exposed to different mycobacterial species in the northern and southern parts of the district. A total of 148 soil samples and 24 water samples were collected from various locations and examined to determine the presence of mycobacteria. The detection method involved semiselective culturing and acid-fast staining, following decontamination of samples to enrich mycobacteria and reduce the numbers of other microorganisms, or PCR with primers specific for the mycobacterial 16S rRNA gene, using DNA extracted directly from soil and water samples. Mycobacteria were detected in the majority of the samples, and subsequent sequence analysis of PCR products amplified directly from soil DNA indicated that most of the products were related to known environmental mycobacteria. For both methods the rates of recovery were consistently higher for dry season samples than for wet season samples. All isolates cultured from soil appeared to be strains of Mycobacterium fortuitum. This study revealed a complex pattern for the environmental mycobacterial flora but identified no clear differences between the northern and southern parts of Karonga District.     

FGF gene family characterization provides insights into its adaptive evolution in Carnivora

Fibroblast growth factors (FGFs) encoded by the FGF gene family can regulate development and physiology in animals. However, their evolutionary characteristics in Carnivora are largely unknown. In this study, we identified 660 sequences of three types of FGF genes from 30 unannotated genomes of Carnivora animals (before 7th May 2020), and the FGF genes from 52 Carnivora species were analyzed through the method of comparative genomics. Phylogenetic and selective pressure analyses were carried out based on the FGF genes of these 52 Carnivora species. The phylogenetic analysis results demonstrated that the FGF gene family was divided into 10 subfamilies and that FGF5 formed one clade rather than belonging to the subfamilies of FGF4 and FGF6. The evolutionary analysis results showed that the FGF genes were prominently subjected to purifying selection and were highly conserved in the process of Carnivora evolution. We also carried out phylogenetic comparative analyses, which indicated that the habitat was one of the factors that shaped the evolution of Carnivora FGF genes. The FGF1 and FGF6 genes were positively selected in the Carnivora animals, and positive selection signals were detected for the FGF19 gene in semiaquatic Carnivora animals. In summary, we clarified the phylogenetic and evolutionary characteristics of Carnivora FGF genes and provided valuable data for future studies on evolutionary characterization of Carnivora animals.     

Characterization of Cultures Enriched from Acidic Polycyclic Aromatic Hydrocarbon-Contaminated Soil for Growth on Pyrene at Low pH

Two polycyclic aromatic hydrocarbon (PAH)-contaminated soils of pH 2 were successfully used as inoculum to enrich cultures growing on phenanthrene and pyrene at different pHs, including pH 3. Selected pyrene-utilizing cultures obtained at pH 3, pH 5, and pH 7 were further characterized. All showed rapid [¹⁴C]pyrene mineralization at pH 3 and pH 5 and grew on pyrene at pH values ranging from 2 to 6. Eubacterial and mycobacterial 16S rRNA gene denaturing gradient gel electrophoresis fingerprinting and sequencing indicated that the cultures were dominated by a single bacterium closely related to Mycobacterium montefiorense, belonging to the slow-growing Mycobacterium sp. In contrast, a culture enriched on pyrene at pH 7 from a slightly alkaline soil sampled at the same site was dominated by Pseudomonas putida and a fast-growing Mycobacterium sp. The M. montefiorense-related species dominating the pyrene-utilizing cultures enriched from the acidic soils was also the dominant Mycobacterium species in the acidic soils. Our data indicate that a slow-growing Mycobacterium species is involved in PAH degradation in that culture and show that bacteria able to degrade high-molecular-weight PAHs at low pH are present in acidic PAH-contaminated soil.     

Mycobacteriophages: From Petri dish to patient

Mycobacteriophages—bacteriophages infecting Mycobacterium hosts—contribute substantially to our understanding of viral diversity and evolution, provide resources for advancing Mycobacterium genetics, are the basis of high-impact science education programs, and show considerable therapeutic potential. Over 10,000 individual mycobacteriophages have been isolated by high school and undergraduate students using the model organism Mycobacterium smegmatis mc²155 and 2,100 have been completely sequenced, giving a high-resolution view of the phages that infect a single common host strain. The phage genomes are revealed to be highly diverse and architecturally mosaic and are replete with genes of unknown function. Mycobacteriophages have provided many widely used tools for Mycobacterium genetics including integration-proficient vectors and recombineering systems, as well as systems for efficient delivery of reporter genes, transposons, and allelic exchange substrates. The genomic insights and engineering tools have facilitated exploration of phages for treatment of Mycobacterium infections, although their full therapeutic potential has yet to be realized.