Arrow Plot: a new graphical tool for selecting up and down regulated genes and genes differentially expressed on samples subgroups

Background: A common task in analyzing microarray data is to determine which genes are differentially expressed across two (or more) kind of tissue samples or samples submitted under experimental conditions. Several statistical methods have been proposed to accomplish this goal, generally based on measures of distance between classes. It is well known that biological samples are heterogeneous because of factors such as molecular subtypes or genetic background that are often unknown to the experimenter. For instance, in experiments which involve molecular classification of tumors it is important to identify significant subtypes of cancer. Bimodal or multimodal distributions often reflect the presence of subsamples mixtures. Consequently, there can be genes differentially expressed on sample subgroups which are missed if usual statistical approaches are used. In this paper we propose a new graphical tool which not only identifies genes with up and down regulations, but also genes with differential expression in different subclasses, that are usually missed if current statistical methods are used. This tool is based on two measures of distance between samples, namely the overlapping coefficient (OVL) between two densities and the area under the receiver operating characteristic (ROC) curve. The methodology proposed here was implemented in the open-source R software. Results: This method was applied to a publicly available dataset, as well as to a simulated dataset. We compared our results with the ones obtained using some of the standard methods for detecting differentially expressed genes, namely Welch t-statistic, fold change (FC), rank products (RP), average difference (AD), weighted average difference (WAD), moderated t-statistic (modT), intensity-based moderated t-statistic (ibmT), significance analysis of microarrays (samT) and area under the ROC curve (AUC). On both datasets all differentially expressed genes with bimodal or multimodal distributions were not selected by all standard selection procedures. We also compared our results with (i) area between ROC curve and rising area (ABCR) and (ii) the test for not proper ROC curves (TNRC). We found our methodology more comprehensive, because it detects both bimodal and multimodal distributions and different variances can be considered on both samples. Another advantage of our method is that we can analyze graphically the behavior of different kinds of differentially expressed genes. Conclusion: Our results indicate that the arrow plot represents a new flexible and useful tool for the analysis of gene expression profiles from microarrays.

Decreased STAMP2 expression in association with visceral adipose tissue dysfunction.

CONTEXT Six-transmembrane protein of prostate 2 (STAMP2) is a counter-regulator of inflammation and insulin resistance according to findings in mice. However, there have been contradictory reports in humans. OBJECTIVE We aimed to explore STAMP2 in association with inflammatory and metabolic status of human obesity. DESIGN, PATIENTS, AND METHODS STAMP2 gene expression was analyzed in adipose tissue samples (171 visceral and 67 sc depots) and during human preadipocyte differentiation. Human adipocytes were treated with macrophage-conditioned medium, TNF-α, and rosiglitazone. RESULTS In visceral adipose tissue, STAMP2 gene expression was significantly decreased in obese subjects, mainly in obese subjects with type 2 diabetes. STAMP2 gene expression and protein were significantly and inversely associated with obesity phenotype measures (body mass index, waist, hip, and fat mass) and obesity-associated metabolic disturbances (systolic blood pressure and fasting glucose). In addition, STAMP2 gene expression was positively associated with lipogenic (FASN, ACC1, SREBP1, THRSP14, TRα, and TRα1), CAV1, IRS1, GLUT4, and CD206 gene expression. In sc adipose tissue, STAMP2 gene expression was not associated with metabolic parameters. In both fat depots, STAMP2 gene expression in stromovascular cells was significantly higher than in mature adipocytes. STAMP2 gene expression was significantly increased during the differentiation process in parallel to adipogenic genes, being increased in preadipocytes derived from lean subjects. Macrophage-conditioned medium (25%) and TNF-α (100 ng/ml) administration increased whereas rosiglitazone (2 μM) decreased significantly STAMP2 gene expression in human differentiated adipocytes. CONCLUSIONS Decreased STAMP2 expression (mRNA and protein) might reflect visceral adipose dysfunction in subjects with obesity and type 2 diabetes.

Breast cancer 1 (BrCa1) may be behind decreased lipogenesis in adipose tissue from obese subjects.

CONTEXT Expression and activity of the main lipogenic enzymes is paradoxically decreased in obesity, but the mechanisms behind these findings are poorly known. Breast Cancer 1 (BrCa1) interacts with acetyl-CoA carboxylase (ACC) reducing the rate of fatty acid biosynthesis. In this study, we aimed to evaluate BrCa1 in human adipose tissue according to obesity and insulin resistance, and in vitro cultured adipocytes. RESEARCH DESIGN AND METHODS BrCa1 gene expression, total and phosphorylated (P-) BrCa1, and ACC were analyzed in adipose tissue samples obtained from a total sample of 133 subjects. BrCa1 expression was also evaluated during in vitro differentiation of human adipocytes and 3T3-L1 cells. RESULTS BrCa1 gene expression was significantly up-regulated in both omental (OM; 1.36-fold, p = 0.002) and subcutaneous (SC; 1.49-fold, p = 0.001) adipose tissue from obese subjects. In parallel with increased BrCa1 mRNA, P-ACC was also up-regulated in SC (p = 0.007) as well as in OM (p = 0.010) fat from obese subjects. Consistent with its role limiting fatty acid biosynthesis, both BrCa1 mRNA (3.5-fold, p<0.0001) and protein (1.2-fold, p = 0.001) were increased in pre-adipocytes, and decreased during in vitro adipogenesis, while P-ACC decreased during differentiation of human adipocytes (p = 0.005) allowing lipid biosynthesis. Interestingly, BrCa1 gene expression in mature adipocytes was restored by inflammatory stimuli (macrophage conditioned medium), whereas lipogenic genes significantly decreased. CONCLUSIONS The specular findings of BrCa1 and lipogenic enzymes in adipose tissue and adipocytes reported here suggest that BrCa1 might help to control fatty acid biosynthesis in adipocytes and adipose tissue from obese subjects.

ESR1 gene promoter region methylation in free circulating DNA and its correlation with estrogen receptor protein expression in tumor tissue in breast cancer patients.

BACKGROUND Tumor expression of estrogen receptor (ER) is an important marker of prognosis, and is predictive of response to endocrine therapy in breast cancer. Several studies have observed that epigenetic events, such methylation of cytosines and deacetylation of histones, are involved in the complex mechanisms that regulate promoter transcription. However, the exact interplay of these factors in transcription activity is not well understood. In this study, we explored the relationship between ER expression status in tumor tissue samples and the methylation of the 5' CpG promoter region of the estrogen receptor gene (ESR1) isolated from free circulating DNA (fcDNA) in plasma samples from breast cancer patients. METHODS Patients (n = 110) with non-metastatic breast cancer had analyses performed of ER expression (luminal phenotype in tumor tissue, by immunohistochemistry method), and the ESR1-DNA methylation status (fcDNA in plasma, by quantitative methylation specific PCR technique). RESULTS Our results showed a significant association between presence of methylated ESR1 in patients with breast cancer and ER negative status in the tumor tissue (p = 0.0179). There was a trend towards a higher probability of ESR1-methylation in those phenotypes with poor prognosis i.e. 80% of triple negative patients, 60% of HER2 patients, compared to 28% and 5.9% of patients with better prognosis such as luminal A and luminal B, respectively. CONCLUSION Silencing, by methylation, of the promoter region of the ESR1 affects the expression of the estrogen receptor protein in tumors of breast cancer patients; high methylation of ESR1-DNA is associated with estrogen receptor negative status which, in turn, may be implicated in the patient's resistance to hormonal treatment in breast cancer. As such, epigenetic markers in plasma may be of interest as new targets for anticancer therapy, especially with respect to endocrine treatment.

Tissue biomarker development in a multicentre trial context: a feasibility study on the PETACC3 stage II and III colon cancer adjuvant treatment trial.

PURPOSE: We evaluated the feasibility of biomarker development in the context of multicenter clinical trials. EXPERIMENTAL DESIGN: Formalin-fixed, paraffin-embedded (FFPE) tissue samples were collected from a prospective adjuvant colon cancer trial (PETACC3). DNA was isolated from tumor as well as normal tissue and used for analysis of microsatellite instability, KRAS and BRAF genotyping, UGT1A1 genotyping, and loss of heterozygosity of 18 q loci. Immunohistochemistry was used to test expression of TERT, SMAD4, p53, and TYMS. Messenger RNA was retrieved and tested for use in expression profiling experiments. RESULTS: Of the 3,278 patients entered in the study, FFPE blocks were obtained from 1,564 patients coming from 368 different centers in 31 countries. In over 95% of the samples, genomic DNA tests yielded a reliable result. Of the immmunohistochemical tests, p53 and SMAD4 staining did best with reliable results in over 85% of the cases. TERT was the most problematic test with 46% of failures, mostly due to insufficient tissue processing quality. Good quality mRNA was obtained, usable in expression profiling experiments. CONCLUSIONS: Prospective clinical trials can be used as framework for biomarker development using routinely processed FFPE tissues. Our results support the notion that as a rule, translational studies based on FFPE should be included in prospective clinical trials.

Insulin degradation by adipose tissue is increased in human obesity

White adipose tissue samples from obese and lean patients were used for the estimation ofinsulin protease and insulin:glutathione transhydrogenase using 1251-labeled insulin. There was no activity detected in the absence of reduced glutathione, which indicates that insulin is cleaved in human adipose "tissue through reduction of the disulfide bridge between the chains. O bese patients showed higher transhydrogenase activity (per U tissue protein wt, per U tissue wt, and in the total adipose tissue mass) than the lean group. There is a significant correlation between the activity per U tissue wt, and protein and total activity in the whole adipose tissue with respect to body mass index, with a higher activity in obese patients. The potential ofinsulin cleavage by adipose tissue in obese patients was a mean 5.6-fold higher than that in controla. The coexistence of high insulinemia and high cleavage capability implies that insulin secretion and turnover are increased in the o bese. Thus, white adipose tissue may be crucial in the control of energy availability through modulation ofinsulin cleavage.

Integrative molecular bioinformatics study of human adrenocortical tumors : microRNA, tissue-specific target prediction, and pathway analysis.

MicroRNAs (miRs) are involved in the pathogenesis of several neoplasms; however, there are no data on their expression patterns and possible roles in adrenocortical tumors. Our objective was to study adrenocortical tumors by an integrative bioinformatics analysis involving miR and transcriptomics profiling, pathway analysis, and a novel, tissue-specific miR target prediction approach. Thirty-six tissue samples including normal adrenocortical tissues, benign adenomas, and adrenocortical carcinomas (ACC) were studied by simultaneous miR and mRNA profiling. A novel data-processing software was used to identify all predicted miR-mRNA interactions retrieved from PicTar, TargetScan, and miRBase. Tissue-specific target prediction was achieved by filtering out mRNAs with undetectable expression and searching for mRNA targets with inverse expression alterations as their regulatory miRs. Target sets and significant microarray data were subjected to Ingenuity Pathway Analysis. Six miRs with significantly different expression were found. miR-184 and miR-503 showed significantly higher, whereas miR-511 and miR-214 showed significantly lower expression in ACCs than in other groups. Expression of miR-210 was significantly lower in cortisol-secreting adenomas than in ACCs. By calculating the difference between dCT(miR-511) and dCT(miR-503) (delta cycle threshold), ACCs could be distinguished from benign adenomas with high sensitivity and specificity. Pathway analysis revealed the possible involvement of G2/M checkpoint damage in ACC pathogenesis. To our knowledge, this is the first report describing miR expression patterns and pathway analysis in sporadic adrenocortical tumors. miR biomarkers may be helpful for the diagnosis of adrenocortical malignancy. This tissue-specific target prediction approach may be used in other tumors too.

Apparent versus true gene expression changes of three hypoxia-related genes in autopsy derived tissue and the importance of normalisation.

The aim of our work was to show how a chosen normal-isation strategy can affect the outcome of quantitative gene expression studies. As an example, we analysed the expression of three genes known to be upregulated under hypoxic conditions: HIF1A, VEGF and SLC2A1 (GLUT1). Raw RT-qPCR data were normalised using two different strategies: a straightforward normalisation against a single reference gene, GAPDH, using the 2(-ΔΔCt) algorithm and a more complex normalisation against a normalisation factor calculated from the quantitative raw data from four previously validated reference genes. We found that the two different normalisation strategies revealed contradicting results: normalising against a validated set of reference genes revealed an upregulation of the three genes of interest in three post-mortem tissue samples (cardiac muscle, skeletal muscle and brain) under hypoxic conditions. Interestingly, we found a statistically significant difference in the relative transcript abundance of VEGF in cardiac muscle between donors who died of asphyxia versus donors who died from cardiac death. Normalisation against GAPDH alone revealed no upregulation but, in some instances, a downregulation of the genes of interest. To further analyse this discrepancy, the stability of all reference genes used were reassessed and the very low expression stability of GAPDH was found to originate from the co-regulation of this gene under hypoxic conditions. We concluded that GAPDH is not a suitable reference gene for the quantitative analysis of gene expression in hypoxia and that validation of reference genes is a crucial step for generating biologically meaningful data.

Nested-PCR for the detection of Mycoplasma hyopneumoniae in bronchial alveolar swabs, frozen tissues and formalin-fixed paraffin-embedded swine lung samples: comparative evaluation with immunohistochemical findings and histological features

The diagnosis of Mycoplasma hyopneumoniae infection is often performed through histopathology, immunohistochemistry (IHC) and polymerase chain reaction (PCR) or a combination of these techniques. PCR can be performed on samples using several conservation methods, including swabs, frozen tissue or formalin-fixed and paraffin-embedded (FFPE) tissue. However, the formalin fixation process often inhibits DNA amplification. To evaluate whether M. hyopneumoniae DNA could be recovered from FFPE tissues, 15 lungs with cranioventral consolidation lesions were collected in a slaughterhouse from swine bred in herds with respiratory disease. Bronchial swabs and fresh lung tissue were collected, and a fragment of the corresponding lung section was placed in neutral buffered formalin for 48 hours. A PCR assay was performed to compare FFPE tissue samples with samples that were only refrigerated (bronchial swabs) or frozen (tissue pieces). M. hyopneumoniae was detected by PCR in all 15 samples of the swab and frozen tissue, while it was detected in only 11 of the 15 FFPE samples. Histological features of M. hyopneumoniae infection were presented in 11 cases and 7 of these samples stained positive in IHC. Concordance between the histological features and detection results was observed in 13 of the FFPE tissue samples. PCR was the most sensitive technique. Comparison of different sample conservation methods indicated that it is possible to detect M. hyopneumoniae from FFPE tissue. It is important to conduct further research using archived material because the efficiency of PCR could be compromised under these conditions.

Interphase cytogenetics using fluorescence in situ hybridization: an overview of its application to diffuse and solid tissue

Interphase cytogenetics, utilizing fluorescence in situ hybridization (FISH) techniques, has been successfully applied to diffuse and solid tissue specimens. Most studies have been performed on isolated cells, such as blood or bone marrow cells; a few have been performed on cells from body fluids, such as amniotic fluid, urine, sperm, and sputum. Mechanically or chemically disaggregated cells from solid tissues have also been used as single cell suspensions for FISH. Additionally, intact organized tissue samples represented by touch preparations or thin tissue sections have been used, especially in cancer studies. Advantages and pitfalls of application of FISH methodology to each type of specimen and some significant biological findings achieved are illustrated in this overview.

Standardization of a fluorimetric assay for the determination of tissue angiotensin-converting enzyme activity in rats

The tripeptide Hip-His-Leu was used to standardize a fluorimetric method to measure tissue angiotensin-converting enzyme (ACE) activity in rats. The fluorescence of the o-phthaldialdehyde-His-Leu adduct was compared in the presence and absence of the homogenate (25 µl) to determine whether the homogenate from different tissues interfered with the fluorimetric determination of the His-Leu product. Only homogenates from lung and renal medulla and cortex showed significantly altered fluorescence intensity. To overcome this problem, the homogenate from these tissues were diluted 10 times with assay buffer. The specificity of the assay was demonstrated by the inhibition of ACE activity with 3 µM enalaprilat (MK-422). There was a linear relationship between product formation and incubation time for up to 90 min for homogenates of renal cortex and medulla and liver, for up to 60 min for ventricles and adrenals and for up to 30 min for the aorta, lung and atrium homogenates. In addition, there was a linear relationship between product formation and the amount of protein in the homogenates within the following range: lung, 30-600 µg; renal cortex and medulla, 40-400 µg; atrium and ventricles, 20-200 µg; adrenal, 20-100 µg; aorta, 5-100 µg; liver, 5-25 µg. No peptidase activity against the His-Leu product (31 nmol), assayed in borate buffer (BB), was detected in the different homogenates except the liver homogenate, which was inhibited by 0.1 mM r-chloromercuribenzoic acid. ACE activity in BB was higher than in phosphate buffer (PB) due, at least in part, to a greater hydrolysis of the His-Leu product in PB. ACE activity of lung increased 20% when BB plus Triton was used. Enzyme activity was stable when the homogenates were stored at -20o or -70oC for at least 30 days. These results indicate a condition whereby ACE activity can be easily and efficiently assayed in rat tissue samples homogenized in BB using a fluorimetric method with Hip-His-Leu as a substrate.

Female pelvic floor disorders – outcome after mesh-augmented procedures and studies on vaginal connective tissue

Pelvic floor disorders, such as urinary incontinence (UI) and pelvic organ prolapse (POP), are common disorders in women. Because of the prolonged life expectancy the prevalence of UI and POP and the probability of ending up in surgery are increasing. However, the pathophysiology behind these disorders is still unsolved. The aim of this thesis is to study possible alterations in the connective tissue in the vaginal wall in patients with and without POP. The long-term outcome and complications of mid-urethral slings (MUS) and mesh-augmented POP surgery were studied in heterogenic patient populations. More elastin and a slight increase in immunostaining of type III and V collagens in tissue samples were obtained from patients with POP compared to controls in whom type I collagen was more prominent. The studies assessing the mesh-augmented procedures revealed good efficacy and high patient satisfaction after a long-term follow-up. Patients operated on because of mixed incontinence and with BMI >30 kg/m² reported significantly more urinary symptoms and a lower quality of life than the patients operated on because of stress urinary incontinence and the ones with BMI ≤30 kg/m². The objective outcome was equal between the groups. Mesh exposure through vaginal mucosa occurred in 23 % of the patients after POP surgery, most of these being asymptomatic. There are alterations in connective tissues in patients with POP. Mid-urethral sling procedures produced good long-term cure rates and patient satisfaction. As to the prolapse surgery, in spite of relatively high exposure rate, mesh-augmented procedure proved to be safe and effective method for the correction of POP.

Influence of birth weight on gene regulators of lipid metabolism and utilization in subcutaneous adipose tissue and skeletal muscle of neonatal pigs.

Epidemiological studies suggest that low-birth weight infants show poor neonatal growth and increased susceptibility to metabolic syndrome, in particular, obesity and diabetes. Adipose tissue development is regulated by many genes, including members of the peroxisome proliferator-activated receptor (PPAR) and the fatty acid-binding protein (FABP) families. The aim of this study was to determine the influence of birth weight on key adipose and skeletal muscle tissue regulating genes. Piglets from 11 litters were ranked according to birth weight and 3 from each litter assigned to small, normal, or large-birth weight groups. Tissue samples were collected on day 7 or 14. Plasma metabolite concentrations and the expression of PPARG2, PPARA, FABP3, and FABP4 genes were determined in subcutaneous adipose tissue and skeletal muscle. Adipocyte number and area were determined histologically. Expression of FABP3 and 4 was significantly reduced in small and large, compared with normal, piglets in adipose tissue on day 7 and in skeletal muscle on day 14. On day 7, PPARA and PPARG2 were significantly reduced in adipose tissue from small and large piglets. Adipose tissue from small piglets contained more adipocytes than normal or large piglets. Birth weight had no effect on adipose tissue and skeletal muscle lipid content. Low-birth weight is associated with tissue-specific and time-dependent effects on lipid-regulating genes as well as morphological changes in adipose tissue. It remains to be seen whether these developmental changes alter an individual's susceptibility to metabolic syndrome.

The quantitation of lipoprotein lipase mRNA in biopsies of human adipose tissue, using the polymerase chain reaction, and the effect of increased consumption of n-3 polyunsaturated fatty acids

Objective: To examine the effects of the consumption of fish oils on the gene expression of lipoprotein lipase (LPL, EC 3.1.1.34) in human adipose tissue. In order to measure LPL mRNA in adipose tissue samples obtained by needle biopsy from human volunteers a competitive, reverse transcriptase PCR (RT-PCR) protocol was developed. Design: A randomised controlled, single blind cross over dietary study which compared the effects of a low level n-3 polyunsaturated fatty acids (PUFA) using normal foods enriched with eicosapentaenoic (EPA) and docosahexaenoic (DHA) (test diet), with non-enriched but otherwise identical foods (control). The diets were consumed for a period of 22 d with a wash out period of 5 months between the diets. Setting: Free-living individuals associated with the University of Surrey. Subjects: Six male subjects with a mean (±sd) age of 51.2±3.6 y were recruited. Major Outcome Measures: Pre-and postprandial blood samples were taken for the measurement of triacylglycerol (TAG), postheparin LPL activity and adipose tissue samples for the measurement of LPL mRNA levels. Results: Mean LPL expression values were 4.12´105 molecules of LPL mRNA per ng total RNA on the control diet and 4.60´105 molecules of LPL mRNA per ng total RNA on the n-3 PUFA enriched (test) diet. There was no significant difference between the levels of LPL expression following each diet, consistent with the lack of change in TAG levels in response to increased dietary n-3 PUFA intake. However, the change in LPL expression (Test-Control diet) correlated significantly with the change in fasting TAG levels (P=0.03, R=-0.87 and R2=0.75) and with the total area under the TAG-time response curve (P=0.003, R=-0.96 and R2=0.92) in individuals. Conclusions: These findings, although based on a small number of subjects, suggest that LPL expression may be a determinant of plasma TAG levels. The development of this methodology should allow further elucidation of the effects of dietary manipulation and disease processes on lipid clearance and regulation in human subjects.

Measurement of paraben concentrations in human breast tissue at serial locations across the breast from axilla to sternum

The concentrations of five esters of p-hydroxybenzoic acid (parabens) were measured using HPLC-MS/MS at four serial locations across the human breast from axilla to sternum using human breast tissue collected from 40 mastectomies for primary breast cancer in England between 2005 and 2008. One or more paraben esters were quantifiable in 158/160 (99%) of the tissue samples and in 96/160 (60%) all five esters were measured. Variation was notable with respect to individual paraben esters, location within one breast and similar locations in different breasts. Overall median values in nanograms per gram tissue for the 160 tissue samples were highest for n-propylparaben [16.8 (range 0–2052.7)] and methylparaben [16.6 (range 0–5102.9)]; levels were lower for n-butylparaben [5.8 (range 0–95.4)], ethylparaben [3.4 (range 0–499.7)] and isobutylparaben 2.1 (range 0–802.9). The overall median value for total paraben was 85.5 ng g−1 tissue (range 0–5134.5). The source of the paraben cannot be identified, but paraben was measured in the 7/40 patients who reported never having used underarm cosmetics in their lifetime. No correlations were found between paraben concentrations and age of patient (37–91 years), length of breast feeding (0–23 months), tumour location or tumour oestrogen receptor content. In view of the disproportionate incidence of breast cancer in the upper outer quadrant, paraben concentrations were compared across the four regions of the breast: n-propylparaben was found at significantly higher levels in the axilla than mid (P = 0.004 Wilcoxon matched pairs) or medial (P = 0.021 Wilcoxon matched pairs) regions (P = 0.010 Friedman ANOVA).