Validation of reference genes for gene expression studies in the honey bee, Apis mellifera, by quantitative real-time RT-PCR

For obtaining accurate and reliable gene expression results it is essential that quantitative real-time RT-PCR (qRT-PCR) data are normalized with appropriate reference genes. The current exponential increase in postgenomic studies on the honey bee, Apis mellifera, makes the standardization of qRT-PCR results an important task for ongoing community efforts. For this aim we selected four candidate reference genes (actin, ribosomal protein 49, elongation factor 1-alpha, tbp-association factor) and used three software-based approaches (geNorm, BestKeeper and NormFinder) to evaluate the suitability of these genes as endogenous controls. Their expression was examined during honey bee development, in different tissues, and after juvenile hormone exposure. Furthermore, the importance of choosing an appropriate reference gene was investigated for two developmentally regulated target genes. The results led us to consider all four candidate genes as suitable genes for normalization in A. mellifera. However, each condition evaluated in this study revealed a specific set of genes as the most appropriated ones.

Sequential analysis of global gene expression profiles in immature and in vitro matured bovine oocytes: potential molecular markers of oocyte maturation

Background: Without intensive selection, the majority of bovine oocytes submitted to in vitro embryo production (IVP) fail to develop to the blastocyst stage. This is attributed partly to their maturation status and competences. Using the Affymetrix GeneChip Bovine Genome Array, global mRNA expression analysis of immature (GV) and in vitro matured (IVM) bovine oocytes was carried out to characterize the transcriptome of bovine oocytes and then use a variety of approaches to determine whether the observed transcriptional changes during IVM was real or an artifact of the techniques used during analysis. Results: 8489 transcripts were detected across the two oocyte groups, of which similar to 25.0% (2117 transcripts) were differentially expressed (p < 0.001); corresponding to 589 over-expressed and 1528 under-expressed transcripts in the IVM oocytes compared to their immature counterparts. Over expression of transcripts by IVM oocytes is particularly interesting, therefore, a variety of approaches were employed to determine whether the observed transcriptional changes during IVM were real or an artifact of the techniques used during analysis, including the analysis of transcript abundance in oocytes in vitro matured in the presence of a-amanitin. Subsets of the differentially expressed genes were also validated by quantitative real-time PCR (qPCR) and the gene expression data was classified according to gene ontology and pathway enrichment. Numerous cell cycle linked (CDC2, CDK5, CDK8, HSPA2, MAPK14, TXNL4B), molecular transport (STX5, STX17, SEC22A, SEC22B), and differentiation (NACA) related genes were found to be among the several over-expressed transcripts in GV oocytes compared to the matured counterparts, while ANXA1, PLAU, STC1and LUM were among the over-expressed genes after oocyte maturation. Conclusion: Using sequential experiments, we have shown and confirmed transcriptional changes during oocyte maturation. This dataset provides a unique reference resource for studies concerned with the molecular mechanisms controlling oocyte meiotic maturation in cattle, addresses the existing conflicting issue of transcription during meiotic maturation and contributes to the global goal of improving assisted reproductive technology.

A novel COL1A1 gene-splicing mutation (c. 1875+1G > C) in a Brazilian patient with osteogenesis imperfecta

Osteogenesis imperfecta is a heterogeneous genetic disorder characterized by bone fragility and deformity, recurrent fractures, blue sclera, short stature, and dentinogenesis imperfecta. Most cases are caused by mutations in COL1A1 and COL1A2 genes. We present a novel splicing mutation in the COL1A1 gene (c. 1875+ 1G>C) in a 16-year-old Brazilian boy diagnosed as a type III osteogenesis imperfecta patient. This splicing mutation and its association with clinical phenotypes will be submitted to the reference database of COL1A1 mutations, which has no other description of this mutation.

Atorvastatin effects on SREBF1a and SCAP gene expression in mononuclear cells and its relation with lowering-lipids response

Background: The transcription factors SREBP1 and SCAP are involved in intracellular cholesterol homeostasis. Polymorphisms of these genes have been associated with variations on serum lipid levels and response to statins that are potent cholesterol-lowering drugs. We evaluated the effects of atorvastatin on SREBF1a and SCAP mRNA expression in peripheral blood mononuclear cells (PBMC) and a possible association with gene polymorphisms and lowering-cholesterol response. Methods: Fifty-nine hypercholesterolemic patients were treated with atorvastatin (10 mg/day for 4 weeks). Serum lipid profile and mRNA expression in PBMC were assessed before and after the treatment. Gene expression was quantified by real-time PCR using GAPD as endogenous reference and mRNA expression in HepG2 cells as calibrator. SREBF1 -36delG and SCAP A2386G polymorphisms were detected by PCR-RFLP. Results: Our results showed that transcription of SREBF1a and SCAP was coordinately regulated by atorvastatin (r=0.595, p<0.001), and that reduction in SCAP transcription was associated with the 2386AA genotype (p=0.019). Individuals who responded to atorvastatin with a downregulation of SCAP had also a lower triglyceride compared to those who responded to atorvastatin with an upregulation of SCAP. Conclusion: Atorvastatin has differential effects on SREBF1a and SCAP mRNA expression in PBMC that are associated with baseline transcription levels, triglycerides response to atorvastatin and SCAP A2386G polymorphism. (c) 2008 Elsevier B.V. All rights reserved.

HER2 Testing in Breast Carcinoma Very Low Concordance Rate Between Reference and Local Laboratories in Brazil

Breast cancer accounts for approximately one quarter of all cancers in females. HER2 gene amplification or HER2 protein overexpression, detected in about 20% of breast carcinomas, predicts a more aggressive clinical course and determines eligibility for targeted therapy with trastuzumab. HER2 testing has become an essential part of the clinical evaluation of all breast carcinoma patients, and accurate HER2 results are critical in identifying patients who may be benefited from targeted therapy. This study investigated the concordance in the results of HER2 immunohistochemistry assays performed in 500 invasive breast carcinomas between a reference laboratory and 149 local laboratories from all geographic regions of Brazil. Our results showed an overall poor concordance (171 of 500 cases, 34.2%) regarding HER2 results between local and reference laboratories, which may be related to the low-volume load of HER2 assays, inexperience with HER2 scoring system, and/or technical issues related to immunohistochemistry in local laboratories. Standardization of HER2 testing with rigorous quality control measures by local laboratories is highly recommended to avoid erroneous treatment of breast cancer patients.

Shared and Unique Gene Expression in Systemic Lupus Erythematosus Depending on Disease Activity

Patients presenting with active Systemic lupus erythematosus (SLE) manifestations may exhibit distinct pathogenetic features in relation to inactive SLE. Also, cDNA microarrays may potentially discriminate the gene expression profile of a disease or disease variant. Therefore, we evaluated the expression profile of 4500 genes in peripheral blood lymphocytes (PBL) of SLE patients. We studied 11 patients with SLE (seven with active SLE and four with inactive SLE) and eight healthy controls. Total RNA was isolated from PBL, reverse transcribed into cDNA, and postlabeled with Cy3 fluorochrome. These probes were then hybridized to a glass slide cDNA microarray containing 4500 human IMAGE cDNA target sequences. An equimolar amount of total RNA from human cell lines served as reference. The microarray images were quantified, normalized, and analyzed using the R environment (ANOVA, significant analysis of microarrays, and cluster-tree view algorithms). Disease activity was assessed by the SLE disease activity index. Compared to the healthy controls, 104 genes in active SLE patients (80 repressed and 24 induced) and 52 genes in nonactive SLE patients (31 induced and 21 repressed) were differentially expressed. The modulation of 12 genes, either induced or repressed, was found in both disease variants; however, each disease variant had differential expression of different genes. Taken together, these results indicate that the two lupus variants studied have common and unique differentially expressed genes. Although the biological significance of the differentially expressed genes discussed above has not been completely understood, they may serve as a platform to further explore the molecular basis of immune deregulation in SLE.

Affected and non-affected skin fibroblasts from systemic sclerosis patients share a gene expression profile deviated from the one observed in healthy individuals

Objectives To evaluate the gene expression profile of fibroblasts from affected and non-affected skin of systemic sclerosis (SSc) patients and from controls. Materials and methods Labeled cDNA from fibroblast cultures from forearm (affected) and axillary (non-affected) skin from six diffuse SSc patients, from three normal controls, and from MOLT-4/HEp-2/normal fibroblasts (reference pool) was probed in microarrays generated with 4193 human cDNAs from the IMAGE Consortium. Microarray images were converted into numerical data and gene expression was calculated as the ratio between fibroblast cDNA (Cy5) and reference pool cDNA (Cy3) data and analyzed by R environment/Aroma, Cluster, Tree View, and SAM softwares. Differential expression was confirmed by real time PCR for a set of selected genes. Results Eighty-eight genes were up- and 241 genes down-regulated in SSc fibroblasts. Gene expression correlation was strong between affected and non-affected fibroblast samples from the same patient (r>0.8), moderate among fibroblasts from all patients (r=0.72) and among fibroblasts from all controls (r=0.70), and modest among fibroblasts from patients and controls (r=0.55). The differential expression was confirmed by real time PCR for all selected genes. Conclusions Fibroblasts from affected and non-affected skin of SSc patients shared a similar abnormal gene expression profile, suggesting that the widespread molecular disturbance in SSc fibroblasts is more sensitive than histological and clinical alterations. Novel molecular elements potentially involved in SSc pathogenesis were identified.

Shared and differential traits in the accessory olfactory bulb of caviomorph rodents with particular reference to the semiaquatic capybara

The vomeronasal system is crucial for social and sexual communication in mammals. Two populations of vomeronasal sensory neurons, each expressing G alpha i2 or G alpha o proteins, send projections to glomeruli of the rostral or caudal accessory olfactory bulb, rAOB and cAOB, respectively. In rodents, the G alpha i2- and G alpha o-expressing vomeronasal pathways have shown differential responses to small/volatile vs. large/non-volatile semiochemicals, respectively. Moreover, early gene expression suggests predominant activation of rAOB and cAOB neurons in sexual vs. aggressive contexts, respectively. We recently described the AOB of Octodon degus, a semiarid-inhabiting diurnal caviomorph. Their AOB has a cell indentation between subdomains and the rAOB is twice the size of the cAOB. Moreover, their AOB receives innervation from the lateral aspect, contrasting with the medial innervation of all other mammals examined to date. Aiming to relate AOB anatomy with lifestyle, we performed a morphometric study on the AOB of the capybara, a semiaquatic caviomorph whose lifestyle differs remarkably from that of O. degus. Capybaras mate in water and scent-mark their surroundings with oily deposits, mostly for male-male communication. We found that, similar to O. degus, the AOB of capybaras shows a lateral innervation of the vomeronasal nerve, a cell indentation between subdomains and heterogeneous subdomains, but in contrast to O. degus the caudal portion is larger than the rostral one. We also observed that four other caviomorph species present a lateral AOB innervation and a cell indentation between AOB subdomains, suggesting that those traits could represent apomorphies of the group. We propose that although some AOB traits may be phylogenetically conserved in caviomorphs, ecological specializations may play an important role in shaping the AOB.

Toward a resolution of Campanulid phylogeny, with special reference to the placement of Dipsacales

Broad-scale phylogenetic analyses of the angiosperms and of the Asteridae have failed to confidently resolve relationships among the major lineages of the campanulid Asteridae (i.e., the euasterid II of APG II, 2003). To address this problem we assembled presently available sequences for a core set of 50 taxa, representing the diversity of the four largest lineages (Apiales, Aquifoliales, Asterales, Dipsacales) as well as the smaller ""unplaced"" groups (e.g., Bruniaceae, Paracryphiaceae, Columelliaceae). We constructed four data matrices for phylogenetic analysis: a chloroplast coding matrix (atpB, matK, ndhF, rbcL), a chloroplast non-coding matrix (rps16 intron, trnT-F region, trnV-atpE IGS), a combined chloroplast dataset (all seven chloroplast regions), and a combined genome matrix (seven chloroplast regions plus 18S and 26S rDNA). Bayesian analyses of these datasets using mixed substitution models produced often well-resolved and supported trees. Consistent with more weakly supported results from previous studies, our analyses support the monophyly of the four major clades and the relationships among them. Most importantly, Asterales are inferred to be sister to a clade containing Apiales and Dipsacales. Paracryphiaceae is consistently placed sister to the Dipsacales. However, the exact relationships of Bruniaceae, Columelliaceae, and an Escallonia clade depended upon the dataset. Areas of poor resolution in combined analyses may be partly explained by conflict between the coding and non-coding data partitions. We discuss the implications of these results for our understanding of campanulid phylogeny and evolution, paying special attention to how our findings bear on character evolution and biogeography in Dipsacales.

The mitochondrial genome of the stingless bee Melipona bicolor (Hymenoptera, Apidae, Meliponini): sequence, gene organization and a unique tRNA translocation event conserved across the tribe Meliponini

At present a complete mtDNA sequence has been reported for only two hymenopterans, the Old World honey bee, Apis mellifera and the sawfly Perga condei. Among the bee group, the tribe Meliponini (stingless bees) has some distinction due to its Pantropical distribution, great number of species and large importance as main pollinators in several ecosystems, including the Brazilian rain forest. However few molecular studies have been conducted on this group of bees and few sequence data from mitochondrial genomes have been described. In this project, we PCR amplified and sequenced 78% of the mitochondrial genome of the stingless bee Melipona bicolor (Apidae, Meliponini). The sequenced region contains all of the 13 mitochondrial protein-coding genes, 18 of 22 tRNA genes, and both rRNA genes (one of them was partially sequenced). We also report the genome organization (gene content and order), gene translation, genetic code, and other molecular features, such as base frequencies, codon usage, gene initiation and termination. We compare these characteristics of M. bicolor to those of the mitochondrial genome of A. mellifera and other insects. A highly biased A+T content is a typical characteristic of the A. mellifera mitochondrial genome and it was even more extreme in that of M. bicolor. Length and compositional differences between M. bicolor and A. mellifera genes were detected and the gene order was compared. Eleven tRNA gene translocations were observed between these two species. This latter finding was surprising, considering the taxonomic proximity of these two bee tribes. The tRNA Lys gene translocation was investigated within Meliponini and showed high conservation across the Pantropical range of the tribe.

A novel missense mutation p.L76P in the GJB2 gene causing nonsyndromic recessive deafness in a Brazilian family

Mutations in the GJB2 gene, encoding connexin 26 (Cx26), are a major cause of nonsyndromic recessive hearing loss in many countries. We report here on a novel point mutation in GJB2, p.L76P (c.227C>T), in compound heterozygosity with a c.35delG mutation, in two Brazilian sibs, one presenting mild and the other profound nonsyndromic neurosensorial hearing impairment. Their father, who carried a wild-type allele and a p.L76P mutation, had normal hearing. The mutation leads to the substitution of leucine (L) by proline (P) at residue 76, an evolutionarily conserved position in Cx26 as well as in other connexins. This mutation is predicted to affect the first extracellular domain (EC1) or the second transmembrane domain (TM2). EC1 is important for connexon-connexon interaction and for the control of channel voltage gating. The segregation of the c.227C>T (p.L76P) mutation together with c.35delG in this family indicates a recessive mode of inheritance. The association between the p.L76P mutation and hearing impairment is further supported by its absence in a normal hearing control group of 100 individuals, 50 European-Brazilians and 50 African-Brazilians.

Essential regulatory elements for NHE3 gene transcription in renal proximal tubule cells

The objectives of the present study were to identify the cis-elements of the promoter absolutely required for the efficient rat NHE3 gene transcription and to locate positive and negative regulatory elements in the 5’-flanking sequence (5’FS), which might modulate the gene expression in proximal tubules, and to compare this result to those reported for intestinal cell lines. We analyzed the promoter activity of different 5’FS segments of the rat NHE3 gene, in the OKP renal proximal tubule cell line by measuring the activity of the reporter gene luciferase. Because the segment spanning the first 157 bp of 5’FS was the most active it was studied in more detail by sequential deletions, point mutations, and gel shift assays. The essential elements for gene transcription are in the region -85 to -33, where we can identify consensual binding sites for Sp1 and EGR-1, which are relevant to NHE3 gene basal transcription. Although a low level of transcription is still possible when the first 25 bp of the 5’FS are used as promoter, efficient transcription only occurs with 44 bp of 5’FS. There are negative regulatory elements in the segments spanning -1196 to -889 and -467 to -152, and positive enhancers between -889 and -479 bp of 5’FS. Transcription factors in the OKP cell nuclear extract efficiently bound to DNA elements of rat NHE3 promoter as demonstrated by gel shift assays, suggesting a high level of similarity between transcription factors of both species, including Sp1 and EGR-1.

Hypothyroidism decreases proinsulin gene expression and the attachment of its mRNA and eEF1A protein to the actin cytoskeleton of INS-1E cells

The actions of thyroid hormone (TH) on pancreatic beta cells have not been thoroughly explored, with current knowledge being limited to the modulation of insulin secretion in response to glucose, and beta cell viability by regulation of pro-mitotic and pro-apoptotic factors. Therefore, the effects of TH on proinsulin gene expression are not known. This led us to measure: a) proinsulin mRNA expression, b) proinsulin transcripts and eEF1A protein binding to the actin cytoskeleton, c) actin cytoskeleton arrangement, and d) proinsulin mRNA poly(A) tail length modulation in INS-1E cells cultured in different media containing: i) normal fetal bovine serum - FBS (control); ii) normal FBS plus 1 µM or 10 nM T3, for 12 h, and iii) FBS depleted of TH for 24 h (Tx). A decrease in proinsulin mRNA content and attachment to the cytoskeleton were observed in hypothyroid (Tx) beta cells. The amount of eEF1A protein anchored to the cytoskeleton was also reduced in hypothyroidism, and it is worth mentioning that eEF1A is essential to attach transcripts to the cytoskeleton, which might modulate their stability and rate of translation. Proinsulin poly(A) tail length and cytoskeleton arrangement remained unchanged in hypothyroidism. T3 treatment of control cells for 12 h did not induce any changes in the parameters studied. The data indicate that TH is important for proinsulin mRNA expression and translation, since its total amount and attachment to the cytoskeleton are decreased in hypothyroid beta cells, providing evidence that effects of TH on carbohydrate metabolism also include the control of proinsulin gene expression.

Pre-processing for noise detection in gene expression classification data

Due to the imprecise nature of biological experiments, biological data is often characterized by the presence of redundant and noisy data. This may be due to errors that occurred during data collection, such as contaminations in laboratorial samples. It is the case of gene expression data, where the equipments and tools currently used frequently produce noisy biological data. Machine Learning algorithms have been successfully used in gene expression data analysis. Although many Machine Learning algorithms can deal with noise, detecting and removing noisy instances from the training data set can help the induction of the target hypothesis. This paper evaluates the use of distance-based pre-processing techniques for noise detection in gene expression data classification problems. This evaluation analyzes the effectiveness of the techniques investigated in removing noisy data, measured by the accuracy obtained by different Machine Learning classifiers over the pre-processed data.

A phylogenetic analysis of Brycon and Henochilus (Characiformes, Characidae, Bryconinae) based on the mitochondrial gene 16S rRNA

The genus Brycon, the largest subunit of the Bryconinae, has 42 valid species distributed from southern Mexico to the La Plata River in Argentina. Henochilus is a monotypic genus, comprising a single species (H. wheatlandii) found in the upper Rio Doce basin. In the present study, partial sequences of the mitochondrial gene 16S were obtained for fifteen species of Brycon and for Henochilus wheatlandii. The results showed that the genus Brycon is paraphyletic, since Henochilus is the sister-group of B. ferox and B. insignis. The most basal species analyzed were the trans-Andean species B. henni, B. petrosus, and B. chagrensis.