Search

Genomic Selection (GS) in tree breeding optimizes genetic gains by leveraging genomic data to enable early selection of seedlings without phenotypic data reducing breeding cycle and increasing selection intensity. Traditional assessments of the potential of GS in forest trees have typically focused on model performance using cross-validation within the same generation but evaluating effectively realized predictive ability (RPA) across generations is crucial. This study estimated RPAs for volume growth (VOL), wood density (WD), and pulp yield (PY) across four generations breeding of Eucalyptus grandis . The training set spanned three generations, including 34,461 trees with three-year growth data, 6,014 trees with wood quality trait data, and 1,918 trees with 12,695 SNPs (single nucleotide polymorphisms) data. Employing single-step genomic BLUP, we compared the genomic predictions of breeding values (GEBVs) for 1,153 fourth-generation full-sib seedlings in the greenhouse with their later-collected phenotypic estimated breeding values (EBVs) at age three years. RPAs were estimated using three GS targets (individual trees, trees within families, and families), two selection criteria (single- and multiple-trait), and training populations of either all 1,918 genotyped trees or the 67 direct ancestors of the selection candidates. RPAs were higher for wood quality traits (0.33 to 0.59) compared to VOL (0.14 to 0.19) and improved for wood traits (0.42 to 0.75) but not for VOL when trained only with direct ancestors, highlighting the challenges in accurately predicting growth traits. GS was more effective at excluding bottom-ranked candidates than selecting top-ranked ones. The between-family GS approach outperformed individual-tree selection for VOL (0.11 to 0.16) and PY (0.72 to 0.75), but not for WD (0.43 vs. 0.42). Furthermore, higher levels of relatedness and lower genotype by environment (G × E) interaction between training and testing populations enhanced RPAs for VOL (0.39). In summary, despite limited effectiveness in ranking top VOL individuals, GS effectively identified low-performing individuals and families. These multi-generational findings underscore GS's potential in tree breeding, stressing the importance of considering relatedness and G × E interaction for optimal performance., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Duarte, Jurcic, Dutour, Villalba, Centurión, Grattapaglia and Cappa.)

Introduction: The coconut tree crop ( Cocos nucifera L.) provides vital resources for millions of people worldwide. Coconut germplasm is largely classified into 'Tall' (Typica) and 'Dwarf' (Nana) types. While Tall coconuts are outcrossing, stress tolerant, and late flowering, Dwarf coconuts are inbred and flower early with a high rate of bunch emission. Precocity determines the earlier production of a plantation and facilitates management and harvest., Methods: A unique outbred F 2 population was used, generated by intercrossing F 1 hybrids between Brazilian Green Dwarf from Jiqui (BGDJ) and West African Tall (WAT) cultivars. Single-nucleotide polymorphism (SNP) markers fixed for alternative alleles in the two varieties, segregating in an F 2 configuration, were used to build a high-density linkage map with ~3,000 SNPs, anchored to the existing chromosome-level genome assemblies, and a quantitative trait locus (QTL) mapping analysis was carried out., Results: The linkage map established the chromosome numbering correspondence between the two reference genome versions and the relationship between recombination rate, physical distance, and gene density in the coconut genomes. Leveraging the strong segregation for precocity inherited from the Dwarf cultivar in the F 2 , a major effect QTL with incomplete dominance was mapped for flowering time. FLOWERING LOCUS T ( FT ) gene homologs of coconut previously described as putatively involved in flowering time by alternative splice variant analysis were colocalized within a ~200-kb window of the major effect QTL [logarithm of the odds (LOD) = 11.86]., Discussion: Our work provides strong phenotype-based evidence for the role of the FT locus as the putative underlying functional variant for the flowering time difference between Dwarf and Tall coconuts. Major effect QTLs were also detected for developmental traits of the palm, plausibly suggesting pleiotropism of the FT locus for other precocity traits. Haplotypes of the two SNPs flanking the flowering time QTL inherited from the Dwarf parent BGDJ caused a reduction in the time to flower of approximately 400 days. These SNPs could be used for high-throughput marker-assisted selection of early-flowering and higher-productivity recombinant lines, providing innovative genetic material to the coconut industry., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Grattapaglia, Alves and Pacheco.)

Eucalyptus dunnii is one of the most important Eucalyptus species for short-fiber pulp production in regions where other species of the genus are affected by poor soil and climatic conditions. In this context, E. dunnii holds promise as a resource to address and adapt to the challenges of climate change. Despite its rapid growth and favorable wood properties for solid wood products, the advancement of its improvement remains in its early stages. In this work, we evaluated the performance of two single nucleotide polymorphism, (SNP), genotyping methods for population genetics analysis and Genomic Selection in E. dunnii . Double digest restriction-site associated DNA sequencing (ddRADseq) was compared with the EUChip60K array in 308 individuals from a provenance-progeny trial. The compared SNP set included 8,011 and 19,008 informative SNPs distributed along the 11 chromosomes, respectively. Although the two datasets differed in the percentage of missing data, genome coverage, minor allele frequency and estimated genetic diversity parameters, they revealed a similar genetic structure, showing two subpopulations with little differentiation between them, and low linkage disequilibrium. GS analyses were performed for eleven traits using Genomic Best Linear Unbiased Prediction (GBLUP) and a conventional pedigree-based model (ABLUP). Regardless of the SNP dataset, the predictive ability (PA) of GBLUP was better than that of ABLUP for six traits (Cellulose content, Total and Ethanolic extractives, Total and Klason lignin content and Syringyl and Guaiacyl lignin monomer ratio). When contrasting the SNP datasets used to estimate PAs, the GBLUP-EUChip60K model gave higher and significant PA values for six traits, meanwhile, the values estimated using ddRADseq gave higher values for three other traits. The PAs correlated positively with narrow sense heritabilities, with the highest correlations shown by the ABLUP and GBLUP-EUChip60K. The two genotyping methods, ddRADseq and EUChip60K, are generally comparable for population genetics and genomic prediction, demonstrating the utility of the former when subjected to rigorous SNP filtering. The results of this study provide a basis for future whole-genome studies using ddRADseq in non-model forest species for which SNP arrays have not yet been developed., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Aguirre, Villalba, García, Filippi, Rivas, Martínez, Acuña, López, López, Pathauer, Palazzini, Harrand, Oberschelp, Marcó, Cisneros, Carreras, Martins Alves, Rodrigues, Hopp, Grattapaglia, Cappa, Paniego and Marcucci Poltri.)

Introduction: Genomic selection (GS) experiments in forest trees have largely reported estimates of predictive abilities from cross-validation among individuals in the same breeding generation. In such conditions, no effects of recombination, selection, drift, and environmental changes are accounted for. Here, we assessed the effectively realized predictive ability (RPA) for volume growth at harvest age by GS across generations in an operational reciprocal recurrent selection (RRS) program of hybrid Eucalyptus ., Methods: Genomic best linear unbiased prediction with additive (GBLUP_G), additive plus dominance (GBLUP_G+D), and additive single-step (HBLUP) models were trained with different combinations of growth data of hybrids and pure species individuals ( N = 17,462) of the G 1 generation, 1,944 of which were genotyped with ~16,000 SNPs from SNP arrays. The hybrid G 2 progeny trial (HPT267) was the GS target, with 1,400 selection candidates, 197 of which were genotyped still at the seedling stage, and genomically predicted for their breeding and genotypic values at the operational harvest age (6 years). Seedlings were then grown to harvest and measured, and their pedigree-based breeding and genotypic values were compared to their originally predicted genomic counterparts., Results: Genomic RPAs ≥0.80 were obtained as the genetic relatedness between G 1 and G 2 increased, especially when the direct parents of selection candidates were used in training. GBLUP_G+D reached RPAs ≥0.70 only when hybrid or pure species data of G 1 were included in training. HBLUP was only marginally better than GBLUP. Correlations ≥0.80 were obtained between pedigree and genomic individual ranks. Rank coincidence of the top 2.5% selections was the highest for GBLUP_G (45% to 60%) compared to GBLUP_G+D. To advance the pure species RRS populations, GS models were best when trained on pure species than hybrid data, and HBLUP yielded ~20% higher predictive abilities than GBLUP, but was not better than ABLUP for ungenotyped trees., Discussion: We demonstrate that genomic data effectively enable accurate ranking of eucalypt hybrid seedlings for their yet-to-be observed volume growth at harvest age. Our results support a two-stage GS approach involving family selection by average genomic breeding value, followed by within-top-families individual GS, significantly increasing selection intensity, optimizing genotyping costs, and accelerating RRS breeding., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Simiqueli, Resende, Takahashi, de Sousa and Grattapaglia.)

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Eucalyptus is an economically important genus comprising more than 890 species in different subgenera and sections. Approximately twenty species of subgenus Symphyomyrtus account for 95% of the world's planted eucalypts. Discrimination of closely related eucalypt taxa is challenging, consistent with their recent phylogenetic divergence and occasional hybridization in nature. Admixture, misclassification or mislabeling of Eucalyptus germplasm resources maintained as exotics have been suggested, although no reports are available. Moreover, hybrids with increased productivity and traits complementarity are planted worldwide, but little is known about their actual genomic ancestry. In this study we examined a set of 440 trees of 16 different Eucalyptus species and 44 interspecific hybrids of multi-species origin conserved in germplasm banks in Brazil. We used genome-wide SNP data to evaluate the agreement between the alleged phylogenetic classification of species and provenances as registered in their historical records, and their observed genetic clustering derived from SNP data. Genetic structure analyses correctly assigned each of the 16 species to a different cluster although the PCA positioning of E. longirostrata was inconsistent with its current taxonomy. Admixture was present for closely related species' materials derived from local germplasm banks, indicating unintended hybridization following germplasm introduction. Provenances could be discriminated for some species, indicating that SNP-based discrimination was directly proportional to geographical distance, consistent with an isolation-by-distance model. SNP-based genomic ancestry analysis showed that the majority of the hybrids displayed realized genomic composition deviating from the expected ones based on their pedigree records, consistent with admixture in their parents and pervasive genome-wide directional selection toward the fast-growing E. grandis genome. SNP data in support of tree breeding provide precise germplasm identity verification, and allow breeders to objectively recognize the actual ancestral origin of superior hybrids to more realistically guide the program toward the development of the desired genetic combinations., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 de Oliveira et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Tropical fruit tree species constitute a yet untapped supply of outstanding diversity of taste and nutritional value, barely developed from the genetics standpoint, with scarce or no genomic resources to tackle the challenges arising in modern breeding practice. We generated a de novo genome assembly of the Psidium guajava, the super fruit "apple of the tropics", and successfully transferred 14,268 SNP probesets from Eucalyptus to Psidium at the nucleotide level, to detect genomic loci linked to resistance to the root knot nematode (RKN) Meloidogyne enterolobii derived from the wild relative P. guineense. Significantly associated loci with resistance across alternative analytical frameworks, were detected at two SNPs on chromosome 3 in a pseudo-assembly of Psidium guajava genome built using a syntenic path approach with the Eucalyptus grandis genome to determine the order and orientation of the contigs. The P. guineense-derived resistance response to RKN and disease onset is conceivably triggered by mineral nutrients and phytohormone homeostasis or signaling with the involvement of the miRNA pathway. Hotspots of mapped resistance quantitative trait loci and functional annotation in the same genomic region of Eucalyptus provide further indirect support to our results, highlighting the evolutionary conservation of genomes across genera of Myrtaceae in the adaptation to pathogens. Marker assisted introgression of the resistance loci mapped should accelerate the development of improved guava cultivars and hybrid rootstocks., Competing Interests: The authors have declared that no competing interests exist.

Cacao is a globally important crop with a long history of domestication and selective breeding. Despite the increased use of elite clones by cacao farmers, worldwide plantations are established mainly using hybrid progeny material derived from heterozygous parents, therefore displaying high tree-to-tree variability. The deliberate development of hybrids from advanced inbred lines produced by successive generations of self-pollination has not yet been fully considered in cacao breeding. This is largely due to the self-incompatibility of the species, the long generation cycles (3-5 years) and the extensive trial areas needed to accomplish the endeavor. We propose a simple and accessible approach to develop inbred lines based on accelerating the buildup of homozygosity based on regular selfing assisted by genome-wide SNP genotyping. In this study we genotyped 90 clones from the Brazilian CEPEC´s germplasm collection and 49 inbred offspring of six S1 or S2 cacao families derived from self-pollinating clones CCN-51, PS-13.19, TSH-1188 and SIAL-169. A set of 3,380 SNPs distributed across the cacao genome were interrogated on the EMBRAPA multi-species 65k Infinium chip. The 90 cacao clones showed considerable variation in genome-wide SNP homozygosity (mean 0.727± 0.182) and 19 of them with homozygosity ≥90%. By assessing the increase in homozygosity across two generations of self-pollinations, SNP data revealed the wide variability in homozygosity within and between S1 and S2 families. Even in small families (<10 sibs), individuals were identified with up to ~1.5 standard deviations above the family mean homozygosity. From baseline homozygosities of 0.476 and 0.454, offspring with homozygosities of 0.862 and 0.879 were recovered for clones TSH-1188 and CCN-51 respectively, in only two generations of selfing (81-93% increase). SNP marker assisted monitoring and selection of inbred individuals can be a practical tool to optimize and accelerate the development of inbred lines of outbred tree species. This approach will allow a faster and more accurate exploitation of hybrid breeding strategies in cacao improvement programs and potentially in other perennial fruit and forest trees., Competing Interests: The authors have declared that no competing interests exist.

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.

Key Message: We propose the application of enviromics to breeding practice, by which the similarity among sites assessed on an "omics" scale of environmental attributes drives the prediction of unobserved genotype performances. Genotype by environment interaction (GEI) studies in plant breeding have focused mainly on estimating genetic parameters over a limited number of experimental trials. However, recent geographic information system (GIS) techniques have opened new frontiers for better understanding and dealing with GEI. These advances allow increasing selection accuracy across all sites of interest, including those where experimental trials have not yet been deployed. Here, we introduce the term enviromics, within an envirotypic-assisted breeding framework. In summary, likewise genotypes at DNA markers, any particular site is characterized by a set of "envirotypes" at multiple "enviromic" markers corresponding to environmental variables that may interact with the genetic background, thus providing informative breeding re-rankings for optimized decisions over different environments. Based on simulated data, we illustrate an index-based enviromics method (the "GIS-GEI") which, due to its higher granular resolution than standard methods, allows for: (1) accurate matching of sites to their most appropriate genotypes; (2) better definition of breeding areas that have high genetic correlation to ensure selection gains across environments; and (3) efficient determination of the best sites to carry out experiments for further analyses. Environmental scenarios can also be optimized for productivity improvement and genetic resources management, especially in the current outlook of dynamic climate change. Envirotyping provides a new class of markers for genetic studies, which are fairly inexpensive, increasingly available and transferable across species. We envision a promising future for the integration of enviromics approaches into plant breeding when coupled with next-generation genotyping/phenotyping and powerful statistical modeling of genetic diversity.

Main Conclusion: Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

High-throughput SNP genotyping has become a precondition to move to higher precision and wider genome coverage genetic analysis of natural and breeding populations of non-model species. We developed a 44,318 annotated SNP catalog for Araucaria angustifolia, a grandiose subtropical conifer tree, one of the only two native Brazilian gymnosperms, critically endangered due to its valuable wood and seeds. Following transcriptome assembly and annotation, SNPs were discovered from RNA-seq and pooled RAD-seq data. From the SNP catalog, an Axiom® SNP array with 3,038 validated SNPs was developed and used to provide a comprehensive look at the genetic diversity and structure of 15 populations across the natural range of the species. RNA-seq was a far superior source of SNPs when compared to RAD-seq in terms of conversion rate to polymorphic markers on the array, likely due to the more efficient complexity reduction of the huge conifer genome. By matching microsatellite and SNP data on the same set of A. angustifolia individuals, we show that SNPs reflect more precisely the actual genome-wide patterns of genetic diversity and structure, challenging previous microsatellite-based assessments. Moreover, SNPs corroborated the known major north-south genetic cline, but allowed a more accurate attribution to regional versus among-population differentiation, indicating the potential to select ancestry-informative markers. The availability of a public, user-friendly 3K SNP array for A. angustifolia and a catalog of 44,318 SNPs predicted to provide ~29,000 informative SNPs across ~20,000 loci across the genome, will allow tackling still unsettled questions on its evolutionary history, toward a more comprehensive picture of the origin, past dynamics and future trend of the species' genetic resources. Additionally, but not less importantly, the SNP array described, unlocks the potential to adopt genomic prediction methods to accelerate the still very timid efforts of systematic tree breeding of A. angustifolia., Competing Interests: The authors have declared that no competing interests exist.

Several studies suggest the relation of DNA methylation to diseases in humans and important phenotypes in plants drawing attention to this epigenetic mark as an important source of variability. In the last decades, several methodologies were developed to assess the methylation state of a genome. However, there is still a lack of affordable and precise methods for genome wide analysis in large sample size studies. Methyl sensitive double digestion MS-DArT sequencing method emerges as a promising alternative for methylation profiling. We developed a computational pipeline for the identification of DNA methylation using MS-DArT-seq data and carried out a pilot study using the Eucalyptus grandis tree sequenced for the species reference genome. Using a statistic framework as in differential expression analysis, 72,515 genomic sites were investigated and 5,846 methylated sites identified, several tissue specific, distributed along the species 11 chromosomes. We highlight a bias towards identification of DNA methylation in genic regions and the identification of 2,783 genes and 842 transposons containing methylated sites. Comparison with WGBS, DNA sequencing after treatment with bisulfite, data demonstrated a precision rate higher than 95% for our approach. The availability of a reference genome is useful for determining the genomic context of methylated sites but not imperative, making this approach suitable for any species. Our approach provides a cost effective, broad and reliable examination of DNA methylation profile on MspI/HpaII restriction sites, is fully reproducible and the source code is available on GitHub (https://github.com/wendelljpereira/ms-dart-seq)., Competing Interests: Authors MRP and DG are employed by EMBRAPA, but this commercial affiliation does not alter our adherence to all PLOS ONE policies on sharing data and materials.

The role of natural selection in shaping patterns of diversity is still poorly understood in the Neotropics. We carried out the first genome-wide population genomics study in a Neotropical tree, Handroanthus impetiginosus (Bignoniaceae), sampling 75,838 SNPs by sequence capture in 128 individuals across 13 populations. We found evidences for local adaptation using Bayesian correlations of allele frequency and environmental variables (32 loci in 27 genes) complemented by an analysis of selective sweeps and genetic hitchhiking events using SweepFinder2 (81 loci in 47 genes). Fifteen genes were identified by both approaches. By accounting for population genetic structure, we also found 14 loci with selection signal in a STRUCTURE-defined lineage comprising individuals from five populations, using Outflank. All approaches pinpointed highly diverse and structurally conserved genes affecting plant development and primary metabolic processes. Spatial interpolation forecasted differences in the expected allele frequencies at loci under selection over time, suggesting that H. impetiginosus may track its habitat during climate changes. However, local adaptation through natural selection may also take place, allowing species persistence due to niche evolution. A high genetic differentiation was seen among the H. impetiginosus populations, which, together with the limited power of the experiment, constrains the improved detection of other types of soft selective forces, such as background, balanced, and purifying selection. Small differences in allele frequency distribution among widespread populations and the low number of loci with detectable adaptive sweeps advocate for a polygenic model of adaptation involving a potentially large number of small genome-wide effects.

Genomic Best Linear Unbiased Prediction (GBLUP) in tree breeding typically only uses information from genotyped trees. However, information from phenotyped but non-genotyped trees can also be highly valuable. The single-step GBLUP approach (ssGBLUP) allows genomic prediction to take into account both genotyped and non-genotyped trees simultaneously in a single evaluation. In this study, we investigated the advantage, in terms of breeding value accuracy and bias, of including phenotypic observation from non-genotyped trees in a standard tree GBLUP evaluation. We compared the efficiency of the conventional pedigree-based (ABLUP), GBLUP and ssGBLUP approaches to evaluate eight growth and wood quality traits in a Eucalyptus hybrid population, genotyped with 33,398 single nucleotide polymorphisms (SNPs) using the EucHIP60k. Theoretical accuracies, predictive ability and bias were calculated by ten-fold cross validation on all traits. The use of additional phenotypic information from non-genotyped trees by means of ssGBLUP provided higher predictive ability (from 37% to 75%) and lower prediction bias (from 21% to 73%) for the genetic component of non-phenotyped but genotyped trees when compared to GBLUP. The increase (decrease) in the prediction accuracy (bias) became stronger as trait heritability decreased. We concluded that ssGBLUP is a promising breeding tool to improve accuracies and bias over classical GBLUP for genomic evaluation in Eucalyptus breeding practice., (Copyright © 2019 Elsevier B.V. All rights reserved.)

A thorough understanding of the heritability, genetic correlations and additive and non-additive variance components of tree growth and wood properties is a requisite for effective tree breeding. This knowledge is essential to maximize genetic gain, that is, the amount of increase in trait performance achieved annually through directional selection. Understanding the genetic attributes of traits targeted by breeding is also important to sustain decade-long genetic progress, that is, the progress made by increasing the average genetic value of the offspring as compared to that of the parental generation. In this study, we report quantitative genetic parameters for fifteen growth, wood chemical and physical traits for the world-famous Eucalyptus urograndis hybrid (E. grandis × E. urophylla). These traits directly impact the optimal use of wood for cellulose pulp, paper, and energy production. A population of 1,000 trees sampled in a progeny trial was phenotyped directly or following the development and use of near-infrared spectroscopy calibration models. Trees were genotyped with 33,398 SNPs and 24,001 DArT-seq genome-wide markers and genomic realized relationship matrices (GRM) were used for parameter estimation with an individual-tree additive-dominant mixed model. Wood chemical properties and wood density showed stronger genetic control than growth, cellulose and fiber traits. Additive effects are the main drivers of genetic variation for all traits, but dominance plays an equally or more important role for growth, singularly in this hybrid. GRM´s with >10,000 markers provided stable relationships estimates and more accurate parameters than pedigrees by capturing the full genetic relationships among individuals and disentangling the non-additive from the additive genetic component. Low correlations between growth and wood properties indicate that simultaneous selection for wood traits can be applied with minor effects on genetic gain for growth. Conversely, moderate to strong correlations between wood density and chemical traits exist, likely due to their interdependency on cell wall structure such that responses to selection will be connected for these traits. Our results illustrate the advantage of using genome-wide marker data to inform tree breeding in general and have important consequences for operational breeding of eucalypt urograndis hybrids., Competing Interests: The authors have declared that no competing interests exist. Author CG was employed by International Paper do Brasil. This commercial affiliation does not alter our adherence to all PLOS ONE policies on sharing data and materials.

Genome-wide association studies (GWAS) in plants typically suffer from limited statistical power. An alternative to the logistical and cost challenge of increasing sample sizes is to gain power by meta-analysis using information from independent studies. We carried out GWAS for growth traits with six single-marker models and regional heritability mapping (RHM) in four Eucalyptus breeding populations independently and by Joint-GWAS, using gene and segment-based models, with data for 3373 individuals genotyped with a communal EUChip60KSNP platform. While single-single nucleotide polymorphism (SNP) GWAS hardly detected significant associations at high-stringency in each population, gene-based Joint-GWAS revealed nine genes significantly associated with tree height. Associations detected using single-SNP GWAS, RHM and Joint-GWAS set-based models explained on average 3-20% of the phenotypic variance. Whole-genome regression, conversely, captured 64-89% of the pedigree-based heritability in all populations. Several associations independently detected for the same SNPs in different populations provided unprecedented GWAS validation results in forest trees. Rare and common associations were discovered in eight genes involved in cell wall biosynthesis and lignification. With the increasing adoption of genomic prediction of complex phenotypes using shared SNPs and much larger tree breeding populations, Joint-GWAS approaches should provide increasing power to pinpoint discrete associations potentially useful toward tree breeding and molecular applications., (© 2018 The Authors. New Phytologist © 2018 New Phytologist Trust.)

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus , references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

Forest tree breeding has been successful at delivering genetically improved material for multiple traits based on recurrent cycles of selection, mating, and testing. However, long breeding cycles, late flowering, variable juvenile-mature correlations, emerging pests and diseases, climate, and market changes, all pose formidable challenges. Genetic dissection approaches such as quantitative trait mapping and association genetics have been fruitless to effectively drive operational marker-assisted selection (MAS) in forest trees, largely because of the complex multifactorial inheritance of most, if not all traits of interest. The convergence of high-throughput genomics and quantitative genetics has established two new paradigms that are changing contemporary tree breeding dogmas. Genomic selection (GS) uses large number of genome-wide markers to predict complex phenotypes. It has the potential to accelerate breeding cycles, increase selection intensity and improve the accuracy of breeding values. Realized genomic relationships matrices, on the other hand, provide innovations in genetic parameters' estimation and breeding approaches by tracking the variation arising from random Mendelian segregation in pedigrees. In light of a recent flow of promising experimental results, here we briefly review the main concepts, analytical tools and remaining challenges that currently underlie the application of genomics data to tree breeding. With easy and cost-effective genotyping, we are now at the brink of extensive adoption of GS in tree breeding. Areas for future GS research include optimizing strategies for updating prediction models, adding validated functional genomics data to improve prediction accuracy, and integrating genomic and multi-environment data for forecasting the performance of genetic material in untested sites or under changing climate scenarios. The buildup of phenotypic and genome-wide data across large-scale breeding populations and advances in computational prediction of discrete genomic features should also provide opportunities to enhance the application of genomics to tree breeding.

Modern genotyping techniques, such as SNP analysis and genotyping by sequencing (GBS), are hampered by poor DNA quality and purity, particularly in challenging plant species, rich in secondary metabolites. We therefore investigated the utility of a pre-wash step using a buffered sorbitol solution, prior to DNA extraction using a high salt CTAB extraction protocol, in a high throughput or miniprep setting. This pre-wash appears to remove interfering metabolites, such as polyphenols and polysaccharides, from tissue macerates. We also investigated the adaptability of the sorbitol pre-wash for RNA extraction using a lithium chloride-based protocol. The method was successfully applied to a variety of tissues, including leaf, cambium and fruit of diverse plant species including annual crops, forest and fruit trees, herbarium leaf material and lyophilized fungal mycelium. We consistently obtained good yields of high purity DNA or RNA in all species tested. The protocol has been validated for thousands of DNA samples by generating high data quality in dense SNP arrays. DNA extracted from Eucalyptus spp. leaf and cambium as well as mycelium from Trichoderma spp. was readily digested with restriction enzymes and performed consistently in AFLP assays. Scaled-up DNA extractions were also suitable for long read sequencing. Successful RNA quality control and good RNA-Seq data for Eucalyptus and cashew confirms the effectiveness of the sorbitol buffer pre-wash for high quality RNA extraction., Competing Interests: The authors have declared that no competing interests exist.

Targeted sequence capture coupled to high-throughput sequencing has become a powerful method for the study of genome-wide sequence variation. Following our recent development of a genome assembly for the Pink Ipê tree (Handroanthus impetiginosus), a widely distributed Neotropical timber species, we now report the development of a set of 24,751 capture probes for single-nucleotide polymorphisms (SNPs) characterization and genotyping across 18,216 distinct loci, sampling more than 10 Mbp of the species genome. This system identifies nearly 200,000 SNPs located inside or in close proximity to almost 14,000 annotated protein-coding genes, generating quality genotypic data in populations spanning wide geographic distances across the species native range. To provide recommendations for future developments of similar systems for highly heterozygous plant genomes we investigated issues such as probe design, sequencing coverage and bioinformatics, including the evaluation of the capture efficiency and a reassessment of the technical reproducibility of the assay for SNPs recall and genotyping precision. Our results highlight the value of a detailed probe screening on a preliminary genome assembly to produce reliable data for downstream genetic studies. This work should inspire and assist the development of similar genomic resources for other orphan crops and forest trees with highly heterozygous genomes.

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.

Non-additive genetic effects can be effectively exploited in control-pollinated families with the availability of genome-wide markers. We used 41,304 SNP markers and compared pedigree vs. marker-based genetic models by analysing height, diameter, basic density and pulp yield for Eucalyptus urophylla × E.grandis control-pollinated families represented by 949 informative individuals. We evaluated models accounting for additive, dominance, and first-order epistatic interactions (additive by additive, dominance by dominance, and additive by dominance). We showed that the models can capture a large proportion of the genetic variance from dominance and epistasis for growth traits as those components are typically not independent. We also showed that we could partition genetic variances more precisely when using relationship matrices derived from markers compared to using only pedigree information. In addition, phenotypic prediction accuracies were only slightly increased by including dominance effects for growth traits since estimates of non-additive variances yielded rather high standard errors. This novel result improves our current understanding of the architecture of quantitative traits and recommends accounting for dominance variance when developing genomic selection strategies in hybrid Eucalyptus., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Background: Handroanthus impetiginosus (Mart. ex DC.) Mattos is a keystone Neotropical hardwood tree widely distributed in seasonally dry tropical forests of South and Mesoamerica. Regarded as the "new mahogany," it is the second most expensive timber, the most logged species in Brazil, and currently under significant illegal trading pressure. The plant produces large amounts of quinoids, specialized metabolites with documented antitumorous and antibiotic effects. The development of genomic resources is needed to better understand and conserve the diversity of the species, to empower forensic identification of the origin of timber, and to identify genes for important metabolic compounds., Findings: The genome assembly covers 503.7 Mb (N50 = 81 316 bp), 90.4% of the 557-Mbp genome, with 13 206 scaffolds. A repeat database with 1508 sequences was developed, allowing masking of ∼31% of the assembly. Depth of coverage indicated that consensus determination adequately removed haplotypes assembled separately due to the extensive heterozygosity of the species. Automatic gene prediction provided 31 688 structures and 35 479 messenger RNA transcripts, while external evidence supported a well-curated set of 28 603 high-confidence models (90% of total). Finally, we used the genomic sequence and the comprehensive gene content annotation to identify genes related to the production of specialized metabolites., Conclusions: This genome assembly is the first well-curated resource for a Neotropical forest tree and the first one for a member of the Bignoniaceae family, opening exceptional opportunities to empower molecular, phytochemical, and breeding studies. This work should inspire the development of similar genomic resources for the largely neglected forest trees of the mega-diverse tropical biomes., (© The Authors 2017. Published by Oxford University Press.)

Streptococcus pyogenes , also known as group A Streptococcus (GAS), is a human pathogen that causes diverse human diseases including streptococcal toxic shock syndrome (STSS). A GAS outbreak occurred in Brasilia, Brazil, during the second half of the year 2011, causing 26 deaths. Whole genome sequencing was performed using Illumina platform. The sequences were assembled and genes were predicted for comparative analysis with emm type 1 strains: MGAS5005 and M1 GAS. Genomics comparison revealed one of the invasive strains that differ from others isolates and from emm 1 reference genomes. Also, the new invasive strain showed differences in the content of virulence factors compared to other isolated in the same outbreak. The evolution of contemporary GAS strains is strongly associated with horizontal gene transfer. This is the first genomic study of a Streptococcal emm 1 outbreak in Brazil, and revealed the rapid bacterial evolution leading to new clones. The emergence of new invasive strains can be a consequence of the injudicious use of antibiotics in Brazil during the past decades.

Background: The advent of high-throughput genotyping technologies coupled to genomic prediction methods established a new paradigm to integrate genomics and breeding. We carried out whole-genome prediction and contrasted it to a genome-wide association study (GWAS) for growth traits in breeding populations of Eucalyptus benthamii (n =505) and Eucalyptus pellita (n =732). Both species are of increasing commercial interest for the development of germplasm adapted to environmental stresses., Results: Predictive ability reached 0.16 in E. benthamii and 0.44 in E. pellita for diameter growth. Predictive abilities using either Genomic BLUP or different Bayesian methods were similar, suggesting that growth adequately fits the infinitesimal model. Genomic prediction models using ~5000-10,000 SNPs provided predictive abilities equivalent to using all 13,787 and 19,506 SNPs genotyped in the E. benthamii and E. pellita populations, respectively. No difference was detected in predictive ability when different sets of SNPs were utilized, based on position (equidistantly genome-wide, inside genes, linkage disequilibrium pruned or on single chromosomes), as long as the total number of SNPs used was above ~5000. Predictive abilities obtained by removing relatedness between training and validation sets fell near zero for E. benthamii and were halved for E. pellita. These results corroborate the current view that relatedness is the main driver of genomic prediction, although some short-range historical linkage disequilibrium (LD) was likely captured for E. pellita. A GWAS identified only one significant association for volume growth in E. pellita, illustrating the fact that while genome-wide regression is able to account for large proportions of the heritability, very little or none of it is captured into significant associations using GWAS in breeding populations of the size evaluated in this study., Conclusions: This study provides further experimental data supporting positive prospects of using genome-wide data to capture large proportions of trait heritability and predict growth traits in trees with accuracies equal or better than those attainable by phenotypic selection. Additionally, our results document the superiority of the whole-genome regression approach in accounting for large proportions of the heritability of complex traits such as growth in contrast to the limited value of the local GWAS approach toward breeding applications in forest trees.

Background: Genomic prediction is a genomics assisted breeding methodology that can increase genetic gains by accelerating the breeding cycle and potentially improving the accuracy of breeding values. In this study, we use 41,304 informative SNPs genotyped in a Eucalyptus breeding population involving 90 E.grandis and 78 E.urophylla parents and their 949 F 1 hybrids to develop genomic prediction models for eight phenotypic traits - basic density and pulp yield, circumference at breast height and height and tree volume scored at age three and six years. We assessed the impact of different genomic prediction methods, the composition and size of the training and validation set and the number and genomic location of SNPs on the predictive ability (PA)., Results: Heritabilities estimated using the realized genomic relationship matrix (GRM) were considerably higher than estimates based on the expected pedigree, mainly due to inconsistencies in the expected pedigree that were readily corrected by the GRM. Moreover, the GRM more precisely capture Mendelian sampling among related individuals, such that the genetic covariance was based on the true proportion of the genome shared between individuals. PA improved considerably when increasing the size of the training set and by enhancing relatedness to the validation set. Prediction models trained on pure species parents could not predict well in F 1 hybrids, indicating that model training has to be carried out in hybrid populations if one is to predict in hybrid selection candidates. The different genomic prediction methods provided similar results for all traits, therefore either GBLUP or rrBLUP represents better compromises between computational time and prediction efficiency. Only slight improvement was observed in PA when more than 5000 SNPs were used for all traits. Using SNPs in intergenic regions provided slightly better PA than using SNPs sampled exclusively in genic regions., Conclusions: The size and composition of the training set and number of SNPs used are the two most important factors for model prediction, compared to the statistical methods and the genomic location of SNPs. Furthermore, training the prediction model based on pure parental species only provide limited ability to predict traits in interspecific hybrids. Our results provide additional promising perspectives for the implementation of genomic prediction in Eucalyptus breeding programs by the selection of interspecific hybrids.

Although genome-wide association studies (GWAS) have provided valuable insights into the decoding of the relationships between sequence variation and complex phenotypes, they have explained little heritability. Regional heritability mapping (RHM) provides heritability estimates for genomic segments containing both common and rare allelic effects that individually contribute too little variance to be detected by GWAS. We carried out GWAS and RHM for seven growth, wood and disease resistance traits in a breeding population of 768 Eucalyptus hybrid trees using EuCHIP60K. Total genomic heritabilities accounted for large proportions (64-89%) of pedigree-based trait heritabilities, providing additional evidence that complex traits in eucalypts are controlled by many sequence variants across the frequency spectrum, each with small contributions to the phenotypic variance. RHM detected 26 quantitative trait loci (QTLs) encompassing 2191 single nucleotide polymorphisms (SNPs), whereas GWAS detected 13 single SNP-trait associations. RHM and GWAS QTLs individually explained 5-15% and 4-6% of the genomic heritability, respectively. RHM was superior to GWAS in capturing larger proportions of genomic heritability. Equated to previously mapped QTLs, our results highlighted genomic regions for further examination towards gene discovery. RHM-QTLs bearing a combination of common and rare variants could be useful enhancements to incorporate prior knowledge of the underlying genetic architecture in genomic prediction models., (© 2016 The Authors. New Phytologist © 2016 New Phytologist Trust.)

The genome-wide heterozygosity at 9590 genes, all heterozygous in a single Eucalyptus grandis parent tree, was examined in a group of 28 S1 offspring. Heterozygosity ranged from 52-79%, averaging 65.5%, much higher than the 50% expected under random segregation, supporting the occurrence of strong (47%) selection against homozygosity. The expected pattern of heterozygosity from theoretical calculations and simulations for recessive detrimentals (pseudo-overdominance) and intrinsic heterozygote advantage was examined and compared with that observed. The observed patterns are consistent with at least several detrimental loci with large effects on both parental chromosomes of the 11 pairs. It is likely that 100 or more genes, many with substantial effects on viability, are contributing to this inbreeding depression. Although our genome-wide analysis of nearly 10 000 genes strongly suggested that pseudo-overdominance was responsible for the observed high inbreeding depression, heterozygote advantage could not be excluded. Finding inconvertible evidence of the cause of inbreeding depression still presents a difficult challenge. This study is the first theoretical examination of the genomic effect of inbreeding in a forest tree and provides an approach to analyze these data to determine the extent and cause of inbreeding depression across other plant genomes., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

Background: Micro RNAs are a class of small non coding RNAs of 20-24 nucleotides transcribed as single stranded precursors from MIR gene loci. Initially described as post-transcriptional regulators involved in development, two decades ago, miRNAs have been proven to regulate a wide range of processes in plants such as germination, morphology and responses to biotic and abiotic stress. Despite wide conservation in plants, a number of miRNAs are lineage specific. We describe the first genome wide survey of Eucalyptus miRNAs based on high throughput sequencing., Results: In addition to discovering small RNA sequences, MIR loci were mapped onto the reference genome and interspecific variability investigated. Sequencing was carried out for the two most world widely planted species, E. grandis and E. globulus. To maximize discovery, E. grandis samples were from BRASUZ1, the same tree whose genome provided the reference sequence. Interspecific analysis reinforces the variability in small RNA repertoire even between closely related species. Characterization of Eucalyptus small RNA sequences showed 95 orthologous to conserved miRNAs and 193 novel miRNAs. In silico target prediction confirmed 163 novel miRNAs and degradome sequencing experimentally confirmed several hundred targets. Experimental evidence based on the exclusive expression of a set of small RNAs across 16 species within Myrtaceae further highlighted variable patterns of conservation and diversity of these regulatory elements., Conclusions: The description of miRNAs in Eucalyptus contributes to scientific knowledge of this vast genre, which is the most widely planted hardwood crop in the tropical and subtropical world, adding another important element to the annotation of Eucalyptus grandis reference genome.

We used high-density single nucleotide polymorphism (SNP) data and whole-genome pooled resequencing to examine the landscape of population recombination (ρ) and nucleotide diversity (ϴw ), assess the extent of linkage disequilibrium (r(2) ) and build the highest density linkage maps for Eucalyptus. At the genome-wide level, linkage disequilibrium (LD) decayed within c. 4-6 kb, slower than previously reported from candidate gene studies, but showing considerable variation from absence to complete LD up to 50 kb. A sharp decrease in the estimate of ρ was seen when going from short to genome-wide inter-SNP distances, highlighting the dependence of this parameter on the scale of observation adopted. Recombination was correlated with nucleotide diversity, gene density and distance from the centromere, with hotspots of recombination enriched for genes involved in chemical reactions and pathways of the normal metabolic processes. The high nucleotide diversity (ϴw = 0.022) of E. grandis revealed that mutation is more important than recombination in shaping its genomic diversity (ρ/ϴw = 0.645). Chromosome-wide ancestral recombination graphs allowed us to date the split of E. grandis (1.7-4.8 million yr ago) and identify a scenario for the recent demographic history of the species. Our results have considerable practical importance to Genome Wide Association Studies (GWAS), while indicating bright prospects for genomic prediction of complex phenotypes in eucalypt breeding., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

Pedigree reconstruction using molecular markers enables efficient management of inbreeding in open-pollinated breeding strategies, replacing expensive and time-consuming controlled pollination. This is particularly useful in preferentially outcrossed, insect pollinated Eucalypts known to suffer considerable inbreeding depression from related matings. A single nucleotide polymorphism (SNP) marker panel consisting of 106 markers was selected for pedigree reconstruction from the recently developed high-density Eucalyptus Infinium SNP chip (EuCHIP60K). The performance of this SNP panel for pedigree reconstruction in open-pollinated progenies of two Eucalyptus nitens seed orchards was compared with that of two microsatellite panels with 13 and 16 markers respectively. The SNP marker panel out-performed one of the microsatellite panels in the resolution power to reconstruct pedigrees and out-performed both panels with respect to data quality. Parentage of all but one offspring in each clonal seed orchard was correctly matched to the expected seed parent using the SNP marker panel, whereas parentage assignment to less than a third of the expected seed parents were supported using the 13-microsatellite panel. The 16-microsatellite panel supported all but one of the recorded seed parents, one better than the SNP panel, although there was still a considerable level of missing and inconsistent data. SNP marker data was considerably superior to microsatellite data in accuracy, reproducibility and robustness. Although microsatellites and SNPs data provide equivalent resolution for pedigree reconstruction, microsatellite analysis requires more time and experience to deal with the uncertainties of allele calling and faces challenges for data transferability across labs and over time. While microsatellite analysis will continue to be useful for some breeding tasks due to the high information content, existing infrastructure and low operating costs, the multi-species SNP resource available with the EuCHIP60k, opens a whole new array of opportunities for high-throughput, genome-wide or targeted genotyping in species of Eucalyptus.

We used whole genome resequencing of pooled individuals to develop a high-density single-nucleotide polymorphism (SNP) chip for Eucalyptus. Genomes of 240 trees of 12 species were sequenced at 3.5× each, and 46 997 586 raw SNP variants were subject to multivariable filtering metrics toward a multispecies, genome-wide distributed chip content. Of the 60 904 SNPs on the chip, 59 222 were genotyped and 51 204 were polymorphic across 14 Eucalyptus species, providing a 96% genome-wide coverage with 1 SNP/12-20 kb, and 47 069 SNPs at ≤ 10 kb from 30 444 of the 33 917 genes in the Eucalyptus genome. Given the EUChip60K multi-species genotyping flexibility, we show that both the sample size and taxonomic composition of cluster files impact heterozygous call specificity and sensitivity by benchmarking against 'gold standard' genotypes derived from deeply sequenced individual tree genomes. Thousands of SNPs were shared across species, likely representing ancient variants arisen before the split of these taxa, hinting to a recent eucalypt radiation. We show that the variable SNP filtering constraints allowed coverage of the entire site frequency spectrum, mitigating SNP ascertainment bias. The EUChip60K represents an outstanding tool with which to address population genomics questions in Eucalyptus and to empower genomic selection, GWAS and the broader study of complex trait variation in eucalypts., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

Keystone species in their native ranges, eucalypts, are ecologically and genetically very diverse, growing naturally along extensive latitudinal and altitudinal ranges and variable environments. Besides their ecological importance, eucalypts are also the most widely planted trees for sustainable forestry in the world. We report the development of a novel collection of 535 microsatellites for species of Eucalyptus, 494 designed from ESTs and 41 from genomic libraries. A selected subset of 223 was evaluated for individual identification, parentage testing, and ancestral information content in the two most extensively studied species, Eucalyptus grandis and Eucalyptus globulus. Microsatellites showed high transferability and overlapping allele size range, suggesting they have arisen still in their common ancestor and confirming the extensive genome conservation between these two species. A consensus linkage map with 437 microsatellites, the most comprehensive microsatellite-only genetic map for Eucalyptus, was built by assembling segregation data from three mapping populations and anchored to the Eucalyptus genome. An overall colinearity between recombination-based and physical positioning of 84% of the mapped microsatellites was observed, with some ordering discrepancies and sporadic locus duplications, consistent with the recently described whole genome duplication events in Eucalyptus. The linkage map covered 95.2% of the 605.8-Mbp assembled genome sequence, placing one microsatellite every 1.55 Mbp on average, and an overall estimate of physical to recombination distance of 618 kbp/cM. The genetic parameters estimates together with linkage and physical position data for this large set of microsatellites should assist marker choice for genome-wide population genetics and comparative mapping in Eucalyptus., (© 2014 John Wiley & Sons Ltd.)

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

The promise of association genetics to identify genes or genomic regions controlling complex traits has generated a flurry of interest. Such phenotype-genotype associations could be useful to accelerate tree breeding cycles, increase precision and selection intensity for late expressing, low heritability traits. However, the prospects of association genetics in highly heterozygous undomesticated forest trees can be severely impacted by the presence of cryptic population and pedigree structure. To investigate how to better account for this, we compared the GLM and five combinations of the Unified Mixed Model ( UMM ) on data of a low-density genome-wide association study for growth and wood property traits carried out in a Eucalyptus globulus population (n = 303) with 7,680 Diversity Array Technology (DArT) markers. Model comparisons were based on the degree of deviation from the uniform distribution and estimates of the mean square differences between the observed and expected p-values of all significant marker-trait associations detected. Our analysis revealed the presence of population and family structure. There was not a single best model for all traits. Striking differences in detection power and accuracy were observed among the different models especially when population structure was not accounted for. The UMM method was the best and produced superior results when compared to GLM for all traits. Following stringent correction for false discoveries, 18 marker-trait associations were detected, 16 for tree diameter growth and two for lignin monomer composition (S:G ratio), a key wood property trait. The two DArT markers associated with S:G ratio on chromosome 10, physically map within 1 Mbp of the ferulate 5-hydroxylase (F5H) gene, providing a putative independent validation of this marker-trait association. This study details the merit of collectively integrate population structure and relatedness in association analyses in undomesticated, highly heterozygous forest trees, and provides additional insights into the nature of complex quantitative traits in 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.

• Genomic selection (GS) is expected to cause a paradigm shift in tree breeding by improving its speed and efficiency. By fitting all the genome-wide markers concurrently, GS can capture most of the 'missing heritability' of complex traits that quantitative trait locus (QTL) and association mapping classically fail to explain. Experimental support of GS is now required. • The effectiveness of GS was assessed in two unrelated Eucalyptus breeding populations with contrasting effective population sizes (N(e) = 11 and 51) genotyped with > 3000 DArT markers. Prediction models were developed for tree circumference and height growth, wood specific gravity and pulp yield using random regression best linear unbiased predictor (BLUP). • Accuracies of GS varied between 0.55 and 0.88, matching the accuracies achieved by conventional phenotypic selection. Substantial proportions (74-97%) of trait heritability were captured by fitting all genome-wide markers simultaneously. Genomic regions explaining trait variation largely coincided between populations, although GS models predicted poorly across populations, likely as a result of variable patterns of linkage disequilibrium, inconsistent allelic effects and genotype × environment interaction. • GS brings a new perspective to the understanding of quantitative trait variation in forest trees and provides a revolutionary tool for applied tree improvement. Nevertheless population-specific predictive models will likely drive the initial applications of GS in forest tree breeding., (© 2012 The Authors. New Phytologist © 2012 New Phytologist Trust.)

Diversity Arrays Technology (DArT) provides a robust, high throughput, cost-effective method to query thousands of sequence polymorphisms in a single assay. Despite the extensive use of this genotyping platform for numerous plant species, little is known regarding the sequence attributes and genome-wide distribution of DArT markers. We investigated the genomic properties of the 7,680 DArT marker probes of a Eucalyptus array, by sequencing them, constructing a high density linkage map and carrying out detailed physical mapping analyses to the Eucalyptus grandis reference genome. A consensus linkage map with 2,274 DArT markers anchored to 210 microsatellites and a framework map, with improved support for ordering, displayed extensive collinearity with the genome sequence. Only 1.4 Mbp of the 75 Mbp of still unplaced scaffold sequence was captured by 45 linkage mapped but physically unaligned markers to the 11 main Eucalyptus pseudochromosomes, providing compelling evidence for the quality and completeness of the current Eucalyptus genome assembly. A highly significant correspondence was found between the locations of DArT markers and predicted gene models, while most of the 89 DArT probes unaligned to the genome correspond to sequences likely absent in E. grandis, consistent with the pan-genomic feature of this multi-Eucalyptus species DArT array. These comprehensive linkage-to-physical mapping analyses provide novel data regarding the genomic attributes of DArT markers in plant genomes in general and for Eucalyptus in particular. DArT markers preferentially target the gene space and display a largely homogeneous distribution across the genome, thereby providing superb coverage for mapping and genome-wide applications in breeding and diversity studies. Data reported on these ubiquitous properties of DArT markers will be particularly valuable to researchers working on less-studied crop species who already count on DArT genotyping arrays but for which no reference genome is yet available to allow such detailed characterization.

The vitamin D receptor (VDR) is an essential protein related to bone metabolism. Some VDR alleles are differentially distributed among ethnic populations and display variable patterns of linkage disequilibrium (LD). In this study, 200 unrelated Brazilians were genotyped using 21 VDR single nucleotide polymorphisms (SNPs) and 28 ancestry informative markers. The patterns of LD and haplotype distribution were compared among Brazilian and the HapMap populations of African (YRI), European (CEU) and Asian (JPT+CHB) origins. Conditional regression and haplotype-specific analysis were performed using estimates of individual genetic ancestry in Brazilians as a quantitative trait. Similar patterns of LD were observed in the 5' and 3' gene regions. However, the frequency distribution of haplotype blocks varied among populations. Conditional regression analysis identified haplotypes associated with European and Amerindian ancestry, but not with the proportion of African ancestry. Individual ancestry estimates were associated with VDR haplotypes. These findings reinforce the need to correct for population stratification when performing genetic association studies in admixed populations.

Background: High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera., Results: We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of Eucalyptus from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for E. grandis. A systematic assessment of in silico SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous in silico constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species.SNP reliability was high across nine Eucalyptus species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased., Conclusions: This study indicates that the GGGT performs well both within and across species of Eucalyptus notwithstanding its nucleotide diversity ≥ 2%. The development of a much larger array of informative SNPs across multiple Eucalyptus species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in Eucalyptus.

Background: Technological advances are progressively increasing the application of genomics to a wider array of economically and ecologically important species. High-density maps enriched for transcribed genes facilitate the discovery of connections between genes and phenotypes. We report the construction of a high-density linkage map of expressed genes for the heterozygous genome of Eucalyptus using Single Feature Polymorphism (SFP) markers., Results: SFP discovery and mapping was achieved using pseudo-testcross screening and selective mapping to simultaneously optimize linkage mapping and microarray costs. SFP genotyping was carried out by hybridizing complementary RNA prepared from 4.5 year-old trees xylem to an SFP array containing 103,000 25-mer oligonucleotide probes representing 20,726 unigenes derived from a modest size expressed sequence tags collection. An SFP-mapping microarray with 43,777 selected candidate SFP probes representing 15,698 genes was subsequently designed and used to genotype SFPs in a larger subset of the segregating population drawn by selective mapping. A total of 1,845 genes were mapped, with 884 of them ordered with high likelihood support on a framework map anchored to 180 microsatellites with average density of 1.2 cM. Using more probes per unigene increased by two-fold the likelihood of detecting segregating SFPs eventually resulting in more genes mapped. In silico validation showed that 87% of the SFPs map to the expected location on the 4.5X draft sequence of the Eucalyptus grandis genome., Conclusions: The Eucalyptus 1,845 gene map is the most highly enriched map for transcriptional information for any forest tree species to date. It represents a major improvement on the number of genes previously positioned on Eucalyptus maps and provides an initial glimpse at the gene space for this global tree genome. A general protocol is proposed to build high-density transcript linkage maps in less characterized plant species by SFP genotyping with a concurrent objective of reducing microarray costs. HIgh-density gene-rich maps represent a powerful resource to assist gene discovery endeavors when used in combination with QTL and association mapping and should be especially valuable to assist the assembly of reference genome sequences soon to come for several plant and animal species.

A set of over 8000 Diversity Arrays Technology (DArT) markers was tested for its utility in high-resolution population and phylogenetic studies across a range of Eucalyptus taxa. Small-scale population studies of Eucalyptus camaldulensis, Eucalyptus cladocalyx, Eucalyptus globulus, Eucalyptus grandis, Eucalyptus nitens, Eucalyptus pilularis and Eucalyptus urophylla demonstrated the potential of genome-wide genotyping with DArT markers to differentiate species, to identify interspecific hybrids and to resolve biogeographic disjunctions within species. The population genetic studies resolved geographically partitioned clusters in E. camaldulensis, E. cladocalyx, E. globulus and E. urophylla that were congruent with previous molecular studies. A phylogenetic study of 94 eucalypt species provided results that were largely congruent with traditional taxonomy and ITS-based phylogenies, but provided more resolution within major clades than had been obtained previously. Ascertainment bias (the bias introduced in a phylogeny from using markers developed in a small sample of the taxa that are being studied) was not detected. DArT offers an unprecedented level of resolution for population genetic, phylogenetic and evolutionary studies across the full range of Eucalyptus species., (Copyright © 2011. Published by Elsevier Inc.)

Background: Eucalyptus species are among the most planted hardwoods in the world because of their rapid growth, adaptability and valuable wood properties. The development and integration of genomic resources into breeding practice will be increasingly important in the decades to come. Bacterial artificial chromosome (BAC) libraries are key genomic tools that enable positional cloning of important traits, synteny evaluation, and the development of genome framework physical maps for genetic linkage and genome sequencing., Results: We describe the construction and characterization of two deep-coverage BAC libraries EG&#95;Ba and EG&#95;Bb obtained from nuclear DNA fragments of E. grandis (clone BRASUZ1) digested with HindIII and BstYI, respectively. Genome coverages of 17 and 15 haploid genome equivalents were estimated for EG&#95;Ba and EG&#95;Bb, respectively. Both libraries contained large inserts, with average sizes ranging from 135 Kb (Eg&#95;Bb) to 157 Kb (Eg&#95;Ba), very low extra-nuclear genome contamination providing a probability of finding a single copy gene ≥ 99.99%. Libraries were screened for the presence of several genes of interest via hybridizations to high-density BAC filters followed by PCR validation. Five selected BAC clones were sequenced and assembled using the Roche GS FLX technology providing the whole sequence of the E. grandis chloroplast genome, and complete genomic sequences of important lignin biosynthesis genes., Conclusions: The two E. grandis BAC libraries described in this study represent an important milestone for the advancement of Eucalyptus genomics and forest tree research. These BAC resources have a highly redundant genome coverage (> 15×), contain large average inserts and have a very low percentage of clones with organellar DNA or empty vectors. These publicly available BAC libraries are thus suitable for a broad range of applications in genetic and genomic research in Eucalyptus and possibly in related species of Myrtaceae, including genome sequencing, gene isolation, functional and comparative genomics. Because they have been constructed using the same tree (E. grandis BRASUZ1) whose full genome is being sequenced, they should prove instrumental for assembly and gap filling of the upcoming Eucalyptus reference genome sequence.

Species of Eucalyptus are keystone species for ecological studies in their natural ranges and are extensively planted in the tropical and subtropical regions of the world to supply high-quality woody biomass for various applications. We report the development of a selected set of 20 dinucleotide and trinucleotide repeat microsatellites derived from Eucalyptus expressed sequence tags (ESTs). These microsatellites were selected for full transferability and homogeneous rate of polymorphism across species. They were evaluated for individual fingerprinting, parentage testing, and intraspecific population structure analyses in 6 of the most extensively studied and planted species worldwide, representing key phylogenetic sections of the largest subgenus Symphyomyrtus. This set of markers provides exceptional resolution for population genetics and molecular breeding applications in the genus Eucalyptus. As they were developed from conserved transcribed regions, the transferability and polymorphism of these microsatellites will most likely extend to the other 300 or more species within the same subgenus.

Background: A number of molecular marker technologies have allowed important advances in the understanding of the genetics and evolution of Eucalyptus, a genus that includes over 700 species, some of which are used worldwide in plantation forestry. Nevertheless, the average marker density achieved with current technologies remains at the level of a few hundred markers per population. Furthermore, the transferability of markers produced with most existing technology across species and pedigrees is usually very limited. High throughput, combined with wide genome coverage and high transferability are necessary to increase the resolution, speed and utility of molecular marker technology in eucalypts. We report the development of a high-density DArT genome profiling resource and demonstrate its potential for genome-wide diversity analysis and linkage mapping in several species of Eucalyptus., Findings: After testing several genome complexity reduction methods we identified the PstI/TaqI method as the most effective for Eucalyptus and developed 18 genomic libraries from PstI/TaqI representations of 64 different Eucalyptus species. A total of 23,808 cloned DNA fragments were screened and 13,300 (56%) were found to be polymorphic among 284 individuals. After a redundancy analysis, 6,528 markers were selected for the operational array and these were supplemented with 1,152 additional clones taken from a library made from the E. grandis tree whose genome has been sequenced. Performance validation for diversity studies revealed 4,752 polymorphic markers among 174 individuals. Additionally, 5,013 markers showed segregation when screened using six inter-specific mapping pedigrees, with an average of 2,211 polymorphic markers per pedigree and a minimum of 859 polymorphic markers that were shared between any two pedigrees., Conclusions: This operational DArT array will deliver 1,000-2,000 polymorphic markers for linkage mapping in most eucalypt pedigrees and thus provide high genome coverage. This array will also provide a high-throughput platform for population genetics and phylogenetics in Eucalyptus. The transferability of DArT across species and pedigrees is particularly valuable for a large genus such as Eucalyptus and will facilitate the transfer of information between different studies. Furthermore, the DArT marker array will provide a high-resolution link between phenotypes in populations and the Eucalyptus reference genome, which will soon be completed.

The rs2476601-T allele at the protein tyrosine phosphatase non-receptor type 22 (PTPN22) gene has been consistently associated with several autoimmune diseases in European-derived populations. However, little is known about the allele and haplotype frequency distributions in PTPN22 among populations derived from other ethnic groups. In the present study, the allele and haplotype frequency distributions of six single nucleotide polymorphisms (SNPs) in the PTPN22 gene were compared among Brazilian populations and HapMap phase 3 dataset. A total of 10 different population samples were evaluated. Additionally, in admixed populations, individual genetic ancestries were estimated for Native American, African, and European contributions. Estimated individual ancestries were used as quantitative traits in a conditional approach for single-marker and haplotype-specific regression analyses. It was shown that several SNPs and haplotypes have different frequencies among different ethnic populations. Individual genetic ancestries were not associated with the rs2476601-T allele, but were associated with PTPN22 haplotypes in Brazilian, Mexican, and African-American admixed populations. Our results suggest caution in the interpretation of results found in association studies involving PTPN22 polymorphisms in admixed populations. Correction for stratification generated by admixture should be mandatory to minimize or avoid chances of spurious association.

Ancestry informative SNPs can be useful to estimate individual and population biogeographical ancestry. Brazilian population is characterized by a genetic background of three parental populations (European, African, and Brazilian Native Amerindians) with a wide degree and diverse patterns of admixture. In this work we analyzed the information content of 28 ancestry-informative SNPs into multiplexed panels using three parental population sources (African, Amerindian, and European) to infer the genetic admixture in an urban sample of the five Brazilian geopolitical regions. The SNPs assigned apart the parental populations from each other and thus can be applied for ancestry estimation in a three hybrid admixed population. Data was used to infer genetic ancestry in Brazilians with an admixture model. Pairwise estimates of F(st) among the five Brazilian geopolitical regions suggested little genetic differentiation only between the South and the remaining regions. Estimates of ancestry results are consistent with the heterogeneous genetic profile of Brazilian population, with a major contribution of European ancestry (0.771) followed by African (0.143) and Amerindian contributions (0.085). The described multiplexed SNP panels can be useful tool for bioanthropological studies but it can be mainly valuable to control for spurious results in genetic association studies in admixed populations., ((c) 2009 Wiley-Liss, Inc.)