Levent Albayrak

Bartonella melophagi in blood of domestic sheep (Ovis aries) and sheep keds (Melophagus ovinus) from the southwestern US: Cultures, genetic characterization, and ecological connections

Bartonella melophagi sp. nov. was isolated from domestic sheep blood and from sheep keds (Melophagus ovinus) from the southwestern United States. The sequence analyses of the reference strain performed by six molecular markers consistently demonstrated that B. melophagi relates to but differ from other Bartonella species isolated from domestic and wild ruminants. Presence of 183 genes specific for B. melophagi, being absent in genomes of other Bartonella species associated with ruminants also supports the separation of this bacterial species from species of other ruminants. Bartonella DNA was detected in all investigated sheep keds; however, culturing of these bacteria from sheep blood rejects a speculation that B. melophagi is an obligatory endosymbiont. Instead, the results support the hypothesis that the domestic sheep is a natural host reservoir for B. melophagi and the sheep ked its main vector. This bacterium was not isolated from the blood of bighorn sheep and domestic goats belonging to the same subfamily Caprinae. B. melophagi has also been shown to be zoonotic and needs to be investigated further.

Genotyping of Bartonella bacteria and their animal hosts: current status and perspectives

Growing evidence demonstrates that bacterial species diversity is substantial, and many of these species are pathogenic in some contexts or hosts. At the same time, laboratories and museums have collected valuable animal tissue and ectoparasite samples that may contain substantial novel information on bacterial prevalence and diversity. However, the identification of bacterial species is challenging, partly due to the difficulty in culturing many microbes and the reliance on molecular data. Although the genomics revolution will surely add to our knowledge of bacterial systematics, these approaches are not accessible to all researchers and rely predominantly on cultured isolates. Thus, there is a need for comprehensive molecular analyses capable of accurately genotyping bacteria from animal tissues or ectoparasites using common methods that will facilitate large-scale comparisons of species diversity and prevalence. To illustrate the challenges of genotyping bacteria, we focus on the genus Bartonella, vector-borne bacteria common in mammals. We highlight the value and limitations of commonly used techniques for genotyping bartonellae and make recommendations for researchers interested in studying the diversity of these bacteria in various samples. Our recommendations could be applicable to many bacterial taxa (with some modifications) and could lead to a more complete understanding of bacterial species diversity.

Microbiome interaction networks and community structure from laboratory-reared and field-collected Aedes aegypti, Aedes albopictus, and Culex quinquefasciatus mosquito vectors

Microbial interactions are an underappreciated force in shaping insect microbiome communities. Although pairwise patterns of symbiont interactions have been identified, we have a poor understanding regarding the scale and the nature of co-occurrence and co-exclusion interactions within the microbiome. To characterize these patterns in mosquitoes, we sequenced the bacterial microbiome of Aedes aegypti, Ae. albopictus, and Culex quinquefasciatus caught in the field or reared in the laboratory and used these data to generate interaction networks. For collections, we used traps that attracted host-seeking or ovipositing female mosquitoes to determine how physiological state affects the microbiome under field conditions. Interestingly, we saw few differences in species richness or microbiome community structure in mosquitoes caught in either trap. Co-occurrence and co-exclusion analysis identified 116 pairwise interactions substantially increasing the list of bacterial interactions observed in mosquitoes. Networks generated from the microbiome of Ae. aegypti often included highly interconnected hub bacteria. There were several instances where co-occurring bacteria co-excluded a third taxa, suggesting the existence of tripartite relationships. Several associations were observed in multiple species or in field and laboratory-reared mosquitoes indicating these associations are robust and not influenced by environmental or host factors. To demonstrate that microbial interactions can influence colonization of the host, we administered symbionts to Ae. aegypti larvae that either possessed or lacked their resident microbiota. We found that the presence of resident microbiota can inhibit colonization of particular bacterial taxa. Our results highlight that microbial interactions in mosquitoes are complex and influence microbiome composition.

The ability of human nuclear DNA to cause false positive low-abundance heteroplasmy calls varies across the mitochondrial genome

BackgroundLow-abundance mutations in mitochondrial populations (mutations with minor allele frequency ≤ 1%), are associated with cancer, aging, and neurodegenerative disorders. While recent progress in high-throughput sequencing technology has significantly improved the heteroplasmy identification process, the ability of this technology to detect low-abundance mutations can be affected by the presence of similar sequences originating from nuclear DNA (nDNA). To determine to what extent nDNA can cause false positive low-abundance heteroplasmy calls, we have identified mitochondrial locations of all subsequences that are common or similar (one mismatch allowed) between nDNA and mitochondrial DNA (mtDNA).ResultsPerformed analysis revealed up to a 25-fold variation in the lengths of longest common and longest similar (one mismatch allowed) subsequences across the mitochondrial genome. The size of the longest subsequences shared between nDNA and mtDNA in several regions of the mitochondrial genome were found to be as low as 11 bases, which not only allows using these regions to design new, very specific PCR primers, but also supports the hypothesis of the non-random introduction of mtDNA into the human nuclear DNA.ConclusionAnalysis of the mitochondrial locations of the subsequences shared between nDNA and mtDNA suggested that even very short (36 bases) single-end sequencing reads can be used to identify low-abundance variation in 20.4% of the mitochondrial genome. For longer (76 and 150 bases) reads, the proportion of the mitochondrial genome where nDNA presence will not interfere found to be 44.5 and 67.9%, when low-abundance mutations at 100% of locations can be identified using 417 bases long single reads. This observation suggests that the analysis of low-abundance variations in mitochondria population can be extended to a variety of large data collections such as NCBI Sequence Read Archive, European Nucleotide Archive, The Cancer Genome Atlas, and International Cancer Genome Consortium.

CoCo: An application to store High-Throughput Sequencing data in compact text and binary file formats

The storage, manipulation, and especially internet transfer of large amounts of data produced by High-Throughput Sequencing (HTS) instruments present major obstacles utilizing the full potential of this promising technology. The current standard is based on storing all data, which are produced in text (FASTQ and FASTA) and often stored in binary (SRA and BAM) formats. To date, significant effort has been devoted to efficiently compressing these cumbersome sequencing data sets in their existing formats. However, given the substantial improvements in the quality of HTS data, we believe that if one can afford to exclude low quality data and read headers, new much more compressed data formats can be used to reduce size of HTS data files by at least two orders of magnitude. Here we present several examples of file formats specifically designed to store only high quality sequencing reads in space efficient text and binary form. The basic principles used to decrease file size include storage of only one copy of a sequence when reads are present in multiple copies; alphabetical sorting of all reads and storage of only the differences (suffixes) between consecutive reads; and optimization of the number of bits/bytes required to store the information in binary formats. While file size reduction depends on properties of the sequencing data, the size of the resulting files can be as low as 0.1 %-5% of the original FASTQ, SRA, or BAM files. The greatest advantage of the proposed formats however, is based on its time and memory efficiency. The time required to convert reads from FASTQ/FAST A files into the proposed formats is up to 10 times faster than gzip and SRA. The conversion of files in the proposed formats back to FAST A is limited only by the time required to read the file from the hard drive. We present the source code of the C++ object (class) implemented to store, sort, and perform I/O operations with equal length subsequences; and two executable LINUX command line applications (CoCo and CoCo-PIus) able to work with all types of sequencing data including paired-end and flexible size reads. Source code, Linux executables, as well as user manual can be downloaded from http://bgl.utmb.edu/publications/34cocoplus.

Secondary Analysis of the NCI-60 Whole Exome Sequencing Data Indicates Significant Presence of Propionibacterium acnes Genomic Material in Leukemia (RPMI-8226) and Central Nervous System (SF-295, SF-539, and SNB-19) Cell Lines.

The NCI-60 human tumor cell line panel has been used in a broad range of cancer research over the last two decades. A landmark 2013 whole exome sequencing study of this panel added an exceptional new resource for cancer biologists. The complementary analysis of the sequencing data produced by this study suggests the presence of Propionibacterium acnes genomic sequences in almost half of the datasets, with the highest abundance in the leukemia (RPMI-8226) and central nervous system (SF-295, SF-539, and SNB-19) cell lines. While the origin of these contaminating bacterial sequences remains to be determined, observed results suggest that computational control for the presence of microbial genomic material is a necessary step in the analysis of the high throughput sequencing (HTS) data.

Using High Throughput DNA Sequencing to Evaluate the Accuracy of Serial Dilution Based Tests of Microbial Activities in Oil Pipelines

Microbial activities have detrimental effects on industrial infrastructure. If not controlled, microbial presence can result in corrosion, biofilm formation, and product degradation. Serial dilution tests are routinely used for evaluating presence and abundance of microorganisms by diluting samples and culturing microbes in specific media designed to support microorganisms with particular properties, such as sulfate reduction. A high-throughput sequencing approach was used to evaluate changes in microbial composition during four standard serial dilution tests. Analysis of 159 isolates revealed significant differences in the microbial compositions of sequential serial dilution titers and identified several cases where: (a) bacteria known to have a detrimental metabolic function (such as acid production) were lost in the serial dilution medium designed to test for this function; (b) bacteria virtually absent in the original sample became dominant in the serial dilution medium. These observations raise concerns regarding the accuracy and overall usefulness of serial dilution tests.

Mosquito bottlenecks alter viral mutant swarm in a tissue and time-dependent manner with contraction and expansion of variant positions and diversity

Abstract Viral diversity is theorized to play a significant role during virus infections, particularly for arthropod-borne viruses (arboviruses) that must infect both vertebrate and invertebrate hosts. To determine how viral diversity influences mosquito infection and dissemination Culex taeniopus mosquitoes were infected with the Venezuelan equine encephalitis virus endemic strain 68U201. Bodies and legs/wings of the mosquitoes were collected individually and subjected to multi-parallel sequencing. Virus sequence diversity was calculated for each tissue. Greater diversity was seen in mosquitoes with successful dissemination versus those with no dissemination. Diversity across time revealed that bottlenecks influence diversity following dissemination to the legs/wings, but levels of diversity are restored by Day 12 post-dissemination. Specific minority variants were repeatedly identified across the mosquito cohort, some in nearly every tissue and time point, suggesting that certain variants are important in mosquito infection and dissemination. This study demonstrates that the interaction between the mosquito and the virus results in changes in diversity and the mutational spectrum and may be essential for successful transition of the bottlenecks associated with arbovirus infection.

Broom: application for non-redundant storage of high throughput sequencing data

Motivation: The data generation capabilities of high throughput sequencing (HTS) instruments have exponentially increased over the last few years, while the cost of sequencing has dramatically decreased allowing this technology to become widely used in biomedical studies. For small labs and individual researchers, however, storage and transfer of large amounts of HTS data present a significant challenge. The recent trends in increased sequencing quality and genome coverage can be used to reconsider HTS data storage strategies. Results: We present Broom, a stand‐alone application designed to select and store only high‐quality sequencing reads at extremely high compression rates. Written in C++, the application accepts single and paired‐end reads in FASTQ and FASTA formats and decompresses data in FASTA format. Availability and implementation: C++ code available at https://scsb.utmb.edu/labgroups/fofanov/broom.asp. Supplementary information: Supplementary data are available at Bioinformatics online.