Pre-vaccination frequency of circulatory Tfh is associated with robust immune response to TV003 dengue vaccine.

It has been estimated that more than 390 million people are infected with Dengue virus every year; around 96 millions of these infections result in clinical pathologies. To date, there is only one licensed viral vector-based Dengue virus vaccine CYD-TDV approved for use in dengue endemic areas. While initially approved for administration independent of serostatus, the current guidance only recommends the use of this vaccine for seropositive individuals. Therefore, there is a critical need for investigating the influence of Dengue virus serostatus and immunological mechanisms that influence vaccine outcome. Here, we provide comprehensive evaluation of sero-status and host immune factors that correlate with robust immune responses to a Dengue virus vector based tetravalent vaccine (TV003) in a Phase II clinical cohort of human participants. We observed that sero-positive individuals demonstrate a much stronger immune response to the TV003 vaccine. Our multi-layered immune profiling revealed that sero-positive subjects have increased baseline/pre-vaccination frequencies of circulating T follicular helper (cTfh) cells and the Tfh related chemokine CXCL13/BLC. Importantly, this baseline/pre-vaccination cTfh profile correlated with the vaccinees' ability to launch neutralizing antibody response against all four sero-types of Dengue virus, an important endpoint for Dengue vaccine clinical trials. Overall, we provide novel insights into the favorable cTfh related immune status that persists in Dengue virus sero-positive individuals that correlate with their ability to mount robust vaccine specific immune responses. Such detailed interrogation of cTfh cell biology in the context of clinical vaccinology will help uncover mechanisms and targets for favorable immuno-modulatory agents.

Long-term Sudan virus Ebola survivors maintain multiple antiviral defense mechanisms.

BACKGROUND: The critical issues of sustained memory immunity following ebolavirus disease among long-term survivors (EVD) are still unclear. METHODS: Here, we examine virus-specific immune and inflammatory responses in 12 Sudan virus (SUDV) long-term survivors from Uganda's 2000-1 Gulu outbreak, 15 years after recovery following in vitro challenge. Total RNA from isolated SUDV-stimulated and unstimulated PBMCs was extracted and analyzed. Matched serum samples were also collected to determine SUDV IgG levels and functionality. RESULTS: We detected persistent humoral (58%, 7 of 12) and cellular (33%, 4 of 12) immune responses in SUDV long-term survivors and identified critical molecular mechanisms of innate and adaptive immunity. Gene expression in immune pathways, the IFN signaling system, antiviral defense response, and activation and regulation of T- and B-cell responses were observed. SUDV long-term survivors also maintained robust virus-specific IgG antibodies capable of polyfunctional responses, including neutralizing and innate Fc effector functions. CONCLUSIONS: Data integration identified significant correlations among humoral and cellular immune responses and pinpointed a specific innate and adaptive gene expression signature associated with long-lasting immunity. This could help identify natural and vaccine correlates of protection against ebolavirus disease.

Phylogenomics of nontuberculous mycobacteria respiratory infections in people with cystic fibrosis.

Nontuberculous mycobacteria (NTM) can cause severe pulmonary disease in people with cystic fibrosis (pwCF). These infections present unique challenges for diagnosis and treatment, prompting a recent interest in understanding NTM transmission and pathogenesis during chronic infection. Major gaps remain in our knowledge regarding basic pathogenesis, immune evasion strategies, population dynamics, recombination potential, and the evolutionary implications of host and antibiotic pressures of long-term NTM infections in pwCF. Phylogenomic techniques have emerged as an important tool for tracking global patterns of transmission and are beginning to be used to ask fundamental biological questions about adaptation to the host during pathogenesis. In this review, we discuss the burden of NTM lung disease (NTM-LD), highlight the use of phylogenomics in NTM research, and address the clinical implications associated with these studies.

Customized optical mapping by CRISPR-Cas9 mediated DNA labeling with multiple sgRNAs.

Whole-genome mapping technologies have been developed as a complementary tool to provide scaffolds for genome assembly and structural variation analysis (1,2). We recently introduced a novel DNA labeling strategy based on a CRISPR-Cas9 genome editing system, which can target any 20bp sequences. The labeling strategy is specifically useful in targeting repetitive sequences, and sequences not accessible to other labeling methods. In this report, we present customized mapping strategies that extend the applications of CRISPR-Cas9 DNA labeling. We first design a CRISPR-Cas9 labeling strategy to interrogate and differentiate the single allele differences in NGG protospacer adjacent motifs (PAM sequence). Combined with sequence motif labeling, we can pinpoint the single-base differences in highly conserved sequences. In the second strategy, we design mapping patterns across a genome by selecting sets of specific single-guide RNAs (sgRNAs) for labeling multiple loci of a genomic region or a whole genome. By developing and optimizing a single tube synthesis of multiple sgRNAs, we demonstrate the utility of CRISPR-Cas9 mapping with 162 sgRNAs targeting the 2Mb Haemophilus influenzae chromosome. These CRISPR-Cas9 mapping approaches could be particularly useful for applications in defining long-distance haplotypes and pinpointing the breakpoints in large structural variants in complex genomes and microbial mixtures.

The bacterial microbiota of Hunner lesion interstitial cystitis/bladder pain syndrome.

OBJECTIVE: To undertake the first comprehensive evaluation of the urinary microbiota associated with Hunner lesion (HL) interstitial cystitis/bladder pain syndrome (IC/BPS). Despite no previous identification of a distinct IC/BPS microbial urotype, HL IC/BPS, an inflammatory subtype of IC/BPS, was hypothesized most likely to be associated with a specific bacterial species or microbial pattern. PARTICIPANTS AND METHODS: The bacterial microbiota of midstream urine specimens from HL IC/BPS and age- and gender-matched IC/BPS patients without HL (non-HL IC/BPS) were examined using the pan-bacterial domain clinical-level molecular diagnostic Pacific Biosciences full-length 16S gene sequencing protocol, informatics pipeline and database. We characterized the differential presence, abundances, and diversity of species, as well as gender-specific differences between and among HL and non-HL IC/BPS patients. RESULTS: A total of 59 patients with IC/BPS were enrolled (29 HL, 30 non-HL; 43 women, 16 men) from a single centre and the microbiota in midstream urine specimens was available for comparison. The species abundance differentiation between the HL and non-HL groups (12 species) was not significantly different after Bonferroni adjustments for multiple comparisons. Similarly, the nine differentiating species noted between female HL and non-HL patients were not significantly different after similar statistical correction. However, four species abundances (out of the 10 species differences identified prior to correction) remained significantly different between male HL and non-HL subjects: Negativicoccus succinivorans, Porphyromonas somerae, Mobiluncus curtisii and Corynebacterium renale. Shannon diversity metrics showed significantly higher diversity among HL male patients than HL female patients (P = 0.045), but no significant diversity differences between HL and non-HL patients overall. CONCLUSIONS: We were not able to identify a unique pathogenic urinary microbiota that differentiates all HL from all non-HL IC/BPS. It is likely that the male-specific differences resulted from colonization/contamination remote from the bladder. We were not able to show that bacteria play an important role in patients with HL IC/BPS.

Learning, visualizing and exploring 16S rRNA structure using an attention-based deep neural network.

Recurrent neural networks with memory and attention mechanisms are widely used in natural language processing because they can capture short and long term sequential information for diverse tasks. We propose an integrated deep learning model for microbial DNA sequence data, which exploits convolutional neural networks, recurrent neural networks, and attention mechanisms to predict taxonomic classifications and sample-associated attributes, such as the relationship between the microbiome and host phenotype, on the read/sequence level. In this paper, we develop this novel deep learning approach and evaluate its application to amplicon sequences. We apply our approach to short DNA reads and full sequences of 16S ribosomal RNA (rRNA) marker genes, which identify the heterogeneity of a microbial community sample. We demonstrate that our implementation of a novel attention-based deep network architecture, Read2Pheno, achieves read-level phenotypic prediction. Training Read2Pheno models will encode sequences (reads) into dense, meaningful representations: learned embedded vectors output from the intermediate layer of the network model, which can provide biological insight when visualized. The attention layer of Read2Pheno models can also automatically identify nucleotide regions in reads/sequences which are particularly informative for classification. As such, this novel approach can avoid pre/post-processing and manual interpretation required with conventional approaches to microbiome sequence classification. We further show, as proof-of-concept, that aggregating read-level information can robustly predict microbial community properties, host phenotype, and taxonomic classification, with performance at least comparable to conventional approaches. An implementation of the attention-based deep learning network is available at https://github.com/EESI/sequence_attention (a python package) and https://github.com/EESI/seq2att (a command line tool).

Lung transcriptional unresponsiveness and loss of early influenza virus control in infected neonates is prevented by intranasal Lactobacillus rhamnosus GG.

Respiratory viral infections contribute substantially to global infant losses and disproportionately affect preterm neonates. Using our previously established neonatal murine model of influenza infection, we demonstrate that three-day old mice are exceptionally sensitive to influenza virus infection and exhibit high mortality and viral load. Intranasal pre- and post-treatment of neonatal mice with Lactobacillus rhamnosus GG (LGG), an immune modulator in respiratory viral infection of adult mice and human preterm neonates, considerably improves neonatal mice survival after influenza virus infection. We determine that both live and heat-killed intranasal LGG are equally efficacious in protection of neonates. Early in influenza infection, neonatal transcriptional responses in the lung are delayed compared to adults. These responses increase by 24 hours post-infection, demonstrating a delay in the kinetics of the neonatal anti-viral response. LGG pretreatment improves immune gene transcriptional responses during early infection and specifically upregulates type I IFN pathways. This is critical for protection, as neonatal mice intranasally pre-treated with IFNß before influenza virus infection are also protected. Using transgenic mice, we demonstrate that the protective effect of LGG is mediated through a MyD88-dependent mechanism, specifically via TLR4. LGG can improve both early control of virus and transcriptional responsiveness and could serve as a simple and safe intervention to protect neonates.

Phenotypic diversity and genotypic flexibility of Burkholderia cenocepacia during long-term chronic infection of cystic fibrosis lungs.

Chronic bacterial infections of the lung are the leading cause of morbidity and mortality in cystic fibrosis patients. Tracking bacterial evolution during chronic infections can provide insights into how host selection pressures-including immune responses and therapeutic interventions-shape bacterial genomes. We carried out genomic and phenotypic analyses of 215 serially collected Burkholderia cenocepacia isolates from 16 cystic fibrosis patients, spanning a period of 2-20 yr and a broad range of epidemic lineages. Systematic phenotypic tests identified longitudinal bacterial series that manifested progressive changes in liquid media growth, motility, biofilm formation, and acute insect virulence, but not in mucoidy. The results suggest that distinct lineages follow distinct evolutionary trajectories during lung infection. Pan-genome analysis identified 10,110 homologous gene clusters present only in a subset of strains, including genes restricted to different molecular types. Our phylogenetic analysis based on 2148 orthologous gene clusters from all isolates is consistent with patient-specific clades. This suggests that initial colonization of patients was likely by individual strains, followed by subsequent diversification. Evidence of clonal lineages shared by some patients was observed, suggesting inter-patient transmission. We observed recurrent gene losses in multiple independent longitudinal series, including complete loss of Chromosome III and deletions on other chromosomes. Recurrently observed loss-of-function mutations were associated with decreases in motility and biofilm formation. Together, our study provides the first comprehensive genome-phenome analyses of B. cenocepacia infection in cystic fibrosis lungs and serves as a valuable resource for understanding the genomic and phenotypic underpinnings of bacterial evolution.

Beyond the pan-genome: current perspectives on the functional and practical outcomes of the distributed genome hypothesis.

The principle of monoclonality with regard to bacterial infections was considered immutable prior to 30 years ago. This view, espoused by Koch for acute infections, has proven inadequate regarding chronic infections as persistence requires multiple forms of heterogeneity among the bacterial population. This understanding of bacterial plurality emerged from a synthesis of what-were-then novel technologies in molecular biology and imaging science. These technologies demonstrated that bacteria have complex life cycles, polymicrobial ecologies, and evolve in situ via the horizontal exchange of genic characters. Thus, there is an ongoing generation of diversity during infection that results in far more highly complex microbial communities than previously envisioned. This perspective is based on the fundamental tenet that the bacteria within an infecting population display genotypic diversity, including gene possession differences, which result from horizontal gene transfer mechanisms including transformation, conjugation, and transduction. This understanding is embodied in the concepts of the supragenome/pan-genome and the distributed genome hypothesis (DGH). These paradigms have fostered multiple researches in diverse areas of bacterial ecology including host-bacterial interactions covering the gamut of symbiotic relationships including mutualism, commensalism, and parasitism. With regard to the human host, within each of these symbiotic relationships all bacterial species possess attributes that contribute to colonization and persistence; those species/strains that are pathogenic also encode traits for invasion and metastases. Herein we provide an update on our understanding of bacterial plurality and discuss potential applications in diagnostics, therapeutics, and vaccinology based on perspectives provided by the DGH with regard to the evolution of pathogenicity.

Genome rearrangements induce biofilm formation in Escherichia coli C - an old model organism with a new application in biofilm research.

BACKGROUND: Escherichia coli C forms more robust biofilms than other laboratory strains. Biofilm formation and cell aggregation under a high shear force depend on temperature and salt concentrations. It is the last of five E. coli strains (C, K12, B, W, Crooks) designated as safe for laboratory purposes whose genome has not been sequenced. RESULTS: Here we present the complete genomic sequence of this strain in which we utilized both long-read PacBio-based sequencing and high resolution optical mapping to confirm a large inversion in comparison to the other laboratory strains. Notably, DNA sequence comparison revealed the absence of several genes thought to be involved in biofilm formation, including antigen 43, waaSBOJYZUL for lipopolysaccharide (LPS) synthesis, and cpsB for curli synthesis. The first main difference we identified that likely affects biofilm formation is the presence of an IS3-like insertion sequence in front of the carbon storage regulator csrA gene. This insertion is located 86 bp upstream of the csrA start codon inside the - 35 region of P4 promoter and blocks the transcription from the sigma32 and sigma70 promoters P1-P3 located further upstream. The second is the presence of an IS5/IS1182 in front of the csgD gene. And finally, E. coli C encodes an additional sigma70 subunit driven by the same IS3-like insertion sequence. Promoter analyses using GFP gene fusions provided insights into understanding this regulatory pathway in E. coli. CONCLUSIONS: Biofilms are crucial for bacterial survival, adaptation, and dissemination in natural, industrial, and medical environments. Most laboratory strains of E. coli grown for decades in vitro have evolved and lost their ability to form biofilm, while environmental isolates that can cause infections and diseases are not safe to work with. Here, we show that the historic laboratory strain of E. coli C produces a robust biofilm and can be used as a model organism for multicellular bacterial research. Furthermore, we ascertained the full genomic sequence of this classic strain, which provides for a base level of characterization and makes it useful for many biofilm-based applications.

Transcriptome-Wide Analysis of Hepatitis B Virus-Mediated Changes to Normal Hepatocyte Gene Expression.

Globally, a chronic hepatitis B virus (HBV) infection remains the leading cause of primary liver cancer. The mechanisms leading to the development of HBV-associated liver cancer remain incompletely understood. In part, this is because studies have been limited by the lack of effective model systems that are both readily available and mimic the cellular environment of a normal hepatocyte. Additionally, many studies have focused on single, specific factors or pathways that may be affected by HBV, without addressing cell physiology as a whole. Here, we apply RNA-seq technology to investigate transcriptome-wide, HBV-mediated changes in gene expression to identify single factors and pathways as well as networks of genes and pathways that are affected in the context of HBV replication. Importantly, these studies were conducted in an ex vivo model of cultured primary hepatocytes, allowing for the transcriptomic characterization of this model system and an investigation of early HBV-mediated effects in a biologically relevant context. We analyzed differential gene expression within the context of time-mediated gene-expression changes and show that in the context of HBV replication a number of genes and cellular pathways are altered, including those associated with metabolism, cell cycle regulation, and lipid biosynthesis. Multiple analysis pipelines, as well as qRT-PCR and an independent, replicate RNA-seq analysis, were used to identify and confirm differentially expressed genes. HBV-mediated alterations to the transcriptome that we identified likely represent early changes to hepatocytes following an HBV infection, suggesting potential targets for early therapeutic intervention. Overall, these studies have produced a valuable resource that can be used to expand our understanding of the complex network of host-virus interactions and the impact of HBV-mediated changes to normal hepatocyte physiology on viral replication.

Transformed Recombinant Enrichment Profiling Rapidly Identifies HMW1 as an Intracellular Invasion Locus in Haemophilus influenza.

Many bacterial species actively take up and recombine homologous DNA into their genomes, called natural competence, a trait that offers a means to identify the genetic basis of naturally occurring phenotypic variation. Here, we describe "transformed recombinant enrichment profiling" (TREP), in which natural transformation is used to generate complex pools of recombinants, phenotypic selection is used to enrich for specific recombinants, and deep sequencing is used to survey for the genetic variation responsible. We applied TREP to investigate the genetic architecture of intracellular invasion by the human pathogen Haemophilus influenzae, a trait implicated in persistence during chronic infection. TREP identified the HMW1 adhesin as a crucial factor. Natural transformation of the hmw1 operon from a clinical isolate (86-028NP) into a laboratory isolate that lacks it (Rd KW20) resulted in ~1,000-fold increased invasion into airway epithelial cells. When a distinct recipient (Hi375, already possessing hmw1 and its paralog hmw2) was transformed by the same donor, allelic replacement of hmw2AHi375 by hmw1A86-028NP resulted in a ~100-fold increased intracellular invasion rate. The specific role of hmw1A86-028NP was confirmed by mutant and western blot analyses. Bacterial self-aggregation and adherence to airway cells were also increased in recombinants, suggesting that the high invasiveness induced by hmw1A86-028NP might be a consequence of these phenotypes. However, immunofluorescence results found that intracellular hmw1A86-028NP bacteria likely invaded as groups, instead of as individual bacterial cells, indicating an emergent invasion-specific consequence of hmw1A-mediated self-aggregation.

Natural competence and the evolution of DNA uptake specificity.

Many bacteria are naturally competent, able to actively transport environmental DNA fragments across their cell envelope and into their cytoplasm. Because incoming DNA fragments can recombine with and replace homologous segments of the chromosome, competence provides cells with a potent mechanism of horizontal gene transfer as well as access to the nutrients in extracellular DNA. This review starts with an introductory overview of competence and continues with a detailed consideration of the DNA uptake specificity of competent proteobacteria in the Pasteurellaceae and Neisseriaceae. Species in these distantly related families exhibit strong preferences for genomic DNA from close relatives, a self-specificity arising from the combined effects of biases in the uptake machinery and genomic overrepresentation of the sequences this machinery prefers. Other competent species tested lack obvious uptake bias or uptake sequences, suggesting that strong convergent evolutionary forces have acted on these two families. Recent results show that uptake sequences have multiple "dialects," with clades within each family preferring distinct sequence variants and having corresponding variants enriched in their genomes. Although the genomic consensus uptake sequences are 12 and 29 to 34 bp, uptake assays have found that only central cores of 3 to 4 bp, conserved across dialects, are crucial for uptake. The other bases, which differ between dialects, make weaker individual contributions but have important cooperative interactions. Together, these results make predictions about the mechanism of DNA uptake across the outer membrane, supporting a model for the evolutionary accumulation and stability of uptake sequences and suggesting that uptake biases may be more widespread than currently thought.

The availability of purine nucleotides regulates natural competence by controlling translation of the competence activator Sxy.

Many bacteria are naturally competent, able to bind and take up DNA from their extracellular environment. This DNA can serve as a significant source of nutrients, in addition to providing genetic material for recombination. The regulation of competence in several model organisms highlights the importance of this nutritional function, although it has often been overlooked. Natural competence is induced by starvation in Haemophilus influenzae, the model for competence regulation in the gamma-proteobacteria. This induction depends on the activation of the global metabolic regulator CRP, which occurs upon depletion of phosphotransferase sugars. In this work, we show that the depletion of purine nucleotides under competence-inducing conditions activates the CRP-dependent competence-specific regulator Sxy. Depletion of extra- or intra-cellular purine nucleotides activates Sxy translation, while high levels inhibit it. This is modulated by the stem structure formed by sxy mRNA. The exact mechanism by which the nucleotide depletion signal is transduced is unclear, but it does not involve direct binding of purine intermediates to the sxy stem, and does not require Hfq or competence proteins. Similar regulation occurs in the relatives of H. influenzae, Actinobacillus pneumoniae and A. suis, confirming the importance of processes enabling competent bacteria to exploit the abundant DNA in their environments.

Defining the DNA uptake specificity of naturally competent Haemophilus influenzae cells.

Some naturally competent bacteria exhibit both a strong preference for DNA fragments containing specific 'uptake sequences' and dramatic overrepresentation of these sequences in their genomes. Uptake sequences are often assumed to directly reflect the specificity of the DNA uptake machinery, but the actual specificity has not been well characterized for any bacterium. We produced a detailed analysis of Haemophilus influenzae's uptake specificity, using Illumina sequencing of degenerate uptake sequences in fragments recovered from competent cells. This identified an uptake motif with the same consensus as the motif overrepresented in the genome, with a 9 bp core (AAGTGCGGT) and two short flanking T-rich tracts. Only four core bases (GCGG) were critical for uptake, suggesting that these make strong specific contacts with the uptake machinery. Other core bases had weaker roles when considered individually, as did the T-tracts, but interaction effects between these were also determinants of uptake. The properties of genomic uptake sequences are also constrained by mutational biases and selective forces acting on USSs with coding and termination functions. Our findings define constraints on gene transfer by natural transformation and suggest how the DNA uptake machinery overcomes the physical constraints imposed by stiff highly charged DNA molecules.

A long-read sequencing strategy with overlapping linkers on adjacent fragments (OLAF-Seq) for targeted resequencing and enrichment.

In this report, we present OLAF-Seq, a novel strategy to construct a long-read sequencing library such that adjacent fragments are linked with end-terminal duplications. We use the CRISPR-Cas9 nickase enzyme and a pool of multiple sgRNAs to perform non-random fragmentation of targeted long DNA molecules (> 300kb) into smaller library-sized fragments (about 20 kbp) in a manner so as to retain physical linkage information (up to 1000 bp) between adjacent fragments. DNA molecules targeted for fragmentation are preferentially ligated with adaptors for sequencing, so this method can enrich targeted regions while taking advantage of the long-read sequencing platforms. This enables the sequencing of target regions with significantly lower total coverage, and the genome sequence within linker regions provides information for assembly and phasing. We demonstrated the validity and efficacy of the method first using phage and then by sequencing a panel of 100 full-length cancer-related genes (including both exons and introns) in the human genome. When the designed linkers contained heterozygous genetic variants, long haplotypes could be established. This sequencing strategy can be readily applied in both PacBio and Oxford Nanopore platforms for both long and short genes with an easy protocol. This economically viable approach is useful for targeted enrichment of hundreds of target genomic regions and where long no-gap contigs need deep sequencing.

Characterization of the family-level Borreliaceae pan-genome and development of an episomal typing protocol.

Background The Borreliaceae family includes many obligate parasitic bacterial species which are etiologically associated with a myriad of zoonotic borrelioses including Lyme disease and vector-borne relapsing fevers. Infections by the Borreliaceae are difficult to detect by both direct and indirect methods, often leading to delayed and missed diagnoses. Efforts to improve diagnoses center around the development of molecular diagnostics (MDx), but due to deep tissue sequestration of the causative spirochaetes and the lack of persistent bacteremias, even MDx assays suffer from a lack of sensitivity. Additionally, the highly extensive genomic heterogeneity among isolates, even within the same species, contributes to the lack of assay sensitivity as single target assays cannot provide universal coverage. This within-species heterogeneity is partly due to differences in replicon repertoires and genomic structures that have likely arisen to support the complex Borreliaceae lifecycle in which these parasites have to survive in multiple hosts each with unique immune responses. Results We constructed a Borreliaceae family-level pangenome and characterized the phylogenetic relationships among the constituent taxa which supports the recent taxonomy of splitting the family into at least two genera. Gene content profiles were created for the majority of the Borreliaceae replicons, providing for the first time their unambiguous molecular typing. Conclusion Our characterization of the Borreliaceae pan-genome supports the splitting of the former Borrelia genus into two genera and provides for the phylogenetic placement of several non-species designated isolates. Mining this family-level pangenome will enable precision diagnostics corresponding to gene content-driven clinical outcomes while also providing targets for interventions.

Genome-wide mapping and assembly of structural variant breakpoints in the mouse genome.

Structural variation (SV) is a rich source of genetic diversity in mammals, but due to the challenges associated with mapping SV in complex genomes, basic questions regarding their genomic distribution and mechanistic origins remain unanswered. We have developed an algorithm (HYDRA) to localize SV breakpoints by paired-end mapping, and a general approach for the genome-wide assembly and interpretation of breakpoint sequences. We applied these methods to two inbred mouse strains: C57BL/6J and DBA/2J. We demonstrate that HYDRA accurately maps diverse classes of SV, including those involving repetitive elements such as transposons and segmental duplications; however, our analysis of the C57BL/6J reference strain shows that incomplete reference genome assemblies are a major source of noise. We report 7196 SVs between the two strains, more than two-thirds of which are due to transposon insertions. Of the remainder, 59% are deletions (relative to the reference), 26% are insertions of unlinked DNA, 9% are tandem duplications, and 6% are inversions. To investigate the origins of SV, we characterized 3316 breakpoint sequences at single-nucleotide resolution. We find that approximately 16% of non-transposon SVs have complex breakpoint patterns consistent with template switching during DNA replication or repair, and that this process appears to preferentially generate certain classes of complex variants. Moreover, we find that SVs are significantly enriched in regions of segmental duplication, but that this effect is largely independent of DNA sequence homology and thus cannot be explained by non-allelic homologous recombination (NAHR) alone. This result suggests that the genetic instability of such regions is often the cause rather than the consequence of duplicated genomic architecture.

Transformation of natural genetic variation into Haemophilus influenzae genomes.

Many bacteria are able to efficiently bind and take up double-stranded DNA fragments, and the resulting natural transformation shapes bacterial genomes, transmits antibiotic resistance, and allows escape from immune surveillance. The genomes of many competent pathogens show evidence of extensive historical recombination between lineages, but the actual recombination events have not been well characterized. We used DNA from a clinical isolate of Haemophilus influenzae to transform competent cells of a laboratory strain. To identify which of the ~40,000 polymorphic differences had recombined into the genomes of four transformed clones, their genomes and their donor and recipient parents were deep sequenced to high coverage. Each clone was found to contain ~1000 donor polymorphisms in 3-6 contiguous runs (8.1±4.5 kb in length) that collectively comprised ~1-3% of each transformed chromosome. Seven donor-specific insertions and deletions were also acquired as parts of larger donor segments, but the presence of other structural variation flanking 12 of 32 recombination breakpoints suggested that these often disrupt the progress of recombination events. This is the first genome-wide analysis of chromosomes directly transformed with DNA from a divergent genotype, connecting experimental studies of transformation with the high levels of natural genetic variation found in isolates of the same species.

Genetic Adaptation and Acquisition of Macrolide Resistance in Haemophilus spp. during Persistent Respiratory Tract Colonization in Chronic Obstructive Pulmonary Disease (COPD) Patients Receiving Long-Term Azithromycin Treatment.

Patients with chronic obstructive pulmonary disease (COPD) benefit from the immunomodulatory effect of azithromycin, but long-term administration may alter colonizing bacteria. Our goal was to identify changes in Haemophilus influenzae and Haemophilus parainfluenzae during azithromycin treatment. Fifteen patients were followed while receiving prolonged azithromycin treatment (Hospital Universitari de Bellvitge, Spain). Four patients (P02, P08, P11, and P13) were persistently colonized by H. influenzae for at least 3 months and two (P04 and P11) by H. parainfluenzae. Isolates from these patients (53 H. influenzae and 18 H. parainfluenzae) were included to identify, by whole-genome sequencing, antimicrobial resistance changes and genetic variation accumulated during persistent colonization. All persistent lineages isolated before treatment were azithromycin-susceptible but developed resistance within the first months, apart from those belonging to P02, who discontinued the treatment. H. influenzae isolates from P08-ST107 acquired mutations in 23S rRNA, and those from P11-ST2480 and P13-ST165 had changes in L4 and L22. In H. parainfluenzae, P04 persistent isolates acquired changes in rlmC, and P11 carried genes encoding MefE/MsrD efflux pumps in an integrative conjugative element, which was also identified in H. influenzae P11-ST147. Other genetic variation occurred in genes associated with cell wall and inorganic ion metabolism. Persistent H. influenzae strains all showed changes in licA and hgpB genes. Other genes (lex1, lic3A, hgpC, and fadL) had variation in multiple lineages. Furthermore, persistent strains showed loss, acquisition, or genetic changes in prophage-associated regions. Long-term azithromycin therapy results in macrolide resistance, as well as genetic changes that likely favor bacterial adaptation during persistent respiratory colonization. IMPORTANCE The immunomodulatory properties of azithromycin reduce the frequency of exacerbations and improve the quality of life of COPD patients. However, long-term administration may alter the respiratory microbiota, such as Haemophilus influenzae, an opportunistic respiratory colonizing bacteria that play an important role in exacerbations. This study contributes to a better understanding of COPD progression by characterizing the clinical evolution of H. influenzae in a cohort of patients with prolonged azithromycin treatment. The emergence of macrolide resistance during the first months, combined with the role of Haemophilus parainfluenzae as a reservoir and source of resistance dissemination, is a cause for concern that may lead to therapeutic failure. Furthermore, genetic variations in cell wall and inorganic ion metabolism coding genes likely favor bacterial adaptation to host selective pressures. Therefore, the bacterial pathoadaptive evolution in these severe COPD patients raise our awareness of the possible spread of macrolide resistance and selection of host-adapted clones.