Differential expression analysis for sequence count data via mixtures of negative binomials

The recent advent of Next-generation sequencing technologies has revolutionized the way of analyzing the genome. This innovation allows to get deeper information at a lower cost and in less time, and provides data that are discrete measurements. One of the most important applications with these data is the differential analysis, that is investigating if one gene exhibit a different expression level in correspondence of two (or more) biological conditions (such as disease states, treatments received and so on). As for the statistical analysis, the final aim will be statistical testing and for modeling these data the Negative Binomial distribution is considered the most adequate one especially because it allows for "over dispersion". However, the estimation of the dispersion parameter is a very delicate issue because few information are usually available for estimating it. Many strategies have been proposed, but they often result in procedures based on plug-in estimates, and in this thesis we show that this discrepancy between the estimation and the testing framework can lead to uncontrolled first-type errors. We propose a mixture model that allows each gene to share information with other genes that exhibit similar variability. Afterwards, three consistent statistical tests are developed for differential expression analysis. We show that the proposed method improves the sensitivity of detecting differentially expressed genes with respect to the common procedures, since it is the best one in reaching the nominal value for the first-type error, while keeping elevate power. The method is finally illustrated on prostate cancer RNA-seq data.

Differential expression of the bone and the liver tissue non-specific alkaline phosphatase isoforms in brain tissues

The enzyme tissue non-specific alkaline phosphatase (TNAP) belongs to the ectophosphatase family. It is present in large amounts in bone in which it plays a role in mineralization but little is known about its function in other tissues. Arguments are accumulating for its involvement in the brain, in particular in view of the neurological symptoms accompanying human TNAP deficiencies. We have previously shown, by histochemistry, alkaline phosphatase (AP) activity in monkey brain vessels and parenchyma in which AP exhibits specific patterns. Here, we clearly attribute this activity to TNAP expression rather than to other APs in primates (human and marmoset) and in rodents (rat and mouse). We have not found any brain-specific transcripts but our data demonstrate that neuronal and endothelial cells exclusively express the bone TNAP transcript in all species tested, except in mouse neurons in which liver TNAP transcripts have also been detected. Moreover, we highlight the developmental regulation of TNAP expression; this also acts during neuronal differentiation. Our study should help to characterize the regulation of the expression of this ectophosphatase in various cell types of the central nervous system.

Differential expression of IGF-I mRNA and peptide in the male and female gonad during early development of a bony fish, the tilapia Oreochromis niloticus

Insulin-like growth factor I (IGF-I) plays a key role in the complex system that regulates bony fish growth, differentiation, and reproduction. The major source of circulating IGF-I is liver, but IGF-I-producing cells also occur in other organs, including the gonads. Because no data are available on the potential production sites of IGF-I in gonad development, developmental stages of monosex breedings of male and female tilapia from 0 day postfertilization (DPF) to 90 DPF were investigated for the production sites of IGF-I at the peptide (immunohistochemistry) and mRNA (in situ hybridization) level. IGF-I mRNA first appeared in somatic cells of the male and female gonad anlage at 7 DPF followed by IGF-I peptide around 9-10 DPF. Gonad anlagen were detected from 7 DPF. Starting at 7 DPF, IGF-I peptide but no IGF-I mRNA was observed in male and female primordial germ cells (PGCs) provided that IGF-I mRNA was not under the detection level, this observation may suggest that IGF-I originates from the somatic cells and is transferred to the PGCs or is of maternal origin. While in female germ cells IGF-I mRNA and peptide appeared at 29 DPF, in male germ cells both were detected as late as at 51-53 DPF. It is assumed that the production of IGF-I in the germ cells is linked to the onset of meiosis that in tilapia ovary starts at around 28 DPF and in testes at around 52-53 DPF. In adult testis, IGF-I mRNA and peptide occurred in the majority of spermatogonia and spermatocytes as well as in Leydig cells, the latter indicating a role of IGF-I in the synthesis of male sex steroids. In adult ovary, IGF-I mRNA and IGF-I peptide were always present in small and previtellogenic oocytes but only IGF-I peptide infrequently occurred in oocytes at the later stages. IGF-I expression appeared in numerous granulosa and some theca cells of follicles at the lipid stage and persisted in follicles with mature oocytes. The results suggest a crucial role of local IGF-I in the formation, differentiation and function of tilapia gonads.

Differential expression of TGFbeta-stimulated clone 22 in normal prostate and prostate cancer

The transforming growth factor-beta (TGFbeta) superfamily and its downstream effector genes are key regulators of epithelial homeostasis. Altered expression of these genes may be associated with malignant transformation of the prostate gland. The cDNA array analysis of differential expression of the TGFbeta superfamily and functionally related genes between patient-matched noncancerous prostate (NP) and prostate cancer (PC) bulk tissue specimens highlighted two genes, namely TGFbeta-stimulated clone-22 (TSC-22) and Id4. Verification of their mRNA expression by real-time PCR in patient-matched NP and PC bulk tissue, in laser-captured pure epithelial and cancer cells and in NP and PC cell lines confirmed TSC-22 underexpression, but not Id4 overexpression, in PC and in human PC cell lines. Immunohistochemical analysis showed that TSC-22 protein expression in NP is restricted to the basal cells and colocalizes with the basal cell marker cytokeratin 5. In contrast, all matched PC samples lack TSC-22 immunoreactivity. Likewise, PC cell lines do not show detectable TSC-22 protein expression as shown by immunoblotting. TSC-22 should be considered as a novel basal cell marker, potentially useful for studying lineage determination within the epithelial compartment of the prostate. Conversely, lack of TSC-22 seems to be a hallmark of malignant transformation of the prostate epithelium. Accordingly, TSC-22 immunohistochemistry may prove to be a diagnostic tool for discriminating benign lesions from malignant ones of the prostate. The suggested tumour suppressor function of TSC-22 warrants further investigation on its role in prostate carcinogenesis and on the TSC-22 pathway as a candidate therapeutic target in PC.

Differential expression of the receptor tyrosine kinase EphB4 and its ligand Ephrin-B2 during human placental development

Normal placentation involves the development of an utero-placental circulation following the migration of the extravillous cytotrophoblasts into the decidua and invasion of the spiral arteries, which are thereby transformed into large vessels of low resistance. Given the documented role of the receptor tyrosine kinase EphB4 and its ligand ephrin-B2 in the establishment of the embryonal vascular network, we hypothesized that these molecules are also instrumental in the development of the human placenta. Monitoring the expression during placental development revealed that in first trimester and term placentae both molecules are equally expressed at the RNA level. In contrast, the protein levels were significantly reduced during gestation. Immunohistochemistry revealed a distinct localization of the EphB4 and ephrin-B2 proteins. EphB4 was predominantly expressed in the villous syncytial trophoblast layer and in a subset of intravillous capillaries. Prominent expression was also observed in the extravillous cytotrophoblast giant cells. In contrast, ephrin-B2 expression was detected in the villous cytotrophoblast and syncytial trophoblast cell layers, as well as initially in all intravillous capillaries. Strong expression was also observed in extravillous anchoring cytotrophoblast cells. Hypoxia is a major inducer of placental development. In vitro studies employing trophoblast-derived cell lines revealed that predominantly ephrin-B2 expression is induced by hypoxia, however, in an Hif-1alpha independent manner. These experiments suggest that EphB4 and ephrin-B2 are instrumental in the establishment of a functional placental structure and of the utero-placental circulation.

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.

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.

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