N-glycosylation of cervicovaginal fluid reflects microbial community, immune activity, and pregnancy status.

Human cervicovaginal fluid (CVF) is a complex, functionally important and glycan rich biological fluid, fundamental in mediating physiological events associated with reproductive health. Using a comprehensive glycomic strategy we reveal an extremely rich and complex N-glycome in CVF of pregnant and non-pregnant women, abundant in paucimannose and high mannose glycans, complex glycans with 2-4 N-Acetyllactosamine (LacNAc) antennae, and Poly-LacNAc glycans decorated with fucosylation and sialylation. N-glycosylation profiles were observed to differ in relation to pregnancy status, microbial composition, immune activation, and pregnancy outcome. Compared to CVF from women experiencing term birth, CVF from women who subsequently experienced preterm birth showed lower sialylation, which correlated to the presence of a diverse microbiome, and higher fucosylation, which correlated positively to pro-inflammatory cytokine concentration. This study is the first step towards better understanding the role of cervicovaginal glycans in reproductive health, their contribution to the mechanism of microbial driven preterm birth, and their potential for preventative therapy.

Functional rewiring of G protein-coupled receptor signaling in human labor.

Current strategies to manage preterm labor center around inhibition of uterine myometrial contractions, yet do not improve neonatal outcomes as they do not address activation of inflammation. Here, we identify that during human labor, activated oxytocin receptor (OTR) reprograms the prostaglandin E2 receptor, EP2, in the pregnant myometrium to suppress relaxatory/Gαs-cAMP signaling and promote pro-labor/inflammatory responses via altered coupling of EP2 from Gαq/11 to Gαi/o. The ability of EP2 to signal via Gαi/o is recapitulated with in vitro OT and only following OTR activation, suggesting direct EP2-OTR crosstalk. Super-resolution imaging with computational modeling reveals OT-dependent reorganization of EP2-OTR complexes to favor conformations for Gαi over Gαs activation. A selective EP2 ligand, PGN9856i, activates the relaxatory/Gαs-cAMP pathway but not the pro-labor/inflammatory responses in term-pregnant myometrium, even following OT. Our study reveals a mechanism, and provides a potential therapeutic solution, whereby EP2-OTR functional associations could be exploited to delay preterm labor.

Microbial-driven preterm labour involves crosstalk between the innate and adaptive immune response.

There has been a surge in studies implicating a role of vaginal microbiota in spontaneous preterm birth (sPTB), but most are associative without mechanistic insight. Here we show a comprehensive approach to understand the causative factors of preterm birth, based on the integration of longitudinal vaginal microbiota and cervicovaginal fluid (CVF) immunophenotype data collected from 133 women at high-risk of sPTB. We show that vaginal depletion of Lactobacillus species and high bacterial diversity leads to increased mannose binding lectin (MBL), IgM, IgG, C3b, C5, IL-8, IL-6 and IL-1ß and to increased risk of sPTB. Cervical shortening, which often precedes preterm birth, is associated with Lactobacillus iners and elevated levels of IgM, C3b, C5, C5a and IL-6. These data demonstrate a role for the complement system in microbial-driven sPTB and provide a scientific rationale for the development of live biotherapeutics and complement therapeutics to prevent sPTB.

Chromosomally normal miscarriage is associated with vaginal dysbiosis and local inflammation.

BACKGROUND: Emerging evidence supports an association between vaginal microbiota composition and risk of miscarriage; however, the underlying mechanisms are poorly understood. We aim to investigate the vaginal microbial composition and the local immune response in chromosomally normal and abnormal miscarriages and compare this to uncomplicated pregnancies delivering at term. METHODS: We used 16S rRNA gene based metataxonomics to interrogate the vaginal microbiota in a cohort of 167 women, 93 miscarriages (54 euploid and 39 aneuploid using molecular cytogenetics) and 74 women who delivered at term and correlate this with the aneuploidy status of the miscarriages. We also measured the concentrations of IL-2, IL-4, IL-6, IL-8, TNF-α, IFN-γ, IL-1ß, IL-18 and IL-10 in cervical vaginal fluid. RESULTS: We show that euploid miscarriage is associated with a significantly higher prevalence of Lactobacillus spp. deplete vaginal microbial communities compared to aneuploid miscarriage (P = 0.01). Integration of matched cervicovaginal fluid immune-profiles showed that Lactobacillus spp. depleted vaginal microbiota associated with pro-inflammatory cytokine levels most strongly in euploid miscarriage compared to viable term pregnancy (IL-1ß; P < 0.001, IL-8; P = 0.01, IL-6; P < 0.001). CONCLUSIONS: Our data suggest the vaginal microbiota plays an important aetiological role in euploid miscarriage and may represent a target to modify risk of pregnancy loss.

Large-scale characterisation of the pregnancy vaginal microbiome and sialidase activity in a low-risk Chinese population.

Vaginal microbiota-host interactions are linked to preterm birth (PTB), which continues to be the primary cause of global childhood mortality. Due to population size, the majority of PTB occurs in Asia, yet there have been few studies of the pregnancy vaginal microbiota in Asian populations. Here, we characterized the vaginal microbiome of 2689 pregnant Chinese women using metataxonomics and in a subset (n = 819), the relationship between vaginal microbiota composition, sialidase activity and leukocyte presence and pregnancy outcomes. Vaginal microbiota were most frequently dominated by Lactobacillus crispatus or L. iners, with the latter associated with vaginal leukocyte presence. Women with high sialidase activity were enriched for bacterial vaginosis-associated genera including Gardnerella, Atopobium and Prevotella. Vaginal microbiota composition, high sialidase activity and/or leukocyte presence was not associated with PTB risk suggesting underlying differences in the vaginal microbiota and/or host immune responses of Chinese women, possibly accounting for low PTB rates in this population.

High Viral Diversity and Mixed Infections in Cerebral Spinal Fluid From Cases of Varicella Zoster Virus Encephalitis.

Background: Varicella zoster virus (VZV) may cause encephalitis, both with and without rash. Here we investigate whether viruses recovered from the central nervous system (CNS; encephalitis or meningitis) differ genetically from those recovered from non-CNS samples. Methods: Enrichment-based deep sequencing of 45 VZV genomes from cerebral spinal fluid (CSF), plasma, bronchoalveolar lavage (BAL), and vesicles was carried out with samples collected from 34 patients with and without VZV infection of the CNS. Results: Viral sequences from multiple sites in the same patient were identical at the consensus level. Virus from vesicle fluid and CSF in cases of meningitis showed low-level diversity. By contrast, plasma, BAL, and encephalitis had higher numbers of variant alleles. Two CSF-encephalitis samples had high genetic diversity, with variant frequency patterns typical of mixed infections with different clades. Conclusions: Low viral genetic diversity in vesicle fluid is compatible with previous observations that VZV skin lesions arise from single or low numbers of virions. A similar result was observed in VZV from cases of VZV meningitis, a generally self-limiting infection. CSF from cases of encephalitis had higher diversity with evidence for mixed clade infections in 2 cases. We hypothesize that reactivation from multiple neurons may contribute to the pathogenesis of VZV encephalitis.

Recombination of Globally Circulating Varicella-Zoster Virus.

UNLABELLED: Varicella-zoster virus (VZV) is a human herpesvirus, which during primary infection typically causes varicella (chicken pox) and establishes lifelong latency in sensory and autonomic ganglia. Later in life, the virus may reactivate to cause herpes zoster (HZ; also known as shingles). To prevent these diseases, a live-attenuated heterogeneous vaccine preparation, vOka, is used routinely in many countries worldwide. Recent studies of another alphaherpes virus, infectious laryngotracheitis virus, demonstrate that live-attenuated vaccine strains can recombine in vivo, creating virulent progeny. These findings raised concerns about using attenuated herpesvirus vaccines under conditions that favor recombination. To investigate whether VZV may undergo recombination, which is a prerequisite for VZV vaccination to create such conditions, we here analyzed 115 complete VZV genomes. Our results demonstrate that recombination occurs frequently for VZV. It thus seems that VZV is fully capable of recombination if given the opportunity, which may have important implications for continued VZV vaccination. Although no interclade vaccine wild-type recombinant strains were found, intraclade recombinants were frequently detected in clade 2, which harbors the vaccine strains, suggesting that the vaccine strains have already been involved in recombination events, either in vivo or in vitro during passages in cell culture. Finally, previous partial and complete genomic studies have described strains that do not cluster phylogenetically to any of the five established clades. The additional VZV strains sequenced here, in combination with those previously published, have enabled us to formally define a novel sixth VZV clade. IMPORTANCE: Although genetic recombination has been demonstrated to frequently occur for other human alphaherpesviruses, herpes simplex viruses 1 and 2, only a few ancient and isolated recent recombination events have hitherto been demonstrated for VZV. In the present study, we demonstrate that VZV also frequently undergoes genetic recombination, including strains belonging to the clade containing the vOKA strain.

Rates of vaccine evolution show strong effects of latency: implications for varicella zoster virus epidemiology.

Varicella-zoster virus (VZV) causes chickenpox and shingles, and is found in human populations worldwide. The lack of temporal signal in the diversity of VZV makes substitution rate estimates unreliable, which is a barrier to understanding the context of its global spread. Here, we estimate rates of evolution by studying live attenuated vaccines, which evolved in 22 vaccinated patients for known periods of time, sometimes, but not always undergoing latency. We show that the attenuated virus evolves rapidly (∼ 10(-6) substitutions/site/day), but that rates decrease dramatically when the virus undergoes latency. These data are best explained by a model in which viral populations evolve for around 13 days before becoming latent, but then undergo no replication during latency. This implies that rates of viral evolution will depend strongly on transmission patterns. Nevertheless, we show that implausibly long latency periods are required to date the most recent common ancestor of extant VZV to an "out-of-Africa" migration with humans, as has been previously suggested.

Whole-genome enrichment and sequencing of Chlamydia trachomatis directly from clinical samples.

BACKGROUND: Chlamydia trachomatis is a pathogen of worldwide importance, causing more than 100 million cases of sexually transmitted infections annually. Whole-genome sequencing is a powerful high resolution tool that can be used to generate accurate data on bacterial population structure, phylogeography and mutations associated with antimicrobial resistance. The objective of this study was to perform whole-genome enrichment and sequencing of C. trachomatis directly from clinical samples. METHODS: C. trachomatis positive samples comprising seven vaginal swabs and three urine samples were sequenced without prior in vitro culture in addition to nine cultured C. trachomatis samples, representing different serovars. A custom capture RNA bait set, that captures all known diversity amongst C. trachomatis genomes, was used in a whole-genome enrichment step during library preparation to enrich for C. trachomatis DNA. All samples were sequenced on the MiSeq platform. RESULTS: Full length C. trachomatis genomes (>95-100% coverage of a reference genome) were successfully generated for eight of ten clinical samples and for all cultured samples. The proportion of reads mapping to C. trachomatis and the mean read depth across each genome were strongly linked to the number of bacterial copies within the original sample. Phylogenetic analysis confirmed the known population structure and the data showed potential for identification of minority variants and mutations associated with antimicrobial resistance. The sensitivity of the method was >10-fold higher than other reported methodologies. CONCLUSIONS: The combination of whole-genome enrichment and deep sequencing has proven to be a non-mutagenic approach, capturing all known variation found within C. trachomatis genomes. The method is a consistent and sensitive tool that enables rapid whole-genome sequencing of C. trachomatis directly from clinical samples and has the potential to be adapted to other pathogens with a similar clonal nature.

Evolution of cocirculating varicella-zoster virus genotypes during a chickenpox outbreak in Guinea-Bissau.

UNLABELLED: Varicella-zoster virus (VZV), a double-stranded DNA alphaherpesvirus, is associated with seasonal outbreaks of varicella in nonimmunized populations. Little is known about whether these outbreaks are associated with a single or multiple viral genotypes and whether new mutations rapidly accumulate during transmission. Here, we take advantage of a well-characterized population cohort in Guinea-Bissau and produce a unique set of 23 full-length genome sequences, collected over 7 months from eight households. Comparative sequence analysis reveals that four distinct genotypes cocirculated among the population, three of which were present during the first week of the outbreak, although no patients were coinfected, which indicates that exposure to infectious virus from multiple sources is common during VZV outbreaks. Transmission of VZV was associated with length polymorphisms in the R1 repeat region and the origin of DNA replication. In two cases, these were associated with the formation of distinct lineages and point to the possible coevolution of these loci, despite the lack of any known functional link in VZV or related herpesviruses. We show that these and all other sequenced clade 5 viruses possess a distinct R1 repeat motif that increases the acidity of an ORF11p protein domain and postulate that this has either arisen or been lost following divergence of the major clades. Thus, sequencing of whole VZV genomes collected during an outbreak has provided novel insights into VZV biology, transmission patterns, and (recent) natural history. IMPORTANCE: VZV is a highly infectious virus and the causative agent of chickenpox and shingles, the latter being particularly associated with the risk of painful complications. Seasonal outbreaks of chickenpox are very common among young children, yet little is known about the dynamics of the virus during person-to-person to transmission or whether multiple distinct viruses seed and/or cocirculate during an outbreak. In this study, we have sequenced chickenpox viruses from an outbreak in Guinea-Bissau that are supported by detailed epidemiological data. Our data show that multiple different virus strains seeded and were maintained throughout the 6-month outbreak period and that viruses transmitted between individuals accumulated new mutations in specific genomic regions. Of particular interest is the potential coevolution of two distinct parts of the genomes and our calculations of the rate of viral mutation, both of which increase our understanding of how VZV evolves over short periods of time in human populations.

Deep sequencing of viral genomes provides insight into the evolution and pathogenesis of varicella zoster virus and its vaccine in humans.

Immunization with the vOka vaccine prevents varicella (chickenpox) in children and susceptible adults. The vOka vaccine strain comprises a mixture of genotypes and, despite attenuation, causes rashes in small numbers of recipients. Like wild-type virus, the vaccine establishes latency in neuronal tissue and can later reactivate to cause Herpes zoster (shingles). Using hybridization-based methodologies, we have purified and sequenced vOka directly from skin lesions. We show that alleles present in the vaccine can be recovered from the lesions and demonstrate the presence of a severe bottleneck between inoculation and lesion formation. Genotypes in any one lesion appear to be descended from one to three vaccine-genotypes with a low frequency of novel mutations. No single vOka haplotype and no novel mutations are consistently present in rashes, indicating that neither new mutations nor recombination with wild type are critical to the evolution of vOka rashes. Instead, alleles arising from attenuation (i.e., not derived from free-living virus) are present at lower frequencies in rash genotypes. We identify 11 loci at which the ancestral allele is selected for in vOka rash formation and show genotypes in rashes that have reactivated from latency cannot be distinguished from rashes occurring immediately after inoculation. We conclude that the vOka vaccine, although heterogeneous, has not evolved to form rashes through positive selection in the mode of a quasispecies, but rather alleles that were essentially neutral during the vaccine production have been selected against in the human subjects, allowing us to identify key loci for rash formation.

Next-generation whole genome sequencing identifies the direction of norovirus transmission in linked patients.

BACKGROUND: Noroviruses are a highly transmissible and major cause of nosocomial gastroenteritis resulting in bed and hospital-ward closures. Where hospital outbreaks are suspected, it is important to determine the routes of spread so that appropriate infection-control procedures can be implemented. To investigate a cluster of norovirus cases occurring in children undergoing bone marrow transplant, we undertook norovirus genome sequencing by next-generation methods. Detailed comparison of sequence data from 2 linked cases enabled us to identify the likely direction of spread. METHODS: Norovirus complementary DNA was amplified by overlapping polymerase chain reaction (PCR) from 13 stool samples from 5 diagnostic real-time PCR-positive patients. The amplicons were sequenced by Roche 454, the genomes assembled by de novo assembly, and the data analyzed phylogenetically. RESULTS: Phylogenetic analysis indicated that patients were infected by viruses similar to 4 distinct GII.4 subtypes and 2 patients were linked by the same virus. Of the 14 sites at which there were differences between the consensus sequences of the 2 linked viral genomes, 9 had minor variants present within one or the other patient. Further analysis confirmed that minor variants at all 9 sites in patient B w ere present as the consensus sequence in patient A. CONCLUSIONS: Phylogenetic analysis excluded a common source of infection in this apparent outbreak. Two of 3 patients on the same ward had closely related viruses, raising the possibility of cross-infection despite protective isolation. Analysis of deep sequencing data enabled us to establish the likely direction of nosocomial transmission.

Tracking viral evolution during a disease outbreak: the rapid and complete selective sweep of a circovirus in the endangered Echo parakeet.

Circoviruses are among the smallest and simplest of all viruses, but they are relatively poorly characterized. Here, we intensively sampled two sympatric parrot populations from Mauritius over a period of 11 years and screened for the circovirus Beak and feather disease virus (BFDV). During the sampling period, a severe outbreak of psittacine beak and feather disease, which is caused by BFDV, occurred in Echo parakeets. Consequently, this data set presents an ideal system for studying the evolution of a pathogen in a natural population and to understand the adaptive changes that cause outbreaks. Unexpectedly, we discovered that the outbreak was most likely caused by changes in functionally important regions of the normally conserved replication-associated protein gene and not the immunogenic capsid. Moreover, these mutations were completely fixed in the Echo parakeet host population very shortly after the outbreak. Several capsid alleles were linked to the replication-associated protein outbreak allele, suggesting that whereas the key changes occurred in the latter, the scope of the outbreak and the selective sweep may have been influenced by positive selection in the capsid. We found evidence for viral transmission between the two host populations though evidence for the invasive species as the source of the outbreak was equivocal. Finally, the high evolutionary rate that we estimated shows how rapidly new variation can arise in BFDV and is consistent with recent results from other small single-stranded DNA viruses.

Patterns of MHC selection in African mole-rats, family Bathyergidae: the effects of sociality and habitat.

African mole-rats are a family of rodents exhibiting an eclectic range of social behaviour and occupying a variety of habitat types. These differences are likely to impact upon the risk of parasite transmission and virulence, with increasing sociality predicted to correspond to an increased risk of transmission. We investigate these factors by analysing the major histocompatibility complex (MHC), a set of genes responsible for encoding highly variable intermediaries of the vertebrate adaptive immune response. To this end we assessed selection at exons 2 and 3 of the MHC class II DQalpha1 gene of four African mole-rat species representing a range of social behaviours. We demonstrate that: (i) the overall pattern of selection at these exons differentiates according to the predicted function of different regions, with the presence of positive selection indicating the likely influence of host-parasite coevolution; and (ii) contrary to the often observed and predicted positive correspondence between sociality and the risk of parasite transmission, two highly social African mole-rat species in fact appear to have comparatively weak positive selection, suggesting diminished host immunity and thus a low overall risk of parasite transmission.