Search

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.

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.

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.

Lymphatic filariasis is a devastating disease caused by filarial nematode roundworms, which contain obligate Wolbachia endosymbionts. Here, we assembled the genome of w Bp, 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.

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.

Summary 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.

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., 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.

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 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, thus improving amplification success from low bacterial bioburden samples.ImportanceAmplicon 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 which 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 the Illumina MiSeq platform, without sacrificing amplicon length. To make this method more flexible to 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.

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.

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., Author summary Since its discovery in 1881, invasive disease caused by Streptococcus pneumoniae, the leading cause of community-acquired pneumonia and a prototypical extracellular pathogen, has been tied exclusively to the planktonic phenotype, i.e. individual diplococci or short chains. Herein, we report that heart-invaded pneumococci can instead replicate intracellularly and transition into an immunoquiescent biofilm. Using dual RNA-seq technology we capture the complete gene expression profile of S. pneumoniae within infected hearts as well as the bloodstream. In doing so, we affirm that pneumococci within the heart are indeed forming biofilms and identify virulence determinants that play important roles in vivo. Accordingly, we identify a novel role for the genomic island Region of Diversity 12 in promotion of biofilm formation, virulence, and dampening of the host response. We subsequently show that biofilm pneumococci prevent immune cell infiltration into the heart not by passive means but instead through enhanced fratricide-mediated release of the toxin pneumolysin that kills resident macrophages, pre-empting their response. Collectively our manuscript describes a novel site, i.e. intracellular; previously unreported growth phenotype during invasive disease, i.e. biofilm formation; and counter-intuitive molecular mechanism, i.e. resident macrophage killing; for how pneumococci establish themselves in the heart without inflammation.

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 recently introduced bacterial species Chlamydia gallinacea is known to occur in domestic poultry and other birds. Its potential as an avian pathogen and zoonotic agent is under investigation. The whole-genome sequence of its type strain, 08-1274/3, consists of a 1,059,583-bp chromosome with 914 protein-coding sequences (CDSs) and a plasmid (p1274) comprising 7,619 bp with 9 CDSs.

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.

Apolipoprotein C-III (apoC-III) inhibits triglyceride hydrolysis and has been implicated in coronary artery disease. Through a genome-wide association study, we have found that about 5% of the Lancaster Amish are heterozygous carriers of a null mutation (R19X) in the gene encoding apoC-III ( APOC3 ) and, as a result, express half the amount of apoC-III present in noncarriers. Mutation carriers compared with noncarriers had lower fasting and postprandial serum triglycerides, higher levels of HDL-cholesterol and lower levels of LDL-cholesterol. Subclinical atherosclerosis, as measured by coronary artery calcification, was less common in carriers than noncarriers, which suggests that lifelong deficiency of apoC-III has a cardioprotective effect.

OBJECTIVE— We sought to identify type 2 diabetes susceptibility genes through a genome-wide association scan (GWAS) in the Amish. RESEARCH DESIGN AND METHODS— DNA from 124 type 2 diabetic case subjects and 295 control subjects with normal glucose tolerance were genotyped on the Affymetrix 100K single nucleotide polymorphism (SNP) array. A total of 82,485 SNPs were tested for association with type 2 diabetes. Type 2 diabetes–associated SNPs were further prioritized by the following: 1) associations with 5 oral glucose tolerance test (OGTT) traits in 427 nondiabetic Amish subjects, and 2) in silico replication from three independent 100L SNP GWASs (Framingham Heart Study Caucasians, Pima Indians, and Mexican Americans) and a 500K GWAS in Scandinavians. RESULTS— The strongest association (P = 1.07 × 10−5) was for rs2237457, which is located in growth factor receptor–bound protein 10 (Grb10), an adaptor protein that regulate insulin receptor signaling. rs2237457 was also strongly associated with OGTT glucose area under the curve in nondiabetic subjects (P = 0.001). Of the 1,093 SNPs associated with type 2 diabetes at P < 0.01, 67 SNPs demonstrated associations with at least one OGTT trait in nondiabetic individuals; 80 SNPs were nominally associated with type 2 diabetes in one of the three independent 100K GWASs, 3 SNPs (rs2540317 in MFSD9, rs10515353 on chromosome 5, and rs2242400 in BCAT1 were associated with type 2 diabetes in more than one population), and 11 SNPs were nominally associated with type 2 diabetes in Scandinavians. One type 2 diabetes–associated SNP (rs3845971, located in FHIT) showed replication with OGTT traits and also in another population. CONCLUSIONS— Our GWAS of type 2 diabetes identified several gene variants associated with type 2 diabetes, some of which are worthy of further study.

Rho guanine nucleotide exchange factor 11 (ARHGEF11), located on chromosome 1q21, is involved in G protein signaling and is a pathway known to play a role in both insulin secretion and action. We genotyped 52 single nucleotide polymorphims (SNPs) in ARHGEF11 and compared the genotype frequencies of subjects with type 2 diabetes (n = 145) or type 2 diabetes/impaired glucose tolerance (IGT) (n = 293) with those of control subjects with normal glucose tolerance (NGT) (n = 358). Thirty SNPs, spanning the entire gene, were significantly associated with type 2 diabetes or type 2 diabetes/IGT. The most significantly associated SNP was rs6427340 (intron 2), in which the less common allele was the risk allele (odds ratio [OR] 1.82 [95% CI 1.20–2.70], P = 0.005 for type 2 diabetes vs. NGT and 1.79 [1.27–2.50], P = 0.0008 for type 2 diabetes/IGT vs. NGT). In an expanded set of nondiabetic subjects (n = 754), most of the type 2 diabetes–and IGT-associated SNPs were significantly associated with glucose levels during an oral glucose tolerance test, with the same SNP (rs6427340) showing the most significant associations (P = 0.007). All type 2 diabetes–and IGT-associated SNPs were in high linkage disequilibrium and constitute a single 133-kb haplotype block. These results, coupled with similar findings in Pima Indians, suggest that sequence variation in ARHGEF11 may influence risk of type 2 diabetes.

Background: Nucleosomes are the building blocks of chromatin where gene regulation takes place. Chromatin landscapes have been profiled for several species, providing insights into the fundamental mechanisms of chromatin-mediated transcriptional regulation of gene expression. However, knowledge is missing for several major and deep-branching eukaryotic groups, such as the Stramenopiles, which include the diatoms. Diatoms are highly diverse and ubiquitous species of phytoplankton that play a key role in global biogeochemical cycles. Dissecting chromatin-mediated regulation of genes in diatoms will help understand the ecological success of these organisms in contemporary oceans. Results: Here, we use high resolution mass spectrometry to identify a full repertoire of post-translational modifications on histones of the marine diatom Phaeodactylum tricornutum, including eight novel modifications. We map five histone marks coupled with expression data and show that P. tricornutum displays both unique and broadly conserved chromatin features, reflecting the chimeric nature of its genome. Combinatorial analysis of histone marks and DNA methylation demonstrates the presence of an epigenetic code defining activating or repressive chromatin states. We further profile three specific histone marks under conditions of nitrate depletion and show that the histone code is dynamic and targets specific sets of genes. Conclusions: This study is the first genome-wide characterization of the histone code from a stramenopile and a marine phytoplankton. The work represents an important initial step for understanding the evolutionary history of chromatin and how epigenetic modifications affect gene expression in response to environmental cues in marine environments. Instituto de Biotecnología Fil: Veluchamy, Alaguraj. Centre National de la Recherche Scientifique. Institut de Biologie de l’École Normale Supérieure. Ecology and Evolutionary Biology Section; Francia. King Abdullah University of Science and Technology. BESE Division. Center for Desert Agriculture; Arabia Saudita Fil: Rastogi, Achal. Centre National de la Recherche Scientifique. Institut de Biologie de l’École Normale Supérieure. Ecology and Evolutionary Biology Section; Francia Fil: Lin, Xin. Centre National de la Recherche Scientifique. Institut de Biologie de l’École Normale Supérieure. Ecology and Evolutionary Biology Section; Francia. Xiamen University. State key lab of Marine Environmental Science; China Fil: Lombard, Bérangère. Institut Curie PSL Research University. Centre de Recherche. Laboratoire de Spectrométrie de Masse Protéomique; Francia Fil: Murik, Omer. Centre National de la Recherche Scientifique. Institut de Biologie de l’École Normale Supérieure. Ecology and Evolutionary Biology Section; Francia Fil: Thomas, Yann. Centre National de la Recherche Scientifique. Institut de Biologie de l’École Normale Supérieure. Ecology and Evolutionary Biology Section; Francia Fil: Dingli, Florent. Institut Curie PSL Research University. Centre de Recherche. Laboratoire de Spectrométrie de Masse Protéomique; Francia Fil: Rivarola, Maximo Lisandro. Instituto Nacional de Tecnología Agropecuaria (INTA). Instituto de Biotecnología; Argentina. University of Maryland School of Medicine. Institute for Genome Sciences; Estados Unidos Fil: Ott, Sandra. University of Maryland School of Medicine. Institute for Genome Sciences; Estados Unidos Fil: Liu, Xinyue. University of Maryland School of Medicine. Institute for Genome Sciences; Estados Unidos Fil: Sun, Yezhou. University of Maryland School of Medicine. Institute for Genome Sciences; Estados Unidos Fil: Rabinowicz, Pablo D. University of Maryland School of Medicine. Institute for Genome Sciences; Estados Unidos Fil: McCarthy, James. J. Craig Venter Institute; Estados Unidos Fil: Allen, Andrew E. J. Craig Venter Institute; Estados Unidos. University of California. Integrative Oceanography Division. Scripps Institution of Oceanography; Estados Unidos Fil: Loew, Damarys. Institut Curie PSL Research University. Centre de Recherche. Laboratoire de Spectrométrie de Masse Protéomique; Francia Fil: Bowler, Chris. Centre National de la Recherche Scientifique. Institut de Biologie de l’École Normale Supérieure. Ecology and Evolutionary Biology Section; Francia Fil: Tirichine, Leïla. Centre National de la Recherche Scientifique. Institut de Biologie de l’École Normale Supérieure. Ecology and Evolutionary Biology Section; Francia

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.

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 μg/ml and deletion homozygosity increasing levels by 3.0 ± 0.8 μg/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.

Calsequestrin (CASQ)1 is involved in intracellular storage and release of calcium, a process that has been shown to mediate glucose transport in muscle. Its gene, CASQ1, is encoded on chromosome 1q21, a region that has been linked to type 2 diabetes in the Amish and several other populations. We screened all 11 exons, exon-intron junctions, and the proximal regulatory region of CASQ1 for mutations. We detected four novel single nucleotide polymorphisms (SNPs) (−1470C→T, −1456delG, −1366insG, and 593C→T). Ten informative SNPs within CASQ1 were genotyped in Amish subjects with type 2 diabetes (n = 145), impaired glucose tolerance (n = 148), and normal glucose tolerance (n = 358). Rs2275703 and rs617698 in introns 4 and 2 were significantly associated with type 2 diabetes (P = 0.008 and 0.04, respectively); three other SNPs showed borderline evidence for association to type 2 diabetes (P = 0.076–0.093). Furthermore, in nondiabetic subjects (n = 754), both rs2275703 and rs617698 were significantly associated with glucose area under the curve during an oral glucose tolerance test (P = 0.035 and 0.013, respectively). Haplotype analysis suggested that no haplotype could explain these associations better than rs2275703. These findings, coupled with similar findings in Utah Caucasians, suggest that sequence variation in CASQ1 may influence risk of type 2 diabetes.

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. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-788) contains supplementary material, which is available to authorized users.

Macrolide resistance, emerging in Streptococcus pneumoniae and other Gram-positive bacteria, is increasingly due to efflux pumps encoded by mef/mel(msr) operons found on discrete mobile genetic elements. The regulation of mef/mel(msr) in these elements is not well understood. We identified the mef(E)/mel transcriptional start, localized the mef(E)/mel promoter, and demonstrated attenuation of transcription as a mechanism of regulation of macrolide-inducible mef-mediated macrolide resistance in S. pneumoniae. The mef(E)/mel transcriptional start site was a guanine 327 bp upstream of mef(E). Consensus pneumococcal promoter -10 (5'-TATACT-3') and -35 (5'-TTGAAC-3') boxes separated by 17 bp were identified 7 bp upstream of the start site. Analysis of the predicted secondary structure of the 327 5' region identified four pairs of inverted repeats R1-R8 predicted to fold into stem-loops, a small leader peptide [MTASMRLR, (Mef(E)L)] required for macrolide induction and a Rho-independent transcription terminator. RNA-seq analyses provided confirmation of transcriptional attenuation. In addition, expression of mef(E)L was also influenced by mef(E)L-dependent mRNA stability. The regulatory region 5' of mef(E) was highly conserved in other mef/mel(msr)-containing elements including Tn1207.1 and the 5612IQ complex in pneumococci and Tn1207.3 in Group A streptococci, indicating a regulatory mechanism common to a wide variety of Gram-positive bacteria containing mef/mel(msr) elements.

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. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-738) contains supplementary material, which is available to authorized users.

Cassava is a major tropical food crop in the Euphorbiaceae family that has high carbohydrate production potential and adaptability to diverse environments. Here we present the draft genome sequences of a wild ancestor and a domesticated variety of cassava and comparative analyses with a partial inbred line. We identify 1,584 and 1,678 gene models specific to the wild and domesticated varieties, respectively, and discover high heterozygosity and millions of single-nucleotide variations. Our analyses reveal that genes involved in photosynthesis, starch accumulation and abiotic stresses have been positively selected, whereas those involved in cell wall biosynthesis and secondary metabolism, including cyanogenic glucoside formation, have been negatively selected in the cultivated varieties, reflecting the result of natural selection and domestication. Differences in microRNA genes and retrotransposon regulation could partly explain an increased carbon flux towards starch accumulation and reduced cyanogenic glucoside accumulation in domesticated cassava. These results may contribute to genetic improvement of cassava through better understanding of its biology., Cassava is a major source of food in tropical and subtropical regions. Here the authors sequence the genomes of wild and domesticated cassava varieties and identify genes that have been selected for and against during the evolution and domestication of this economically important crop.

Background To take advantage of affordable high-throughput next-generation sequencing technologies to characterize microbial community composition often requires the development of improved methods to overcome technical limitations inherent to the sequencing platforms. Sequencing low sequence diversity libraries such as 16S rRNA amplicons has been problematic on the Illumina MiSeq platform and often generates sequences of suboptimal quality. Results Here we present an improved dual-indexing amplification and sequencing approach to assess the composition of microbial communities from clinical samples using the V3-V4 region of the 16S rRNA gene on the Illumina MiSeq platform. We introduced a 0 to 7 bp “heterogeneity spacer” to the index sequence that allows an equal proportion of samples to be sequenced out of phase. Conclusions Our approach yields high quality sequence data from 16S rRNA gene amplicons using both 250 bp and 300 bp paired-end MiSeq protocols and provides a flexible and cost-effective sequencing option.

We report the draft genome sequence of Mortierella alpina isolate CDC-B6842. M. alpina is a nonpathogenic member of the Mucoromycotina subphylum of fungi that is an important model for understanding the molecular mechanisms of lipid production and metabolism.

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 w Bm genome has been transferred to its nematode host and retains high levels of similarity, including 227 w Bm 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.

TAS1R- and TAS2R-type taste receptors are expressed in the gustatory system, where they detect sweet- and bitter-tasting stimuli, respectively. These receptors are also expressed in subsets of cells within the mammalian gastrointestinal tract, where they mediate nutrient assimilation and endocrine responses. For example, sweeteners stimulate taste receptors on the surface of gut enteroendocrine L cells to elicit an increase in intracellular Ca(2+) and secretion of the incretin hormone glucagon-like peptide-1 (GLP-1), an important modulator of insulin biosynthesis and secretion. Because of the importance of taste receptors in the regulation of food intake and the alimentary responses to chemostimuli, we hypothesized that differences in taste receptor efficacy may impact glucose homeostasis. To address this issue, we initiated a candidate gene study within the Amish Family Diabetes Study and assessed the association of taste receptor variants with indicators of glucose dysregulation, including a diagnosis of type 2 diabetes mellitus and high levels of blood glucose and insulin during an oral glucose tolerance test. We report that a TAS2R haplotype is associated with altered glucose and insulin homeostasis. We also found that one SNP within this haplotype disrupts normal responses of a single receptor, TAS2R9, to its cognate ligands ofloxacin, procainamide and pirenzapine. Together, these findings suggest that a functionally compromised TAS2R receptor negatively impacts glucose homeostasis, providing an important link between alimentary chemosensation and metabolic disease.

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 ∼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 ≤ 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.

Hepatocyte nuclear factor 4-α (HNF4A) is a transcription factor located on chromosome 20q13 that regulates expression of genes involved in glucose metabolism and homeostasis. Recently, two groups independently identified single nucleotide polymorphism (SNPs) in an alternate upstream promoter (P2) of HNF4A that were associated with type 2 diabetes in Ashkenazi Jews and Finns. We genotyped haplotype-tagging SNPs (htSNPs) across the two promoter regions and the coding region of HNF4A in individuals with type 2 diabetes (n = 137), impaired glucose tolerance (IGT) (n = 139), and normal glucose tolerance (n = 342) from the Amish Family Diabetes Study (AFDS) to test for association with type 2 diabetes. In the P1 promoter region, we observed a significant association between the A allele of rs2425640 and type 2 diabetes (odds ratio [OR] 1.60, P = 0.03). Furthermore, the mean age of type 2 diabetes onset was, on average, 5.1 years earlier in those with the AA or GA genotype at SNP rs2425640 than in those with the GG genotype (57.8 vs. 62.9 years, P = 0.011). In the P2 promoter, the htSNP rs1884614 showed borderline association with both type 2 diabetes (OR 1.40, P = 0.09) and the combined type 2 diabetes/IGT trait (1.35, P = 0.07). In an expanded set of 698 nondiabetic AFDS subjects, we found association between rs1884614 and glucose area under the curve during an oral glucose tolerance test (additive model, P = 0.022; dominant model, P = 0.010). The results of this study provide evidence that variants in both the P1 and P2 promoters of HNF4A increase risk for typical type 2 diabetes.

In July 1940 the Armistice Commission sent a German delegation to Pau in south-western France. The delegation included Andre Muller, a Nazi who also worked for the German Security Police. Muller twenty-nine letters to a female collaborator, alongside the testimonies of those who knew the two friends, provide an unusual opportunity to recapture a sense of day-to-day relations between the occupied and occupiers. The Nazi officer's letters reveal his reflections on the war, his hopes and, above all, his material needs. The cycle of epistolary exchange that linked the two friends entailed a unidirectional flow of gifts from the female collaborator to Muller. The letters, goods and services they exchanged raise important questions about the location of power, indebtedness and gratitude in Franco-German partnerships; they also provide an opportunity to test assumptions about the inequality of such relationships in occupied France. The court testimonies in Pau also offer rare insights into female involvement in collaborationist movements, which were primarily a male sphere of activity and influence.