Homopolymeric Tracts Research Articles

Frequently arising ESX-1-associated phase variants influence Mycobacterium tuberculosis fitness in the presence of host and antibiotic pressures.

Mycobacterium tuberculosis (Mtb) is responsible for more deaths worldwide than any other single infectious agent. Understanding how this pathogen adapts to the varied environmental pressures imposed by host immunity and antibiotics has important implications for the design of more effective therapies. In this work, we show that the genome of Mtb contains multiple contingency loci that control the activity of the ESX-1 secretion system, which is critical for interactions with the host. These loci consist of homopolymeric DNA tracts in gene regulatory regions that are subject to high-frequency reversible variation and act to tune the activity of ESX-1. We find that variation at these sites increases the fitness of Mtb in the presence of antibiotic and/or during infection. These findings indicate that Mtb has the ability to diversify its genome in specific sites to create subpopulations of cells that are preadapted to new conditions.

An additional proofreader contributes to DNA replication fidelity in mycobacteria

The removal of mis-incorporated nucleotides by proofreading activity ensures DNA replication fidelity. Whereas the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, it has been shown that proofreading in a majority of bacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase, despite the presence of a DnaQ homolog that is structurally and functionally distinct from E. coli DnaQ. However, the biological functions of this type of noncanonical DnaQ remain unclear. Here, we provide independent evidence that noncanonical DnaQ functions as an additional proofreader for mycobacteria. Using the mutation accumulation assay in combination with whole-genome sequencing, we showed that depletion of DnaQ in Mycolicibacterium smegmatis leads to an increased mutation rate, resulting in AT-biased mutagenesis and increased insertions/deletions in the homopolymer tract. Our results showed that mycobacterial DnaQ binds to the β clamp and functions synergistically with the PHP domain proofreader to correct replication errors. Furthermore, the loss of dnaQ results in replication fork dysfunction, leading to attenuated growth and increased mutagenesis on subinhibitory fluoroquinolones potentially due to increased vulnerability to fork collapse. By analyzing the sequence polymorphism of dnaQ in clinical isolates of Mycobacterium tuberculosis (Mtb), we demonstrated that a naturally evolved DnaQ variant prevalent in Mtb lineage 4.3 may enable hypermutability and is associated with drug resistance. These results establish a coproofreading model and suggest a division of labor between DnaQ and PHP domain proofreader. This study also provides real-world evidence that a mutator-driven evolutionary pathway may exist during the adaptation of Mtb.

Phase variable colony variants are conserved across Gardnerella spp. and exhibit different virulence-associated phenotypes.

The Gardnerella genus, comprising at least 13 species, is associated with the polymicrobial disorder bacterial vaginosis (BV). However, the details of BV pathogenesis are poorly defined, and the contributions made by individual species, including Gardnerella spp., are largely unknown. We report here that colony phenotypes characterized by size (large and small) and opacity (opaque and translucent) are phase variable and are conserved among all tested Gardnerella strains, representing at least 10 different species. With the hypothesis that these different variants could be an important missing piece to the enigma of how BV develops in vivo, we characterized their phenotypic, proteomic, and genomic differences. Beyond increased colony size, large colony variants showed reduced vaginolysin secretion and faster growth rate relative to small colony variants. The ability to inhibit the growth of Neisseria gonorrhoeae and commensal Lactobacillus species varied by strain and, in some instances, differed between variants. Proteomics analyses indicated that 127-173 proteins were differentially expressed between variants. Proteins with increased expression in large variants of both strains were associated with amino acid and protein synthesis and protein folding, whereas those increased in small variants were related to nucleotide synthesis, phosphate transport, ABC transport, and glycogen breakdown. Furthermore, whole genome sequencing analyses revealed an abundance of genes associated with variable homopolymer tracts, implicating slipped strand mispairing in Gardnerella phase variation and illuminating the potential for previously unrecognized heterogeneity within clonal populations. Collectively, these results suggest that phase variants may be primed to serve different roles in BV pathogenesis.IMPORTANCEBacterial vaginosis is the most common gynecological disorder in women of childbearing age. Gardnerella species are crucial to the development of this dysbiosis, but the mechanisms involved in the infection are not understood. We discovered that Gardnerella species vary between two different forms, reflected in bacterial colony size. A slow-growing form makes large amounts of the toxin vaginolysin and is better able to survive in human cervix tissue. A fast-growing form is likely the one that proliferates to high numbers just prior to symptom onset and forms the biofilm that serves as a scaffold for multiple BV-associated anaerobic bacteria. Identification of the proteins that vary between different forms of the bacteria as well as those that vary randomly provides insight into the factors important for Gardnerella infection and immune avoidance.

The first genomic characterization of a stable, hemin-dependent small colony variant strain of Staphylococcus epidermidis isolated from a prosthetic-joint infection

Phenotype switching from a wild type (WT) to a slow-growing subpopulation, referred to as small colony variants (SCVs), supports an infectious lifestyle of Staphylococcus epidermidis, the leading cause of medical device-related infections. Specific mechanisms underlying formation of SCVs and involved in the shaping of their pathogenic potential are of particular interest for stable strains as they have been only rarely cultured from clinical specimens. As the SCV phenotype stability implies the existence of genetic changes, the whole genome sequence of a stable, hemin-dependent S. epidermidis SCV strain (named 49SCV) involved in a late prosthetic joint infection was analyzed. The strain was isolated in a monoculture without a corresponding WT clone, therefore, its genome was compared against five reference S. epidermidis strains (ATCC12228, ATCC14990, NBRC113846, O47, and RP62A), both at the level of the genome structure and coding sequences. According to the Multilocus Sequence Typing analysis, the 49SCV strain represented the sequence type 2 (ST2) regarded as the most prominent infection-causing lineage with a worldwide dissemination. Genomic features unique to 49SCV included the absence of the Staphylococcal Cassette Chromosome (SCC), ~12 kb deletion with the loss of genes involved in the arginine deiminase pathway, and frameshift-generating mutations within the poly(A) and poly(T) homopolymeric tracts. Indels were identified in loci associated with adherence, metabolism, stress response, virulence, and cell wall synthesis. Of note, deletion in the poly(A) of the hemA gene has been considered a possible trigger factor for the phenotype transition and hemin auxotrophy in the strain. To our knowledge, the study represents the first genomic characterization of a clinical, stable and hemin-dependent S. epidermidis SCV strain. We propose that previously unreported indels in the homopolymeric tracts can constitute a background of the SCV phenotype due to a resulting truncation of the corresponding proteins and their possible biological dysfunction. Streamline of genetic content evidenced by the loss of the SCC and a large genomic deletion can represent a possible strategy associated both with the SCV phenotype and its adaptation to chronicity.

Whole genome sequencing in ROHHAD trios proved inconclusive: what's beyond?

Rapid-onset Obesity with Hypothalamic dysfunction, Hypoventilation and Autonomic Dysregulation (ROHHAD) is a rare, life-threatening, pediatric disorder of unknown etiology, whose diagnosis is made difficult by poor knowledge of clinical manifestation, and lack of any confirmatory tests. Children with ROHHAD usually present with rapid onset weight gain which may be followed, over months or years, by hypothalamic dysfunction, hypoventilation, autonomic dysfunction, including impaired bowel motility, and tumors of neural crest origin. Despite the lack of evidence of inheritance in ROHHAD, several studies have been conducted in recent years that have explored possible genetic origins, with unsuccessful results. In order to broaden the search for possible genetic risk factors, an attempt was made to analyse the non-coding variants in two trios (proband with parents), recruited in the Gaslini Children's Hospital in Genoa (Italy). Both patients were females, with a typical history of ROHHAD. Gene variants (single nucleotide variants, short insertions/deletions, splice variants or in tandem expansion of homopolymeric tracts) or altered genomic regions (copy number variations or structural variants) shared between the two probands were searched. Currently, we have not found any potentially pathogenic changes, consistent with the ROHHAD clinical phenotype, and involving genes, regions or pathways shared between the two trios. To definitively rule out the genetic etiology, third-generation sequencing technologies (e.g., long-reads sequencing, optical mapping) should be applied, as well as other pathways, including those associated with immunological and autoimmune disorders, should be explored, making use not only of genomics but also of different -omic datasets.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA).

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

Phase variation as a major mechanism of adaptation in Mycobacterium tuberculosis complex

Phase variation induced by insertions and deletions (INDELs) in genomic homopolymeric tracts (HT) can silence and regulate genes in pathogenic bacteria, but this process is not characterized in MTBC (Mycobacterium tuberculosis complex) adaptation. We leverage 31,428 diverse clinical isolates to identify genomic regions including phase-variants under positive selection. Of 87,651 INDEL events that emerge repeatedly across the phylogeny, 12.4% are phase-variants within HTs (0.02% of the genome by length). We estimated the in-vitro frameshift rate in a neutral HT at 100× the neutral substitution rate at [Formula: see text] frameshifts/HT/year. Using neutral evolution simulations, we identified 4,098 substitutions and 45 phase-variants to be putatively adaptive to MTBC (P < 0.002). We experimentally confirm that a putatively adaptive phase-variant alters the expression of espA, a critical mediator of ESX-1-dependent virulence. Our evidence supports the hypothesis that phase variation in the ESX-1 system of MTBC can act as a toggle between antigenicity and survival in the host.

Variations in the poly-histidine repeat motif of HOXA1 contribute to bicuspid aortic valve in mouse and zebrafish

Bicuspid aortic valve (BAV), the most common cardiovascular malformation occurs in 0.5–1.2% of the population. Although highly heritable, few causal mutations have been identified in BAV patients. Here, we report the targeted sequencing of HOXA1 in a cohort of BAV patients and the identification of rare indel variants in the homopolymeric histidine tract of HOXA1. In vitro analysis shows that disruption of this motif leads to a significant reduction in protein half-life and defective transcriptional activity of HOXA1. In zebrafish, targeting hoxa1a ortholog results in aortic valve defects. In vivo assays indicates that these variants behave as dominant negatives leading abnormal valve development. In mice, deletion of Hoxa1 leads to BAV with a very small, rudimentary non-coronary leaflet. We also show that 17% of homozygous Hoxa1−1His knock-in mice present similar phenotype. Genetic lineage tracing in Hoxa1−/− mutant mice reveals an abnormal reduction of neural crest-derived cells in the valve leaflet, which is caused by a failure of early migration of these cells.

Non-canonical interactions tune condensation and kinetically arrested material states of short RNA homopolymers.

Phase separation related processes have been frequently linked to biological compartmentalization and organization, and associated with cellular responses to stress, development, and disease. As a result, understanding biomolecular condensates and the parameters that modulate material and chemical properties has been an area of interest. Although much initial work focused on the role of proteins in phase separation, increasingly RNA-RNA interactions are being recognized for their potential to tune the condensate properties, influencing material properties, morphology, and dynamic organization. While many of these studies have focused on Watson-Crick base pairing as the primary networking interaction, previous work has also suggested a potential for non-canonical RNA-RNA interactions in long homopolymers to promote phase separation and modulate overall material properties. Here, we examine the phase separation of short RNA homopolymers in the presence of divalent cations and disordered peptides. Our studies reveal that particular short homopolymers form internally kinetically arrested droplets with limited monomer mixing despite continued fusion. We also probe the combined cooperative effects of various ions in modulating phase behavior and investigate the potential for homopolymer tracts in slightly more complex sequences to tune the material and phase separation propensities of different oligomers. Collectively, our results suggest a potential for tracts predominantly interacting through non-canonical interactions to tune the phase behavior and mesoscale properties of simplified condensates, as well as offer insight into the fundamental biophysics of RNA-based condensate systems and short homopolymer tracts in more complex biological systems.

Diversity of glpK Gene and Its Effect on Drug Sensitivity in Mycobacterium bovis.

BackgroundGlycerol kinase (glpK) is essential for the first step of glycerol catabolism in Mycobacterium tuberculosis. However, Mycobacterium bovis has been known to grow poorly in glycerol media because of a base insertion in the glpK gene.MethodsWe analyzed the glpK gene sequences of 60 clinical M. bovis isolates, and determined the minimum inhibitory concentration of 14 drugs by microdilution method to evaluate the effect of frameshift mutations on drug sensitivity. The effect of M. bovis growth rate on its drug sensitivity was investigated using bacteria grown on glycerol or pyruvate.ResultsA total of 44 (73.33%) clinical M. bovis isolates have frameshift mutations in a homopolymeric tract of 7 cytosines in the glpK gene. 15.00% M. bovis isolates showed phenotypic drug resistance. Glycerol metabolism-deficient M. bovis showed reduced susceptibility to 9 out of 14 tested drugs. Mutations in the glpK gene can lead to impaired growth in glycerol-based media, while the minimal inhibitory concentration values of slow-growing M. bovis were higher.ConclusionMutations in the glpK gene can lead to slowed growth and reduced susceptibility to drugs in M. bovis, which may contribute to the emergence of drug-resistant M. bovis and pose a threat to human health owing to the zoonotic capacity of M. bovis.

In-field genetic stock identification of overwintering coho salmon in the Gulf of Alaska: Evaluation of Nanopore sequencing for remote real-time deployment.

Genetic stock identification (GSI) from genotyping‐by‐sequencing of single nucleotide polymorphism (SNP) loci has become the gold standard for stock of origin identification in Pacific salmon. The sequencing platforms currently applied require large batch sizes and multiday processing in specialized facilities to perform genotyping by the thousands. However, recent advances in third‐generation single‐molecule sequencing platforms, such as the Oxford Nanopore minION, provide base calling on portable, pocket‐sized sequencers and promise real‐time, in‐field stock identification of variable batch sizes. Here we evaluate utility and comparability to established GSI platforms of at‐sea stock identification of coho salmon (Oncorhynchus kisutch) using targeted SNP amplicon sequencing on the minION platform during a high‐sea winter expedition to the Gulf of Alaska. As long read sequencers are not optimized for short amplicons, we concatenate amplicons to increase coverage and throughput. Nanopore sequencing at‐sea yielded data sufficient for stock assignment for 50 out of 80 individuals. Nanopore‐based SNP calls agreed with Ion Torrent‐based genotypes in 83.25%, but assignment of individuals to stock of origin only agreed in 61.5% of individuals, highlighting inherent challenges of Nanopore sequencing, such as resolution of homopolymer tracts and indels. However, poor representation of assayed salmon in the queried baseline data set contributed to poor assignment confidence on both platforms. Future improvements will focus on lowering turnaround time and cost, increasing accuracy and throughput, as well as augmentation of the existing baselines. If successfully implemented, Nanopore sequencing will provide an alternative method to the large‐scale laboratory approach by providing mobile small batch genotyping to diverse stakeholders.

KLF3 and PAX6 are candidate driver genes in late-stage, MSI-hypermutated endometrioid endometrial carcinomas.

Endometrioid endometrial carcinomas (EECs) are the most common histological subtype of uterine cancer. Late-stage disease is an adverse prognosticator for EEC. The purpose of this study was to analyze EEC exome mutation data to identify late-stage-specific statistically significantly mutated genes (SMGs), which represent candidate driver genes potentially associated with disease progression. We exome sequenced 15 late-stage (stage III or IV) non-ultramutated EECs and paired non-tumor DNAs; somatic variants were called using Strelka, Shimmer, SomaticSniper and MuTect. Additionally, somatic mutation calls were extracted from The Cancer Genome Atlas (TCGA) data for 66 late-stage and 270 early-stage (stage I or II) non-ultramutated EECs. MutSigCV (v1.4) was used to annotate SMGs in the two late-stage cohorts and to derive p-values for all mutated genes in the early-stage cohort. To test whether late-stage SMGs are statistically significantly mutated in early-stage tumors, q-values for late-stage SMGs were re-calculated from the MutSigCV (v1.4) early-stage p-values, adjusting for the number of late-stage SMGs tested. We identified 14 SMGs in the combined late-stage EEC cohorts. When the 14 late-stage SMGs were examined in the TCGA early-stage data, only Krüppel-like factor 3 (KLF3) and Paired box 6 (PAX6) failed to reach significance as early-stage SMGs, despite the inclusion of enough early-stage cases to ensure adequate statistical power. Within TCGA, nonsynonymous mutations in KLF3 and PAX6 were, respectively, exclusive or nearly exclusive to the microsatellite instability (MSI)-hypermutated molecular subgroup and were dominated by insertions-deletions at homopolymer tracts. In conclusion, our findings are hypothesis-generating and suggest that KLF3 and PAX6, which encode transcription factors, are MSI target genes and late-stage-specific SMGs in EEC.

Comparative Genome Analysis Reveals Accumulation of Single-Nucleotide Repeats in Pathogenic Escherichia Lineages.

Homopolymeric tracts (HPTs) can lead to phase variation and DNA replication slippage, driving adaptation to environmental changes and evolution of genes and genomes. However, there is limited information on HPTs in Escherichia; therefore, we conducted a comprehensive cross-strain search for HPTs in Escherichia genomes. We determined the HPT genomic distribution and identified a pattern of high-frequency HPT localization in pathogenic Escherichia lineages. Notably, HPTs localized near transcriptional regulatory genes. Additionally, excessive repeats accumulated in toxin-coding genes. Moreover, the genomic localization of some HPTs might be derived from exogenous DNA, such as that of bacteriophages. Altogether, our findings may prove useful for understanding the role of HPTs in Escherichia genomes.

Tatajuba: exploring the distribution of homopolymer tracts.

Length variation of homopolymeric tracts, which induces phase variation, is known to regulate gene expression leading to phenotypic variation in a wide range of bacterial species. There is no specialized bioinformatics software which can, at scale, exhaustively explore and describe these features from sequencing data. Identifying these is non-trivial as sequencing and bioinformatics methods are prone to introducing artefacts when presented with homopolymeric tracts due to the decreased base diversity. We present tatajuba, which can automatically identify potential homopolymeric tracts and help predict their putative phenotypic impact, allowing for rapid investigation. We use it to detect all tracts in two separate datasets, one of Campylobacter jejuni and one of three Bordetella species, and to highlight those tracts that are polymorphic across samples. With this we confirm homopolymer tract variation with phenotypic impact found in previous studies and additionally find many more with potential variability. The software is written in C and is available under the open source licence GNU GPLv3.

Phase Variation During Host Colonization and Invasion by Campylobacter jejuni and Other Campylobacter Species.

Phase variation (PV) is a phenomenon common to a variety of bacterial species for niche adaption and survival in challenging environments. Among Campylobacter species, PV depends on the presence of intergenic and intragenic hypermutable G/C homopolymeric tracts. The presence of phase-variable genes is of especial interest for species that cause foodborne or zoonotic infections in humans. PV influences the formation and the structure of the lipooligosaccharide, flagella, and capsule in Campylobacter species. PV of components of these molecules is potentially important during invasion of host tissues, spread within hosts and transmission between hosts. Motility is a critical phenotype that is potentially modulated by PV. Variation in the status of the phase-variable genes has been observed to occur during colonization in chickens and mouse infection models. Interestingly, PV is also involved in bacterial survival of attack by bacteriophages even during chicken colonization. This review aims to explore and discuss observations of PV during model and natural infections by Campylobacter species and how PV may affect strategies for fighting infections by this foodborne pathogen.

Полиморфизм ДНК собак (Canis familiaris) и его применение. IV. мтДНК

The structural organization of the mitochondrial genome of dogs is briefly considered, paying more attention to its variable areas, including the control region and its hypervariable parts HV1 and HV2, as well as a block of tandem decameric repeats and a homopolymer tract of cytosines and thymines. It is noted that the detected polymorphism of mitogenomes forming the clades with haplotypes, with few exceptions, are practically unrelated to dog breeds and their geographical place of residence. In general, the polymorphism of the mitochondrial DNA of dogs does not allow unambiguously identifying a specific individual and in criminalistics, in most cases, it is possible with one or another probability, depending on the variable areas taken into analysis, only to exclude suspected dogs that could lead to their owners from further investigation. Moreover, this requires the creation of population databases of canine mitogenomes, which can be used to calculate the probability of coincidence of haplotypes. Examples of individual cases of investigation of crimes in which polymorphism of DNA of dogs were used, including those that ended with a conviction, are given. Attention is drawn to the need for wider introduction into forensic practice of usage of canine DNA, which can obtained from traces left by dogs mainly in the form of saliva or wool (individual hairs).

Experimental Evolution of Campylobacter jejuni Leads to Loss of Motility, rpoN (σ54) Deletion and Genome Reduction.

Evolution experiments in the laboratory have focused heavily on model organisms, often to the exclusion of clinically relevant pathogens. The foodborne bacterial pathogen Campylobacter jejuni belongs to a genus whose genomes are small compared to those of its closest genomic relative, the free-living genus Sulfurospirillum, suggesting genome reduction during the course of evolution to host association. In an in vitro experiment, C. jejuni serially passaged in rich medium in the laboratory exhibited loss of flagellar motility–an essential function for host colonization. At early time points the motility defect was often reversible, but after 35 days of serial culture, motility was irreversibly lost in most cells in 5 independently evolved populations. Population re-sequencing revealed disruptive mutations to genes in the flagellar transcriptional cascade, rpoN (σ54)—therefore disrupting the expression of the genes σ54 regulates—coupled with deletion of rpoN in all evolved lines. Additional mutations were detected in virulence-related loci. In separate in vivo experiments, we demonstrate that a phase variable (reversible) motility mutant carrying an adenine deletion within a homopolymeric tract resulting in truncation of the flagellar biosynthesis gene fliR was deficient for colonization in a C57BL/6 IL-10–/– mouse disease model. Re-insertion of an adenine residue partially restored motility and ability to colonize mice. Thus, a pathogenic C. jejuni strain was rapidly attenuated by experimental laboratory evolution and demonstrated genomic instability during this evolutionary process. The changes observed suggest C. jejuni is able to evolve in a novel environment through genome reduction as well as transition, transversion, and slip-strand mutations.

Genetic Basis Underlying the Hyperhemolytic Phenotype of Streptococcus agalactiae Strain CNCTC10/84.

Streptococcus agalactiae (group B streptococcus [GBS]) is a major cause of infections in newborns, pregnant women, and immunocompromised patients. GBS strain CNCTC10/84 is a clinical isolate that has high virulence in animal models of infection and has been used extensively to study GBS pathogenesis. Two unusual features of this strain are hyperhemolytic activity and hypo-CAMP factor activity. These two phenotypes are typical of GBS strains that are functionally deficient in the CovR-CovS two-component regulatory system. A previous whole-genome sequencing study found that strain CNCTC10/84 has intact covR and covS regulatory genes. We investigated CovR-CovS regulation in CNCTC10/84 and discovered that a single-nucleotide insertion in a homopolymeric tract in the covR promoter region underlies the strong hemolytic activity and weak CAMP activity of this strain. Using isogenic mutant strains, we demonstrate that this single-nucleotide insertion confers significantly decreased expression of covR and covS and altered expression of CovR-CovS-regulated genes, including that of genes encoding β-hemolysin and CAMP factor. This single-nucleotide insertion also confers significantly increased GBS survival in human whole blood ex vivoIMPORTANCE Group B streptococcus (GBS) is the leading cause of neonatal sepsis, pneumonia, and meningitis. GBS strain CNCTC10/84 is a highly virulent blood isolate that has been used extensively to study GBS pathogenesis for over 20 years. Strain CNCTC10/84 has an unusually strong hemolytic activity, but the genetic basis is unknown. In this study, we discovered that a single-nucleotide insertion in an intergenic homopolymeric tract is responsible for the elevated hemolytic activity of CNCTC10/84.

Campylobacter jejuni and Campylobacter coli autotransporter genes exhibit lineage-associated distribution and decay

BackgroundCampylobacter jejuni and Campylobacter coli are major global causes of bacterial gastroenteritis. Whilst several individual colonisation and virulence factors have been identified, our understanding of their role in the transmission, pathogenesis and ecology of Campylobacter has been hampered by the genotypic and phenotypic diversity within C. jejuni and C. coli. Autotransporter proteins are a family of outer membrane or secreted proteins in Gram-negative bacteria such as Campylobacter, which are associated with virulence functions. In this study we have examined the distribution and predicted functionality of the previously described capC and the newly identified, related capD autotransporter gene families in Campylobacter.ResultsTwo capC-like autotransporter families, designated capC and capD, were identified by homology searches of genomes of the genus Campylobacter. Each family contained four distinct orthologs of CapC and CapD. The distribution of these autotransporter genes was determined in 5829 C. jejuni and 1347 C. coli genomes. Autotransporter genes were found as intact, complete copies and inactive formats due to premature stop codons and frameshift mutations. Presence of inactive and intact autotransporter genes was associated with C. jejuni and C. coli multi-locus sequence types, but for capC, inactivation was independent from the length of homopolymeric tracts in the region upstream of the capC gene. Inactivation of capC or capD genes appears to represent lineage-specific gene decay of autotransporter genes. Intact capC genes were predominantly associated with the C. jejuni ST-45 and C. coli ST-828 generalist lineages. The capD3 gene was only found in the environmental C. coli Clade 3 lineage. These combined data support a scenario of inter-lineage and interspecies exchange of capC and subsets of capD autotransporters.ConclusionsIn this study we have identified two novel, related autotransporter gene families in the genus Campylobacter, which are not uniformly present and exhibit lineage-specific associations and gene decay. The distribution and decay of the capC and capD genes exemplifies the erosion of species barriers between certain lineages of C. jejuni and C. coli, probably arising through co-habitation. This may have implications for the phenotypic variability of these two pathogens and provide opportunity for new, hybrid genotypes to emerge.

Genetic heterogeneity of the Spy1336/R28-Spy1337 virulence axis in Streptococcus pyogenes and effect on gene transcript levels and pathogenesis.

Streptococcus pyogenes is a strict human pathogen responsible for more than 700 million infections annually worldwide. Strains of serotype M28 S. pyogenes are typically among the five more abundant types causing invasive infections and pharyngitis in adults and children. Type M28 strains also have an unusual propensity to cause puerperal sepsis and neonatal disease. We recently discovered that a one-nucleotide indel in an intergenic homopolymeric tract located between genes Spy1336/R28 and Spy1337 altered virulence in a mouse model of infection. In the present study, we analyzed size variation in this homopolymeric tract and determined the extent of heterogeneity in the number of tandemly-repeated 79-amino acid domains in the coding region of Spy1336/R28 in large samples of strains recovered from humans with invasive infections. Both repeat sequence elements are highly polymorphic in natural populations of M28 strains. Variation in the homopolymeric tract results in (i) changes in transcript levels of Spy1336/R28 and Spy1337 in vitro, (ii) differences in virulence in a mouse model of necrotizing myositis, and (iii) global transcriptome changes as shown by RNAseq analysis of isogenic mutant strains. Variation in the number of tandem repeats in the coding sequence of Spy1336/R28 is responsible for size variation of R28 protein in natural populations. Isogenic mutant strains in which genes encoding R28 or transcriptional regulator Spy1337 are inactivated are significantly less virulent in a nonhuman primate model of necrotizing myositis. Our findings provide impetus for additional studies addressing the role of R28 and Spy1337 variation in pathogen-host interactions.