Assessing the expected response to genomic selection of individuals and families in Eucalyptus breeding with an additive-dominant model.

We report a genomic selection (GS) study of growth and wood quality traits in an outbred F2 hybrid Eucalyptus population (n=768) using high-density single-nucleotide polymorphism (SNP) genotyping. Going beyond previous reports in forest trees, models were developed for different selection targets, namely, families, individuals within families and individuals across the entire population using a genomic model including dominance. To provide a more breeder-intelligible assessment of the performance of GS we calculated the expected response as the percentage gain over the population average expected genetic value (EGV) for different proportions of genomically selected individuals, using a rigorous cross-validation (CV) scheme that removed relatedness between training and validation sets. Predictive abilities (PAs) were 0.40-0.57 for individual selection and 0.56-0.75 for family selection. PAs under an additive+dominance model improved predictions by 5 to 14% for growth depending on the selection target, but no improvement was seen for wood traits. The good performance of GS with no relatedness in CV suggested that our average SNP density (~25 kb) captured some short-range linkage disequilibrium. Truncation GS successfully selected individuals with an average EGV significantly higher than the population average. Response to GS on a per year basis was ~100% more efficient than by phenotypic selection and more so with higher selection intensities. These results contribute further experimental data supporting the positive prospects of GS in forest trees. Because generation times are long, traits are complex and costs of DNA genotyping are plummeting, genomic prediction has good perspectives of adoption in tree breeding practice.

Genetic diversity and assessment of 23 microsatellite markers for parentage testing of Santa Inês hair sheep in Brazil.

Santa Inês is the most common hair sheep breed in Brazil and probably has the highest genetic diversity among sheep breeds in this country. Successful breeding programs for Brazilian sheep breeds are not common for various reasons, including a lack of control of parentage in the flocks. We developed an allele frequency database for 23 STR loci for the Santa Inês breed based on 285 animals sampled from five populations distributed across the central-western and north-eastern regions of Brazil. The marker set included seven microsatellites used in the 2011 International Society for Animal Genetics sheep genotyping comparison tests and all eight microsatellites currently approved by the Brazilian Agricultural Ministry laboratory accreditation guidelines for sheep identification. The microsatellites had an average of 10 alleles and a mean expected heterozygosity of 0.745. Combined paternity exclusion probabilities when no parent or one parent was known were >99.99%. A small proportion (5.8%) of the existing genetic variation was found to be among the Santa Inês populations, possibly derived from genetic drift and selection. We found that the marker panel proposed by the Agricultural Ministry, although generally useful, should be enhanced by including more markers for improved exclusionary power in parentage testing. This database provides a useful tool for parentage testing of this major Brazilian breed, contributing to improved management and breeding of existing herds.

Accelerating the domestication of trees using genomic selection: accuracy of prediction models across ages and environments.

• Genomic selection is increasingly considered vital to accelerate genetic improvement. However, it is unknown how accurate genomic selection prediction models remain when used across environments and ages. This knowledge is critical for breeders to apply this strategy in genetic improvement. • Here, we evaluated the utility of genomic selection in a Pinus taeda population of c. 800 individuals clonally replicated and grown on four sites, and genotyped for 4825 single-nucleotide polymorphism (SNP) markers. Prediction models were estimated for diameter and height at multiple ages using genomic random regression best linear unbiased predictor (BLUP). • Accuracies of prediction models ranged from 0.65 to 0.75 for diameter, and 0.63 to 0.74 for height. The selection efficiency per unit time was estimated as 53-112% higher using genomic selection compared with phenotypic selection, assuming a reduction of 50% in the breeding cycle. Accuracies remained high across environments as long as they were used within the same breeding zone. However, models generated at early ages did not perform well to predict phenotypes at age 6 yr. • These results demonstrate the feasibility and remarkable gain that can be achieved by incorporating genomic selection in breeding programs, as long as models are used at the relevant selection age and within the breeding zone in which they were estimated.

An international parentage and identification panel for the domestic cat (Felis catus).

Seventeen commercial and research laboratories participated in two comparison tests under the auspices of the International Society for Animal Genetics to develop an internationally tested, microsatellite-based parentage and identification panel for the domestic cat (Felis catus). Genetic marker selection was based on the polymorphism information content and allele ranges from seven random-bred populations (n = 261) from the USA, Europe and Brazil and eight breeds (n = 200) from the USA. Nineteen microsatellite markers were included in the comparison test and genotyped across the samples. Based on robustness and efficiency, nine autosomal microsatellite markers were ultimately selected as a single multiplex 'core' panel for cat identification and parentage testing. Most markers contained dinucleotide repeats. In addition to the autosomal markers, the panel included two gender-specific markers, amelogenin and zinc-finger XY, which produced genotypes for both the X and Y chromosomes. This international cat parentage and identification panel has a power of exclusion comparable to panels used in other species, ranging from 90.08% to 99.79% across breeds and 99.47% to 99.87% in random-bred cat populations.

Plant microsatellite genotyping with 4-color fluorescent detection using multiple-tailed primers

We extended the concept of fluorescent microsatellite genotyping with a single-universal tailed primer to the simultaneous use of three different tailed primers to allow multiplexed 4-color detection for medium throughput genotyping of plant species. The method was tested on Eucalyptus DNA samples using three forward primer sequences of human microsatellite markers labeled with different fluorescent dyes. The robustness of the method was tested for the simultaneous detection and genetic analysis of microsatellites in a genetic mapping experiment. This method allows reliable and cost-effective genotyping with the same level of multiplexing attained in regular microsatellite fluorescent detection assays. Besides the enhanced quality of the genotypic data provided by the fluorescent detection method when compared to colorimetric ones, the economy brought about by this method becomes greater with an increasing number of microsatellite markers. This method has been particularly useful for genotyping populations of several tropical tree species addressing community-wide population genetics and conservation questions.

A microsatellite-based, gene-rich linkage map for the AA genome of Arachis (Fabaceae).

Cultivated peanut (Arachis hypogaea) is an important crop, widely grown in tropical and subtropical regions of the world. It is highly susceptible to several biotic and abiotic stresses to which wild species are resistant. As a first step towards the introgression of these resistance genes into cultivated peanut, a linkage map based on microsatellite markers was constructed, using an F(2) population obtained from a cross between two diploid wild species with AA genome (A. duranensis and A. stenosperma). A total of 271 new microsatellite markers were developed in the present study from SSR-enriched genomic libraries, expressed sequence tags (ESTs), and by "data-mining" sequences available in GenBank. Of these, 66 were polymorphic for cultivated peanut. The 271 new markers plus another 162 published for peanut were screened against both progenitors and 204 of these (47.1%) were polymorphic, with 170 codominant and 34 dominant markers. The 80 codominant markers segregating 1:2:1 (P<0.05) were initially used to establish the linkage groups. Distorted and dominant markers were subsequently included in the map. The resulting linkage map consists of 11 linkage groups covering 1,230.89 cM of total map distance, with an average distance of 7.24 cM between markers. This is the first microsatellite-based map published for Arachis, and the first map based on sequences that are all currently publicly available. Because most markers used were derived from ESTs and genomic libraries made using methylation-sensitive restriction enzymes, about one-third of the mapped markers are genic. Linkage group ordering is being validated in other mapping populations, with the aim of constructing a transferable reference map for Arachis.

Power of microsatellite markers for fingerprinting and parentage analysis in Eucalyptus grandis breeding populations.

We report the genetic analysis of 192 unrelated individuals of an elite breeding population of Eucalyptus grandis (Hill ex Maiden) with a selected set of six highly polymorphic microsatellite markers developed for species of the genus Eucalyptus. A full characterization of this set of six loci was carried out generating allele frequency distributions that were used to estimate parameters of genetic information content of these loci, including expected heterozygosity, polymorphism information content (PIC), power of exclusion, and probability of identity. The number of detected alleles per locus ranged from 6 to 33, with an average of 19.8 +/- 9.2. The average expected heterozygosity was 0.86 +/- 0.11 and the average PIC was 0.83 +/- 0.16. Using only three loci, it was possible to discriminate all 192 individuals. The overall probability of identity considering all six EMBRA microsatellite markers combined was lower than 1 in 2 billion. An analysis of the sample size necessary to estimate expected heterozygosity with minimum variance indicated that at least 64 individuals have to be genotyped to characterize this parameter with adequate accuracy for most microsatellites in Eucalyptus. The high degree of multiallelism and the clear and simple codominant Mendelian inheritance of the set of microsatellites used provide an extremely powerful system for the unique identification of Eucalyptus individuals for fingerprinting purposes and parentage testing.

Retrospective selection of elite parent trees using paternity testing with microsatellite markers: an alternative short term breeding tactic for Eucalyptus.

The conventional way to drive modifications in old forest tree seed orchards is to establish progeny trials involving each parent tree and then evaluate its contribution to the performance of the progeny by estimating its general and specific combining ability (GCA and SCA). In this work, we successfully applied an alternative parent selection tactic based on paternity testing of superior offspring derived from a hybrid seed orchard established with a single Eucalyptus grandis seed parents and six E. urophylla pollen parents. A battery of 14 microsatellite markers was used to carry out parentage tests of 256 progeny individuals including two independent samples of selected trees and one control unselected sample, all derived from 6-year-old forest stands in eastern Brazil. Paternity determination was carried out for all progeny individuals by a sequential paternity exclusion procedure. Exclusion was declared only when the obligatory paternal allele in the progeny tree was not present in the alleged parent tree for at least four independent markers to avoid false exclusions due to mutation or null alleles. After maternity checks to identify seed mixtures and selfed individuals, the paternity tests revealed that approximately 29% of the offspring was sired by pollen parents outside the orchard. No selfed progeny were found in the selected samples. Three pollen parents were found to have sired essentially all of the offspring in the samples of selected and non-selected progeny individuals. One of these three parents sired significantly more selected progeny than unselected ones ( P< or =0.0002 in a Fisher exact test). Based on these results, low-reproductive-successful parents were culled from the orchard, and management procedures were adopted to minimize external pollen contamination. A significant difference ( P<0.01) in mean annual increment was observed between forest stands produced with seed from the orchard before and after selection of parents and revitalization of the orchard. An average realized gain of 24.3% in volume growth was obtained from the selection of parents as measured in forest stands at age 2-4 years. The marker-assisted tree-breeding tactic presented herein efficiently identified top parents in a seed orchard and resulted in an improved seed variety. It should be applicable for rapidly improving the output quality of seed orchards, especially when an emergency demand for improved seed is faced by the breeder.

Integrating genomics into Eucalyptus breeding

The advent of high throughput genomic technologies has opened new perspectives in the speed, scale and detail with which one can investigate genes, genomes and complex traits in Eucalyptus species. A genomic approach to a more detailed understanding of important metabolic and physiological processes, which affect tree growth and stress resistance, and the identification of genes and their allelic variants, which determine the major chemical and physical features of wood properties, should eventually lead to new opportunities for directed genetic modifications of far-reaching economic impact in forest industry. It should be kept in mind, however, that basic breeding strategies, coupled with sophisticated quantitative methods, breeder's experience and breeder's intuition, will continue to generate significant genetic gains and have a clear measurable impact on production forestry. Even with a much more global view of genetic processes, genomics will only succeed in contributing to the development of improved industrial forests if it is strongly interconnected with intensive fieldwork and creative breeding. Integrated genomic projects involving multi-species expressed sequence tag sequencing and quantitative trait locus detection, single nucleotide polymorphism discovery for association mapping, and the development of a gene-rich physical map for the Eucalyptus genome will quickly move toward linking phenotypes to genes that control the wood formation processes that define industrial-level traits. Exploiting the full power of the superior natural phenotypic variation in wood properties found in Eucalyptus genetic resources will undoubtedly be a key factor to reach this goal.

Genetic structure, mating system, and long-distance gene flow in heart of palm (Euterpe edulis Mart.).

We report a detailed analysis of the population genetic structure, mating system, and gene flow of heart of palm (Euterpe edulis Mart.-Arecaceae) in central Brazil. This palm is considered a keystone species because it supplies fruits for birds and rodents all year and is intensively harvested for culinary purposes. Two populations of this palm tree were examined, using 18 microsatellite loci. The species displays a predominantly outcrossed mating system (tm = 0.94), with a probability of full sibship greater than 70% within open-pollinated families. The following estimates of interpopulation genetic variation were calculated and found significant: FIT = 0.17, FIS = 0.12, FST = 0.06, and RST = 0.07. This low but significant level of interpopulation genetic variation indicates high levels of gene flow. Two adult trees were identified as likely seed parents (P > 99.9%) of juveniles located at a distance of 22 km. Gene flow over such distances has not been reported before for tropical tree species. The establishment and management of in situ genetic reserves or ex situ conservation and breeding populations for E. edulis should contemplate the collection of several hundreds open-pollinated maternal families from relatively few distant populations to maximize the genetic sampling of a larger number of pollen parents.

Resistance to rust ( Puccinia psidii Winter) in eucalyptus: mode of inheritance and mapping of a major gene with RAPD markers.

Rust is one of the most-damaging eucalypt diseases in Brazil and is considered a potential threat to eucalypt plantations worldwide. To determine the mode of inheritance of resistance in the Eucalyptus grandis- Puccinia psidii pathosystem, ten full-sib families, generated from crosses between susceptible and resistant trees, were inoculated with a single-pustule isolate of the pathogen and rust severity was scored. The observed segregation ratios in segregating families suggested major gene control of rust resistance, although clearly incomplete penetrance, variable expressivity and minor genes are also involved in the global rust-resistance response. To identify markers linked to the resistance locus, screening of RAPD polymorphisms was conducted using bulked segregant analysis in a large full-sib family. A linkage group was built around the Ppr1 gene ( P. psidii resistance gene 1) encompassing six RAPD markers, with a genetic window spanning 5 cM with the two most-closely linked flanking markers. Besides these two flanking markers, RAPD marker AT9/917 co-segregated with Ppr1 without a single recombinant in 994 meioses. This tightly linked marker should prove useful for marker-assisted introgression and will provide an initial lead for a positional cloning effort of this resistance allele. This is the first report of a disease resistance gene identified in Eucalyptus, and one of the few examples of the involvement of a major gene in a non-coevolved pathosystem.

Multiplexed systems of microsatellite markers for genetic analysis of mahogany, Swietenia macrophylla King (Meliaceae), a threatened neotropical timber species.

Mahogany (Swietenia macrophylla King [Meliaceae]) is the most valuable hardwood species in the neotropics. Its conservation status has been the subject of increasing concern due to overexploitation and habitat destruction. In this work we report the development and characterization of 10 highly variable microsatellite loci for S. macrophylla. Twenty-nine percent of the 126 sequenced mahogany clones yielded useful microsatellite loci. Three high-throughput genotyping systems were developed based on polymerase chain reaction (PCR) multiplexing of these mahogany loci. We identified a total of 158 alleles in 121 adult individuals of S. macrophylla, with an average of 15.8 alleles (range 11-25) per locus. All loci showed Mendelian inheritance in open-pollinated half-sib families. The mean expected heterozygosity was 0.84 and the mean observed heterozygosity was 0.73. The combined probability of identity-the probability that two individuals selected at random from a population would have identical genotypes--was 7.0 x 10(-15), and combined probability of paternity exclusion was 0.999998 overall loci. These microsatellite loci permit precise estimates of parameters such as gene flow, mating system, and paternity, thus providing important insights into the population genetics and conservation of S. macrophylla.

Towards a genus-wide reference linkage map for Eucalyptus based exclusively on highly informative microsatellite markers.

A novel set of 50 highly polymorphic microsatellite markers were developed and mapped on existing RAPD framework maps of Eucalyptus grandis and E. urophylla. Together with the twenty previously developed microsatellite markers, these were used to align the existing maps for the two most commercially important Eucalyptus species in the tropics. Sixty-three microsatellite markers were placed on the E. grandis map in 11 linkage groups, and 53 on the E. urophylla map distributed in 10 linkage groups. Approximately 66% of the microsatellite markers segregated in a fully informative fashion, allowing the establishment of colinear syntenic linkage groups between the two maps. The 50 new microsatellite markers were highly informative, with an average of 14 alleles per locus, and average expected heterozygosity between 0.82 and 0.87. Furthermore, within the subgenus Symphyomyrtus, to which the vast majority of commercially important Eucalyptus species belong, these markers display on average 90% transportability. This set of 70 mapped microsatellite markers represents a significant step toward the development of a genus-wide reference linkage map for Eucalyptus. These highly multiallelic and transportable markers constitute a powerful tool for QTL discovery and validation, and can be used in directed searches for QTL allele variation across Eucalyptus pedigrees.

Conservation and synteny of SSR loci and QTLs for vegetative propagation in four Eucalyptus species.

Conservation of microsatellite loci, heterozygous in Eucalyptus grandis, Eucalyptus urophylla, Eucalyptus tereticornis and Eucalyptus globulus, allowed us to propose homeologies among genetic linkage groups in these species, supported by at least three SSR loci in two different linkage groups. Marker-trait associations for sprouting and adventitious rooting ability were also compared in the four species. Putative quantitative trait loci (QTLs) influencing vegetative propagation traits were located on homeologous linkage groups. Our findings indicate high transferability of microsatellite markers between Eucalyptus species of the Symphyomyrtus subgenus and establish foundations for the use of synteny in the genetic analysis of this genus. Microsatellite markers should help integrate eucalypt genetic linkage maps from various sources. The availability of comparative linkage maps provides a basis of more-efficient use of genetic information for molecular breeding and evolutionary studies in Eucalyptus.

Cost-effective method to synthesize a fluorescent internal DNA standard for automated fragment sizing.

We describe a simple and cost-effective method for the synthesis of an internal fluorescently labeled DNA standard for fragment sizing using an automated DNA sequencer. A set of primer pairs labeled with ROX was developed to amplify 12 DNA fragments, 58-417 bp, derived from a conserved region of plant chloroplast DNA. These amplified fragments were mixed together, constituting a fluorescent internal DNA size marker. The precision of the size standard was evaluated by estimating the size of 20 alleles that were amplified at four dinucleotide microsatellite loci with the synthesized size standard and the commercial internal sizing standard, GeneScan Rox500. A number of intra-gel and inter-gel comparisons were run, and an analysis of variance was carried out. No significant difference was observed between the size estimates obtained with the synthesized DNA standard and the commercial standard. This facile and general PCR-based method for the synthesis of internal standards allows for significant savings in the implementation of large genotyping experiments using microsatellite or AFLP markers.

Brazilian population database for the 13 STR loci of the AmpFlSTR Profiler Plus and Cofiler multiplex kits.

Allele frequencies for the 13 STR core loci (D3S1358, VWA, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, D7S820, CSF1PO, TPOX, THO1 and D16S539) included in the AmpFlSTR((R)) Profiler Plus and AmpFlSTR((R)) Cofiler kits were obtained for a sample of 700-800 genetically unrelated Brazilians. The expected performance of these loci for personal identification and paternity testing in the Brazilian population was estimated.

High resolution microsatellite based analysis of the mating system allows the detection of significant biparental inbreeding in Caryocar brasiliense, an endangered tropical tree species.

In this work we investigate the mating system of four populations of the endangered tropical tree species Caryocar brasiliense, using genetic data from 10 microsatellite loci. Eight to 10 open-pollinated progeny arrays of 16 individuals, together with their mother tree, were sampled per population. Mating system parameters were estimated under the mixed mating model, implemented by the software MLTR. The single-locus outcrossing rate (ts) varied among loci and populations, but multilocus outcrossing rates (tm) were equal to one for all four populations. Nevertheless, biparental inbreeding (tm - ts) was different from zero for all populations, indicating that outcrossing events may occur between relatives. Our results also indicate that the high polymorphism of microsatellite markers provide an extraordinary resolution to discriminate precisely selfing events from outcrossing events between close relatives. Our results indicate that, although highly outcrossed, C. brasiliense shows high levels of biparental inbreeding, most likely due to the limited flight range of pollinators and restriction in seed dispersal. Furthermore, these results suggest that Cerrado fragmentation could limit gene flow by isolating seed dispersers and territorial small sized bat pollinators inside fragments, increasing the rate of mating between close relatives. The conservation of nonisolated populations in large preserved areas may be necessary to foster outcrossing events between unrelated individuals and thus maintain species viability.

Population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at microsatellite loci.

We report the population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at 10 microsatellite loci. Additionally, we compare heterozygosity and inbreeding estimates for continuous and fragmented populations and discuss the consequences for conservation. For a total of 314 individuals over 10 populations, the number of alleles per locus ranged from 20 to 27 and expected and observed heterozygosity varied from 0.129 to 0.924 and 0.067 to 1.000, respectively. Significant values of theta and R(ST) showed important genetic differentiation among populations. theta was much lower than R(ST), suggesting that identity by state and identity by descent have diverged in these populations. Although a significant amount of inbreeding was found under the identity by descent model (f = 0.11), an estimate of inbreeding for microsatellite markers based on a more adequate stepwise mutation model showed no evidence of nonrandom mating (R(IS) = 0.04). Differentiation (pairwise F(ST)) was positively correlated with geographical distance, as expected under the isolation by distance model. No effect of fragmentation on heterozygosity or inbreeding could be detected. This is most likely due to the fact that Cerrado fragmentation is a relatively recent event (approximately 60 years) compared to the species life cycle. Also, the populations surveyed from both fragmented and disturbed areas were composed mainly of adult individuals, already present prior to ecosystem fragmentation. Adequate hypothesis testing of the effect of habitat fragmentation will require the recurrent analysis of juveniles across generations in both fragmented and nonfragmented areas.

Development and characterization of microsatellite markers for genetic analysis of a Brazilian endangered tree species Caryocar brasiliense.

In this work we report the development and characterization of 10 microsatellite loci for the endangered tree species Caryocar brasiliense. Using genomic library enrichment, the efficiency of SSR marker development was 14.4% from sequencing data to operationally useful loci. Primer sequences for this set of 10 loci are made available together with their estimates of expected heterozygosity, probability of paternity exclusion and probability of identity. Mendelian inheritance and segregation was confirmed for all 10 loci in open-pollinated half-sib families as well as the absolute transferability of these 10 loci to five other species of the same genus. Number of alleles per locus ranged from 10 to 22 with a mean value of 16 and expected heterozygosity varying from 0.84 to 0.94. The combined probability of genetic identity was on the order of 10-17 clearly demonstrating that SSR multilocus genotypes are likely to be unique and capable of readily discriminating individuals of C. brasiliense. The very high combined probability of paternity exclusion (0.99999995) also indicates that these markers will permit detailed parentage studies in natural populations even in situations where both maternity and paternity are unknown. The battery of microsatellite markers developed and characterized in this study opens a new perspective for the generation of fundamental population genetic data for devising sound collection and conservation procedures for C. brasiliense and related species of the genus.

Genetic mapping of quantitative trait loci controlling growth and wood quality traits in Eucalyptus grandis using a maternal half-sib family and RAPD markers.

Quantitative trait loci (QTL) mapping of forest productivity traits was performed using an open pollinated half-sib family of Eucalyptus grandis. For volume growth, a sequential QTL mapping approach was applied using bulk segregant analysis (BSA), selective genotyping (SG) and cosegregation analysis (CSA). Despite the low heritability of this trait and the heterogeneous genetic background employed for mapping, BSA detected one putative QTL and SG two out of the three later found by CSA. The three putative QTL for volume growth were found to control 13.7% of the phenotypic variation, corresponding to an estimated 43.7% of the genetic variation. For wood specific gravity five QTL were identified controlling 24.7% of the phenotypic variation corresponding to 49% of the genetic variation. Overlapping QTL for CBH, WSG and percentage dry weight of bark were observed. A significant case of digenic epistasis was found, involving unlinked QTL for volume. Our results demonstrate the applicability of the within half-sib design for QTL mapping in forest trees and indicate the existence of major genes involved in the expression of economically important traits related to forest productivity in Eucalyptus grandis. These findings have important implications for marker-assisted tree breeding.