Influence of sex on basal and dickkopf-1 regulated gene expression in the bovine morula

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudo de associação genômica e perfil de expressão gênica de folículos antrais em novilhas das raças nelore e angus, associadas ao fenótipo de alta e baixa população folicular

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Gene expression profiles displayed by peripheral blood mononuclear cells from patients with type 2 diabetes mellitus focusing on biological processes implicated on the pathogenesis of the disease

Patients with type 2 diabetes mellitus (T2DM) exhibit insulin resistance associated with obesity and inflammatory response, besides an increased level of oxidative DNA damage as a consequence of the hyperglycemic condition and the generation of reactive oxygen species (ROS). In order to provide information on the mechanisms involved in the pathophysiology of T2DM, we analyzed the transcriptional expression patterns exhibited by peripheral blood mononuclear cells (PBMCs) from patients with T2DM compared to non-diabetic subjects, by investigating several biological processes: inflammatory and immune responses, responses to oxidative stress and hypoxia, fatty acid processing, and DNA repair. PBMCs were obtained from 20 T2DM patients and eight non-diabetic subjects. Total RNA was hybridized to Agilent whole human genome 4x44K one-color oligo-microarray. Microarray data were analyzed using the GeneSpring GX 11.0 software (Agilent). We used BRB-ArrayTools software (gene set analysis - GSA) to investigate significant gene sets and the Genomica tool to study a possible influence of clinical features on gene expression profiles. We showed that PBMCs from T2DM patients presented significant changes in gene expression, exhibiting 1320 differentially expressed genes compared to the control group. A great number of genes were involved in biological processes implicated in the pathogenesis of T2DM. Among the genes with high fold-change values, the up-regulated ones were associated with fatty acid metabolism and protection against lipid-induced oxidative stress, while the down-regulated ones were implicated in the suppression of pro-inflammatory cytokines production and DNA repair. Moreover, we identified two significant signaling pathways: adipocytokine, related to insulin resistance; and ceramide, related to oxidative stress and induction of apoptosis. In addition, expression profiles were not influenced by patient features, such as age, gender, obesity, pre/post-menopause age, neuropathy, glycemia, and HbA(1c) percentage. Hence, by studying expression profiles of PBMCs, we provided quantitative and qualitative differences and similarities between T2DM patients and non-diabetic individuals, contributing with new perspectives for a better understanding of the disease. (C) 2012 Elsevier B.V. All rights reserved.

Differential Gene Expression Profiles May Differentiate Responder and Nonresponder Patients with Rheumatoid Arthritis for Methotrexate (MTX) Monotherapy and MTX plus Tumor Necrosis Factor Inhibitor Combined Therapy

Objective. We aimed to evaluate whether the differential gene expression profiles of patients with rheumatoid arthritis (RA) could distinguish responders from nonresponders to methotrexate (MTX) and, in the case of MTX nonresponders, responsiveness to MTX plus anti-tumor necrosis factor-alpha (anti-TNF) combined therapy. Methods. We evaluated 25 patients with RA taking MTX 15-20 mg/week as a monotherapy (8 responders and 17 nonresponders). All MTX nonresponders received intliximab and were reassessed after 20 weeks to evaluate their anti-TNF responsiveness using the European League Against Rheumatism response criteria. A differential gene expression analysis from peripheral blood mononuclear cells was performed in terms of hierarchical gene clustering, and an evaluation of differentially expressed genes was performed using the significance analysis of microarrays program. Results. Hierarchical gene expression clustering discriminated MTX responders from nonresponders, and MTX plus anti-TNF responders from nonresponders. The evaluation of only highly modulated genes (fold change > 1.3 or < 0.7) yielded 5 induced (4 antiapoptotic and CCL4) and 4 repressed (4 proapoptotic) genes in MTX nonresponders compared to responders. In MTX plus anti-TNF nonresponders, the CCL4, CD83, and BCL2A1 genes were induced in relation to responders. Conclusion. Study of the gene expression profiles of RA peripheral blood cells permitted differentiation of responders from nonresponders to MTX and anti-TNF. Several candidate genes in MTX non-responders (CCL4, HTRA2, PRKCD, BCL2A1, CAV1, TNIP1 CASP8AP2, MXD1, and BTG2) and 3 genes in MTX plus anti-TNF nonresponders (CCL4, CD83, and BCL2A1) were identified for further study. (First Release July 1 2012; J Rheumatol 2012;39:1524-32; doi:10.3899/jrheum.120092)

Targeting the p53–MDM2 interaction by the small-molecule MDM2 antagonist Nutlin-3a: a new challenged target therapy in adult Philadelphia positive acute lymphoblastic leukemia patients

The human p53 tumor suppressor, known as the “guardian of the genome”, is one of the most important molecules in human cancers. One mechanism for suppressing p53 uses its negative regulator, MDM2, which modulates p53 by binding directly to and decreasing p53 stability. In testing novel therapeutic approaches activating p53, we investigated the preclinical activity of the MDM2 antagonist, Nutlin-3a, in Philadelphia positive (Ph+) and negative (Ph-) leukemic cell line models, and primary B-Acute lymphoblastic leukemia (ALL) patient samples. In this study we demonstrated that treatment with Nutlin-3a induced grow arrest and apoptosis mediated by p53 pathway in ALL cells with wild-type p53, in time and dose-dependent manner. Consequently, MDM2 inhibitor caused an increase of pro-apoptotic proteins and key regulators of cell cycle arrest. The dose-dependent reduction in cell viability was confirmed in primary blast cells from Ph+ ALL patients with the T315I Bcr-Abl kinase domain mutation. In order to better elucidate the implications of p53 activation and to identify biomarkers of clinical activity, gene expression profiling analysis in sensitive cell lines was performed. A total of 621 genes were differentially expressed (p < 0.05). We found a strong down-regulation of GAS41 (growth-arrest specific 1 gene) and BMI1 (a polycomb ring-finger oncogene) (fold-change -1.35 and -1.11, respectively; p-value 0.02 and 0.03, respectively) after in vitro treatment as compared to control cells. Both genes are repressors of INK4/ARF and p21. Given the importance of BMI in the control of apoptosis, we investigated its pattern in treated and untreated cells, confirming a marked decrease after exposure to MDM2 inhibitor in ALL cells. Noteworthy, the BMI-1 levels remained constant in resistant cells. Therefore, BMI-1 may be used as a biomarker of response. Our findings provide a strong rational for further clinical investigation of Nutlin-3a in Ph+ and Ph-ALL.

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.

DYSREGULATION OF MEOX2 FOLLOWING WT1 MUTATION IN KIDNEY DEVELOPMENT AND WILMS TUMORIGENESIS

Wilms tumor (WT) is a childhood tumor of the kidney and a productive model for understanding the role of genetic alteration and interactions in tumorigenesis. The Wilms tumor gene 1 (WT1) is a transcriptional factor and one of the few genes known to have genetic alterations in WT and has been shown be inactivated in 20% of WTs. However, the mechanisms of how WT1 mutations lead to Wilms tumorigenesis and its influence on downstream genes are unknown. Since it has been established that WT1 is a transcriptional regulator, it has been hypothesized that the loss of WT1 leads to the dysregulation of downstream genes, in turn result in the formation of WTs. To identify the dysregulated downstream genes following WT1 mutations, an Affymetrix GeneChip Human Genome Array was previously conducted to assess the differentially expressed genes in the WT1-wildtype human and WT1-mutant human WTs. Approximately 700 genes were identified as being significantly dysregulated. These genes were further prioritized based on their statistical significance, fold change, chromosomal region, spatial pattern of gene expression and known or putative cellular functions. Mesenchyme homeobox 2 (MEOX2) was one of the most significantly upregulated genes in WT1-mutant WT. MEOX2 is known to play a role in cell proliferation, apoptosis, and differentiation. In addition to its biological roles, it is expressed during early kidney development in the condensed mesenchyme similar to WT1. Furthermore, the use of the Match® web-based tool from the BIOBASE Biological Data base identified a significant predicted WT1 binding site within the first intron of MEOX2. The similarity in spatial gene expression in the developing kidney and the significant predicted WT1 binding site found in the first intron of MEOX2 lead to the development of my hypothesis that MEOX2 is upregulated via a WT1-dependent manner. Here as a part of my master’s work, I have validated the Affymetrix GeneChip Human Genome Array data using an independent set of Wilms tumors. MEOX2 remained upregulated in the mutant WT1 Wilms tumor by 41-fold. Wt1 and Meox2 gene expression were assessed in murine newborn kidney; both Wt1 and Meox2 were expressed in the condensed, undifferentiated metanephric mesenchyme. I have shown that the in vivo ablation of Wt1 during embryonic development at embryonic day (E) 13.5 resulted in the slight increase of Meox2 gene expression by two fold. In order to functionally demonstrate the effect of the loss of Wt1 on Meox2 gene expression in undifferentiated metanephric mesenchyme, I have generated a kidney mesenchymal cell line to genetically ablate Wt1 in vitro by adenoviral infection. The ablation of Wt1 in the kidney mesenchymal cell line resulted in the upregulation of Meox2 by 61-fold. Moreover, the upregulation of Meox2 resulted in the significant induction of p21 and Itgb5. In addition to the dysregulation of these genes the ablation of Wt1 in the kidney mesenchymal cells resulted in decrease in cell growth and loss of cellular adherence. However, it is uncertain whether the upregulation of Meox2 caused this particular cellular phenotype. Overall, I have demonstrated that the upregulation of Meox2 is Wt1-dependent during early kidney development.

Application of stellar photometry to the analysis of microarray images

Improvements in the analysis of microarray images are critical for accurately quantifying gene expression levels. The acquisition of accurate spot intensities directly influences the results and interpretation of statistical analyses. This dissertation discusses the implementation of a novel approach to the analysis of cDNA microarray images. We use a stellar photometric model, the Moffat function, to quantify microarray spots from nylon microarray images. The inherent flexibility of the Moffat shape model makes it ideal for quantifying microarray spots. We apply our novel approach to a Wilms' tumor microarray study and compare our results with a fixed-circle segmentation approach for spot quantification. Our results suggest that different spot feature extraction methods can have an impact on the ability of statistical methods to identify differentially expressed genes. We also used the Moffat function to simulate a series of microarray images under various experimental conditions. These simulations were used to validate the performance of various statistical methods for identifying differentially expressed genes. Our simulation results indicate that tests taking into account the dependency between mean spot intensity and variance estimation, such as the smoothened t-test, can better identify differentially expressed genes, especially when the number of replicates and mean fold change are low. The analysis of the simulations also showed that overall, a rank sum test (Mann-Whitney) performed well at identifying differentially expressed genes. Previous work has suggested the strengths of nonparametric approaches for identifying differentially expressed genes. We also show that multivariate approaches, such as hierarchical and k-means cluster analysis along with principal components analysis, are only effective at classifying samples when replicate numbers and mean fold change are high. Finally, we show how our stellar shape model approach can be extended to the analysis of 2D-gel images by adapting the Moffat function to take into account the elliptical nature of spots in such images. Our results indicate that stellar shape models offer a previously unexplored approach for the quantification of 2D-gel spots. ^

Molecular and biochemical analyses of wounding in mesocarp of peach (Prunus persica L. Batsch) ripe fruits

The physiological and molecular responses of ripe fruit to wounding were evaluated in two peach (Prunus persica) varieties ('Glohaven', GH, melting and 'BigTop', BT, slow melting nectarine) by comparing mesocarp samples from wedges (as in minimal processing) and whole fruit as the control. Slight differences between the two varieties were detected in terms of ethylene production, whereas total phenol and flavonoid concentrations, and PPO and POD enzyme activities showed a general increase in wounded GH but not in BT. This was associated with the better appearance of the BT wedges at the end of the experimental period (72 h). Microarray (genome-wide ?PEACH3.0) analysis revealed that a total number of 2218 genes were differentially expressed (p < 0.01, log2 fold change expression ratio >1 or <-1) in GH 24 h after wounding compared to the control. This number was much lower (1208) in BT. According to the enrichment analysis, cell wall, plasma membrane, response to stress, secondary metabolic processes, oxygen binding were the GO categories over-represented among the GH up-regulated genes, whereas plasma membrane and response to endogenous stimulus were the categories over-represented among the down-regulated genes. Only 32 genes showed a common expression trend in the two varieties 24 h after wounding, whereas a total of 512 genes (with highly represented transcription factors), displayed opposite behavior. Quantitative RT-PCR analysis confirmed the microarray data for 18 out of a total of 20 genes selected. Specific WRKY, AP2/ERF and HSP20 genes were markedly up-regulated in wounded GH, indicating the activation of regulatory and signaling mechanisms probably related to different hormone categories. Compared to BT, the expression of specific genes involved in phenylpropanoid and triterpenoid biosynthetic pathways showed a more pronounced induction in GH, highlighting the difference between the two peach varieties in terms of molecular responses to wounding in the mesocarp tissue.

Quaternary carbonate record of ODP Hole 115-709A

The CaCO3 content in Quaternary deep-sea sediments from Pacific and Atlantic oceans have been suggested to respond differently to glacial/interglacial cycles; CaCO3 contents are highest during glacials in the Pacific but highest during interglacials in the Atlantic Ocean. It is not yet clear as to whether a Pacific or an Atlantic pattern of CaCO3 fluctuations dominates the Indian Ocean. We have analyzed the Ocean Drilling Program (ODP) Site 709A from the western equatorial Indian Ocean for the last 1370 ka to determine the relationships between percentages and fluxes of CaCO3 and Quaternary paleoclimatic changes. We also analyzed the coarse (>25 µm) and fine (<25 µm) fractions of CaCO3 in an attempt at estimating the influence of differences in productivity of foraminifera and calcareous nannofossils in shaping the CaCO3 record. Carbon isotopes and Ba/Al ratios were used as indices of productivity. Percentages and fluxes of CaCO3 in the total sediment and <25 µm fraction do not show any clear relationships to glacial/interglacial cycles derived from d18O of the planktonic foraminifera Globigerinoides ruber. This indicates that CaCO3 fluctuations at this site do not show either a Pacific or an Atlantic pattern of CaCO3 fluctuations. Fluxes of CaCO3 (0.38 to 2.46 g/cm**2/ ka) in total sediment and Ba/Al ratios (0.58 to 3.93 g/cm**2/ka) show six-fold variability through the last 1370 ka, which points out that productivity changes are significant at this site. Fluxes of the fine CaCO3 component demonstrate a 26-fold change (0.02 to 0.52 g/cm**2/ka), whereas the coarse CaCO3 component exhibit eight-fold change (0.13 to 1.07 g/cm**2/ka). This suggests that productivity variations of calcareous nannofossils are greater in comparison with the foraminifera. On the other hand, mean values of coarse CaCO3 fluxes are higher compared to those of fine CaCO3, which reveals that the foraminifera contribute more to the bulk CaCO3 flux than the calcareous nannofossils in the equatorial Indian Ocean.

Major cellular and physiological impacts of ocean acidification on a reef building coral

As atmospheric levels of CO2 increase, reef-building corals are under greater stress from both increased sea surface temperatures and declining sea water pH. To date, most studies have focused on either coral bleaching due to warming oceans or declining calcification due to decreasing oceanic carbonate ion concentrations. Here, through the use of physiology measurements and cDNA microarrays, we show that changes in pH and ocean chemistry consistent with two scenarios put forward by the Intergovernmental Panel on Climate Change (IPCC) drive major changes in gene expression, respiration, photosynthesis and symbiosis of the coral, Acropora millepora, before affects on biomineralisation are apparent at the phenotype level. Under high CO2 conditions corals at the phenotype level lost over half their Symbiodinium populations, and had a decrease in both photosynthesis and respiration. Changes in gene expression were consistent with metabolic suppression, an increase in oxidative stress, apoptosis and symbiont loss. Other expression patterns demonstrate upregulation of membrane transporters, as well as the regulation of genes involved in membrane cytoskeletal interactions and cytoskeletal remodeling. These widespread changes in gene expression emphasize the need to expand future studies of ocean acidification to include a wider spectrum of cellular processes, many of which may occur before impacts on calcification.

Evolution of extreme stomach pH in bilateria inferred from gastric alkalization mechanisms in basal deuterostomes

The stomachs of most vertebrates operate at an acidic pH of 2 generated by the gastric H+/K+-ATPase located in parietal cells. The acidic pH in stomachs of vertebrates is believed to aid digestion and to protect against environmental pathogens. Little attention has been placed on whether acidic gastric pH regulation is a vertebrate character or a deuterostome ancestral trait. Here, we report alkaline conditions up to pH 10.5 in the larval digestive systems of ambulacraria (echinoderm + hemichordate), the closest relative of the chordate. Microelectrode measurements in combination with specific inhibitors for acid-base transporters and ion pumps demonstrated that the gastric alkalization machinery in sea urchin larvae is mainly based on direct H+ secretion from the stomach lumen and involves a conserved set of ion pumps and transporters. Hemichordate larvae additionally utilized HCO3- transport pathways to generate even more alkaline digestive conditions. Molecular analyses in combination with acidification experiments supported these findings and identified genes coding for ion pumps energizing gastric alkalization. Given that insect larval guts were also reported to be alkaline, our discovery raises the hypothesis that the bilaterian ancestor utilized alkaline digestive system while the vertebrate lineage has evolved a strategy to strongly acidify their stomachs.

Whole transcriptome analysis of the coral Acropora millepora reveals complex responses to CO2-driven acidification during the initiation of calcification

The impact of ocean acidification (OA) on coral calcification, a subject of intense current interest, is poorly understood in part because of the presence of symbionts in adult corals. Early life history stages of Acropora spp. provide an opportunity to study the effects of elevated CO(2) on coral calcification without the complication of symbiont metabolism. Therefore, we used the Illumina RNAseq approach to study the effects of acute exposure to elevated CO(2) on gene expression in primary polyps of Acropora millepora, using as reference a novel comprehensive transcriptome assembly developed for this study. Gene ontology analysis of this whole transcriptome data set indicated that CO(2) -driven acidification strongly suppressed metabolism but enhanced extracellular organic matrix synthesis, whereas targeted analyses revealed complex effects on genes implicated in calcification. Unexpectedly, expression of most ion transport proteins was unaffected, while many membrane-associated or secreted carbonic anhydrases were expressed at lower levels. The most dramatic effect of CO(2) -driven acidification, however, was on genes encoding candidate and known components of the skeletal organic matrix that controls CaCO(3) deposition. The skeletal organic matrix effects included elevated expression of adult-type galaxins and some secreted acidic proteins, but down-regulation of other galaxins, secreted acidic proteins, SCRiPs and other coral-specific genes, suggesting specialized roles for the members of these protein families and complex impacts of OA on mineral deposition. This study is the first exhaustive exploration of the transcriptomic response of a scleractinian coral to acidification and provides an unbiased perspective on its effects during the early stages of calcification.

Invariant size–frequency distributions along a latitudinal gradient in marine bivalves

In the most extensive analysis of body size in marine invertebrates to date, we show that the size–frequency distributions of northeastern Pacific bivalves at the provincial level are surprisingly invariant in modal and median size as well as size range, despite a 4-fold change in species richness from the tropics to the Arctic. The modal sizes and shapes of these size–frequency distributions are consistent with the predictions of an energetic model previously applied to terrestrial mammals and birds. However, analyses of the Miocene–Recent history of body sizes within 82 molluscan genera show little support for the expectation that the modal size is an evolutionary attractor over geological time.