Christopher P. Kolbert

Ribosomal DNA sequencing as a tool for identification of bacterial pathogens.

Conventional methods for the identification and characterization of clinical isolates of bacterial pathogens sometimes fall short when such isolates exhibit unusual phenotypic profiles. Recent advances in DNA sequencing technology have greatly enhanced the ability of the microbiologist to determine the identity of a bacterial isolate. Given the relative objectivity of DNA sequence information and growing availability of sequence information databases, a significant movement is now afoot to use molecular methods for the identification of clinical pathogens.

Brain Expression Genome-Wide Association Study (eGWAS) Identifies Human Disease-Associated Variants

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimers disease (AD, cerebellar n = 197, temporal cortex n = 202) and with other brain pathologies (non–AD, cerebellar n = 177, temporal cortex n = 197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ±100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non–ADs (q<0.05, p = 7.70×10−5–1.67×10−82). Of these, 2,089 were also significant in the temporal cortex (p = 1.85×10−5–1.70×10−141). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10−6). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinsons disease (PD) MMRN1/rs6532197, Pagets disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9–3.3 fold enrichment (p<10−6) of significant cisSNPs with suggestive AD–risk association (p<10−3) in the Alzheimers Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non–CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.

Bordetella holmesii-Like Organisms Associated with Septicemia, Endocarditis, and Respiratory Failure

We recovered an unusual bacterial strain from blood or sputum of three patients with septicemia, endocarditis, and/or respiratory failure. The three isolates were thin, curved, gram-negative, light brown, pigment-producing bacilli with variable catalase activity. They were asaccharolytic, oxidase-negative, nonmotile, and fastidious. Identification was not possible on the basis of these characteristics alone or in combination with cellular fatty acid profiles. Nucleic acid amplification and sequence analysis of the 16S rRNA gene revealed that all three isolates were identical and most closely related to the emerging pathogen Bordetella holmesii, diverging from the published sequence at three nucleotide positions (99.8% similarity). Isolation of a B. holmesii-like pathogen from sputum suggests that, in addition to producing septicemia, the organism may inhabit the respiratory tract like other Bordetella species.

Novel late-onset Alzheimer disease loci variants associate with brain gene expression

Objective: Recent genome-wide association studies (GWAS) of late-onset Alzheimer disease (LOAD) identified 9 novel risk loci. Discovery of functional variants within genes at these loci is required to confirm their role in Alzheimer disease (AD). Single nucleotide polymorphisms that influence gene expression (eSNPs) constitute an important class of functional variants. We therefore investigated the influence of the novel LOAD risk loci on human brain gene expression. Methods: We measured gene expression levels in the cerebellum and temporal cortex of autopsied AD subjects and those with other brain pathologies (∼400 total subjects). To determine whether any of the novel LOAD risk variants are eSNPs, we tested their cis-association with expression of 6 nearby LOAD candidate genes detectable in human brain (ABCA7, BIN1, CLU, MS4A4A, MS4A6A, PICALM) and an additional 13 genes ±100 kb of these SNPs. To identify additional eSNPs that influence brain gene expression levels of the novel candidate LOAD genes, we identified SNPs ±100 kb of their location and tested for cis-associations. Results: CLU rs11136000 (p = 7.81 × 10−4) and MS4A4A rs2304933/rs2304935 (p = 1.48 × 10−4–1.86 × 10−4) significantly influence temporal cortex expression levels of these genes. The LOAD-protective CLU and risky MS4A4A locus alleles associate with higher brain levels of these genes. There are other cis-variants that significantly influence brain expression of CLU and ABCA7 (p = 4.01 × 10−5–9.09 × 10−9), some of which also associate with AD risk (p = 2.64 × 10−2–6.25 × 10−5). Conclusions: CLU and MS4A4A eSNPs may at least partly explain the LOAD risk association at these loci. CLU and ABCA7 may harbor additional strong eSNPs. These results have implications in the search for functional variants at the novel LOAD risk loci.

Quantitative miRNA expression analysis using fluidigm microfluidics dynamic arrays.

BackgroundMicroRNAs (miRNAs) represent a growing class of small non-coding RNAs that are important regulators of gene expression in both plants and animals. Studies have shown that miRNAs play a critical role in human cancer and they can influence the level of cell proliferation and apoptosis by modulating gene expression. Currently, methods for the detection and measurement of miRNA expression include small and moderate-throughput technologies, such as standard quantitative PCR and microarray based analysis. However, these methods have several limitations when used in large clinical studies where a high-throughput and highly quantitative technology needed for the efficient characterization of a large number of miRNA transcripts in clinical samples. Furthermore, archival formalin fixed, paraffin embedded (FFPE) samples are increasingly becoming the primary resource for gene expression studies because fresh frozen (FF) samples are often difficult to obtain and requires special storage conditions. In this study, we evaluated the miRNA expression levels in FFPE and FF samples as well as several lung cancer cell lines employing a high throughput qPCR-based microfluidic technology. The results were compared to standard qPCR and hybridization-based microarray platforms using the same samples.ResultsWe demonstrated highly correlated Ct values between multiplex and singleplex RT reactions in standard qPCR assays for miRNA expression using total RNA from A549 (R = 0.98; p < 0.0001) and H1299 (R = 0.95; p < 0.0001) lung cancer cell lines. The Ct values generated by the microfluidic technology (Fluidigm 48.48 dynamic array systems) resulted in a left-shift toward lower Ct values compared to those observed by ABI 7900 HT (mean difference, 3.79), suggesting that the microfluidic technology exhibited a greater sensitivity. In addition, we show that as little as 10 ng total RNA can be used to reliably detect all 48 or 96 tested miRNAs using a 96-multiplexing RT reaction in both FFPE and FF samples. Finally, we compared miRNA expression measurements in both FFPE and FF samples by qPCR using the 96.96 dynamic array and Affymetrix microarrays. Fold change comparisons for comparable genes between the two platforms indicated that the overall correlation was R = 0.60. The maximum fold change detected by the Affymetrix microarray was 3.5 compared to 13 by the 96.96 dynamic array.ConclusionThe qPCR-array based microfluidic dynamic array platform can be used in conjunction with multiplexed RT reactions for miRNA gene expression profiling. We showed that this approach is highly reproducible and the results correlate closely with the existing singleplex qPCR platform at a throughput that is 5 to 20 times higher and a sample and reagent usage that was approximately 50-100 times lower than conventional assays. We established optimal conditions for using the Fluidigm microfluidic technology for rapid, cost effective, and customizable arrays for miRNA expression profiling and validation.

Direct Identification of Bacteria from Positive Blood Cultures by Amplification and Sequencing of the 16S rRNA Gene: Evaluation of BACTEC 9240 Instrument True- Positive and False-Positive Results

ABSTRACT In a previous study which evaluated the BACTEC 9240 automated blood culture system (Becton Dickinson Diagnostic Instrument Systems, Sparks, Md.), we noted a 1.3% “instrument false-positive” rate. That is, the BACTEC system signaled that a bottle (BACTEC Plus Aerobic/F bottle or BACTEC Anaerobic Lytic/10 bottle) culture was positive but a Gram stain was negative and there was no growth of bacteria or yeasts on subculture to chocolate agar. Furthermore, from the same sample of blood, cultures for fungi using the Isolator blood culture system (Wampole Laboratories, Cranbury, N.J.) were negative for growth. For the present study, we evaluated 76 instrument false-positive samples for the presence of 16S ribosomal DNA using the MicroSeq 500 kit (PE Biosystems, Foster City, Calif.). These samples also were negative for fungi by the Isolator method. This kit has a PCR module and sequencing module for the amplification and sequencing of the 16S RNA gene and provides a database for sequence alignment and identification of bacteria. To optimize the assay, we evaluated the effect of adding 0.5% bovine serum albumin to the sample from blood culture bottles and found that it decreased the effects of inhibitors on the PCR. Two control groups of blood culture specimens were also evaluated. One group (n = 45) were “instrument true positives”; the instrument signaled positive, and subsequent Gram stains were positive and subcultures on chocolate agar grew bacteria. The other group (n = 20) were “instrument true negatives”; the instrument signaled negative, the Gram stain was negative, and subcultures on chocolate agar and from the Isolator tube on fungal media showed no growth. None of the 76 instrument false-positive samples had evidence for 16S rRNA gene sequences. All of the instrument true-positive samples and all of the instrument true-negative specimens were positive and negative, respectively, using the MicroSeq 500 kit. Total peripheral white blood cell counts were statistically significantly higher for patients who had instrument false-positive results than for patients who had instrument true-positive or true-negative results (P = 0.001). We conclude that instrument false positives signaled by the BACTEC 9240 system are not due to bacteria in the blood culture samples.

Multi-Platform Analysis of MicroRNA Expression Measurements in RNA from Fresh Frozen and FFPE Tissues

MicroRNAs play a role in regulating diverse biological processes and have considerable utility as molecular markers for diagnosis and monitoring of human disease. Several technologies are available commercially for measuring microRNA expression. However, cross-platform comparisons do not necessarily correlate well, making it difficult to determine which platform most closely represents the true microRNA expression level in a tissue. To address this issue, we have analyzed RNA derived from cell lines, as well as fresh frozen and formalin-fixed paraffin embedded tissues, using Affymetrix, Agilent, and Illumina microRNA arrays, NanoString counting, and Illumina Next Generation Sequencing. We compared the performance within- and between the different platforms, and then verified these results with those of quantitative PCR data. Our results demonstrate that the within-platform reproducibility for each method is consistently high and although the gene expression profiles from each platform show unique traits, comparison of genes that were commonly detectable showed that detection of microRNA transcripts was similar across multiple platforms.

Increased miR-708 Expression in NSCLC and Its Association with Poor Survival in Lung Adenocarcinoma from Never Smokers

Purpose: miRNA plays an important role in human disease and cancer. We seek to investigate the expression status, clinical relevance, and functional role of miRNA in non–small cell lung cancer. Experimental Design: We conducted miRNA expression profiling in matched lung adenocarcinoma and uninvolved lung using 56 pairs of fresh-frozen (FF) and 47 pairs of formalin-fixed, paraffin-embedded (FFPE) samples from never smokers. The most differentially expressed miRNA genes were evaluated by Cox analysis and log-rank test. Among the best candidate, miR-708 was further examined for differential expression in two independent cohorts. Functional significance of miR-708 expression in lung cancer was examined by identifying its candidate mRNA target and through manipulating its expression levels in cultured cells. Results: Among the 20 miRNAs most differentially expressed between tested tumor and normal samples, high expression level of miR-708 in the tumors was most strongly associated with an increased risk of death after adjustments for all clinically significant factors including age, sex, and tumor stage (FF cohort: HR, 1.90; 95% CI, 1.08–3.35; P = 0.025 and FFPE cohort: HR, 1.93; 95% CI, 1.02–3.63; P = 0.042). The transcript for TMEM88 gene has a miR-708 binding site in its 3′ UTR and was significantly reduced in tumors high of miR-708. Forced miR-708 expression reduced TMEM88 transcript levels and increased the rate of cell proliferation, invasion, and migration in culture. Conclusions: miRNA-708 acts as an oncogene contributing to tumor growth and disease progression by directly downregulating TMEM88, a negative regulator of the Wnt signaling pathway in lung cancer. Clin Cancer Res; 18(13); 3658–67. ©2012 AACR.

Key inflammatory signaling pathways are regulated by the proteasome.

ABSTRACT Lipopolysaccharide (LPS) is a major structural component of all Gram-negative organisms and has been implicated in Gram-negative sepsis and septic shock. In the present study, Affymetrix microarray analysis of RNA derived from murine macrophages treated with LPS in the absence or presence of the proteasome inhibitor lactacystin revealed that the vast majority of genes regulated by LPS is under control of the proteasome. Analysis of the data has revealed that the products of these genes participate in 14 distinct signaling pathways. This represents a novel approach to the identification of signaling pathways that are both toll-like receptor 4- and proteasome-dependent and may lead to the development of new drug targets in Gram-negative sepsis and septic shock.

Association of MAPT haplotypes with Alzheimer’s disease risk and MAPT brain gene expression levels

IntroductionMAPT encodes for tau, the predominant component of neurofibrillary tangles that are neuropathological hallmarks of Alzheimer’s disease (AD). Genetic association of MAPT variants with late-onset AD (LOAD) risk has been inconsistent, although insufficient power and incomplete assessment of MAPT haplotypes may account for this.MethodsWe examined the association of MAPT haplotypes with LOAD risk in more than 20,000 subjects (n-cases = 9,814, n-controls = 11,550) from Mayo Clinic (n-cases = 2,052, n-controls = 3,406) and the Alzheimer’s Disease Genetics Consortium (ADGC, n-cases = 7,762, n-controls = 8,144). We also assessed associations with brain MAPT gene expression levels measured in the cerebellum (n = 197) and temporal cortex (n = 202) of LOAD subjects. Six single nucleotide polymorphisms (SNPs) which tag MAPT haplotypes with frequencies greater than 1% were evaluated.ResultsH2-haplotype tagging rs8070723-G allele associated with reduced risk of LOAD (odds ratio, OR = 0.90, 95% confidence interval, CI = 0.85-0.95, p = 5.2E-05) with consistent results in the Mayo (OR = 0.81, p = 7.0E-04) and ADGC (OR = 0.89, p = 1.26E-04) cohorts. rs3785883-A allele was also nominally significantly associated with LOAD risk (OR = 1.06, 95% CI = 1.01-1.13, p = 0.034). Haplotype analysis revealed significant global association with LOAD risk in the combined cohort (p = 0.033), with significant association of the H2 haplotype with reduced risk of LOAD as expected (p = 1.53E-04) and suggestive association with additional haplotypes. MAPT SNPs and haplotypes also associated with brain MAPT levels in the cerebellum and temporal cortex of AD subjects with the strongest associations observed for the H2 haplotype and reduced brain MAPT levels (β = -0.16 to -0.20, p = 1.0E-03 to 3.0E-03).ConclusionsThese results confirm the previously reported MAPT H2 associations with LOAD risk in two large series, that this haplotype has the strongest effect on brain MAPT expression amongst those tested and identify additional haplotypes with suggestive associations, which require replication in independent series. These biologically congruent results provide compelling evidence to screen the MAPT region for regulatory variants which confer LOAD risk by influencing its brain gene expression.