Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiologic agent of the coronavirus disease 2019 (COVID-19) pandemic. Clinical manifestations of the disease range from an asymptomatic condition to life-threatening events and death, with more severe courses being associated with age, male sex, and comorbidities. Besides these risk factors, intrinsic characteristics of the virus as well as genetic factors of the host are expected to account for COVID-19 clinical heterogeneity. Genetic studies have long been recognized as fundamental to identify biological mechanisms underlying congenital diseases, to pinpoint genes/proteins responsible for the susceptibility to different inherited conditions, to highlight targets of therapeutic relevance, to suggest drug repurposing, and even to clarify causal relationships that make modifiable some environmental risk factors. Though these studies usually take long time to be concluded and, above all, to translate their discoveries to patients' bedside, the scientific community moved really fast to deliver genetic signals underlying different COVID-19 phenotypes. In this Review, besides a concise description of COVID-19 symptomatology and of SARS-CoV-2 mechanism of infection, we aimed to recapitulate the current literature in terms of host genetic factors that specifically associate with an increased severity of the disease.

Biofilm-related infections are resistant forms of pathogens that are regarded as a medical problem, particularly due to the spread of multiple drug resistance. One of the factors associated with biofilm drug resistance is the presence of various types of efflux pumps in bacteria. Efflux pumps also play a role in biofilm formation by influencing Physical-chemical interactions, mobility, gene regulation, quorum sensing (QS), extracellular polymeric substances (EPS), and toxic compound extrusion. According to the findings of studies based on efflux pump expression analysis, their role in the anatomical position within the biofilm will differ depending on the biofilm formation stage, encoding gene expression level, the type and concentration of substrate. In some cases, the function of the efflux pumps can overlap with each other, so it seems necessary to accurate identify the efflux pumps of biofilm-forming bacteria along with their function in this process. Such studies will help to choose treatment strategy, at least in combination with antibiotics. Furthermore, if the goal of treatment is an efflux pump manipulation, we should not limit it to inhibition.

CRISPR-Cas systems are the only RNA- guided adaptive immunity pathways that trigger the detection and destruction of invasive phages and plasmids in bacteria and archaea. Due to its prevalence and mystery, the Class 1 CRISPR-Cas system has lately been the subject of several studies. This review highlights the specificity of CRISPR-Cas system III-A in Mycobacterium tuberculosis, the tuberculosis-causing pathogen, for over twenty years. We discuss the difference between the several subtypes of Type III and their defence mechanisms. The anti-CRISPRs (Acrs) recently described, the critical role of Reverse transcriptase (RT) and housekeeping nuclease for type III CRISPR-Cas systems, and the use of this cutting-edge technology, its impact on the search for novel anti-tuberculosis drugs.

The COVID-19 pandemic has brought significant challenges for genomic surveillance strategies in public health systems worldwide. During the past thirty-four months, many countries faced several epidemic waves of SARS-CoV-2 infections, driven mainly by the emergence and spread of novel variants. In that line, genomic surveillance has been a crucial toolkit to study the real-time SARS-CoV-2 evolution, for the assessment and optimization of novel diagnostic assays, and to improve the efficacy of existing vaccines. During the pandemic, the identification of emerging lineages carrying lineage-specific mutations (particularly those in the Receptor Binding domain) showed how these mutations might significantly impact viral transmissibility, protection from reinfection and vaccination. So far, an unprecedented number of SARS-CoV-2 viral genomes has been released in public databases (i.e., GISAID, and NCBI), achieving 14 million genome sequences available as of early-November 2022. In the present review, we summarise the global landscape of SARS-CoV-2 during the first thirty-four months of viral circulation and evolution. It demonstrates the urgency and importance of sustained investment in genomic surveillance strategies to timely identify the emergence of any potential viral pathogen or associated variants, which in turn is key to epidemic and pandemic preparedness.

Giardia duodenalis, Giardia enterica, Giardia intestinalis and Giardia lamblia are the synonyms for a species complex of 8-11 phylogenetically distinct species of Giardia infecting a broad range of animals including humans. Retrospective alignment of 8409 gene sequences from 3 loci confirmed host associations of Assemblages and sub-Assemblages within this species complex and molecular species delimitation testing confirmed that the Assemblages and sub-Assemblages AI and AII should be recognised as distinct species. It is recommended to synonymise the Assemblages with historic species descriptions based on host associations and consider descriptions for new species where no corresponding description exists. Synonyms, Giardia duodenalis, Giardia intestinalis and Giardia enterica, to be removed from synonymy: synonymise "Giardia duodenalis-Assemblage AI" syn. n. to Giardia duodenalis (Davaine, 1875), Kofoid and Christansen, 1915, synonymise "Giardia duodenalis-Assemblage AII" syn. n. to Giardia intestinalis (Lambl, 1859; Blanchard, 1885), Alexeieff, 1914 and synonymise "Giardia duodenalis-Assemblage B" syn. n. to Giardia enterica (Grassi, 1881), Kofoid, 1920. Host specific Assemblages synonymised: synonymise canid-associated "Giardia duodenalis-Assemblage C" syn. n. to Giardia canisHegner, 1922; synonymise artiodactyl-associated "Giardia duodenalis-Assemblage E" syn. n. to Giardia bovisFantham, 1921; synonymise feline-associated "Giardia duodenalis-Assemblage F" syn. n. to Giardia catiDeschiens, 1925; and synonymise rodent-associated "Giardia duodenalis-Assemblage G" syn. n. to Giardia simoniLavier, 1924. New description for parasite type infecting specific host: canid-associated "Giardia duodenalis-Assemblage D" named Giardia lupus, sp. n. (LSID: urn:lsid:zoobank.org:act:1651A8CB-CBA8-40D9-AB59-D4AB11AC18A3). New proposed names and descriptions for consideration for parasite types infecting specific hosts: cervid-associated "Giardia duodenalis-sub-Assemblage AIII" for consideration "cervus" and Pinnipedia-associated "Giardia duodenalis-Assemblage H" for consideration "pinnipedis".

Giardia is a genus of flagellated protozoan parasites that infect the small intestine of humans and animals, causing the diarrheal illness known as giardiasis. Giardia exhibits significant genetic diversity among its isolates, which can have important implications for disease transmission and clinical presentation. This diversity is influenced by the coevolution of Giardia with its host, resulting in the development of unique genetic assemblages with distinct phenotypic characteristics. Although panmixia has not been observed, some assemblages appear to have a broader host range and exhibit higher transmission rates. Molecular diagnostic methods enable researchers to examine the genetic diversity of Giardia populations, enhancing our understanding of the genetic diversity, population structure, and transmission patterns of this pathogen and providing insights into clinical presentations of giardiasis.

Malaria still poses a major burden on human health around the world, especially in endemic areas. Plasmodium resistance to several antimalarial drugs has been one of the major hindrances in control of malaria. Thus, the World Health Organization recommended artemisinin-based combination therapy (ACT) as a front-line treatment for malaria. The emergence of parasites resistant to artemisinin, along with resistant to ACT partner drugs, has led to ACT treatment failure. The artemisinin resistance is mostly related to the mutations in the propeller domain of the kelch13 (k13) gene that encodes protein Kelch13 (K13). The K13 protein has an important role in parasite reaction to oxidative stress. The most widely spread mutation in K13, with the highest degree of resistance, is a C580Y mutation. Other mutations, which are already identified as markers of artemisinin resistance, are R539T, I543T, and Y493H. The objective of this review is to provide current molecular insights into artemisinin resistance in Plasmodium falciparum. The trending use of artemisinin beyond its antimalarial effect is described. Immediate challenges and future research directions are discussed. Better understanding of the molecular mechanisms underlying artemisinin resistance will accelerate implementation of scientific findings to solve problems with malarial infection.

Aboriginal and Torres Strait Islander People (respectfully referred to as Indigenous Australians herein) are disparately burdened by many infectious and chronic diseases relative to Australians with European genetic ancestry. Some of these diseases are described in other populations to be influenced by the inherited profile of complement genes. These include complement factor B, H, I and complement factor H-related (CFHR) genes that can contribute to a polygenic complotype. Here the focus is on the combined deletion of CFHR1 and 3 to form a common haplotype (CFHR3-1Δ). The prevalence of CFHR3-1Δ is high in people with Nigerian and African American genetic ancestry and correlates to a higher frequency and severity of systemic lupus erythematosus (SLE) but a lower prevalence of age-related macular degeneration (AMD) and IgA-nephropathy (IgAN). This pattern of disease is similarly observed among Indigenous Australian communities. Additionally, the CFHR3-1Δ complotype is also associated with increased susceptibility to infection with pathogens, such as Neisseria meningitidis and Streptococcus pyogenes, which also have high incidences in Indigenous Australian communities. The prevalence of these diseases, while likely influenced by social, political, environmental and biological factors, including variants in other components of the complement system, may also be suggestive of the CFHR3-1Δ haplotype in Indigenous Australians. These data highlight a need to define the Indigenous Australian complotypes, which may lead to the discovery of new risk factors for common diseases and progress towards precision medicines for treating complement-associated diseases in Indigenous and non-Indigenous populations. Herein, the disease profiles suggestive of a common complement CFHR3-1Δ control haplotype are examined.

To systematically investigate the prophages carrying in Porphyromonas gingivalis (P. gingivalis) strains, analyze potential antibiotic resistance genes (ARGs) and virulence genes in these prophages. We collected 90 whole genome sequences of P. gingivalis from NCBI and utilized the Prophage Hunter online software to predict prophages; Comprehensive antibiotic research database (CARD) and virulence factors database (VFDB) were adopted to analyze the ARGs and virulence factors (VFs) carried by the prophages. Sixty-nine prophages were identified among 24/90 P. gingivalis strains, including 17 active prophages (18.9%) and 52 ambiguous prophages (57.8%). The proportion of prophages carried by each P. gingivalis genome ranged from 0.5% to 6.7%. A total of 188 antibiotic resistance genes belonging to 25 phenotypes and 46 different families with six mechanisms of antibiotic resistance were identified in the 17 active prophages. Three active prophages encoded 4 virulence genes belonging to type III and type VI secretion systems. The potential hosts of these virulence genes included Escherichia coli, Shigella sonnei, Salmonella typhi, and Klebsiella pneumoniae. In conclusion, 26.7% P. gingivalis strains carry prophages, while the proportion of prophage genes in the P. gingivalis genome is relatively low. In addition, approximately 39.7% of the P. gingivalis prophage genes have ARGs identified, mainly against streptogramin, peptides, and aminoglycosides. Only a few prophages carry virulence genes. Prophages may play an important role in the acquisition, dissemination of antibiotic resistance genes, and pathogenicity evolution in P. gingivalis.

Although the occurrence of three fiber genes in monkey adenoviruses had already been described, the relatedness of the "extra" fibers have not yet been discussed. Here we report the genome analysis of two simian adenovirus (SAdV) serotypes from Old World monkeys and the phylogenetic analysis of the multiple fiber genes found in these and related AdVs. One of the newly sequenced serotypes (SAdV-2), isolated from a rhesus macaque (Macaca mulatta), was classified into species Human mastadenovirus G (HAdV-G), while the other serotype (SAdV-17), originating from a grivet (Chlorocebus aethiops), classified to Simian mastadenovirus F (SAdV-F). We identified unique features in the gene content of these SAdVs compared to those typical for other members of the genus Mastadenovirus. Namely, in the E1B region of SAdV-2, the 19K gene was replaced by an ITR repetition and a copy of the E4 ORF1 gene. Among the 37 genes in both SAdVs, three genes of different lengths, predicted to code for the cellular attachment proteins (the fibers), were found. These proteins exhibit high diversity. Yet, phylogenetic calculations of their conserved parts could reveal the probable evolutionary steps leading to the multiple-fibered contemporary HAdV and SAdV species. Seemingly, there existed (a) common ancestor(s) with two fiber genes for the lineages of the AdVs in species SAdV-B, -E, -F and HAdV-F, alongside a double-fibered ancestor for today's SAdV-C and HAdV-G, which later diverged into descendants forming today's species. Additionally, some HAdV-G members picked up a third fiber gene either to the left-hand or to the in-between position from the existing two. A SAdV-F progenitor also obtained a third copy to the middle, as observed in SAdV-17. The existence of three fiber genes in these contemporary AdVs brings novel possibilities for the design of optimised AdV-based vectors with potential multiple target binding abilities.

Blastocystis sp., is an intestinal protist with a broad host range and a high prevalence in human populations worldwide, even in developed Western countries. The publication of conflicting evidence has divided the scientific community about the pathogenic role of this parasite. Even though, genetic studies on Blastocystis sp. revealed associations between genotypes and different pathogenic profiles. Conventionally, the detection of this parasite is based on microscopic or PCR methods, which offer meager or null performance in detecting mixed infections. In this work, we applied a metataxonomic NGS approach targeting the V4 region of the eukaryotic SSU-rRNA gene and classical phylogenetic methods. This approach allowed us to detect Blastocystis sp. in stool samples from infected children living in an urban setting in the city of Medellin attending the same daycare center. Phylogenetic analysis identified the subtypes present in the children as ST1, ST2, and ST3. Besides, mixed infections of subtypes ST1 + ST3 were spotted in 16% of the analyzed stool samples.

Bats remains as reservoirs for highly contagious and pathogenic viral families including the Coronaviridae, Filoviridae, Paramyxoviruses, and Rhabdoviridae. Spill over of viral species (SARS-CoV, MERS-CoV & SARS-CoV2) from bats (as a possible potential reservoirs) have recently caused worst outbreaks. Early detection of viral species of pandemic potential in bats is of great importance. We detected beta coronaviruses in the studied bats population (positive samples from Rousettus leschenaultia) and performed the evolutionary analysis, amino acid sequence alignment, and analysed the 3-Dimentional protein structure. We detected the coronaviruses for the first time in bats from Pakistan. Our analysis based on RdRp partial gene sequencing suggest that the studied viral strains are closely related to MERS-CoV-like viruses as they exhibit close structure similarities (with few substitutions) and also observed a substitution in highly conserved SDD in the palm subdomain of motif C to ADD, when compared with earlier reported viral strains. It could be concluded from our study that coronaviruses are circulating among the bat's population in Pakistan. Based on the current findings, we suggest large scale screening procedures of bat virome across the country to detect potential pathogenic viral species.

Over the past 20 years, the Seneca Valley virus (SVV) has emerged in various countries and regions around the world. Infected pigs display symptoms similar to foot-and-mouth disease and other vesicular diseases, causing severe economic losses to affected countries. In recent years, the number of SVV infections has been increasing in Brazil, China, and the United States. In this study, we comprehensively analyzed SVV genomic sequence data from the perspectives of evolutionary dynamics, phylogeography, and codon usage bias. We aimed to gain further insights into SVV's genetic diversity, spatiotemporal distribution patterns, and evolutionary adaptations. Phylogenetic analysis revealed that SVV has evolved into eight distinct lineages. Based on the results of phylogeographic analysis, it is speculated that the United States might have been the source of SVV, from where it subsequently spread to different countries and regions. Moreover, our analysis of positive selection sites in SVV capsid proteins suggests their potential importance in the process of receptor recognition. Finally, codon preference analysis indicates that natural selection has been a primary evolutionary driver influencing SVV codon usage bias. In conclusion, our in-depth investigation into SVV's origin, dissemination, evolution, and adaptation emphasizes the significance of SVV surveillance and control measures.

Picobirnaviruses are double-stranded RNA viruses known from a wide range of host species and locations but with unknown pathogenicity and host relationships. Here, we examined the diversity of picobirnaviruses from cattle and gorillas within and around Bwindi Impenetrable Forest National Park (BIFNP), Uganda, where wild and domesticated animals and humans live in relatively close contact. We use metagenomic sequencing with bioinformatic analyses to examine genetic diversity. We compared our findings to global Picobirnavirus diversity using clustering-based analyses. Picobirnavirus diversity at Bwindi was high, with 14 near-complete RdRp and 15 capsid protein sequences, and 497 new partial viral sequences recovered from 44 gorilla samples and 664 from 16 cattle samples. Sequences were distributed throughout a phylogenetic tree of globally derived picobirnaviruses. The relationship with Picobirnavirus diversity and host taxonomy follows a similar pattern to the global dataset, generally lacking pattern with either host or geography.

Malaria as vector-borne disease remains important health concern with over 200 million cases globally. Novel antimalarial medicines and more effective vaccines must be developed to eliminate and eradicate malaria. Appraisal of preceding genome editing approaches confirmed the CRISPR/Cas nuclease system as a novel proficient genome editing system and a tool for species-specific diagnosis, and drug resistance researches for Plasmodium species, and gene drive to control Anopheles population. CRISPR/Cas technology, as a handy tool for genome editing can be justified for the production of transgenic malaria parasites like Plasmodium transgenic lines expressing Cas9, chimeric Plasmodium transgenic lines, knockdown and knockout transgenic parasites, and transgenic parasites expressing alternative alleles, and also mutant strains of Anopheles such as only male mosquito populations, generation of wingless mosquitoes, and creation of knock-out/ knock-in mutants. Though, the incorporation of traditional methods and novel molecular techniques could noticeably enhance the quality of results. The striking development of a CRISPR/Cas-based diagnostic kit that can specifically diagnose the Plasmodium species or drug resistance markers is highly required in malaria settings with affordable cost and high-speed detection. Furthermore, the advancement of genome modifications by CRISPR/Cas technologies resolves contemporary restrictions to culturing, maintaining, and analyzing these parasites, and the aptitude to investigate parasite genome functions opens up new vistas in the better understanding of pathogenesis.

Genetic variation in the major histocompatibility complex (MHC) may be associated with resistance to the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). The pathogen originated in Asia, then spread worldwide, causing amphibian population declines and species extinctions. We compared the expressed MHC IIβ1 alleles of a Bd-resistant species, Bufo gargarizans, from South Korea with those of a Bd-susceptible Australasian species, Litoria caerulea. We found at least six expressed MHC IIβ1 loci in each of the two species. Amino acid diversity encoded by these MHC alleles was similar between species, but the genetic distance of those alleles known for potential broader pathogen-derived peptide binding was greater in the Bd-resistant species. In addition, we found a potentially rare allele in one resistant individual from the Bd-susceptible species. Deep next-generation sequencing recovered approximately triple the genetic resolution accessible from traditional cloning-based genotyping. Targeting the full MHC IIβ1 enables us to better understand how host MHC may adapt to emerging infectious diseases.

Mammalian orthoreoviruses (reoviruses) are currently classified based on properties of the attachment protein, σ1. Four reovirus serotypes have been identified, three of which are represented by well-studied prototype human reovirus strains. Reoviruses contain ten segments of double-stranded RNA that encode 12 proteins and can reassort during coinfection. To understand the breadth of reovirus genetic diversity and its potential influence on reassortment, the sequence of the entire genome should be considered. While much is known about the prototype strains, a thorough analysis of all ten reovirus genome segment sequences has not previously been conducted. We analyzed phylogenetic relationships and nucleotide sequence conservation for each of the ten segments of more than 60 complete or nearly complete reovirus genome sequences, including those of the prototype strains. Using these relationships, we defined genotypes for each segment, with minimum nucleotide identities of 77-88% for most genotypes that contain several representative sequences. We applied segment genotypes to determine reovirus genome constellations, and we propose implementation of an updated reovirus genome classification system that incorporates genotype information for each segment. For most sequenced reoviruses, segments other than S1, which encodes σ1, cluster into a small number of genotypes and a limited array of genome constellations that do not differ greatly over time or based on animal host. However, a small number of reoviruses, including prototype strain Jones, have constellations in which segment genotypes differ from those of most other sequenced reoviruses. For these reoviruses, there is little evidence of reassortment with the major genotype. Future basic research studies that focus on the most genetically divergent reoviruses may provide new insights into reovirus biology. Analysis of available partial sequences and additional complete reovirus genome sequencing may also reveal reassortment biases, host preferences, or infection outcomes that are based on reovirus genotype.

Plant pathogenic Pseudomonas species use multiple classes of toxins and virulence factors during host infection. The genes encoding these pathogenicity factors are often located on plasmids and other mobile genetic elements, suggesting that they are acquired through horizontal gene transfer to confer an evolutionary advantage for successful adaptation to host infection. However, the genetic rearrangements that have led to mobilization of the pathogenicity genes are not fully understood. In this study, we have sequenced and analyzed the complete genome sequences of four Pseudomonas amygdali pv. aesculi (Pae), which infect European horse chestnut trees (Aesculus hippocastanum) and belong to phylogroup 3 of the P. syringae species complex. The four investigated genomes contain six groups of plasmids that all encode pathogenicity factors. Effector genes were found to be mostly associated with insertion sequence elements, suggesting that virulence genes are generally mobilized and potentially undergo horizontal gene transfer after transfer to a conjugative plasmid. We show that the biosynthetic gene cluster encoding the phytotoxin coronatine was recently transferred from a chromosomal location to a mobilizable plasmid that subsequently formed a co-integrate with a conjugative plasmid.

Phage therapy is revitalized as an alternative to antibiotics therapy against antimicrobials resistant pathogens. Mycobacteriophages are genetically diverse viruses that can specifically infect Mycobacterium genus including Mycobacterium tuberculosis and Mycobacterium smegmatis. Here, we isolated and annotated the genome of a mycobacteriophage Lang, a temperate mycobacteriophage isolated from the soil of Hohhot, Inner Mongolia, China, by using Mycolicibacterium smegmatis mc2 155 as the host. It belongs to the Siphoviridae family of Caudovirales as determined by transmission electron microscopy. The morphological characteristics and certain biological properties of the phage were considered in detail. Phage Lang genomes is 41,487 bp in length with 66.85% GC content and encodes 60 putative open reading frames and belongs to the G1 sub-cluster. Genome annotation indicated that genes for structure proteins, assembly proteins, replications/transcription and lysis of the host are present in function clucters. The genome sequence of phage Lang is more than 95% similar to that of mycobacteriophage Grizzly and Sweets, differing in substitutions, insertions and deletions in Lang. One-step growth curve revealed that Lang has a latent period of 30 min and a outbreak period of 90 min. The short latent period and rapid outbreak mark the unique properties of phage Lang, which can be another potential source for combating M. tuberculosis.

Epstein-Barr virus (EBV) infection is extremely common worldwide, with approximately 90% of adults testing positive for EBV antibodies. Human are susceptible to EBV infection, and primary EBV infection typically occurs early in life. EBV infection can cause infectious mononucleosis (IM) as well as some severe non-neoplastic diseases, such as chronic active EBV infection (CAEBV) and EBV-associated hemophagocytic lymphohistiocytosis (EBV-HLH), which can have a heavy disease burden. After primary EBV infection, individuals develop robust EBV-specific T cell immune responses, with EBV-specific CD8+ and part of CD4+ T cells functioning as cytotoxic T cells, defending against virus. Different proteins expressed during EBV's lytic replication and latent proliferation can cause varying degrees of cellular immune responses. Strong T cell immunity plays a key role in controlling infection by decreasing viral load and eliminating infected cells. However, the virus persists as latent infection in EBV healthy carriers even with robust T cell immune response. When reactivated, it undergoes lytic replication and then transmits virions to a new host. Currently, the relationship between the pathogenesis of lymphoproliferative diseases and the adaptive immune system is still not fully clarified and needs to be explored in the future. Investigating the T cell immune responses evoked by EBV and utilizing this knowledge to design promising prophylactic vaccines are urgent issues for future research due to the importance of T cell immunity.

Japanese encephalitis (JE) is a major global public health threat. Using Japanese encephalitis incidence data from 2004 to 2010 in Guangxi Province, China, this study comprehensively explored the driving forces and the interactive effects between environmental and social factors of Japanese encephalitis using the Geo-detector method. The results indicated that the incidence of Japanese encephalitis showed a fluctuating downward trend from 2004 to 2010. The onset of JE was seasonal, mainly concentrated in June-July, and highly aggregated in northwestern Guangxi. Among the factors associated with Japanese encephalitis, days with temperatures >30 °C, accumulated temperatures >25 °C, slope, the normalized difference vegetation index, the gross domestic product of tertiary industries, the gross domestic product of primary industries and the number of pigs slaughtered showed higher contributions to Japanese encephalitis incidence. An enhanced interactive effect was found between environmental and social factors, and the interaction between days with humidity levels >80% and the gross domestic product of tertiary industries had the greatest combined effect on JE. These findings enhanced the understanding of the combined effect of social and environmental factors on the incidence of Japanese encephalitis and could help improve Japanese encephalitis transmission control and prevention strategies.

Rhipicephalus microplus, a hematophagous vector prevalent in the tropics and subtropics, is responsible for huge economic losses throughout the globe. However, the taxonomy of the tick species, especially prevalent in north India and south China has been challenged in the recent past. The present study attempted to assess the cryptic status of R. microplus ticks of north India based on two mitochondrial markers; the 16S rRNA and cox1 gene sequences. The phylogenetic tree corresponding to both markers demonstrated the presence of three distinct genetic assemblages/ clades of R. microplus. The present study isolates (n = five and seven for the cox1 and 16S rRNA gene sequences, respectively) from north India along with other isolates from India assorted in the R. microplus clade C sensu. Based on the median joining network analysis corresponding to the 16S rRNA gene sequence, 18 haplotypes were recorded, exhibiting a stellate shape, which was indicative of rapid population expansion. For the cox1 gene, the haplotypes corresponding to clades A, B and C were distantly placed with two exceptions. While performing the population structure analysis, low nucleotide (0.04745 ± 0.00416 and 0.01021 ± 0.00146) and high haplotype diversities (0.913 ± 0.032 and 0.794 ± 0.058) were recorded for the different clades of R. microplus based on the cox1 and 16S rRNA mitochondrial markers, respectively. Eventually, high genetic differentiation and low gene flow were recorded among the different clades. A negative value for the neutrality indices (Tajima's D = -1.44125, Fu's Fs = -4.879, Fu and Li's D = -2.78031 and Fu and Li's F = -2.75229) corresponding to the 16S rRNA gene for the overall dataset evinced an expansion of population size. Based on the detailed studies, it was inferred that the R. microplus tick species circulating in north India belonged to clade C sensu, similar to that of the species prevalent in the other parts of the country as well as in the Indian subcontinent.

Lagovirus europaeus/GI.1 is the virus that causes severe and dangerous rabbit haemorrhagic disease (RHD) in rabbits. Recombination formation in RHD viruses is common. Recombination is thought to play a key role in the evolution of lagoviruses and therefore most likely influences the pathogenicity of L. europaeus/GI strains. Immunological events also play a key role in the control of RHD, and an in-depth knowledge of these phenomena provides insights into the characteristics of the infection, which can help implement appropriate infection control measures. To obtain a more complete picture of RHD caused by different GI.1 strains, it is necessary to correlate the genetic diversity within L. europaeus/GI.1 strains and the immune picture in response to infection. We performed a phylogenetic analysis of the L. europaeus/GI strains and compared the recombinant L. europaeus/GI.1 strain with the GI.1a strain on the basis of a thorough statistical analysis of immunological traits performed previously. Our phylogenetic analysis based on the sequence of the gene encoding the VP60 capsid protein of 34 strains of Lagovirus europaeus showed that the Hartmannsdorf strain forms a separate clade from the other GI.1a strains and is separate from the GI.1b-d strains. Next, we showed significant differences in the levels of individual parameters for non-specific cellular and humoral immunity in infection with the GI.1a strain and the Hartmannsdorf recombinant strain. Against the background of this study, our results indicate that the characteristics of each recombinant should be considered individually.

Contrary to the global trend, between 2010 and 2020, an increase of 43% new HIV infections was recorded in Eastern Europe and Central Asia. Analyses of phylogenetic relationship, and routes and modes of transmission of the HIV-1 subtype B across the former Soviet Union (FSU) region are currently lacking. The objective of this analysis was to investigate the origin and transmission routes of HIV subtype B in FSU countries. We performed phylogenetic and phylodynamic analyses using 21,007 publicly available subtype B sequences from Europe and Asia, including thirteen FSU countries. Our study suggests that BFSU strain evolved more recently in FSU countries (Russia, Kyrgyzstan, Uzbekistan, Tajikistan, Belarus, Ukraine, Lithuania, Latvia, Georgia, Armenia, Azerbaijan) compared to the Western B variant in Western Europe (Austria, Belgium, Germany, Luxembourg, Netherlands, Switzerland). The primary high-risk group responsible for the transmission of subtype B was found to be MSM/homosexual. Intermixing of phylogenetic clusters among high-risk groups and bridging with the general population indicated that the HIV epidemic is no longer confined to distinct key populations - emphasizing an urgent need to improve the HIV harm-reduction efforts among high risk as well as general populations.

Contagious ecthyma in small ruminants is a zoonotic disease caused by Orf virus (ORFV) in the genus Parapoxvirus that can be deadly to its natural hosts. It causes significant losses worldwide, and commonly infects humans. However, the literature about its comparative severity in sheep and goat hosts is misleading; and while contagious ecthyma has been shown to occur in camels and transmit to humans, there is confusion as to whether ORFV is responsible. Camels are important from a 'One Health' perspective as they have been implicated as a reservoir host for the virus causing Middle East Respiratory Syndrome (MERS), which has a case fatality rate of 35% in humans. We compared ORFV gene sequences and mortality data from the West Bank in Palestine, where ORFV has not been reported previously, with data from the region. Surprisingly, we found that infections of camels that had been attributed to ORFV were more closely related to a different member of the genus Parapoxvirus. Two Middle East ORFVs isolated from humans were unrelated and sat alongside sheep and goat derived sequences on two distinct ORFV lineages of a maximum likelihood B2L gene tree. One of the viral lineages bifurcated to produce a monophyletic group of goat-derived ORFVs characterized uniquely by a glycine at amino acid reside 249. We found that serine is the ancestral allele shared between ORFV infections of sheep and also two closely related Parapoxviruses (PCPV and CCEV), indicating that the glycine allele represents a more recent shift in virus host range adaptation to goats. Furthermore, and contrary to some reports that ORFV is more severe in goats than in sheep, we observed median mortality of up to 24.5% in sheep, but none in goats. We also identified trans-boundary spread of ORFV between the West Bank and Israel.

Trypanosoma cruzi is the parasite responsible for Chagas disease. The parasite has been classified into six taxonomic assemblages: TcI-TcVI and TcBat (aka Discrete Typing Units or Near-Clades). No studies have focused on describing the genetic diversity of T. cruzi in the northwestern region of Mexico. Within the Baja California peninsula lives Dipetalogaster maxima, the largest vector species for CD. The study aimed to describe the genetic diversity of T. cruzi within D. maxima. A total of three Discrete Typing Units (DTUs) were found (TcI, TcIV, and TcIV-USA). TcI was the predominant DTU found (∼75% of samples), in concordance with studies from the southern USA, one sample was described as TcIV while the other ∼20% pertained to TcIV-USA, which has recently been proposed to have enough genetic divergence from TcIV, to merit its own DTU. Potential phenotype differences between TcIV and TcIV-USA should be assessed in future studies.

We analyzed the prevalence and genotypes of Parechovirus A (PeV-A) in children with diarrhea in Beijing, China, 2017-2019. A total of 1734 stool samples collected from children <5 years of age with diarrhea were tested for the presence of PeV-A. Viral RNA was detected by real-time RT-PCR, and then genotyped by nested RT-PCR. We detected PeV-A in 93 (5.4%, 93/1734) samples, of which 87 could be genotyped by amplification of either the complete or partial VP1 region or the VP3/VP1 junction region. The median age of PeV-A infected children was 10 months. Most PeV-A infections were observed between August and November, with a peak in September. Seven known genotypes of PeV-A1A, -A1B, -A3, -A4, -A6, -A8 and -A11 were detected and PeV-A1B was the most prevalent genotype. Coinfection with other diarrheal viruses was observed in 30.1% (28/93) of PeV-A positive samples. All strains of PeV-A1A, -A1B, -A4 and -A6 obtained in this study contained the arginine-glycine-aspartic acid (RGD) motif, while all strains of PeV-A3, -A8 and -A11 lacked it. This study revealed a high genetic diversity of PeV-A circulating in Beijing and PeV-A11 was reported for the first time in children with diarrhea in China.

Human adenoviruses (HAdVs) are the viral agents responsible for a wide spectrum of acute and chronic diseases. HAdVs are the most important etiological agents of acute gastroenteritis (AGE) and are identified as the major contributor to the deaths of diarrheal children globally. The significant rise in HAdV infections in rotavirus-vaccinated children documented in multiple studies demands continuous monitoring of HAdV strains. After the inclusion of rotavirus vaccines in the immunization schedule of India, public health research regarding prevalence, etiology, and risk factors is highly necessary for evidence-based policies and their implementation to sustain diarrhea prevention programs. In the present study, children admitted for AGE between 2013 and 2016 in seven different hospitals in Maharashtra and Gujrat states of Western India were subjected for investigation. HAdVs were found in 5.2% of the fecal specimens with the dominance of species-F (52.4%) strains, followed by the occurrence of non-enteric adenoviruses of species A (17.4%), C (11.4%), B (8.2%), and D (3.2%). The species-F strains were predominant in Ahmadabad (78.5%), Mumbai (61.5%), and Surat (57.1%) cities, followed by species-A strains. In Pune city, species B strains were detected in all HAdV patients, with none of the species A strains. Clinically, patients infected with enteric and non-enteric HAdV strains were indistinguishable. However, a high viral load was observed in species-F specimens as compared to non-species-F. The present study on fecal specimens collected in the pre-rotavirus vaccination era from hospitalized AGE patients will be important for future comparative analysis to know the exact impact of vaccination in children of Western India.

Outbreaks of HFMD in children aged <5 years have been reported worldwide and the major causative agents are Coxsackievirus (CV) A16, enterovirus (EV)-A71 and recently CVA6. In India, HFMD is a disease that is not commonly reported. The purpose of the study was to identify the enterovirus type(s) associated with large outbreak of Hand, foot, and mouth disease during COVID-19 pandemic in 2022. Four hundred and twenty five clinical samples from 196-suspected cases were collected from different parts of the country. This finding indicated the emergence of CVA6 in HFMD along with CVA16, soon after the gradual easing of non-pharmaceutical interventions during-pandemic COVID-19 and the relevance of continued surveillance of circulating enterovirus types in the post-COVID pandemic era.

The beak and feather disease virus (BFDV), causative agent of Psittacine beak and feather disease (PBFD), is a highly fatal and widespread virus that infects both the wild and captive Psittaciformes around the world. The BFDV genome is a ssDNA of approximately 2 kb in size, making it among the smallest known pathogenic viruses. Though, the virus is placed in Circoviridae family of the Circovirus genus, there is no classification system on clade and sub-clade level according to the International Committee on Taxonomy of Viruses and the strains are grouped on the bases of geographic locations. Thus, we provide the latest and robust phylogenetic classification of BFDVs in this study based on full-length genomic sequences, grouping all the available 454 strains detected during 1996-2022 into two distinct clades, e.g., GI and GII. The GI clade is further divided into six sub-clades (GI a-f), while GII into two sub-clades (GII a and b). In addition, the phylogeographic network identified high variability among the BFDV strains, showing several branches, where all the branches are connected to four strains, e.g., BFDV-ZA-PGM-70A(GenBank ID: HM748921.1, 2008-South Africa), BFDV-ZA-PGM-81A(GenBank ID: JX221009.1, 2008-South Africa), BFDV14(GenBank ID: GU015021.1, 2010-Thailand) and BFDV-isolate-9IT11(GenBank ID: KF723390.1, 2014-Italy). Furthermore, we identified 27 recombination events in the rep (replication-associated protein) and cap (capsid protein) coding regions using the complete genomes of BFDVs. Similarly, the amino acids variability analysis indicated that both the rep and cap regions are highly variable with values exceeding the variability coefficient estimation limit of 1.00, speculating the possible amino acids drift with the emergence of new strains. The findings provided in this study may offer the latest phylogenetic, phylogeographic and evolutionary landscape of the BFDVs.

Amphibian chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), emerged from Asia and spread globally. By comparing functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) with those of an Australasian Bd-susceptible species (Litoria caerulea), we identified MHC genotypes associated with Bd resistance. These alleles encode a glycine deletion (G90β1) and adjacent motifs in the deepest pathogen-derived peptide-binding groove. Every Bd-resistant individual, but no susceptible individuals, possessed at least one allele encoding the variant. We detected trans-species polymorphism at the end of the MHC IIβ1 sequences. The G90β1 deletion was encoded by different alleles in the two species, suggesting it may have evolved independently in each species rather than having been derived from a common ancestor. These results are consistent with a scenario by which MHC adaptations that confer resistance to the pathogen have evolved by convergent evolution. Immunogenetic studies such as this are critical to ongoing conservation efforts.

Hypervirulent carbapenem-resistant Klebsiella pneumoniae (hv-CRKP) has been widely reported and poses a global threat. However, the comprehensive genetic structure of ST11-KL64 hv-CRKP and the possible evolutionary mechanisms from a genetic structure perspective of this high-risk clone remain unclear. Here, a blaKPC-2-blaNDM-1-positive ST11-KL64 hv-CRKP isolate was obtained from a human bloodstream infection (BSI). Whole-genome sequencing and bioinformatics analyses revealed that it contained a fusion plasmid, pKPTCM2-1. pKPTCM2-1 is a conjugative plasmid composed of an oriT-positive pLVPK-like virulence plasmid and a type IV secretion system-produced blaNDM-1-bearing IncX3 plasmid mediated by IS26-based co-integration. This progress generated 8-bp target site duplications (TGAAAACC) on both sides. The fusion plasmid possessed self-transferability and could be transferred to blaKPC-2-harboring ST11-KL64 CRKP to form the ST11-KL64 hv-CRKP clone. The pLVPK-like-positive ST11-KL64 strain exhibited virulence levels similar to those of the typical hypervirulent K. pneumoniae NTUH-2044. The mutation, Tet(A) (A276S), which was believed to lead to tigecycline resistance was observed. Overall, this high-risk clone has emerged as a tremendous threat in fatal BSIs and thus, targeted surveillance is an urgent need to contain the hv-CRKP clones.

Canine circovirus (CanineCV) is a DNA virus affecting domestic dogs and other wild carnivore species. Despite the potential implications for dogs' health and wildlife conservation, data on CanineCV presence, epidemiology and genetic features from Africa is still poor. In the present study, biological specimens collected between 2020 and 2022 from a total of 32 jackals and 575 domestic dogs were tested for the presence of CanineCV DNA to evaluate its frequency. Furthermore, sequencing was conducted on positive samples to characterize the strains and compare them with publicly available sequences through phylogenetic analysis. A high CanineCV prevalence was observed both in jackals (43.75%; 95 CI: 28.17% - 60.67%) and domestic dogs (27.13%; 95 CI: 23.66% - 30.91%). All aside from one Namibian strain formed an independent clade, suggestive of extremely rare introduction events, followed by local persistence, circulation, and evolution. Remarkably, different recombination events were observed involving strains from both jackals and domestic dogs, which testify to the likely strain exchange between these populations. Distinctive amino acid residues were also observed in jackals. The limitations of the considered host populations however prevent a definitive conclusion on host adaptation, biological, and clinical features. Further studies should be performed to expand our current knowledge of the CanineCV disease scenario in Namibia, other African regions, and associated host species in Africa.

Polyomaviruses (PyVs) are known to infect a diverse range of vertebrate host species. We report the discovery of PyVs in vesper bats (family Vespertilionidae) from sampling in Central Europe. Seven partial VP1 sequences from different PyVs were detected in samples originating from six distinct vesper bat species. Using a methodology based on conserved segments within the major capsid virus protein 1 (VP1) among known PyVs, the complete genomes of two different novel bat PyVs were determined. The genetic distances of the large T antigen coding sequences from these PyVs compared to previously-described bat PyVs exceeded 15% meriting classification as representatives of two novel PyV species: Alphapolyomavirus epserotinus and Alphapolyomavirus myodaubentonii. Phylogenetic analysis revealed that both belong to the genus Alphapolyomavirus and clustered together with high confidence in clades including other bat alphapolyomaviruses reported from China, South America and Africa. In silico protein modeling of the VP1 subunits and capsid pentamers, and electrostatic surface potential comparison of the pentamers showed significant differences between the reference template (murine polyomavirus) and the novel bat PyVs. An electrostatic potential difference pattern between the two bat VP1 pentamers was also revealed. Disaccharide molecular docking studies showed that the reference template and both bat PyVs possess the typical shallow sialic acid-binding site located between two VP1 subunits, with relevant oligosaccharide-binding affinities. The characterisation of these novel bat PyVs and the reported properties of their capsid proteins will potentially contribute in the elucidation of the conditions creating the host-pathogen restrictions associated with these viruses.

Diarrheal cases caused by non-toxigenic Vibrio cholerae have been reported globally. Lineages L3b and L9, characterized as ctxAB-negative and tcpA-positive (CNTP), pose the highest risk and have caused long-term epidemics in different regions worldwide. From 2001 to 2018, two waves (2001-2012 and 2013-2018) of epidemic caused by non-toxigenic V. cholerae occurred in the developed city of Hangzhou, China. In this study, through the integrated analysis of 207 genomes of Hangzhou isolates from these two waves (119 and 88) and 1573 publicly available genomes, we showed that L3b and L9 lineages together caused the second wave as had happened in the first wave, but the dominant lineage shifted from L3b (first wave: 69%) to L9 (second wave: 50%). We further found that the genotype of a key virulence gene, tcpF, in the L9 lineage during the second wave shifted to type I, which may have enhanced bacterial colonization in humans and potentially promoted the pathogenic lineage shift. Moreover, we found that 21% of L3b and L9 isolates had changed to predicted cholera toxin producers, providing evidence that gain of complete CTXφ-carrying ctxAB genes, rather than ctxAB gain in pre-CTXφ-carrying isolates, led to the transition. Taken together, our findings highlight the possible public health risk associated with L3b and L9 lineages due to their potential to cause long-term epidemics and turn into high-virulent cholera toxin producers, which necessitates a more comprehensive and unbiased sampling in further disease control efforts.

Monkeypox virus (MPXV) belongs to the Orthopoxvirus genus. The worldwide outbreak of MPXV in 2022 has caused widespread concerns. Cross-reactive antibodies induced by vaccinia-inoculation can provide protection against reinfection by MPXV. The vaccinia Tian Tan (VTT) strain, which was widely inoculated in the Chinese population before the 1980s, has genomic differences from other vaccinia strains, although they all belong to the orthopoxviruses family. The current seroprevalence of VTT-vaccinated populations remains unclear more than four decades after the termination of vaccination campaigns in China. Our results showed that cross-reactive IgG antibodies against MPXV were present in 31.8% (75/236) of vaccinees four decades after VTT-vaccination, suggesting that vaccination with VTT may provide long-term protection against MPXV infection in some individuals.

Salmonella enterica serovar Dublin is highly adapted to cattle and a relatively rare cause of human infections. In Denmark S. Dublin has been endemic in the cattle population for many years. A national surveillance program in the cattle population was established at herd-level to reduce the occurrence of S. Dublin. In this study, we analyzed 421 S. Dublin genomes from cattle and food in order to determine the trend of S. Dublin's population size over time in Denmark and the impact of intervention in the cattle industry on the bacterial population size. A phylogenetic tree based on SNPs exhibited two major clades and one small cluster. All isolates were ST10. The temporal phylogenetic tree for the S. Dublin isolates showed that the most recent common ancestor was estimated to be in ∼1980 for the two major clades. An effective population size over time based on a Bayesian skyline plot showed that the population size of S. Dublin decreased significantly between 2014 and 2019 in both major clades. This result was concordant with the decrease of infected human cases by S. Dublin in Denmark. The strengthening of a surveillance program in Denmark could be the cause for the reduction of S. Dublin's effective population size. This study showed that whole genome sequencing combined with computer intensive phylogenetic analysis estimating the effective size of the S. Dublin's population over time is a strongly relevant measure with respect to assessing the impact of control measures aiming to reduce the bacterial population in the reservoir and the risk for human infection.

Mycobacteroides abscessus complex (MAB), a fast-growing nontuberculous mycobacterium, is emerging as a significant infectious disease threat, due to both intrinsic and acquired resistance mechanisms to antibiotics and disinfectants and the need for extensive and multidrug regimens for treatment. Despite the prolonged regimens, outcomes are poor and persistence cases have been reported. Here, we describe clinical, microbiologic and genomic features of a M. abscessus subsp. bolletii (M. bolletii) strain consecutively isolated from a patient within an eight-year infection period. From April 2014 to September 2021, the National Reference Laboratory for Mycobacteria received eight strains isolated from a male patient. Species identification, molecular resistance profile and phenotypic drug susceptibility were determined. Five of these isolates were recovered for further in-depth genomic analysis. Genomic analysis confirmed the multidrug resistant pattern of the strain and also other genetic changes associated with adaptation to environment and defence mechanisms. We highlight the identification of new mutations in locus MAB_1881c and in locus MAB_4099c (mps1 gene), already described as associated with macrolides resistance and morphotype switching, respectively. Additionally, we also observed the emergence and fixation of a mutation in locus MAB_0364c that appeared at a frequency of 36% for the 2014 isolate, 57% for the 2015 isolate and 100% for the 2017 and 2021 isolates, clearly illustrating a fixation process underlying a microevolution of the MAB strain within the patient. Altogether these results suggest that the observed genetic alterations are a reflection of the bacterial population's continuous adaptation and survival to the host environment during infection, contributing to persistence and treatment failure.

Influenza A viruses (IAV) are fast-evolving pathogens with a very high mutation rate (2.0 × 10-6 to 2.0 × 10-4) compared to the influenza B (IBV) and influenza C (ICV) viruses. Generally, the tropical regions are considered as the reservoir for the IAV's genetic and antigenic evolutionary modification to be reintroduced into the temperate region. Therefore, in connection to the above facts, the present study emphasized on the evolutionary dynamic of the pandemic-2009 H1N1 (pdmH1N1) influenza virus in India. A total of Ninety-two whole genome sequences of pdmH1N1 viruses circulating in India during the 2009 post-pandemic era were analysed. The temporal signal of the study, indicating strict molecular clock evolutionary process and the overall substitution rate is 2.21 × 10-3/site/year. We are using the nonparametric Bayesian Skygrid coalescent model to estimates the effective past population dynamic or size over time. The study exhibits a strong relation between the genetic distances and collection dates of the Indian pdmH1N1 strain. The skygrid plot represents the highest exponential growth of IAV in rainy and winter seasons. All the genes of Indian pdmH1N1 were under purifying selective pressure. The Bayesian time-imprinted phylogenetic tree represents the following clade distributions in the country within the last 10 years; I) clade 6, 6C, and 7 were co-circulating between the 2011 to 2012 flu season; II) the clade 6B was introduced into circulation in the late seasons of 2012; III) lastly, the clade 6B remain existing in the circulation and segregated into subclade 6B.1 with five different subgroup (6B.1A, 6B.1A.1, 6B.1A.5a, 6B.1A.5a.2, 6B.1A.7). The recent circulating strain of Indian H1N1 strain represent the insertion of basic-amino acid arginine (R) in the cleavage site (325/K-R) of HA protein and amino acid mutation (314/I-M) on the lateral head surface domain of NA protein. Moreover, the study indicates the sporadic presence of the oseltamivir-resistant (275/H-Y) H1N1 variant in circulation. The present study suggests the purifying selective pressure and stochastic ecological factors for the existence and adaptation of a certain clade 6B in the host populations and additional information on the emergence of mutated strains in the circulation.

OBJECTIVES: Hepatitis A (HAV) virus causes asymptomatic to life-treating fulminant hepatitis. During infection, patients show large viral excretion in their stools. Resistance of HAV to environmental conditions, allows us to recover viral nucleotide sequences from wastewater and trace its evolutionary history. METHODS: We characterize twelve years of HAV circulation in wastewater from Santiago, Chile, and conducted phylogenetic analyses to decipher the dynamics of circulating lineages. RESULTS: We observed the exclusive circulation of the HAV IA genotype. The molecular epidemiologic analyses showed a steady circulation of a dominant lineage with low genetic diversity (d = 0,007) between 2010 and 2017. An outbreak of Hepatitis A associated with men who have sex with men, in 2017 was associated with the irruption of a new lineage. Remarkably, a dramatic change in the dynamic of HAV circulation was observed in the period post-outbreak; between 2017 and 2021 when 4 different lineages were transiently detected. Exhaustive phylogenetic analyses indicate that these lineages were introduced and possibly derived from isolates from other Latin American countries. CONCLUSION: The HAV circulation in recent years in Chile is rapidly changing and suggests that this phenomenon could be a consequence of massive population migrations in Latin America caused by political instability and natural disasters.

Newcastle disease (ND) is a highly contagious viral disease that affects many bird species worldwide. This study presents the results of the molecular characterization and phylogenetic analysis of 15 virulent ND viruses (NDV) isolated from chickens during outbreaks reported in 2016 and 2018, in the provinces of Namibe and Huíla, in southern Angola. A 561-nucleotide fragment of the F gene was amplified by RT-PCR and sequenced for molecular characterization. Results showed that in all isolates the amino acid sequence comprising the cleavage site of fusion protein is characteristic of virulent viruses (RRQKR/F). Blast analysis revealed high similarity (99.2%) between two isolates from Huíla province, HLA4 and HLA6, and strain 5620 (GenBank accession number KY747479) isolated from chickens in the neighboring country Namibia, in 2016. The other isolates investigated are more related (97.0%) with strain 6195 (GenBank accession number KY747480), also isolated in Namibia in 2016. Phylogenetic analysis performed by Maximum Likelihood, Neighbor-joining and Bayesian methods revealed that like the strains isolated in Namibia, the isolates from southern Angola also belong to subgenotype 2 of genotype VII (VII.2). The network analysis revealed that NBA1 isolate from Angola is closer to a common ancestor than the isolates from Namibia, suggesting that transmission of ND viruses occurred from Angola to Namibia.

Coquillettidia mosquitoes are important nuisance-biting pests and a vector of brugian filariasis in Thailand. However, comprehensive information about these mosquitoes remains unavailable such as molecular and morphometric differences among species. The lack of vector knowledge on Coquillettidia species could affect future disease control. This study aims to investigate differences in molecular variations based on mitochondrial cytochrome oxidase subunit I (COI) gene and wing geometric traits of three Coquillettidia species, namely Cq. crassipes, Cq. nigrosignata, and Cq. ochracea in Thailand. The results of molecular analyses revealed the differences among three Coquillettidia species. The genetic difference measure based on the Kimura two-parameter model among three Coquillettidia species showed low intraspecific distances (0%-3.05%) and large interspecific distances (10.10%-12.41%). The values of intra- and inter-genetic differences of three Coquillettidia species did not overlap which showed the existence of a barcoding gap indicating the efficiency of the identification based on the COI gene. As with molecular analysis, the landmark-based geometric morphometrics approach based on wing shape analysis indicated three distinct species groups which were supported by the high total performance score of cross-validated classification (97.16%). These results provide the first evidence of taxonomic signal based on molecular and wing geometric differences to support species identification and biological variations of Coquillettidia mosquitoes in Thailand for understanding these rare vector mosquitoes in depth and leading to effective further mosquito control.

Bovine Parainfluenza Type 3 virus (BPIV-3) is an enveloped, non-segmented single-stranded, negative-sense RNA virus belonging to the Paramyxoviridae family (genus Respirovirus) with a well-known role in Bovine Respiratory Disease (BRD) onset. Being isolated for the first time in 1959, BPIV-3 currently circulates worldwide in cattle herds and is routinely tested in suspected BRD cases. Different commercial vaccines are available to prevent infection and/or to reduce the clinical signs associated with BPIV-3 infection, which are essential to prevent secondary infections. Despite years of molecular surveillance, a very limited number of complete genome sequences were made publicly available, preventing thus the understanding of the genetic diversity of the circulating strains in the field. In addition, no data about the genetic identity between field and vaccine strains is currently available. In this study, we sequenced the full-genome and genetically characterized BPIV-3 strains isolated from animals displaying respiratory illness in France and Sweden, as well as the vaccine strains contained in three different commercialized vaccines. Our results show that the sequences from France and Sweden belong to genotype C. However, a third sequence from Sweden from 2017 clustered within genotype A. The sequencing of vaccine strains revealed that two of the vaccine strains clustered within genotype C, whereas the third vaccine strain belonged to genotype A. Altogether, our findings suggest that both genotypes A and C circulate in Europe and that BPIV-3 field and vaccine strains are genetically divergent. Our sequencing results could be useful to better understand the genetic differences between the circulating field and vaccine BPIV-3 strains. This is crucial for a correct interpretation of diagnostic findings and for the assessment of BPIV-3 prevalence in cattle population.

OBJECTIVES: Clustering pathogen sequence data is a common practice in epidemiology to gain insights into the genetic diversity and evolutionary relationships among pathogens. We can find groups of cases with a shared transmission history and common origin, as well as identifying transmission hotspots. Motivated by the experience of clustering SARS-CoV-2 cases using whole genome sequence data during the COVID-19 pandemic to aid with public health investigation, we investigated how differences in epidemiology and sampling can influence the composition of clusters that are identified. METHODS: We performed genomic clustering on simulated SARS-CoV-2 outbreaks produced with different transmission rates and levels of genomic diversity, along with varying the proportion of cases sampled. RESULTS: In single outbreaks with a low transmission rate, decreasing the sampling fraction resulted in multiple, separate clusters being identified where intermediate cases in transmission chains are missed. Outbreaks simulated with a high transmission rate were more robust to changes in the sampling fraction and largely resulted in a single cluster that included all sampled outbreak cases. When considering multiple outbreaks in a sampled jurisdiction seeded by different introductions, low genomic diversity between introduced cases caused outbreaks to be merged into large clusters. If the transmission and sampling fraction, and diversity between introductions was low, a combination of the spurious break-up of outbreaks and the linking of closely related cases in different outbreaks resulted in clusters that may appear informative, but these did not reflect the true underlying population structure. Conversely, genomic clusters matched the true population structure when there was relatively high diversity between introductions and a high transmission rate. CONCLUSION: Differences in epidemiology and sampling can impact our ability to identify genomic clusters that describe the underlying population structure. These findings can help to guide recommendations for the use of pathogen clustering in public health investigations.

Since the late 1930s, resistance to sulfonamides has been accumulating across bacterial species including Acinetobacter baumannii, an opportunistic pathogen increasingly implicated the spread of antimicrobial resistance worldwide. Our study aimed to explore events involved in the acquisition of sulfonamide resistance genes, particularly sul2, among the earliest available isolates of A. baumannii. The study utilized the genomic data of 19 strains of A. baumannii isolated before 1985. The whole genomes of 5 clinical isolates obtained from the Culture Collection University of Göteborg (CCUG), Sweden, were sequenced using the Illumina MiSeq system. Acquired resistance genes, insertion sequence elements and plasmids were detected using ResFinder, ISfinder and Plasmidseeker, respectively, while sequence types (STs) were assigned using the PubMLST Pasteur scheme. BLASTn was used to verify the occurrence of sul genes and to map their genetic surroundings. The sul1 and sul2 genes were detected in 4 and 9 isolates, respectively. Interestingly, sul2 appeared thirty years earlier than sul1. The sul2 gene was first located in the genomic island GIsul2 located on a plasmid, hereafter called NCTC7364p. With the emergence of international clone 1, the genetic context of sul2 evolved toward transposon Tn6172, which was also plasmid-mediated. Sulfonamide resistance in A. baumannii was efficiently acquired and transferred vertically, e.g., among the ST52 and ST1 isolates, as well as horizontally among non-related strains by means of a few efficient transposons and plasmids. Timely acquisition of the sul genes has probably contributed to the survival skill of A. baumannii under the high antimicrobial stress of hospital settings.

Outer inflammatory protein A (OipA), which is encoded by the oipA gene, can induce interleukin-8 secretion in gastric epithelial cells. The functional status of the oipA gene is regulated by the slipped-strand mispairing mechanism based on the CT dinucleotide repeat number in the 5' region. This study aimed to investigate the oipA functional status ("on/off") of Helicobacter pylori (H. pylori) and its association with gastroduodenal diseases in southwestern Vietnam. The cross-sectional study was conducted on 173H. pylori isolates from 173 patients with gastroduodenal diseases. Sanger sequencing was used to determine the functional status of oipA. Multivariable logistic regression analysis was performed to identify the association between oipA status and gastroduodenal diseases. The oipA "on" status accounted for 96% of H. pylori isolates. Twenty-five CT repeat patterns of the oipA 5' signal region were observed, five of which were novel CT repeat patterns. The oipA "on" status was found in 100%, 97.8%, and 86.8% of H. pylori isolates from patients with peptic ulcer, precancerous lesions, and chronic gastritis, respectively (p < 0.01). The oipA "on" status was related to gastric precancerous lesions versus chronic gastritis (adjusted OR = 7.39, 95% CI: 1.35-40.59, p = 0.021) and peptic ulcers versus chronic gastritis (adjusted OR = 12.79, 95% CI: 1.19-1760.32, p = 0.033). Our data show a high prevalence of the oipA "on" status, which was associated with precancerous gastric lesions and peptic ulcers. Moreover, genetic diversity in the number and pattern of CT dinucleotide repeat of oipA among Vietnamese H. pylori strains was identified.

Leptospirosis, a major zoonotic disease caused by pathogenic Leptospira spp. is recognized globally as an emerging zoonotic disease. Whole-genome sequencing reveals hidden messages about Leptospira's pathogenesis. We used Single Molecule Real-Time (SMRT) sequencing to obtain complete genome sequences of twelve L. interrogans isolates from febrile patients from Sri Lanka for a comparative whole genome sequencing study. The sequence data generated 12 genomes with a coverage greater than X600 with sizes ranging from 4.62 Mb to 5.16 Mb, and a G + C content ranging from 35.00% to 35.42%. The total number of coding sequences predicted by the NCBI (National Center for Biotechnology Information) genome assembly platform ranged from 3845 to 4621 for the twelve strains. Leptospira serogroup with similar-sized LPS biosynthetic loci that belonged to the same clade had a close relationship in the phylogenetic analysis. Nonetheless, variations in the genes encoding sugar biosynthesis were found in the serovar determinant region (rfb locus). Type I and Type III CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) systems were found in all of the strains. Genome BLAST Distance Phylogeny of these sequences allowed for detailed genomic strain typing. These findings may help us better understand the pathogenesis, develop a tools for early diagnosis, comparative genomic analysis and evolution of Leptospira.

Genotyping and virulence studies of Toxoplasma gondii are essential to investigate the pathogenesis of strains circulating worldwide. In this study, eight T. gondii isolates obtained from a congenitally infected newborn, a calf, two cats, three dogs, and a wallaby from five states of México were genotyped by Mn-PCR-RFLP with 11 typing markers (SAG1, SAG2 5'3', alt. SAG2, SAG3, BTUB, GRA6, c22-8, c29-2, L358, PK1 and Apico), five virulence markers (CS3, ROP16, ROP17, ROP18 and ROP5), 15 microsatellite markers (TUB-2, W35, TgM-A, B18, B17, M33, IV.1, XI.1, M48, M102, N60, N82, AA, N61, N83), and sequencing. A phylogenetic network was built to determine the relationship between Mexican isolates and those reported worldwide. Six different genotypes were identified by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), ToxoDB #8, #10, #28 (n = 3), #48, #116, and #282. Genotyping by microsatellite analysis differentiated the three PCR-RFLP genotype #28 isolates into two strains, revealing a total of seven microsatellite genotypes. Three different allele combinations of ROP18/ROP5 virulence markers were also found, 3/3, 1/1, and 4/1. The last two combinations are predicted to be highly virulent in the murine model. According to the phylogenetic network, the T. gondii strains studied here are related to archetypal strains I and III, but none are related to the strains previously reported in México. The genotypes identified in this study in different species of animals demonstrate the great genetic diversity of T. gondii in México. The ToxoDB-PCR-RFLP #28 genotype was found in three isolates from different hosts and states. Additionally, four of the isolates are predicted to be highly virulent in mice. The next step will be to perform in vitro and in vivo assays to determine the phenotype of these T. gondii isolates in murine models.

We retrospectively investigated CRKP isolates among 92 pediatric patients (32 neonates and 60 non‑neonates) in 2019 and 2020 (59 and 33 isolates, respectively) to investigate the molecular characteristics and virulence factors of carbapenem-resistant Klebsiella pneumoniae (CRKP) isolated from pediatric patients,. All the CRKP isolates were subjected to antimicrobial susceptibility testing, string testing, molecular typing of virulence and carbapenemase genes, and multilocus sequence typing. Hypervirulent K. pneumoniae (Hvkp) was defined based on the detection of the regulator of mucoid phenotype A (rmpA).Sequence type 11 (ST11) accounted for the majority of infections in both neonates (37.5%) and non‑neonates (43.3%) (P > 0.05), whereas it increased from 30.5% (18/59) in 2019 to 60.6% (20/33) in 2020 (P < 0.05). Carbapenemase gene KPC-2 was predominant in both neonates and non‑neonates (46.9% vs. 51.7%, respectively), followed by New Delhi metallo-beta-lactamase 1 (NDM-1) (34.4% vs. 28.3%, respectively) (all P > 0.05). Compared to 2019, the proportion of blaNDM-1 decreased (44.1% vs. 6.1%) (P < 0.001), while that of blaKPC-2 increased (40.7% vs. 66.7%) (P = 0.017) in 2020. ybtS and iutA had a higher positivity rate in KPC-2 and ST11 producers (all P < 0.05); the KPC-2-, ybtS-, and iutA-positive isolates showed relatively higher resistance to fluoroquinolones and aminoglycosides, nitrofurantoin, and piperacillin/tazobactam, respectively. Furthermore, the combined expression (95.7%, 88/92) of carbapenemase and virulence-associated genes was detected, with the carbapenemase genes blaKPC-2 and blaTEM-1 combined with virulence-associated genes entB, mrkD, and ybtS accounting for the highest percentage (20.7%).Carbapenemase gene mutations in the CRKP strain from 2019 to 2020 highlight the importance of dynamic monitoring. The spread of hypervirulence-associated genes in CRKP strains and the high positivity rates of ybtS and iutA in KPC-2- and ST11-producing ones signify their high virulence potential in pediatric patients.

The wide-spread of drug-resistant Acinetobacter baumannii is a global health problem. This study investigated the clonal distribution and antimicrobial resistance of 167 A. baumannii isolates from two Korean university hospitals from 2009 to 2019 by analyzing the sequence types (STs), antimicrobial resistance, and resistance determinants of carbapenems and aminoglycosides. Twenty STs, including 16 pre-existing STs and four unassigned STs, were identified in A. baumannii isolates using the Oxford multilocus sequence typing scheme. Two STs, ST191 (n = 77) and ST451 (n = 40), were prevalent, and majority (n = 153) of the isolates belonged to clonal complex 208. The ST191 isolates were detected during the study period, whereas ST451 isolates were detected after 2016. One hundred forty-seven (87%) of 167 A. baumannii isolates were non-susceptible to carbapenems. The ST191 and ST451 isolates exhibited higher resistance to antimicrobial agents than that of the sporadic ST isolates. Interestingly, ST451 isolates exhibited lower susceptibility to minocycline and tigecycline than the other ST isolates. All carbapenem-non-susceptible A. baumannii isolates, except four, carried the ISAbaI-blaOXA-23 structure. armA was detected in all amikacin-non-susceptible isolates (n = 128) except for one isolate. Five aminoglycoside-modifying enzyme (AME) genes were detected, but their carriage varied between STs; ant(3″)-Ia and aac(6')-Ib were more common in ST191 than in ST451, while aph(3')-Ia was more common in ST451 than in ST191. This study demonstrated the clonal evolution related to antimicrobial resistance in A. baumannii.

DENV-2 was the main responsible for a 70% increase in dengue incidence in Brazil during 2019. That year, our metagenomic study by Illumina NextSeq on serum samples from acute febrile patients (n = 92) with suspected arbovirus infection, sampled in 22 cities of the state of Mato Grosso (MT), in the middle west of Brazil, revealed eight complete genomes and two near-complete sequences of DENV-2 genotype III, one Human parvovirus B19 genotype I (5,391 nt) and one Coxsackievirus A6 lineage D (4,514 nt). These DENV-2 sequences share the aminoacidic identities of BR4 lineage on E protein domains I, II and III, and were included in a clade with sequences of the same lineage circulating in the southeast of Brazil in the same year. Nevertheless, 11/34 non-synonymous mutations are unique to three strains inthis study, distributed in the E (n = 6), NS3 (n = 2) and NS5 (n = 3) proteins. Other 14 aa changes on C (n = 1), E (n = 3), NS1 (n = 2), NS2A (n = 1) and NS5 (n = 7) were first reported in a genotype III lineage, having been already reported only in other DENV-2 genotypes. All 10 sequences have mutations in the NS5 protein (14 different aa changes). Nine E protein aa changes found in two sequences, six of which are unique, are in the ectodomain; where the E:M272T change is on the hinge of the E protein at domain II, in a region critical for the anchoring to the host cell receptor. The NS5:G81R mutation, in the methyltransferase domain, was found in one strain of this study. Altogether, these data points to an important evolution of DENV-2 genotype III lineage BR4 during this outbreak in 2019 in MT. Genomic surveillance is essential to detect virus etiology and evolution, possibly related to immune evasion and viral fitness changes leading to future novel outbreaks.

Human movement may be an important driver of transmission dynamics for enteric pathogens but has largely been underappreciated except for international 'travelers' diarrhea or cholera. Phylodynamic methods, which combine genomic and epidemiological data, are used to examine rates and dynamics of disease matching underlying evolutionary history and biogeographic distributions, but these methods often are not applied to enteric bacterial pathogens. We used phylodynamics to explore the phylogeographic and evolutionary patterns of diarrheagenic E. coli in northern Ecuador to investigate the role of human travel in the geographic distribution of strains across the country. Using whole genome sequences of diarrheagenic E. coli isolates, we built a core genome phylogeny, reconstructed discrete ancestral states across urban and rural sites, and estimated migration rates between E. coli populations. We found minimal structuring based on site locations, urban vs. rural locality, pathotype, or clinical status. Ancestral states of phylogenomic nodes and tips were inferred to have 51% urban ancestry and 49% rural ancestry. Lack of structuring by location or pathotype E. coli isolates imply highly connected communities and extensive sharing of genomic characteristics across isolates. Using an approximate structured coalescent model, we estimated rates of migration among circulating isolates were 6.7 times larger for urban towards rural populations compared to rural towards urban populations. This suggests increased inferred migration rates of diarrheagenic E. coli from urban populations towards rural populations. Our results indicate that investments in water and sanitation prevention in urban areas could limit the spread of enteric bacterial pathogens among rural populations.

Malaria in India is declining, in part due to the use of long-lasting insecticide-treated nets (LLINs) and vector control. Historically, the north-eastern region of India has contributed ~10%-12% of the nation's malaria burden. The important mosquito vectors in northeast India have long been considered to be Anopheles baimaii and An. minimus, both associated with forest habitats. Local deforestation and increased rice cultivation, along with widespread LLIN use, may be changing vector species composition. Understanding if and how vector species composition is changing is critical to successful malaria control. In Meghalaya state, malaria is now at a low level of endemicity with occasional seasonal outbreaks. In a biodiverse setting like Meghalaya, where >24 Anopheles mosquito species have been recorded, accurate morphological identification of all species is logistically challenging. To accurately determine Anopheles species richness in the West Khasi Hills (WKH) and West Jaintia Hills (WJH) districts, adult and larval mosquitoes were collected and identified using molecular methods of allele-specific PCR and cytochrome oxidase I DNA barcoding. In 14 villages across both districts, we identified high species richness, 19 species in total. Molecular findings indicated that An. minimus and An. baimaii were rare, while four other species (An. maculatus, An. pseudowillmori, An. jeyporiensis and An. nitidus) were abundant. Anopheles maculatus was highly prevalent in WKH (39% of light trap collections) and An. pseudowillmori in WJH (45%). Larvae of these four species were found in rice fields, suggesting that land cover change is influencing species composition change. Our results suggest that rice fields might be contributing to the observed abundance of An. maculatus and An. pseudowillmori, which could be playing a role in malaria transmission, either independently due to their high abundance, or in combination with An. baimaii and/or An. minimus.

The NADC34-like porcine reproductive and respiratory syndrome virus 2 (PRRSV-2) first emerged in China in 2017 and has the potential to become the dominant PRRSV strain in China. Here, a novel PRRSV-2, SCcd2020, was isolated from diseased piglets in Sichuan province, southwest China in 2020. The complete viral genome was determined and analyzed. An ORF5-based phylogenetic analysis showed that SCcd2020 clustered with NADC34-like strains, whereas the genome sequence clustered the isolate with NADC30-like viruses and it contains a discontinuous 131-aa deletion in NSP2 when compared to NADC30 strain. Notably, recombination analyses indicated that SCcd2020 is a multiple recombinant virus from NADC30-like, NADC34-like and JXA1-like strains, which is the first description of Chinese domestic HP-PRRSV involving the recombination event of an NADC34-like strain. Importantly, an animal challenge study in 4-week-old piglets showed that SCcd2020 causes high fever and severe hemorrhagic pneumonia with pulmonary consolidation and edema, and it has a high mortality rate (60%), which indicated that SCcd2020 is a highly pathogenic PRRSV strain. The study reports the emergence of a novel highly pathogenic NADC34-like recombinant strain, and it highlights the importance of monitoring newly emerging PRRSV strains in China.

Equine ocular setariasis is mainly caused by Setaria digitata, and the identification of this filarial nematode is based on morphology. However, morphological characterization alone is insufficient for the detection and differentiation of S. digitata from its congeners. In Thailand, the molecular detection of S. digitata is lacking and its genetic diversity is still unknown. This study aimed to phylogenetically characterize equine S. digitata from Thailand based on sequences derived from the mitochondrial cytochrome c oxidase subunit 1 (COI), the mitochondrial small subunit ribosomal DNA (12S rDNA), the nuclear internal transcribed spacer 1 (ITS1) and Wolbachia surface protein (wsp). Five samples of S. digitata were characterized, submitted to the NCBI database, and used for phylogenetic analysis as well as the assessment of similarity, entropy, and haplotype diversity. Phylogenetic analyses revealed that the S. digitata Thai strain was similar to S. digitata from China and Sri Lanka, with 99 to 100% similarity. The entropy and haplotype diversity indicated that the S. digitata Thai isolate was conserved and closely related to S. digitata worldwide. This is the first report on the molecular detection of equine ocular setariasis caused by S. digitata in Thailand.

Highly pathogenic avian influenza (HPAI) has caused widespread mortality in both wild and domestic birds in Europe during 2020-2022. Virus types H5N8 and H5N1 have dominated the epidemic. Isolated spill-over infections in mammals started to emerge as the epidemic continued. In autumn 2021, HPAI H5N1 caused a series of mass mortality events in farmed and released pheasants (Phasianus colchicus) in a restricted area in southern Finland. Later, in the same area, an otter (Lutra lutra), two red foxes (Vulpes vulpes) and a lynx (Lynx lynx) were found moribund or dead and infected with H5N1 HPAI virus. Phylogenetically, H5N1 strains from pheasants and mammals clustered together. Molecular analyses of the four mammalian virus strains revealed mutations in the PB2 gene segment (PB2-E627K and PB2-D701N) that are known to facilitate viral replication in mammals. This study revealed that avian influenza cases in mammals were spatially and temporally connected with avian mass mortalities suggesting increased infection pressure from birds to mammals.

In the recent survey, Japanese encephalitis (JE) is one of the most common mosquito-borne diseases, accounting for ∼30% of fatalities. The outbreaks of the JE virus (JEV) suggests that exhaustive study is essential for the prevention and management of the disease. The disease mainly spreads in humans and pigs by the vector: mosquito; as this is a major concern, this study had employed various bioinformatics tools to investigate the codon usage bias, evolutionary inference and selection pressure analysis of the Japanese encephalitis virus disease. The results indicated that the JE virus was biased and natural selection was the main factor shaping the codon usage that was determined and confirmed with the Nc, neutrality, PR2 plots and correlation analysis. The evolutionary analysis revealed that the virus had a substitution rate of 1.54 × 10-4 substitution/site/year and the tMRCA was found to be in 1723. The transmission of the virus in the map found transmissions mostly from China and transmitted across Asia and Africa. The selection pressure analysis employed three methods which had 969th codon site as diversifying site and had many purifying sites that shows the virus had evolved rapidly. The inferences from this study would aid people to employ this methodology on various diseases and also perform insilico studies in the field of vaccinology and immunoinformatics.

In early 2020, the emergence and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the human population quickly developed into a global pandemic. SARS-CoV-2 is the etiological agent of coronavirus disease 2019 (COVID-19) which has a broad range of respiratory illnesses. As the virus circulates, it acquires nucleotide changes. These mutations are potentially due to the inherent differences in the selection pressures within the human population compared to the original zoonotic reservoir of SARS-CoV-2 and formerly naïve humans. The acquired mutations will most likely be neutral, but some may have implications for viral transmission, disease severity, and resistance to therapies or vaccines. This is a follow-up study from our early report (Hartley et al. J Genet Genomics. 01202021;48(1):40-51) which detected a rare variant (nsp12, RdRp P323F) circulating within Nevada in mid 2020 at high frequency. The primary goals of the current study were to determine the phylogenetic relationship of the SARS-CoV-2 genomes within Nevada and to determine if there are any unusual variants within Nevada compared to the current database of SARS-CoV-2 sequences. Whole genome sequencing and analysis of SARS-CoV-2 from 425 positively identified nasopharyngeal/nasal swab specimens were performed from October 2020 to August 2021 to determine any variants that could result in potential escape from current therapeutics. Our analysis focused on nucleotide mutations that generated amino acid variations in the viral Spike (S) protein, Receptor binding domain (RBD), and the RNA-dependent RNA-polymerase (RdRp) complex. The data indicate that SARS-CoV-2 sequences from Nevada did not contain any unusual variants that had not been previously reported. Additionally, we did not detect the previously identified the RdRp P323F variant in any of the samples. This suggests that the rare variant we detected before was only able to circulate because of the stay-at-home orders and semi-isolation experience during the early months of the pandemic. IMPORTANCE: SARS-COV-2 continues to circulate in the human population. In this study, SARS-CoV-2 positive nasopharyngeal/nasal swab samples were used for whole genome sequencing to determine the phylogenetic relationship of SARS-CoV-2 sequences within Nevada from October 2020 to August 2021. The resulting data is being added to a continually growing database of SARS-CoV-2 sequences that will be important for understanding the transmission and evolution of the virus as it spreads around the globe.

Despite the high incidence of urogenital carcinoma (UGC) in California sea lions stranded along California, no UGC has been reported in other areas of their distribution; however, cell morphologies typical of premalignant states have been found. Risk factors for UGC include high of organochlorines and infection with a gammaherpesvirus, OtHV-1, but the importance of the bacteriome for epithelial status remains unknown. We characterized the genital bacteriome of adult female California sea lions along their distribution in the Gulf of California and examined whether the diversity and abundance of the bacteriome varied spatially, whether there were detectable differences in the bacteriome between healthy and altered epithelia, and whether the bacteriome was different in California sea lions infected with OtHV-1 or papillomavirus. We detected 2270 ASVs in the genital samples, of which 35 met the criteria for inclusion in the core bacteriome. Fusobacteriia and Clostridia were present in all samples, at high abundances, and Actinobacteria, Alphaproteobacteria, and Campylobacteria were also well-represented. Alpha diversity and abundance of the California sea lion genital bacteriome varied geographically. The abundance of bacterial ASVs varied depending on the genital epithelial status and inflammation, with differences driven by classes Fusobacteriia, Clostridia, Campylobacteria and Alphaproteobacteria. Alpha diversity and abundance were lowest in samples in which OtHV-1 was detected, and highest those with papillomavirus. Our study is the first investigation of how the bacteriome is related to epithelial status in a wild marine species prone to developing cancer.

BACKGROUND: Group A rotaviruses (RVA) are the primary pathogens of acute gastroenteritis. Currently, two live attenuated RVA vaccines, LLR and RotaTeq, have been introduced into mainland China but are not included in the national immunization program. Because of the unknown genetic evolution of group A rotavirus in an all-age population in Ningxia, China, we monitored the epidemiological characteristics and circulating genotypes of RVA as a reference for developing vaccine strategies. METHODS: We conducted seven years of consecutive surveillance of RVA based on stool samples from patients with acute gastroenteritis in sentinel hospitals in Ningxia, China, from 2015 to 2021. Reverse transcription quantitative polymerase chain reaction(RT-qPCR) was used to detect RVA in stool samples. Genotyping and phylogenetic analysis of VP7, VP4 and NSP4 genes were performed by reverse transcription-polymerase chain reaction(RT-PCR) and nucleotide sequence determination. RESULTS: RVA was detected in 16.58% (1436/8662) of 8662 stool samples. The positive rates were 7.17% (201/2805) and 21.09% (1235/5857) in adults and children, respectively. The most affected age group was infants and children aged 12-23 months, with a positive rate of 29.53% (p < 0.05). A significant winter/spring seasonality was observed. 23.29% positive rate in 2020 was the highest in 7 years (p < 0.05). The region with the highest positive rate in the adult group was Yinchuan, and the children's group was Guyuan. A total of 9 genotype combinations were found to be distributed in Ningxia. The dominant genotype combinations in this region gradually changed from G9P[8]-E1, G3P[8]-E1, G1P[8]-E1 to G9P[8]-E1, G9P[8]-E2, and G3P[8]-E2 during these seven years. Rare strains (e.g., G9P[4]-E1, G3P[9]-E3 and G1P[8]-E2) were occasionally detected during the study. CONCLUSIONS: During the study period, changes in the significant RVA circulating genotype combinations and the emergence of reassortment strains were observed, particularly the emergence and prevalence of G9P[8]-E2, G3P[8]-E2 reassortants in the region. These results indicate the importance of continuous monitoring of the molecular evolution and recombination characteristics of RVA, and should not be limited to G/P genotyping but should consider multi-gene fragment co-analysis and whole genome sequencing.

South-East Asia (SEA) and South Asia (SA) are two important geographic regions with the most severe enzootic rabies in the world. In these regions, phylogenetic analysis of rabies virus (RABV) has been conducted only at a country level; the results obtained from different countries are scattered and unequal, with a non-uniform system to name RABV genotypes. Therefore, it is difficult to undertake origin-tracking and compare inter-country RABV evolution and transmission. To avoid the confusion in understanding and to generate a panoramic picture of RABV genetic diversity, distribution, and transmission in SEA and SA, the present study conducted a systematic phylogenetic analysis by combining all sequences representing 2368 RABV strains submitted to GenBank by 14 rabies endemic SEA and SA countries. The results showed that RABVs circulating in two regions were classified into four major clades and many subclades: the Asia clade is circulating only in SEA, the Indian subcontinent, and Arctic-like clades only in SA, while the Cosmopolitan clade has been detected in both regions. The results also showed a wide range of hosts were infected by divergent RABV subclades, with dogs being the major transmission source. However, wildlife rabies was also found to be an important issue with 6 wild carnivore species identified as potential sources of spillover risk for sylvatic rabies to humans, domestic animals, and other wild animals. Current findings indicate that the two regions have separate virus clades circulating thus indicating the absence of cross-transmission between the regions. The study emphasizes the importance of phylogenetic analysis in the regions using uniform genotyping and naming systems for rabies surveillance, to coordinate actions of member countries to eliminate dog-mediated human rabies by 2030.

The local diversity and population structure of malaria parasites vary across different regions of the world, reflecting variations in transmission intensity, host immunity, and vector species. This study aimed to use amplicon sequencing to investigate the genotypic patterns and population structure of P. vivax isolates from a highly endemic province of Thailand in recent years. Amplicon deep sequencing was performed on 70 samples for the 42-kDa region of pvmsp1 and domain II of pvdbp. Unique haplotypes were identified and a network constructed to illustrate genetic relatedness in northwestern Thailand. Based on this dataset of 70 samples collected between 2015 and 2021, 16 and 40 unique haplotypes were identified in pvdbpII and pvmsp142kDa, respectively. Nucleotide diversity was higher in pvmsp142kDa than in pvdbpII (π = 0.027 and 0.012), as was haplotype diversity (Hd = 0.962 and 0.849). pvmsp142kDa also showed a higher recombination rate and higher levels of genetic differentiation (Fst) in northwestern Thailand versus other regions (0.2761-0.4881). These data together suggested that the genetic diversity of P. vivax in northwestern Thailand at these two studied loci evolved under a balancing selection, most likely host immunity. The lower genetic diversity of pvdbpII may reflect its stronger functional constrain. In addition, despite the balancing selection, a decrease in genetic diversity was observed. Hd of pvdbpII decreased from 0.874 in 2015-2016 to 0.778 in 2018-2021; π of pvmsp142kDa decreased from 0.030 to 0.022 over the same period. Thus, the control activities must have had a strong impact on the parasite population size. The findings from this study provide an understanding of P. vivax population structure and the evolutionary force on vaccine candidates. They also established a new baseline for tracking future changes in P. vivax diversity in the most malarious area of Thailand.

The wide dissemination of plasmids carrying antibiotic resistance determinants among bacteria is a severe threat to global public health. Here, we characterized an extensively drug-resistant (XDR) Klebsiella pneumoniae NTU107224 by whole genome sequencing (WGS) in combination with phenotypic tests. Broth dilution method was used to determine the minimal inhibitory concentrations (MICs) of NTU107224 to 24 antibiotics. The whole genome sequence of NTU107224 was determined by Nanopore/Illumina hybrid genome sequencing. Conjugation assay was performed to determine the transferability of plasmids in NTU107224 to recipient K. pneumoniae 1706. Larvae infection model was used to determine the effect(s) of conjugative plasmid pNTU107224-1 on bacterial virulence. Among the 24 antibiotics tested, XDR K. pneumoniae NTU107224 had low MICs only for amikacin (≤1 μg/mL), polymyxin B (0.25 μg/mL), colistin (0.25 μg/mL), eravacycline (0.25 μg/mL), cefepime/zidebactam (1 μg/mL), omadacycline (4 μg/mL), and tigecycline (0.5 μg/mL). Whole genome sequencing showed that the closed NTU107224 genome comprises a 5,076,795-bp chromosome, a 301,404-bp plasmid named pNTU107224-1, and a 78,479-bp plasmid named pNTU107224-2. IncHI1B plasmid pNTU107224-1 contained three class 1 integrons accumulated various antimicrobial resistance genes (including carbapenemase genes blaVIM-1, blaIMP-23, and truncated blaOXA-256) and the blast results suggested the dissemination of IncHI1B plasmids in China. By day 7 after infection, larvae infected with K. pneumoniae 1706 and transconjugant had 70% and 15% survival rates, respectively. We found that the conjugative plasmid pNTU107224-1 is closely related to IncHI1B plasmids disseminated in China and contributes to the virulence and antibiotic resistance of pathogens.

Infection with the hepatitis C virus (HCV) remains a considerable public health concern in the Middle East and North Africa (MENA). The objectives of this study were to analyze the HCV genotype (GT) and sub-genotype (SGT) distribution in the MENA region and to assess the temporal change in the number of sequences within the MENA region. All HCV molecular sequences collected in the MENA region had been retrieved from GenBank as of 1 August 2022. The number of HCV sequences retrieved was 6740 representing sequences from a total of 17 MENA countries with a majority from Iran (n = 1969, 29.2%), Egypt (n = 1591, 23.6%), Tunisia (n = 1305, 19.4%) and Saudi Arabia (n = 1085, 16.1%). The determination of GT/SGT was based on the NCBI genotyping and Blast tool. Genotype 1 (GT1) dominated infections in the MENA (n = 2777, 41.2%), followed by GT4 (n = 2566, 39.0%). Additionally, SGT4a (1515/6393, 23.7%) was the most common SGT in the MENA, and SGT4a was dominant in Egypt and Saudi Arabia, followed by SGT1b (n = 1308, 20.5%), which was dominant in Morocco and Tunisia, while SGT1a (n = 1275, 19.9%) was common in Iran, Iraq and Palestine. Furthermore, significant temporal increase in the number of HCV MENA sequences was observed. On the SGT level, specific patterns of HCV genetic diversity were seen in the MENA region, with the most common SGT being 4a, in addition to increasing the availability of HCV sequences in the MENA region.

In the worst future pandemic, effective vaccines and medicines could be unavailable for a long time. In such circumstances, it is necessary to evaluate whether a periodic screening can protect isolated communities and critical facilities and avoid a complete shutdown. In this study, we introduced an epidemiological model that included the essential parameters of infection transmission and screening. With varying parameters, we studied the dynamics of viral infection in the isolated communities. In the scenario with a periodic infection screening once per 3 days and a viral basic reproduction number 3.0, >85% of the infection waves have a duration <7 days and the infection size in each of the waves is generally <4 individuals when the efficiency of infection discovery is 0.9 in the screening. When the period of screening was elongated to once per 7 days, the cases of infection dramatically increased to 5 folds of that mentioned previously. Further, with a weak discovery efficiency of 0.7 and the aforementioned low screening frequency, the spread of infection would be out of control. Our study suggests that frequent periodic screening is capable of controlling a future epidemic in isolated communities without other measures.

In October 2021, the world's first malaria vaccine RTS,S was endorsed by WHO for broad use in children, despite its low efficacy. This study examined polyclonal infections and the associations of parasite genetic variations with binding affinity to human leukocyte antigen (HLA). Multiplicity of infection was determined by amplicon deep sequencing of PfMSP1. Genetic variations in PfCSP were examined across 88 samples from Ghana and analyzed together with 1655 PfCSP sequences from other African and non-African isolates. Binding interactions of PfCSP peptide variants and HLA were predicted using NetChop and HADDOCK. High polyclonality was detected among infections, with each infection harboring multiple non-3D7 PfCSP variants. Twenty-seven PfCSP haplotypes were detected in the Ghanaian samples, and they broadly represented PfCSP diversity across Africa. The number of genetic differences between 3D7 and non-3D7 PfCSP variants does not influence binding to HLA. However, CSP peptide length after proteolytic degradation significantly affects its molecular weight and binding affinity to HLA. Despite the high diversity of HLA, the majority of the HLAI and II alleles interacted/bound with all Ghana CSP peptides. Multiple non-3D7 strains among P. falciparum infections could impact the effectiveness of RTS,S. Longer peptides of the Th2R/Th3R CSP regions should be considered in future versions of RTS,S.

The Insertion Sequence 711 (IS711) is linked to the Brucella genus. Mapping the genomic distribution of IS711 can help understand this insertion element's biological and evolutionary role. This work aimed to delineate the genomic distribution of the IS711 element and to study its association with Brucella evolution. A total of 124 genomes representing 9 Brucella species were searched using BLASTn sequence alignment tool to identify complete and truncated copies of IS711. Based on the genomic context, each IS711 locus was assigned a code using the initial letters of its neighboring genes. Various tools were used to annotate the neighboring genes and determine the shared synteny around orthologous IS711 loci. The tool Islandviewer 4 was used to scan for genomic islands. The Codon Tree method was used to build phylogenetic trees of B. melitensis, B. abortus, and B. suis genomes. The phylogenetic trees of the three species were analyzed, taking into account the genomic distribution patterns of IS711. The result of IS711 frequency analysis showed a relatively conserved number of copies/genome for the different species and for some biovars. The analysis showed that Brucella species with a relatively low IS711 copy number (4-8 copies/genome) are linked to domestic animals as primary hosts and have potential for zoonotic transmission. However, species with a relatively higher copy number (12-30 copies/genome) are less zoonotic and tend to be linked with wild animals as primary hosts. Analyzing the genomic distribution map of IS711 loci showed several unique patterns of IS711 distribution that are correlated with the evolution of Brucella species and biovars. The results also showed that 46.2% of the conserved IS711 elements are located within genomic islands. Based on our results and previous data, we postulate a model explaining the IS711 role in Brucella evolution. We assume that during the transition from a free-living to an intracellular lifestyle, a descendant of the Brucella genus had acquired a progenitor sequence of the IS711. Subsequently, a burst in IS711 transposition occurred. This parasitic expansion can be deleterious and has to be counteracted by evolutionary forces to prevent lineage extension and to promote adaptation to host. Similar to other plasmid-free pathogenic α-Proteobacteria bacteria, the balance of expansion and reduction of insertion elements could be one of the mechanisms to control genome reduction and streamlining. We hypothesize that the IS711-mediated genomic changes and other small sequence nucleotide changes in specific orthologous genes could significantly contribute to Brucella's evolution and adaptation to different animal hosts.

Recent reports on identification of canine coronavirus (CCoV) in humans have emphasized the urgency to strengthen surveillance of animal CoVs. The fact that recombinations between CCoV with feline, porcine CoVs brought about new types of CoVs indicated that more attention should be paid to domestic animals like dogs, cats and pigs, and the CoVs they carried. However, there are about ten kinds of CoVs that infect above animals, and thus representative CoVs with zoonotic potentials were considered in this study. Multiplex RT-PCR against CCoV, Feline coronavirus (FCoV), porcine deltacoronavirus and porcine acute diarrhea syndrome coronavirus was developed to investigate the prevalence of CoVs from domestic dogs in Chengdu, Southwest China. Samples from a total of 117 dogs were collected from a veterinary hospital, and only CCoV (34.2%, 40/117) was detected. Therefore, this study focused on CCoV and its characteristics of S, E, M, N and ORF3abc genes. Compared with CoVs that are capable of infecting humans, CCoV strains showed highest nucleotide identity with the novel canine-feline recombinant detected from humans (CCoV-Hupn-2018). Phylogenetic analysis based on S gene, CCoV strains were not only clustered with CCoV-II strains, but also closely related to FCoV-II strains ZJU1617 and SMU-CD59/2018. As for assembled ORF3abc, E, M, N sequences, CCoV strains had the closest relationship with CCoV-II (B203_GZ_2019, B135_JS_2018 and JS2103). What's more, specific amino acid variations were found, especially in S and N proteins, and some mutations were consistent with FCoV, TGEV strains. Altogether, this study provided a novel insight into the identification, diversification and evolution of CoVs from domestic dogs. It is of top priority to recognize zoonotic potential of CoVs, and continued comprehensive surveillance will help better understand the emergence, spreading, and ecology of animal CoVs.

Data obtained from new sequencing technologies are evolving rapidly, leading to the development of specific bioinformatic tools, pipelines and softwares. Several algorithms and tools are today available allowing a better identification and description of Mycobacterium tuberculosis complex (MTBC) isolates worldwide. Our approach consists in applying existing methods to analyze DNA sequencing data (from FASTA or FASTQ files), and tentatively extract meaningful information that would facilitate identification as well as a better understanding and management of MTBC isolates (taking into account whole genome sequencing and classical genotyping data). The aim of this study is to propose a pipeline analysis allowing to potentially simplify MTBC data analysis by providing different ways to interpret genomic or genotyping information based on existing tools. Furthermore, we propose a "reconciledTB" list making a link with results directly obtained from whole genome sequencing (WGS) data and results obtained from classical genotyping analysis (data inferred from SpoTyping and MIRUReader). Data visualization graphics and trees generated provide additional elements to better understand and confer associations among information overlap analyses. Additionally, comparison between data entered in an international genotyping database (SITVITEXTEND) and ensuing data obtained from the pipeline not only provide meaningful information, but further suggest that simpiTB could also be suitable for new data integration in specific TB genotyping databases.

Gonorrhea is an urgent antimicrobial resistance threat and its therapeutic options are continuously getting restricted. Moreover, no vaccine has been approved against it so far. Hence, the present study aimed to introduce novel immunogenic and drug targets against antibiotic-resistant Neisseria gonorrhoeae strains. In the first step, the core proteins of 79 complete genomes of N. gonorrhoeae were retrieved. Next, the surface-exposed proteins were evaluated from different aspects such as antigenicity, allergenicity, conservancy, and B-cell and T-cell epitopes to introduce promising immunogenic candidates. Then, the interactions with human Toll-like receptors (TLR-1, 2, and 4), and immunoreactivity to elicit humoral and cellular immune responses were simulated. On the other hand, to identify novel broad-spectrum drug targets, the cytoplasmic and essential proteins were detected. Then, the N. gonorrhoeae metabolome-specific proteins were compared to the drug targets of the DrugBank, and novel drug targets were retrieved. Finally, the protein data bank (PDB) file availability and prevalence among the ESKAPE group and common sexually transmitted infection (STI) agents were assessed. Our analyses resulted in the recognition of ten novel and putative immunogenic targets including murein transglycosylase A, PBP1A, Opa, NlpD, Azurin, MtrE, RmpM, LptD, NspA, and TamA. Moreover, four potential and broad-spectrum drug targets were identified including UMP kinase, GlyQ, HU family DNA-binding protein, and IF-1. Some of the shortlisted immunogenic and drug targets have confirmed roles in adhesion, immune evasion, and antibiotic resistance that can induce bactericidal antibodies. Other immunogenic and drug targets might be associated with the virulence of N. gonorrhoeae as well. Thus, further experimental studies and site-directed mutations are recommended to investigate the role of potential vaccine and drug targets in the pathogenesis of N. gonorrhoeae. It seems that the efforts for proposing novel vaccines and drug targets appear to be paving the way for a prevention-treatment strategy against this bacterium. Additionally, a combination of bactericidal monoclonal antibodies and antibiotics is a promising approach to curing N. gonorrhoeae.

BACKGROUND: Emergence and predominance of hepatitis B virus (HBV) variants carrying S gene mutations frequently occur in HBV-infected individuals. Here, coexistent serum anti-HBsAg antibody (HBsAb) and HBV surface antigen (HBsAg) were detected in a chronic HBV patient. The patient's HBsAg proteins possessed amino acid substitutions sK122R and sV96A. We reported this case and conducted relevant studies to investigate differences in expression levels and antibody neutralization of HBsAg proteins bearing sK122R and sV96A amino acid substitutions to explore causes of antigen-antibody coexistence in a chronic hepatitis B patient. STUDY DESIGN: We first sequenced the S gene from HBV present within the patient's serum. Based on the S gene sequence, we cloned wild-type and mutated S gene sequences via site-directed mutagenesis to construct expression plasmids pJW4303-WT (wild-type), pJW4303-sV96A, pJW4303-sK122R, and pJW4303-sV96A-sK122R. Plasmids were transfected into HEK 293 T cells then culture supernatants and cells were collected. Collected cells and supernatants were next subjected to a series of quantitative and functional tests to assess expression and neutralization characteristics of wild-type and mutant HBsAg proteins. RESULTS: Based on quantification of HBsAg expression in cells transfected with the four plasmids, HBsAg-sK122R-sV96A was more intracellularly retained and less secreted than HBsAg-sV96A single-mutant protein and WT. Neutralization ability of serum from chronic HBV patient against culture supernatants containing recombinant HBsAg proteins were ranked from highest to lowest as HBsAg-sV96A, HBsAg-sV96A-sK122R, and HBsAg-sK122R. However, no significant differences of neutralization efficiency by high-potency antibodies from HBV-vaccinees against these three mutant proteins were observed. CONCLUSIONS: The levels of HBsAg proteins with amino acid substitutions sV96A-sK122R were greatly reduced in culture supernatants but were apparently increased in the intracellular fraction. This may account for the higher levels of HBV replication in patients. HBsAg neutralization by HBsAb in this patient may have been compromised by the HBsAg sK122R amino acid substitution, suggesting that antibodies produced by the patient had lost their HBV-neutralizing effect.

Stomoxys flies (Diptera: Muscidae) are hematophagous ectoparasites of medical and veterinary importance. In this study, three Stomoxys species, i.e. S. bengalensis, S. calcitrans, and S. sitiens, were collected from three provinces in Central Thailand with the aim of estimating the genetic divergence between species, for species identification, as well as within species, for a genetic diversity study based on the cytochrome c oxidase subunit I (COI) gene. Our results showed that the average intraspecific genetic divergences of Stomoxys flies ranged from 0.11% in S. sitiens to 0.98% in S. calcitrans, whereas the average interspecific genetic divergences ranged from 5.24% between S. sitiens and S. bengalensis to 6.69% between S. calcitrans and S. bengalensis. In addition, there was no overlap between the intraspecific and interspecific genetic divergences. The COI sequence analysis revealed a high haplotype diversity and low nucleotide diversity, reflecting a rapid population expansion after past bottlenecks. Moreover, there was no significant difference (P > 0.05) in the pairwise population differentiation (Fst) among Stomoxys flies in Central Thailand, because of the lack of natural barriers, thus allowing genetic exchange between them. The monitoring of the haplotype network revealed that two lineages of S. calcitrans in Central Thailand were distributed in all study areas, including the Nakhon Pathom, Pathum Thani, and Saraburi Provinces. These findings may improve our understanding of the genetic patterns of these three Stomoxys flies, as well as the underlying biological mechanisms, which is knowledge that can be used for further effective control of these flies.

OBJECT: High-risk human papillomavirus (HR-HPV) is a main reason for cervical cancer. The long control region (LCR) of the genome plays a variety of roles in the transcription of the virus. METHODS: LCR sequences were amplified by polymerase chain reaction (PCR) and confirmed by DNA sequencing. MEGA 11.0 software and NCBI blast were used to analyze the sequences and construct the Neighbor-Joining tree. In addition, the JASPAR database was used to predict the potential transcription factor binding sites (TFBS). RESULTS: For HPV-52 LCR, 68 single nucleotide polymorphisms (SNPs), 8 deletions, and 1 insertion were found, 17 of which were novel variations. Most of the variants were clustered in B2 sub-lineage (96.22%). For HPV-58 LCR, 25.43% of samples were prototype. 49 SNPs, 2 deletions, and 1 insertion were observed in the remaining samples. A1 sub-lineage was the most frequent (64.16%). For HPV-16 LCR, 75 SNPs and 2 deletions were identified, 13 of which were newly identified. A total of 55.68% of the variants were distributed in A4 sub-lineage. The JASPAR results suggested that multiple variations occurred in TFBSs, which might affect the function of transcription factors. CONCLUSIONS: This study provides experimental data for further studies on the epidemiology and biological function of LCR. Various LCR mutational data may prove useful for exploring the carcinogenic mechanism of HPV.

BACKGROUND: Aedes albopictus is an important vector of many mosquito-borne viral diseases, including dengue fever and Zika. In recent years, it has spread and colonized tropical, subtropical and temperate regions worldwide. Monitoring of Ae. albopictus population dynamics is an important tool for early warning of mosquito-borne infections. Because the genetic diversity and genetic structure of natural populations are the genetic bases of population dynamics, studies of population genetics can reveal the origin, differentiation and dispersal characteristics of Ae. albopictus populations. Then, their evolutionary potential and environmental adaptability can be analyzed, providing a theoretical basis for the formulation of accurate Ae. albopictus surveillance and integrated control programs. METHODS: In 2018, 552 Ae. albopictus larvae were collected during an invasive mosquito species surveillance project in China's Yangtze River Basin. Morphological analysis was performed to assign the adult mosquitoes to species, and then the genetic marker ITS2 was amplified and sequenced. RESULTS: There were 179 haplotypes among 552 ITS2 sequences. In total, 155/179 (86.59%) haplotypes were specific to individual populations, and 24/179 (13.41%) haplotypes were shared by populations. Hap4 (126), Hap7 (43), and Hap16 (34) were the most numerous haplotypes and the most widely distributed. The overall Hd was 0.928, π was 0.031, the mean nucleotide difference number (K) was 7.255, and the number of segregating sites was 169. TCS network maps mainly showed a single star-like scattered distribution. According to geographical location, there were no obvious haplotype groups, and the haplotypes were intricately connected. The genetic diversity of Ae. albopictus populations in the Yangtze River Basin was high. The molecular variance observed in the populations of Ae. albopictus mainly occurred among individuals within populations, accounting for 98.79% of the total, while that among populations accounted for only 1.21% of the total. Only the populations of Ae. albopictus in the Chongqing and Sichuan regions showed a moderate degree of population genetic differentiation, while genetic differentiation between the other regions were small, gene exchange was very common, and genetic differentiation within populations was minimal. CONCLUSIONS: According to this study, the genetic diversity of Ae. albopictus populations in the Yangtze River Basin is high, the genetic differentiation among populations is small, and gene exchange is common. In addition, frequent interregional exchange exacerbates the abnormal spread of vectors. This study highlighted the potential spread route of the vector Ae. albopictus in the Yangtze River Basin. There are three potential dispersal routes for Ae. albopictus populations in the Yangtze River Basin. The findings could be helpful for effective surveillance and early warning of Ae. albopictus vectors.

BACKGROUND: The COVID-19 pandemic has put many medical systems on the verge of collapse in the last two years. Virus mutation was one of the important factors affecting the COVID-19 infection severity and hospitalizations. Although over ten thousand SARS-CoV-2 mutations being reported since the beginning of the COVID-19 pandemic, only a small percentage of mutations are likely to affect the virus phenotype and change its severity. Finding out which amino acids have the greatest impact on COVID-19 hospitalization rate is an important research question. METHODS: This observational study used the COVID-19 case hospitalization ratio (CHR) to represent the virus severity related with hospitalization. The database is based on the daily state-level epidemiological and genomic sequential data in the United States from the Alpha wave to the first Omicron wave. The critical amino acids that mostly affected the CHR were determined by using four types of models including extreme gradient boosting decision trees (XGBoost), artificial neural networks (ANNs), logistic regression and Lasso regression models. RESULTS: The XGBoost, ANN, logistic regression, and Lasso regression models all produce excellent results (mean square error for all state-level models does not exceed 0.0008 using the testing dataset). Based on the rank of importance of all covariates, the critical amino acids most affecting the CHR were identified, including T19, L24, P25, P26, A27, A67, H69, V70, T95, G142, V143, Y145, E156, F157, N211, L212, V213, R214, D215, G339, R346, S373, L452, S477, T478, E484, N501, A570, P681, and T716. CONCLUSION: This study identified critical amino acids that are most likely to affect the hospitalization rate, allowing public health workers to monitor these highly risky amino acids and raise an alarm immediately when more severe mutations occur. Furthermore, the methodology and results may be extended to other regions.

Eliminating schistosomiasis as a public health problem by 2030 requires a better understanding of the disease transmission, especially the asymmetric distribution of worm burden in individuals living and sharing the same environment. It is in this light that this study was designed to identify human genetic determinants associated with high burden of S. mansoni and also with the plasma concentrations of IgE and four cytokines in children from two schistosomiasis endemic areas of Cameroon. In school-aged children of schistosomiasis endemic areas of Makenene and Nom-Kandi of Cameroon, S. mansoni infections and their infection intensities were evaluated in urine and stool samples using respectively the Point-of-care Circulating Cathodic Antigen test (POC-CCA) and the Kato Katz (KK) test. Thereafter, blood samples were collected in children harbouring high burden of schistosome infections as well as in their parents and siblings. DNA extracts and plasma were obtained from blood. Polymorphisms at 14 loci of five genes were assessed using PCR-restriction fragment length polymorphism and amplification-refractory mutation system. The ELISA test enabled to determine the plasma concentrations of IgE, IL-13, IL-10, IL-4 and IFN-γ. The prevalence of S. mansoni infections was significantly higher (P < 0.0001 for POC-CCA; P = 0.001 for KK) in Makenene (48.6% for POC-CCA and 7.9% for KK) compared to Nom-Kandi (31% for POC-CCA and 4.3% for KK). The infection intensities were also higher (P < 0.0001 for POC-CCA; P = 0.001 for KK) in children from Makenene than those from Nom-Kandi. The allele C of SNP rs3024974 of STAT6 was associated with an increased risk of bearing high burden of S. mansoni both in the additive (p = 0.009) and recessive model (p = 0.01) while the allele C of SNP rs1800871 of IL10 was protective (p = 0.0009) against high burden of S. mansoni. The alleles A of SNP rs2069739 of IL13 and G of SNP rs2243283 of IL4 were associated with an increased risk of having low plasma concentrations of IL-13 (P = 0.04) and IL-10 (P = 0.04), respectively. This study showed that host genetic polymorphisms may influence the outcome (high or low worm burden) of S. mansoni infections and also the plasma concentrations of some cytokines.

A decreased susceptibility of Fulani populations to malaria infections has been shown in Africa. A previous longitudinal cohort study conducted in the Atacora region of northern Benin showed a high merozoite-phagocytosis capacity in young Fulani. Here, we explored the combined polymorphisms in the constant region of the IgG3 heavy chain (presence/absence of the G3m6 allotype) and in Fc gamma receptors (FcγRs) as potentially involved in the natural protection against malaria of young Fulani in Benin. An active malaria follow-up was conducted among individuals from Fulani, Bariba, Otamari and Gando ethnic groups living in sympatry in Atacora, over the full malaria transmission season. FcγRIIA 131R/H (rs1801274), FcγRIIC C/T (rs3933769) and FcγRIIIA 176F/V (rs396991) were determined using the TaqMan method; FcγRIIIB NA1/NA2 was assessed by polymerase chain reaction (PCR) using allele-specific primers and G3m6 using allotype by PCR-RFLP. Individual carriage of G3m6 (+) was associated with an increased risk of Pf malaria infection (logistic multivariate regression model (lmrm), OR = 2.25, 95% CI = 1.06;4.74, P = 0.034). Combined haplotype G3m6 (+) - FcγRIIA 131H - FcγRIIC T - FcγRIIIA 176F - FcγRIIIB NA2 was also associated with an increased risk of Pf malaria infection (lmrm, OR = 13.01, 95% CI = 1.69;99.76, P = 0.014). G3m6 (-), FcγRIIA 131R and FcγRIIIB NA1 were more prevalent in young Fulani (P = 0.002, P < 0.001 and P = 0.049, respectively), while no Fulani presented the combined G3m6 (+) - FcγRIIA 131H - FcγRIIC T - FcγRIIIA 176F - FcγRIIIB NA2 haplotype that was carried by a majority of infected children. Our results highlight the combined factors G3m6 - FcγR as potentially involved in the merozoite-phagocytosis capacity and in the natural protection of young Fulani individuals against P. falciparum malaria in Benin.

Studies on antimicrobial resistance (AMR) profiles and epidemiological affirmation for AMR transmission are limited in fisheries and aquaculture settings. Since 2015, based on Global Action Plan on AMR by World Health Organization (WHO) and World Organization for Animal Health (OIE), several initiatives have been under taken to enhance the knowledge, skills and capacity to establish AMR trends through surveillance and strengthening of epidemiological evidence. The focus of this study was to determine the prevalence of antimicrobial resistance (AMR), its resistance profiles and molecular characterization with respect to phylogroups, antimicrobial resistance genes (ARGs), virulence genes (VGs), quaternary ammonium compounds resistance (QAC) genes and plasmid typing in retail market fishes. Pulse field gel electrophoresis (PFGE) to understand the genetic lineage of the two most important Enterobacteriaceae members, E. coli and Klebsiella sp. was performed. 94 fish samples were collected from three different sites viz., Silagrant (S1), Garchuk (S2) and North Guwahati Town Committee (NGTC) Region (S3) in Guwahati, Assam. Out of the 113 microbial isolates from the fish samples, 45 (39.82%) were E. coli; 23 (20.35%) belonged to Klebsiella genus. Among E. coli, 48.88% (n = 22) of the isolates were alerted by the BD Phoenix M50 instrument as ESBL, 15.55% (n = 7) as PCP and 35.55% (n = 16) as non-ESBL. E. coli (39.82%) was the most prevalent pathogen among the Enterobacteriaceae members screened and showed resistance to ampicillin (69%) followed by cefazoline (64%), cefotaxime (49%) and piperacillin (49%). In the present study, 66.66% of E. coli and 30.43% of Klebsiella sp. were categorized as multi drug resistance (MDR) bacteria. CTX-M-gp-1, with CTX-M-15 variant (47%), was the most widely circulating beta-lactamase gene, while other ESBL genes blaTEM (7%), blaSHV (2%) and blaOXA-1-like (2%) were also identified in E. coli. Out of the 23 isolates of Klebsiella, 14(60.86%) were ampicillin (AM)-resistant (11(47.82%) K. oxytoca, 3(13.04%) K. aerogenes), whereas 8(34.78%) isolates of K. oxytoca showed intermediate resistance to AM. All Klebsiella isolates were susceptible to AN, SCP, MEM and TZP, although two K. aerogenes were resistant to imipenem. DHA and LAT genes were detected, respectively, in 7(16%) and 1(2%) of the E. coli strains while a single K. oxytoca (4.34%) isolate carried MOX, DHA and blaCMY-2 genes. The fluoroquinolone resistance genes identified in E. coli included qnrB (71%), qnrS (84%), oqxB (73%) and aac(6)-Ib-cr (27%); however, in Klebsiella, these genes, respectively, had a prevalence of 87%, 26%, 74% and 9%. The E. coli isolates belonged to phylogroup A(47%), B1(33%) and D(14%). All of the 22(100%) ESBL E. coli had chromosome-mediated disinfectant resistance genes viz., ydgE, ydgF, sugE(c), mdfA while 82% of ESBL E. coli had emrE. Among the non-ESBL E. coli isolates, 87% of them showed the presence of ydgE, ydgF and sugE(c) genes, while 78% of the isolates had mdfA and 39% had emrE genes respectively. 59% of the ESBL and 26% of the non-ESBL E. coli had showed the presence of qacEΔ1. The sugE(p) was present in 27% of the ESBL-producing E. coli and in 9% of non-ESBL isolates. Out of the 3 ESBL-producing Klebsiella isolates, 2(66.66%) K. oxytoca isolates were found harboring plasmid-mediated qacEΔ1 gene while one (33.33%) K. oxytoca isolate had sugE(p) gene. IncFI was the most prevalent plasmid type detected in the isolates studied, with A/C (18%), P (14%), X, Y (9% each) and I1-Iγ (14%, 4%). 50% (n = 11) of the ESBL and 17% (n = 4) of the non-ESBL E. coli isolates harboured IncFIB and 45% (n = 10) ESBL and one (4.34%) non-ESBL E. coli isolates harboured IncFIA. Dominance of E. coli over other Enterobacterales and diverse phylogenetic profiles of E. coli and Klebsiella sp. suggests the possibility of contamination and this may be due to compromised hygienic practices along the supply chain and contamination of aquatic ecosystem. Continuous surveillance in domestic markets must be a priority in addressing antimicrobial resistance in fishery settings and to identify any unwarranted epidemic clones of E. coli and Klebsiella that can challenge public health sector.