Genomic Epidemiology of Large Blastomycosis Outbreak, Ontario, Canada, 2021.

Using phylogenomic analysis, we provide genomic epidemiology analysis of a large blastomycosis outbreak in Ontario, Canada, caused by Blastomyces gilchristii. The outbreak occurred in a locale where blastomycosis is rarely diagnosed, signaling a possible shift in geographically associated incidence patterns. Results elucidated fungal population genetic structure, enhancing understanding of the outbreak.

Detection of the novel SARS-CoV-2 European lineage B.1.177 in Ontario, Canada.

Background: Travel-related dissemination of SARS-CoV-2 continues to contribute to the global pandemic. A novel SARS-CoV-2 lineage (B.1.177) reportedly arose in Spain in the summer of 2020, with subsequent spread across Europe linked to travel by infected individuals. Surveillance and monitoring through the use of whole genome sequencing (WGS) offers insights into the global and local movement of pathogens such as SARS-CoV-2 and can detect introductions of novel variants. Methods: We analysed the genomes of SARS-CoV-2 sequenced for surveillance purposes from specimens received by Public Health Ontario (Sept 6 - Oct 10, 2020), collected from individuals in eastern Ontario, which comprised the study sample. Taxonomic lineages were identified using pangolin (v2.08) and phylogenetic analysis incorporated publicly available genomes covering the same time period as the study sample. Epidemiological data collected from laboratory requisitions and standard reportable disease case investigation was integrated into the analysis. Results: Genomic surveillance identified a COVID-19 case with SARS-CoV-2 lineage B.1.177 from an individual in eastern Ontario in late September, 2020. The individual had recently returned from Europe. Genomic analysis with publicly available data indicate the most closely related genomes to this specimen were from Southern Europe. Genomic surveillance did not identify further cases with this lineage. Conclusions: Genomic surveillance allowed for early detection of a novel SARS-CoV-2 lineage in Ontario which was deemed to be travel related. This type of genomic-based surveillance is a key tool to measure the effectiveness of public health measures such as mandatory self-isolation for returned travellers, aimed at preventing onward transmission of newly introduced lineages of SARS-CoV-2.

Are shed hair genomes the most effective noninvasive resource for estimating relationships in the wild?

Knowledge of relationships in wild populations is critical for better understanding mating systems and inbreeding scenarios to inform conservation strategies for endangered species. To delineate pedigrees in wild populations, study genetic connectivity, study genotype-phenotype associations, trace individuals, or track wildlife trade, many identified individuals need to be genotyped at thousands of loci, mostly from noninvasive samples. This requires us to (a) identify the most common noninvasive sample available from identified individuals, (b) assess the ability to acquire genome-wide data from such samples, and (c) evaluate the quality of such genome-wide data, and its ability to reconstruct relationships between animals within a population.We followed identified individuals from a wild endangered tiger population and found that shed hair samples were the most common compared to scat samples, opportunistically found carcasses, and opportunistic invasive samples. We extracted DNA from these samples, prepared whole genome sequencing libraries, and sequenced genomes from these.Whole genome sequencing methods resulted in between 25%-98% of the genome sequenced for five such samples. Exploratory population genetic analyses revealed that these data were free of holistic biases and could recover expected population structure and relatedness. Mitochondrial genomes recovered matrilineages in accordance with long-term monitoring data. Even with just five samples, we were able to uncover the matrilineage for three individuals with unknown ancestry.In summary, we demonstrated that noninvasive shed hair samples yield adequate quality and quantity of DNA in conjunction with sensitive library preparation methods, and provide reliable data from hundreds of thousands of SNPs across the genome. This makes shed hair an ideal noninvasive resource for studying individual-based genetics of elusive endangered species in the wild.

Non-coding and Coding Transcriptional Profiles Are Significantly Altered in Pediatric Retinoblastoma Tumors.

Retinoblastoma is a rare pediatric tumor of the retina, caused by the homozygous loss of the Retinoblastoma 1 (RB1) tumor suppressor gene. Previous microarray studies have identified changes in the expression profiles of coding genes; however, our understanding of how non-coding genes change in this tumor is absent. This is an important area of research, as in many adult malignancies, non-coding genes including LNC-RNAs are used as biomarkers to predict outcome and/or relapse. To establish a complete and in-depth RNA profile, of both coding and non-coding genes, in Retinoblastoma tumors, we conducted RNA-seq from a cohort of tumors and normal retina controls. This analysis identified widespread transcriptional changes in the levels of both coding and non-coding genes. Unexpectedly, we also found rare RNA fusion products resulting from genomic alterations, specific to Retinoblastoma tumor samples. We then determined whether these gene expression changes, of both coding and non-coding genes, were also found in a completely independent Retinoblastoma cohort. Using our dataset, we then profiled the potential effects of deregulated LNC-RNAs on the expression of neighboring genes, the entire genome, and on mRNAs that contain a putative area of homology. This analysis showed that most deregulated LNC-RNAs do not act locally to change the transcriptional environment, but potentially function to modulate genes at distant sites. From this analysis, we selected a strongly down-regulated LNC-RNA in Retinoblastoma, DRAIC, and found that restoring DRAIC RNA levels significantly slowed the growth of the Y79 Retinoblastoma cell line. Collectively, our work has generated the first non-coding RNA profile of Retinoblastoma tumors and has found that these tumors show widespread transcriptional deregulation.

Metabolite systems profiling identifies exploitable weaknesses in retinoblastoma.

Retinoblastoma (RB) is a childhood eye cancer. Currently, chemotherapy, local therapy, and enucleation are the main ways in which these tumors are managed. The present work is the first study that uses constraint-based reconstruction and analysis approaches to identify and explain RB-specific survival strategies, which are RB tumor specific. Importantly, our model-specific secretion profile is also found in RB1-depleted human retinal cells in vitro and suggests that novel biomarkers involved in lipid metabolism may be important. Finally, RB-specific synthetic lethals have been predicted as lipid and nucleoside transport proteins that can aid in novel drug target development.

Somatic Variations in Cervical Cancers in Indian Patients.

There are very few reports that describe the mutational landscape of cervical cancer, one of the leading cancers in Indian women. The aim of the present study was to investigate the somatic mutations that occur in cervical cancer. Whole exome sequencing of 10 treatment naïve tumour biopsies with matched blood samples, from a cohort of Indian patients with locally advanced disease, was performed. The data revealed missense mutations across 1282 genes, out of 1831 genes harbouring somatic mutations. These missense mutations (nonsynonymous + stop-gained) when compared with pre-existing mutations in the COSMIC database showed that 272 mutations in 250 genes were already reported although from cancers other than cervical cancer. More than 1000 novel somatic variations were obtained in matched tumour samples. Pathways / genes that are frequently mutated in various other cancers were found to be mutated in cervical cancers. A significant enrichment of somatic mutations in the MAPK pathway was observed, some of which could be potentially targetable. This is the first report of whole exome sequencing of well annotated cervical cancer samples from Indian women and helps identify trends in mutation profiles that are found in an Indian cohort of cervical cancer.

A genomewide mutagenesis screen identifies multiple genes contributing to Vi capsular expression in Salmonella enterica serovar Typhi.

A transposon-based, genomewide mutagenesis screen exploiting the killing activity of a lytic ViII bacteriophage was used to identify Salmonella enterica serovar Typhi genes that contribute to Vi polysaccharide capsule expression. Genes enriched in the screen included those within the viaB locus (tviABCDE and vexABCDE) as well as oxyR, barA/sirA, and yrfF, which have not previously been associated with Vi expression. The role of these genes in Vi expression was confirmed by constructing defined null mutant derivatives of S. Typhi, and these were negative for Vi expression as determined by agglutination assays with Vi-specific sera or susceptibility to Vi-targeting bacteriophages. Transcriptome analysis confirmed a reduction in expression from the viaB locus in these S. Typhi mutant derivatives and defined regulatory networks associated with Vi expression.

Characterization of the yehUT two-component regulatory system of Salmonella enterica Serovar Typhi and Typhimurium.

Proteins exhibiting hyper-variable sequences within a bacterial pathogen may be associated with host adaptation. Several lineages of the monophyletic pathogen Salmonella enterica serovar Typhi (S. Typhi) have accumulated non-synonymous mutations in the putative two-component regulatory system yehUT. Consequently we evaluated the function of yehUT in S. Typhi BRD948 and S. Typhimurium ST4/74. Transcriptome analysis identified the cstA gene, encoding a carbon starvation protein as the predominantly yehUT regulated gene in both these serovars. Deletion of yehUT had no detectable effect on the ability of these mutant Salmonella to invade cultured epithelial cells (S. Typhi and S. Typhimurium) or induce colitis in a murine model (S. Typhimurium only). Growth, metabolic and antimicrobial susceptibility tests identified no obvious influences of yehUT on these phenotypes.

An overview of prokaryotic transcription factors : a summary of function and occurrence in bacterial genomes.

Transcriptional initiation is arguably the most important control point for gene expression. It is regulated by a combination of factors, including DNA sequence and its three-dimensional topology, proteins and small molecules. In this chapter, we focus on the trans-acting factors of bacterial regulation. Initiation begins with the recruitment of the RNA polymerase holoenzyme to a specific locus upstream of the gene known as its promoter. The sigma factor, which is a component of the holoenzyme, provides the most fundamental mechanisms for orchestrating broad changes in gene expression state. It is responsible for promoter recognition as well as recruiting the holoenzyme to the promoter. Distinct sigma factors compete with for binding to a common pool of RNA polymerases, thus achieving condition-dependent differential expression. Another important class of bacterial regulators is transcription factors, which activate or repress transcription of target genes typically in response to an environmental or cellular trigger. These factors may be global or local depending on the number of genes and range of cellular functions that they target. The activities of both global and local transcription factors may be regulated either at a post-transcriptional level via signal-sensing protein domains or at the level of their own expression. In addition to modulating polymerase recruitment to promoters, several global factors are considered as "nucleoid-associated proteins" that impose structural constraints on the chromosome by altering the conformation of the bound DNA, thus influencing other processes involving DNA such as replication and recombination. This chapter concludes with a discussion of how regulatory interactions between transcription factors and their target genes can be represented as a network.

Staphylococcus aureus nasal carriage and its contributing factors.

Staphylococcus aureus is a medically important pathogen that is often acquired from hospital settings (nosocomial) as well as from the community (community acquired). Bacteria that reside in anterior nares of hosts serve as reservoirs for both the spread of the pathogen and predispose the host to subsequent infections. Here, we will review the extent and variability of nasal carriage, and the possible causative factors--both from the host and the bacterium. We also discuss the existing molecular typing techniques used for studying variations among strains of S. aureus. Finally, we discuss the possible areas of studies that are open in this field. Given the pathogen's importance in healthcare setting, such areas of study vary vastly, from fundamental research to applied medical care and use of alternative medical regimes for control of S. aureus nasal carriage. Unsurprisingly, our conclusions also underscore the importance of making policy decisions based on local ethnic and socioeconomic population structure.

Pathogenesis gene families in the common minimal genome of Staphylococcus aureus are hypervariable.

Staphylococcus aureus is a versatile pathogen that shows high levels of inter-strain genetic variability and positive evolution in certain pathogenesis-related genes. Apart from gene content differences, variability in shared genes may affect pathogenicity. Studying such variability requires that the common minimal genome (CMG) be identified. In this study, we have surveyed the CMG of S. aureus with respect to variability amongst orthologous family members, and determined that genes involved in pathogenesis preferentially accumulate variations. A negative correlation between variability of genes and their evolution was found, suggesting a preservation of host-specific function while exhibiting sequence diversity. Variation in key pathogenesis genes in S. aureus might predispose them to functional modulation, thereby playing an important role in evasion of host immunity.

Genome sequencing and analysis reveals possible determinants of Staphylococcus aureus nasal carriage.

BACKGROUND: Nasal carriage of Staphylococcus aureus is a major risk factor in clinical and community settings due to the range of etiologies caused by the organism. We have identified unique immunological and ultrastructural properties associated with nasal carriage isolates denoting a role for bacterial factors in nasal carriage. However, despite extensive molecular level characterizations by several groups suggesting factors necessary for colonization on nasal epithelium, genetic determinants of nasal carriage are unknown. Herein, we have set a genomic foundation for unraveling the bacterial determinants of nasal carriage in S. aureus. RESULTS: MLST analysis revealed no lineage specific differences between carrier and non-carrier strains suggesting a role for mobile genetic elements. We completely sequenced a model carrier isolate (D30) and a model non-carrier strain (930918-3) to identify differential gene content. Comparison revealed the presence of 84 genes unique to the carrier strain and strongly suggests a role for Type VII secretion systems in nasal carriage. These genes, along with a putative pathogenicity island (SaPIBov) present uniquely in the carrier strains are likely important in affecting carriage. Further, PCR-based genotyping of other clinical isolates for a specific subset of these 84 genes raise the possibility of nasal carriage being caused by multiple gene sets. CONCLUSION: Our data suggest that carriage is likely a heterogeneic phenotypic trait and implies a role for nucleotide level polymorphism in carriage. Complete genome level analyses of multiple carriage strains of S. aureus will be important in clarifying molecular determinants of S. aureus nasal carriage.

Codon choice in genes depends on flanking sequence information--implications for theoretical reverse translation.

Algorithms for theoretical reverse translation have direct applications in degenerate PCR. The conventional practice is to create several degenerate primers each of which variably encode the peptide region of interest. In the current work, for each codon we have analyzed the flanking residues in proteins and determined their influence on codon choice. From this, we created a method for theoretical reverse translation that includes information from flanking residues of the protein in question. Our method, named the neighbor correlation method (NCM) and its enhancement, the consensus-NCM (c-NCM) performed significantly better than the conventional codon-usage statistic method (CSM). Using the methods NCM and c-NCM, we were able to increase the average sequence identity from 77% up to 81%. Furthermore, we revealed a significant increase in coverage, at 80% identity, from < 20% (CSM) to > 75% (c-NCM). The algorithms, their applications and implications are discussed herein.

Role of hydrophobic interactions and salt-bridges in beta-hairpin folding.

Beta-hairpins are the simplest form of beta-sheets which, due to the presence of long-range interactions, can be considered as tertiary structures. Molecular dynamics simulation is a powerful tool that can unravel whole pathways of protein folding/unfolding at atomic resolution. We have performed several molecular dynamics simulations, to a total of over 250 ns, of a beta-hairpin peptide in water using GROMACS. We show that hydrophobic interactions are necessary for initiating the folding of the peptide. Once formed, the peptide is stabilized by hydrogen bonds and disruption of hydrophobic interactions in the folded peptide does not denature the structure. In the absence of hydrophobic interactions, the peptide fails to fold. However, the introduction of a salt-bridge compensates for the loss of hydrophobic interactions to a certain extent.

Promoter addresses: revelations from oligonucleotide profiling applied to the Escherichia coli genome.

BACKGROUND: Transcription is the first step in cellular information processing. It is regulated by cis-acting elements such as promoters and operators in the DNA, and trans-acting elements such as transcription factors and sigma factors. Identification of cis-acting regulatory elements on a genomic scale requires computational analysis. RESULTS: We have used oligonucleotide profiling to predict regulatory regions in a bacterial genome. The method has been applied to the Escherichia coli K12 genome and the results analyzed. The information content of the putative regulatory oligonucleotides so predicted is validated through intra-genomic analyses, correlations with experimental data and inter-genome comparisons. Based on the results we have proposed a model for the bacterial promoter. The results show that the method is capable of identifying, in the E.coli genome, cis-acting elements such as TATAAT (sigma70 binding site), CCCTAT (1 base relative of sigma32 binding site), CTATNN (LexA binding site), AGGA-containing hexanucleotides (Shine Dalgarno consensus) and CTAG-containing hexanucleotides (core binding sites for Trp and Met repressors). CONCLUSION: The method adopted is simple yet effective in predicting upstream regulatory elements in bacteria. It does not need any prior experimental data except the sequence itself. This method should be applicable to most known genomes. Profiling, as applied to the E.coli genome, picks up known cis-acting and regulatory elements. Based on the profile results, we propose a model for the bacterial promoter that is extensible even to eukaryotes. The model is that the core promoter lies within a plateau of bent AT-rich DNA. This bent DNA acts as a homing segment for the sigma factor to recognize the promoter. The model thus suggests an important role for local landscapes in prokaryotic and eukaryotic gene regulation.