Low but contrasting neutral genetic differentiation shaped by winter temperature in European great tits

Gene flow is usually thought to reduce genetic divergence and impede local adaptation by homogenising gene pools between populations. However, evidence for local adaptation and phenotypic differentiation in highly mobile species, experiencing high levels of gene flow, is emerging. Assessing population genetic structure at different spatial scales is thus a crucial step towards understanding mechanisms underlying intraspecific differentiation and diversification. Here, we studied the population genetic structure of a highly mobile species – the great tit Parus major – at different spatial scales. We analysed 884 individuals from 30 sites across Europe including 10 close-by sites (< 50 km), using 22 microsatellite markers. Overall we found a low but significant genetic differentiation among sites (FST = 0.008). Genetic differentiation was higher, and genetic diversity lower, in south-western Europe. These regional differences were statistically best explained by winter temperature. Overall, our results suggest that great tits form a single patchy metapopulation across Europe, in which genetic differentiation is independent of geographical distance and gene flow may be regulated by environmental factors via movements related to winter severity. This might have important implications for the evolutionary trajectories of sub-populations, especially in the context of climate change, and calls for future investigations of local differences in costs and benefits of philopatry at large scales.

Contemporary genetic, historical genetic and morphological data sets for the Yellow-tufted Honeyeater Lichenostomus melanops

Understanding the evolutionary history of threatened populations can improve their conservation management. Re-establishment of past but recent gene flow could re-invigorate threatened populations and replenish genetic diversity, necessary for population persistence. One of the four nominal subspecies of the common yellow-tufted honeyeater, Lichenostomus melanops cassidix, is critically endangered despite substantial conservation efforts over 55 years. Using a combination of morphometric, genetic and modelling approaches we tested for its evolutionary distinctiveness and conservation merit. We confirmed that cassidix has at least one morphometric distinction. It also differs genetically from the other subspecies in allele frequencies but not phylogenetically, implying that its evolution was recent. Modelling historical distribution supported the lack of vicariance and suggested a possibility of gene flow among subspecies at least since the late Pleistocene. Multi-locus coalescent analyses indicated that cassidix diverged from its common ancestor with neighbouring subspecies gippslandicus sometime from the mid-Pleistocene to the Holocene, and that it has the smallest historical effective population size of all subspecies. It appears that cassidix diverged from its ancestor with gippslandicus through a combination of drift and local selection. From patterns of genetic subdivision on two spatial scales and morphological variation we concluded that cassidix, gippslandicus and (melanops + meltoni) are diagnosable as subspecies. Low genetic diversity and effective population size of cassidix may translate to low genetic fitness and evolutionary potential, thus managed gene flow from gippslandicus is recommended for its recovery.

Combined use of gliadins and SSRs to analyse the genetic variability of the Spanish collection of cultivated diploid wheat (Triticum monococcum L. ssp. monococcum)

This work studied the combined use of gliadins and SSRs to analyse inter- and intra-accession variability of the Spanish collection of cultivated einkorn (Triticum monococcum L. ssp. monococcum) maintained at the CRF-INIA. In general, gliadin loci presented higher discrimination power than SSRs, reflecting the high variability of the gliadins. The loci on chromosome 6A were the most polymorphic with similar PIC values for both marker systems, showing that these markers are very useful for genetic variability studies in wheat. The gliadin results indicated that the Spanish einkorn collection possessed high genetic diversity, being the differentiation large between varieties and small within them. Some associations between gliadin alleles and geographical and agro-morphological data were found. Agro-morphological relations were also observed in the clusters of the SSRs dendrogram. A high concordance was found between gliadins and SSRs for genotype identification. In addition, both systems provide complementary information to resolve the different cases of intra-accession variability not detected at the agro-morphological level, and to identify separately all the genotypes analysed. The combined use of both genetic markers is an excellent tool for genetic resource evaluation in addition to agro-morphological evaluation.

Analysis of genetic variability in a sample of the durum wheat (Triticum durum Desf.) Spanish collection based on gliadin markers

In this work gliadin proteins were used to analyse the genetic variability in a sample of the durum wheat Spanish collection conserved at the CRF-INIA. In total 38 different alleles were identified at the loci Gli-A1, Gli-A3, Gli-B5, Gli-B1, Gli-A2 and Gli-B2. All the gliadin loci were polymorphic, possessed large genetic diversity and small and large differentiation within and between varieties, respectively. The Gli-A2 and Gli-B2 loci were the most polymorphic, the most fixed within varieties and the most useful to distinguish among varieties. Alternatively, Gli-B1 locus presented the least genetic variability out of the four main loci Gli-A1, Gli-B1, Gli-A2 and Gli-B2. The Gli-B1 alleles coding for the gliadin γ-45, associated with good quality, had an accumulated frequency of 69.7%, showing that the Spanish germplasm could be a good source for breeding quality. The Spanish landraces studied showed new gliadin alleles not catalogued so far. These new alleles might be associated with specific Spanish environment factors. The large number of new alleles identified also indicates that durum wheat Spanish germplasm is rather unique.

Molecular tools to improve chestnut management: El Bierzo as a case study

The European chestnut (Castanea sativa Mill.) is a multipurpose species that has been widely cultivated around the Mediterranean basin since ancient times. New varieties were brought to the Iberian Peninsula during the Roman Empire, which coexist since then with native populations that survived the last glaciation. The relevance of chestnut cultivation has being steadily growing since the Middle Ages, until the rural decline of the past century put a stop to this trend. Forest fires and diseases were also major factors. Chestnut cultivation is gaining momentum again due to its economic (wood, fruits) and ecologic relevance, and represents currently an important asset in many rural areas of Europe. In this Thesis we apply different molecular tools to help improve current management strategies. For this study we have chosen El Bierzo (Castile and Leon, NW Spain), which has a centenary tradition of chestnut cultivation and management, and also presents several unique features from a genetic perspective (next paragraph). Moreover, its nuts are widely appreciated in Spain and abroad for their organoleptic properties. We have focused our experimental work on two major problems faced by breeders and the industry: the lack of a fine-grained genetic characterization and the need for new strategies to control blight disease. To characterize with sufficient detail the genetic diversity and structure of El Bierzo orchards, we analyzed DNA from 169 trees grafted for nut production covering the entire region. We also analyzed 62 nuts from all traditional varieties. El Bierzo constitutes an outstanding scenario to study chestnut genetics and the influence of human management because: (i) it is located at one extreme of the distribution area; (ii) it is a major glacial refuge for the native species; (iii) it has a long tradition of human management (since Roman times, at least); and (iv) its geographical setting ensures an unusual degree of genetic isolation. Thirteen microsatellite markers provided enough informativeness and discrimination power to genotype at the individual level. Together with an unexpected level of genetic variability, we found evidence of genetic structure, with three major gene pools giving rise to the current population. High levels of genetic differentiation between groups supported this organization. Interestingly, genetic structure does not match with spatial boundaries, suggesting that the exchange of material and cultivation practices have strongly influenced natural gene flow. The microsatellite markers selected for this study were also used to classify a set of 62 samples belonging to all traditional varieties. We identified several cases of synonymies and homonymies, evidencing the need to substitute traditional classification systems with new tools for genetic profiling. Management and conservation strategies should also benefit from these tools. The avenue of high-throughput sequencing technologies, combined with the development of bioinformatics tools, have paved the way to study transcriptomes without the need for a reference genome. We took advantage of RNA sequencing and de novo assembly tools to determine the transcriptional landscape of chestnut in response to blight disease. In addition, we have selected a set of candidate genes with high potential for developing resistant varieties via genetic engineering. Our results evidenced a deep transcriptional reprogramming upon fungal infection. The plant hormones ET and JA appear to orchestrate the defensive response. Interestingly, our results also suggest a role for auxins in modulating such response. Many transcription factors were identified in this work that interact with promoters of genes involved in disease resistance. Among these genes, we have conducted a functional characterization of a two major thaumatin-like proteins (TLP) that belongs to the PR5 family. Two genes encoding chestnut cotyledon TLPs have been previously characterized, termed CsTL1 and CsTL2. We substantiate here their protective role against blight disease for the first time, including in silico, in vitro and in vivo evidence. The synergy between TLPs and other antifungal proteins, particularly endo-p-1,3-glucanases, bolsters their interest for future control strategies based on biotechnological approaches.

Caracterización de poblaciones de vid silvestre de la Península Ibérica

La vid silvestre se considera como el ancestro autóctono de las vides cultivadas y una enorme reserva genética en peligro de extinción. La prospección llevada a cabo entre 2003 y 2004 permitió catalogar 51 localizaciones de vides silvestres españolas, la mayoría de ellas ubicadas en riberas de ríos. Estos ejemplares se incluyeron en el Banco de Germoplasma de la Finca "El Encín" (BGVCAM - Alcalá de Henares, Madrid, España). En primer lugar, se caracterizó la cantidad y la distribución de su diversidad genética utilizando 25 loci empleando microsatélites nucleares (SSR). Hemos analizado también la posible coexistencia en el hábitat natural de vides silvestres con vides cultivadas naturalizadas y portainjertos. De este modo, los análisis fenotípicos y genéticos identificaron el 19% de las muestras recogidas como derivadas de genotipos cultivados, siendo, o bien vides cultivadas naturalizadas o genotipos híbridos derivados de cruces espontáneos entre vides silvestres y cultivadas. La diversidad genética de las poblaciones de vides silvestres fue similar a la observada en el grupo de las cultivadas. El análisis molecular mostró que el germoplasma de cultivadas y silvestres es genéticamente divergente con bajo nivel de introgresión. Se ha identificado cuatro grupos genéticos, con dos de ellos fundamentalmente representados por los genotipos de vides cultivadas y dos por las accesiones silvestres. El análisis de los vínculos genéticos entre las vides silvestres y cultivadas podría sugerir una contribución genética de las accesiones silvestres españolas a las actuales variedades occidentales. En segundo lugar, se realizó un profundo estudio morfológico "ex situ " y se contrastaron con los resultados de la caracterización realizada en 182 variedades comerciales españolas de la misma colección. Todos los individuos silvestres mostraron diferencias morfológicas con Vitis vinifera L subsp. vinifera, pero no se encontraron diferencias significativas dentro Vitis vinifera L. subsp. sylvestris, ni por localización geográfica ni por sexo. Los resultados de este estudio describen las principales características morfológicas de las vides silvestres españolas y sus rasgos diferenciales con su pariente cultivada. Por último, se analizó la composición antociánica presente en 21 accesiones de vides silvestres de la Península Ibérica conservadas en el BGVCAM de la Finca "El Encín" y seleccionadas basándose en diferencias ampelográficas y caracterización molecular. La concentración de antocianinas es similar la encontrad en vides cultivadas con destino a la vinificación. Las accesiones estudiadas mostraron una variabilidad considerable en su perfil antociánico y fue posible distinguir varios grupos. Sin embargo, la presencia de material silvestre con perfiles antociánicos poco comunes o inexistentes en variedades españolas, sugiere que la variabilidad genética relacionada con antocianinas en poblaciones españolas de vides silvestres podría ser más alta que la de variedades cultivadas comúnmente consideradas de origen español. ABSTRACT The wild grapevine is considered an autochthonous relative of cultivated vines and a huge gene pool endangered in Europe. Prospecting carried out between 2003 and 2004 enabled to inventory 51 Spanish sites with wild grapevines, most of them located near rivers. These individuals were grafted in the collection of "El Encín" (BGVCAM - Alcalá de Henares, Madrid, Spain). Firstly, werw characterized the amount and distribution of their genetic diversity using 25 nuclear SSR loci. We have also analysed the possible coexistence in the natural habitat of wild grapevines with naturalized grapevine cultivars and rootstocks. In this way, phenotypic and genetic analyses identified 19% of the collected samples as derived from cultivated genotypes, being either naturalized cultivars or hybrid genotypes derived from spontaneous crosses between wild and cultivated grapevines. The genetic diversity of wild grapevine populations was similar than that observed in the cultivated group. The molecular analysis showed that cultivated germplasm and wild germplasm are genetically divergent with low level of introgression. We identified four genetic groups, with two of them fundamentally represented among cultivated genotypes and two among wild accessions. The analyses of genetic relationships between wild and cultivated grapevines could suggest a genetic contribution of wild accessions from Spain to current Western cultivars. Secondly, a morphological study was done "ex situ" and were compared with data from 182 Spanish commercial cultivars grown in the same collection. All wild individuals showed morphological differences with Vitis vinifera L. ssp. vinifera but no significant differences were found within Vitis vinifera L subsp. sylvestris neither by geographic origin nor by sex. A pattern with the main characteristics of Spanish wild grapevines is suggested. Ultimately, were investigated the anthocyanin composition of 21 mostly Spanish wild grapevine accessions preserved at BGVCAM "El Encín" and selected in consideration of observed ampelographic differences and molecular characterization. Total anthocyanin concentration was similar to that found in winegrape cultivars. The accessions studied showed considerable variability in their anthocyanin fingerprints and it was possible to distinguish several groups, similar to previous reports on the anthocyanin fingerprint of winegrapes. The anthocyanin composition of wild grapevine accessions was similar to that of cultivated grapes. Nevertheless, the presence of wild accessions with anthocyanin fingerprints uncommon or nonexistent in Spanish cultivated varieties suggests that the genetic variability related to anthocyanins in Spanish wild grapevine populations may be higher than that of cultivated varieties commonly considered of Spanish origin.

Complete congruence between gene diversity estimates derived from genotypic data at enzyme and random amplified polymorphic DNA loci in black spruce.

Controversy still exists over the adaptive nature of variation of enzyme loci. In conifers, random amplified polymorphic DNAs (RAPDs) represent a class of marker loci that is unlikely to fall within or be strongly linked to coding DNA. We have compared the genetic diversity in natural populations of black spruce [Picea mariana (Mill.) B.S.P.] using genotypic data at allozyme loci and RAPD loci as well as phenotypic data from inferred RAPD fingerprints. The genotypic data for both allozymes and RAPDs were obtained from at least six haploid megagametophytes for each of 75 sexually mature individuals distributed in five populations. Heterozygosities and population fixation indices were in complete agreement between allozyme loci and RAPD loci. In black spruce, it is more likely that the similar levels of variation detected at both enzyme and RAPD loci are due to such evolutionary forces as migration and the mating system, rather than to balancing selection and overdominance. Furthermore, we show that biased estimates of expected heterozygosity and among-population differentiation are obtained when using allele frequencies derived from dominant RAPD phenotypes.

Variabilidade populacional em manguezais: análises moleculares e morfológicas em caranguejos Brachyura (Crustacea: Decapoda)

A análise da estruturação populacional de espécies codistribuídas permite a comparação de padrões de estruturação, fornecendo informações acerca dos fatores que influenciam a diferenciação em espécies pertencentes ao mesmo ecossistema. Este projeto teve como objetivo principal analisar a variabilidade intraespecífica em caranguejos habitantes de manguezais, codistribuídos ao longo do Oceano Atlântico Ocidental, por meio de ferramentas moleculares e morfológicas, visando testar a hipótese de elevada estruturação populacional em manguezais. Para este fim, cinco espécies foram utilizadas como modelo (Aratus pisonii, Goniopsis cruentata, Sesarma rectum, Uca thayeri e Ucides cordatus) e avaliadas por meio de marcadores mitocondriais COI e 16S e nuclear H3 e análises morfológicas comparativas e de morfometria. Os dados moleculares revelaram dois padrões, indicando elevada estruturação populacional para as espécies A. pisonii e U. thayeri e ausência de estruturação para G. cruentata, S. rectum e U. cordatus. Os dados morfológicos, no entanto, não acompanham esses padrões, já que não foram encontradas diferenças morfológicas ou morfométricas associadas aos grupos evidenciados pelas análises moleculares. A ausência de fluxo gênico entre regiões para algumas espécies deve-se, muito provavelmente, à existência de fatores que não se limitam ao isolamento por distância, mas também devido a diferenças na duração do estágio larval e a diferenças bruscas em alguns fatores abióticos, como a salinidade, por exemplo, que, associados às diferentes características do desenvolvimento larval de cada espécie, culminam na existência de estruturas populacionais diferentes. Além disso, os padrões de diferenciação genética observados concordam com os cenários biogeográficos propostos para o Atlântico Ocidental, no qual mudanças e flutuações geológicas, climáticas e oceanográficas, resultantes do fechamento do Istmo do Panamá e de ciclos glaciais na América do Norte, promoveram divergência genética.

Genetic effects of chronic habitat fragmentation on tree species: the case of Sorbus aucuparia in a deforested Scottish landscape

Sustainable forest restoration and management practices require a thorough understanding of the influence that habitat fragmentation has on the processes shaping genetic variation and its distribution in tree populations. We quantified genetic variation at isozyme markers and chloroplast DNA (cpDNA), analysed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in severely fragmented populations of Sorbus aucuparia (Rosaceae) in a single catchment (Moffat) in southern Scotland. Remnants maintain surprisingly high levels of gene diversity (H-E) for isozymes (H-E = 0.195) and cpDNA markers (H-E = 0.490). Estimates are very similar to those from non-fragmented populations in continental Europe, even though the latter were sampled over a much larger spatial scale. Overall, no genetic bottleneck or departures from random mating were detected in the Moffat fragments. However, genetic differentiation among remnants was detected for both types of marker (isozymes Theta(n) = 0.043, cpDNA Theta(c) = 0.131; G-test, P-value < 0.001). In this self-incompatible, insect-pollinated, bird-dispersed tree species, the estimated ratio of pollen flow to seed flow between fragments is close to 1 (r = 1.36). Reduced pollen-mediated gene flow is a likely consequence of habitat fragmentation, but effective seed dispersal by birds is probably helping to maintain high levels of genetic diversity within remnants and reduce genetic differentiation between them.

Limited effect of anthropogenic habitat fragmentation on molecular diversity in a rain forest skink, Gnypetoscincus queenslandiae

To examine the effects of recent habitat fragmentation, we assayed genetic diversity in a rain forest endemic lizard, the prickly forest skink (Gnypetoscincus queenslandiae), from seven forest fragments and five sites in continuous forest on the Atherton tableland of northeastern Queensland, Australia. The rain forest in this region was fragmented by logging and clearing for dairy farms in the early 1900s and most forest fragments studied have been isolated for 50-80 years or nine to 12 skink generations. We genotyped 411 individuals at nine microsatellite DNA loci and found fewer alleles per locus in prickly forest skinks from small rain forest fragments and a lower ratio of allele number to allele size range in forest fragments than in continuous forest, indicative of a decrease in effective population size. In contrast, and as expected for populations with small neighbourhood sizes, neither heterozygosity nor variance in allele size differed between fragments and sites in continuous forests. Considering measures of among population differentiation, there was no increase in F-ST among fragments and a significant isolation by distance pattern was identified across all 12 sites. However, the relationship between genetic (F-ST) and geographical distance was significantly stronger for continuous forest sites than for fragments, consistent with disruption of gene flow among the latter. The observed changes in genetic diversity within and among populations are small, but in the direction predicted by the theory of genetic erosion in recently fragmented populations. The results also illustrate the inherent difficulty in detecting genetic consequences of recent habitat fragmentation, even in genetically variable species, and especially when effective population size and dispersal rates are low.

Optimal sampling strategy for estimation of spatial genetic structure in tree populations

Fine-scale spatial genetic structure (SGS) in natural tree populations is largely a result of restricted pollen and seed dispersal. Understanding the link between limitations to dispersal in gene vectors and SGS is of key interest to biologists and the availability of highly variable molecular markers has facilitated fine-scale analysis of populations. However, estimation of SGS may depend strongly on the type of genetic marker and sampling strategy (of both loci and individuals). To explore sampling limits, we created a model population with simulated distributions of dominant and codominant alleles, resulting from natural regeneration with restricted gene flow. SGS estimates from subsamples (simulating collection and analysis with amplified fragment length polymorphism (AFLP) and microsatellite markers) were correlated with the 'real' estimate (from the full model population). For both marker types, sampling ranges were evident, with lower limits below which estimation was poorly correlated and upper limits above which sampling became inefficient. Lower limits (correlation of 0.9) were 100 individuals, 10 loci for microsatellites and 150 individuals, 100 loci for AFLPs. Upper limits were 200 individuals, five loci for microsatellites and 200 individuals, 100 loci for AFLPs. The limits indicated by simulation were compared with data sets from real species. Instances where sampling effort had been either insufficient or inefficient were identified. The model results should form practical boundaries for studies aiming to detect SGS. However, greater sample sizes will be required in cases where SGS is weaker than for our simulated population, for example, in species with effective pollen/seed dispersal mechanisms.

Contrasting quantitative traits and neutral genetic markers for genetic resource assessment of Mesoamerican Cedrela odorata

We compared within-population variability and degree of population differentiation for neutral genetic markers (RAPDS) and eight quantitative traits in Central American populations of the endangered tree, Cedrela odorata. Whilst population genetic diversity for neutral markers (Shannon index) and quantitative traits (heritability, coefficient of additive genetic variation) were uncorrelated, both marker types revealed strong differentiation between populations from the Atlantic coast of Costa Rica and the rest of the species' distribution. The degree of interpopulation differentiation was higher for RAPD markers (F-ST 0.67 for the sampled Mesoamerican range) than for quantitative traits (Q(ST) = 0.30). Hence, the divergence in quantitative traits was lower than could have been achieved by genetic drift alone, suggesting that balancing selection for similar phenotypes in different populations of this species. Nevertheless, a comparison of pair-wise estimates of population differentiation in neutral genetic markers and quantitative traits revealed a strong positive correlation (r = 0.66) suggesting that, for C. odorata, neutral marker divergence could be used as a surrogate for adaptive gene divergence for conservation planning. The utility of this finding and suggested further work are discussed.

Genetic structure of a reef-building coral from thermally distinct environments on the Great Barrier Reef

Adaptation to localised thermal regimes is facilitated by restricted gene flow, ultimately leading to genetic divergence among populations and differences in their physiological tolerances. Allozyme analysis of six polymorphic loci was used to assess genetic differentiation between nine populations of the reef-building coral Acropora millepora over a latitudinal temperature gradient on the inshore regions of the Great Barrier Reef (GBR). Small but significant genetic differentiation indicative of moderate levels of gene flow (pairwise F-ST 0.023 to 0.077) was found between southern populations of A. millepora in cooler regions of the GBR and the warmer, central or northern GBR populations. Patterns of genetic differentiation at these putatively neutral allozyme loci broadly matched experimental variation in thermal tolerance and were consistent with local thermal regimes (warmest monthly-averages) for the A. millepora populations examined. It is therefore hypothesized that natural selection has influenced the thermal tolerance of the A. millepora populations examined and greater genetic divergence is likely to be revealed by examination of genetic markers under the direct effects of natural selection.

Spatial and temporal variation in algal-assemblage structure in isolated dryland river waterholes, Cooper Creek and Warrego River, Australia

The scale at which algal biodiversity is partitioned across the landscape, and the biophysical processes and biotic interactions which shape these communities in dryland river refugia was studied on two occasions from 30 sites in two Australian dryland rivers. Despite the waterholes studied having characteristically high levels of abiogenic turbidity, a total of 186 planktonic microalgae, 253 benthic diatom and 62 macroalgal species were recorded. The phytoplankton communities were dominated by flagellated cryptophytes, euglenophytes and chlorophytes, the diatom communities by cosmopolitan taxa known to tolerate wide environmental conditions, and the macroalgal communities by filamentous cyanobacteria. All algal communities showed significant differences between catchments and sampling times, with a suite of between 5 and 12 taxa responsible for similar to 50% of the observed change. In general, algal assemblage patterns were poorly correlated with the measured environmental variables. Phytoplankton and diatom assemblage patterns were weakly correlated with several waterhole geomorphic measures, whereas macroalgal assemblage patterns showed some association with variability in ionic concentration.

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.