SARS-CoV-2 genomes from Saudi Arabia implicate nucleocapsid mutations in host response and increased viral load.

Monitoring SARS-CoV-2 spread and evolution through genome sequencing is essential in handling the COVID-19 pandemic. Here, we sequenced 892 SARS-CoV-2 genomes collected from patients in Saudi Arabia from March to August 2020. We show that two consecutive mutations (R203K/G204R) in the nucleocapsid (N) protein are associated with higher viral loads in COVID-19 patients. Our comparative biochemical analysis reveals that the mutant N protein displays enhanced viral RNA binding and differential interaction with key host proteins. We found increased interaction of GSK3A kinase simultaneously with hyper-phosphorylation of the adjacent serine site (S206) in the mutant N protein. Furthermore, the host cell transcriptome analysis suggests that the mutant N protein produces dysregulated interferon response genes. Here, we provide crucial information in linking the R203K/G204R mutations in the N protein to modulations of host-virus interactions and underline the potential of the nucleocapsid protein as a drug target during infection.

Genome-wide transposon mutagenesis analysis of Burkholderia pseudomallei reveals essential genes for in vitro and in vivo survival.

Introduction: Burkholderia pseudomallei, a soil-dwelling microbe that infects humans and animals is the cause of the fatal disease melioidosis. The molecular mechanisms that underlie B. pseudomallei's versatility to survive within a broad range of environments are still not well defined. Methods: We used the genome-wide screening tool TraDIS (Transposon Directed Insertion-site Sequencing) to identify B. pseudomallei essential genes. Transposon-flanking regions were sequenced and gene essentiality was assessed based on the frequency of transposon insertions within each gene. Transposon mutants were grown in LB and M9 minimal medium to determine conditionally essential genes required for growth under laboratory conditions. The Caenorhabditis elegans infection model was used to assess genes associated with in vivo B. pseudomallei survival. Transposon mutants were fed to the worms, recovered from worm intestines, and sequenced. Two selected mutants were constructed and evaluated for the bacteria's ability to survive and proliferate in the nematode intestinal lumen. Results: Approximately 500,000 transposon-insertion mutants of B. pseudomallei strain R15 were generated. A total of 848,811 unique transposon insertion sites were identified in the B. pseudomallei R15 genome and 492 genes carrying low insertion frequencies were predicted to be essential. A total of 96 genes specifically required to support growth under nutrient-depleted conditions were identified. Genes most likely to be involved in B. pseudomallei survival and adaptation in the C. elegans intestinal lumen, were identified. When compared to wild type B. pseudomallei, a Tn5 mutant of bpsl2988 exhibited reduced survival in the worm intestine, was attenuated in C. elegans killing and showed decreased colonization in the organs of infected mice. Discussion: The B. pseudomallei conditional essential proteins should provide further insights into the bacteria's niche adaptation, pathogenesis, and virulence.

A genetic barcode of SARS-CoV-2 for monitoring global distribution of different clades during the COVID-19 pandemic.

OBJECTIVE: The SARS-CoV-2 pathogen has established endemicity in humans. This necessitates the development of rapid genetic surveillance methodologies to serve as an adjunct with existing comprehensive, albeit though slower, genome sequencing-driven approaches. METHODS: A total of 21,789 complete genomes were downloaded from GISAID on May 28, 2020 for analyses. We have defined the major clades and subclades of circulating SARS-CoV-2 genomes. A rapid sequencing-based genotyping protocol was developed and tested on SARS-CoV-2-positive RNA samples by next-generation sequencing. RESULTS: We describe 11 major mutations which defined five major clades (G614, S84, V251, I378 and D392) of globally circulating viral populations. The clades can specifically identify using an 11-nucleotide genetic barcode. An analysis of amino acid variation in SARS-CoV-2 proteins provided evidence of substitution events in the viral proteins involved in both host entry and genome replication. CONCLUSION: Globally circulating SARS-CoV-2 genomes could be classified into 5 major clades based on mutational profiles defined by an 11-nucleotide barcode. We have successfully developed a multiplexed sequencing-based, rapid genotyping protocol for high-throughput classification of major clade types of SARS-CoV-2 in clinical samples. This barcoding strategy will be required to monitor decreases in genetic diversity as treatment and vaccine approaches become widely available.

Influx of diverse, drug resistant and transmissible Plasmodium falciparum into a malaria-free setting in Qatar.

BACKGROUND: Successful control programs have impeded local malaria transmission in almost all Gulf Cooperation Council (GCC) countries: Qatar, Bahrain, Kuwait, Oman, the United Arab Emirates (UAE) and Saudi Arabia. Nevertheless, a prodigious influx of imported malaria via migrant workers sustains the threat of local transmission. Here we examine the origin of imported malaria in Qatar, assess genetic diversity and the prevalence of drug resistance genes in imported Plasmodium falciparum, and finally, address the potential for the reintroduction of local transmission. METHODS: This study examined imported malaria cases reported in Qatar, between 2013 and 2016. We focused on P. falciparum infections and estimated both total parasite and gametocyte density, using qPCR and qRT-PCR, respectively. We also examined ten neutral microsatellites and four genes associated with drug resistance, Pfmrp1, Pfcrt, Pfmdr1, and Pfkelch13, to assess the genetic diversity of imported P. falciparum strains, and the potential for propagating drug resistance genotypes respectively. RESULTS: The majority of imported malaria cases were P. vivax, while P. falciparum and mixed species infections (P. falciparum / P. vivax) were less frequent. The primary origin of P. vivax infection was the Indian subcontinent, while P. falciparum was mostly presented by African expatriates. Imported P. falciparum strains were highly diverse, carrying multiple genotypes, and infections also presented with early- and late-stage gametocytes. We observed a high prevalence of mutations implicated in drug resistance among these strains, including novel SNPs in Pfkelch13. CONCLUSIONS: The influx of genetically diverse P. falciparum, with multiple drug resistance markers and a high capacity for gametocyte production, represents a threat for the reestablishment of drug-resistant malaria into GCC countries. This scenario highlights the impact of mass international migration on the reintroduction of malaria to areas with absent or limited local transmission.

Integrated transcriptomic and proteomic analysis of pathogenic mycobacteria and their esx-1 mutants reveal secretion-dependent regulation of ESX-1 substrates and WhiB6 as a transcriptional regulator.

The mycobacterial type VII secretion system ESX-1 is responsible for the secretion of a number of proteins that play important roles during host infection. The regulation of the expression of secreted proteins is often essential to establish successful infection. Using transcriptome sequencing, we found that the abrogation of ESX-1 function in Mycobacterium marinum leads to a pronounced increase in gene expression levels of the espA operon during the infection of macrophages. In addition, the disruption of ESX-1-mediated protein secretion also leads to a specific down-regulation of the ESX-1 substrates, but not of the structural components of this system, during growth in culture medium. This effect is observed in both M. marinum and M. tuberculosis. We established that down-regulation of ESX-1 substrates is the result of a regulatory process that is influenced by the putative transcriptional regulator whib6, which is located adjacent to the esx-1 locus. In addition, the overexpression of the ESX-1-associated PE35/PPE68 protein pair resulted in a significantly increased secretion of the ESX-1 substrate EsxA, demonstrating a functional link between these proteins. Taken together, these data show that WhiB6 is required for the secretion-dependent regulation of ESX-1 substrates and that ESX-1 substrates are regulated independently from the structural components, both during infection and as a result of active secretion.

Enteric Infections Circulating during Hajj Seasons, 2011-2013.

Hajj, the annual Muslim pilgrimage to Mecca, Saudi Arabia, is a unique mass gathering event that raises public health concerns in the host country and globally. Although gastroenteritis and diarrhea are common among Hajj pilgrims, the microbial etiologies of these infections are unknown. We collected 544 fecal samples from pilgrims with medically attended diarrheal illness from 40 countries during the 2011-2013 Hajj seasons and screened the samples for 16 pathogens commonly associated with diarrheal infections. Bacteria were the main agents detected, in 82.9% of the 228 positive samples, followed by viral (6.1%) and parasitic (5.3%) agents. Salmonella spp., Shigella/enteroinvasive Escherichia coli, and enterotoxigenic E. coli were the main pathogens associated with severe symptoms. We identified genes associated with resistance to third-generation cephalosporins ≈40% of Salmonella- and E. coli-positive samples. Hajj-associated foodborne infections pose a major public health risk through the emergence and transmission of antimicrobial drug-resistant bacteria.

Candidate Essential Genes in Burkholderia cenocepacia J2315 Identified by Genome-Wide TraDIS.

Burkholderia cenocepacia infection often leads to fatal cepacia syndrome in cystic fibrosis patients. However, antibiotic therapy rarely results in complete eradication of the pathogen due to its intrinsic resistance to many clinically available antibiotics. Recent attention has turned to the identification of essential genes as the proteins encoded by these genes may serve as potential targets for development of novel antimicrobials. In this study, we utilized TraDIS (Transposon Directed Insertion-site Sequencing) as a genome-wide screening tool to facilitate the identification of B. cenocepacia genes essential for its growth and viability. A transposon mutant pool consisting of approximately 500,000 mutants was successfully constructed, with more than 400,000 unique transposon insertion sites identified by computational analysis of TraDIS datasets. The saturated library allowed for the identification of 383 genes that were predicted to be essential in B. cenocepacia. We extended the application of TraDIS to identify conditionally essential genes required for in vitro growth and revealed an additional repertoire of 439 genes to be crucial for B. cenocepacia growth under nutrient-depleted conditions. The library of B. cenocepacia mutants can subsequently be subjected to various biologically related conditions to facilitate the discovery of genes involved in niche adaptation as well as pathogenicity and virulence.

Genome-scale comparison of expanded gene families in Plasmodium ovale wallikeri and Plasmodium ovale curtisi with Plasmodium malariae and with other Plasmodium species.

Malaria in humans is caused by six species of Plasmodium parasites, of which the nuclear genome sequences for the two Plasmodium ovale spp., P. ovale curtisi and P. ovale wallikeri, and Plasmodium malariae have not yet been analyzed. Here we present an analysis of the nuclear genome sequences of these three parasites, and describe gene family expansions therein. Plasmodium ovale curtisi and P. ovale wallikeri are genetically distinct but morphologically indistinguishable and have sympatric ranges through the tropics of Africa, Asia and Oceania. Both P. ovale spp. show expansion of the surfin variant gene family, and an amplification of the Plasmodium interspersed repeat (pir) superfamily which results in an approximately 30% increase in genome size. For comparison, we have also analyzed the draft nuclear genome of P. malariae, a malaria parasite causing mild malaria symptoms with a quartan life cycle, long-term chronic infections, and wide geographic distribution. Plasmodium malariae shows only a moderate level of expansion of pir genes, and unique expansions of a highly diverged transmembrane protein family with over 550 members and the gamete P25/27 gene family. The observed diversity in the P. ovale wallikeri and P. ovale curtisi surface antigens, combined with their phylogenetic separation, supports consideration that the two parasites be given species status.

Genomic and Phenotypic Analyses Reveal the Emergence of an Atypical Salmonella enterica Serovar Senftenberg Variant in China.

Human infections with Salmonella enterica subspecies enterica serovar Senftenberg are often associated with exposure to poultry flocks, farm environments, or contaminated food. The recent emergence of multidrug-resistant isolates has raised public health concerns. In this study, comparative genomics and phenotypic analysis were used to characterize 14 Salmonella Senftenberg clinical isolates recovered from multiple outbreaks in Shenzhen and Shanghai, China, between 2002 and 2011. Single-nucleotide polymorphism analyses identified two phylogenetically distinct clades of S Senftenberg, designated SC1 and SC2, harboring variations in Salmonella pathogenicity island 1 (SPI-1) and SPI-2 and exhibiting distinct biochemical and phenotypic signatures. Although the two variants shared the same serotype, the SC2 isolates of sequence type 14 (ST14) harbored intact SPI-1 and -2 and hence were characterized by possessing efficient invasion capabilities. In contrast, the SC1 isolates had structural deletion patterns in both SPI-1 and -2 that correlated with an impaired capacity to invade cultured human cells and also the year of their isolation. These atypical SC1 isolates also lacked the capacity to produce hydrogen sulfide. These findings highlight the emergence of atypical Salmonella Senftenberg variants in China and provide genetic validation that variants lacking SPI-1 and regions of SPI-2, which leads to impaired invasion capacity, can still cause clinical disease. These data have identified an emerging public health concern and highlight the need to strengthen surveillance to detect the prevalence and transmission of nontyphoidal Salmonella species.

Genomic expression catalogue of a global collection of BCG vaccine strains show evidence for highly diverged metabolic and cell-wall adaptations.

Although Bacillus Calmette-Guérin (BCG) vaccines against tuberculosis have been available for more than 90 years, their effectiveness has been hindered by variable protective efficacy and a lack of lasting memory responses. One factor contributing to this variability may be the diversity of the BCG strains that are used around the world, in part from genomic changes accumulated during vaccine production and their resulting differences in gene expression. We have compared the genomes and transcriptomes of a global collection of fourteen of the most widely used BCG strains at single base-pair resolution. We have also used quantitative proteomics to identify key differences in expression of proteins across five representative BCG strains of the four tandem duplication (DU) groups. We provide a comprehensive map of single nucleotide polymorphisms (SNPs), copy number variation and insertions and deletions (indels) across fourteen BCG strains. Genome-wide SNP characterization allowed the construction of a new and robust phylogenic genealogy of BCG strains. Transcriptional and proteomic profiling revealed a metabolic remodeling in BCG strains that may be reflected by altered immunogenicity and possibly vaccine efficacy. Together, these integrated-omic data represent the most comprehensive catalogue of genetic variation across a global collection of BCG strains.

Genome Sequence of a Multidrug-Resistant Strain of Stenotrophomonas maltophilia with Carbapenem Resistance, Isolated from King Abdullah Medical City, Makkah, Saudi Arabia.

The emergence and spread of multidrug-resistant (MDR) bacteria have been regarded as major challenges among health care-associated infections worldwide. Here, we report the draft genome sequence of an MDR Stenotrophomonas maltophilia strain isolated in 2014 from King Abdulla Medical City, Makkah, Saudi Arabia.

Chromerid genomes reveal the evolutionary path from photosynthetic algae to obligate intracellular parasites.

The eukaryotic phylum Apicomplexa encompasses thousands of obligate intracellular parasites of humans and animals with immense socio-economic and health impacts. We sequenced nuclear genomes of Chromera velia and Vitrella brassicaformis, free-living non-parasitic photosynthetic algae closely related to apicomplexans. Proteins from key metabolic pathways and from the endomembrane trafficking systems associated with a free-living lifestyle have been progressively and non-randomly lost during adaptation to parasitism. The free-living ancestor contained a broad repertoire of genes many of which were repurposed for parasitic processes, such as extracellular proteins, components of a motility apparatus, and DNA- and RNA-binding protein families. Based on transcriptome analyses across 36 environmental conditions, Chromera orthologs of apicomplexan invasion-related motility genes were co-regulated with genes encoding the flagellar apparatus, supporting the functional contribution of flagella to the evolution of invasion machinery. This study provides insights into how obligate parasites with diverse life strategies arose from a once free-living phototrophic marine alga.

Comprehensive evaluation of Toxoplasma gondii VEG and Neospora caninum LIV genomes with tachyzoite stage transcriptome and proteome defines novel transcript features.

Toxoplasma gondii is an important protozoan parasite that infects all warm-blooded animals and causes opportunistic infections in immuno-compromised humans. Its closest relative, Neospora caninum, is an important veterinary pathogen that causes spontaneous abortion in livestock. Comparative genomics of these two closely related coccidians has been of particular interest to identify genes that contribute to varied host cell specificity and disease. Here, we describe a manual evaluation of these genomes based on strand-specific RNA sequencing and shotgun proteomics from the invasive tachyzoite stages of these two parasites. We have corrected predicted structures of over one third of the previously annotated gene models and have annotated untranslated regions (UTRs) in over half of the predicted protein-coding genes. We observe distinctly long UTRs in both the organisms, almost four times longer than other model eukaryotes. We have also identified a putative set of cis-natural antisense transcripts (cis-NATs) and long intergenic non-coding RNAs (lincRNAs). We have significantly improved the annotation quality in these genomes that would serve as a manually curated dataset for Toxoplasma and Neospora research communities.

Identification of a novel conjugative plasmid in mycobacteria that requires both type IV and type VII secretion.

UNLABELLED: Conjugative plasmids have been identified in a wide variety of different bacteria, ranging from proteobacteria to firmicutes, and conjugation is one of the most efficient routes for horizontal gene transfer. The most widespread mechanism of plasmid conjugation relies on different variants of the type IV secretion pathway. Here, we describe the identification of a novel type of conjugative plasmid that seems to be unique for mycobacteria. Interestingly, while this plasmid is efficiently exchanged between different species of slow-growing mycobacteria, including Mycobacterium tuberculosis, it could not be transferred to any of the fast-growing mycobacteria tested. Genetic analysis of the conjugative plasmid showed the presence of a locus containing homologues of three type IV secretion system components and a relaxase. In addition, a new type VII secretion locus was present. Using transposon insertion mutagenesis, we show that in fact both these secretion systems are essential for conjugation, indicating that this plasmid represents a new class of conjugative plasmids requiring two secretion machineries. This plasmid could form a useful new tool to exchange or introduce DNA in slow-growing mycobacteria. IMPORTANCE: Conjugative plasmids play an important role in horizontal gene transfer between different bacteria and, as such, in their adaptation and evolution. This effect is most obvious in the spread of antibiotic resistance genes. Thus far, conjugation of natural plasmids has been described only rarely for mycobacterial species. In fact, it is generally accepted that M. tuberculosis does not show any recent sign of horizontal gene transfer. In this study, we describe the identification of a new widespread conjugative plasmid that can also be efficiently transferred to M. tuberculosis. This plasmid therefore poses both a threat and an opportunity. The threat is that, through the acquisition of antibiotic resistance markers, this plasmid could start a rapid spread of antibiotic resistance genes between pathogenic mycobacteria. The opportunity is that we could use this plasmid to generate new tools for the efficient introduction of foreign DNA in slow-growing mycobacteria.

The population structure of Vibrio cholerae from the Chandigarh Region of Northern India.

BACKGROUND: Cholera infection continues to be a threat to global public health. The current cholera pandemic associated with Vibrio cholerae El Tor has now been ongoing for over half a century. METHODOLOGY/PRINCIPAL FINDINGS: Thirty-eight V. cholerae El Tor isolates associated with a cholera outbreak in 2009 from the Chandigarh region of India were characterised by a combination of microbiology, molecular typing and whole-genome sequencing. The genomic analysis indicated that two clones of V. cholera circulated in the region and caused disease during this time. These clones fell into two distinct sub-clades that map independently onto wave 3 of the phylogenetic tree of seventh pandemic V. cholerae El Tor. Sequence analyses of the cholera toxin gene, the Vibrio seventh Pandemic Island II (VSPII) and SXT element correlated with this phylogenetic position of the two clades on the El Tor tree. The clade 2 isolates, characterized by a drug-resistant profile and the expression of a distinct cholera toxin, are closely related to the recent V. cholerae isolated elsewhere, including Haiti, but fell on a distinct branch of the tree, showing they were independent outbreaks. Multi-Locus Sequence Typing (MLST) distinguishes two sequence types among the 38 isolates, that did not correspond to the clades defined by whole-genome sequencing. Multi-Locus Variable-length tandem-nucleotide repeat Analysis (MLVA) identified 16 distinct clusters. CONCLUSIONS/SIGNIFICANCE: The use of whole-genome sequencing enabled the identification of two clones of V. cholerae that circulated during the 2009 Chandigarh outbreak. These clones harboured a similar structure of ICEVchHai1 but differed mainly in the structure of CTX phage and VSPII. The limited capacity of MLST and MLVA to discriminate between the clones that circulated in the 2009 Chandigarh outbreak highlights the value of whole-genome sequencing as a route to the identification of further genetic markers to subtype V. cholerae isolates.

SVAMP: sequence variation analysis, maps and phylogeny.

SUMMARY: SVAMP is a stand-alone desktop application to visualize genomic variants (in variant call format) in the context of geographical metadata. Users of SVAMP are able to generate phylogenetic trees and perform principal coordinate analysis in real time from variant call format (VCF) and associated metadata files. Allele frequency map, geographical map of isolates, Tajima's D metric, single nucleotide polymorphism density, GC and variation density are also available for visualization in real time. We demonstrate the utility of SVAMP in tracking a methicillin-resistant Staphylococcus aureus outbreak from published next-generation sequencing data across 15 countries. We also demonstrate the scalability and accuracy of our software on 245 Plasmodium falciparum malaria isolates from three continents. AVAILABILITY AND IMPLEMENTATION: The Qt/C++ software code, binaries, user manual and example datasets are available at http://cbrc.kaust.edu.sa/svamp CONTACT: arnab.pain@kaust.edu.sa or arnab.pain@cantab.net SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

READSCAN: a fast and scalable pathogen discovery program with accurate genome relative abundance estimation.

UNLABELLED: READSCAN is a highly scalable parallel program to identify non-host sequences (of potential pathogen origin) and estimate their genome relative abundance in high-throughput sequence datasets. READSCAN accurately classified human and viral sequences on a 20.1 million reads simulated dataset in <27 min using a small Beowulf compute cluster with 16 nodes (Supplementary Material). AVAILABILITY: http://cbrc.kaust.edu.sa/readscan.