Capturing clinically relevant Campylobacter attributes through direct whole genome sequencing of stool.

Campylobacter is the leading bacterial cause of infectious intestinal disease, but the pathogen typically accounts for a very small proportion of the overall stool microbiome in each patient. Diagnosis is even more difficult due to the fastidious nature of Campylobacter in the laboratory setting. This has, in part, driven a change in recent years, from culture-based to rapid PCR-based diagnostic assays which have improved diagnostic detection, whilst creating a knowledge gap in our clinical and epidemiological understanding of Campylobacter genotypes - no isolates to sequence. In this study, direct metagenomic sequencing approaches were used to assess the possibility of replacing genome sequences with metagenome sequences; metagenomic sequencing outputs were used to describe clinically relevant attributes of Campylobacter genotypes. A total of 37 diarrhoeal stool samples with Campylobacter and five samples with an unknown pathogen result were collected and processed with and without filtration, DNA was extracted, and metagenomes were sequenced by short-read sequencing. Culture-based methods were used to validate Campylobacter metagenome-derived genome (MDG) results. Sequence output metrics were assessed for Campylobacter genome quality and accuracy of characterization. Of the 42 samples passing quality checks for analysis, identification of Campylobacter to the genus and species level was dependent on Campylobacter genome read count, coverage and genome completeness. A total of 65% (24/37) of samples were reliably identified to the genus level through Campylobacter MDG, 73% (27/37) by culture and 97% (36/37) by qPCR. The Campylobacter genomes with a genome completeness of over 60% (n=21) were all accurately identified at the species level (100%). Of those, 72% (15/21) were identified to sequence types (STs), and 95% (20/21) accurately identified antimicrobial resistance (AMR) gene determinants. Filtration of stool samples enhanced Campylobacter MDG recovery and genome quality metrics compared to the corresponding unfiltered samples, which improved the identification of STs and AMR profiles. The phylogenetic analysis in this study demonstrated the clustering of the metagenome-derived with culture-derived genomes and revealed the reliability of genomes from direct stool sequencing. Furthermore, Campylobacter genome spiking percentages ranging from 0 to 2% total metagenome abundance in the ONT MinION sequencer, configured to adaptive sequencing, exhibited better assembly quality and accurate identification of STs, particularly in the analysis of metagenomes containing 2 and 1% of Campylobacter jejuni genomes. Direct sequencing of Campylobacter from stool samples provides clinically relevant and epidemiologically important genomic information without the reliance on cultured genomes.

Codon pairs of the HIV-1 vif gene correlate with CD4+ T cell count.

BACKGROUND: The human APOBEC3G (A3G) protein activity is associated with innate immunity against HIV-1 by inducing high rates of guanosines to adenosines (G-to-A) mutations (viz., hypermutation) in the viral DNA. If hypermutation is not enough to disrupt the reading frames of viral genes, it may likely increase the HIV-1 diversity. To counteract host innate immunity HIV-1 encodes the Vif protein that binds A3G protein and form complexes to be degraded by cellular proteolysis. METHODS: Here we studied the pattern of substitutions in the vif gene and its association with clinical status of HIV-1 infected individuals. To perform the study, unique vif gene sequences were generated from 400 antiretroviral-naïve individuals. RESULTS: The codon pairs: 78-154, 85-154, 101-157, 105-157, and 105-176 of vif gene were associated with CD4+ T cell count lower than 500 cells per mm(3). Some of these codons were located in the (81)LGQGVSIEW(89) region and within the BC-Box. We also identified codons under positive selection clustered in the N-terminal region of Vif protein, between (21)WKSLVK(26) and (40)YRHHY(44) regions (i.e., 31, 33, 37, 39), within the BC-Box (i.e., 155, 159) and the Cullin5-Box (i.e., 168) of vif gene. All these regions are involved in the Vif-induced degradation of A3G/F complexes and the N-terminal of Vif protein binds to viral and cellular RNA. CONCLUSIONS: Adaptive evolution of vif gene was mostly to optimize viral RNA binding and A3G/F recognition. Additionally, since there is not a fully resolved structure of the Vif protein, codon pairs associated with CD4+ T cell count may elucidate key regions that interact with host cell factors. Here we identified and discriminated codons under positive selection and codons under functional constraint in the vif gene of HIV-1.

Exposure of Salmonella biofilms to antibiotic concentrations rapidly selects resistance with collateral tradeoffs.

Most bacteria in nature exist in biofilms, which are inherently tolerant to antibiotics. There is currently very limited understanding of how biofilms evolve in response to sub-lethal concentrations of antimicrobials. In this study, we use a biofilm evolution model to study the effects of sub-inhibitory concentrations of three antibiotics on Salmonella Typhimurium biofilms. We show that biofilms rapidly evolve resistance to each antibiotic they are exposed to, demonstrating a strong selective pressure on biofilms from low antibiotic concentrations. While all antibiotics tested select for clinical resistance, there is no common mechanism. Adaptation to antimicrobials, however, has a marked cost for other clinically important phenotypes, including biofilm formation and virulence. Cefotaxime selects mutants with the greatest deficit in biofilm formation followed by azithromycin and then ciprofloxacin. Understanding the impacts of exposure of biofilms to antibiotics will help understand evolutionary trajectories and may help guide how best to use antibiotics in a biofilm context. Experimental evolution in combination with whole-genome sequencing is a powerful tool for the prediction of evolution trajectories associated with antibiotic resistance in biofilms.

Phylogeography of loaches of the genus lefua (balitoridae, cypriniformes) inferred from mitochondrial DNA sequences.

In order to elucidate phylogenetic relationships and intraspecific variations and to infer the evolutionary process of loaches of the genus Lefua, we analyzed nucleotide sequences of the mitochondrial D-loop region of 100 specimens obtained from 97 localities in Japan and Korea. The genus Lefua includes three described species, L. nikkonis, L. echigonia, and L. costata and an undescribed species, Lefua sp. Our results showed that each species of Lefua formed a monophyletic group, indicating clearly that Lefua species can be genetically distinguished from one another. Lefua nikkonis was the most closely related to L. costata, while L. sp. was the most closely related to L. echigonia. Specimens of L. sp. were grouped into two intraspecific populations and specimens of L. echigonia were grouped into six populations. These populations were well separated geographically from one another by mountain ranges and highlands. We estimated the evolutionary time for splitting of the species and intraspecific populations, and speculated on the evolutionary process of the genus Lefua. Species of Lefua are severely threatened. Fundamental genetic information is indispensable for conservation. We presented genetic background in order to protect these threatened loaches.

Proving universal common ancestry with similar sequences.

Douglas Theobald recently developed an interesting test putatively capable of quantifying the evidence for a Universal Common Ancestry uniting the three domains of life (Eukarya, Archaea and Bacteria) against hypotheses of Independent Origins for some of these domains. We review here his model, in particular in relation to the treatment of Horizontal Gene Transfer (HGT) and to the quality of sequence alignment.

Phylogenetic detection of recombination with a Bayesian prior on the distance between trees.

Genomic regions participating in recombination events may support distinct topologies, and phylogenetic analyses should incorporate this heterogeneity. Existing phylogenetic methods for recombination detection are challenged by the enormous number of possible topologies, even for a moderate number of taxa. If, however, the detection analysis is conducted independently between each putative recombinant sequence and a set of reference parentals, potential recombinations between the recombinants are neglected. In this context, a recombination hotspot can be inferred in phylogenetic analyses if we observe several consecutive breakpoints. We developed a distance measure between unrooted topologies that closely resembles the number of recombinations. By introducing a prior distribution on these recombination distances, a Bayesian hierarchical model was devised to detect phylogenetic inconsistencies occurring due to recombinations. This model relaxes the assumption of known parental sequences, still common in HIV analysis, allowing the entire dataset to be analyzed at once. On simulated datasets with up to 16 taxa, our method correctly detected recombination breakpoints and the number of recombination events for each breakpoint. The procedure is robust to rate and transitionratiotransversion heterogeneities for simulations with and without recombination. This recombination distance is related to recombination hotspots. Applying this procedure to a genomic HIV-1 dataset, we found evidence for hotspots and de novo recombination.