Comparative Analysis of Genome of Ehrlichia sp. HF, a Model Bacterium to Study Fatal Human Ehrlichiosis.

BACKGROUND: The genus Ehrlichia consists of tick-borne obligatory intracellular bacteria that can cause deadly diseases of medical and agricultural importance. Ehrlichia sp. HF, isolated from Ixodes ovatus ticks in Japan [also referred to as I. ovatus Ehrlichia (IOE) agent], causes acute fatal infection in laboratory mice that resembles acute fatal human monocytic ehrlichiosis caused by Ehrlichia chaffeensis. As there is no small laboratory animal model to study fatal human ehrlichiosis, Ehrlichia sp. HF provides a needed disease model. However, the inability to culture Ehrlichia sp. HF and the lack of genomic information have been a barrier to advance this animal model. In addition, Ehrlichia sp. HF has several designations in the literature as it lacks a taxonomically recognized name. RESULTS: We stably cultured Ehrlichia sp. HF in canine histiocytic leukemia DH82 cells from the HF strain-infected mice, and determined its complete genome sequence. Ehrlichia sp. HF has a single double-stranded circular chromosome of 1,148,904 bp, which encodes 866 proteins with a similar metabolic potential as E. chaffeensis. Ehrlichia sp. HF encodes homologs of all virulence factors identified in E. chaffeensis, including 23 paralogs of P28/OMP-1 family outer membrane proteins, type IV secretion system apparatus and effector proteins, two-component systems, ankyrin-repeat proteins, and tandem repeat proteins. Ehrlichia sp. HF is a novel species in the genus Ehrlichia, as demonstrated through whole genome comparisons with six representative Ehrlichia species, subspecies, and strains, using average nucleotide identity, digital DNA-DNA hybridization, and core genome alignment sequence identity. CONCLUSIONS: The genome of Ehrlichia sp. HF encodes all known virulence factors found in E. chaffeensis, substantiating it as a model Ehrlichia species to study fatal human ehrlichiosis. Comparisons between Ehrlichia sp. HF and E. chaffeensis will enable identification of in vivo virulence factors that are related to host specificity, disease severity, and host inflammatory responses. We propose to name Ehrlichia sp. HF as Ehrlichia japonica sp. nov. (type strain HF), to denote the geographic region where this bacterium was initially isolated.

Streptococcus pneumoniae in the heart subvert the host response through biofilm-mediated resident macrophage killing.

For over 130 years, invasive pneumococcal disease has been associated with the presence of extracellular planktonic pneumococci, i.e. diplococci or short chains in affected tissues. Herein, we show that Streptococcus pneumoniae that invade the myocardium instead replicate within cellular vesicles and transition into non-purulent biofilms. Pneumococci within mature cardiac microlesions exhibited salient biofilm features including intrinsic resistance to antibiotic killing and the presence of an extracellular matrix. Dual RNA-seq and subsequent principal component analyses of heart- and blood-isolated pneumococci confirmed the biofilm phenotype in vivo and revealed stark anatomical site-specific differences in virulence gene expression; the latter having major implications on future vaccine antigen selection. Our RNA-seq approach also identified three genomic islands as exclusively expressed in vivo. Deletion of one such island, Region of Diversity 12, resulted in a biofilm-deficient and highly inflammogenic phenotype within the heart; indicating a possible link between the biofilm phenotype and a dampened host-response. We subsequently determined that biofilm pneumococci released greater amounts of the toxin pneumolysin than did planktonic or RD12 deficient pneumococci. This allowed heart-invaded wildtype pneumococci to kill resident cardiac macrophages and subsequently subvert cytokine/chemokine production and neutrophil infiltration into the myocardium. This is the first report for pneumococcal biofilm formation in an invasive disease setting. We show that biofilm pneumococci actively suppress the host response through pneumolysin-mediated immune cell killing. As such, our findings contradict the emerging notion that biofilm pneumococci are passively immunoquiescent.

The Streptococcus agalactiae Stringent Response Enhances Virulence and Persistence in Human Blood.

Streptococcus agalactiae (group B Streptococcus [GBS]) causes serious infections in neonates. We previously reported a transposon sequencing (Tn-seq) system for performing genomewide assessment of gene fitness in GBS. In order to identify molecular mechanisms required for GBS to transition from a mucosal commensal lifestyle to bloodstream invasion, we performed Tn-seq on GBS strain A909 with human whole blood. Our analysis identified 16 genes conditionally essential for GBS survival in blood, of which 75% were members of the capsular polysaccharide (cps) operon. Among the non-cps genes identified as conditionally essential was relA, which encodes an enzyme whose activity is central to the bacterial stringent response-a conserved adaptation to environmental stress. We used blood coincubation studies of targeted knockout strains to confirm the expected growth defects of GBS deficient in capsule or stringent response activation. Unexpectedly, we found that the relA knockout strains demonstrated decreased expression of ß-hemolysin/cytolysin, an important cytotoxin implicated in facilitating GBS invasion. Furthermore, chemical activation of the stringent response with serine hydroxamate increased ß-hemolysin/cytolysin expression. To establish a mechanism by which the stringent response leads to increased cytotoxicity, we performed transcriptome sequencing (RNA-seq) on two GBS strains grown under stringent response or control conditions. This revealed a conserved decrease in the expression of genes in the arginine deiminase pathway during stringent response activation. Through coincubation with supplemental arginine and the arginine antagonist canavanine, we show that arginine availability is a determinant of GBS cytotoxicity and that the pathway between stringent response activation and increased virulence is arginine dependent.

A comparative analysis of library prep approaches for sequencing low input translatome samples.

BACKGROUND: Cell type-specific ribosome-pulldown has become an increasingly popular method for analysis of gene expression. It allows for expression analysis from intact tissues and monitoring of protein synthesis in vivo. However, while its utility has been assessed, technical aspects related to sequencing of these samples, often starting with a smaller amount of RNA, have not been reported. In this study, we evaluated the performance of five library prep protocols for ribosome-associated mRNAs when only 250 pg-4 ng of total RNA are used. RESULTS: We obtained total and RiboTag-IP RNA, in three biological replicates. We compared 5 methods of library preparation for Illumina Next Generation sequencing: NuGEN Ovation RNA-Seq system V2 Kit, TaKaRa SMARTer Stranded Total RNA-Seq Kit, TaKaRa SMART-Seq v4 Ultra Low Input RNA Kit, Illumina TruSeq RNA Library Prep Kit v2 and NEBNext® Ultra™ Directional RNA Library Prep Kit using slightly modified protocols each with 4 ng of total RNA. An additional set of samples was processed using the TruSeq kit with 70 ng, as a 'gold standard' control and the SMART-Seq v4 with 250 pg of total RNA. TruSeq-processed samples had the best metrics overall, with similar results for the 4 ng and 70 ng samples. The results of the SMART-Seq v4 processed samples were similar to TruSeq (Spearman correlation > 0.8) despite using lower amount of input RNA. All RiboTag-IP samples had an increase in the intronic reads compared with the corresponding whole tissue, suggesting that the IP captures some immature mRNAs. The SMARTer-processed samples had a higher representation of ribosomal and non-coding RNAs leading to lower representation of protein coding mRNA. The enrichment or depletion of IP samples compared to corresponding input RNA was similar across all kits except for SMARTer kit. CONCLUSION: RiboTag-seq can be performed successfully with as little as 250 pg of total RNA when using the SMART-Seq v4 kit and 4 ng when using the modified protocols of other library preparation kits. The SMART-Seq v4 and TruSeq kits resulted in the highest quality libraries. RiboTag IP RNA contains some immature transcripts.

Neuraminidase A-Exposed Galactose Promotes Streptococcus pneumoniae Biofilm Formation during Colonization.

Streptococcus pneumoniae is an opportunistic pathogen that colonizes the nasopharynx. Herein we show that carbon availability is distinct between the nasopharynx and bloodstream of adult humans: glucose is absent from the nasopharynx, whereas galactose is abundant. We demonstrate that pneumococcal neuraminidase A (NanA), which cleaves terminal sialic acid residues from host glycoproteins, exposed galactose on the surface of septal epithelial cells, thereby increasing its availability during colonization. We observed that S. pneumoniae mutants deficient in NanA and ß-galactosidase A (BgaA) failed to form biofilms in vivo despite normal biofilm-forming abilities in vitro Subsequently, we observed that glucose, sucrose, and fructose were inhibitory for biofilm formation, whereas galactose, lactose, and low concentrations of sialic acid were permissive. Together these findings suggested that the genes involved in biofilm formation were under some form of carbon catabolite repression (CCR), a regulatory network in which genes involved in the uptake and metabolism of less-preferred sugars are silenced during growth with preferred sugars. Supporting this notion, we observed that a mutant deficient in pyruvate oxidase, which converts pyruvate to acetyl-phosphate under non-CCR-inducing growth conditions, was unable to form biofilms. Subsequent comparative transcriptome sequencing (RNA-seq) analyses of planktonic and biofilm-grown pneumococci showed that metabolic pathways involving the conversion of pyruvate to acetyl-phosphate and subsequently leading to fatty acid biosynthesis were consistently upregulated during diverse biofilm growth conditions. We conclude that carbon availability in the nasopharynx impacts pneumococcal biofilm formation in vivo Additionally, biofilm formation involves metabolic pathways not previously appreciated to play an important role.

The essential genome of Streptococcus agalactiae.

BACKGROUND: Next-generation sequencing of transposon-genome junctions from a saturated bacterial mutant library (Tn-seq) is a powerful tool that permits genome-wide determination of the contribution of genes to fitness of the organism under a wide range of experimental conditions. We report development, testing, and results from a Tn-seq system for use in Streptococcus agalactiae (group B Streptococcus; GBS), an important cause of neonatal sepsis. METHODS: Our method uses a Himar1 mini-transposon that inserts at genomic TA dinucleotide sites, delivered to GBS on a temperature-sensitive plasmid that is subsequently cured from the bacterial population. In order to establish the GBS essential genome, we performed Tn-seq on DNA collected from three independent mutant libraries-with at least 135,000 mutants per library-at serial 24 h time points after outgrowth in rich media. RESULTS: After statistical analysis of transposon insertion density and distribution, we identified 13.5 % of genes as essential and 1.2 % as critical, with high levels of reproducibility. Essential and critical genes are enriched for fundamental cellular housekeeping functions, such as acyl-tRNA biosynthesis, nucleotide metabolism, and glycolysis. We further validated our system by comparing fitness assignments of homologous genes in GBS and a close bacterial relative, Streptococcus pyogenes, which demonstrated 93 % concordance. Finally, we used our fitness assignments to identify signal transduction pathway components predicted to be essential or critical in GBS. CONCLUSIONS: We believe that our baseline fitness assignments will be a valuable tool for GBS researchers and that our system has the potential to reveal key pathogenesis gene networks and potential therapeutic/preventative targets.

Extensive duplication of the Wolbachia DNA in chromosome four of Drosophila ananassae.

BACKGROUND: Lateral gene transfer (LGT) from bacterial Wolbachia endosymbionts has been detected in ~20% of arthropod and nematode genome sequencing projects. Many of these transfers are large and contain a substantial part of the Wolbachia genome. RESULTS: Here, we re-sequenced three D. ananassae genomes from Asia and the Pacific that contain large LGTs from Wolbachia. We find that multiple copies of the Wolbachia genome are transferred to the Drosophila nuclear genome in all three lines. In the D. ananassae line from Indonesia, the copies of Wolbachia DNA in the nuclear genome are nearly identical in size and sequence yielding an even coverage of mapped reads over the Wolbachia genome. In contrast, the D. ananassae lines from Hawaii and India show an uneven coverage of mapped reads over the Wolbachia genome suggesting that different parts of these LGTs are present in different copy numbers. In the Hawaii line, we find that this LGT is underrepresented in third instar larvae indicative of being heterochromatic. Fluorescence in situ hybridization of mitotic chromosomes confirms that the LGT in the Hawaii line is heterochromatic and represents ~20% of the sequence on chromosome 4 (dot chromosome, Muller element F). CONCLUSIONS: This collection of related lines contain large lateral gene transfers composed of multiple Wolbachia genomes that constitute >2% of the D. ananassae genome (~5 Mbp) and partially explain the abnormally large size of chromosome 4 in D. ananassae.

Single molecule sequencing and genome assembly of a clinical specimen of Loa loa, the causative agent of loiasis.

BACKGROUND: More than 20% of the world's population is at risk for infection by filarial nematodes and >180 million people worldwide are already infected. Along with infection comes significant morbidity that has a socioeconomic impact. The eight filarial nematodes that infect humans are Wuchereria bancrofti, Brugia malayi, Brugia timori, Onchocerca volvulus, Loa loa, Mansonella perstans, Mansonella streptocerca, and Mansonella ozzardi, of which three have published draft genome sequences. Since all have humans as the definitive host, standard avenues of research that rely on culturing and genetics have often not been possible. Therefore, genome sequencing provides an important window into understanding the biology of these parasites. The need for large amounts of high quality genomic DNA from homozygous, inbred lines; the availability of only short sequence reads from next-generation sequencing platforms at a reasonable expense; and the lack of random large insert libraries has limited our ability to generate high quality genome sequences for these parasites. However, the Pacific Biosciences single molecule, real-time sequencing platform holds great promise in reducing input amounts and generating sufficiently long sequences that bypass the need for large insert paired libraries. RESULTS: Here, we report on efforts to generate a more complete genome assembly for L. loa using genetically heterogeneous DNA isolated from a single clinical sample and sequenced on the Pacific Biosciences platform. To obtain the best assembly, numerous assemblers and sequencing datasets were analyzed, combined, and compared. Quiver-informed trimming of an assembly of only Pacific Biosciences reads by HGAP2 was selected as the final assembly of 96.4 Mbp in 2,250 contigs. This results in ~9% more of the genome in ~85% fewer contigs from ~80% less starting material at a fraction of the cost of previous Roche 454-based sequencing efforts. CONCLUSIONS: The result is the most complete filarial nematode assembly produced thus far and demonstrates the utility of single molecule sequencing on the Pacific Biosciences platform for genetically heterogeneous metazoan genomes.

Rapid transcriptome sequencing of an invasive pest, the brown marmorated stink bug Halyomorpha halys.

BACKGROUND: Halyomorpha halys (Stål) (Insecta:Hemiptera;Pentatomidae), commonly known as the Brown Marmorated Stink Bug (BMSB), is an invasive pest of the mid-Atlantic region of the United States, causing economically important damage to a wide range of crops. Native to Asia, BMSB was first observed in Allentown, PA, USA, in 1996, and this pest is now well-established throughout the US mid-Atlantic region and beyond. In addition to the serious threat BMSB poses to agriculture, BMSB has become a nuisance to homeowners, invading home gardens and congregating in large numbers in human-made structures, including homes, to overwinter. Despite its significance as an agricultural pest with limited control options, only 100 bp of BMSB sequence data was available in public databases when this project began. RESULTS: Transcriptome sequencing was undertaken to provide a molecular resource to the research community to inform the development of pest control strategies and to provide molecular data for population genetics studies of BMSB. Using normalized, strand-specific libraries, we sequenced pools of all BMSB life stages on the Illumina HiSeq. Trinity was used to assemble 200,000 putative transcripts in >100,000 components. A novel bioinformatic method that analyzed the strand-specificity of the data reduced this to 53,071 putative transcripts from 18,573 components. By integrating multiple other data types, we narrowed this further to 13,211 representative transcripts. CONCLUSIONS: Bacterial endosymbiont genes were identified in this dataset, some of which have a copy number consistent with being lateral gene transfers between endosymbiont genomes and Hemiptera, including ankyrin-repeat related proteins, lysozyme, and mannanase. Such genes and endosymbionts may provide novel targets for BMSB-specific biocontrol. This study demonstrates the utility of strand-specific sequencing in generating shotgun transcriptomes and that rapid sequencing shotgun transcriptomes is possible without the need for extensive inbreeding to generate homozygous lines. Such sequencing can provide a rapid response to pest invasions similar to that already described for disease epidemiology.

Extensively duplicated and transcriptionally active recent lateral gene transfer from a bacterial Wolbachia endosymbiont to its host filarial nematode Brugia malayi.

BACKGROUND: Lymphatic filariasis is a neglected tropical disease afflicting more than 120 million people, while another 1.3 billion people are at risk of infection. The nematode worm Brugia malayi is one of the causative agents of the disease and exists in a mutualistic symbiosis with Wolbachia bacteria. Since extensive lateral gene transfer occurs frequently between Wolbachia and its hosts, we sought to measure the extent of such LGT in B. malayi by whole genome sequencing of Wolbachia-depleted worms. RESULTS: A considerable fraction (at least 115.4-kbp, or 10.6%) of the 1.08-Mbp Wolbachia wBm genome has been transferred to its nematode host and retains high levels of similarity, including 227 wBm genes and gene fragments. Complete open reading frames were transferred for 32 of these genes, meaning they have the potential to produce functional proteins. Moreover, four transfers have evidence of life stage-specific regulation of transcription at levels similar to other nematode transcripts, strengthening the possibility that they are functional. CONCLUSIONS: There is extensive and ongoing transfer of Wolbachia DNA to the worm genome and some transfers are transcribed in a stage-specific manner at biologically relevant levels.

Genetic deletion of trkB.T1 increases neuromuscular function.

Neurotrophin-dependent activation of the tyrosine kinase receptor trkB.FL modulates neuromuscular synapse maintenance and function; however, it is unclear what role the alternative splice variant, truncated trkB (trkB.T1), may have in the peripheral neuromuscular axis. We examined this question in trkB.T1 null mice and demonstrate that in vivo neuromuscular performance and nerve-evoked muscle tension are significantly increased. In vitro assays indicated that the gain-in-function in trkB.T1(-/-) animals resulted specifically from an increased muscle contractility, and increased electrically evoked calcium release. In the trkB.T1 null muscle, we identified an increase in Akt activation in resting muscle as well as a significant increase in trkB.FL and Akt activation in response to contractile activity. On the basis of these findings, we conclude that the trkB signaling pathway might represent a novel target for intervention across diseases characterized by deficits in neuromuscular function.

Successful Profiling of Plasmodium falciparum var Gene Expression in Clinical Samples via a Custom Capture Array.

var genes encode Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) antigens. These highly diverse antigens are displayed on the surface of infected erythrocytes and play a critical role in immune evasion and sequestration of infected erythrocytes. Studies of var expression using non-leukocyte-depleted blood are challenging because of the predominance of host genetic material and lack of conserved var segments. Our goal was to enrich for parasite RNA, allowing de novo assembly of var genes and detection of expressed novel variants. We used two overall approaches: (i) enriching for total mRNA in the sequencing library preparations and (ii) enriching for parasite RNA with a custom capture array based on Roche's SeqCap EZ enrichment system. The capture array was designed with probes based on the whole 3D7 reference genome and an additional >4,000 full-length var gene sequences from other P. falciparum strains. We tested each method on RNA samples from Malian children with severe or uncomplicated malaria infections. All reads mapping to the human genome were removed, the remaining reads were assembled de novo into transcripts, and from these, var-like transcripts were identified and annotated. The capture array produced the longest maximum length and largest numbers of var gene transcripts in each sample, particularly in samples with low parasitemia. Identifying the most-expressed var gene sequences in whole-blood clinical samples without the need for extensive processing or generating sample-specific reference genome data is critical for understanding the role of PfEMP1s in malaria pathogenesis. IMPORTANCE Malaria parasites display antigens on the surface of infected red blood cells in the human host that facilitate attachment to blood vessels, contributing to the severity of infection. These antigens are highly variable, allowing the parasite to evade the immune system. Identifying these expressed antigens is critical to understanding the development of severe malarial disease. However, clinical samples contain limited amounts of parasite genetic material, a challenge for sequencing efforts further compounded by the extreme diversity of the parasite surface antigens. We present a method that enriches for these antigen sequences in clinical samples using a custom capture array, requiring minimal processing in the field. While our results are focused on the malaria parasite Plasmodium falciparum, this approach has broad applicability to other highly diverse antigens from other parasites and pathogens such as those that cause giardiasis and leishmaniasis.

Nearly Complete Genome Sequence of Brugia pahangi FR3.

Brugia pahangi is a zoonotic parasite that is closely related to human-infecting filarial nematodes. Here, we report the nearly complete genome of Brugia pahangi, including assemblies of four autosomes and an X chromosome, with only seven gaps. The Y chromosome is still not completely assembled.

Complete Genome Sequence of wBp, the Wolbachia Endosymbiont of Brugia pahangi FR3.

Lymphatic filariasis is a devastating disease caused by filarial nematode roundworms, which contain obligate Wolbachia endosymbionts. Here, we assembled the genome of wBp, the Wolbachia endosymbiont of the filarial nematode Brugia pahangi, from Illumina, Pacific Biosciences, and Oxford Nanopore data. The complete, circular genome is 1,072,967 bp.

Ultrahigh-Throughput Multiplexing and Sequencing of >500-Base-Pair Amplicon Regions on the Illumina HiSeq 2500 Platform.

Amplification, sequencing, and analysis of the 16S rRNA gene affords characterization of microbial community composition. As this tool has become more popular and amplicon-sequencing applications have grown in the total number of samples, growth in sample multiplexing is becoming necessary while maintaining high sequence quality and sequencing depth. Here, modifications to the Illumina HiSeq 2500 platform are described which produce greater multiplexing capabilities and 300-bp paired-end reads of higher quality than those produced by the current Illumina MiSeq platform. To improve the feasibility and flexibility of this method, a 2-step PCR amplification protocol is also described that allows for targeting of different amplicon regions, and enhances amplification success from samples with low bacterial bioburden. IMPORTANCE Amplicon sequencing has become a popular and widespread tool for surveying microbial communities. Lower overall costs associated with high-throughput sequencing have made it a widely adopted approach, especially for projects that necessitate sample multiplexing to eliminate batch effect and reduced time to acquire data. The method for amplicon sequencing on the Illumina HiSeq 2500 platform described here provides improved multiplexing capabilities while simultaneously producing greater quality sequence data and lower per-sample cost relative to those of the Illumina MiSeq platform without sacrificing amplicon length. To make this method more flexible for various amplicon-targeted regions as well as improve amplification from low-biomass samples, we also present and validate a 2-step PCR library preparation method.

Polymorphisms in the transcription factor 7-like 2 (TCF7L2) gene are associated with type 2 diabetes in the Amish: replication and evidence for a role in both insulin secretion and insulin resistance.

Transcription factor 7-like 2 (TCF7L2) regulates genes involved in cell proliferation and differentiation. The TCF7L2 gene is located on chromosome 10q25 in a region of replicated linkage to type 2 diabetes. Recently, a microsatellite marker in intron 3 (DG10S478) and five correlated single nucleotide polymorphisms (SNPs) were identified in Icelandic individuals that showed strong association with type 2 diabetes, which was replicated in Danish and European-American cohorts. We genotyped four of the SNPs (rs7901695, rs7903146, rs11196205, and rs12255372) in Amish subjects with type 2 diabetes (n = 137), impaired glucose tolerance (IGT; n = 139), and normal glucose tolerance (NGT; n = 342). We compared genotype frequencies in subjects with type 2 diabetes with those with NGT and found marginal association for rs7901695 (P = 0.05; odds ratio [OR] 1.51); comparison between NGT control subjects and the combined type 2 diabetes/IGT case group showed strong association with rs7901695 and rs7903146 (P = 0.008-0.01; OR 1.53-1.57) and marginal association with rs11196205 and rs12255372 (P = 0.07 and P = 0.04, respectively). In an expanded set of 698 Amish subjects without diabetes, we found no association with insulin and glucose levels during a 3-h oral glucose tolerance test. We also genotyped these SNPs in nondiabetic, non-Amish subjects (n = 48), in whom intravenous glucose tolerance tests were performed, and found an association between rs7901695 and rs7903146 and insulin sensitivity (P = 0.003 and P = 0.005, respectively) and disposition index (P = 0.04 and P = 0.007, respectively). These data provide replicating evidence that variants in TCF7L2 increase the risk for type 2 diabetes and novel evidence that the variants likely influence both insulin secretion and insulin sensitivity.

Analysis of complete genome sequence and major surface antigens of Neorickettsia helminthoeca, causative agent of salmon poisoning disease.

Neorickettsia helminthoeca, a type species of the genus Neorickettsia, is an endosymbiont of digenetic trematodes of veterinary importance. Upon ingestion of salmonid fish parasitized with infected trematodes, canids develop salmon poisoning disease (SPD), an acute febrile illness that is particularly severe and often fatal in dogs without adequate treatment. We determined and analysed the complete genome sequence of N. helminthoeca: a single small circular chromosome of 884 232 bp encoding 774 potential proteins. N. helminthoeca is unable to synthesize lipopolysaccharides and most amino acids, but is capable of synthesizing vitamins, cofactors, nucleotides and bacterioferritin. N. helminthoeca is, however, distinct from majority of the family Anaplasmataceae to which it belongs, as it encodes nearly all enzymes required for peptidoglycan biosynthesis, suggesting its structural hardiness and inflammatory potential. Using sera from dogs that were experimentally infected by feeding with parasitized fish or naturally infected in southern California, Western blot analysis revealed that among five predicted N. helminthoeca outer membrane proteins, P51 and strain-variable surface antigen were uniformly recognized. Our finding will help understanding pathogenesis, prevalence of N. helminthoeca infection among trematodes, canids and potentially other animals in nature to develop effective SPD diagnostic and preventive measures. Recent progresses in large-scale genome sequencing have been uncovering broad distribution of Neorickettsia spp., the comparative genomics will facilitate understanding of biology and the natural history of these elusive environmental bacteria.

Genetic variation in adiponectin receptor 1 and adiponectin receptor 2 is associated with type 2 diabetes in the Old Order Amish.

Adiponectin receptor 1 (ADIPOR1) and adiponectin receptor 2 (ADIPOR2) are newly identified receptors for adiponectin, an adipocytokine with anti-inflammatory and insulin-sensitizing properties. We screened for polymorphisms by performing sequence analysis on all eight exons, splice junctions, and approximately 2 kb of the 5' flanking regions of each receptor. We detected 5 single nucleotide polymorphisms (SNPs) in ADIPOR1 and 16 SNPs in ADIPOR2. We genotyped these SNPs in Amish subjects with type 2 diabetes (n = 137), impaired glucose tolerance (IGT) (n = 139), and normal glucose tolerance (n = 342) to test for association with type 2 diabetes. Three intronic SNPs in ADIPOR1 were significantly associated with type 2 diabetes (P = 0.014-0.007; odds ratio [OR] 1.61-1.65) and in high linkage disequilibrium (r2 = 0.97-1.0). In ADIPOR2, we found that five SNPs delineated one large haplotype block (r2= 0.9-1.0) spanning >98 kb of the gene and promoter region, which was strongly associated with the combined type 2 diabetes/IGT trait (P < or = 0.001; OR 1.64-1.71). To our knowledge, these data provide the first evidence for association between variation in the adiponectin receptors and type 2 diabetes.

Linkage of plasma adiponectin levels to 3q27 explained by association with variation in the APM1 gene.

We performed a genome-wide linkage scan of plasma adiponectin levels in 569 nondiabetic participants in the Amish Family Diabetes Study. The highest logarithm of odds (LOD) score (2.13; P = 0.0009) occurred on chromosome 3q27 between markers D3S1602 and D3S1580, which flank APM1/ACDC, the adiponectin gene. The APM1 +2019 A/- insertion/deletion polymorphism in the 3' untranslated region (single nucleotide polymorphism [SNP] +2019; deletion allele frequency 0.30 in Amish) showed strong association with adiponectin levels in a dosage-dependent manner in a direction consistent with that reported in previous studies, with deletion heterozygosity increasing adiponectin levels by 1.3 +/- 0.5 microg/ml and deletion homozygosity increasing levels by 3.0 +/- 0.8 microg/ml (P < 0.0001). Two other SNPs, rs2241766 and rs1501299, showed moderate association. In a subset of 523 subjects genotyped for both SNP +2019 and rs2241766, including the APM1 SNP +2019 genotype as a covariate reduced the linkage signal at 3q27 by 1.26 LOD units (from 2.22 to 0.96) and including both SNPs reduced the signal by 1.51 LOD units (to 0.71). These findings, combined with a two-point LOD score of 2.35 for SNP +2019, provide evidence that variation in APM1 is responsible for linkage of adiponectin levels to 3q27 in the Old Order Amish.

Efficient Enrichment of Bacterial mRNA from Host-Bacteria Total RNA Samples.

Despite numerous advances in genomics and bioinformatics, technological hurdles remain to examine host-microbe transcriptomics. Sometimes the transcriptome of either or both can be ascertained merely by generating more sequencing reads. However, many cases exist where bacterial mRNA needs to be enriched further to enable cost-effective sequencing of the pathogen or endosymbiont. While a suitable method is commercially available for mammalian samples of this type, development of such methods has languished for invertebrate samples. Furthermore, a common method across multiple taxa would facilitate comparisons between bacteria in invertebrate vectors and their vertebrate hosts. Here, a method is described to concurrently remove polyadenylated transcripts, prokaryotic rRNA, and eukaryotic rRNA, including those with low amounts of starting material (e.g. 100 ng). In a Wolbachia-Drosophila system, this bacterial mRNA enrichment yielded a 3-fold increase in Wolbachia mRNA abundance and a concomitant 3.3-fold increase in the percentage of transcripts detected. More specifically, 70% of the genome could be recovered by transcriptome sequencing compared to 21% in the total RNA. Sequencing of similar bacterial mRNA-enriched samples generated from Ehrlichia-infected canine cells covers 93% of the Ehrlichia genome, suggesting ubiquitous transcription across the entire Ehrlichia chaffeensis genome. This technique can potentially be used to enrich bacterial mRNA in many studies of host-microbe interactions.