Single Nucleotide Polymorphism Analysis Research Articles

Genomic Epidemiology of Carbapenem-Resistant Klebsiella in Qatar: Emergence and Dissemination of Hypervirulent Klebsiella pneumoniae Sequence Type 383 Strains.

The emergence of carbapenem-resistant, hypervirulent Klebsiella pneumoniae is a new threat to health care. We studied the molecular epidemiology of carbapenem-resistant Klebsiella pneumoniae isolates in Qatar using whole-genome sequence data. We also characterized the prevalence and genetic basis of hypervirulent phenotypes and established the virulence potential using a Galleria mellonella model. Of 100 Klebsiella isolates studied, NDM and OXA-48 were the most common carbapenemases. Core genome single-nucleotide polymorphism (SNP) analysis indicated the presence of diverse sequence types and clonal lineages; isolates belonging to Klebsiella quasipneumoniae subsp. quasipneumoniae sequence type 196 (ST196) and ST1416 may be disseminated among several health care centers. Ten K. pneumoniae isolates carried rmpA and/or truncated rmpA2, and 2 isolates belonged to KL2, indicating low prevalence of classical hypervirulent isolates. Isolates carrying both carbapenem resistance and hypervirulence genes were confined mainly to ST231 and ST383 isolates. One ST383 isolate was further investigated by MinION sequencing, and the assembled genome indicated that blaNDM was located on an IncHI1B-type plasmid (pFQ61_ST383_NDM-5) which coharbored several virulence factors, including the regulator of the mucoid phenotype (rmpA), the regulator of mucoid phenotype 2 (rmpA2), and aerobactin (iucABCD and iutA), likely resulting from recombination events. Comparative genomics indicated that this hybrid plasmid may be present in two additional Qatari ST383 isolates. Carbapenem-resistant, hypervirulent K. pneumoniae ST383 isolates pose an emerging threat to global health due to their simultaneous hypervirulence and multidrug resistance.

Genomic Surveillance of SARS-CoV-2 Variants That Emerged in South and Southeast Asia during Early 2022.

The continuously emerging new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have made the global coronavirus disease 2019 (COVID-19) pandemic unpredictable. Since the beginning of the pandemic, densely populated South and Southeast Asia have suffered great losses due to multiple COVID-19 surges because of vaccine and other medical resource shortages. Therefore, it is crucial to closely monitor the SARS-CoV-2 epidemic and to understand the evolutionary and transmission characteristics of SARS-CoV-2 in these regions. Here, we document the evolution of epidemic strains in the Philippines, Pakistan, and Malaysia from late 2021 to early 2022. Our results confirmed the circulation of at least five SARS-CoV-2 genotypes in these countries in January 2022, when Omicron BA.2, with a detection rate of 69.11%, replaced Delta B.1.617 as the dominant strain. Single-nucleotide polymorphism analysis indicated the distinct evolutionary directions of the Omicron and Delta isolates, with S, Nsp1, and Nsp6 genes potentially playing a significant role in the host adaptation of the Omicron strain. These findings are able to provide insights for predicting the evolutionary direction of SARS-CoV-2 in terms of variant competition, developing multi-part vaccines, and to support the evaluation and adjustment of current surveillance, prevention, and control strategies in South and Southeast Asia.

Association analysis of PMEL gene expression and single nucleotide polymorphism with plumage color in quail

To explore the relationship between PMEL gene and quail plumage color, to provide a reference for subsequent quail plumage color breeding. In this experiment, RT-qPCR technology was used to analyze the relative mRNA expression levels of Korean quail (maroon) and Beijing white quail embryos at different developmental stages. Two SNPs in PMEL gene were screened based on the RNA-Seq data of skin tissues of Korean quail and Beijing white quail during embryonic stage. The KASP technology was used for genotyping in the resource population and correlation analysis was carried out with the plumage color traits of quail. Finally, the bioinformatics technology was used to predict the effects of these two SNPs on the structure and function of the encoded protein. The results showed that the expression levels of PMEL gene during the embryonic development of Beijing white quail were extremely significantly higher than that of Korean quail (p < 0.01). The frequency distribution of the three genotypes (AA, AB, and BB) of the Beijing white quail at the c. 1030C > T and c. 1374A > G mutation sites were extremely significantly different from that of the Korean quail (p < 0.01). And there was a significant correlation between the c. 1374A > G mutation site with white plumage phenotype. Bioinformatics analysis showed that SNP1 (c. c1030t) located in exon 6 was a harmful mutation site, and SNP2 (c. a1374g) located in exon 7 was a neutral mutation site. Protein conservation prediction showed that the coding protein P344S site caused by SNP1 (c. c1030t) site and the coding protein I458M site caused by SNP2 (c. g2129a) site were non-conservative sites. The results of this experiment showed that the PMEL gene was associated with the plumage color traits of quail and could be used as a candidate gene for studying the plumage color of quail.

In host evolution of Exophiala dermatitidis in cystic fibrosis lung micro-environment.

Individuals with cystic fibrosis (CF) are susceptible to chronic lung infections that lead to inflammation and irreversible lung damage. While most respiratory infections that occur in CF are caused by bacteria, some are dominated by fungi such as the slow-growing black yeast Exophiala dermatitidis. Here, we analyze isolates of E. dermatitidis cultured from two samples, collected from a single subject 2 years apart. One isolate genome was sequenced using long-read Nanopore technology as an in-population reference to use in comparative single nucleotide polymorphism and insertion-deletion variant analyses of 23 isolates. We then used population genomics and phylo-genomics to compare the isolates to each other as well as the reference genome strain E. dermatitidis NIH/UT8656. Within the CF lung population, three E. dermatitidis clades were detected, each with varying mutation rates. Overall, the isolates were highly similar suggesting that they were recently diverged. All isolates were MAT 1-1, which was consistent with their high relatedness and the absence of evidence for mating or recombination between isolates. Phylogenetic analysis grouped sets of isolates into clades that contained isolates from both early and late time points indicating there are multiple persistent lineages. Functional assessment of variants unique to each clade identified alleles in genes that encode transporters, cytochrome P450 oxidoreductases, iron acquisition, and DNA repair processes. Consistent with the genomic heterogeneity, isolates showed some stable phenotype heterogeneity in melanin production, subtle differences in antifungal minimum inhibitory concentrations, and growth on different substrates. The persistent population heterogeneity identified in lung-derived isolates is an important factor to consider in the study of chronic fungal infections, and the analysis of changes in fungal pathogens over time may provide important insights into the physiology of black yeasts and other slow-growing fungi in vivo.

Whole-genome sequencing analysis of molecular epidemiology and silent transmissions causing meticillin-resistant Staphylococcus aureus bloodstream infections in a university hospital

The emergence of novel genomic-type clones, such as community-associated meticillin-resistant Staphylococcus aureus (MRSA) and livestock-associated MRSA, and their invasion into hospitals have become major concerns worldwide; however, little information is available regarding the prevalence of MRSA in Japan. Whole-genome sequencing (WGS) has been conducted to analyse various pathogens worldwide. Therefore, it is important to establish a genome database of clinical MRSA isolates available in Japan. A molecular epidemiological analysis of MRSA strains isolated from bloodstream-infected patients in a Japanese university hospital was conducted using WGS and single-nucleotide polymorphism (SNP) analysis. Additionally, through a review of patients' clinical characteristics, the effectiveness of SNP analysis as a tool for detecting silent nosocomial transmission that may be missed by other methods was evaluated in diverse settings and various time points of detection. Polymerase-chain-reaction-based staphylococcal cassette chromosome mec (SCCmec) typing was performed using 135 isolates obtained between 2014 and 2018, and WGS was performed using 88 isolates obtained between 2015 and 2017. SCCmec type II strains, prevalent in 2014, became rare in 2018, whereas the prevalence of SCCmec type IV strains increased from 18.75% to 83.87% of the population, and became the dominant clones. Clonal complex (CC) 5 CC8 and CC1 were detected between 2015 and 2017, with CC1 being dominant. In 88 cases, SNP analyses revealed nosocomial transmissions among 20 patients which involved highly homologous strains. Routine monitoring of MRSA by whole-genome analysis is effective not only for gaining knowledge regarding molecular epidemiology, but also for detecting silent nosocomial transmission.

Fine-scale genomic tracking of Ross River virus using nanopore sequencing

BackgroundRoss River virus (RRV) is Australia’s most common and widespread mosquito-transmitted arbovirus and is of significant public health concern. With increasing anthropogenic impacts on wildlife and mosquito populations, it is important that we understand how RRV circulates in its endemic hotspots to determine where public health efforts should be directed. Current surveillance methods are effective in locating the virus but do not provide data on the circulation of the virus and its strains within the environment. This study examined the ability to identify single nucleotide polymorphisms (SNPs) within the variable E2/E3 region by generating full-length haplotypes from a range of mosquito trap-derived samples.MethodsA novel tiled primer amplification workflow for amplifying RRV was developed with analysis using Oxford Nanopore Technology’s MinION and a custom ARTIC/InterARTIC bioinformatic protocol. By creating a range of amplicons across the whole genome, fine-scale SNP analysis was enabled by specifically targeting the variable region that was amplified as a single fragment and established haplotypes that informed spatial-temporal variation of RRV in the study site in Victoria.ResultsA bioinformatic and laboratory pipeline was successfully designed and implemented on mosquito whole trap homogenates. Resulting data showed that genotyping could be conducted in real time and that whole trap consensus of the viruses (with major SNPs) could be determined in a timely manner. Minor variants were successfully detected from the variable E2/E3 region of RRV, which allowed haplotype determination within complex mosquito homogenate samples.ConclusionsThe novel bioinformatic and wet laboratory methods developed here will enable fast detection and characterisation of RRV isolates. The concepts presented in this body of work are transferable to other viruses that exist as quasispecies in samples. The ability to detect minor SNPs, and thus haplotype strains, is critically important for understanding the epidemiology of viruses their natural environment.Graphical

Increased genetic diversity and immigration after West Nile virus emergence in American crows: No evidence for a genetic bottleneck.

Infectious diseases can cause steep declines in wildlife populations, leading to changes in genetic diversity that may affect the susceptibility of individuals to infection and the overall resilience of populations to pathogen outbreaks. Here, we examine evidence for a genetic bottleneck in a population of American crows (Corvus brachyrhynchos) before and after the emergence of West Nile virus (WNV). More than 50% of marked birds in this population were lost over the 2-year period of the epizootic, representing a 10-fold increase in adult mortality. Using analyses of single-nucleotide polymorphisms (SNPs) and microsatellite markers, we tested for evidence of a genetic bottleneck and compared levels of inbreeding and immigration in the pre- and post-WNV populations. Counter to expectations, genetic diversity (allelic diversity and the number of new alleles) increased after WNV emergence. This was likely due to increases in immigration, as the estimated membership coefficients were lower in the post-WNV population. Simultaneously, however, the frequency of inbreeding appeared to increase: Mean inbreeding coefficients were higher among SNP markers, and heterozygosity-heterozygosity correlations were stronger among microsatellite markers, in the post-WNV population. These results indicate that loss of genetic diversity at the population level is not an inevitable consequence of a population decline, particularly in the presence of gene flow. The changes observed in post-WNV crows could have very different implications for their response to future pathogen risks, potentially making the population as a whole more resilient to a changing pathogen community, while increasing the frequency of inbred individuals with elevated susceptibility to disease.

Structural and genomic analysis of single nucleotide polymorphisms in human host factor endothelial protein C receptor (EPCR) reveals complex interplay with malaria parasites.

Plasmodium parasites responsible for malaria follow a complex life cycle of which half takes place inside the human host. Parasites present diverse antigens at different stages of their life cycle and interact with many surface molecules to attach to and enter host cells. The CIDRα1 domain of Plasmodium falciparum Erythrocyte Membrane Protein 1 (PfEMP1) in infected erythrocytes adheres to one such vascular receptor endothelial protein C receptor (EPCR). EPCR is implicated in the pathogenesis of severe malaria as preferential binding of CIDRα1 to endothelium results in widespread sequestration of infected erythrocytes leading to endothelium inflammation and severe disease. A single EPCR variant S219G is clinically reported to provide protection from severe malaria. In this work, we have collated all single nucleotide polymorphisms (SNPs) in EPCR from dbSNP. We structurally mapped the SNPs on the three-dimensional complex of EPCR and PfEMP1 CIDRα1. Analysis shows that most EPCR mutations lie on the receptor surface and are non-conservative. Of the 11 mutations in the CIDRα1-interaction region of EPCR, S88P, L96V/I, and R98L/H/P/C are seen with comparably higher occurrences in diverse populations. Our structural analysis details a framework of the interactions between the parasite ligand and host factor EPCR. These structural glimpses provide a blueprint for designing both field-based variant sequencing studies and vaccine development.

Challenges Associated with Investigating Salmonella Enteritidis with Low Genomic Diversity in New York State: The Impact of Adjusting Analytical Methods and Correlation with Epidemiological Data.

Defining investigation-worthy genomic clusters among strains of Salmonella Enteritidis is challenging because of their highly clonal nature. We investigated a cluster identified by core genome multilocus sequence typing (cgMLST) consisting of 265 isolates with isolation dates spanning two and a half years. This cluster experienced chaining, growing to a range of 14 alleles. The volume of isolates and broad allele range of this cluster made it difficult to ascertain whether it represented a common-source outbreak. We explored laboratory-based methods to subdivide and refine this cluster. These methods included using cgMLST with a narrower allele range, whole genome multilocus sequence typing (wgMLST) and high-quality single-nucleotide polymorphism (hqSNP) analysis. At each analysis level, epidemiologists retroactively reviewed exposures, geography, and temporality for potential commonalities. Lowering the threshold to 0 alleles using cgMLST proved an effective method to refine this analysis, resulting in this large cluster being subdivided into 34 smaller clusters. Additional analysis by wgMLST and hqSNP provided enhanced cluster resolution, with the majority of clusters being further refined. These analysis methods combined with more stringent allele thresholds and layering of epidemiologic data proved useful in helping to subdivide this large cluster into actionable subclusters.

Immunity gene IFITM3 variant: relation to cognition and Alzheimer pathology

AbstractBackgroundIFITM3, an innate immune response protein and inhibitor of viral infection, was reported to modulate amyloid‐β production in Alzheimer’s disease (AD). We aimed to identify single‐nucleotide polymorphisms (SNPs) in IFITM3 associated with cognition and AD biomarkers.MethodWe used genetic, longitudinal cognition and AD biomarker data from Alzheimer’s Disease Neuroimaging Initiative (ADNI; N = 1,565) and AddNeuroMed (N = 633) as discovery and replication samples, respectively. First, we performed gene‐based association analysis of SNPs in IFITM3 with cognitive performance. Second, we performed SNP‐based association analysis in IFITM3 with cognitive decline and AD biomarkers from amyloid positron emission tomography (PET), cerebrospinal fluid (CSF), and magnetic resonance imaging (MRI).ResultGene‐based association analysis showed that IFITM3 was significantly associated with cognitive performance (permutation‐corrected p = 1.25×10−3). Particularly, among two SNPs (rs10751647, rs2091850) in IFITM3 significantly associated with cognitive performance, rs10751647 was associated with cognitive decline in ADNI, which was replicated in AddNeuroMed. In addition, rs10751647 was significantly associated with amyloid‐β deposition measured by amyloid PET scan, CSF phosphorylated tau levels, and entorhinal cortical thickness measured by MRI scan in ADNI. The association of rs10751647 with entorhinal cortical thickness was replicated in AddNeuroMed. Participants with minor alleles (C) of rs10751647 have less cognitive decline, less amyloid and tau burden, and less brain atrophy. eQTL analysis showed that rs10751647 is associated with IFITM3 expression levels in blood and brain.ConclusionThis suggests that rs10751647 in IFITM3 is associated with less vulnerability for cognitive decline and AD biomarkers, providing mechanistic insight regarding involvement of immune activity and infection in AD.

Genomic Characterization of Arcobacter butzleri Strains Isolated from Various Sources in Lithuania.

Arcobacter (A.) butzleri, the most widespread species within the genus Arcobacter, is considered as an emerging pathogen causing gastroenteritis in humans. Here, we performed a comparative genome-wide analysis of 40 A. butzleri strains from Lithuania to determine the genetic relationship, pangenome structure, putative virulence, and potential antimicrobial- and heavy-metal-resistance genes. Core genome single nucleotide polymorphism (cgSNP) analysis revealed low within-group variability (≤4 SNPs) between three milk strains (RCM42, RCM65, RCM80) and one human strain (H19). Regardless of the type of input (i.e., cgSNPs, accessory genome, virulome, resistome), these strains showed a recurrent phylogenetic and hierarchical grouping pattern. A. butzleri demonstrated a relatively large and highly variable accessory genome (comprising of 6284 genes with around 50% of them identified as singletons) that only partially correlated to the isolation source. Downstream analysis of the genomes resulted in the detection of 115 putative antimicrobial- and heavy-metal-resistance genes and 136 potential virulence factors that are associated with the induction of infection in host (e.g., cadF, degP, iamA), survival and environmental adaptation (e.g., flagellar genes, CheA-CheY chemotaxis system, urease cluster). This study provides additional knowledge for a better A. butzleri-related risk assessment and highlights the need for further genomic epidemiology studies in Lithuania and other countries.

Computational analysis of non-synonymous single nucleotide polymorphism in the bovine PKLR gene

Pyruvate kinase (PKLR) is a potential candidate gene for milk production traits in cows. The main aim of this work is to investigate the potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the PKLR gene by using several computational tools. In silico tools including SIFT, Polyphen-2, SNAP2 and Panther indicated only 18 nsSNPs out of 170 were considered deleterious. The analysis of proteins’ stability change due to amino acid substitution performed by the use of the I-mutant, MUpro, CUPSTAT, SDM and Dynamut confirmed that 9 nsSNPs decreased protein stability. ConSurf analysis predicted that all 18 nsSNPs were evolutionary moderately or highly conserved. Two different domains of PKLR protein were revealed by the InterPro tool with 12 nsSNPs positioned in the Pyruvate Kinase barrel domain and 6 nsSNP present in the Pyruvate Kinase C Terminal. The PKLR 3D model was predicted by MODELLER software and validated via Ramachandran plot and Prosa which indicated a good quality model. The analysis of energy minimizations for the native and mutated structures was performed by SWISS PDB viewer with GROMOS 96 program and showed that 3 structural and 4 functional residues had total energy higher than the native model. These findings indicate that these mutant structures (rs441424814, rs449326723, rs476805413, rs472263384, rs474320860, rs475521477, rs441633284) were less stable than the native model. Molecular Dynamics simulations were performed to confirm the impact of nsSNPs on the protein structure and function. The present study provides useful information about functional SNPs that have an impact on PKLR protein in cattle. Communicated by Ramaswamy H. Sarma

FKBP5 genetic variants are associated with respiratory- and sleep-related parameters in Chinese patients with obstructive sleep apnea.

Obstructive sleep apnea (OSA) has been associated with psychiatric disorders, especially depression and posttraumatic stress disorder (PTSD). FKBP5 genetic variants have been previously reported to confer the risk of depression and PTSD. This study aimed to investigate the association of single nucleotide polymorphisms (SNPs) in the FKBP5 gene with OSA and OSA-related quantitative traits. Four SNPs within the FKBP5 gene (rs1360780, rs3800373, rs9296158, rs9470080) were genotyped in 5773 participants with anthropometric and polysomnography data. Linear regression and logistic regression analyses were performed to evaluate the relationship between FKBP5 SNPs and OSA-related traits. Binary logistic regression was used to assess the effect of SNPs on OSA susceptibility. Interacting genes of SNPs were assessed based on the 3DSNP database, and expression quantitative trait loci (eQTL) analysis for SNPs was adopted to examine the correlation of SNPs with gene expression. Gene expression analyses in human brains were performed with the aid of Brain Atlas. In moderate-to-severe OSA patients, all four SNPs were positively associated with AHIREM, and rs9296158 showed the strongest association (ß = 1.724, p = 0.001). Further stratified analyses showed that in men with moderate OSA, rs1360780, rs3800373 and rs9470080 were positively associated with wake time (p = 0.0267, p = 0.0254 and p = 0.0043, respectively). Rs1360780 and rs3800373 were 28 and 29.4%more likely to rate a higher ordered MAI category (OR (95% CI) = 1.280 (1.042 - 1.575), p = 0.019; OR (95% CI) = 1.294 (1.052 - 1.592), p = 0.015, respectively). Rs9296158 and rs9470080 increased the risk of low sleep efficiency by 25.7 and 28.1% (OR (95% CI) = 1.257 (1.003 - 1.575), p = 0.047; OR (95% CI) = 1.281 (1.026-1.6), p = 0.029, respectively). Integrated analysis of eQTL and gene expression patterns revealed that four SNPs may exert their effects by regulating FKBP5, TULP1, and ARMC12. Single nucleotide polymorphisms in the FKBP5 gene were associated with sleep respiratory events in moderate-to-severe OSA patients during REM sleep and associated with sleep architecture variables in men with moderate OSA. FKBP5 variants may be a potential predisposing factor for sleep disorders, especially in REM sleep.

Burkholderia contaminans Bacteriophage CSP3 Requires O-Antigen Polysaccharides for Infection.

The Burkholderia cepacia complex is a group of opportunistic pathogens that cause both severe acute and chronic respiratory infections. Due to their large genomes containing multiple intrinsic and acquired antimicrobial resistance mechanisms, treatment is often difficult and prolonged. One alternative to traditional antibiotics for treatment of bacterial infections is bacteriophages. Therefore, the characterization of bacteriophages infective for the Burkholderia cepacia complex is critical to determine their suitability for any future use. Here, we describe the isolation and characterization of novel phage, CSP3, infective against a clinical isolate of Burkholderia contaminans. CSP3 is a new member of the Lessievirus genus that targets various Burkholderia cepacia complex organisms. Single nucleotide polymorphism (SNP) analysis of CSP3-resistant B. contaminans showed that mutations to the O-antigen ligase gene, waaL, consequently inhibited CSP3 infection. This mutant phenotype is predicted to result in the loss of cell surface O-antigen, contrary to a related phage that requires the inner core of the lipopolysaccharide for infection. Additionally, liquid infection assays showed that CSP3 provides suppression of B. contaminans growth for up to 14 h. Despite the inclusion of genes that are typical of the phage lysogenic life cycle, we saw no evidence of CSP3's ability to lysogenize. Continuation of phage isolation and characterization is crucial in developing large and diverse phage banks for global usage in cases of antibiotic-resistant bacterial infections. IMPORTANCE Amid the global antibiotic resistance crisis, novel antimicrobials are needed to treat problematic bacterial infections, including those from the Burkholderia cepacia complex. One such alternative is the use of bacteriophages; however, a lot is still unknown about their biology. Bacteriophage characterization studies are of high importance for building phage banks, as future work in developing treatments such as phage cocktails should require well-characterized phages. Here, we report the isolation and characterization of a novel Burkholderia contaminans phage that requires the O-antigen for infection, a distinct phenotype seen among other related phages. Our findings presented in this article expand on the ever-evolving phage biology field, uncovering unique phage-host relationships and mechanisms of infection.

Whole-Genome Sequencing of Shiga Toxin-Producing Escherichia coli for Characterization and Outbreak Investigation.

Shiga toxin-producing Escherichia coli (STEC) causes high frequencies of foodborne infections worldwide and has been linked to numerous outbreaks each year. Pulsed-field gel electrophoresis (PFGE) has been the gold standard for surveillance until the recent transition to whole-genome sequencing (WGS). To further understand the genetic diversity and relatedness of outbreak isolates, a retrospective analysis of 510 clinical STEC isolates was conducted. Among the 34 STEC serogroups represented, most (59.6%) belonged to the predominant six non-O157 serogroups. Core genome single nucleotide polymorphism (SNP) analysis differentiated clusters of isolates with similar PFGE patterns and multilocus sequence types (STs). One serogroup O26 outbreak strain and another non-typeable (NT) strain, for instance, were identical by PFGE and clustered together by MLST; however, both were distantly related in the SNP analysis. In contrast, six outbreak-associated serogroup O5 strains clustered with five ST-175 serogroup O5 isolates, which were not part of the same outbreak as determined by PFGE. The use of high-quality SNP analyses enhanced the discrimination of these O5 outbreak strains into a single cluster. In all, this study demonstrates how public health laboratories can more rapidly use WGS and phylogenetics to identify related strains during outbreak investigations while simultaneously uncovering important genetic attributes that can inform treatment practices.

Development and testing of a polygenic risk score for breast cancer aggressiveness

Aggressive breast cancers portend a poor prognosis, but current polygenic risk scores (PRSs) for breast cancer do not reliably predict aggressive cancers. Aggressiveness can be effectively recapitulated using tumor gene expression profiling. Thus, we sought to develop a PRS for the risk of recurrence score weighted on proliferation (ROR-P), an established prognostic signature. Using 2363 breast cancers with tumor gene expression data and single nucleotide polymorphism (SNP) genotypes, we examined the associations between ROR-P and known breast cancer susceptibility SNPs using linear regression models. We constructed PRSs based on varying p-value thresholds and selected the optimal PRS based on model r2 in 5-fold cross-validation. We then used Cox proportional hazards regression to test the ROR-P PRS’s association with breast cancer-specific survival in two independent cohorts totaling 10,196 breast cancers and 785 events. In meta-analysis of these cohorts, higher ROR-P PRS was associated with worse survival, HR per SD = 1.13 (95% CI 1.06–1.21, p = 4.0 × 10–4). The ROR-P PRS had a similar magnitude of effect on survival as a comparator PRS for estrogen receptor (ER)-negative versus positive cancer risk (PRSER-/ER+). Furthermore, its effect was minimally attenuated when adjusted for PRSER-/ER+, suggesting that the ROR-P PRS provides additional prognostic information beyond ER status. In summary, we used integrated analysis of germline SNP and tumor gene expression data to construct a PRS associated with aggressive tumor biology and worse survival. These findings could potentially enhance risk stratification for breast cancer screening and prevention.

Genetic Characterization of the Norwegian Apple Collection

Commercial fruit production in Norway is located at around latitude 60° north, demanding a careful choice of adapted cultivars. The most comprehensive collection of apple genetic resources in Norway is being kept in the Norwegian Apple Collection (NAC) at the Njøs Fruit and Berry Centre (NJØS). The collection contains around 350 accessions and was recently genotyped with a single nucleotide polymorphism (SNP) array. Curated SNP data were used for the assessment of structure and diversity, pedigree confirmation, and core collection development. In the following SNP analysis, we identified several duplicates and parent-child relationships. Across the geographic regions represented, the collection was equally diverse. Different methods for analyzing population structure were applied. K-means clustering and a Bayesian modeling approach with prior assumptions of the data revealed five subpopulations associated with geographic breeding centers. The collection has a distinct genetic structure and low relatedness among the accessions; hence, two core collections with 100 accessions in each were created. These new core collections will allow breeders and researchers to use the NAC efficiently. The results from this study suggest that several of the accessions in the Norwegian Apple Collection could be of high importance for breeding purposes.

Investigation of a Staphylococcus aureus sequence type 72 food poisoning outbreak associated with food-handler contamination in Italy.

On August 2019 a staphylococcal food poisoning outbreak occurred in an elderly home in Piedmont, Italy. The epidemiological investigation performed among the persons that consumed the meal identified chicken salad as the most likely source of the outbreak. Staphylococcus aureus was isolated from a total of seven samples, namely one vomit sample from a guest of the nursing home, two food samples (chicken salad with and without mayonnaise) and nasal swabs collected from a total of four persons working in the kitchen of the nursing home. The maximum likelihood tree obtained using single nucleotide polymorphisms analysis revealed that the isolates from the aforementioned samples clustered together. Multilocus sequence typing revealed that they belonged to Sequence Type 72. Fourier transform infrared spectroscopy (FTIR) was used in parallel to single nucleotide polymorphisms and whole genome sequencing for the determination of the degree of relatedness of the isolates. The results of the FTIR showed the same clustering obtained with single nucleotide polymorphisms and whole genome sequencing and revealed the source of infection. This study underlines the importance of both laboratory evidence and epidemiological data for outbreak investigation and further confirms that FTIR is a suitable support for the short-term epidemiological investigation on source attribution in case of a S. aureus infection.

Analysis of single nucleotide polymorphisms associated with the vitamin D pathway in the placentas of women with gestational diabetes mellitus: a laboratory study.

The aim of this study was to analyze the single nucleotide polymorphisms (SNPs) of genes known to be involved in vitamin D metabolism in the placenta using the placental tissue of mothers diagnosed with gestational diabetes mellitus (GDM) to determine whether the SNPs and occurrence of GDM are related. We enrolled 80 women of the same gestational age, 40 with and 40 without GDM. The placenta was obtained from each woman after delivery and SNP genotyping was performed on seven SNPs in the CYP27B1 (rs10877012), CYP24A1 (rs2248359, rs6013897, and rs2209314), and GC (rs2282679, rs16847024, and rs3733359) genes. Maternal serum 25-hydroxyvitamin D levels were measured during the first trimester of pregnancy and before delivery. At the time of delivery, vitamin D levels were lower (21.05±12.05 mg/dL vs. 31.31±20.72 mg/dL, p=0.012) and the frequency of vitamin D deficiency was higher (60.7% vs. 32.5%, p=0.040) in the GDM group. In women with GDM, the G allele of rs10877012 was more common (86.3% vs. 65.0%, p=0.002). The rs10877012 GG genotype was more common in the GDM group (72.5% vs. 42.5%, p=0.007) and the rs10877012 TT genotype was more common in the control group (12.5% vs. 0%, p=0.007). Mothers with GDM have lower serum concentrations of vitamin D before delivery than healthy controls and vitamin D deficiency is common. A polymorphism in CYP27B1 (rs10877012), is considered to be a cause of GDM pathogenesis.

Comparative Genomic Analyses of Lactococcus garvieae Isolated from Bovine Mastitis in China.

Lactococcus garvieae is an emerging zoonotic pathogen, but there are few reports regarding bovine mastitis. The prevalence of L. garvieae represents an increasing disease threat and global public health risk. Thirty-nine L. garvieae isolates were obtained from 2,899 bovine clinical mastitis milk samples in 6 provinces of China from 2017 to 2021. Five clonal complexes were determined from 32 multilocus sequence types (MLSTs) of L. garvieae: sequence type 46 (ST46) was the predominant sequence type, and 13 novel MLSTs were identified. All isolates were resistant to chloramphenicol and clindamycin, but susceptible to penicillin, ampicillin, amoxicillin-clavulanic acid, imipenem, ceftiofur, enrofloxacin, and marbofloxacin. Based on genomic analyses, L. garvieae had 6,310 genes, including 1,015 core, 3,641 accessory, and 1,654 unique genes. All isolates had virulence genes coding for collagenase, fibronectin-binding protein, glyceraldehyde-3-phosphate dehydrogenase, superoxide dismutase, and NADH oxidase. Most isolates had lsaD and mdtA antimicrobial resistance (AMR) genes. Based on COG (Clusters of Orthologous Genes database) results, the functions of defense, transcription and replication, and recombination and repair were enhanced in unique genes, whereas functions of translation, ribosomal structure, and biogenesis were enhanced in core genes. The KEGG functional categories enriched in unique genes included human disease and membrane transport, whereas COG functional categories enriched in core genes included energy metabolism, nucleotide metabolism, and translation. No gene was significantly associated with host specificity. In addition, analysis of core genome single nucleotide polymorphisms (SNPs) implied potential host adaptation of some isolates in several sequence types. In conclusion, this study characterized L. garvieae isolated from mastitis and detected potential adaptations of L. garvieae to various hosts. IMPORTANCE This study provides important genomic insights into a bovine mastitis pathogen, Lactococcus garvieae. Comprehensive genomic analyses of L. garvieae from dairy farms have not been reported. This study is a detailed and comprehensive report of novel features of isolates of L. garvieae, an important but poorly characterized bacterium, recovered in the past 5 years in 6 Chinese provinces. We documented diverse genetic features, including predominant sequence type ST46 and 13 novel MLSTs. Lactococcus garvieae had 6,310 genes, including 1,015 core, 3,641 accessory, and 1,654 unique genes. All isolates had virulence genes coding for collagenase, fibronectin-binding protein, glyceraldehyde-3-phosphate dehydrogenase, superoxide dismutase, and NADH oxidase and resistance to chloramphenicol and clindamycin. Most isolates had lsaD and mdtA antimicrobial resistance genes. However, no gene was significantly associated with host specificity. This is the first report that characterized L. garvieae isolates from bovine mastitis and revealed potential host adaptations of L. garvieae to various hosts.