A simulation study comparing advanced marker-assisted selection with genomic selection in tree breeding programs.

Advances in DNA sequencing technologies allow the sequencing of whole genomes of thousands of individuals and provide several million single nucleotide polymorphisms (SNPs) per individual. These data combined with precise and high-throughput phenotyping enable genome-wide association studies (GWAS) and the identification of SNPs underlying traits with complex genetic architectures. The identified causal SNPs and estimated allelic effects could then be used for advanced marker-assisted selection (MAS) in breeding programs. But could such MAS compete with the broadly used genomic selection (GS)? This question is of particular interest for the lengthy tree breeding strategies. Here, with our new software "SNPscan breeder," we simulated a simple tree breeding program and compared the impact of different selection criteria on genetic gain and inbreeding. Further, we assessed different genetic architectures and different levels of kinship among individuals of the breeding population. Interestingly, apart from progeny testing, GS using gBLUP performed best under almost all simulated scenarios. MAS based on GWAS results outperformed GS only if the allelic effects were estimated in large populations (ca. 10,000 individuals) of unrelated individuals. Notably, GWAS using 3,000 extreme phenotypes performed as good as the use of 10,000 phenotypes. GS increased inbreeding and thus reduced genetic diversity more strongly compared to progeny testing and GWAS-based selection. We discuss the practical implications for tree breeding programs. In conclusion, our analyses further support the potential of GS for forest tree breeding and improvement, although MAS may gain relevance with decreasing sequencing costs in the future.

SNP Markers as a Successful Molecular Tool for Assessing Species Identity and Geographic Origin of Trees in the Economically Important South American Legume Genus Dipteryx.

Dipteryx timber has been heavily exploited in South America since 2000s due to the increasing international demand for hardwood. Developing tools for the genetic identification of Dipteryx species and their geographical origin can help to promote legal trading of timber. A collection of 800 individual trees, belonging to 6 different Dipteryx species, was genotyped based on 171 molecular markers. After the exclusion of markers out of Hardy-Weinberg equilibrium or with no polymorphism or low amplification, 83 nuclear, 29 chloroplast, 13 mitochondrial single nucleotide polymorphisms (SNPs), and 2 chloroplast and 5 mitochondrial INDELS remained. Six genetic groups were identified using Bayesian Structure analyses of the nuclear SNPs, which corresponded to the different Dipteryx species collected in the field. Seventeen highly informative markers were identified as suitable for species identification and obtained self-assignment success rates to species level of 78-96%. An additional set of 15 molecular markers was selected to determine the different genetic clusters found in Dipteryx odorata and Dipteryx ferrea, obtaining self-assignment success rates of 91-100%. The success to assign samples to the correct country of origin using all or only the informative markers improved when using the nearest neighbor approach (69-92%) compared to the Bayesian approach (33-80%). While nuclear and chloroplast SNPs were more suitable for differentiating the different Dipteryx species, mitochondrial SNPs were ideal for determining the genetic clusters of D. odorata and D. ferrea. These 32 selected SNPs will be invaluable genetic tools for the accurate identification of species and country of origin of Dipteryx timber.

Molecular evidence for three genetic species of Dipteryx in the Peruvian Amazon.

There is a high international demand for timber from the genus Dipteryx, or "shihuahuaco" as it is known in Peru. Developing tools that allow the identification and discrimination of Dipteryx species is therefore important for supporting management of natural populations and to underpin legal trade of its timber. The objective of this study was the molecular characterization of Dipteryx species in the Peruvian Amazonia. Two plastid regions (cpDNA: trnH-psbA and matK) were sequenced and 11 microsatellite markers (nDNA) were genotyped for 32 individuals identified as Dipteryx charapilla, D. micrantha morphotype 1 and D. micrantha morphotype 2. Using the concatenated sequences of the plastid genes, we identified ten haplotypes that were not shared between the species or between the D. micrantha morphotypes. Haplotypic diversity was greater in D. micrantha morphotype 2 and D. charapilla than in D. micrantha morphotype 1, which presented only one haplotype with a wide distribution in Peru. The microsatellites allowed the discrimination of the same three clades and identified diagnostic alleles for each clade. These results allowed us to demonstrate that the two morphotypes of D. micrantha are different at both the plastid and nuclear markers, which supports the existence of three genetically distinct species in Peru. This study provides information for the genetic discrimination of Dipteryx species and emphasises the importance of conserving the genetic variability of this genus in the Peruvian Amazonia.

Assessing the Ability of Chloroplast and Nuclear DNA Gene Markers to Verify the Geographic Origin of Jatoba (Hymenaea courbaril L.) Timber.

Deforestation-reinforced by illegal logging-is a serious problem in many tropical regions and causes pervasive environmental and economic damage. Existing laws that intend to reduce illegal logging need efficient, fraud resistant control methods. We developed a genetic reference database for Jatoba (Hymenaea courbaril), an important, high value timber species from the Neotropics. The data set can be used for controls on declarations of wood origin. Samples from 308 Hymenaea trees from 12 locations in Brazil, Bolivia, Peru, and French Guiana have been collected and genotyped on 10 nuclear microsatellites (nSSRs), 13 chloroplast SNPs (cpSNP), and 1 chloroplast indel marker. The chloroplast gene markers have been developed using Illumina DNA sequencing. Bayesian cluster analysis divided the individuals based on the nSSRs into 8 genetic groups. Using self-assignment tests, the power of the genetic reference database to judge on declarations on the location has been tested for 3 different assignment methods. We observed a strong genetic differentiation among locations leading to high and reliable self-assignment rates for the locations between 50% to 100% (average of 88%). Although all 3 assignment methods came up with similar mean self-assignment rates, there were differences for some locations linked to the level of genetic diversity, differentiation, and heterozygosity. Our results show that the nuclear and chloroplast gene markers are effective to be used for a genetic certification system and can provide national and international authorities with a robust tool to confirm legality of timber.

Complete Chloroplast Genome Sequences of Four Meliaceae Species and Comparative Analyses.

The Meliaceae family mainly consists of trees and shrubs with a pantropical distribution. In this study, the complete chloroplast genomes of four Meliaceae species were sequenced and compared with each other and with the previously published Azadirachta indica plastome. The five plastomes are circular and exhibit a quadripartite structure with high conservation of gene content and order. They include 130 genes encoding 85 proteins, 37 tRNAs and 8 rRNAs. Inverted repeat expansion resulted in a duplication of rps19 in the five Meliaceae species, which is consistent with that in many other Sapindales, but different from many other rosids. Compared to Azadirachta indica, the four newly sequenced Meliaceae individuals share several large deletions, which mainly contribute to the decreased genome sizes. A whole-plastome phylogeny supports previous findings that the four species form a monophyletic sister clade to Azadirachta indica within the Meliaceae. SNPs and indels identified in all complete Meliaceae plastomes might be suitable targets for the future development of genetic markers at different taxonomic levels. The extended analysis of SNPs in the matK gene led to the identification of four potential Meliaceae-specific SNPs as a basis for future validation and marker development.

A nearest neighbour approach by genetic distance to the assignment of individual trees to geographic origin.

During the past decade, the use of DNA for forensic applications has been extensively implemented for plant and animal species, as well as in humans. Tracing back the geographical origin of an individual usually requires genetic assignment analysis. These approaches are based on reference samples that are grouped into populations or other aggregates and intend to identify the most likely group of origin. Often this grouping does not have a biological but rather a historical or political justification, such as "country of origin". In this paper, we present a new nearest neighbour approach to individual assignment or classification within a given but potentially imperfect grouping of reference samples. This method, which is based on the genetic distance between individuals, functions better in many cases than commonly used methods. We demonstrate the operation of our assignment method using two data sets. One set is simulated for a large number of trees distributed in a 120km by 120km landscape with individual genotypes at 150 SNPs, and the other set comprises experimental data of 1221 individuals of the African tropical tree species Entandrophragma cylindricum (Sapelli) genotyped at 61 SNPs. Judging by the level of correct self-assignment, our approach outperformed the commonly used frequency and Bayesian approaches by 15% for the simulated data set and by 5-7% for the Sapelli data set. Our new approach is less sensitive to overlapping sources of genetic differentiation, such as genetic differences among closely-related species, phylogeographic lineages and isolation by distance, and thus operates better even for suboptimal grouping of individuals.

Revealing hidden species diversity in closely related species using nuclear SNPs, SSRs and DNA sequences - a case study in the tree genus Milicia.

BACKGROUND: Species delimitation in closely related plant taxa can be challenging because (i) reproductive barriers are not always congruent with morphological differentiation, (ii) use of plastid sequences might lead to misinterpretation, (iii) rare species might not be sampled. We revisited molecular-based species delimitation in the African genus Milicia, currently divided into M. regia (West Africa) and M. excelsa (from West to East Africa). We used 435 samples collected in West, Central and East Africa. We genotyped SNP and SSR loci to identify genetic clusters, and sequenced two plastid regions (psbA-trnH, trnC-ycf6) and a nuclear gene (At103) to confirm species' divergence and compare species delimitation methods. We also examined whether ecological niche differentiation was congruent with sampled genetic structure. RESULTS: West African M. regia, West African and East African M. excelsa samples constituted three well distinct genetic clusters according to SNPs and SSRs. In Central Africa, two genetic clusters were consistently inferred by both types of markers, while a few scattered samples, sympatric with the preceding clusters but exhibiting leaf traits of M. regia, were grouped with the West African M. regia cluster based on SNPs or formed a distinct cluster based on SSRs. SSR results were confirmed by sequence data from the nuclear region At103 which revealed three distinct 'Fields For Recombination' corresponding to (i) West African M. regia, (ii) Central African samples with leaf traits of M. regia, and (iii) all M. excelsa samples. None of the plastid sequences provide indication of distinct clades of the three species-like units. Niche modelling techniques yielded a significant correlation between niche overlap and genetic distance. CONCLUSIONS: Our genetic data suggest that three species of Milicia could be recognized. It is surprising that the occurrence of two species in Central Africa was not reported for this well-known timber tree. Globally, our work highlights the importance of collecting samples in a systematic way and the need for combining different nuclear markers when dealing with species complexes. Recognizing cryptic species is particularly crucial for economically exploited species because some hidden taxa might actually be endangered as they are merged with more abundant species.

Development of Molecular Markers for Determining Continental Origin of Wood from White Oaks (Quercus L. sect. Quercus).

To detect and avoid illegal logging of valuable tree species, identification methods for the origin of timber are necessary. We used next-generation sequencing to identify chloroplast genome regions that differentiate the origin of white oaks from the three continents; Asia, Europe, and North America. By using the chloroplast genome of Asian Q. mongolica as a reference, we identified 861 variant sites (672 single nucleotide polymorphisms (SNPs); 189 insertion/deletion (indel) polymorphism) from representative species of three continents (Q. mongolica from Asia; Q. petraea and Q. robur from Europe; Q. alba from North America), and we identified additional chloroplast polymorphisms in pools of 20 individuals each from Q. mongolica (789 variant sites) and Q. robur (346 variant sites). Genome sequences were screened for indels to develop markers that identify continental origin of oak species, and that can be easily evaluated using a variety of detection methods. We identified five indels and one SNP that reliably identify continent-of-origin, based on evaluations of up to 1078 individuals representing 13 white oak species and three continents. Due to the size of length polymorphisms revealed, this marker set can be visualized using capillary electrophoresis or high resolution gel (acrylamide or agarose) electrophoresis. With these markers, we provide the wood trading market with an instrument to comply with the U.S. and European laws that require timber companies to avoid the trade of illegally harvested timber.

Efficient long-distance gene flow into an isolated relict oak stand.

Geographically isolated and small populations outside a species' central distribution range are likely to be of major importance to a species' ability to quickly adjust its distribution range to global change dynamics. Gene flow from the outside plays a pivotal role in the fate of these marginal populations. It has been proposed that spatial fragmentation and perceived geographic isolation do not necessarily reflect a loss of genetic connectivity in tree species. However, the spatial limits of long-distance gene flow, as well as its magnitude and impact, are still generally unknown. In the present study, we analyzed long-distance pollen-mediated gene flow into an isolated relict stand consisting of 7 individuals of Quercus robur L. based on a total sample of 177 trees and 9 microsatellite loci. We show that pollen-mediated gene flow across more than 80 km in this wind-pollinated tree species contributed at least 35% of all successful pollinations in the investigated isolated and small oak stand at the eastern limit of the species' distribution. The observed pollen immigration shaped the genetic diversity of acorn progenies in the stand and might explain the comparably high genetic diversity in the persisting adult population.

A very small and isolated population of the green Oak Leaf Roller, Tortrix viridana L., with high genetic diversity--how does this work?

A geographically isolated population of an herbivorous insect (Tortrix viridana, Lepidoptera: Tortricidae) in the Bashkir Transural region and 5 further populations were investigated for genetic variation using eight microsatellite markers. The sample size per population was between 48 and 62 individuals. The genetic variation was higher within the isolated population than within populations in the center of the distribution area. No bottleneck effects were discovered during analyses that could have formatted the gene pool of this population. Balancing or directed selection toward preservation of specific alleles or promotion of heterozygous individuals could be an explanation for the unexpected high genetic diversity within this small and isolated population.

Fine-scale genetic structure and gene dispersal inferences in 10 neotropical tree species.

The extent of gene dispersal is a fundamental factor of the population and evolutionary dynamics of tropical tree species, but directly monitoring seed and pollen movement is a difficult task. However, indirect estimates of historical gene dispersal can be obtained from the fine-scale spatial genetic structure of populations at drift-dispersal equilibrium. Using an approach that is based on the slope of the regression of pairwise kinship coefficients on spatial distance and estimates of the effective population density, we compare indirect gene dispersal estimates of sympatric populations of 10 tropical tree species. We re-analysed 26 data sets consisting of mapped allozyme, SSR (simple sequence repeat), RAPD (random amplified polymorphic DNA) or AFLP (amplified fragment length polymorphism) genotypes from two rainforest sites in French Guiana. Gene dispersal estimates were obtained for at least one marker in each species, although the estimation procedure failed under insufficient marker polymorphism, limited sample size, or inappropriate sampling area. Estimates generally suffered low precision and were affected by assumptions regarding the effective population density. Averaging estimates over data sets, the extent of gene dispersal ranged from 150 m to 1200 m according to species. Smaller gene dispersal estimates were obtained in species with heavy diaspores, which are presumably not well dispersed, and in populations with high local adult density. We suggest that limited seed dispersal could indirectly limit effective pollen dispersal by creating higher local tree densities, thereby increasing the positive correlation between pollen and seed dispersal distances. We discuss the potential and limitations of our indirect estimation procedure and suggest guidelines for future studies.