Overexpression of DR-nm23, a protein encoded by a member of the nm23 gene family, inhibits granulocyte differentiation and induces apoptosis in 32Dc13 myeloid cells.

Chronic myelogenous leukemia evolves in two clinically distinct stages: a chronic and a blast crisis phase. The molecular changes associated with chronic phase to blast crisis transition are largely unknown. We have identified a cDNA clone, DR-nm23, differentially expressed in a blast-crisis cDNA library, which has approximately 70% sequence similarity to the putative metastatic suppressor genes, nm23-H1 and nm23-H2. The deduced amino acid sequence similarity to the proteins encoded by these two latter genes is approximately 65% and includes domains and amino acid residues (the leucine zipper-like and the RGD domain, a serine and a histidine residue in the NH2- and in the COOH-terminal portion of the protein, respectively) postulated to be important for nm23 function. DR-nm23 mRNA is preferentially expressed at early stages of myeloid differentiation of highly purified CD34+ cells. Its constitutive expression in the myeloid precursor 32Dc13 cell line, which is growth-factor dependent for both proliferation and differentiation, results in inhibition of granulocytic differentiation induced by granulocyte colony-stimulating factor and causes apoptotic cell death. These results are consistent with a role for DR-nm23 in normal hematopoiesis and raise the possibility that its overexpression contributes to differentiation arrest, a feature of blastic transformation in chronic myelogenous leukemia.

Perfil da expressão gênica de larvas de Tetrapedia diversipes (Hymenoptera: Apidae) em diapausa

A diapausa é um fenômeno amplamente presente nos artrópodes e é considerada como primordial para o sucesso evolutivo da Classe Insecta, pois possibilita a sobrevivência em condições adversas, como estações frias e secas. Sabe-se que durante a diapausa ocorre o silenciamento de muitos genes e que outros são unicamente expressos nesta fase. Embora existam evidências de que o processo da diapausa tenha se mantido conservado durante a evolução das espécies, ainda há lacunas no conhecimento sobre o nível de conservação dos padrões metabólicos. Um bom modelo para se estudar a diapausa é Tetrapedia diversipes, uma espécie bivoltina de abelha solitária. Os indivíduos que nascem na primeira geração seguem o desenvolvimento desde ovo até adulto em tempo bem menor do que aqueles que nascem na segunda geração; estes retardam o desenvolvimento na fase larval. Além disso, essa espécie é de fácil obtenção no seu ambiente natural, pois apresenta alta taxa de nidificação em ninhos-armadilha. O objetivo deste trabalho foi comparar o perfil de expressão de genes entre as larvas da 1ª geração (que não entram em diapausa), larvas da 2ª geração (que entrariam em diapausa) e das larvas em diapausa. Foram identificados 196 genes diferencialmente expressos, destes 87 foram anotados. Muitos destes genes já foram descritos na literatura como relacionados à diapausa em outras espécies, no entanto, o padrão de expressão não é conservado. Os genes aqui identificados foram divididos em cinco grupos: relacionados à desintoxicação celular, cutícula e citoesqueleto, metabolismo de lipídeos e esteróis, ciclo celular e outros genes relacionados à diapausa

Estudos sobre a organização genômica, evolução e expressão de microRNAs

Os microRNAs (miRNAs) são pequenos RNAs não codificadores de proteínas presentes na maioria dos eucariotos. Esses RNAs regulam a expressão gênica em nível pós-transcricional através do silenciamento de mRNAs-alvo que possuem sítios complementares às suas sequências, atuando em praticamente todos os processos celulares. Embora a estrutura e função dos miRNAs estejam bem caracterizadas, aspectos relacionados à sua organização genômica, evolução e atuação em doenças são tópicos que apresentam enormes lacunas. Nesta tese, utilizamos abordagens computacionais para investigar estes temas em três trabalhos. No primeiro, processamos e integramos um vasto volume de dados publicamente disponíveis referentes aos miRNAs e genes codificadores de proteínas para cinco espécies de vertebrados. Com isso, construimos uma ferramenta web que permite a fácil inspeção da organização genômica dos miRNAs em regiões inter e intragênicas, o acesso a dados de expressão de miRNAs e de genes codificadores de proteínas (classificados em genes hospedeiros e não hospedeiros de miRNAs), além de outras informações pertinentes. Verificamos que a ferramenta tem sido amplamente utilizada pela comunidade científica e acreditamos que ela possa facilitar a geração de hipóteses associadas à regulação dos miRNAs, principalmente quando estão inseridos em genes hospedeiros. No segundo estudo, buscamos compreender como o contexto genômico e a origem evolutiva dos genes hospedeiros influenciam a expressão e evolução dos miRNAs humanos. Nossos achados mostraram que os miRNAs intragênicos surgem preferencialmente em genes antigos (origem anterior à divergência de vertebrados). Observamos que os miRNAs inseridos em genes antigos têm maior abrangência de expressão do que os inseridos em genes novos. Surpreendentemente, miRNAs jovens localizados em genes antigos são expressos em um maior número de tecidos do que os intergênicos de mesma idade, sugerindo uma vantagem adaptativa inicial que pode estar relacionada com o controle da expressão dos genes hospedeiros, e como consequência, expondo-os a contextos celulares e conjuntos de alvos diversos. Na evolução a longo prazo, vimos que genes antigos conferem maior restrição nos padrões de expressão (menor divergência de expressão) para miRNAs intragênicos, quando comparados aos intergênicos. Também mostramos possíveis associações funcionais relacionadas ao contexto genômico, tais como o enriquecimento da expressão de miRNAs intergênicos em testículo e dos intragênicos em tecidos neurais. Propomos que o contexto genômico e a idade dos genes hospedeiros são fatores-chave para a evolução e expressão dos miRNAs. Por fim, buscamos estabelecer associações entre a expressão diferencial de miRNAs e a quimioresistência em câncer colorretal utilizando linhagens celulares sensíveis e resistentes às drogas 5-Fluoruracil e Oxaliplatina. Dentre os miRNAs identificados, o miR-342 apresentou níveis elevados de expressão nas linhagens sensíveis à Oxaliplatina. Com base na análise dos alvos preditos, detectamos uma significativa associação de miR-342 com a apoptose. A superexpressão de miR-342 na linhagem resistente SW620 evidenciou alterações na expressão de genes da via apoptótica, notavelmente a diminuição da expressão do fator de crescimento PDGFB, um alvo predito possivelmente sujeito à regulação direta pelo miR-342.

Investigação de mecanismos genéticos e epigenéticos em distúrbios do crescimento humano

Parcela considerável de pacientes com distúrbios de crescimento não têm a causa de seus quadros clínicos estabelecida, incluindo aproximadamente 50% dos pacientes com diagnóstico clínico de síndrome de Silver−Russell (SRS) e 10-20% dos pacientes com síndrome de Beckwith-Wiedemann (BWS). O objetivo deste estudo foi investigar as causas genéticas e epigenéticas de distúrbios de crescimento, de etiologia desconhecida, numa contribuição para o entendimento de mecanismos que regulam o crescimento. O estudo compreendeu: (1) a investigação de microdesequilíbrios cromossômicos, por aCGH; (2) a análise do perfil de expressão alelo-específica de genes sujeitos a imprinting (IG), por pirossequenciamento (PSQ) ou sequenciamento de Sanger; (3) a investigação do padrão de metilação global em pacientes com restrição de crescimento, utilizando microarray de metilação. A casuística constituiu-se de 41 pacientes não aparentados, com distúrbios de crescimento, de etiologia desconhecida: (1) 25, com hipótese diagnóstica de SRS; (2) seis, com restrição de crescimento intrauterino e peso ao nascimento abaixo do 10º percentil, associados a outros sinais clínicos; (3) sete, com hipótese diagnóstica de BWS; e (4) três, com macrossomia pré-natal ou pós-natal, associada a outros sinais. A investigação de microdesequilíbrios cromossômicos foi realizada em 40 pacientes. Foram detectadas 58 variantes raras em 30/40 pacientes (75%): 40 foram consideradas provavelmente benignas (18 pacientes, 45%), 12, com efeito patogênico desconhecido (11 pacientes, 27,5%), duas, provavelmente patogênicas (um paciente, 2,5%) e quatro, patogênicas (três pacientes, 7,5%). Essas frequências são comparáveis àquelas descritas em estudos que investigaram CNV em grupos de pacientes com distúrbios de crescimento e outras alterações congênitas, incluindo SRS, e mostram a importância da investigação de microdesequilíbrios cromossômicos nesses pacientes. A diversidade dos microdesequilíbrios cromossômicos identificados é reflexo da heterogeneidade clínica das casuísticas. Neste estudo, muitos dos pacientes com hipótese diagnóstica de SRS e BWS apresentavam sinais clínicos atípicos, explicando a ausência neles das alterações (epi)genéticas que causam essas síndromes. A identificação de CNV características de outras síndromes reflete a sobreposição de sinais clínicos com BWS e SRS. A análise do perfil de expressão alelo-específica de IG foi realizada em um subgrupo de 18 pacientes com restrição de crescimento. Trinta IG com função em proliferação celular, crescimento fetal ou neurodesenvolvimento foram inicialmente selecionados. Após seleção de SNP transcritos com alta frequência na população, genotipagem de pacientes, genitores e indivíduos controle, determinação da expressão dos IG em sangue periférico e seu padrão de expressão (mono ou bialélico), 13 IG, expressos no sangue, tiveram a expressão alelo-específica avaliada, sete deles por PSQ e seis por sequenciamento de Sanger. Alterações no perfil de expressão de dois genes, de expressão normalmente paterna, foram detectadas em 4/18 pacientes (22%). Este estudo é o primeiro a utilizar pirossequenciamento e sequenciamento de Sanger na avaliação do perfil de expressão alelo-específica de IG, em pacientes com restrição de crescimento. Apesar de terem limitações, ambas as técnicas mostraram-se robustas e revelaram alterações de expressão alélica interessantes; entretanto, a relação dessas alterações com o quadro clínico dos pacientes permanece por esclarecer. A investigação da metilação global do DNA foi realizada em subgrupo de 21 pacientes com restrição de crescimento e em 24 indivíduos controle. Dois tipos de análise foram realizados: (1) análise diferencial de grupo e (2) análise diferencial individual. Na primeira análise, em que foi comparado o padrão de metilação do grupo de pacientes com quadro clínico sugestivo de SRS (n=16) com o do grupo controle (n=24), não houve indicação de hipo ou hipermetilação global no grupo SRS. Na segunda análise, foi comparado o padrão de metilação de cada um dos 21 pacientes com restrição de crescimento e dos 24 indivíduos controle, com o padrão de metilação do grupo controle. O número médio de CpG hipermetilados e de segmentos diferencialmente metilados (SDM) foi significativamente maior nos pacientes. Foram identificados 82 SDM hipermetilados, estando 57 associados a gene(s) (69,5%), em 16 pacientes, e 51 SDM hipometilados, 41 deles associados a gene(s) (80,4%), em 10 pacientes. A análise de ontologia genética dos 61 genes associados aos SDM hipo ou hipermetilados nos pacientes destacou genes que atuam no desenvolvimento e na morfogênese do sistema esquelético e de órgãos fetais, e na regulação da transcrição gênica e de processos metabólicos. Alterações de metilação em genes que atuam em processos de proliferação e diferenciação celulares e crescimento foram identificadas em 9/20 dos pacientes (45%), sugerindo implicação clínica. Não foi detectada alteração epigenética comum aos pacientes com diagnóstico clínico de SRS, explicável provavelmente pela heterogeneidade clínica. A investigação de metilação global, utilizando microarray, produziu novos dados que podem contribuir para a compreensão de mecanismos moleculares que influenciam o crescimento pré- e pós-natal. Na translocação aparentemente equilibrada - t(5;6)(q35.2;p22.3)dn, detectada em paciente com suspeita clínica de SRS, a interrupção de um gene, pela quebra no cromossomo 6, pode ser a causa do quadro clínico; alternativamente, a translocação pode ter impactado a regulação de genes de desenvolvimento localizados próximos aos pontos de quebra. A análise de expressão em sangue periférico mostrou que os níveis de cDNA do gene, interrompido pelo ponto de quebra da translocação, estavam reduzidos à metade. Além de sinais típicos da SRS, a paciente apresentava algumas características clínicas sugestivas de displasia cleidocraniana. Assim, a translocação t(5;6) pode ter alterado a interação de genes de desenvolvimento e seus elementos reguladores, levando à desregulação de sua expressão espaço-temporal, e resultando num fenótipo atípico, com características sobrepostas de mais de uma síndrome genética

Developmental control of H+/amino acid permease gene expression during seed development of Arabidopsis

Long distance transport of amino acids is mediated by several families of differentially expressed amino acid transporters. The two genes AAP1 and AAP2 encode broad specificity H+-amino acid co-transporters and are expressed to high levels in siliques of Arabidopsis, indicating a potential role in supplying the seeds with organic nitrogen. The expression of both genes is developmentally controlled and is strongly induced in siliques at heart stage of embryogenesis, shortly before induction of storage protein genes. Histochemical analysis of transgenic plants expressing promoter-GUS fusions shows that the genes have non-overlapping expression patterns in siliques. AAP1 is expressed in the endosperm and the cotyledons whereas AAP2 is expressed in the vascular strands of siliques and in funiculi. The endosperm expression of AAP1 during early stages of seed development indicates that the endosperm serves as a transient storage tissue for organic nitrogen. Amino acids are transported in both xylem and phloem but during seed filling are imported only via the phloem. AAP2, which is expressed in the phloem of stems and in the veins supplying seeds, may function in uptake of amino acids assimilated in the green silique tissue, in the retrieval of amino acids leaking passively out of the phloem and in xylem-to-phloem transfer along the path. The promoters provide excellent tools to study developmental, hormonal and metabolic control of nitrogen nutrition during development and may help to manipulate the timing and composition of amino acid import into seeds.

Transgenerational effects persist down the maternal line in marine sticklebacks: gene expression matches physiology in a warming ocean, supplementary material

Transgenerational effects can buffer populations against environmental change, yet little is known about underlying mechanisms, their persistence, or the influence of environmental cue timing. We investigated mitochondrial respiratory capacity (MRC) and gene expression of marine sticklebacks that experienced acute or developmental acclimation to simulated ocean warming (21°C) across three generations. Previous work showed that acute acclimation of grandmothers to 21°C led to lower (optimised) offspring MRCs. Here, developmental acclimation of mothers to 21°C led to higher, but more efficient offspring MRCs. Offspring with a 21°Cx17°C grandmother-mother environment mismatch showed metabolic compensation: their MRCs were as low as offspring with a 17°C thermal history across generations. Transcriptional analyses showed primarily maternal but also grandmaternal environment effects: genes involved in metabolism and mitochondrial protein biosynthesis were differentially expressed when mothers developed at 21°C, whereas 21°C grandmothers influenced genes involved in hemostasis and apoptosis. Genes involved in mitochondrial respiration all showed higher expression when mothers developed at 21° and lower expression in the 21°Cx17°C group, matching the phenotypic pattern for MRCs. Our study links transcriptomics to physiology under climate change, and demonstrates that mechanisms underlying transgenerational effects persist across multiple generations with specific outcomes depending on acclimation type and environmental mismatch between generations.

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.

Continued discovery of transcriptional units expressed in cells of the mouse mononuclear phagocyte lineage

The current RIKEN transcript set represents a significant proportion of the mouse transcriptome but transcripts expressed in the innate and acquired immune systems are poorly represented. In the present study we have assessed the complexity of the transcriptome expressed in mouse macrophages before and after treatment with lipopolysaccharide, a global regulator of macrophage gene expression, using existing RIKEN 19K arrays. By comparison to array profiles of other cells and tissues, we identify a large set of macrophage-enriched genes, many of which have obvious functions in endocytosis and phagocytosis. In addition, a significant number of LPS-inducible genes were identified. The data suggest that macrophages are a complex source of mRNA for transcriptome studies. To assess complexity and identify additional macrophage expressed genes, cDNA libraries were created from purified populations of macrophage and dendritic cells, a functionally related cell type. Sequence analysis revealed a high incidence of novel mRNAs within these cDNA libraries. These studies provide insights into the depths of transcriptional complexity still untapped amongst products of inducible genes, and identify macrophage and dendritic cell populations as a starting point for sampling the inducible mammalian transcriptome.

Identifying the molecular phenotype of renal progenitor cells

Although many of the molecular interactions in kidney development are now well understood, the molecules involved in the specification of the metanephric mesenchyme from surrounding intermediate mesoderm and, hence, the formation of the renal progenitor population are poorly characterized. In this study, cDNA microarrays were used to identify genes enriched in the murine embryonic day 10.5 (E10.5) uninduced metanephric mesenchyme, the renal progenitor population, in comparison with more rostral derivatives of the intermediate mesoderm. Microarray data were analyzed using R statistical software to determine accurately genes differentially expressed between these populations. Microarray outliers were biologically verified, and the spatial expression pattern of these genes at E10.5 and subsequent stages of early kidney development was determined by RNA in situ hybridization. This approach identified 21 genes preferentially expressed by the E10.5 metanephric mesenchyme, including Ewing sarcoma homolog, 14-3-3 theta, retinoic acid receptor-alpha, stearoyl-CoA desaturase 2, CD24, and cadherin-11, that may be important in formation of renal progenitor cells. Cell surface proteins such as CD24 and cadherin-11 that were strongly and specifically expressed in the uninduced metanephric mesenchyme and mark the renal progenitor population may prove useful in the purification of renal progenitor cells by FACS. These findings may assist in the isolation and characterization of potential renal stem cells for use in cellular therapies for kidney disease.

Correlating gene expression with larval competence and the effect of age and parentage on metamorphosis in the tropical abalone Haliotis asinina

The non-geniculate crustose coralline alga (CCA) Mastophora pacifica can induce the metamorphosis of competent Haliotis asinina (Vetigastropoda) larvae. The ability to respond to this natural cue varies considerably with larval age, with a higher proportion of older larvae (e.g. 90 h) able to metamorphose in response to M. pacifica than younger larvae (e.g. 66 h). Here we document the variation in time to acquisition of competence within a larval age class. For example, after 18 h of exposure to M. pacifica, approximately 15 and 36% of 84 and 90-h-old H. asinina larvae had initiated metamorphosis, respectively. This age-dependent response to M. pacifica is also observed when different aged larvae are exposed to CCA for varying periods. A higher proportion of older larvae require shorter periods of exposure to CCA than younger larvae in order to initiate metamorphosis. In this experiment, as in the previous, a small proportion of young larvae were able to respond to brief periods of CCA exposure, suggesting that they had developed the same state of competency as the majority of their older counterparts. Comparisons of the proportions of larvae undergoing metamorphosis between families reveals that parentage also has a significant (P < 0.05) affect on whether an individual will initiate metamorphosis at a given age. These familial differences are more pronounced when younger, largely pre-competent larvae (i.e. 66 h old) are exposed to M. pacifica, with proportions of larvae undergoing metamorphosis differing by as much as 10 fold between families. As these data suggest that variation in the rate of development of the competent state has a genetic basis, and as a first step towards identifying the molecular basis to this variation, we have identified numerous genes that are differentially expressed later in larval development using a differential display approach. Spatial expression analysis of these genes suggests that they may be directly involved in the acquisition of competence, or may play a functional role in the postlarva following metamorphosis.

Expression analysis of a human hepatic cell line in response to palmitate

Saturated fat plays a role in common debilitating diseases such as obesity, type 2 diabetes, and coronary heart disease. It is also clear that certain fatty acids act as regulators of metabolism via both direct and indirect signalling of target tissues. As the molecular mechanisms of saturated fatty acid signalling in the liver are poorly defined, hepatic gene expression analysis was undertaken in a human hepatocyte cell line after incubation with palmitate. Profiling of mRNA expression using cDNA microarray analysis revealed that 162 of approximately 18,000 genes tested were differentially expressed after incubation with palmitate for 48 h. Altered transcription profiles were observed in a wide variety of genes, including genes involved in lipid and cholesterol transport, cholesterol catabolism, cell growth and proliferation, cell signalling, P-oxidation, and oxidative stress response. While palinitate signalling has been examined in pancreatic beta-cells, this is the first report showing that palmitate regulates expression of numerous genes via direct molecular signalling mechanisms in liver cells. (C) 2005 Elsevier Inc. All rights reserved.

A simple implementation of a normal mixture approach to differential gene expression in multiclass microarrays

Motivation: An important problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. We provide a straightforward and easily implemented method for estimating the posterior probability that an individual gene is null. The problem can be expressed in a two-component mixture framework, using an empirical Bayes approach. Current methods of implementing this approach either have some limitations due to the minimal assumptions made or with more specific assumptions are computationally intensive. Results: By converting to a z-score the value of the test statistic used to test the significance of each gene, we propose a simple two-component normal mixture that models adequately the distribution of this score. The usefulness of our approach is demonstrated on three real datasets.

Parasitic castration by the digenian trematode Allopodocotyle sp alters gene expression in the brain of the host mollusc Haliotis asinina

Infection of molluscs by digenean trematode parasites typically results in the repression of reproduction - the so-called parasitic castration. This is known to occur by altering the expression of a range of host neuropeptide genes. Here we analyse the expression levels of 10 members of POU, Pax, Sox and Hox transcription factor gene families, along with genes encoding FNIRFamide, prohormone convertase and P-tubulin, in the brain ganglia of actively reproducing (summer), non-reproducing (winter) and infected Haliotis asinina (a vetigastropod mollusc). A number of the regulatory genes are differentially expressed in parasitised H. asinina, but in only a few cases do expression patterns in infected animals match those occurring in animals where reproduction is normally repressed. (c) 2006 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.

The potential role of environmental exposures and genomic signaling in development of central nervous system tumors

The etiology of central nervous system tumors (CNSTs) is mainly unknown. Aside from extremely rare genetic conditions, such as neurofibromatosis and tuberous sclerosis, the only unequivocally identified risk factor is exposure to ionizing radiation, and this explains only a very small fraction of cases. Using meta-analysis, gene networking and bioinformatics methods, this dissertation explored the hypothesis that environmental exposures produce genetic and epigenetic alterations that may be involved in the etiology of CNSTs. A meta-analysis of epidemiological studies of pesticides and pediatric brain tumors revealed a significantly increased risk of brain tumors among children whose mothers had farm-related exposures during pregnancy. A dose response was recognized when this risk estimate was compared to those for risk of brain tumors from maternal exposure to non-agricultural pesticides during pregnancy, and risk of brain tumors among children exposed to agricultural activities. Through meta-analysis of several microarray studies which compared normal tissue to astrocytomas, we were able to identify a list of 554 genes which were differentially expressed in the majority of astrocytomas. Many of these genes have in fact been implicated in development of astrocytoma, including EGFR, HIF-1α, c-Myc, WNT5A, and IDH3A. Reverse engineering of these 554 genes using Bayesian network analysis produced a gene network for each grade of astrocytoma (Grade I-IV), and ‘key genes’ within each grade were identified. Genes found to be most influential to development of the highest grade of astrocytoma, Glioblastoma multiforme (GBM) were: COL4A1, EGFR, BTF3, MPP2, RAB31, CDK4, CD99, ANXA2, TOP2A, and SERBP1. Lastly, bioinformatics analysis of environmental databases and curated published results on GBM was able to identify numerous potential pathways and geneenvironment interactions that may play key roles in astrocytoma development. Findings from this research have strong potential to advance our understanding of the etiology and susceptibility to CNSTs. Validation of our ‘key genes’ and pathways could potentially lead to useful tools for early detection and novel therapeutic options for these tumors.

High natural gene expression variation in the reef-building coral Acropora millepora: potential for acclimative and adaptive plasticity

Background: Ecosystems worldwide are suffering the consequences of anthropogenic impact. The diverse ecosystem of coral reefs, for example, are globally threatened by increases in sea surface temperatures due to global warming. Studies to date have focused on determining genetic diversity, the sequence variability of genes in a species, as a proxy to estimate and predict the potential adaptive response of coral populations to environmental changes linked to climate changes. However, the examination of natural gene expression variation has received less attention. This variation has been implicated as an important factor in evolutionary processes, upon which natural selection can act. Results: We acclimatized coral nubbins from six colonies of the reef-building coral Acropora millepora to a common garden in Heron Island (Great Barrier Reef, GBR) for a period of four weeks to remove any site-specific environmental effects on the physiology of the coral nubbins. By using a cDNA microarray platform, we detected a high level of gene expression variation, with 17% (488) of the unigenes differentially expressed across coral nubbins of the six colonies (jsFDR-corrected, p < 0.01). Among the main categories of biological processes found differentially expressed were transport, translation, response to stimulus, oxidation-reduction processes, and apoptosis. We found that the transcriptional profiles did not correspond to the genotype of the colony characterized using either an intron of the carbonic anhydrase gene or microsatellite loci markers. Conclusion: Our results provide evidence of the high inter-colony variation in A. millepora at the transcriptomic level grown under a common garden and without a correspondence with genotypic identity. This finding brings to our attention the importance of taking into account natural variation between reef corals when assessing experimental gene expression differences. The high transcriptional variation detected in this study is interpreted and discussed within the context of adaptive potential and phenotypic plasticity of reef corals. Whether this variation will allow coral reefs to survive to current challenges remains unknown.