Consideration of within-patient diversity highlights transmission pathways and antimicrobial resistance gene variability in vancomycin-resistant Enterococcus faecium.

BACKGROUND: WGS is increasingly being applied to healthcare-associated vancomycin-resistant Enterococcus faecium (VREfm) outbreaks. Within-patient diversity could complicate transmission resolution if single colonies are sequenced from identified cases. OBJECTIVES: Determine the impact of within-patient diversity on transmission resolution of VREfm. MATERIALS AND METHODS: Fourteen colonies were collected from VREfm positive rectal screens, single colonies were collected from clinical samples and Illumina WGS was performed. Two isolates were selected for Oxford Nanopore sequencing and hybrid genome assembly to generate lineage-specific reference genomes. Mapping to closely related references was used to identify genetic variations and closely related genomes. A transmission network was inferred for the entire genome set using Phyloscanner. RESULTS AND DISCUSSION: In total, 229 isolates from 11 patients were sequenced. Carriage of two or three sequence types was detected in 27% of patients. Presence of antimicrobial resistance genes and plasmids was variable within genomes from the same patient and sequence type. We identified two dominant sequence types (ST80 and ST1424), with two putative transmission clusters of two patients within ST80, and a single cluster of six patients within ST1424. We found transmission resolution was impaired using fewer than 14 colonies. CONCLUSIONS: Patients can carry multiple sequence types of VREfm, and even within related lineages the presence of mobile genetic elements and antimicrobial resistance genes can vary. VREfm within-patient diversity could be considered in future to aid accurate resolution of transmission networks.

The genus Serratia revisited by genomics.

The genus Serratia has been studied for over a century and includes clinically-important and diverse environmental members. Despite this, there is a paucity of genomic information across the genus and a robust whole genome-based phylogenetic framework is lacking. Here, we have assembled and analysed a representative set of 664 genomes from across the genus, including 215 historic isolates originally used in defining the genus. Phylogenomic analysis of the genus reveals a clearly-defined population structure which displays deep divisions and aligns with ecological niche, as well as striking congruence between historical biochemical phenotyping data and contemporary genomics data. We highlight the genomic, phenotypic and plasmid diversity of Serratia, and provide evidence of different patterns of gene flow across the genus. Our work provides a framework for understanding the emergence of clinical and other lineages of Serratia.

Presence of optrA-mediated linezolid resistance in multiple lineages and plasmids of Enterococcus faecalis revealed by long read sequencing.

Transferable linezolid resistance due to optrA, poxtA, cfr and cfr-like genes is increasingly detected in enterococci associated with animals and humans globally. We aimed to characterize the genetic environment of optrA in linezolid-resistant Enterococcus faecalis isolates from Scotland. Six linezolid-resistant E. faecalis isolated from urogenital samples were confirmed to carry the optrA gene by PCR. Short read (Illumina) sequencing showed the isolates were genetically distinct (>13900 core SNPs) and belonged to different MLST sequence types. Plasmid contents were examined using hybrid assembly of short and long read (Oxford Nanopore MinION) sequencing technologies. The optrA gene was located on distinct plasmids in each isolate, suggesting that transfer of a single plasmid did not contribute to optrA dissemination in this collection. pTM6294-2, BX5936-1 and pWE0438-1 were similar to optrA-positive plasmids from China and Japan, while the remaining three plasmids had limited similarity to other published examples. We identified the novel Tn6993 transposon in pWE0254-1 carrying linezolid (optrA), macrolide (ermB) and spectinomycin [ANT(9)-Ia] resistance genes. OptrA amino acid sequences differed by 0-20 residues. We report multiple variants of optrA on distinct plasmids in diverse strains of E. faecalis. It is important to identify the selection pressures driving the emergence and maintenance of resistance against linezolid to retain the clinical utility of this antibiotic.

Pan-Resistome Characterization of Uropathogenic Escherichia coli and Klebsiella pneumoniae Strains Circulating in Uganda and Kenya, Isolated from 2017-2018.

Urinary tract infection (UTI) develops after a pathogen adheres to the inner lining of the urinary tract. Cases of UTIs are predominantly caused by several Gram-negative bacteria and account for high morbidity in the clinical and community settings. Of greater concern are the strains carrying antimicrobial resistance (AMR)-conferring genes. The gravity of a UTI is also determined by a spectrum of other virulence factors. This study represents a pilot project to investigate the burden of AMR among uropathogens in East Africa. We examined bacterial samples isolated in 2017-2018 from in- and out-patients in Kenya (KY) and Uganda (UG) that presented with clinical symptoms of UTI. We reconstructed the evolutionary history of the strains, investigated their population structure, and performed comparative analysis their pangenome contents. We found 55 Escherichia coli and 19 Klebsiella pneumoniae strains confirmed uropathogenic following screening for the prevalence of UTI virulence genes including fimH, iutA, feoA/B/C, mrkD, and foc. We identified 18 different sequence types in E. coli population while all K. pneumoniae strains belong to ST11. The most prevalent E. coli sequence types were ST131 (26%), ST335/1193 (10%), and ST10 (6%). Diverse plasmid types were observed in both collections such as Incompatibility (IncF/IncH/IncQ1/IncX4) and Col groups. Pangenome analysis of each set revealed a total of 2862 and 3464 genes comprised the core genome of E. coli and K. pneumoniae population, respectively. Among these are acquired AMR determinants including fluoroquinolone resistance-conferring genes aac(3)-Ib-cr and other significant genes: aad, tet, sul1, sul2, and cat, which are associated with aminoglycoside, tetracycline, sulfonamide, and chloramphenicol resistance, respectively. Accessory genomes of both species collections were detected several ß-lactamase genes, blaCTX-M, blaTEM and blaOXA, or blaNDM. Overall, 93% are multi-drug resistant in the E. coli collection while 100% of the K. pneumoniae strains contained genes that are associated with resistance to three or more antibiotic classes. Our findings illustrate the abundant acquired resistome and virulome repertoire in uropathogenic E. coli and K. pneumoniae, which are mainly disseminated via clonal and horizontal transfer, circulating in the East African region. We further demonstrate here that routine genomic surveillance is necessary for high-resolution bacterial epidemiology of these important AMR pathogens.

The Microevolution and Epidemiology of Staphylococcus aureus Colonization during Atopic Eczema Disease Flare.

Staphylococcus aureus is an opportunistic pathogen and variable component of the human microbiota. A characteristic of atopic eczema (AE) is colonization by S. aureus, with exacerbations associated with an increased bacterial burden of the organism. Despite this, the origins and genetic diversity of S. aureus colonizing individual patients during AE disease flares is poorly understood. To examine the microevolution of S. aureus colonization, we deep sequenced S. aureus populations from nine children with moderate to severe AE and 18 non-atopic children asymptomatically carrying S. aureus nasally. Colonization by clonal S. aureus populations was observed in both AE patients and control participants, with all but one of the individuals carrying colonies belonging to a single sequence type. Phylogenetic analysis showed that disease flares were associated with the clonal expansion of the S. aureus population, occurring over a period of weeks to months. There was a significant difference in the genetic backgrounds of S. aureus colonizing AE cases versus controls (Fisher exact test, P = 0.03). Examination of intra-host genetic heterogeneity of the colonizing S. aureus populations identified evidence of within-host selection in the AE patients, with AE variants being potentially selectively advantageous for intracellular persistence and treatment resistance.

MRSA Transmission Dynamics Among Interconnected Acute, Intermediate-Term, and Long-Term Healthcare Facilities in Singapore.

BACKGROUND: Methicillin-resistant Staphylococcus aureus (MRSA) is the most common healthcare-associated multidrug-resistant organism. Despite the interconnectedness between acute care hospitals (ACHs) and intermediate- and long-term care facilities (ILTCFs), the transmission dynamics of MRSA between healthcare settings is not well understood. METHODS: We conducted a cross-sectional study in a network comprising an ACH and 5 closely affiliated ILTCFs in Singapore. A total of 1700 inpatients were screened for MRSA over a 6-week period in 2014. MRSA isolates underwent whole-genome sequencing, with a pairwise single-nucleotide polymorphism (Hamming distance) cutoff of 60 core genome single-nucleotide polymorphisms used to define recent transmission clusters (clades) for the 3 major clones. RESULTS: MRSA prevalence was significantly higher in intermediate-term (29.9%) and long-term (20.4%) care facilities than in the ACH (11.8%) (P < .001). The predominant clones were sequence type [ST] 22 (n = 183; 47.8%), ST45 (n = 129; 33.7%), and ST239 (n = 26; 6.8%), with greater diversity of STs in ILTCFs relative to the ACH. A large proportion of the clades in ST22 (14 of 21 clades; 67%) and ST45 (7 of 13; 54%) included inpatients from the ACH and ILTCFs. The most frequent source of the interfacility transmissions was the ACH (n = 28 transmission events; 36.4%). CONCLUSIONS: MRSA transmission dynamics between the ACH and ILTCFs were complex. The greater diversity of STs in ILTCFs suggests that the ecosystem in such settings might be more conducive for intrafacility transmission events. ST22 and ST45 have successfully established themselves in ILTCFs. The importance of interconnected infection prevention and control measures and strategies cannot be overemphasized.

Further evidence for a parent-of-origin effect at the NOP9 locus on language-related phenotypes.

BACKGROUND: Specific language impairment (SLI) is a common neurodevelopmental disorder, observed in 5-10 % of children. Family and twin studies suggest a strong genetic component, but relatively few candidate genes have been reported to date. A recent genome-wide association study (GWAS) described the first statistically significant association specifically for a SLI cohort between a missense variant (rs4280164) in the NOP9 gene and language-related phenotypes under a parent-of-origin model. Replications of these findings are particularly challenging because the availability of parental DNA is required. METHODS: We used two independent family-based cohorts characterised with reading- and language-related traits: a longitudinal cohort (n = 106 informative families) including children with language and reading difficulties and a nuclear family cohort (n = 264 families) selected for dyslexia. RESULTS: We observed association with language-related measures when modelling for parent-of-origin effects at the NOP9 locus in both cohorts: minimum P = 0.001 for phonological awareness with a paternal effect in the first cohort and minimum P = 0.0004 for irregular word reading with a maternal effect in the second cohort. Allelic and parental trends were not consistent when compared to the original study. CONCLUSIONS: A parent-of-origin effect at this locus was detected in both cohorts, albeit with different trends. These findings contribute in interpreting the original GWAS report and support further investigations of the NOP9 locus and its role in language-related traits. A systematic evaluation of parent-of-origin effects in genetic association studies has the potential to reveal novel mechanisms underlying complex traits.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts.

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10(-8)) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left-right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

Copy Number Variation Screen Identifies a Rare De Novo Deletion at Chromosome 15q13.1-13.3 in a Child with Language Impairment.

A significant proportion of children (up to 7% in the UK) present with pronounced language difficulties that cannot be explained by obvious causes like other neurological and medical conditions. A substantial genetic component is predicted to underlie such language problems. Copy number variants (CNVs) have been implicated in neurodevelopmental and psychiatric conditions, such as autism and schizophrenia, but it is not fully established to what extent they might contribute to language disorders. We conducted a CNV screen in a longitudinal cohort of young children with language-related difficulties (n = 85), focusing on single events at candidate loci. We detected a de novo deletion on chromosome 15q13.1-13.3. The adjacent 15q11-13.1 locus is disrupted in Prader-Willi and Angelman syndromes, while disruptions across the breakpoints (BP1-BP6) have previously been implicated in different neurodevelopmental phenotypes including autism, intellectual disability (ID), seizures and developmental delay (DD). This is the first report of a deletion at BP3-BP5 being linked to a deficit confined to language impairment, in the absence of ID, expanding the range of phenotypes that implicate the chromosome 15q13 locus.

Differential TERT promoter methylation and response to 5-aza-2'-deoxycytidine in acute myeloid leukemia cell lines: TERT expression, telomerase activity, telomere length, and cell death.

The catalytic subunit of human telomerase (TERT) is highly expressed in cancer cells, and correlates with complex cytogenetics and disease severity in acute myeloid leukemia (AML). The TERT promoter is situated within a large CpG island, suggesting that expression is methylation-sensitive. Studies suggest a correlation between hypermethylation and TERT overexpression. We investigated the relationship between TERT promoter methylation and expression and telomerase activity in human leukemia and lymphoma cell lines. DAC-induced demethylation and cell death were observed in all three cell lines, as well as telomere shortening in HL-60 cells. DAC treatment reduced TERT expression and telomerase activity in OCI/AML3 and HL-60 cells, but not in U937 cells. Control U937 cells expressed lower levels of TERT mRNA, carried a highly methylated TERT core promoter, and proved more resistant to DAC-induced repression of TERT expression and cell death. AML patients had significantly lower methylation levels at several CpGs than "well elderly" individuals. This study, the first to investigate the relationship between TERT methylation and telomerase activity in leukemia cells, demonstrated a differential methylation pattern and response to DAC in three AML cell lines. We suggest that, although DAC treatment reduces TERT expression and telomerase activity, this is unlikely to occur via direct demethylation of the TERT promoter. However, further investigations on the regions spanning CpGs 7-12 and 14-16 may reveal valuable information regarding transcriptional regulation of TERT.

Implications of discoveries from genome-wide association studies in current cardiovascular practice.

Genome-wide association studies (GWAS) have identified several genetic variants associated with coronary heart disease (CHD), and variations in plasma lipoproteins and blood pressure (BP). Loci corresponding to CDKN2A/CDKN2B/ANRIL, MTHFD1L, CELSR2, PSRC1 and SORT1 genes have been associated with CHD, and TMEM57, DOCK7, CELSR2, APOB, ABCG5, HMGCR, TRIB1, FADS2/S3, LDLR, NCAN and TOMM40-APOE with total cholesterol. Similarly, CELSR2-PSRC1-SORT1, PCSK9, APOB, HMGCR, NCAN-CILP2-PBX4, LDLR, TOMM40-APOE, and APOC1-APOE are associated with variations in low-density lipoprotein cholesterol levels. Altogether, forty, forty three and twenty loci have been associated with high-density lipoprotein cholesterol, triglycerides and BP phenotypes, respectively. Some of these identified loci are common for all the traits, some do not map to functional genes, and some are located in genes that encode for proteins not previously known to be involved in the biological pathway of the trait. GWAS have been successful at identifying new and unexpected genetic loci common to diseases and traits, thus rapidly providing key novel insights into disease biology. Since genotype information is fixed, with minimum biological variability, it is useful in early life risk prediction. However, these variants explain only a small proportion of the observed variance of these traits. Therefore, the utility of genetic determinants in assessing risk at later stages of life has limited immediate clinical impact. The future application of genetic screening will be in identifying risk groups early in life to direct targeted preventive measures.

Fine mapping the KLK3 locus on chromosome 19q13.33 associated with prostate cancer susceptibility and PSA levels.

Measurements of serum prostate-specific antigen (PSA) protein levels form the basis for a widely used test to screen men for prostate cancer. Germline variants in the gene that encodes the PSA protein (KLK3) have been shown to be associated with both serum PSA levels and prostate cancer. Based on a resequencing analysis of a 56 kb region on chromosome 19q13.33, centered on the KLK3 gene, we fine mapped this locus by genotyping tag SNPs in 3,522 prostate cancer cases and 3,338 controls from five case-control studies. We did not observe a strong association with the KLK3 variant, reported in previous studies to confer risk for prostate cancer (rs2735839; P = 0.20) but did observe three highly correlated SNPs (rs17632542, rs62113212 and rs62113214) associated with prostate cancer [P = 3.41 × 10(-4), per-allele trend odds ratio (OR) = 0.77, 95% CI = 0.67-0.89]. The signal was apparent only for nonaggressive prostate cancer cases with Gleason score <7 and disease stage 8 or stage ≥III (P = 0.31, per-allele trend OR = 1.12, 95% CI = 0.90-1.40). One of the three highly correlated SNPs, rs17632542, introduces a non-synonymous amino acid change in the KLK3 protein with a predicted benign or neutral functional impact. Baseline PSA levels were 43.7% higher in control subjects with no minor alleles (1.61 ng/ml, 95% CI = 1.49-1.72) than in those with one or more minor alleles at any one of the three SNPs (1.12 ng/ml, 95% CI = 0.96-1.28) (P = 9.70 × 10(-5)). Together our results suggest that germline KLK3 variants could influence the diagnosis of nonaggressive prostate cancer by influencing the likelihood of biopsy.

Resequencing of the CCL5 and CCR5 genes and investigation of variants for association with diabetic nephropathy.

Chemokine (C-C motif) ligand 5 (CCL5) and chemokine (C-C motif) receptor 5 are implicated in the pathogenesis of diabetic nephropathy (DN). We hypothesize that variants in these genes may be associated with DN. The CCL5 and chemokine receptor type 5 (CCR5) genes were resequenced, variants identified (n=58), allele frequencies determined in 46 individuals (92 chromosomes) and efficient haplotype tag single-nucleotide polymorphisms (htSNPs) selected to effectively evaluate the common variation in these genes. One reportedly functional gene variant and eight htSNPs were genotyped in a case-control association study involving Caucasian individuals with type 1 diabetes (267 cases with DN and 442 non-nephropathic diabetic controls). Genotyping was performed using MassARRAY iPLEX, TaqMan, gel electrophoresis and direct capillary sequencing. After correction for multiple testing, there were no statistically significant associations between variants in the CCL5 and CCR5 genes and DN.

A GREM1 gene variant associates with diabetic nephropathy.

Gremlin, a cell growth and differentiation factor, promotes the development of diabetic nephropathy in animal models, but whether GREM1 gene variants associate with diabetic nephropathy is unknown. We comprehensively screened the 5' upstream region (including the predicted promoter), all exons, intron-exon boundaries, complete untranslated regions, and the 3' region downstream of the GREM1 gene. We identified 31 unique variants, including 24 with a minor allele frequency exceeding 5%, and 9 haplotype-tagging single nucleotide polymorphisms (htSNPs). We selected one additional variant that we predicted to alter transcription factor binding. We genotyped 709 individuals with type 1 diabetes of whom 267 had nephropathy (cases) and 442 had no evidence of kidney disease (controls). Three individual SNPs significantly associated with nephropathy at the 5% level, and two remained significant after adjustment for multiple testing. Subsequently, we genotyped a replicate population comprising 597 cases and 502 controls: this population supported an association with one of the SNPs (rs1129456; P = 0.0003). Combined analysis, adjusted for recruitment center (n = 8), suggested that the T allele conferred greater odds of nephropathy (OR 1.69; 95% CI 1.36 to 2.11). In summary, the GREM1 variant rs1129456 associates with diabetic nephropathy, perhaps explaining some of the genetic susceptibility to this condition.

Vascular endothelial growth factor promotes physical wound repair and is anti-apoptotic in primary distal lung epithelial and A549 cells.

OBJECTIVE: There is evidence to suggest a beneficial role for growth factors, including vascular endothelial growth factor (VEGF), in tissue repair and proliferation after injury within the lung. Whether this effect is mediated predominantly by actions on endothelial cells or epithelial cells is unknown. This study tested the hypothesis that VEGF acts as an autocrine trophic factor for human adult alveolar epithelial cells and that under situations of pro-apoptotic stress, VEGF reduces cell death. DESIGN: In vitro cell culture study looking at the effects of 0.03% H2O2 on both A549 and primary distal lung epithelial cells. MEASUREMENT AND MAIN RESULTS: Primary adult human distal lung epithelial cells express both the soluble and membrane-associated VEGF isoforms and VEGF receptors 1 and 2. At physiologically relevant doses, soluble VEGF isoforms stimulate wound repair and have a proliferative action. Specific receptor ligands confirmed that this effect was mediated by VEGF receptor 1. In addition to proliferation, we demonstrate that VEGF reduces A549 and distal lung epithelial cell apoptosis when administered after 0.03% H2O2 injury. This effect occurs due to reduced caspase-3 activation and is phosphatidylinositol 3'-kinase dependent. CONCLUSION: In addition to its known effects on endothelial cells, VEGF acts as a growth and anti-apoptotic factor on alveolar epithelial cells. VEGF treatment may have potential as a rescue therapy for diseases associated with alveolar epithelial damage such as acute respiratory distress syndrome.