Highly parallel genome-wide expression analysis of single mammalian cells

Background

We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs.

Methodology/Principal Findings

The method uses tagged priming and template-switching, resulting in the incorporation of universal PCR priming sites at both ends of the synthesized cDNA for global PCR amplification. Coupled with a whole-genome gene expression microarray platform, we routinely obtain expression correlation values of R²∼0.76–0.80 between individual cells and R²∼0.69 between 50 pg total RNA replicates. Expression profiles generated from single cells or 50 pg total RNA correlate well with that generated with higher input (1 ng total RNA) (R²∼0.80). Also, the assay is sufficiently sensitive to detect, in a single cell, approximately 63% of the number of genes detected with 1 ng input, with approximately 97% of the genes detected in the single-cell input also detected in the higher input.

Conclusions/Significance

In summary, our method facilitates whole-genome gene expression profiling in contexts where starting material is extremely limiting, particularly in areas such as the study of progenitor cells in early development and tumor stem cell biology.     

Activated TLR signaling in atherosclerosis among women with lower Framingham risk score: the multi-ethnic study of atherosclerosis

Background

Atherosclerosis is the leading cause of cardiovascular disease (CVD). Traditional risk factors can be used to identify individuals at high risk for developing CVD and are generally associated with the extent of atherosclerosis; however, substantial numbers of individuals at low or intermediate risk still develop atherosclerosis.

Results

A case-control study was performed using microarray gene expression profiling of peripheral blood from 119 healthy women in the Multi-Ethnic Study of Atherosclerosis cohort aged 50 or above. All participants had low (<10%) to intermediate (10% to 20%) predicted Framingham risk; cases (N = 48) had coronary artery calcium (CAC) score >100 and carotid intima-media thickness (IMT) >1.0 mm, whereas controls (N = 71) had CAC<10 and IMT <0.65 mm. We identified two major expression profiles significantly associated with significant atherosclerosis (odds ratio 4.85; P<0.001); among those with Framingham risk score <10%, the odds ratio was 5.30 (P<0.001). Ontology analysis of the gene signature reveals activation of a major innate immune pathway, toll-like receptors and IL-1R signaling, in individuals with significant atherosclerosis.

Conclusion

Gene expression profiles of peripheral blood may be a useful tool to identify individuals with significant burden of atherosclerosis, even among those with low predicted risk by clinical factors. Furthermore, our data suggest an intimate connection between atherosclerosis and the innate immune system and inflammation via TLR signaling in lower risk individuals.     

Estimation of relevant variables on high-dimensional biological patterns using iterated weighted kernel functions

Background

The analysis of complex proteomic and genomic profiles involves the identification of significant markers within a set of hundreds or even thousands of variables that represent a high-dimensional problem space. The occurrence of noise, redundancy or combinatorial interactions in the profile makes the selection of relevant variables harder.

Methodology/Principal Findings

Here we propose a method to select variables based on estimated relevance to hidden patterns. Our method combines a weighted-kernel discriminant with an iterative stochastic probability estimation algorithm to discover the relevance distribution over the set of variables. We verified the ability of our method to select predefined relevant variables in synthetic proteome-like data and then assessed its performance on biological high-dimensional problems. Experiments were run on serum proteomic datasets of infectious diseases. The resulting variable subsets achieved classification accuracies of 99% on Human African Trypanosomiasis, 91% on Tuberculosis, and 91% on Malaria serum proteomic profiles with fewer than 20% of variables selected. Our method scaled-up to dimensionalities of much higher orders of magnitude as shown with gene expression microarray datasets in which we obtained classification accuracies close to 90% with fewer than 1% of the total number of variables.

Conclusions

Our method consistently found relevant variables attaining high classification accuracies across synthetic and biological datasets. Notably, it yielded very compact subsets compared to the original number of variables, which should simplify downstream biological experimentation.     

Adipose gene expression prior to weight loss can differentiate and weakly predict dietary responders

Background

The ability to identify obese individuals who will successfully lose weight in response to dietary intervention will revolutionize disease management. Therefore, we asked whether it is possible to identify subjects who will lose weight during dietary intervention using only a single gene expression snapshot.

Methodology/Principal Findings

The present study involved 54 female subjects from the Nutrient-Gene Interactions in Human Obesity-Implications for Dietary Guidelines (NUGENOB) trial to determine whether subcutaneous adipose tissue gene expression could be used to predict weight loss prior to the 10-week consumption of a low-fat hypocaloric diet. Using several statistical tests revealed that the gene expression profiles of responders (8–12 kgs weight loss) could always be differentiated from non-responders (<4 kgs weight loss). We also assessed whether this differentiation was sufficient for prediction. Using a bottom-up (i.e. black-box) approach, standard class prediction algorithms were able to predict dietary responders with up to 61.1%±8.1% accuracy. Using a top-down approach (i.e. using differentially expressed genes to build a classifier) improved prediction accuracy to 80.9%±2.2%.

Conclusion

Adipose gene expression profiling prior to the consumption of a low-fat diet is able to differentiate responders from non-responders as well as serve as a weak predictor of subjects destined to lose weight. While the degree of prediction accuracy currently achieved with a gene expression snapshot is perhaps insufficient for clinical use, this work reveals that the comprehensive molecular signature of adipose tissue paves the way for the future of personalized nutrition.     

Gene expression signatures of radiation response are specific, durable and accurate in mice and humans

Background

Previous work has demonstrated the potential for peripheral blood (PB) gene expression profiling for the detection of disease or environmental exposures.

Methods and Findings

We have sought to determine the impact of several variables on the PB gene expression profile of an environmental exposure, ionizing radiation, and to determine the specificity of the PB signature of radiation versus other genotoxic stresses. Neither genotype differences nor the time of PB sampling caused any lessening of the accuracy of PB signatures to predict radiation exposure, but sex difference did influence the accuracy of the prediction of radiation exposure at the lowest level (50 cGy). A PB signature of sepsis was also generated and both the PB signature of radiation and the PB signature of sepsis were found to be 100% specific at distinguishing irradiated from septic animals. We also identified human PB signatures of radiation exposure and chemotherapy treatment which distinguished irradiated patients and chemotherapy-treated individuals within a heterogeneous population with accuracies of 90% and 81%, respectively.

Conclusions

We conclude that PB gene expression profiles can be identified in mice and humans that are accurate in predicting medical conditions, are specific to each condition and remain highly accurate over time.     

Identification of factors contributing to variability in a blood-based gene expression test

Background

Corus CAD is a clinically validated test based on age, sex, and expression levels of 23 genes in whole blood that provides a score (1–40 points) proportional to the likelihood of obstructive coronary disease. Clinical laboratory process variability was examined using whole blood controls across a 24 month period: Intra-batch variability was assessed using sample replicates; inter-batch variability examined as a function of laboratory personnel, equipment, and reagent lots.

Methods/Results

To assess intra-batch variability, five batches of 132 whole blood controls were processed; inter-batch variability was estimated using 895 whole blood control samples. ANOVA was used to examine inter-batch variability at 4 process steps: RNA extraction, cDNA synthesis, cDNA addition to assay plates, and qRT-PCR. Operator, machine, and reagent lots were assessed as variables for all stages if possible, for a total of 11 variables. Intra- and inter-batch variations were estimated to be 0.092 and 0.059 Cp units respectively (SD); total laboratory variation was estimated to be 0.11 Cp units (SD). In a regression model including all 11 laboratory variables, assay plate lot and cDNA kit lot contributed the most to variability (p = 0.045; 0.009 respectively). Overall, reagent lots for RNA extraction, cDNA synthesis, and qRT-PCR contributed the most to inter-batch variance (52.3%), followed by operators and machines (18.9% and 9.2% respectively), leaving 19.6% of the variance unexplained.

Conclusion

Intra-batch variability inherent to the PCR process contributed the most to the overall variability in the study while reagent lot showed the largest contribution to inter-batch variability.     

The burden of selected chronic non-communicable diseases and their risk factors in Malawi: nationwide STEPS survey

Background

Chronic non-communicable diseases (NCDs) are becoming significant causes of morbidity and mortality, particularly in sub-Saharan African countries, although local, high-quality data to inform evidence-based policies are lacking.

Objectives

To determine the magnitude of NCDs and their risk factors in Malawi.

Methods

Using the WHO STEPwise approach to chronic disease risk factor surveillance, a population-based, nationwide cross-sectional survey was conducted between July and September 2009 on participants aged 25–64 years. Socio-demographic and behaviour risk factors were collected in Step 1. Physical anthropometric measurements and blood pressure were documented in Step 2. Blood cholesterol and fasting blood glucose were measured in Step 3.

Results and Conclusion

A total of 5,206 adults (67% females) were surveyed. Tobacco smoking, alcohol drinking and raised blood pressure (BP) were more frequent in males than females, 25% vs 3%, 30% vs 4% and 37% vs 29%. Overweight, physical inactivity and raised cholesterol were more common in females than males, 28% vs 16%, 13% vs 6% and 11% vs 6%. Tobacco smoking was more common in rural than urban areas 11% vs 7%, and overweight and physical inactivity more common in urban than rural areas 39% vs 22% and 24% vs 9%, all with p<0.05. Overall (both sexes) prevalence of tobacco smoking, alcohol consumption, overweight and physical inactivity was 14%, 17%, 22%, 10% and prevalence of raised BP, fasting blood sugar and cholesterol was 33%, 6% and 9% respectively. These data could be useful in the formulation and advocacy of NCD policy and action plan in Malawi.     

Genome-wide DNA methylation indicates silencing of tumor suppressor genes in uterine leiomyoma

Background

Uterine leiomyomas, or fibroids, represent the most common benign tumor of the female reproductive tract. Fibroids become symptomatic in 30% of all women and up to 70% of African American women of reproductive age. Epigenetic dysregulation of individual genes has been demonstrated in leiomyoma cells; however, the in vivo genome-wide distribution of such epigenetic abnormalities remains unknown.

Principal Findings

We characterized and compared genome-wide DNA methylation and mRNA expression profiles in uterine leiomyoma and matched adjacent normal myometrial tissues from 18 African American women. We found 55 genes with differential promoter methylation and concominant differences in mRNA expression in uterine leiomyoma versus normal myometrium. Eighty percent of the identified genes showed an inverse relationship between DNA methylation status and mRNA expression in uterine leiomyoma tissues, and the majority of genes (62%) displayed hypermethylation associated with gene silencing. We selected three genes, the known tumor suppressors KLF11, DLEC1, and KRT19 and verified promoter hypermethylation, mRNA repression and protein expression using bisulfite sequencing, real-time PCR and western blot. Incubation of primary leiomyoma smooth muscle cells with a DNA methyltransferase inhibitor restored KLF11, DLEC1 and KRT19 mRNA levels.

Conclusions

These results suggest a possible functional role of promoter DNA methylation-mediated gene silencing in the pathogenesis of uterine leiomyoma in African American women.     

Circulating micro-RNAs as potential blood-based markers for early stage breast cancer detection

Introduction

MicroRNAs (miRNAs, miRs) are a class of small, non-coding RNA molecules with relevance as regulators of gene expression thereby affecting crucial processes in cancer development. MiRNAs offer great potential as biomarkers for cancer detection due to their remarkable stability in blood and their characteristic expression in many different diseases. We investigated whether microarray-based miRNA profiling on whole blood could discriminate between early stage breast cancer patients and healthy controls.

Methods

We performed microarray-based miRNA profiling on whole blood of 48 early stage breast cancer patients at diagnosis along with 57 healthy individuals as controls. This was followed by a real-time semi-quantitative Polymerase Chain Reaction (RT-qPCR) validation in a separate cohort of 24 early stage breast cancer patients from a breast cancer screening unit and 24 age matched controls using two differentially expressed miRNAs (miR-202, miR-718).

Results

Using the significance level of p<0.05, we found that 59 miRNAs were differentially expressed in whole blood of early stage breast cancer patients compared to healthy controls. 13 significantly up-regulated miRNAs and 46 significantly down-regulated miRNAs in our microarray panel of 1100 miRNAs and miRNA star sequences could be detected. A set of 240 miRNAs that was evaluated by radial basis function kernel support vector machines and 10-fold cross validation yielded a specificity of 78.8%, and a sensitivity of 92.5%, as well as an accuracy of 85.6%. Two miRNAs were validated by RT-qPCR in an independent cohort. The relative fold changes of the RT-qPCR validation were in line with the microarray data for both miRNAs, and statistically significant differences in miRNA-expression were found for miR-202.

Conclusions

MiRNA profiling in whole blood has potential as a novel method for early stage breast cancer detection, but there are still challenges that need to be addressed to establish these new biomarkers in clinical use.     

Association of low level viremia with inflammation and mortality in HIV-infected adults

Background

Whether HIV viremia, particularly at low levels is associated with inflammation, increased coagulation, and all-cause mortality is unclear.

Methods

The associations of HIV RNA level with C-reactive protein (CRP), fibrinogen, interleukin (IL)-6 and mortality were evaluated in 1116 HIV-infected participants from the Study of Fat Redistribution and Metabolic Change in HIV infection. HIV RNA level was categorized as undetectable (i.e., “target not detected”), 1–19, 20–399, 400–9999, and ≥10,000 copies/ml. Covariates included demographics, lifestyle, adipose tissue, and HIV-related factors.

Results

HIV RNA level had little association with CRP. Categories of HIV RNA below 10,000 copies/ml had similar levels of IL-6 compared with an undetectable HIV RNA level, while HIV RNA ≥10,000 copies/ml was associated with 89% higher IL-6 (p<0.001). This association was attenuated by ∼50% after adjustment for CD4+ cell count. Higher HIV RNA was associated with higher fibrinogen. Compared to an undetectable HIV RNA level, fibrinogen was 0.6%, 1.9%, 4.5%, 4.6%, and 9.4% higher across HIV RNA categories, respectively, and statistically significant at the highest level (p = 0.0002 for HIV RNA ≥10,000 copies/ml). Higher HIV RNA was associated with mortality during follow-up in unadjusted analysis, but showed little association after adjustment for CD4+ cell count and inflammation.

Conclusion

HIV RNA ≥10,000 copies/ml was associated with higher IL-6 and fibrinogen, but lower levels of viremia appeared similar, and there was little association with CRP. The relationship of HIV RNA with IL-6 was strongly affected by CD4 cell depletion. After adjustment for CD4+ cell count and inflammation, viremia did not appear to be substantially associated with mortality risk over 5 years.     

Gene Expression Analysis of Peripheral Cells for Subclassification of Pediatric Inflammatory Bowel Disease in Remission

Objective

In current clinical practice, optimal treatment of inflammatory bowel disease (IBD) aims at the induction and maintenance of clinical remission. Clinical remission is apparent when laboratory markers of inflammation are normal and clinical symptoms are absent. However, sub-clinical inflammation can still be present. A detailed analysis of the immune status during this inactive state of disease may provide a useful tool to categorize patients with clinical remission into subsets with variable states of immune activation.

Design

By using Affymetrix GeneChips, we analysed RNA gene expression profiles of peripheral blood leukocytes from pediatric IBD patients in clinical remission and controls. We performed (un)supervised clustering analysis of IBD-associated genes and applied Ingenuity® pathway software to identify specific molecular profiles between patients.

Results

Pediatric IBD patients with disease in clinical remission display heterogeneously distributed gene expression profiles that are significantly distinct from controls. We identified three clusters of IBD patients, each displaying specific expression profiles of IBD-associated genes.

Conclusion

The expression of immune- and IBD-associated genes in peripheral blood leukocytes from pediatric IBD patients in clinical remission was different from healthy controls, indicating that sub-clinical immune mechanisms are still active during remission. As such, RNA profiling of peripheral blood may allow for non-invasive patient subclassification and new perspectives in treatment regimes of IBD patients in the future.     

Hepatitis E virus seroprevalence among blood donors in southwest Switzerland

Aim

The aim of this study was to determine the seroprevalence of Hepatitis E virus (HEV) among blood donors in southwest Switzerland.

Background

HEV is recognized as a food-borne disease in industrialized countries, transmitted mainly through pork meat. Cases of transmission through blood transfusion have also been reported. Recent studies have revealed seroprevalence rates of 13.5%, 16.6% and 20.6% among blood donors in England, France and Denmark, respectively.

Methods

We analyzed 550 consecutive blood donor samples collected in the region of Lausanne, canton of Vaud, Switzerland, for the presence of anti-HEV IgG, using the MP Diagnostics HEV ELISA kit. For each donor, we documented age, sex and alanine aminotransferase (ALT) value.

Results

The study panel was composed of 332 men (60.4%) and 218 women (39.6%). Overall, anti-HEV IgG was found in 27 of 550 samples (4.9%). The seroprevalence was 5.4% (18/332) in men and 4.1% (9/218) in women. The presence of anti-HEV IgG was not correlated with age, gender or ALT values. However, we observed a peak in seroprevalence of 5.3% in individuals aged 51 to 70 years old.

Conclusions

Compared with other European countries, HEV seroprevalence among blood donors in southwest Switzerland is low. The low seroprevalence may be explained by the sensitivity of commercial tests used and/or the strict regulation of animal and meat imports. Data regarding HEV prevalence in Swiss livestock are lacking and merit exploration.     

Family history and breast cancer hormone receptor status in a Spanish cohort

Background

Breast cancer is a heterogenous disease that impacts racial/ethnic groups differently. Differences in genetic composition, lifestyles, reproductive factors, or environmental exposures may contribute to the differential presentation of breast cancer among Hispanic women.

Materials and Methods

A population-based study was conducted in the city of Santiago de Compostela, Spain. A total of 645 women diagnosed with operable invasive breast cancer between 1992 and 2005 participated in the study. Data on demographics, breast cancer risk factors, and clinico-pathological characteristics of the tumors were collected. Hormone receptor negative tumors were compared with hormone receptor postive tumors on their clinico-pathological characteristics as well as risk factor profiles.

Results

Among the 645 breast cancer patients, 78% were estrogen receptor-positive (ER+) or progesterone receptor-positive (PR+), and 22% were ER−&PR−. Women with a family history of breast cancer were more likely to have ER−&PR− tumors than women without a family history (Odds ratio, 1.43; 95% confidence interval, 0.91–2.26). This association was limited to cancers diagnosed before age 50 (Odds ratio, 2.79; 95% confidence interval, 1.34–5.81).

Conclusions

An increased proportion of ER−&PR− breast cancer was observed among younger Spanish women with a family history of the disease.     

Diagnosis of Indian visceral leishmaniasis by nucleic acid detection using PCR

Background

PCR based diagnosis for Visceral Leishmaniasis (VL), despite numerous published primers, remains far from being applied in the field. The present study was planned to design a Leishmania specific diagnostic assay and to evaluate its sensitivity and specificity on a sample size, which to the best of our knowledge is the largest ever screened in one study.

Methods

Leishmania specific primers were developed using 18S rRNA gene and their sensitivity was evaluated on 500 parasitologically confirmed patients with VL and 25 Post Kala-azar Dermal Leishmaniasis (PKDL) patients. Specificity was calculated on 250 healthy endemic controls, 250 healthy non endemic controls and 250 non leishmanial diseases like malaria.

Results

Our PCR assay had a sensitivity of 87.8% (95%CI: 84.1–89.8) using 200 µL of patient''s peripheral-blood. Specificity was absolute in non-endemic healthy controls and in subjects with different diseases while in endemic controls it was 84% (95%CI: 78.9–88.0). Its overall specificity was 94.6% (95%CI-92.8–96.1).

Conclusions

The PCR assay developed is sensitive enough to detect the 18S rRNA gene in an amount equivalent to a single parasite or less in a one million human cell environment. The high sensitivity of this PCR diagnostic test with relatively non-invasive peripheral blood sampling method opens up the possibility of its deployment in field for the routine diagnosis of VL.     

Diagnostic accuracy of PCR in gambiense sleeping sickness diagnosis, staging and post-treatment follow-up: a 2-year longitudinal study

Background

The polymerase chain reaction (PCR) has been proposed for diagnosis, staging and post-treatment follow-up of sleeping sickness but no large-scale clinical evaluations of its diagnostic accuracy have taken place yet.

Methodology/Principal Findings

An 18S ribosomal RNA gene targeting PCR was performed on blood and cerebrospinal fluid (CSF) of 360 T. brucei gambiense sleeping sickness patients and on blood of 129 endemic controls from the Democratic Republic of Congo. Sensitivity and specificity (with 95% confidence intervals) of PCR for diagnosis, disease staging and treatment failure over 2 years follow-up post-treatment were determined. Reference standard tests were trypanosome detection for diagnosis and trypanosome detection and/or increased white blood cell concentration in CSF for staging and detection of treatment failure. PCR on blood showed a sensitivity of 88.4% (84.4–92.5%) and a specificity of 99.2% (97.7–100%) for diagnosis, while for disease staging the sensitivity and specificity of PCR on cerebrospinal fluid were 88.4% (84.8–91.9%) and 82.9% (71.2–94.6%), respectively. During follow-up after treatment, PCR on blood had low sensitivity to detect treatment failure. In cerebrospinal fluid, PCR positivity vanished slowly and was observed until the end of the 2 year follow-up in around 20% of successfully treated patients.

Conclusions/Significance

For T.b. gambiense sleeping sickness diagnosis and staging, PCR performed better than, or similar to, the current parasite detection techniques but it cannot be used for post-treatment follow-up. Continued PCR positivity in one out of five cured patients points to persistence of living or dead parasites or their DNA after successful treatment and may necessitate the revision of some paradigms about the pathophysiology of sleeping sickness.