Transcriptomic Prediction of Breeding Values in Loblolly Pine

Phenotypic variation in forest trees can be partitioned into subsets controlled by genetic variation and by environmental factors, and heritability expressed as the proportion of total phenotypic variation attributed to genetic variation. Applied tree breeding programs can use matrices of relationships, based either on recorded pedigrees in structured breeding populations or on genotypes of molecular genetic markers, to model genetic covariation among related individuals and predict genetic values for individuals for whom no phenotypic measurements are available. This study tests the hypothesis that genetic covariation among individuals of similar genetic value will be reflected in shared patterns of gene expression. We collected gene expression data by high-throughput sequencing of RNA isolated from pooled seedlings from parents of known genetic value, and compared alternative approaches to data analysis to test this hypothesis. Selection of specific sets of transcripts increased the predictive power of models over that observed using all transcripts. Using information on presence of putative mutations in protein-coding sequences increased predictive accuracy for some traits but not for others. Known pedigree relationships are not required for this approach to modeling genetic variation, so it has potential to allow broader application of genetic covariance modeling to natural populations of forest trees.