Analysis of the transcriptome of the needles and bark of Pinus radiata induced by bark stripping and methyl jasmonate.

BACKGROUND: Plants are attacked by diverse insect and mammalian herbivores and respond with different physical and chemical defences. Transcriptional changes underlie these phenotypic changes. Simulated herbivory has been used to study the transcriptional and other early regulation events of these plant responses. In this study, constitutive and induced transcriptional responses to artificial bark stripping are compared in the needles and the bark of Pinus radiata to the responses from application of the plant stressor, methyl jasmonate. The time progression of the responses was assessed over a 4-week period. RESULTS: Of the 6312 unique transcripts studied, 86.6% were differentially expressed between the needles and the bark prior to treatment. The most abundant constitutive transcripts were related to defence and photosynthesis and their expression did not differ between the needles and the bark. While no differential expression of transcripts were detected in the needles following bark stripping, in the bark this treatment caused an up-regulation and down-regulation of genes associated with primary and secondary metabolism. Methyl jasmonate treatment caused differential expression of transcripts in both the bark and the needles, with individual genes related to primary metabolism more responsive than those associated with secondary metabolism. The up-regulation of genes related to sugar break-down and the repression of genes related with photosynthesis, following both treatments was consistent with the strong down-regulation of sugars that has been observed in the same population. Relative to the control, the treatments caused a differential expression of genes involved in signalling, photosynthesis, carbohydrate and lipid metabolism as well as defence and water stress. However, non-overlapping transcripts were detected between the needles and the bark, between treatments and at different times of assessment. Methyl jasmonate induced more transcriptional responses in the bark than bark stripping, although the peak of expression following both treatments was detected 7 days post treatment application. The effects of bark stripping were localised, and no systemic changes were detected in the needles. CONCLUSION: There are constitutive and induced differences in the needle and bark transcriptome of Pinus radiata. Some expression responses to bark stripping may differ from other biotic and abiotic stresses, which contributes to the understanding of plant molecular responses to diverse stresses. Whether the gene expression changes are heritable and how they differ between resistant and susceptible families identified in earlier studies needs further investigation.

Expansion of the rare Eucalyptus risdonii under climate change through hybridization with a closely related species despite hybrid inferiority.

BACKGROUND AND AIMS: Hybridization is increasingly recognized as an integral part of the dynamics of species range expansion and contraction. Thus, it is important to understand the reproductive barriers between co-occurring species. Extending previous studies that argued that the rare Eucalyptus risdonii was expanding into the range of the surrounding E. amygdalina by both seed and pollen dispersal, we here investigate the long-term fitness of both species and their hybrids and whether expansion is continuing. METHODS: We assessed the survival of phenotypes representing a continuum between the two pure species in a natural hybrid swarm after 29 years, along with seedling recruitment. The performance of pure species as well as of artificial and natural hybrids was also assessed over 28 years in a common garden trial. KEY RESULTS: In the hybrid zone, E. amygdalina adults showed greater mortality than E. risdonii, and the current seedling cohort is still dominated by E. risdonii phenotypes. Morphologically intermediate individuals appeared to be the least fit. Similar results were observed after growing artificial first-generation and natural hybrids alongside pure species families in a common garden trial. Here, the survival, reproduction, health and growth of the intermediate hybrids were significantly less than those of either pure species, consistent with hybrid inferiority, although this did not manifest until later reproductive ages. Among the variable progeny of natural intermediate hybrids, the most E. risdonii-like phenotypes were the most fit. CONCLUSIONS: This study contributes to the increasing number of reports of hybrid inferiority in Eucalyptus, suggesting that post-zygotic barriers contribute to the maintenance of species integrity even between closely related species. However, with fitness rapidly recovered following backcrossing, it is argued that hybridization can still be an important evolutionary process, in the present case appearing to contribute to the range expansion of the rare E. risdonii in response to climate change.

Genetic-based interactions among tree neighbors: identification of the most influential neighbors, and estimation of correlations among direct and indirect genetic effects for leaf disease and growth in Eucalyptus globulus.

An individual's genes may influence the phenotype of neighboring conspecifics. Such indirect genetic effects (IGEs) are important as they can affect the apparent total heritable variance in a population, and thus the response to selection. We studied these effects in a large, pedigreed population of Eucalyptus globulus using variance component analyses of Mycosphearella leaf disease, diameter growth at age 2 years, and post-infection diameter growth at ages 4 and 8 years. In a novel approach, we initially modeled IGEs using a factor analytic (FA) structure to identify the most influential neighbor positions, with the FA loadings being position-specific regressions on the IGEs. This involved sequentially comparing FA models for the variance-covariance matrices of the direct and indirect effects of each neighbor. We then modeled IGEs as a distance-based, combined effect of the most influential neighbors. This often increased the magnitude and significance of indirect genetic variance estimates relative to using all neighbors. The extension of a univariate IGEs model to bivariate analyses also provided insights into the genetic architecture of this population, revealing that: (1) IGEs arising from increased probability of neighbor infection were not associated with reduced growth of neighbors, despite adverse fitness effects being evident at the direct genetic level; and (2) the strong, genetic-based competitive interactions for growth, established early in stand development, were highly positively correlated over time. Our results highlight the complexities of genetic-based interactions at the multi-trait level due to (co)variances associated with IGEs, and the marked discrepancy occurring between direct and total heritable variances.

Comparative genomics of Eucalyptus and Corymbia reveals low rates of genome structural rearrangement.

BACKGROUND: Previous studies suggest genome structure is largely conserved between Eucalyptus species. However, it is unknown if this conservation extends to more divergent eucalypt taxa. We performed comparative genomics between the eucalypt genera Eucalyptus and Corymbia. Our results will facilitate transfer of genomic information between these important taxa and provide further insights into the rate of structural change in tree genomes. RESULTS: We constructed three high density linkage maps for two Corymbia species (Corymbia citriodora subsp. variegata and Corymbia torelliana) which were used to compare genome structure between both species and Eucalyptus grandis. Genome structure was highly conserved between the Corymbia species. However, the comparison of Corymbia and E. grandis suggests large (from 1-13 MB) intra-chromosomal rearrangements have occurred on seven of the 11 chromosomes. Most rearrangements were supported through comparisons of the three independent Corymbia maps to the E. grandis genome sequence, and to other independently constructed Eucalyptus linkage maps. CONCLUSIONS: These are the first large scale chromosomal rearrangements discovered between eucalypts. Nonetheless, in the general context of plants, the genomic structure of the two genera was remarkably conserved; adding to a growing body of evidence that conservation of genome structure is common amongst woody angiosperms.

A latitudinal cline in disease resistance of a host tree.

The possible drivers and implications of an observed latitudinal cline in disease resistance of a host tree were examined. Mycosphaerella leaf disease (MLD) damage, caused by Teratosphaeria species, was assessed in five Eucalyptus globulus (Tasmanian blue gum) common garden trials containing open-pollinated progeny from 13 native-forest populations. Significant population and family within population variation in MLD resistance was detected, which was relatively stable across different combinations of trial sites, ages, seasons and epidemics. A distinct genetic-based latitudinal cline in MLD damage among host populations was evident. Two lines of evidence argue that the observed genetic-based latitudinal trend was the result of direct pathogen-imposed selection for MLD resistance. First, MLD damage was positively associated with temperature and negatively associated with a prediction of disease risk in the native environment of these populations; and, second, the quantitative inbreeding coefficient (QST) significantly exceeded neutral marker FST at the trial that exhibited the greatest MLD damage, suggesting that diversifying selection contributed to differentiation in MLD resistance among populations. This study highlights the potential for spatial variation in pathogen risk to drive adaptive differentiation across the geographic range of a foundation host tree species.

The effects of age and environment on the expression of inbreeding depression in Eucalyptus globulus.

Inbreeding adversely affects fitness traits in many plant and animal species, and the magnitude, stability and genetic basis of inbreeding depression (ID) will have short- and long-term evolutionary consequences. The effects of four degrees of inbreeding (selfing, f=50%; full- and half-sib matings, f=25 and 12.5%; and unrelated outcrosses, f=0%) on survival and growth of an island population of Eucalyptus globulus were studied at two sites for over 14 years. For selfs, ID in survival increased over time, reaching a maximum of 49% by age 14 years. However, their inbreeding depression for stem diameter remained relatively stable with age, and ranged from 28 to 36% across years and sites. ID for survival was markedly greater on the more productive site, possibly due to greater and earlier onset of inter-tree competition, but was similar on both sites for the diameter of survivors. The deleterious trait response to increasing inbreeding coefficients was linear for survival and diameter. Non-significant quadratic effects suggested that epistasis did not contribute considerably to the observed ID at the population level. Among- and within-family coefficients of variation for diameter increased with inbreeding degree, and the variance among the outcrossed families was significant only on the more productive site. The performance of self-families for diameter was highly stable between sites. This suggests that, for species with mixed mating systems, environmentally stable inbreeding effects in open-pollinated progenies may tend to mask the additive genotype-by-environment interaction for fitness traits and the adaptive response to the environment.

Recurrent nuclear DNA introgression accompanies chloroplast DNA exchange between two eucalypt species.

Numerous studies within plant genera have found geographically structured sharing of chloroplast (cp) DNA among sympatric species, consistent with introgressive hybridization. Current research is aimed at understanding the extent, direction and significance of nuclear (nr) DNA exchange that accompanies putative cpDNA exchange. Eucalyptus is a complex tree genus for which cpDNA sharing has been established between multiple species. Prior phylogeographic analysis has indicated cpDNA introgression into the widespread forest species Eucalyptus globulus from its rare congener E. cordata. In this study, we use AFLP markers to characterize corresponding nrDNA introgression, on both a broad and fine spatial scale. Using 388 samples we examine (i) the fine-scale spatial structure of cp and nrDNA introgression from E. cordata into E. globulus at a site in natural forest and (ii) broad-scale patterns of AFLP marker introgression at six additional mixed populations. We show that while E. globulus and E. cordata retain strongly differentiated nuclear gene pools overall, leakage of nrDNA occurs at mixed populations, with some AFLP markers being transferred to E. globulus recurrently at different sites. On the fine scale, different AFLP fragments show varying distances of introgression into E. globulus, while introgression of cpDNA is extensive. The frequency of E. cordata markers in E. globulus is correlated with spatial proximity to E. cordata, but departs from expectations based on AFLP marker frequency in E. cordata, indicating that selection may be governing the persistence of introgressed fragments in E. globulus.

The relationship of the female reproductive success of Eucalyptus globulus to the endogenous properties of the flower.

Low capsule and seed set is a major factor limiting seed production in Eucalyptus globulus seed orchards. Controlled pollination studies showed that the reproductive success (number of seeds produced per flower pollinated) was primarily determined by the female. We aimed to identify the factors contributing to the differences in reproductive success between female genotypes in terms of the physical and anatomical properties of the flower. We studied pairs of genotypes of high and low reproductive success from each of three races (Furneaux Group, Strzelecki Ranges and Western Otways) growing in a seed orchard. Controlled pollinations were performed on six females and along with flower physical measurements, pollen tube growth and seed set were assessed. Overall tree reproductive success was positively correlated with flower size, ovule numbers, style size, cross-sectional area of conductive tissue within the style (all of which were inter-correlated) and the proportion of pollen tubes reaching the bottom of the style. Significant positive correlations of reproductive success and flower physical properties between different ramets of the same genotypes across seasons suggests a genetic basis to the variation observed. The majority of pollen tube attrition occurred within the first millimetre of the cut style and appeared to be associated with differences in style physiology. When examined as pairs within races the difference in reproductive success for the Western Otways pair was simply explained by differences in flower size and the number of ovules per flower. Physical features did not differ significantly for the Strzelecki Ranges pair, but the proportion of pollen tubes reaching the bottom of the style was lower in the less reproductively successful genotype, suggesting an endogenous physiological constraint to pollen tube growth. The difference in reproductive success between the females from the Furneaux Group was associated with a combination of these factors.

Quantitative trait loci for key defensive compounds affecting herbivory of eucalypts in Australia.

* Formylated phloroglucinols (FPCs) are key defensive compounds that influence herbivory by mammals and arthropods in eucalypts. However, the genetic architecture underlying variation in their levels remains poorly understood. * Quantitative trait loci (QTL) analysis for the concentrations of two major FPCs, sideroxylonal A and macrocarpal G, was conducted using juvenile leaves from 112 clonally duplicated progenies from an outcross F2 of Eucalyptus globulus. * Two unlinked QTL were located for macrocarpal, while another unlinked QTL was located for sideroxylonal. The sideroxylonal QTL collocated with one for total sideroxylonal previously reported using adult Eucalyptus nitens foliage, providing independent validation in a different evolutionary lineage and a different ontogenetic stage. * Given the potential widespread occurrence of these QTL, their ontogenetic stability, and their impact on a range of dependent herbivores, it is possible that they have extended phenotypic effects in the Australian forest landscape.

Early ovule development following self- and cross-pollinations in Eucalyptus globulus Labill. ssp. globulus.

The study was conducted to identify the self-incompatibility mechanism in Eucalyptus globulus ssp. globulus. Controlled self- and cross-pollinations were conducted on individual flowers from three mature trees that had self-incompatibility levels of 76, 99.6 and 100%. Flowers were harvested at 4, 6 and 8 weeks after pollination. Embryology was investigated by bright field microscopy on material harvested at 4 and 6 weeks after pollination. Fertilization had taken place at 4 weeks after pollination with zygotes and free nuclear endosperm visible. There was a greater proportion of healthy, fertilized ovules in the cross- compared with the self-pollination treatment, and approx. half the ovules examined from both pollen treatments were not fertilized or were degenerating. By 6 weeks after pollination a few zygotes were starting to divide. The number of healthy, fertilized ovules was still greater in the cross-pollination treatment, but the number of healthy fertilized ovules was lower in both treatments compared with 4 weeks after pollination, and many ovules were degenerating. Fertilized ovules were significantly larger than non-fertilized or degenerating ovules and this difference was detectable by eye at 6 and 8 weeks after pollination. The mechanism of self-incompatibility appears to have both late pre- and post-zygotic components.

Maternal inheritance of the chloroplast genome in Eucalyptus globulus and interspecific hybrids.

The utility of chloroplast DNA (cpDNA) in Eucalyptus, either as a molecular marker for genetic studies or as a potential vehicle for genetic manipulation, is based on knowledge of its mode of inheritance. Chloroplast inheritance in angiosperms can vary among and within species, and anomalous inheritance has been reported in some interspecific-hybrid combinations. In Eucalyptus, abnormalities of pollen-tube growth occur in a number of interspecific-hybrid combinations, and this might increase the likelihood of anomalous chloroplast transmission. We used a rapid PCR technique to determine chloroplast heritability in 425 progeny of Eucalyptus, comprising 194 progeny of the premier pulpwood species E. globulus and 231 interspecific hybrids between E. globulus and E. nitens (F1, F2, and backcrosses). At this sampling intensity, no pollen-mediated transmission of cpDNA was found in any of the 40 families tested. The results are discussed with reference to chloroplast engineering and the use of cpDNA as a seed-specific marker in phylogeographic studies of Eucalyptus.

Chloroplast sharing in the Tasmanian eucalypts.

The biogeographic pattern of chloroplast DNA (cpDNA) haplotypes in Eucalyptus on the island of Tasmania is consistent with reticulate evolution, involving at least 12 Tasmanian species from the subgenus Symphyomyrtus. Intraspecific cpDNA polymorphism in 14 of 17 species is coupled with extensive sharing of identical haplotypes across populations of different species in the same geographic area. Haplotype diversity is lowest in central regions of Tasmania formerly occupied by alpine vegetation during glacial intervals and in northern regions that were periodically linked to continental Australia by land bridges. The observed distribution of several cpDNA haplotypes unique to Tasmania coincides with modeled locations of glacial refugia in coastal areas of Tasmania and shows the power of cpDNA in unraveling the complex history of past distributions of Eucalyptus. The results suggest that the model of evolution of the eucalypts should be reassessed to allow for the anastomosing effects of interspecific hybridization and introgression.

F1 hybrid inviability in eucalyptus: the case of E. ovata x E. globulus.

The impact of inbreeding and hybridization on fitness was compared in the two co-occurring forest tree species, Eucalyptus ovata and E. globulus, aimed at explaining the rarity of their hybrids in nature. The success of selfing, open-pollination and outcrossing of both species and interspecific hybridization was monitored from seed-set to 10-year's growth in a field trial. There was a unilateral barrier to hybridization with seed-set obtained only with E. ovata females. The F1 hybrids exhibited reduced viability compared to intraspecific cross-types at virtually all stages of the life cycle and are clearly at a selective disadvantage compared with their open-pollinated E. ovata half-sibs with which they would directly compete in nature. Eucalyptus ovata and E. globulus overlap in their flowering time but the F1 hybrids flowered later with virtually no overlap with either species. The asynchronous flowering and reduced reproductive fitness of F1 hybrids would markedly limit the opportunity for advanced generation hybridization. Inbreeding similarly had a deleterious effect on the fitness of both species, and the F1 hybrids were most competitive with the E. ovata selfs. It is argued that changes in inbreeding levels of parental populations may be a key factor affecting the relative fitness of hybrids and their potential to impact on the pure species gene pool. Reduced fitness of the pure species through inbreeding may result in hybridization having its greatest evolutionary impact in small founder or relict populations.

Plant genetics affects arthropod community richness and composition: evidence from a synthetic eucalypt hybrid population.

To examine how genetic variation in a plant population affects arthropod community richness and composition, we quantified the arthropod communities on a synthetic population of Eucalyptus amygdalina, E. risdonii, and their F1 and advanced-generation hybrids. Five major patterns emerged. First, the pure species and hybrid populations supported significantly different communities. Second, species richness was significantly greatest on hybrids (F1 > F2 > E. amygdalina > E. risdonii). These results are similar to those from a wild population of the same species and represent the first case in which both synthetic and wild population studies confirm a genetic component to community structure. Hybrids also acted as centers of biodiversity by accumulating both the common and specialist taxa of both parental species (100% in the wild and 80% in the synthetic population). Third, species richness was significantly greater on F1s than the single F2 family, suggesting that the increased insect abundance on hybrids may not be caused by the breakup of coadapted gene complexes. Fourth, specialist arthropod taxa were most likely to show a dominance response to F1 hybrids, whereas generalist taxa exhibited a susceptible response. Fifth, in an analysis of 31 leaf terpenoids that are thought to play a role in plant defense, hybrids were generally intermediate to the parental chemotypes. Within the single F2 family, we found significant associations between the communities of individual trees and five individual oil components, including oil yield, demonstrating that there is a genetic effect on plant defensive chemistry that, in turn, may affect community structure. These studies argue that hybridization has important community-level consequences and that the genetic variation present in hybrid zones can be used to explore the genetic-based mechanisms that structure communities.

Incongruence between chloroplast and species phylogenies in Eucalyptus subgenus Monocalyptus (Myrtaceae).

Seventy-eight polymorphic cpDNA (chloroplast DNA) characters were found in 13 closely related taxa from Eucalyptus series Amygdalinae (subgenus Monocalyptus) and seven potential outgroup taxa. The strict consensus of six cladograms generated from cpDNA data confirmed monophyly of Monocalyptus. However, cpDNA phylogeny within Monocalyptus was incongruent with taxonomic classification, being more related to geography, even when accessions were from divergent series. Monocalyptus cpDNA formed two major clades. On the island of Tasmania cpDNA was restricted to a single clade, exhibited very little variation, and was phylogenetically related to cpDNA found in central and western Victoria. In contrast, cpDNA of mainland monocalypt taxa was more variable, even within the Amygdalinae. Four out of six Tasmanian Amygdalinae species were polymorphic. The difference between cpDNA of replicates was often greater than differences between species from different series. The low level of cpDNA variation and extensive morphological intergradation between the Tasmanian endemics suggest recent speciation. However, the transfer of cpDNA through hybridization between lineages is the most likely explanation for the observed sharing of cpDNA across series. This study highlights that the geographical pattern to cpDNA variation in Eucalyptus may be an important source of information on past plant distributions in Australia.

ITS sequence data resolve higher level relationships among the eucalypts.

Sequences of the internal transcribed spacer (ITS) region of the nuclear ribosomal DNA were obtained for 35 species of Eucalyptus s.s. and seven taxa representing five outgroup genera (Allosyncarpia, Angophora, Arillastrum, Corymbia, and Stockwellia). The sequences were analyzed cladistically. The data distinguished clearly between the two major subgenera of Eucalyptus s.s. (Symphyomyrtus and Monocalyptus) but indicated that subgenus Eudesmia may be paraphyletic. ITS sequence data demonstrated the potential to resolve relationships between sections within subgenus Symphyomyrtus. Within sections, however, taxa were poorly differentiated. At the generic level, Corymbia appeared to be paraphyletic due to the exclusion of Angophora. The positions of Allosyncarpia and Arillastrum relative to the ingroup remain unresolved. ITS sequence data may prove valuable for resolving other phylogenetic relationships at higher taxonomic levels within Eucalyptus.

A comparison of genetic information from open-pollinated and control-pollinated progeny tests in two eucalypt species.

Genetic-parameter estimates and parental breeding-value predictions were compared from open-pollinated and control-pollinated progeny populations of Eucalyptus globulus and two populations of E. nitens. For E. globulus there were two types of open-pollinated populations (native stand open-pollinated and seed orchard open-pollinated) and two types of control-pollinated populations (intra-provenance and interprovenance full-sib families). For E. nitens there were two populations, a seed orchard open-pollinated population and intra-provenance full-sib families. Progeny tests were established across multiple sites and 2-year height and diameter were measured and volume calculated. Genetic parameters from native stand open-pollinated E. globulus were unlike the parameters from the other three E. globulus populations; heritability estimates were severely inflated, presumably due to high levels, and possibly differential levels, of inbreeding depression relative to the other populations. Estimates of dominance variance in the E. globulus full-sib populations were high, but were zero in the E. nitens population. Correlations among parental breeding values, predicted using data from the different populations, were generally low and non-significant, with two exceptions: predictions from the two E. globulus full-sib populations were significantly correlated (r=0.54, P = 0.001), as were predictions from the E. nitens seed orchard OP and full-sib population (r = 0.61, P = 0.08). There was some indication that superior parents of E. globulus native stand open-pollinated families also tended to have above-average breeding values based on the performance of intra-provenance full-sib offspring. The consequences of these results for exploitation of base-population collections from native stands are discussed.

Plant hybrid zones as centers of biodiversity: the herbivore community of two endemic Tasmanian eucalypts.

We found the hybrid zone between Eucalyptus amygdalina and Eucalyptus risdonii to be a center of insect and fungal species richness and abundance. Of 40 taxa examined, 73% were significantly more abundant in the hybrid zone than in pure zones, 25% showed on significant differences, and 2% were most abundant on a pure host species. The average hybrid tree supported 53% more insect and fungal species, and relative abundances were, on average, 4 times greater on hybrids than on either eucalypt species growing in pure stands. Hybrids may act as refugia for rare species: 5 of 40 species were largely restricted to the hybrid zone. Also, 50% of the species coexisted only in the hybrid zone, making for mique species assemblages. Although hybrids support more species and greater abundances, all hybrids are not equal: 68% of the 40 taxa examined were significantly more abundant on one hybrid phenotype than another. While herbivore concentrations on F1 type intermediates were rare, concentrations were common on phenotypes resembling backcrosses either to E. amygdalina or E. risdonii. For specialist herbivores, the hybrid phenotype most heavily utilized appears to be determined by its phenotypic affinity to its host species. Generalists exhibit an overall greater abundance on hybrids, but are less likely to utilize one hybrid phenotype over another. Mechanistic explanations for these distributions are numerous and probably species specific, but are likely to include: increased genetic susceptibility of hybrids due to hybrid breakdown; increased stress in the hybrid zone resulting in greater plant susceptibility; and a greater diversity of resources in the hybrid zone which could support more species. Seed capsule production by hybrids and their parental species is negatively correlated with herbivory. However, it is difficult to determine whether herbivores cause this pattern as hybrids may have inherently lower sexual reproduction. Laws enacted to protect rare and endangered species do not include hybrids. We argue that a re-examination of our current "hybrid policy" is warranted. Plant hybrid zones are centers of plant evolution and speciation, sources of economically important plants and potential biocontrol agents, and, as our study suggests, also provide essential habitats for phytophagous communities.

Inheritance of freezing resistance in interspecific F1 hybrids of Eucalyptus.

The inheritance of freezing resistance in interspecific F1 hybrid families of Eucalyptus encompassing 27 different species combinations and a range of levels of hardiness was examined. Freezing resistance was assessed by determining the temperatures required to cause either 30% (T30), 40% (T40), or 50% (T50) leakage of electrolytes from excised leaf discs subjected to artificial freezing. Highly significant variation in freezing resistance occurred between species; the maximum difference between parents in any specific combination was over 9°C (E. gunnii x E. globulus). Freezing resistance was inherited in a predominantly additive manner in interspecific hybrids, although there was a tendency towards partial dominance toward the more sensitive species in some combinations (e.g., E. nitens x E. Globulus, E. nitens x E. camaldulensis, E. gunnii x E. globulus). The full expression of this genetic variation appeared to increase with hardiness and in some cases appeared to vary with ontogeny. Estimates of individual narrow-sense heritability of freezing resistance for pure E. nitens families were h (2) = 0.66±0.44 and 0.46±0.44. Across all species combinations examined, the heritability of F1 family means estimated from midparent regression was h (2) = 0.76±0.06 and h (2) = 0.89±0.06 for T40 and T50 values, respectively. The advantage of using selected parents for interspecific hybridization is demonstrated and the implications of these results for breeding for freezing resistance in Eucalyptus are discussed.