Search

Background: The Covid-19 pandemic has been characterized by the emergence of novel SARS-CoV-2 variants, each with distinct properties influencing transmission dynamics, immune escape, and virulence, which, in turn, influence their impact on local populations. Swift analysis of the properties of newly emerged variants is essential in the initial days and weeks to enhance readiness and facilitate the scaling of clinical and public health system responses., Methods: This paper introduces a two-variant metapopulation compartmental model of disease transmission to simulate the dynamics of disease transmission during a period of transition to a newly dominant strain. Leveraging novel S-gene dropout analysis data and genomic sequencing data, combined with confirmed Covid-19 case data, we estimate the epidemiological characteristics of the Omicron variant, which replaced the Delta variant in late 2021 in Philadelphia, PA. We utilized a grid-search method to identify plausible combinations of model parameters, followed by an ensemble adjustment Kalman filter for parameter inference., Results: The model successfully estimated key epidemiological parameters; we estimated the ascertainment rate of 0.22 (95% credible interval 0.15-0.29) and transmission rate of 5.0 (95% CI 2.4-6.6) for the Omicron variant., Conclusions: The study demonstrates the potential for this model-inference framework to provide real-time insights during the emergence of novel variants, aiding in timely public health responses., (© 2024. The Author(s).)

Background: The critical issues of sustained memory immunity following ebolavirus disease among long-term survivors are still unclear., Methods: Here, we examine virus-specific immune and inflammatory responses following in vitro challengd in 12 Sudan virus (SUDV) long-term survivors from Uganda's 2000-2001 Gulu outbreak, 15 years after recovery. Total RNA from isolated SUDV-stimulated and unstimulated peripheral blood mononuclear cells 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 interferon 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., Competing Interests: Potential conflicts of interest. All authors: No reported conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2023. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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., (© 2024. The Author(s).)

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., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Ltd.. All rights reserved.)

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.

In our community-based prospective cohort study in young children, we observed a significant increase in pneumococcal serotype 35B nasopharyngeal (NP) commensal colonization during the 2011-2014 timeframe, but these strains were not associated with disease. Beginning in 2015 and continuing through to the present, the serotype 35B virulence changed, and it became the dominant bacteria isolated and associated with pneumococcal acute otitis-media (AOM) in our cohort. We performed comparative analyses of 250 35B isolates obtained from 140 children collected between 2006 and 2019. Changes in prevalence, clonal-complex composition, and antibiotic resistance were analyzed. Seventy-two (29%) of 35B isolates underwent whole-genome sequencing to investigate genomic changes associated with the shift in virulence that resulted in increased rates of 35B-associated AOM disease. 35B strains that were commensals and AOM disease-causing were mainly associated with sequence type (ST) 558. Antibiotic concentrations of β-lactams and ofloxacin necessary to inhibit growth of 35B strains rose significantly (2006-2019) (p<0.005). However, only isolates from the 35B/ST558 showed significant increases in MIC 50 of penicillin and ofloxacin between the years 2006-2014 and 2015-2019 (p=0.007 and p<0.0001). One hundred thirty-eight SNPs located in 34 different genes were significantly associated with post-2015 strains. SNPs were found in nrdG (metal binding, 10%); metP and metN (ABC transporter, 9%); corA (Mg 2+ transporter, 6%); priA (DNA replication, 5%); and on the enzymic gene ldcB (LD-carboxypeptidase, 3%). Pneumococcal serotype 35B strains was a common NP commensal during 2010-2014. In 2015, a shift in increasing number of AOM cases occurred in young children caused by 35B, that was associated with changes in genetic composition and antibiotic susceptibility., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Fuji, Pichichero, Ehrlich, Mell, Ehrlich and Kaur.)

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&#95;attention (a python package) and https://github.com/EESI/seq2att (a command line tool)., Competing Interests: The authors have declared that no competing interests exist.

Background and Aims: Outbreaks of severe and chronic tick-borne diseases (TBDs) are on the rise. This is through the transmission of infectious disease agents to humans during tick feeding. The transmission rate and extent of microbial exchange, however, vary based on the tick microbiome composition. While select microbes are determined to be members of the normal tick microbiome and others are clearly recognized mammalian and/or avian pathogens, the status of many other members of the tick microbiota with respect to human and alternate host pathogenesis remains unclear. Moreover, the species-level 16S microbiome of prominent TBD vectors, including Ixodes pacificus , have not been extensively studied. To elucidate the I. pacificus microbiome composition, we performed a pan-domain species-specific characterization of the bacterial microbiome on adult I. pacificus ticks collected from two regional parks within Western California. Our methods provide for characterizing nuances within cohort microbiomes and their relationships to geo-locale of origin, surrounding fauna, and prevalences of known and suspected pathogens in relation to current TBD epidemiological zones. Methods: Ninety-two adult I. pacificus bacterial microbiomes were characterized using a high-fidelity, pan-domain, species-specific, full-length 16S rRNA amplification method using circular consensus sequencing performed on the Pacific Biosciences Sequel platform. Data analyses were performed with the MCSMRT data analysis package and database. Results: The species-specific I. pacificus microbiome composition illustrates a complex assortment of microflora, including over 900 eubacterial species with high taxonomic diversity, which was revealed to vary by sex and geo-locale, though the use of full-length 16S gene sequencing. The TBD-associated pathogens, such as Borrelia burgdorferi , Anaplasma phagocytophilum , and Rickettsia monacensis , were identified along with a host of bacteria previously unassociated with ticks. Conclusion: Species-level taxonomic classification of the I. pacificus microbiome revealed that full-length bacterial 16S gene sequencing is required for the granularity to elucidate the microbial diversity within and among ticks based on geo-locale.

A key to longevity assurance is the nutrient-sensing mTOR pathway. Inhibition of mTOR extends lifespan in a variety of organisms. However, the downstream effectors of the mTOR pathway for lifespan regulation are elusive. In a recent report, we described the role of Maf1 as a critical lifespan regulator downstream of the mTOR pathway in fission yeast. Maf1 is the master negative regulator of RNA polymerase III-directed transcription (e.g. tRNAs and 5S rRNAs) and is regulated by mTOR-mediated phosphorylation. We demonstrated that Maf1 is required for lifespan extension under calorie restriction or when mTOR is inhibited. We also showed that Maf1 prevents DNA damage at tRNA genes, which appears to contribute to lifespan maintenance by Maf1. Here we highlight these observations and present additional results to discuss the role of the mTOR-Maf1-Pol III axis in promoting genomic integrity in the face of DNA replication-transcription conflicts in order to maintain normal lifespan.

The genomes of naturally competent Pasteurellaceae and Neisseriaceae have many short uptake sequences (USS), which allow them to distinguish self-DNA from foreign DNA. To fully characterize this preference we developed genome-wide maps of DNA uptake using both a sequence-based computational model and genomic DNA that had been sequenced after uptake by and recovery from competent Haemophilus influenzae cells. When DNA fragments were shorter than the average USS spacing of ∼1,000 bp, sharp peaks of uptake were centered at USS and separated by valleys with 1000-fold lower uptake. Long DNA fragments (1.5-17 kb) gave much less variation, with 90% of positions having uptake within 2-fold of the mean. All detectable uptake biases arose from sequences that fit the USS uptake motif. Simulated competition predicted that, in its respiratory tract environment, H. influenzae will efficiently take up its own DNA even when human DNA is present in 100-fold excess., (© 2020 The Author(s).)

Maf1 is the master repressor of RNA polymerase III responsible for transcription of tRNAs and 5S rRNAs. Maf1 is negatively regulated via phosphorylation by the mTOR pathway, which governs protein synthesis, growth control, and lifespan regulation in response to nutrient availability. Inhibiting the mTOR pathway extends lifespan in various organisms. However, the downstream effectors for the regulation of cell homeostasis that are critical to lifespan extension remain elusive. Here we show that fission yeast Maf1 is required for lifespan extension. Maf1's function in tRNA repression is inhibited by mTOR-dependent phosphorylation, whereas Maf1 is activated via dephosphorylation by protein phosphatase complexes, PP4 and PP2A. Mutational analysis reveals that Maf1 phosphorylation status influences lifespan, which is correlated with elevated tRNA and protein synthesis levels in maf1∆ cells. However, mTOR downregulation, which negates protein synthesis, fails to rescue the short lifespan of maf1∆ cells, suggesting that elevated protein synthesis is not a cause of lifespan shortening in maf1∆ cells. Interestingly, maf1∆ cells accumulate DNA damage represented by formation of Rad52 DNA damage foci and Rad52 recruitment at tRNA genes. Loss of the Rad52 DNA repair protein further exacerbates the shortened lifespan of maf1∆ cells. Strikingly, PP4 deletion alleviates DNA damage and rescues the short lifespan of maf1∆ cells even though tRNA synthesis is increased in this condition, suggesting that elevated DNA damage is the major cause of lifespan shortening in maf1∆ cells. We propose that Maf1-dependent inhibition of tRNA synthesis controls fission yeast lifespan by preventing genomic instability that arises at tRNA genes., (© 2019 The Authors. Aging Cell published by the Anatomical Society and John Wiley & Sons Ltd.)

Analysis of microbiome data involves identifying co-occurring groups of taxa associated with sample features of interest (e.g., disease state). Elucidating such relations is often difficult as microbiome data are compositional, sparse, and have high dimensionality. Also, the configuration of co-occurring taxa may represent overlapping subcommunities that contribute to sample characteristics such as host status. Preserving the configuration of co-occurring microbes rather than detecting specific indicator species is more likely to facilitate biologically meaningful interpretations. Additionally, analyses that use taxonomic relative abundances to predict the abundances of different gene functions aggregate predicted functional profiles across taxa. This precludes straightforward identification of predicted functional components associated with subsets of co-occurring taxa. We provide an approach to explore co-occurring taxa using "topics" generated via a topic model and link these topics to specific sample features (e.g., disease state). Rather than inferring predicted functional content based on overall taxonomic relative abundances, we instead focus on inference of functional content within topics, which we parse by estimating interactions between topics and pathways through a multilevel, fully Bayesian regression model. We apply our methods to three publicly available 16S amplicon sequencing datasets: an inflammatory bowel disease dataset, an oral cancer dataset, and a time-series dataset. Using our topic model approach to uncover latent structure in 16S rRNA amplicon surveys, investigators can (1) capture groups of co-occurring taxa termed topics; (2) uncover within-topic functional potential; (3) link taxa co-occurrence, gene function, and environmental/host features; and (4) explore the way in which sets of co-occurring taxa behave and evolve over time. These methods have been implemented in a freely available R package: https://cran.r-project.org/package=themetagenomics, https://github.com/EESI/themetagenomics., Competing Interests: The authors have declared that no competing interests exist.

The CRISPR/Cas9 system has been proposed as a cure strategy for HIV. However, few published guide RNAs (gRNAs) are predicted to cleave the majority of HIV-1 viral quasispecies (vQS) observed within and among patients. We report the design of a novel pipeline to identify gRNAs that target HIV across a large number of infected individuals. Next generation sequencing (NGS) of LTRs from 269 HIV-1-infected samples in the Drexel CARES Cohort was used to select gRNAs with predicted broad-spectrum activity. In silico, D-LTR-P4-227913 (package of the top 4 gRNAs) accounted for all detectable genetic variation within the vQS of the 269 samples and the Los Alamos National Laboratory HIV database. In silico secondary structure analyses from NGS indicated extensive TAR stem-loop malformations predicted to inactivate proviral transcription, which was confirmed by reduced viral gene expression in TZM-bl or P4R5 cells. Similarly, a high sensitivity in vitro CRISPR/Cas9 cleavage assay showed that the top-ranked gRNA was the most effective at cleaving patient-derived HIV-1 LTRs from five patients. Furthermore, the D-LTR-P4-227913 was predicted to cleave a median of 96.1% of patient-derived sequences from other HIV subtypes. These results demonstrate that the gRNAs possess broad-spectrum cutting activity and could contribute to an HIV cure.

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., Competing Interests: The authors have declared that no competing interests exist.

The clustered regularly interspaced short palindromic repeats (CRISPR)-associated Cas9 system has been used to excise the HIV-1 proviral genome from latently infected cells, potentially offering a cure for HIV-infected patients. Recent studies have shown that most published HIV-1 guide RNAs (gRNAs) do not account for the diverse viral quasispecies within or among patients, which continue to diversify with time even in long-term antiretroviral therapy (ART)-suppressed patients. Given this observation, proviral genomes were deep sequenced from 23 HIV-1-infected patients in the Drexel Medicine CNS AIDS Research and Eradication Study cohort at two different visits. Based on the spectrum of integrated proviral DNA polymorphisms observed, three gRNA design strategies were explored: based on the patient's own HIV-1 sequences (personalized), based on consensus sequences from a large sample of patients [broad-spectrum (BS)], or a combination of both approaches. Using a bioinformatic algorithm, the personalized gRNA design was predicted to cut 46 of 48 patient samples at 90% efficiency, whereas the top 4 BS gRNAs (BS4) were predicted to excise provirus from 44 of 48 patient samples with 90% efficiency. Using a mixed design with the top three BS gRNAs plus one personalized gRNA (BS3 + PS1) resulted in predicted excision of provirus from 45 of 48 patient samples with 90% efficiency. In summary, these studies used an algorithmic design strategy to identify potential BS gRNAs to target a spectrum of HIV-1 long teriminal repeat (LTR) quasispecies for use with a small HIV-1-infected population. This approach should advance CRISPR/Cas9 excision technology taking into account the extensive molecular heterogeneity of HIV-1 that persists in situ after prolonged ART.

Background: Pan-bacterial 16S rRNA microbiome surveys performed with massively parallel DNA sequencing technologies have transformed community microbiological studies. Current 16S profiling methods, however, fail to provide sufficient taxonomic resolution and accuracy to adequately perform species-level associative studies for specific conditions. This is due to the amplification and sequencing of only short 16S rRNA gene regions, typically providing for only family- or genus-level taxonomy. Moreover, sequencing errors often inflate the number of taxa present. Pacific Biosciences' (PacBio's) long-read technology in particular suffers from high error rates per base. Herein, we present a microbiome analysis pipeline that takes advantage of PacBio circular consensus sequencing (CCS) technology to sequence and error correct full-length bacterial 16S rRNA genes, which provides high-fidelity species-level microbiome data., Results: Analysis of a mock community with 20 bacterial species demonstrated 100% specificity and sensitivity with regard to taxonomic classification. Examination of a 250-plus species mock community demonstrated correct species-level classification of > 90% of taxa, and relative abundances were accurately captured. The majority of the remaining taxa were demonstrated to be multiply, incorrectly, or incompletely classified. Using this methodology, we examined the microgeographic variation present among the microbiomes of six sinonasal sites, by both swab and biopsy, from the anterior nasal cavity to the sphenoid sinus from 12 subjects undergoing trans-sphenoidal hypophysectomy. We found greater variation among subjects than among sites within a subject, although significant within-individual differences were also observed. Propiniobacterium acnes (recently renamed Cutibacterium acnes) was the predominant species throughout, but was found at distinct relative abundances by site., Conclusions: Our microbial composition analysis pipeline for single-molecule real-time 16S rRNA gene sequencing (MCSMRT, https://github.com/jpearl01/mcsmrt ) overcomes deficits of standard marker gene-based microbiome analyses by using CCS of entire 16S rRNA genes to provide increased taxonomic and phylogenetic resolution. Extensions of this approach to other marker genes could help refine taxonomic assignments of microbial species and improve reference databases, as well as strengthen the specificity of associations between microbial communities and dysbiotic states.

Tracking bacterial evolution during chronic infection provides insights into how host selection pressures shape bacterial genomes. The human-restricted opportunistic pathogen nontypeable Haemophilus influenzae (NTHi) infects the lower airways of patients suffering chronic obstructive pulmonary disease (COPD) and contributes to disease progression. To identify bacterial genetic variation associated with bacterial adaptation to the COPD lung, we sequenced the genomes of 92 isolates collected from the sputum of 13 COPD patients over 1 to 9 years. Individuals were colonized by distinct clonal types (CTs) over time, but the same CT was often reisolated at a later time or found in different patients. Although genomes from the same CT were nearly identical, intra-CT variation due to mutation and recombination occurred. Recurrent mutations in several genes were likely involved in COPD lung adaptation. Notably, nearly a third of CTs were polymorphic for null alleles of ompP1 (also called fadL ), which encodes a bifunctional membrane protein that both binds the human carcinoembryonic antigen-related cell adhesion molecule 1 (hCEACAM1) receptor and imports long-chain fatty acids (LCFAs). Our computational studies provide plausible three-dimensional models for FadL's interaction with hCEACAM1 and LCFA binding. We show that recurrent fadL mutations are likely a case of antagonistic pleiotropy, since loss of FadL reduces NTHi's ability to infect epithelia but also increases its resistance to bactericidal LCFAs enriched within the COPD lung. Supporting this interpretation, truncated fadL alleles are common in publicly available NTHi genomes isolated from the lower airway tract but rare in others. These results shed light on molecular mechanisms of bacterial pathoadaptation and guide future research toward developing novel COPD therapeutics. IMPORTANCE Nontypeable Haemophilus influenzae is an important pathogen in patients with chronic obstructive pulmonary disease (COPD). To elucidate the bacterial pathways undergoing in vivo evolutionary adaptation, we compared bacterial genomes collected over time from 13 COPD patients and identified recurrent genetic changes arising in independent bacterial lineages colonizing different patients. Besides finding changes in phase-variable genes, we found recurrent loss-of-function mutations in the ompP1 ( fadL ) gene. We show that loss of OmpP1/FadL function reduces this bacterium's ability to infect cells via the hCEACAM1 epithelial receptor but also increases its resistance to bactericidal fatty acids enriched within the COPD lung, suggesting a case of antagonistic pleiotropy that restricts Δ fadL strains' niche. These results show how H. influenzae adapts to host-generated inflammatory mediators in the COPD airways., (Copyright © 2018 Moleres et al.)

Clustered regularly interspaced short palindromic repeats (CRISPR) CRISPR-associated protein 9 (Cas9), including specific guide RNAs (gRNAs), can excise integrated human immunodeficiency virus type 1 (HIV-1) provirus from host chromosomes. To date, anti-HIV-1 gRNAs have been designed to account for off-target activity, however, they seldom account for genetic variation in the HIV-1 genome within and between patients, which will be crucial for therapeutic application of this technology. This analysis tests the ability of published anti-HIV-1 gRNAs to cleave publicly available patient-derived HIV-1 sequences to inform gRNA design and provides basic computational tools to researchers in the field.

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., (© 2017 Lee et al.; Published by Cold Spring Harbor Laboratory Press.)

We report here the complete genomic sequence and methylome of Aggregatibacter actinomycetemcomitans strain IDH781. This rough strain is used extensively as a model organism to characterize localized aggressive periodontitis pathogenesis, the basic biology and oral cavity colonization of A. actinomycetemcomitans, and its interactions with other members of the oral microbiome., (Copyright © 2016 May et al.)

Evolutionary divergence of the mitochondrial genome has given rise to distinct haplogroups. These haplogroups have arisen in specific geographical locations and are responsible for subtle functional changes in the mitochondria that may provide an evolutionary advantage in a given environment. Based on these functional differences, haplogroups could define disease susceptibility in chronic settings. In this study, we undertook a detailed neuropsychological analysis of a cohort of long-term HIV-1-infected individuals in conjunction with sequencing of their mitochondrial genomes. Stepwise regression analysis showed that the best model for predicting both working memory and declarative memory were age and years since diagnosis. In contrast, years since diagnosis and sub-haplogroup were significantly predictive of psychomotor speed. Consistent with this, patients with haplogroup L3e obtained better scores on psychomotor speed and dexterity tasks when compared to the remainder of the cohort, suggesting that this haplogroup provides a protective advantage when faced with the combined stress of HIV-1 infection and long-term antiretroviral therapies. Differential performance on declarative memory tasks was noted for individuals with other sub-L haplogroups, but these differences were not as robust as the association between L3e and psychomotor speed and dexterity tasks. This work provides evidence that mitochondrial haplogroup is related to neuropsychological test performance among patients in chronic disease settings such as HIV-1 infection., Competing Interests: The authors have declared that no competing interests exist.

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.

Originally isolated from a pediatric patient with otitis media, Haemophilus influenzae strain 375 (Hi375) has been extensively studied as a model system for intracellular invasion of airway epithelial cells and other pathogenesis traits. Here, we report its complete genome sequence and methylome., (Copyright © 2014 Mell et al.)

Nontypable Haemophilus influenzae (NTHi) has emerged as an important opportunistic pathogen causing infection in adults suffering obstructive lung diseases. Existing evidence associates chronic infection by NTHi to the progression of the chronic respiratory disease, but specific features of NTHi associated with persistence have not been comprehensively addressed. To provide clues about adaptive strategies adopted by NTHi during persistent infection, we compared sequential persistent isolates with newly acquired isolates in sputa from six patients with chronic obstructive lung disease. Pulse field gel electrophoresis (PFGE) identified three patients with consecutive persistent strains and three with new strains. Phenotypic characterisation included infection of respiratory epithelial cells, bacterial self-aggregation, biofilm formation and resistance to antimicrobial peptides (AMP). Persistent isolates differed from new strains in showing low epithelial adhesion and inability to form biofilms when grown under continuous-flow culture conditions in microfermenters. Self-aggregation clustered the strains by patient, not by persistence. Increasing resistance to AMPs was observed for each series of persistent isolates; this was not associated with lipooligosaccharide decoration with phosphorylcholine or with lipid A acylation. Variation was further analyzed for the series of three persistent isolates recovered from patient 1. These isolates displayed comparable growth rate, natural transformation frequency and murine pulmonary infection. Genome sequencing of these three isolates revealed sequential acquisition of single-nucleotide variants in the AMP permease sapC, the heme acquisition systems hgpB, hgpC, hup and hxuC, the 3-deoxy-D-manno-octulosonic acid kinase kdkA, the long-chain fatty acid transporter ompP1, and the phosphoribosylamine glycine ligase purD. Collectively, we frame a range of pathogenic traits and a repertoire of genetic variants in the context of persistent infection by NTHi.

Naturally competent bacterial species actively take up environmental DNA and can incorporate it into their chromosomes by homologous recombination. This can bring genetic variation from environmental DNA to recipient chromosomes, often in multiple long "donor" segments. Here, we report the results of genome sequencing 96 colonies of a laboratory Haemophilus influenzae strain, which had been experimentally transformed by DNA from a diverged clinical isolate. Donor segments averaged 6.9 kb (spanning several genes) and were clustered into recombination tracts of ~19.5 kb. Individual colonies had replaced from 0.1 to 3.2% of their chromosomes, and ~1/3 of all donor-specific single-nucleotide variants were present in at least one recombinant. We found that nucleotide divergence did not obviously limit the locations of recombination tracts, although there were small but significant reductions in divergence at recombination breakpoints. Although indels occasionally transformed as parts of longer recombination tracts, they were common at breakpoints, suggesting that indels typically block progression of strand exchange. Some colonies had recombination tracts in which variant positions contained mixtures of both donor and recipient alleles. These tracts were clustered around the origin of replication and were interpreted as the result of heteroduplex segregation in the original transformed cell. Finally, a pilot experiment demonstrated the utility of natural transformation for genetically dissecting natural phenotypic variation. We discuss our results in the context of the potential to merge experimental and population genetic approaches, giving a more holistic understanding of bacterial gene transfer.

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.

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.

Background: Genome architecture is profoundly influenced by transposable elements (TEs), and natural selection against their harmful effects is a critical factor limiting their spread. Genome defense by the piRNA silencing pathway also plays a crucial role in limiting TE proliferation. How these two forces jointly determine TE abundance is not well understood. To shed light on the nature of factors that predict TE success, we test three distinct hypotheses in the Drosophila genus. First, we determine whether TE abundance and relaxed genome-wide purifying selection on protein sequences are positively correlated. This serves to test the hypothesis that variation in TE abundance in the Drosophila genus can be explained by the strength of natural selection, relative to drift, acting in parallel against mildly deleterious non-synonymous mutations. Second, we test whether increasing TE abundance is correlated with an increased rate of amino-acid evolution in genes encoding the piRNA machinery, as might be predicted by an evolutionary arms race model. Third, we test whether increasing TE abundance is correlated with greater codon bias in genes of the piRNA machinery. This is predicted if increasing TE abundance selects for increased efficiency in the machinery of genome defense., Results: Surprisingly, we find neither of the first two hypotheses to be true. Specifically, we found that genome-wide levels of purifying selection, measured by the ratio of non-synonymous to synonymous substitution rates (ω), were greater in species with greater TE abundance. In addition, species with greater TE abundance have greater levels of purifying selection in the piRNA machinery. In contrast, it appears that increasing TE abundance has primarily driven adaptation in the piRNA machinery by increasing codon bias., Conclusions: These results indicate that within the Drosophila genus, a historically reduced strength of selection relative to drift is unlikely to explain patterns of increased TE success across species. Other factors, such as ecological exposure, are likely to contribute to variation in TE abundances within species. Furthermore, constraints on the piRNA machinery may temper the evolutionary arms race that would drive increasing rates of evolution at the amino acid level. In the face of these constraints, selection may act primarily by improving the translational efficiency of the machinery of genome defense through efficient codon usage.

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.

Trans-acting factors involved in the early meiotic recombination pathway play a major role in promoting homolog pairing during meiosis in many plants, fungi, and mammals. Here we address whether or not allelic sites have higher levels of interaction when in cis to meiotic recombination events in the budding yeast Saccharomyces cerevisiae. We used Cre/loxP site-specific recombination to genetically measure the magnitude of physical interaction between loxP sites located at allelic positions on homologous chromosomes during meiosis. We observed nonrandom coincidence of Cre-mediated loxP recombination events and meiotic recombination events when the two occurred at linked positions. Further experiments showed that a subset of recombination events destined to become crossover products increased the frequency of nearby Cre-mediated loxP recombination. Our results support a simple physical model of homolog pairing in budding yeast, where recombination at numerous genomic positions generally serves to loosely coalign homologous chromosomes, while crossover-bound recombination intermediates locally stabilize interactions between allelic sites.

We report a novel instance of negative interference during Saccharomyces cerevisiae meiosis, where Cre-mediated recombination between pairs of allelic loxP sites is more frequent than expected. We suggest that endogenous crossover recombination mediates cooperative pairing interactions between all four chromatids of a meiotic bivalent.

A unique aspect of meiosis is the segregation of homologous chromosomes at the meiosis I division. The pairing of homologous chromosomes is a critical aspect of meiotic prophase I that aids proper disjunction at anaphase I. We have used a site-specific recombination assay in Saccharomyces cerevisiae to examine allelic interaction levels during meiosis in a series of mutants defective in recombination, chromatin structure, or intracellular movement. Red1, a component of the chromosome axis, and Mnd1, a chromosome-binding protein that facilitates interhomolog interaction, are critical for achieving high levels of allelic interaction. Homologous recombination factors (Sae2, Rdh54, Rad54, Rad55, Rad51, Sgs1) aid in varying degrees in promoting allelic interactions, while the Srs2 helicase appears to play no appreciable role. Ris1 (a SWI2/SNF2 related protein) and Dot1 (a histone methyltransferase) appear to play minor roles. Surprisingly, factors involved in microtubule-mediated intracellular movement (Tub3, Dhc1, and Mlp2) appear to play no appreciable role in homolog juxtaposition, unlike their counterparts in fission yeast. Taken together, these results support the notion that meiotic recombination plays a major role in the high levels of homolog interaction observed during budding yeast meiosis.

A talented artist developed a progressive aphasia syndrome associated with frontotemporal dementia (FTD). As her disease progressed, language and executive skills declined, but her paintings became freer and more original. She demonstrates that artistic development can occur in the setting of language-dominant types of FTD. The study of artistic development in the setting of FTD suggests that language is not required for, and may even inhibit, certain types of visual creativity.