Exploring the evolution and epidemiology of European CC1-MRSA-IV: tracking a multidrug-resistant community-associated meticillin-resistant Staphylococcus aureus clone.

This study investigated the evolution and epidemiology of the community-associated and multidrug-resistant Staphylococcus aureus clone European CC1-MRSA-IV. Whole-genome sequences were obtained for 194 European CC1-MRSA-IV isolates (189 of human and 5 of animal origin) from 12 countries, and 10 meticillin-susceptible precursors (from North-Eastern Romania; all of human origin) of the clone. Phylogenetic analysis was performed using a maximum-likelihood approach, a time-measured phylogeny was reconstructed using Bayesian analysis, and in silico microarray genotyping was performed to identify resistance, virulence-associated and SCCmec (staphylococcal cassette chromosome mec) genes. Isolates were typically sequence type 1 (190/204) and spa type t127 (183/204). Bayesian analysis indicated that European CC1-MRSA-IV emerged in approximately 1995 before undergoing rapid expansion in the late 1990s and 2000s, while spreading throughout Europe and into the Middle East. Phylogenetic analysis revealed an unstructured meticillin-resistant S. aureus (MRSA) population, lacking significant geographical or temporal clusters. The MRSA were genotypically multidrug-resistant, consistently encoded seh, and intermittently (34/194) encoded an undisrupted hlb gene with concomitant absence of the lysogenic phage-encoded genes sak and scn. All MRSA also harboured a characteristic ~5350 nt insertion in SCCmec adjacent to orfX. Detailed demographic data from Denmark showed that there, the clone is typically (25/35) found in the community, and often (10/35) among individuals with links to South-Eastern Europe. This study elucidated the evolution and epidemiology of European CC1-MRSA-IV, which emerged from a meticillin-susceptible lineage prevalent in North-Eastern Romania before disseminating rapidly throughout Europe.

Taking the next-gen step: Comprehensive antimicrobial resistance detection from Burkholderia pseudomallei.

BACKGROUND: Antimicrobial resistance (AMR) poses a major threat to human health. Whole-genome sequencing holds great potential for AMR identification; however, there remain major gaps in accurately and comprehensively detecting AMR across the spectrum of AMR-conferring determinants and pathogens. METHODS: Using 16 wild-type Burkholderia pseudomallei and 25 with acquired AMR, we first assessed the performance of existing AMR software (ARIBA, CARD, ResFinder, and AMRFinderPlus) for detecting clinically relevant AMR in this pathogen. B. pseudomallei was chosen due to limited treatment options, high fatality rate, and AMR caused exclusively by chromosomal mutation (i.e. single-nucleotide polymorphisms [SNPs], insertions-deletions [indels], copy-number variations [CNVs], inversions, and functional gene loss). Due to poor performance with existing tools, we developed ARDaP (Antimicrobial Resistance Detection and Prediction) to identify the spectrum of AMR-conferring determinants in B. pseudomallei. FINDINGS: CARD, ResFinder, and AMRFinderPlus failed to identify any clinically-relevant AMR in B. pseudomallei; ARIBA identified AMR encoded by SNPs and indels that were manually added to its database. However, none of these tools identified CNV, inversion, or gene loss determinants, and ARIBA could not differentiate AMR determinants from natural genetic variation. In contrast, ARDaP accurately detected all SNP, indel, CNV, inversion, and gene loss AMR determinants described in B. pseudomallei (n≈50). Additionally, ARDaP accurately predicted three previously undescribed determinants. In mixed strain data, ARDaP identified AMR to as low as ~5% allelic frequency. INTERPRETATION: Existing AMR software packages are inadequate for chromosomal AMR detection due to an inability to detect resistance conferred by CNVs, inversions, and functional gene loss. ARDaP overcomes these major shortcomings. Further, ARDaP enables AMR prediction from mixed sequence data down to 5% allelic frequency, and can differentiate natural genetic variation from AMR determinants. ARDaP databases can be constructed for any microbial species of interest for comprehensive AMR detection. FUNDING: National Health and Medical Research Council (BJC, EPP, DSS); Australian Government (DEM, ES); Advance Queensland (EPP, DSS).

Development and validation of a RAD-Seq target-capture based genotyping assay for routine application in advanced black tiger shrimp (Penaeus monodon) breeding programs.

BACKGROUND: The development of genome-wide genotyping resources has provided terrestrial livestock and crop industries with the unique ability to accurately assess genomic relationships between individuals, uncover the genetic architecture of commercial traits, as well as identify superior individuals for selection based on their specific genetic profile. Utilising recent advancements in de-novo genome-wide genotyping technologies, it is now possible to provide aquaculture industries with these same important genotyping resources, even in the absence of existing genome assemblies. Here, we present the development of a genome-wide SNP assay for the Black Tiger shrimp (Penaeus monodon) through utilisation of a reduced-representation whole-genome genotyping approach (DArTseq). RESULTS: Based on a single reduced-representation library, 31,262 polymorphic SNPs were identified across 650 individuals obtained from Australian wild stocks and commercial aquaculture populations. After filtering to remove SNPs with low read depth, low MAF, low call rate, deviation from HWE, and non-Mendelian inheritance, 7542 high-quality SNPs were retained. From these, 4236 high-quality genome-wide loci were selected for baits-probe development and 4194 SNPs were included within a finalized target-capture genotype-by-sequence assay (DArTcap). This assay was designed for routine and cost effective commercial application in large scale breeding programs, and demonstrates higher confidence in genotype calls through increased call rate (from 80.2 ± 14.7 to 93.0% ± 3.5%), increased read depth (from 20.4 ± 15.6 to 80.0 ± 88.7), as well as a 3-fold reduction in cost over traditional genotype-by-sequencing approaches. CONCLUSION: Importantly, this assay equips the P. monodon industry with the ability to simultaneously assign parentage of communally reared animals, undertake genomic relationship analysis, manage mate pairings between cryptic family lines, as well as undertake advance studies of genome and trait architecture. Critically this assay can be cost effectively applied as P. monodon breeding programs transition to undertaking genomic selection.

Evolution and Global Transmission of a Multidrug-Resistant, Community-Associated Methicillin-Resistant Staphylococcus aureus Lineage from the Indian Subcontinent.

The evolution and global transmission of antimicrobial resistance have been well documented for Gram-negative bacteria and health care-associated epidemic pathogens, often emerging from regions with heavy antimicrobial use. However, the degree to which similar processes occur with Gram-positive bacteria in the community setting is less well understood. In this study, we traced the recent origins and global spread of a multidrug-resistant, community-associated Staphylococcus aureus lineage from the Indian subcontinent, the Bengal Bay clone (ST772). We generated whole-genome sequence data of 340 isolates from 14 countries, including the first isolates from Bangladesh and India, to reconstruct the evolutionary history and genomic epidemiology of the lineage. Our data show that the clone emerged on the Indian subcontinent in the early 1960s and disseminated rapidly in the 1990s. Short-term outbreaks in community and health care settings occurred following intercontinental transmission, typically associated with travel and family contacts on the subcontinent, but ongoing endemic transmission was uncommon. Acquisition of a multidrug resistance integrated plasmid was instrumental in the emergence of a single dominant and globally disseminated clade in the early 1990s. Phenotypic data on biofilm, growth, and toxicity point to antimicrobial resistance as the driving force in the evolution of ST772. The Bengal Bay clone therefore combines the multidrug resistance of traditional health care-associated clones with the epidemiological transmission of community-associated methicillin-resistant S. aureus (MRSA). Our study demonstrates the importance of whole-genome sequencing for tracking the evolution of emerging and resistant pathogens. It provides a critical framework for ongoing surveillance of the clone on the Indian subcontinent and elsewhere.IMPORTANCE The Bengal Bay clone (ST772) is a community-associated and multidrug-resistant Staphylococcus aureus lineage first isolated from Bangladesh and India in 2004. In this study, we showed that the Bengal Bay clone emerged from a virulent progenitor circulating on the Indian subcontinent. Its subsequent global transmission was associated with travel or family contact in the region. ST772 progressively acquired specific resistance elements at limited cost to its fitness and continues to be exported globally, resulting in small-scale community and health care outbreaks. The Bengal Bay clone therefore combines the virulence potential and epidemiology of community-associated clones with the multidrug resistance of health care-associated S. aureus lineages. This study demonstrates the importance of whole-genome sequencing for the surveillance of highly antibiotic-resistant pathogens, which may emerge in the community setting of regions with poor antibiotic stewardship and rapidly spread into hospitals and communities across the world.

Global Scale Dissemination of ST93: A Divergent Staphylococcus aureus Epidemic Lineage That Has Recently Emerged From Remote Northern Australia.

Background: In Australia, community-associated methicillin-resistant Staphylococcus aureus (MRSA) lineage sequence type (ST) 93 has rapidly risen to dominance since being described in the early 1990s. We examined 459 ST93 genome sequences from Australia, New Zealand, Samoa, and Europe to investigate the evolutionary history of ST93, its emergence in Australia and subsequent spread overseas. Results: Comparisons with other S. aureus genomes indicate that ST93 is an early diverging and recombinant lineage, comprising of segments from the ST59/ST121 lineage and from a divergent but currently unsampled Staphylococcal population. However, within extant ST93 strains limited genetic diversity was apparent with the most recent common ancestor dated to 1977 (95% highest posterior density 1973-1981). An epidemic ST93 population arose from a methicillin-susceptible progenitor in remote Northern Australia, which has a proportionally large Indigenous population, with documented overcrowded housing and a high burden of skin infection. Methicillin-resistance was acquired three times in these regions, with a clade harboring a staphylococcal cassette chromosome mec (SCCmec) IVa expanding and spreading to Australia's east coast by 2000. We observed sporadic and non-sustained introductions of ST93-MRSA-IVa to the United Kingdom. In contrast, in New Zealand, ST93-MRSA-IVa was sustainably transmitted with clonal expansion within the Pacific Islander population, who experience similar disadvantages as Australian Indigenous populations. Conclusion: ST93 has a highly recombinant genome including portions derived from an early diverging S. aureus population. Our findings highlight the need to understand host population factors in the emergence and spread of antimicrobial resistant community pathogens.

Identification of key contributors in complex population structures.

Evaluating the genetic contribution of individuals to population structure is essential to select informative individuals for genome sequencing, genotype imputation and to ascertain complex population structures. Existing methods for the selection of informative individuals for genomic imputation solely focus on the identification of key ancestors, which can lead to a loss of phasing accuracy of the reference population. Currently many methods are independently applied to investigate complex population structures. Based on the Eigenvalue Decomposition (EVD) of a genomic relationship matrix we describe a novel approach to evaluate the genetic contribution of individuals to population structure. We combined the identification of key contributors with model-based clustering and population network visualization into an integrated three-step approach, which allows identification of high-resolution population structures and substructures around such key contributors. The approach was applied and validated in four disparate datasets including a simulated population (5,100 individuals and 10,000 SNPs), a highly structured experimental sheep population (1,421 individuals and 44,693 SNPs) and two large complex pedigree populations namely horse (1,077 individuals and 38,124 SNPs) and cattle (2,457 individuals and 45,765 SNPs). In the simulated and experimental sheep dataset, our method, which is unsupervised, successfully identified all known key contributors. Applying our three-step approach to the horse and cattle populations, we observed high-resolution population substructures including the absence of obvious important key contributors. Furthermore, we show that compared to commonly applied strategies to select informative individuals for genotype imputation including the computation of marginal gene contributions (Pedig) and the optimization of genetic relatedness (Rel), the selection of key contributors provided the highest phasing accuracies within the selected reference populations. The presented approach opens new perspectives in the characterization and informed management of populations in general, and in areas such as conservation genetics and selective animal breeding in particular, where assessing the genetic contribution of influential and admixed individuals is crucial for research and management applications. As such, this method provides a valuable complement to common applied tools to visualize complex population structures and to select individuals for re-sequencing.

netview p: a network visualization tool to unravel complex population structure using genome-wide SNPs.

Network-based approaches are emerging as valuable tools for the analysis of complex genetic structure in wild and captive populations. netview p combines data quality control with the construction of population networks through mutual k-nearest neighbours thresholds applied to genome-wide SNPs. The program is cross-platform compatible, open-source and efficiently operates on data ranging from hundreds to hundreds of thousands of SNPs. The pipeline was used for the analysis of pedigree data from simulated (n = 750, SNPs = 1279) and captive silver-lipped pearl oysters (n = 415, SNPs = 1107), wild populations of the European hake from the Atlantic and Mediterranean (n = 834, SNPs = 380) and grey wolves from North America (n = 239, SNPs = 78 255). The population networks effectively visualize large- and fine-scale genetic structure within and between populations, including family-level structure and relationships. netview p comprises a network-based addition to other population analysis tools and provides user-friendly access to a complex network analysis pipeline through implementation in python.

Single-molecule sequencing reveals the molecular basis of multidrug-resistance in ST772 methicillin-resistant Staphylococcus aureus.

BACKGROUND: Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of hospital-associated infection, but there is growing awareness of the emergence of multidrug-resistant lineages in community settings around the world. One such lineage is ST772-MRSA-V, which has disseminated globally and is increasingly prevalent in India. Here, we present the complete genome sequence of DAR4145, a strain of the ST772-MRSA-V lineage from India, and investigate its genomic characteristics in regards to antibiotic resistance and virulence factors. RESULTS: Sequencing using single-molecule real-time technology resulted in the assembly of a single continuous chromosomal sequence, which was error-corrected, annotated and compared to nine draft genome assemblies of ST772-MRSA-V from Australia, Malaysia and India. We discovered numerous and redundant resistance genes associated with mobile genetic elements (MGEs) and known core genome mutations that explain the highly antibiotic resistant phenotype of DAR4145. Staphylococcal toxins and superantigens, including the leukotoxin Panton-Valentinin Leukocidin, were predominantly associated with genomic islands and the phage φ-IND772PVL. Some of these mobile resistance and virulence factors were variably present in other strains of the ST772-MRSA-V lineage. CONCLUSIONS: The genomic characteristics presented here emphasize the contribution of MGEs to the emergence of multidrug-resistant and highly virulent strains of community-associated MRSA. Antibiotic resistance was further augmented by chromosomal mutations and redundancy of resistance genes. The complete genome of DAR4145 provides a valuable resource for future investigations into the global dissemination and phylogeography of ST772-MRSA-V.

Effects of growth phase and nitrogen starvation on expression of fatty acid desaturases and fatty acid composition of Isochrysis aff. galbana (TISO).

Very long-chain polyunsaturated fatty acids (VLC-PUFAs) are important dietary requirements for maintaining human health. Many marine microalgae are naturally high in ω-3 VLC-PUFAs, however, the molecular mechanisms underpinning fatty acid (FA) desaturation and elongation in algae are poorly understood. An advanced molecular understanding would facilitate improvements of this nascent industry. We aimed to investigate expression responses of four front-end fatty acid desaturase genes and downstream effects on FA profiles to nitrogen limitation and cultivation growth stage in Isochrysis aff. galbana (TISO). Cultures were grown in nitrogen-replete and -deplete medium; samples were harvested during logarithmic, late logarithmic and stationary growth phases to analyse FA content/composition and gene expression of ∆(6)-, ∆(8)-, ∆(5)- and ∆(4)-desaturases (d6FAD (putative), d8FAD, d5FAD and d4FAD, respectively). d6FAD (putative) exhibited no differential expression, while d8FAD, d5FAD and d4FAD were significantly upregulated during logarithmic growth of nutrient-replete cultures, coinciding with rapid cell division. In conclusion, it is demonstrated that expression of some FADs in I. aff. galbana varies with culture age and nitrogen status which has downstream consequences on FA desaturation levels. This has implications for the commercial production of VLC-PUFAs where a trade-off between total lipid yield and VLC-PUFAs has to be made.

The effect of nitrogen limitation on acetyl-CoA carboxylase expression and fatty acid content in Chromera velia and Isochrysis aff. galbana (TISO).

Lipids from microalgae have become a valuable product with applications ranging from biofuels to human nutrition. While changes in fatty acid (FA) content and composition under nitrogen limitation are well documented, the involved molecular mechanisms are poorly understood. Acetyl-CoA carboxylase (ACCase) is a key enzyme in the FA synthesis and elongation pathway. Plastidial and cytosolic ACCases provide malonyl-CoA for de novo FA synthesis in the plastid and FA elongation in the endoplasmic reticulum, respectively. The present study aimed at investigating the expression of plastidial and cytosolic ACCase in Chromera velia and Isochrysis aff. galbana (TISO) and their impact on FA content and elongation level when grown under nitrogen-deplete conditions. In C. velia, plastidial ACCase was significantly upregulated during nitrogen starvation and with culture age, strongly correlating with increased FA content. Conversely, plastidial ACCase of I. aff. galbana was not differentially expressed in nitrogen-deplete cultures, but upregulated during the logarithmic phase of nitrogen-replete cultures. In contrast to plastidial ACCase, the cytosolic ACCase of C. velia was downregulated with culture age and nitrogen-starvation, strongly correlating with an increase in medium-chain FAs. In conclusion, the expression of plastidial and cytosolic ACCase changed with growth phase and nutrient status in a species-specific manner and nitrogen limitation did not always result in FA accumulation.