Detecting genetic gain and loss events in terms of protein domain: Method and implementation.

Continuous gain and loss of genes are the primary driving forces of bacterial evolution and environmental adaptation. Studying bacterial evolution in terms of protein domain, which is the fundamental function and evolutionary unit of proteins, can provide a more comprehensive understanding of bacterial differentiation and phenotypic adaptation processes. Therefore, we proposed a phylogenetic tree-based method for detecting genetic gain and loss events in terms of protein domains. Specifically, the method focuses on a single domain to trace its evolution process or on multiple domains to investigate their co-evolution principles. This novel method was validated using 122 Shigella isolates. We found that the loss of a significant number of domains was likely the main driving force behind the evolution of Shigella, which could reduce energy expenditure and preserve only the most essential functions. Additionally, we observed that simultaneously gained and lost domains were often functionally related, which can facilitate and accelerate phenotypic evolutionary adaptation to the environment. All results obtained using our method agree with those of previous studies, which validates our proposed method.

Comparative genomic analysis provides insights into the genetic diversity and pathogenicity of the genus Brucella.

Some Brucella spp. are important pathogens. According to the latest prokaryotic taxonomy, the Brucella genus consists of facultative intracellular parasitic Brucella species and extracellular opportunistic or environmental Brucella species. Intracellular Brucella species include classical and nonclassical types, with different species generally exhibiting host preferences. Some classical intracellular Brucella species can cause zoonotic brucellosis, including B. melitensis, B. abortus, B. suis, and B. canis. Extracellular Brucella species comprise opportunistic or environmental species which belonged formerly to the genus Ochrobactrum and thus nowadays renamed as for example Brucella intermedia or Brucella anthropi, which are the most frequent opportunistic human pathogens within the recently expanded genus Brucella. The cause of the diverse phenotypic characteristics of different Brucella species is still unclear. To further investigate the genetic evolutionary characteristics of the Brucella genus and elucidate the relationship between its genomic composition and prediction of phenotypic traits, we collected the genomic data of Brucella from the NCBI Genome database and conducted a comparative genomics study. We found that classical and nonclassical intracellular Brucella species and extracellular Brucella species exhibited differences in phylogenetic relationships, horizontal gene transfer and distribution patterns of mobile genetic elements, virulence factor genes, and antibiotic resistance genes, showing the close relationship between the genetic variations and prediction of phenotypic traits of different Brucella species. Furthermore, we found significant differences in horizontal gene transfer and the distribution patterns of mobile genetic elements, virulence factor genes, and antibiotic resistance genes between the two chromosomes of Brucella, indicating that the two chromosomes had distinct dynamics and plasticity and played different roles in the survival and evolution of Brucella. These findings provide new directions for exploring the genetic evolutionary characteristics of the Brucella genus and could offer new clues to elucidate the factors influencing the phenotypic diversity of the Brucella genus.

Spatio-temporal spread and evolution of Lassa virus in West Africa.

BACKGROUND: Lassa fever is a hemorrhagic disease caused by Lassa virus (LASV), which has been classified by the World Health Organization as one of the top infectious diseases requiring prioritized research. Previous studies have provided insights into the classification and geographic characteristics of LASV lineages. However, the factor of the distribution and evolution characteristics and phylodynamics of the virus was still limited. METHODS: To enhance comprehensive understanding of LASV, we employed phylogenetic analysis, reassortment and recombination detection, and variation evaluation utilizing publicly available viral genome sequences. RESULTS: The results showed the estimated the root of time of the most recent common ancestor (TMRCA) for large (L) segment was approximately 634 (95% HPD: [385879]), whereas the TMRCA for small (S) segment was around 1224 (95% HPD: [10301401]). LASV primarily spread from east to west in West Africa through two routes, and in route 2, the virus independently spread to surrounding countries through Liberia, resulting in a wider spread of LASV. From 1969 to 2018, the effective population size experienced two significant increased, indicating the enhanced genetic diversity of LASV. We also found the evolution rate of L segment was faster than S segment, further results showed zinc-binding protein had the fastest evolution rate. Reassortment events were detected in multiple lineages including sub-lineage IIg, while recombination events were observed within lineage V. Significant amino acid changes in the glycoprotein precursor of LASV were identified, demonstrating sequence diversity among lineages in LASV. CONCLUSION: This study comprehensively elucidated the transmission and evolution of LASV in West Africa, providing detailed insights into reassortment events, recombination events, and amino acid variations.

Monkeypox: Can we count on the current smallpox immunization?

Two vaccines ACAM 2000 and JYNNEOS have obtained approval from the Food and Drug Administration as preventive measures against monkeypox, contributing significantly to the management of the monkeypox epidemic. Nonetheless, research has demonstrated that smallpox vaccination offers approximately 88.8% protection against monkeypox, while immunization with these vaccines generates relatively low levels of neutralizing antibodies. In this work, we performed a comprehensive comparison of antigens between the 2022-2023 monkeypox strains and the smallpox vaccine strains. Our analysis has revealed considerable amino acid changes in all 27 antigens, including core and envelope proteins. Amino acid substitutions within B cell epitopes were observed in 26 of these antigens, with at least half of the antigen substitutions occurring within B cell epitopes in 20 out of the 26 antigens analyzed. These findings may raise potential concerns regarding the efficacy of these vaccines.

Evolution of the ability to evade host innate immune defense by Talaromyces marneffei.

Talaromyces (Penicillium) marneffei is an intracellular pathogenic fungus. Some strains of this fungus have been misidentified due to the similarity between Talaromyces and Penicillium. T. marneffei has mainly been found to afflict immunocompromised individuals, causing respiratory, skin, and systemic mycosis. Mp1p is a key virulence factor that can help T. marneffei evade clearance by the normally functioning immune system. Understanding how novel functions arise is an intriguing question in many fields of biology. Mp1p has two homologous domains (Mp1p-LBD1 and Mp1p-LBD2). Sequence similarity searches with Mp1p-LBD sequences revealed Mp1p homologs in many other pathogenic fungi. Integrated information on the taxonomic distribution, phylogenetic relationships, and sequence similarity of Mp1p domains revealed that the ancestor of Mp1p-LBDs was acquired through horizontal gene transfer (HGT). Additional evidence revealed that Mp1p homologs have undergone extensive gene duplications in T. marneffei. Mp1p might be a result of gene fusion following gene duplication. Furthermore, we propose a new method for identifying Talaromyces and identify 4 strains with misclassification errors. Our results characterize the evolutionary mechanism of T. marneffei evasion of host innate immune defense and clearly demonstrate the role of gene duplication and HGT in the evolution of host immune escape by T. marneffei.

Genetic characteristics of classical astroviruses in Shenzhen, China, 2016-2019.

Human astrovirus (HAstV) is a single-stranded, positive-sense RNA virus and is the leading cause of viral gastroenteritis. However, despite its prevalence, astroviruses still remain one of the least studied enteroviruses. In this study, we sequenced 11 classical astrovirus strains from clinical samples collected in Shenzhen, China from 2016 to 2019, analyzed their genetic characteristics, and deposited them into GenBank. We conducted phylogenetic analysis using IQ-TREE software, with references to astrovirus sequences worldwide. The phylogeographic analysis was performed using the Bayesian Evolutionary Analysis Sampling Trees program, through Bayesian Markov Chain Monte Carlo sampling. We also conducted recombination analysis with the Recombination Detection Program. The newly sequenced strains were categorized as HAstV genotype 1, which is the predominant genotype in Shenzhen. Phylogeographic reconstruction indicated that HAstV-1 may have migrated from the United States to China, followed by frequent transmission between China and Japan. The recombination analysis revealed recombination events within and across genotypes, and identified a recombination-prone region that produced relatively uniform recombination breakpoints and fragment lengths. The genetic analysis of HAstV strains in Shenzhen addresses the current lack of astrovirus data in the region of Shenzhen and provides key insights to the evolution and transmission of astroviruses worldwide. These findings highlight the importance of improving surveillance of astroviruses.

Characterization and Diversity of Klebsiella pneumoniae Prophages.

Klebsiella pneumoniae is a common human commensal and opportunistic pathogen. In recent years, the clinical isolation and resistance rates of K. pneumoniae have shown a yearly increase, leading to a special interest in mobile genetic elements. Prophages are a representative class of mobile genetic elements that can carry host-friendly genes, transfer horizontally between strains, and coevolve with the host's genome. In this study, we identified 15,946 prophages from the genomes of 1437 fully assembled K. pneumoniae deposited in the NCBI database, with 9755 prophages on chromosomes and 6191 prophages on plasmids. We found prophages to be notably diverse and widely disseminated in the K. pneumoniae genomes. The K. pneumoniae prophages encoded multiple putative virulence factors and antibiotic resistance genes. The comparison of strain types with prophage types suggests that the two may be related. The differences in GC content between the same type of prophages and the genomic region in which they were located indicates the alien properties of the prophages. The overall distribution of GC content suggests that prophages integrated on chromosomes and plasmids may have different evolutionary characteristics. These results suggest a high prevalence of prophages in the K. pneumoniae genome and highlight the effect of prophages on strain characterization.

Evaluating the effect of SARS-CoV-2 spike mutations with a linear doubly robust learner.

Driven by various mutations on the viral Spike protein, diverse variants of SARS-CoV-2 have emerged and prevailed repeatedly, significantly prolonging the pandemic. This phenomenon necessitates the identification of key Spike mutations for fitness enhancement. To address the need, this manuscript formulates a well-defined framework of causal inference methods for evaluating and identifying key Spike mutations to the viral fitness of SARS-CoV-2. In the context of large-scale genomes of SARS-CoV-2, it estimates the statistical contribution of mutations to viral fitness across lineages and therefore identifies important mutations. Further, identified key mutations are validated by computational methods to possess functional effects, including Spike stability, receptor-binding affinity, and potential for immune escape. Based on the effect score of each mutation, individual key fitness-enhancing mutations such as D614G and T478K are identified and studied. From individual mutations to protein domains, this paper recognizes key protein regions on the Spike protein, including the receptor-binding domain and the N-terminal domain. This research even makes further efforts to investigate viral fitness via mutational effect scores, allowing us to compute the fitness score of different SARS-CoV-2 strains and predict their transmission capacity based solely on their viral sequence. This prediction of viral fitness has been validated using BA.2.12.1, which is not used for regression training but well fits the prediction. To the best of our knowledge, this is the first research to apply causal inference models to mutational analysis on large-scale genomes of SARS-CoV-2. Our findings produce innovative and systematic insights into SARS-CoV-2 and promotes functional studies of its key mutations, serving as reliable guidance about mutations of interest.

Comparative genomics reveals new insights into the evolution of the IncA and IncC family of plasmids.

Incompatibility groups IncA and IncC plasmids are of great concern due to their ability to disseminate antibiotic resistance in bacteria via conjugative transfer. A deep understanding of their genomic structures and evolutionary characteristics is of great significance for improving our knowledge about its multidrug-resistance evolution and dissemination. However, current knowledge of their backbone structure, features of core functional modules and the characteristics of variable regions is based on a few plasmids, which highlights the need for a comprehensive systematic study. The present study thoroughly compared and analysed 678 IncA and IncC plasmid genomes. We found that their core functional genes were occasionally deficient and sometimes existed as multiple functional copies/multiple families, which resulted in much diversity. The phylogeny of 13 core functional genes corresponded well to the plasmid subtypes. The conjugative transfer system gained diverse complexity and exhibited many previously unnoticed types with multiple combinations. The insertion of mobile genetic elements (MGEs) in plasmids varied between types and was present in 4 insertion spots in different types of plasmids with certain types of transposons, integrons and insertion sequences. The impact of gene duplication, deletion, the insertion of MGEs, genome rearrangement and recombination resulted in the complex dynamic variable backbone of IncA and IncC plasmids. And IncA and IncC plasmids were more complex than their closest relative SXT/R391 integrative conjugative elements (ICEs), which included nearly all of the diversity of SXT/R391 in key systems. Our work demonstrated a global and systematic view of the IncA and IncC plasmids and provides many new insights into their genome evolution.

A Standardized Framework for Better Understanding of Phenotypic Differences within Bacterial Phyla Based on Protein Domain.

We propose a standardized framework to classify target species based on their protein domains, which can be utilized in different contexts, like eukaryotes and prokaryotes. In this study, by applying the framework to the bacterial kingdom as an implementation example and comparing the results with the current taxonomy standards at the phylum level, we came to the conclusion that the sequence of domains rather than the content of domains in a protein and the presence of one domain rather than the number of occurrences of one domain play more important roles in deciding bacterial phenotypes as well as matching the current taxonomy. In addition, the comparison also helps us to better focus on the species that conflict with the current phylum category, as well as to further investigate their phenotypic or genotypic differences. IMPORTANCE A 3-step framework was designed which can be applied to clustering species based on their protein domains, and different candidate models are proposed in each step for better adaptation of various scenarios. We show its implementation for the bacterial kingdom as an example, which helps us to find the most appropriate model combination that will best reflect the relationship between domains and phenotypes in this context. In addition, identifying species that are distant in the results but should be closely related phylogenetically can help us to focus on the mismatch for better understanding of their key phenotypic or genotypic differences.

A Vaccine Based on the Receptor-Binding Domain of the Spike Protein Expressed in Glycoengineered Pichia pastoris Targeting SARS-CoV-2 Stimulates Neutralizing and Protective Antibody Responses.

In 2020 and 2021, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel coronavirus, caused a global pandemic. Vaccines are expected to reduce the pressure of prevention and control, and have become the most effective strategy to solve the pandemic crisis. SARS-CoV-2 infects the host by binding to the cellular receptor angiotensin converting enzyme 2 (ACE2) via the receptor-binding domain (RBD) of the surface spike (S) glycoprotein. In this study, a candidate vaccine based on a RBD recombinant subunit was prepared by means of a novel glycoengineered yeast Pichia pastoris expression system with characteristics of glycosylation modification similar to those of mammalian cells. The candidate vaccine effectively stimulated mice to produce high-titer anti-RBD specific antibody. Furthermore, the specific antibody titer and virus-neutralizing antibody (NAb) titer induced by the vaccine were increased significantly by the combination of the double adjuvants Al(OH)3 and CpG. Our results showed that the virus-NAb lasted for more than six months in mice. To summarize, we have obtained a SARS-CoV-2 vaccine based on the RBD of the S glycoprotein expressed in glycoengineered Pichia pastoris, which stimulates neutralizing and protective antibody responses. A technical route for fucose-free complex-type N-glycosylation modified recombinant subunit vaccine preparation has been established.

In vivo structure and dynamics of the SARS-CoV-2 RNA genome.

The dynamics of SARS-CoV-2 RNA structure and their functional relevance are largely unknown. Here we develop a simplified SPLASH assay and comprehensively map the in vivo RNA-RNA interactome of SARS-CoV-2 genome across viral life cycle. We report canonical and alternative structures including 5'-UTR and 3'-UTR, frameshifting element (FSE) pseudoknot and genome cyclization in both cells and virions. We provide direct evidence of interactions between Transcription Regulating Sequences, which facilitate discontinuous transcription. In addition, we reveal alternative short and long distance arches around FSE. More importantly, we find that within virions, while SARS-CoV-2 genome RNA undergoes intensive compaction, genome domains remain stable but with strengthened demarcation of local domains and weakened global cyclization. Taken together, our analysis reveals the structural basis for the regulation of replication, discontinuous transcription and translational frameshifting, the alternative conformations and the maintenance of global genome organization during the whole life cycle of SARS-CoV-2, which we anticipate will help develop better antiviral strategies.

Characterizing the Core Internal Gene Pool of H9N2 Responsible for Continuous Reassortment With Other Influenza A Viruses.

Reassortment among avian influenza viruses is the main source of novel avian influenza virus subtypes. Studies have shown that the H9N2 virus often donates internal segments to generate novel reassortant avian influenza viruses, acting as a reassortment template. However, the characteristics of the internal pattern of reassortment remain unclear. In this article, we first defined the core gene pool of the internal segments of the H9N2 virus that provide templates for reassortment. We used genetic distance and sequence similarity to define typical clusters in the core gene pool. Then, we analyzed the phylogenetic relationships, feature vector distances, geographic distributions and mutation sites of strains related to the core gene pool. Strains in the same typical clusters have close phylogenetic relationships and feature vector distances. We also found that these typical clusters can be divided into three categories according to their main geographic distribution area. Furthermore, typical clusters in the same geographic area contain some common mutation patterns. Our results suggest that typical clusters in the core gene pool affect the reassortment events of the H9N2 virus in many respects, such as geographic distribution and amino acid mutation sites.

Reassortment Network of Influenza A Virus.

Influenza A virus (IAV) genomes are composed of eight single-stranded RNA segments. Genetic exchange through reassortment of the segmented genomes often endows IAVs with new genetic characteristics, which may affect transmissibility and pathogenicity of the viruses. However, a comprehensive understanding of the reassortment history of IAVs remains lacking. To this end, we assembled 40,296 whole-genome sequences of IAVs for analysis. Using a new clustering method based on Mean Pairwise Distances in the phylogenetic trees, we classified each segment of IAVs into clades. Correspondingly, reassortment events among IAVs were detected by checking the segment clade compositions of related genomes under specific environment factors and time period. We systematically identified 1,927 possible reassortment events of IAVs and constructed their reassortment network. Interestingly, minimum spanning tree of the reassortment network reproved that swine act as an intermediate host in the reassortment history of IAVs between avian species and humans. Moreover, reassortment patterns among related subtypes constructed in this study are consistent with previous studies. Taken together, our genome-wide reassortment analysis of all the IAVs offers an overview of the leaping evolution of the virus and a comprehensive network representing the relationships of IAVs.

Characterization of Burkholderia cepacia Complex Core Genome and the Underlying Recombination and Positive Selection.

Recombination and positive selection are two key factors that play a vital role in pathogenic microorganisms' population adaptation and diversification. The Burkholderia cepacia complex (Bcc) represents bacterial species with high similarity, which can cause severe infections among cases suffering from the chronic granulomatous disorder and cystic fibrosis (CF). At present, no genome-wide study has been carried out focusing on investigating the core genome of Bcc associated with the two evolutionary forces. The general characteristics of the core genome of Bcc species remain scarce as well. In this study, we explored the core orthologous genes of 116 Bcc strains using comparative genomic analysis and studied the two adaptive evolutionary forces: recombination and positive selection. We estimated 1005 orthogroups consisting entirely of single copy genes. These single copy orthologous genes in some Cluster of Orthologous Groups (COG) categories showed significant differences in the comparison of several evolutionary properties, and the encoding proteins were relatively simple and compact. Our findings showed that 5.8% of the core orthologous genes strongly supported recombination; in the meantime, 1.1% supported positive selection. We found that genes involved in protein synthesis as well as material transport and metabolism are favored by selection pressure. More importantly, homologous recombination contributed more genetic variation to a large number of genes and largely maintained the genetic cohesion in Bcc. This high level of recombination between Bcc species blurs their taxonomic boundaries, which leads Bcc species to be difficult or impossible to distinguish phenotypically and genotypically.

Nifurtimox Hampered the Progression of Astroglioma In vivo Via Manipulating the AKT-GSK3ß axis.

BACKGROUND: Astroglioma, one major form of brain tumors, has remained principally tough to handle for decades, due to the complexity of tumor pathology and the poor response to chemo- and radio-therapies. METHODS: Our previous study demonstrated that nifurtimox could regulate the signaling axis of AKT-GSK3ß in various tumor types including the astroglioma U251 cells. Intriguingly, earlier case studies suggested that nifurtimox could possibly permeate the blood brain barrier and arrest neuroblastoma in the brain. These observations jointly encouraged us to explore whether nifurtimox would hinder the growth of astroglioma in vivo. RESULTS: Our results exhibited that nifurtimox could competently hinder the development of astroglioma in the mouse brain as compared to temozolomide, the first line of drug for brain tumors. Meanwhile the surviving rate, as well as the body-weight was dramatically upregulated upon nifurtimox treatment, as compared to that of temozolomide. These findings offered nifurtimox as a better alternative drug in treating astroglioma in vivo. CONCLUSION: Persistently, the manipulation of the signaling axis of AKT-GSK3ß in astroglioma was found in line with earlier findings in neuroblastoma when treated with nifurtimox.