Differential Expression with the Bioconductor Project

A basic, yet challenging task in the analysis of microarray gene expression data is the identification of changes in gene expression that are associated with particular biological conditions. We discuss different approaches to this task and illustrate how they can be applied using software from the Bioconductor Project. A central problem is the high dimensionality of gene expression space, which prohibits a comprehensive statistical analysis without focusing on particular aspects of the joint distribution of the genes expression levels. Possible strategies are to do univariate gene-by-gene analysis, and to perform data-driven nonspecific filtering of genes before the actual statistical analysis. However, more focused strategies that make use of biologically relevant knowledge are more likely to increase our understanding of the data.

Differential expression of CD45 on canine microglial cells

CD45, also called leucocyte common antigen is a transmembrane protein tyrosine phosphatase on the surface of nearly all white blood cells and has a functional role in signal transduction. In the brain, the expression of CD45 can be used to distinguish microglial cells with a characteristic phenotype of CD11b/c+ and CD45(low) from other central nervous system (CNS) macrophages which show an expression of CD11b/c+ and CD45(high). In the course of pathological changes in the CNS, microglia in rodents is known to readily upregulate expression of various surface molecules, such as CD45. Understanding the mechanisms that regulate expression of surface molecules is essential to study the pathogenesis of CNS diseases. In the present study, the expression of CD45 on microglia of 42 dogs was examined ex vivo by means of flow cytometry. The dogs were classified in two groups according to the histopathological diagnosis in the CNS. All dogs without changes in the CNS (group I; n = 22) only showed low percentages of CD45+ microglial cells. In group II consisting of 20 dogs with different intracranial diseases varying results were obtained. Thirteen dogs showed a low percentage of CD45+ microglial cells whereas seven dogs exhibited high percentages of microglial cells expressing CD45. Evaluation of expression intensity in these seven dogs revealed two subpopulations of CD45+ microglial cells: a large subpopulation with CD45(low) and a small subpopulation with CD45(high). The expression intensity of CD45(high) was comparable with that of canine monocytes. It was attempted to correlate these findings to age of the animals, underlying disease, duration of clinical signs, medical treatment, occurrence of seizure activity and the expression of other surface molecules. It appeared that dogs with high percentages of CD45+ suffered from long-lasting CNS disease with seizures. In future studies, the reason and consequences for upregulated CD45 in long-lasting CNS diseases has to be further evaluated.

Vitamin C deficiency in weanling guinea pigs: differential expression of oxidative stress and DNA repair in liver and brain

Neonates are particularly susceptible to malnutrition due to their limited reserves of micronutrients and their rapid growth. In the present study, we examined the effect of vitamin C deficiency on markers of oxidative stress in plasma, liver and brain of weanling guinea pigs. Vitamin C deficiency caused rapid and significant depletion of ascorbate (P < 0.001), tocopherols (P < 0.001) and glutathione (P < 0.001), and a decrease in superoxide dismutase activity (P = 0.005) in the liver, while protein oxidation was significantly increased (P = 0.011). No changes in lipid oxidation or oxidatively damaged DNA were observed in this tissue. In the brain, the pattern was markedly different. Of the measured antioxidants, only ascorbate was significantly depleted (P < 0.001), but in contrast to the liver, ascorbate oxidation (P = 0.034), lipid oxidation (P < 0.001), DNA oxidation (P = 0.13) and DNA incision repair (P = 0.014) were all increased, while protein oxidation decreased (P = 0.003). The results show that the selective preservation of brain ascorbate and induction of DNA repair in vitamin C-deficient weanling guinea pigs is not sufficient to prevent oxidative damage. Vitamin C deficiency may therefore be particularly adverse during the neonatal period.

Differential expression of stx2 variants in Shiga toxin-producing Escherichia coli belonging to seropathotypes A and C

Only a subset of Shiga toxin (Stx)-producing Escherichia coli (STEC) are human pathogens, but the characteristics that account for differences in pathogenicity are not well understood. In this study, we investigated the distribution of the stx variants coding for Stx2 and its variants in highly virulent STEC of seropathotype A and low-pathogenic STEC of seropathotype C. We analysed and compared transcription of the corresponding genes, production of Shiga toxins, and stx-phage release in basal as well as in induced conditions. We found that the stx(2) variant was mainly associated with strains of seropathotype A, whereas most of the strains of seropathotype C possessed the stx(2-vhb) variant, which was frequently associated with stx(2), stx(2-vha) or stx(2c). Levels of stx(2) and stx(2)-related mRNA were higher in strains belonging to seropathotype A and in those strains of seropathotype C that express the stx(2) variant than in the remaining strains of seropathotype C. The stx(2-vhb) genes were the least expressed, in basal as well as in induced conditions, and in many cases did not seem to be carried by an inducible prophage. A clear correlation was observed between stx mRNA levels and stx-phage DNA in the culture supernatants, suggesting that most stx(2)-related genes are expressed only when they are carried by a phage. In conclusion, some relationship between stx(2)-related gene expression in vitro and the seropathotype of the STEC strains was observed. A higher expression of the stx(2) gene and a higher release of its product, in basal as well as in induced conditions, was observed in pathogenic strains of seropathotype A. A subset of strains of seropathotype C shows the same characteristics and could be a high risk to human health.

Differential expression of ABC transporters and their regulatory genes during lactation and dry period in bovine mammary tissue

ATP-binding cassette (ABC) transporters play a pivotal role in human physiology, and mutations in these genes often result in severe hereditary diseases. ABC transporters are expressed in the bovine mammary gland but their physiological role in this organ remains elusive. Based on findings in the context of human disorders we speculated that candidate ABC transporters are implicated in lipid and cholesterol transport in the mammary gland. Therefore we investigated the expression pattern of selected genes that are associated with sterol transport in lactating and nonlactating mammary glands of dairy cows. mRNA levels from mammary gland biopsies taken during lactation and in the first and second week of the dry period were analysed using quantitative PCR. Five ABC transporter genes, namely ABCA1, ABCA7, ABCG1, ABCG2 and ABCG5, their regulating genes LXRalpha, PPARgamma, SREBP1 and the milk proteins lactoferrin and alpha-lactalbumin were assessed. A significantly enhanced expression in the dry period was observed for ABCA1 while a significant decrease of expression in this period was detected for ABCA7, ABCG2, SREBP1 and alpha-lactalbumin. ABCG1, ABCG5, LXRalpha, PPARgamma and lactoferrin expression was not altered between lactation and dry period. These results indicate that candidate ABC transporters involved in lipid and cholesterol transport show differential mRNA expression between lactation and the dry period. This may be due to physiological changes in the mammary gland such as immigration of macrophages or the accumulation of fat due to the loss of liquid in the involuting mammary gland. The current mRNA expression analysis of transporters in the mammary gland is the prerequisite for elucidating novel molecular mechanisms underlying cholesterol and lipid transfer into milk.

Metalloprotease meprin beta in rat kidney: glomerular localization and differential expression in glomerulonephritis

Meprin (EC 3.4.24.18) is an oligomeric metalloendopeptidase found in microvillar membranes of kidney proximal tubular epithelial cells. Here, we present the first report on the expression of meprin beta in rat glomerular epithelial cells and suggest a potential involvement in experimental glomerular disease. We detected meprin beta in glomeruli of immunostained rat kidney sections on the protein level and by quantitative RT-PCR of laser-capture microdissected glomeruli on the mRNA level. Using immuno-gold staining we identified the membrane of podocyte foot processes as the main site of meprin beta expression. The glomerular meprin beta expression pattern was altered in anti-Thy 1.1 and passive Heymann nephritis (PHN). In addition, the meprin beta staining pattern in the latter was reminiscent of immunostaining with the sheep anti-Fx1A antiserum, commonly used in PHN induction. Using Western blot and immunoprecipitation assays we demonstrated that meprin beta is recognized by Fx1A antiserum and may therefore represent an auto-antigen in PHN. In anti-Thy 1.1 glomerulonephritis we observed a striking redistribution of meprin beta in tubular epithelial cells from the apical to the basolateral side and the cytosol. This might point to an involvement of meprin beta in this form of glomerulonephritis.

Expression profiling with RNA from formalin-fixed, paraffin-embeeded material

Background Molecular characterization of breast and other cancers by gene expression profiling has corroborated existing classifications and revealed novel subtypes. Most profiling studies are based on fresh frozen (FF) tumor material which is available only for a limited number of samples while thousands of tumor samples exist as formalin-fixed, paraffin-embedded (FFPE) blocks. Unfortunately, RNA derived of FFPE material is fragmented and chemically modified impairing expression measurements by standard procedures. Robust protocols for isolation of RNA from FFPE material suitable for stable and reproducible measurement of gene expression (e.g. by quantitative reverse transcriptase PCR, QPCR) remain a major challenge. Results We present a simple procedure for RNA isolation from FFPE material of diagnostic samples. The RNA is suitable for expression measurement by QPCR when used in combination with an optimized cDNA synthesis protocol and TaqMan assays specific for short amplicons. The FFPE derived RNA was compared to intact RNA isolated from the same tumors. Preliminary scores were computed from genes related to the ER response, HER2 signaling and proliferation. Correlation coefficients between intact and partially fragmented RNA from FFPE material were 0.83 to 0.97. Conclusion We developed a simple and robust method for isolating RNA from FFPE material. The RNA can be used for gene expression profiling. Expression measurements from several genes can be combined to robust scores representing the hormonal or the proliferation status of the tumor.

Differential expression and activity of 11beta-hydroxysteroid dehydrogenase in human placenta and fetal membranes from pregnancies with intrauterine growth restriction

OBJECTIVES: To study the expression and the function of the 11beta-hydroxysteroid dehydrogenase enzyme 1 (11beta-HSD1) and 2 (11beta-HSD2) in placenta and the fetal membranes from pregnancies with intrauterine growth restriction (IUGR) and from controls. METHODS: Amnion, chorion, decidua and cotyledon were separated from placenta; mRNA was analyzed by TaqMan real-time technology and proteins by Western blot; enzyme activities were measured by the conversion of 3H-cortisol to 3H-cortisone and vice versa. RESULTS: Predominant mRNA expression (p < 0.001) was found for 11beta-HSD1 in chorion and for 11beta-HSD2 in decidua and cotyledon. In pregnancies with IUGR, 11beta-HSD1 was upregulated in chorion (mean DeltaCt 11beta-HSD:18S mRNA 193.5 vs. 103.0 in controls respectively, p < 0.05) and 11beta-HSD2 was downregulated in decidua (mean DeltaCt 11beta-HSD2:18S mRNA 0.18 vs. 15.88 in controls respectively, p < 0.05). 11beta-HSD1 protein levels were reduced in amnion and 11beta-HSD1 and 11beta-HSD2 oxidase activity in decidua and cotyledon were reduced from pregnancies with IUGR. CONCLUSION: Reduced synthesis or activity of 11beta-HSD1 or 2 in cases of IUGR is shown in some but not in all tissues. The local mRNA expression of 11beta-HSD1 in chorion may reflect a mechanism on the post-transcriptional gene regulation to stimulate the formation of cortisone in IUGR. To provoke increasing activity with oxidase stimulators could be a future therapy in cases of IUGR.

Differential expression and localization of lipid transporters in the bovine mammary gland during the pregnancy-lactation cycle

The transport of lipids across mammary gland epithelial cells (MEC) determines milk lipid content and composition. We investigated the expression of lipid transporters and their regulators in comparison to blood metabolites during lactation and dry period (DP) in dairy cows. Repeated mammary gland biopsies and blood samples were taken from 10 animals at 7 stages of the pregnancy-lactation cycle. Expression levels of the specific mRNAs were determined by quantitative reverse transcription-PCR, whereas ABCA1 was localized by immunohistochemistry. Blood serum metabolites were determined by common enzymatic chemistries. Elevated mRNA profiles of ABCA1 and ABCA7 were found during DP as compared with lactation and were inversely associated with blood cholesterol levels. Elevated levels of ABCG2, NPC1, SREBP1, SREBP2, LXR alpha, and PPAR gamma were found postpartum, whereas ABCG1 did not differ between the functional stages of the mammary gland. The ABCA1 protein was localized in MEC and showed differential activity between DP and lactation suggesting a role of ABCA1 in the removal of excess cellular cholesterol from MEC during the DP. The expression profiles of ABCA7 and NPC1 may reflect a role of these transporters in the clearance of apoptotic cells and the intracellular redistribution of cholesterol, respectively. Regulation of lipid transporters in the mammary gland is partially associated with transcription factors that control lipid homeostasis.

Differential expression of chemokines in normal pancreas and in chronic pancreatitis

Cellular infiltrates are present already in early stages of chronic pancreatitis. The mechanisms responsible for their recruitment are unknown. Hence, we determined the differential expression of chemokine genes and their cellular sources in normal and affected pancreatic tissues.

Gene expression profiling in anti-CD11d mAb-treated spinal cord-injured rats

Acute administration of a mononclonal antibody (mAb) raised against the CD11d subunit of the leukocyte CD11d/CD18 integrin after spinal cord injury (SCI) in the rat greatly improves neurological outcomes. We have profiled gene expression in anti-CD11d and isotyped-matched control mAb-treated rats after SCI. Microarray analysis demonstrated reduced expression of pro-inflammatory cytokines and increased expression of inflammatory mediators that promote wound healing and the expression of scar proteins predicted to improve nerve growth. These changes in gene expression may reflect changes in the types of macrophages that populate the lesions in anti-CD11d mAb-treated rats.

Comprehensive MicroRNA expression profiling identifies novel markers in follicular variant of papillary thyroid carcinoma

BACKGROUND Follicular variant of papillary thyroid carcinoma (FVPTC) shares features of papillary (PTC) and follicular (FTC) thyroid carcinomas on a clinical, morphological, and genetic level. MicroRNA (miRNA) deregulation was extensively studied in PTCs and FTCs. However, very limited information is available for FVPTC. The aim of this study was to assess miRNA expression in FVPTC with the most comprehensive miRNA array panel and to correlate it with the clinicopathological data. METHODS Forty-four papillary thyroid carcinomas (17 FVPTC, 27 classic PTC) and eight normal thyroid tissue samples were analyzed for expression of 748 miRNAs using Human Microarray Assays on the ABI 7900 platform (Life Technologies, Carlsbad, CA). In addition, an independent set of 61 tumor and normal samples was studied for expression of novel miRNA markers detected in this study. RESULTS Overall, the miRNA expression profile demonstrated similar trends between FVPTC and classic PTC. Fourteen miRNAs were deregulated in FVPTC with a fold change of more than five (up/down), including miRNAs known to be upregulated in PTC (miR-146b-3p, -146-5p, -221, -222 and miR-222-5p) and novel miRNAs (miR-375, -551b, 181-2-3p, 99b-3p). However, the levels of miRNA expression were different between these tumor types and some miRNAs were uniquely dysregulated in FVPTC allowing separation of these tumors on the unsupervised hierarchical clustering analysis. Upregulation of novel miR-375 was confirmed in a large independent set of follicular cell derived neoplasms and benign nodules and demonstrated specific upregulation for PTC. Two miRNAs (miR-181a-2-3p, miR-99b-3p) were associated with an adverse outcome in FVPTC patients by a Kaplan-Meier (p < 0.05) and multivariate Cox regression analysis (p < 0.05). CONCLUSIONS Despite high similarity in miRNA expression between FVPTC and classic PTC, several miRNAs were uniquely expressed in each tumor type, supporting their histopathologic differences. Highly upregulated miRNA identified in this study (miR-375) can serve as a novel marker of papillary thyroid carcinoma, and miR-181a-2-3p and miR-99b-3p can predict relapse-free survival in patients with FVPTC thus potentially providing important diagnostic and predictive value.

Differential expression of urate oxidase in rat liver

An affinity-purified monospecific antibody was prepared to study the differential expression of the peroxisomal enzyme urate oxidase in rat liver during development and in various metabolic states. Monospecific antibody for urate oxidase was affinity purified from a pool of antibodies initially produced against a mixture of proteins from a Percoll density gradient fraction. Immunogold staining of samples of the gradient fraction and rat liver tissue with the affinity-purified antibody demonstrated labelling of peroxisomal core structures. Screening of liver homogenates from rats at different developmental stages using immunoblot analysis demonstrated low levels of urate oxidase prior to 20 days of age; at 20 days of age, urate oxidase levels are 2-fold greater than the 15-day old levels and approximate adult levels. Catalase expression during rat development mimicked the differential expression pattern of urate oxidase. The increase between days 15 and 20 was determined to be independent of the process of weaning. Administration of exogenous glucocorticoid hormone to 10-day old rats resulted in a precocious rise (2.5-fold) in urate oxidase levels, but adrenalectomy at 10 days of age did not cause decreased expression in the fourth week of life. In adult animals, exogenous glucocorticoid did not influence urate oxidase levels, but adrenalectomized rats had urate oxidase levels that were 40 percent of control expression 4 days post-surgery. Catalase expression was not influenced by glucocorticoid status in these studies. Glucocorticoid regulation of urate oxidase expression appears to be one part of a more complex mechanism controlling levels of the enzyme. Exogenous glucocorticoid administration influenced urate oxidase levels in an age-dependent manner; in addition, it is possible that the control mechanism for urate oxidase may include factors which can modulate expression in the absence of glucocorticoids. The effect of glucocorticoids on urate oxidase expression can not be extended to include all peroxisomal proteins, since catalase is unaffected. Glucocorticoids appear to participate in the complex regulation of urate oxidase expression; glucocorticoids influence urate oxidase specifically and do not modulate all peroxisomal proteins. ^

The differential expression of the $\beta\sb2$-adrenergic receptor modifies agonist-dependent adenylylcyclase activation

There have been multiple reports which indicate that variations in $\beta$AR expression affect the V$\sb{\rm max}$ observed for the agonist-dependent activation of adenylylcyclase. This observation has been ignored by most researchers when V$\sb{\rm max}$ values obtained for wild type and mutant receptors are compared. Such an imprecise analysis may lead to erroneous conclusions concerning the ability of a receptor to activate adenylylcyclase. Equations were derived from the Cassel-Selinger model of GTPase activity and Tolkovsky and Levitzki's Collision Coupling model which predict that the EC$\sb{50}$ and V$\sb{\rm max}$ for the activation of adenylylcyclase are a function of receptor number. Experimental results for L cell clones in which either hamster or human $\beta$AR were transfected at varying levels showed that EC$\sb{50}$ decreases and V$\sb{\rm max}$ increases as receptor number increases. Comparison of these results with simulations obtained from the equations describing EC$\sb{50}$ and V$\sb{\rm max}$ showed a close correlation. This documents that the kinetic parameters of adenylylcyclase activation change with the level of receptor expression and relates this phenomenon to a theoretical framework concerning the mechanisms involved in $\beta$AR signal transduction.^ One of the terms used in the equations which expressed the EC$\sb{50}$ and V$\sb{\rm max}$ as a function of receptor number is coupling efficiency, defined as $\rm k\sb1/k\sb{-1}$. Calculation of $\rm k\sb1/k\sb{-1}$ can be accomplished for wild type receptors with the easily measured experimental values of agonist K$\sb{\rm d}$, EC$\sb{50}$ and receptor number. This was demonstrated for hamster $\beta$AR which yielded a coupling efficiency of 0.15 $\pm$ 0.003 and human $\beta$AR which yielded a coupling efficiency of 0.90 $\pm$ 0.031. $\rm k\sb1/k\sb{-1}$ replaces the traditional qualitative evaluation of the ability to activate adenylylcyclase, which utilizes V$\sb{\rm max}$ without correction for variation in receptor number, with a quantitative definition that more accurately describes the ability of $\beta$AR to couple to G$\sb{\rm s}$.^ The equations which express the EC$\sb{50}$ and V$\sb{\rm max}$ for adenylylcyclase activation as a function of receptor number and coupling efficiency were tested to determine whether they could accurately simulate the changes seen in these parameters during desensitization. Data from original desensitization experiments and data from the literature (24,25,52,54,83) were compared to simulated changes in EC$\sb{50}$ and V$\sb{\rm max}$. In a variety of systems the predictions of the equations were consistent with the changes observed in EC$\sb{50}$ and V$\sb{\rm max}$. In addition reductions in the calculated value of $\rm k\sb1/k\sb{-1}$ was shown to correlate well with $\beta$AR phosphorylation and to be minimally affected by sequestration and down-regulation. ^

The AtProT family. Compatible solute transporters with similar substrate specificity but differential expression patterns

Proline transporters (ProTs) mediate transport of the compatible solutes Pro, glycine betaine, and the stress-induced compound gamma-aminobutyric acid. A new member of this gene family, AtProT3, was isolated from Arabidopsis (Arabidopsis thaliana), and its properties were compared to AtProT1 and AtProT2. Transient expression of fusions of AtProT and the green fluorescent protein in tobacco (Nicotiana tabacum) protoplasts revealed that all three AtProTs were localized at the plasma membrane. Expression in a yeast (Saccharomyces cerevisiae) mutant demonstrated that the affinity of all three AtProTs was highest for glycine betaine (K-m = 0.1-0.3 mM), lower for Pro (K-m = 0.4-1 mM), and lowest for gamma-aminobutyric acid (K-m = 4-5 mM). Relative quantification of the mRNA level using real-time PCR and analyses of transgenic plants expressing the beta-glucuronidase (uidA) gene under control of individual AtProT promoters showed that the expression pattern of AtProTs are complementary. AtProT1 expression was found in the phloem or phloem parenchyma cells throughout the whole plant, indicative of a role in long-distance transport of compatible solutes. beta-Glucuronidase activity under the control of the AtProT2 promoter was restricted to the epidermis and the cortex cells in roots, whereas in leaves, staining could be demonstrated only after wounding. In contrast, AtProT3 expression was restricted to the above-ground parts of the plant and could be localized to the epidermal cells in leaves. These results showed that, although intracellular localization, substrate specificity, and affinity are very similar, the transporters fulfill different roles in planta.