High density linkage maps, genetic architecture, and genomic prediction of growth and wood properties in Pinus radiata.

BACKGROUND: The growing availability of genomic resources in radiata pine paves the way for significant advances in fundamental and applied genomic research. We constructed robust high-density linkage maps based on exome-capture genotyping in two F1 populations, and used these populations to perform quantitative trait locus (QTL) scans, genomic prediction and quantitative analyses of genetic architecture for key traits targeted by tree improvement programmes. RESULTS: Our mapping approach used probabilistic error correction of the marker data, followed by an iterative approach based on stringent parameters. This approach proved highly effective in producing high-density maps with robust marker orders and realistic map lengths (1285-4674 markers per map, with sizes ranging from c. 1643-2292 cM, and mean marker intervals of 0.7-2.1 cM). Colinearity was high between parental linkage maps, although there was evidence for a large chromosomal rearrangement (affecting ~ 90 cM) in one of the parental maps. In total, 28 QTL were detected for growth (stem diameter) and wood properties (wood density and fibre properties measured by Silviscan) in the QTL discovery population, with 1-3 QTL of small to moderate effect size detected per trait in each parental map. Four of these QTL were validated in a second, unrelated F1 population. Results from genomic prediction and analyses of genetic architecture were consistent with those from QTL scans, with wood properties generally having moderate to high genomic heritabilities and predictive abilities, as well as somewhat less complex genetic architectures, compared to growth traits. CONCLUSIONS: Despite the economic importance of radiata pine as a plantation forest tree, robust high-density linkage maps constructed from reproducible, sequence-anchored markers have not been published to date. The maps produced in this study will be a valuable resource for several applications, including the selection of marker panels for genomic prediction and anchoring a recently completed de novo whole genome assembly. We also provide the first map-based evidence for a large genomic rearrangement in radiata pine. Finally, results from our QTL scans, genomic prediction, and genetic architecture analyses are informative about the genomic basis of variation in important phenotypic traits.

Genetic analysis of early phenology in lentil identifies distinct loci controlling component traits.

Modern-day domesticated lentil germplasm is generally considered to form three broad adaptation groups: Mediterranean, South Asian, and northern temperate, which correspond to the major global production environments. Reproductive phenology plays a key role in lentil adaptation to this diverse ecogeographic variation. Here, we dissect the characteristic earliness of the pilosae ecotype, suited to the typically short cropping season of South Asian environments. We identified two loci, DTF6a and DTF6b, at which dominant alleles confer early flowering, and we show that DTF6a alone is sufficient to confer early flowering under extremely short photoperiods. Genomic synteny confirmed the presence of a conserved cluster of three florigen (FT) gene orthologues among potential candidate genes, and expression analysis in near-isogenic material showed that the early allele is associated with a strong derepression of the FTa1 gene in particular. Sequence analysis revealed a 7.4 kb deletion in the FTa1-FTa2 intergenic region in the pilosae parent, and a wide survey of >350 accessions with diverse origin showed that the dtf6a allele is predominant in South Asian material. Collectively, these results contribute to understanding the molecular basis of global adaptation in lentil, and further emphasize the importance of this conserved genomic region for adaptation in temperate legumes generally.

Pests, diseases, and aridity have shaped the genome of Corymbia citriodora.

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent.

Independent genetic control of drought resistance, recovery, and growth of Eucalyptus globulus seedlings.

Drought is a major stress impacting forest ecosystems worldwide. We utilized quantitative trait loci (QTL) analysis to study the genetic basis of variation in (a) drought resistance and recovery and (b) candidate traits that may be associated with this variation in the forest tree Eucalyptus globulus. QTL analysis was performed using a large outcrossed F2 mapping population from which 300 trees were phenotyped based on the mean performance of their open-pollinated F3 progeny. Progenies were grown in a glasshouse in a randomized complete block design. A subset of seedlings was subjected to a drought treatment after which they were rewatered and scored for damage and growth postdrought. Nondroughted seedlings were assessed for growth traits as well as lignotuber size and resprouting following severe damage to the main stem. QTL were detected for most traits. Importantly, independent QTL were detected for (a) drought damage and plant size, (b) drought damage and growth recovery, and (c) lignotuber size and resprouting capacity. Such independence argues that trade-offs are unlikely to be a major limitation to the response to selection and at the early life history stage studied; there are opportunities to improve resilience to drought without adverse effects on productivity.

Altered Expression of an FT Cluster Underlies a Major Locus Controlling Domestication-Related Changes to Chickpea Phenology and Growth Habit.

Flowering time is a key trait in breeding and crop evolution, due to its importance for adaptation to different environments and for yield. In the particular case of chickpea, selection for early phenology was essential for the successful transition of this species from a winter to a summer crop. Here, we used genetic and expression analyses in two different inbred populations to examine the genetic control of domestication-related differences in flowering time and growth habit between domesticated chickpea and its wild progenitor Cicer reticulatum. A single major quantitative trait locus for flowering time under short-day conditions [Days To Flower (DTF)3A] was mapped to a 59-gene interval on chromosome three containing a cluster of three FT genes, which collectively showed upregulated expression in domesticated relative to wild parent lines. An equally strong association with growth habit suggests a pleiotropic effect of the region on both traits. These results indicate the likely molecular explanation for the characteristic early flowering of domesticated chickpea, and the previously described growth habit locus Hg. More generally, they point to de-repression of this specific gene cluster as a conserved mechanism for achieving adaptive early phenology in temperate legumes.

Comparison of host susceptibilities to native and exotic pathogens provides evidence for pathogen-imposed selection in forest trees.

The extent to which spatial structuring of host resistance in wild plant populations reflects direct pathogen-imposed selection is a subject of debate. To examine this issue, genetic susceptibilities to an exotic and a coevolved native fungal pathogen were compared using two Australian host tree species. Damage to common host germplasm of Corymbia citriodora ssp. variegata (CCV) and Eucalyptus globulus, caused by recently introduced (Austropuccinia psidii) and native (Quambalaria pitereka and Teratosphaeria sp.) pathogens was evaluated in common-garden experiments. There was significant additive genetic variation within host species for susceptibility to both the exotic and native pathogens. However, susceptibility to A. psidii was not genetically correlated with susceptibility to either native pathogen, providing support for pathogen-specific rather than general mechanisms of resistance. Population differentiation (QST ) for susceptibility to the native pathogens was greater than neutral expectations (molecular FST ), arguing for divergent selection. Coupled with lower native, but not exotic, pathogen susceptibility in host populations from areas climatically more prone to fungal proliferation, these findings suggest that pathogen-imposed selection has contributed directly to a geographic mosaic of host resistance to native pathogens.

Annotation of the Corymbia terpene synthase gene family shows broad conservation but dynamic evolution of physical clusters relative to Eucalyptus.

Terpenes are economically and ecologically important phytochemicals. Their synthesis is controlled by the terpene synthase (TPS) gene family, which is highly diversified throughout the plant kingdom. The plant family Myrtaceae are characterised by especially high terpene concentrations, and considerable variation in terpene profiles. Many Myrtaceae are grown commercially for terpene products including the eucalypts Corymbia and Eucalyptus. Eucalyptus grandis has the largest TPS gene family of plants currently sequenced, which is largely conserved in the closely related E. globulus. However, the TPS gene family has been well studied only in these two eucalypt species. The recent assembly of two Corymbia citriodora subsp. variegata genomes presents an opportunity to examine the conservation of this important gene family across more divergent eucalypt lineages. Manual annotation of the TPS gene family in C. citriodora subsp. variegata revealed a similar overall number, and relative subfamily representation, to that previously reported in E. grandis and E. globulus. Many of the TPS genes were in physical clusters that varied considerably between Eucalyptus and Corymbia, with several instances of translocation, expansion/contraction and loss. Notably, there was greater conservation in the subfamilies involved in primary metabolism than those involved in secondary metabolism, likely reflecting different selective constraints. The variation in cluster size within subfamilies and the broad conservation between the eucalypts in the face of this variation are discussed, highlighting the potential contribution of selection, concerted evolution and stochastic processes. These findings provide the foundation to better understand terpene evolution within the ecologically and economically important Myrtaceae.

Quantitative Trait Loci (QTLs) for Intumescence Severity in Eucalyptus globulus and Validation of QTL Detection Based on Phenotyping Using Open-Pollinated Families of a Mapping Population.

Intumescence is a nonpathogenic physiological disorder characterized by leaf blistering. This disorder can affect growth and development in glasshouses and growth chambers and may be confused with pathogenic diseases. We used quantitative trait loci (QTL) analysis to examine the genetic basis of variation in intumescence severity in Eucalyptus globulus, and test for colocation with previously detected QTLs for pathogen susceptibility. QTL analysis used the phenotype means of open-pollinated (OP) families of an outcrossed F2 mapping family (OP F3; n = 300) of E. globulus and the linkage map constructed in the F2. We validate this phenotyping approach for QTL analysis by assessing a trait previously used for QTL discovery in the F2 and showing the same major QTL was detected with the OP F3. For intumescence severity, five putative QTLs were detected across four linkage groups. Four of these did not colocate with previously reported QTLs for fungal pathogen susceptibility in Eucalyptus, suggesting the mechanisms underlying susceptibility to intumescence and to the two fungal pathogens are largely independent. This study demonstrates there is a genetic basis for variation in intumescence severity, reports the first QTL for intumescence severity in plants, and provides a robust framework for investigating the potential mechanisms involved.

Genome-wide variation in recombination rate in Eucalyptus.

BACKGROUND: Meiotic recombination is a fundamental evolutionary process. It not only generates diversity, but influences the efficacy of natural selection and genome evolution. There can be significant heterogeneity in recombination rates within and between species, however this variation is not well understood outside of a few model taxa, particularly in forest trees. Eucalypts are forest trees of global economic importance, and dominate many Australian ecosystems. We studied recombination rate in Eucalyptus globulus using genetic linkage maps constructed in 10 unrelated individuals, and markers anchored to the Eucalyptus reference genome. This experimental design provided the replication to study whether recombination rate varied between individuals and chromosomes, and allowed us to study the genomic attributes and population genetic parameters correlated with this variation. RESULTS: Recombination rate varied significantly between individuals (range = 2.71 to 3.51 centimorgans/megabase [cM/Mb]), but was not significantly influenced by sex or cross type (F1 vs. F2). Significant differences in recombination rate between chromosomes were also evident (range = 1.98 to 3.81 cM/Mb), beyond those which were due to variation in chromosome size. Variation in chromosomal recombination rate was significantly correlated with gene density (r = 0.94), GC content (r = 0.90), and the number of tandem duplicated genes (r = -0.72) per chromosome. Notably, chromosome level recombination rate was also negatively correlated with the average genetic diversity across six species from an independent set of samples (r = -0.75). CONCLUSIONS: The correlations with genomic attributes are consistent with findings in other taxa, however, the direction of the correlation between diversity and recombination rate is opposite to that commonly observed. We argue this is likely to reflect the interaction of selection and specific genome architecture of Eucalyptus. Interestingly, the differences amongst chromosomes in recombination rates appear stable across Eucalyptus species. Together with the strong correlations between recombination rate and features of the Eucalyptus reference genome, we maintain these findings provide further evidence for a broad conservation of genome architecture across the globally significant lineages of Eucalyptus.

Genetic control of cuticular wax compounds in Eucalyptus globulus.

Plant cuticular wax compounds perform functions that are essential for the survival of terrestrial plants. Despite their importance, the genetic control of these compounds is poorly understood outside of model taxa. Here we investigate the genetic basis of variation in cuticular compounds in Eucalyptus globulus using quantitative genetic and quantitative trait loci (QTL) analyses. Quantitative genetic analysis was conducted using 246 open-pollinated progeny from 13 native sub-races throughout the geographic range. QTL analysis was conducted using 112 clonally replicated progeny from an outcross F2 population. Nine compounds exhibited significant genetic variation among sub-races with three exhibiting signals of diversifying selection. Fifty-two QTL were found with co-location of QTL for related compounds commonly observed. Notable among these was the QTL for five wax esters, which co-located with a gene from the KCS family, previously implicated in the biosynthesis of cuticular waxes in Arabidopsis. In combination, the QTL and quantitative genetic analyses suggest the variation and differentiation in cuticular wax compounds within E. globulus has a complex genetic origin. Sub-races exhibited independent latitudinal and longitudinal differentiation in cuticular wax compounds, likely reflecting processes such as historic gene flow and diversifying selection acting upon genes that have diverse functions in distinct biochemical pathways.

Genomic patterns of species diversity and divergence in Eucalyptus.

We examined genome-wide patterns of DNA sequence diversity and divergence among six species of the important tree genus Eucalyptus and investigated their relationship with genomic architecture. Using c. 90 range-wide individuals of each Eucalyptus species (E. grandis, E. urophylla, E. globulus, E. nitens, E. dunnii and E. camaldulensis), genetic diversity and divergence were estimated from 2840 polymorphic diversity arrays technology markers covering the 11 chromosomes. Species differentiating markers (SDMs) identified in each of 15 pairwise species comparisons, along with species diversity (HHW ) and divergence (FST ), were projected onto the E. grandis reference genome. Across all species comparisons, SDMs totalled 1.1-5.3% of markers and were widely distributed throughout the genome. Marker divergence (FST and SDMs) and diversity differed among and within chromosomes. Patterns of diversity and divergence were broadly conserved across species and significantly associated with genomic features, including the proximity of markers to genes, the relative number of clusters of tandem duplications, and gene density within or among chromosomes. These results suggest that genomic architecture influences patterns of species diversity and divergence in the genus. This influence is evident across the six species, encompassing diverse phylogenetic lineages, geography and ecology.

Genetic control of heterochrony in Eucalyptus globulus.

A change in the timing or rate of developmental events throughout ontogeny is referred to as heterochrony, and it is a major evolutionary process in plants and animals. We investigated the genetic basis for natural variation in the timing of vegetative phase change in the tree Eucalyptus globulus, which undergoes a dramatic change in vegetative morphology during the juvenile-to-adult transition. Quantitative trait loci analysis in an outcross F2 family derived from crosses between individuals from a coastal population of E. globulus with precocious vegetative phase change and individuals from populations in which vegetative phase change occurs several years later implicated the microRNA EglMIR156.5 as a potential contributor to this heterochronic difference. Additional evidence for the involvement of EglMIR156.5 was provided by its differential expression in trees with early and late phase change. Our findings suggest that changes in the expression of miR156 underlie natural variation in vegetative phase change in E. globulus, and may also explain interspecific differences in the timing of this developmental transition.

Quantitative trait loci in hop (Humulus lupulus L.) reveal complex genetic architecture underlying variation in sex, yield and cone chemistry.

BACKGROUND: Hop (Humulus lupulus L.) is cultivated for its cones, the secondary metabolites of which contribute bitterness, flavour and aroma to beer. Molecular breeding methods, such as marker assisted selection (MAS), have great potential for improving the efficiency of hop breeding. The success of MAS is reliant on the identification of reliable marker-trait associations. This study used quantitative trait loci (QTL) analysis to identify marker-trait associations for hop, focusing on traits related to expediting plant sex identification, increasing yield capacity and improving bittering, flavour and aroma chemistry. RESULTS: QTL analysis was performed on two new linkage maps incorporating transferable Diversity Arrays Technology (DArT) markers. Sixty-three QTL were identified, influencing 36 of the 50 traits examined. A putative sex-linked marker was validated in a different pedigree, confirming the potential of this marker as a screening tool in hop breeding programs. An ontogenetically stable QTL was identified for the yield trait dry cone weight; and a QTL was identified for essential oil content, which verified the genetic basis for variation in secondary metabolite accumulation in hop cones. A total of 60 QTL were identified for 33 secondary metabolite traits. Of these, 51 were pleiotropic/linked, affecting a substantial number of secondary metabolites; nine were specific to individual secondary metabolites. CONCLUSIONS: Pleiotropy and linkage, found for the first time to influence multiple hop secondary metabolites, have important implications for molecular selection methods. The selection of particular secondary metabolite profiles using pleiotropic/linked QTL will be challenging because of the difficulty of selecting for specific traits without adversely changing others. QTL specific to individual secondary metabolites, however, offer unequalled value to selection programs. In addition to their potential for selection, the QTL identified in this study advance our understanding of the genetic control of traits of current economic and breeding significance in hop and demonstrate the complex genetic architecture underlying variation in these traits. The linkage information obtained in this study, based on transferable markers, can be used to facilitate the validation of QTL, crucial to the success of MAS.

Stability of quantitative trait loci for growth and wood properties across multiple pedigrees and environments in Eucalyptus globulus.

· Eucalypts are one of the most planted tree genera worldwide, and there is increasing interest in marker-assisted selection for tree improvement. Implementation of marker-assisted selection requires a knowledge of the stability of quantitative trait loci (QTLs). This study aims to investigate the stability of QTLs for wood properties and growth across contrasting sites and multiple pedigrees of Eucalyptus globulus. · Saturated linkage maps were constructed using 663 genotypes from four separate families, grown at three widely separated sites, and were employed to construct a consensus map. This map was used for QTL analysis of growth, wood density and wood chemical traits, including pulp yield. · Ninety-eight QTLs were identified across families and sites: 87 for wood properties and 11 for growth. These QTLs mapped to 38 discrete regions, some of which co-located with candidate genes. Although 16% of QTLs were verified across different families, 24% of wood property QTLs and 38% of growth QTLs exhibited significant genotype-by-environment interaction. · This study provides the most detailed assessment of the effect of environment and pedigree on QTL detection in the genus. Despite markedly different environments and pedigrees, many QTLs were stable, providing promising targets for the application of marker-assisted selection.

A reference linkage map for Eucalyptus.

BACKGROUND: Genetic linkage maps are invaluable resources in plant research. They provide a key tool for many genetic applications including: mapping quantitative trait loci (QTL); comparative mapping; identifying unlinked (i.e. independent) DNA markers for fingerprinting, population genetics and phylogenetics; assisting genome sequence assembly; relating physical and recombination distances along the genome and map-based cloning of genes. Eucalypts are the dominant tree species in most Australian ecosystems and of economic importance globally as plantation trees. The genome sequence of E. grandis has recently been released providing unprecedented opportunities for genetic and genomic research in the genus. A robust reference linkage map containing sequence-based molecular markers is needed to capitalise on this resource. Several high density linkage maps have recently been constructed for the main commercial forestry species in the genus (E. grandis, E. urophylla and E. globulus) using sequenced Diversity Arrays Technology (DArT) and microsatellite markers. To provide a single reference linkage map for eucalypts a composite map was produced through the integration of data from seven independent mapping experiments (1950 individuals) using a marker-merging method. RESULTS: The composite map totalled 1107 cM and contained 4101 markers; comprising 3880 DArT, 213 microsatellite and eight candidate genes. Eighty-one DArT markers were mapped to two or more linkage groups, resulting in the 4101 markers being mapped to 4191 map positions. Approximately 13% of DArT markers mapped to identical map positions, thus the composite map contained 3634 unique loci at an average interval of 0.31 cM. CONCLUSION: The composite map represents the most saturated linkage map yet produced in Eucalyptus. As the majority of DArT markers contained on the map have been sequenced, the map provides a direct link to the E. grandis genome sequence and will serve as an important reference for progressing eucalypt research.

Few Mendelian genes underlie the quantitative response of a forest tree, Eucalyptus globulus, to a natural fungal epidemic.

Foliar fungal pathogens from the genus Mycosphaerella affect eucalyptus in natural forests and plantations worldwide. QTL analysis was conducted to dissect the genetic control of resistance in Eucalyptus globulus to a natural infection by Mycosphaerella leaf disease, using a clonally replicated outbred F2 family (112 genotypes) planted in a field trial. Two major QTL, with high LOD support (20.2 and 10.9) and high genomewide significance, explained a large proportion (52%) of the phenotypic variance in the severity of damage by Mycosphaerella cryptica, which may be indicative of oligogenic control. Both QTL were validated in a second F2 family and one was validated in a third F2 family. The mean values of different genotype classes at both major QTL argue for Mendelian inheritance with resistance dominant over susceptibility. There were strong correlations between the levels of Mycosphaerella damage in related genetic material planted in three widely separated locations in Tasmania. These findings together provide evidence that the genes controlling resistance to Mycosphaerella damage are stable in different genetic backgrounds and across different environments.