Analysis of ripening-related gene expression in papaya using an Arabidopsis-based microarray

Abstract Background Papaya (Carica papaya L.) is a commercially important crop that produces climacteric fruits with a soft and sweet pulp that contain a wide range of health promoting phytochemicals. Despite its importance, little is known about transcriptional modifications during papaya fruit ripening and their control. In this study we report the analysis of ripe papaya transcriptome by using a cross-species (XSpecies) microarray technique based on the phylogenetic proximity between papaya and Arabidopsis thaliana. Results Papaya transcriptome analyses resulted in the identification of 414 ripening-related genes with some having their expression validated by qPCR. The transcription profile was compared with that from ripening tomato and grape. There were many similarities between papaya and tomato especially with respect to the expression of genes encoding proteins involved in primary metabolism, regulation of transcription, biotic and abiotic stress and cell wall metabolism. XSpecies microarray data indicated that transcription factors (TFs) of the MADS-box, NAC and AP2/ERF gene families were involved in the control of papaya ripening and revealed that cell wall-related gene expression in papaya had similarities to the expression profiles seen in Arabidopsis during hypocotyl development. Conclusion The cross-species array experiment identified a ripening-related set of genes in papaya allowing the comparison of transcription control between papaya and other fruit bearing taxa during the ripening process.

Deregulation of apoptosis-related genes is associated with PRV1 overexpression and JAK2 V617F allele burden in Essential Thrombocythemia and Myelofibrosis

Background: Essential Thrombocythemia (ET) and Primary Myelofibrosis (PMF) are Chronic Myeloproliferative Neoplasms (MPN) characterized by clonal myeloproliferation/myeloaccumulation without cell maturation impairment. The JAK2 V617F mutation and PRV1 gene overexpression may contribute to MPN physiopathology. We hypothesized that deregulation of the apoptotic machinery may also play a role in the pathogenesis of ET and PMF. In this study we evaluated the apoptosis-related gene and protein expression of BCL2 family members in bone marrow CD34(+) hematopoietic stem cells (HSC) and peripheral blood leukocytes from ET and PMF patients. We also tested whether the gene expression results were correlated with JAK2 V617F allele burden percentage, PRV1 overexpression, and clinical and laboratory parameters. Results: By real time PCR assay, we observed that A1, MCL1, BIK and BID, as well as A1, BCLW and BAK gene expression were increased in ET and PMF CD34(+) cells respectively, while pro-apoptotic BAX and anti-apoptotic BCL2 mRNA levels were found to be lower in ET and PMF CD34(+) cells respectively, in relation to controls. In patients' leukocytes, we detected an upregulation of anti-apoptotic genes A1, BCL2, BCL-XL and BCLW. In contrast, pro-apoptotic BID and BIMEL expression were downregulated in ET leukocytes. Increased BCL-XL protein expression in PMF leukocytes and decreased BID protein expression in ET leukocytes were observed by Western Blot. In ET leukocytes, we found a correlation between JAK2 V617F allele burden and BAX, BIK and BAD gene expression and between A1, BAX and BIK and PRV1 gene expression. A negative correlation between PRV1 gene expression and platelet count was observed, as well as a positive correlation between PRV1 gene expression and splenomegaly. Conclusions: Our results suggest the participation of intrinsic apoptosis pathway in the MPN physiopathology. In addition, PRV1 and JAK2 V617F allele burden were linked to deregulation of the apoptotic machinery.

Expression profile of apoptosis-related genes in childhood adrenocortical tumors: low level of expression of BCL2 and TNF genes suggests a poor prognosis

Background: Impaired apoptosis has been implicated in the development of childhood adrenocortical tumors (ACT), although the expression of apoptosis-related gene expression in such tumors has not been reported. Methods: The mRNA expression levels of the genes CASP3, CASP8, CASP9, FAS, TNF, NFKB, and BCL2 were analyzed by quantitative real-time PCR in consecutive tumor samples obtained at diagnosis from 60 children with a diagnosis of ACT and in 11 non-neoplastic adrenal samples. BCL2 and TNF protein expression was analyzed by immunohistochemistry. Results: A significant association was observed between tumor size >= 100 g and lower expression levels of the BCL2 (P=0.03) and TNF (P=0.05) genes; between stage IV and lower expression levels of CASP3 (P=0.008), CASP9 (P=0.02), BCL2 (P=0.002), TNF (P=0.05), and NFKB (P=0.03); Weiss score >= 3 and lower expression of TNF (P=0.01); unfavorable event and higher expression values of CASP9 (P=0.01) and lower values of TNF (P=0.02); and death and lower expression of BCL2 (P=0.04). Underexpression of TNF was associated with lower event-free survival in uni- and multivariate analyses (P<0.01). Similar results were observed when patients with Weiss score <3 were excluded. Conclusion: This study supports the participation of apoptosis-related genes in the biology and prognosis of childhood ACT and suggests the complex role of these genes in the pathogenesis of this tumor.

Elevated tissue omega-3 fatty acid status prevents age-related glucose intolerance in fat-1 transgenic mice

The objective of this study was to investigate the impact of elevated tissue omega-3 (n-3) polyunsaturated fatty acids (PUFA) status on age-related glucose intolerance utilizing the fat-1 transgenic mouse model, which can endogenously synthesize n-3 PUFA from omega-6 (n-6) PUFA. Fat-1 and wild-type mice, maintained on the same dietary regime of a 10% corn oil diet, were tested at two different ages (2months old and 8months old) for various glucose homeostasis parameters and related gene expression. The older wild-type mice exhibited significantly increased levels of blood insulin, fasting blood glucose, liver triglycerides, and glucose intolerance, compared to the younger mice, indicating an age-related impairment of glucose homeostasis. In contrast, these age-related changes in glucose metabolism were largely prevented in the older fat-1 mice. Compared to the older wild-type mice, the older fat-1 mice also displayed a lower capacity for gluconeogenesis, as measured by pyruvate tolerance testing (PTT) and hepatic gene expression of phosphoenolpyruvate carboxykinase (PEPCK) and glucose 6 phosphatase (G6Pase). Furthermore, the older fat-1 mice showed a significant decrease in body weight, epididymal fat mass, inflammatory activity (NFκ-B and p-IκB expression), and hepatic lipogenesis (acetyl-CoA carboxylase (ACC) and fatty acid synthase (FAS) expression), as well as increased peroxisomal activity (70-kDa peroxisomal membrane protein (PMP70) and acyl-CoA oxidase1 (ACOX1) expression). Altogether, the older fat-1 mice exhibit improved glucose homeostasis in comparison to the older wild-type mice. These findings support the beneficial effects of elevated tissue n-3 fatty acid status in the prevention and treatment of age-related chronic metabolic diseases

Up-regulation of fas and fasL pro-apoptotic genes expression in type 1 diabetes patients after autologous haematopoietic stem cell transplantation

Type 1 diabetes (T1D) is a chronic autoimmune disease characterized by T cell-mediated destruction of pancreatic beta cells, resulting in insulin deficiency and hyperglycaemia. Recent studies have described that apoptosis impairment during central and peripheral tolerance is involved in T1D pathogenesis. In this study, the apoptosis-related gene expression in T1D patients was evaluated before and after treatment with high-dose immunosuppression followed by autologous haematopoietic stem cell transplantation (HDI-AHSCT). We also correlated gene expression results with clinical response to HDI-AHSCT. We observed a decreased expression of bad, bax and fasL pro-apoptotic genes and an increased expression of a1, bcl-xL and cIAP-2 anti-apoptotic genes in patients' peripheral blood mononuclear cells (PBMCs) compared to controls. After HDI-AHSCT, we found an up-regulation of fas and fasL and a down-regulation of anti-apoptotic bcl-xL genes expression in post-HDI-AHSCT periods compared to pre-transplantation. Additionally, the levels of bad, bax, bok, fasL, bcl-xL and cIAP-1 genes expression were found similar to controls 2 years after HDI-AHSCT. Furthermore, over-expression of pro-apoptotic noxa at 540 days post-HDI-AHSCT correlated positively with insulin-free patients and conversely with glutamic acid decarboxylase autoantibodies (GAD65) autoantibody levels. Taken together, the results suggest that apoptosis-related genes deregulation in patients' PBMCs might be involved in breakdown of immune tolerance and consequently contribute to T1D pathogenesis. Furthermore, HDI-AHSCT modulated the expression of some apoptotic genes towards the levels similar to controls. Possibly, the expression of these apoptotic molecules could be applied as biomarkers of clinical remission of T1D patients treated with HDI-AHSCT therapy.

Prognostic significance of NDRG1 expression in oral and oropharyngeal squamous cell carcinoma

Human N-myc downstream-regulated gene 1 (NDRG1) is a metastasis suppressor gene with several potential functions, including cell differentiation, cell cycle regulation and response to hormones, nickel and stress. The purpose of this study was to investigate the immunoexpression of NDRG1 in oral and oropharyngeal squamous cell carcinomas searching for its role in the clinical course of these tumors. We investigated immunohistochemical expression of NDRG1 protein in 412 tissue microarray cores of tumor samples from 103 patients with oral and oropharyngeal squamous cell carcinomas and in 110 paraffin-embedded surgical margin sections. The results showed NDRG1 up-regulation in 101/103 (98.1 %) tumor samples, but no expression in any normal tissue sample. Western blot assays confirmed the immunohistochemical findings, suggesting that lower levels of NDRG1 are associated with a high mortality rate. NDRG1 overexpression was related to long-term specific survival (HR = 0.38; p = 0.009), whereas the presence of lymph-node metastasis showed the opposite association with survival (HR = 2.45; p = 0.013). Our findings reinforce the idea that NDRG1 plays a metastasis suppressor role in oral and oropharyngeal squamous cell carcinomas and may be a useful marker for these tumors.

Liposuction induces a compensatory increase of visceral fat which is effectively counteracted by physical activity: a randomized trial

Context: Liposuction is suggested to result in long-term body fat regain that could lead to increased cardiometabolic risk. We hypothesized that physical activity could prevent this effect. Objective: Our objective was to investigate the effects of liposuction on body fat distribution and cardiometabolic risk factors in women who were either exercise trained or not after surgery. Design, Setting, and Participants: Thirty-six healthy normal-weight women participated in this 6-month randomized controlled trial at the University of Sao Paulo, Sao Paulo, Brazil. Interventions: Patients underwent a small-volume abdominal liposuction. Two months after surgery, the subjects were randomly allocated into two groups: trained (TR, n = 18, 4-month exercise program) and nontrained (NT, n = 18). Main Outcome Measures: Body fat distribution (assessed by computed tomography) was assessed before the intervention (PRE) and 2 months (POST2), and 6 months (POST6) after surgery. Secondary outcome measures included body composition, metabolic parameters and dietary intake, assessed at PRE, POST2, and POST6, and total energy expenditure, physical capacity, and sc adipocyte size and lipid metabolism-related gene expression, assessed at PRE and POST6. Results: Liposuction was effective in reducing sc abdominal fat (PRE vs. POST2, P = 0.0001). Despite the sustained sc abdominal fat decrement at POST6 (P = 0.0001), the NT group showed a significant 10% increase in visceral fat from PRE to POST6 (P = 0.04; effect size = -0.72) and decreased energy expenditure (P = 0.01; effect size = 0.95) when compared with TR. Dietary intake, adipocyte size, and gene expression were unchanged over time. Conclusion: Abdominal liposuction does not induce regrowth of fat, but it does trigger a compensatory increase of visceral fat, which is effectively counteracted by physical activity. (J Clin Endocrinol Metab 97: 2388-2395, 2012)

Genome analysis of DNA repair genes in the alpha proteobacterium Caulobacter crescentus

Abstract Background The integrity of DNA molecules is fundamental for maintaining life. The DNA repair proteins protect organisms against genetic damage, by removal of DNA lesions or helping to tolerate them. DNA repair genes are best known from the gamma-proteobacterium Escherichia coli, which is the most understood bacterial model. However, genome sequencing raises questions regarding uniformity and ubiquity of these DNA repair genes and pathways, reinforcing the need for identifying genes and proteins, which may respond to DNA damage in other bacteria. Results In this study, we employed a bioinformatic approach, to analyse and describe the open reading frames potentially related to DNA repair from the genome of the alpha-proteobacterium Caulobacter crescentus. This was performed by comparison with known DNA repair related genes found in public databases. As expected, although C. crescentus and E. coli bacteria belong to separate phylogenetic groups, many of their DNA repair genes are very similar. However, some important DNA repair genes are absent in the C. crescentus genome and other interesting functionally related gene duplications are present, which do not occur in E. coli. These include DNA ligases, exonuclease III (xthA), endonuclease III (nth), O6-methylguanine-DNA methyltransferase (ada gene), photolyase-like genes, and uracil-DNA-glycosylases. On the other hand, the genes imuA and imuB, which are involved in DNA damage induced mutagenesis, have recently been described in C. crescentus, but are absent in E. coli. Particularly interesting are the potential atypical phylogeny of one of the photolyase genes in alpha-proteobacteria, indicating an origin by horizontal transfer, and the duplication of the Ada orthologs, which have diverse structural configurations, including one that is still unique for C. crescentus. Conclusion The absence and the presence of certain genes are discussed and predictions are made considering the particular aspects of the C. crescentus among other known DNA repair pathways. The observed differences enlarge what is known for DNA repair in the Bacterial world, and provide a useful framework for further experimental studies in this organism.

Phenotypic alterations of neuropeptide Y and calcitonin gene-related peptide-containing neurons innervating the rat temporomandibular joint during carrageenan-induced arthritis

The aim of this study was to identify immunoreactive neuropeptide Y (NPY) and calcitonin gene-related peptide (CGRP) neurons in the autonomic and sensory ganglia, specifically neurons that innervate the rat temporomandibular joint (TMJ). A possible variation between the percentages of these neurons in acute and chronic phases of carrageenan-induced arthritis was examined. Retrograde neuronal tracing was combined with indirect immunofluorescence to identify NPY-immuno-reactive (NPY-IR) and CGRP-immunoreactive (CGRP-IR) neurons that send nerve fibers to the normal and arthritic temporomandibular joint. In normal joints, NPY-IR neurons constitute 78 +/- 3%, 77 +/- 6% and 10 +/- 4% of double-labeled nucleated neuronal profile originated from the superior cervical, stellate and otic ganglia, respectively. These percentages in the sympathetic ganglia were significantly decreased in acute (58 +/- 2% for superior cervical ganglion and 58 +/- 8% for stellate ganglion) and chronic (60 +/- 2% for superior cervical ganglion and 59 +/- 15% for stellate ganglion) phases of arthritis, while in the otic ganglion these percentages were significantly increased to 19 +/- 5% and 13 +/- 3%, respectively. In the trigeminal ganglion, CGRP-IR neurons innervating the joint significantly increased from 31 +/- 3% in normal animals to 54 +/- 2% and 49 +/- 3% in the acute and chronic phases of arthritis, respectively. It can be concluded that NPY neurons that send nerve fibers to the rat temporomandibular joint are located mainly in the superior cervical, stellate and otic ganglia. Acute and chronic phases of carrageenan-induced arthritis lead to an increase in the percentage of NPY-IR parasympathetic and CGRP-IR sensory neurons and to a decrease in the percentage of NPY-IR sympathetic neurons related to TMJ innervation.

The clinical impact of chromosomal rearrangements with breakpoints upstream of the SOX9 gene: two novel de novo balanced translocations associated with acampomelic campomelic dysplasia

Abstract Background The association of balanced rearrangements with breakpoints near SOX9 [SRY (sex determining region Y)-box 9] with skeletal abnormalities has been ascribed to the presumptive altering of SOX9 expression by the direct disruption of regulatory elements, their separation from SOX9 or the effect of juxtaposed sequences. Case presentation We report on two sporadic apparently balanced translocations, t(7;17)(p13;q24) and t(17;20)(q24.3;q11.2), whose carriers have skeletal abnormalities that led to the diagnosis of acampomelic campomelic dysplasia (ACD; MIM 114290). No pathogenic chromosomal imbalances were detected by a-CGH. The chromosome 17 breakpoints were mapped, respectively, 917–855 kb and 601–585 kb upstream of the SOX9 gene. A distal cluster of balanced rearrangements breakpoints on chromosome 17 associated with SOX9-related skeletal disorders has been mapped to a segment 932–789 kb upstream of SOX9. In this cluster, the breakpoint of the herein described t(17;20) is the most telomeric to SOX9, thus allowing the redefining of the telomeric boundary of the distal breakpoint cluster region related to skeletal disorders to 601–585 kb upstream of SOX9. Although both patients have skeletal abnormalities, the t(7;17) carrier presents with relatively mild clinical features, whereas the t(17;20) was detected in a boy with severe broncheomalacia, depending on mechanical ventilation. Balanced and unbalanced rearrangements associated with disorders of sex determination led to the mapping of a regulatory region of SOX9 function on testicular differentiation to a 517–595 kb interval upstream of SOX9, in addition to TESCO (Testis-specific enhancer of SOX9 core). As the carrier of t(17;20) has an XY sex-chromosome constitution and normal male development for his age, the segment of chromosome 17 distal to the translocation breakpoint should contain the regulatory elements for normal testis development. Conclusions These two novel translocations illustrate the clinical variability in carriers of balanced translocations with breakpoints near SOX9. The translocation t(17;20) breakpoint provides further evidence for an additional testis-specific SOX9 enhancer 517 to 595 kb upstream of the SOX9 gene.

Allelic variations in the vitamin D receptor gene, insulin secretion and parents' heights are independently associated with height in obese children and adolescents

Polymorphisms in the VDR gene were reported to be associated with variations in intrauterine and postnatal growth and with adult height, but also with other traits that are strongly correlated such as the BMI, insulin sensitivity, insulin secretion and hyperglycemia. Here, we assessed the impact of VDR polymorphisms on body height and its interactions with obesity- and glucose tolerance-related traits in obese children and adolescents. We studied 173 prepubertal (Tanner's stage 1) and 146 pubertal (Tanner's stages 2-5) obese children who were referred for a weight-loss program. Three single nucleotide polymorphisms were genotyped: rs1544410 (BsmI), rs7975232 (ApaI) and rs731236 (TaqI). BsmI and TaqI genotypes were significantly associated with height in pubertal children, but the associations did not reach statistical significance in prepubertal children. In stepwise regression analyses, the lean body mass, insulin secretion, BsmI or TaqI genotypes and the father's and the mother's height were independently and positively associated with height in pubertal children. These covariables accounted for 46% of the trait variance. The height of homozygous carriers of the minor allele of BsmI was 0.65 z-scores (4 cm) higher than the height of homozygous carriers of the major allele (P=.0006). Haplotype analyses confirmed the associations of the minor alleles of BsmI and TaqI with increased height. In conclusion, VDR genotypes were significantly associated with height in pubertal obese children. The associations were independent from the effects of confounding traits, such as the body fat mass, insulin secretion, insulin sensitivity and glucose tolerance. (C) 2012 Elsevier Inc. All rights reserved.

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.

Functional single-nucleotide polymorphisms in the DEFB1 gene are associated with systemic lupus erythematosus in Southern Brazilians

Systemic lupus erythematosus (SLE) is an autoimmune disease that results in inflammation and tissue damage. The etiology of SLE remains unknown, but recent studies have shown that the innate immune system may have a role in SLE pathogenesis through the secretion of small cationic peptides named defensins. The aim of the study was to determine the possible involvement in SLE of three functional single nucleotide polymorphisms (SNPs) (c.-52G>A, c.-44C>G and c.-20G>A) in the 5'UTR region of DEFB1 gene, by analyzing them in a population of 139 SLE patients and 288 healthy controls. The c.-52G>A SNP showed significant differences in allele and genotype frequency distribution between SLE patients and controls (p = 0.01 and p = 0.02 respectively) indicating protection against SLE (A allele, OR = 0.68, AA genotype OR = 0.51). Significant differences were also observed for c.-44C>G SNP, the C/G genotype being associated with susceptibility to SLE (OR = 1.60, p = 0.04). Moreover, statistically significant differences between patients and controls were found for two DEFB1 haplotypes (GCA and GGG, p = 0.01 and p = 0.02 respectively). When considering DEFB1 SNPs and SLE clinical and laboratory manifestations, significant association was found with neuropsychiatric disorders, immunological alterations and anti-DNA antibodies. In conclusion, our results evidence a possible role for the c.-52G>A and c.-44C>G DEFB1 polymorphisms in SLE pathogenesis, that can be considered as possible risk factors for development of disease and disease-related clinical manifestations. Additional studies are needed, to corroborate these results as well as functional studies to understand the biological role of these SNPs in the pathogenesis of SLE. Lupus (2012) 21, 625-631.

Gene Expression in Bovine Oocytes and Cumulus Cells After Meiotic Inhibition with the Cyclin-Dependent Kinase Inhibitor Butyrolactone I

Contents The aim of this study was to determine the effect of temporary inhibition of meiosis using the cyclin-dependent kinase inhibitor butyrolactone I (BLI) on gene expression in bovine oocytes and cumulus cells. Immature bovine cumulusoocyte complexes (COCs) were assigned to groups: (i) Control COCs collected immediately after recovery from the ovary or (ii) after in vitro maturation (IVM) for 24 h, (iii) Inhibited COCs collected 24 h after incubation with 100 mu m BLI or (iv) after meiotic inhibition for 24 h followed by IVM for a further 22 h. For mRNA relative abundance analysis, pools of 10 denuded oocytes and respective cumulus cells were collected. Transcripts related to cell cycle regulation and oocyte competence were evaluated in oocytes and cumulus cells by quantitative real-time PCR (qPCR). Most of the examined transcripts were downregulated (p < 0.05) after IVM in control and inhibited oocytes (19 of 35). Nine transcripts remained stable (p > 0.05) after IVM in control oocytes; only INHBA did not show this pattern in inhibited oocytes. Seven genes were upregulated after IVM in control oocytes (p < 0.05), and only PLAT, RBP1 and INHBB were not upregulated in inhibited oocytes after IVM. In cumulus cells, six genes were upregulated (p < 0.05) after IVM and eight were downregulated (p < 0.05). Cells from inhibited oocytes showed the same pattern of expression regarding maturation profile, but were affected by the temporary meiosis inhibition of the oocyte when the same maturation stages were compared between inhibited and control groups. In conclusion, changes in transcript abundance in oocytes and cumulus cells during maturation in vitro were mostly mirrored after meiotic inhibition followed by maturation.

Gene copy number and differential gene expression in haploid and diploid males of the stingless bee, Melipona quadrifasciata

Complementary sex determination in Hymenoptera implies that heterozygosity at the sex locus leads to the development of diploid females, whereas hemizygosity results in haploid males. Diploid males can arise through inbreeding. In social species, these pose a double burden on colony fitness, from significant reduction in its worker force and through being less viable and fertile than haploid males. Apart from being "misfits", diploid males are of interest to assess molecular correlates for possibly ploidy-related bionomic differences. Herein, we generated suppression subtractive cDNA libraries from newly emerged haploid and diploid males of the stingless bee Melipona quadrifasciata to enrich for differentially expressed genes. Gene Ontology classification revealed that in haploid males more DEGs were related to stress responsiveness, biosynthetic processes, reproductive processes and spermatogenesis, whereas in diploid ones differentially expressed genes were associated with cellular organization, nervous system development and amino acid transport were prevalent. Furthermore, both libraries contained over 40 % ESTs representing possibly novel transcripts. Quantitative RT-PCR analyses confirmed the differential expression of a representative DEG set in newly emerged males. Several muscle formation and energy metabolism-related genes were under-expressed in diploid males. On including 5-day-old males in the analysis, changes in transcript abundance during sexual maturation were revealed.