Search

Since their initiation in the 1950s, worldwide selective tree breeding programs followed the recurrent selection scheme of repeated cycles of selection, breeding (mating), and testing phases and essentially remained unchanged to accelerate this process or address environmental contingencies and concerns. Here, we introduce an “end-to-end” selective tree breeding framework that: (1) leverages strategically preselected GWAS-based sequence data capturing trait architecture information, (2) generates unprecedented resolution of genealogical relationships among tested individuals, and (3) leads to the elimination of the breeding phase through the utilization of readily available wind-pollinated (OP) families. Individuals’ breeding values generated from multi-trait multi-site analysis were also used in an optimum contribution selection protocol to effectively manage genetic gain/co-ancestry trade-offs and traits’ correlated response to selection. The proof-of-concept study involved a 40-year-old spruce OP testing population growing on three sites in British Columbia, Canada, clearly demonstrating our method’s superiority in capturing most of the available genetic gains in a substantially reduced timeline relative to the traditional approach. The proposed framework is expected to increase the efficiency of existing selective breeding programs, accelerate the start of new programs for ecologically and environmentally important tree species, and address climate-change caused biotic and abiotic stress concerns more effectively.

Crucial to variety improvement programs is the reliable and accurate prediction of genotype’s performance across environments. However, due to the impactful presence of genotype by environment (G×E) interaction that dictates how changes in expression and function of genes influence target traits in different environments, prediction performance of genomic selection (GS) using single-environment models often falls short. Furthermore, despite the successes of genome-wide association studies (GWAS), the genetic insights derived from genome-to-phenome mapping have not yet been incorporated in predictive analytics, making GS models that use Gaussian kernel primarily an estimator of genomic similarity, instead of the underlying genetics characteristics of the populations. Here, we developed a GS framework that, in addition to capturing the overall genomic relationship, can capitalize on the signal of genetic associations of the phenotypic variation as well as the genetic characteristics of the populations. The capacity of predicting the performance of populations across environments was demonstrated by an overall gain in predictability up to 31% for the winter wheat DH population. Compared to Gaussian kernels, we showed that our multi-environment weighted kernels could better leverage the significance of genetic associations and yielded a marked improvement of 4–33% in prediction accuracy for half-sib families. Furthermore, the flexibility incorporated in our Bayesian implementation provides the generalizable capacity required for predicting multiple highly genetic heterogeneous populations across environments, allowing reliable GS for genetic improvement programs that have no access to genetically uniform material.

Modeling environmental spatial heterogeneity can improve the efficiency of forest tree genomic evaluation. Furthermore, genotyping costs can be lowered by reducing the number of markers needed. We investigated the impact on variance components, breeding value accuracy, and bias of two phenotypic data adjustments (experimental design and autoregressive spatial models), and a relationship matrix calculated from a subset of markers selected for their ability to infer ancestry. Using a multiple-trait multiple-site single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) approach, four scenarios (2 phenotype adjustments × 2 marker sets) were applied to diameter at breast height (DBH), height (HT), and resistance to western gall rust (WGR) in four open-pollinated progeny trials of lodgepole pine, with 1490 (out of 11,188) trees genotyped with 25,099 SNPs. As a control, we fitted the conventional ABLUP model using pedigree information. The highest heritability estimates were achieved for the ABLUP followed closely by the ssGBLUP with the full marker set and using the spatial phenotype adjustments. The highest predictive ability was obtained by using a reduced marker subset (8000 SNPs) when either the spatial (DBH: 0.429, and WGR: 0.513) or design (HT: 0.467) phenotype corrections were used. No significant difference was detected in prediction bias among the six fitted models, and all values were close to 1 (0.918–1.014). Results demonstrated that selecting informative markers, such as those capturing ancestry, can improve the predictive ability. The use of spatial correlation structure increased traits’ heritability and reduced prediction bias, while increases in predictive ability were trait-dependent.

Comparing gene expression patterns and genetic polymorphisms between populations is of central importance for understanding the origin and maintenance of biodiversity. Based on population-specific gene expression levels and allele frequency differences, we sought to identify population divergence (PD) genes across the introgression-resistant genomic regions of Populus trichocarpa. Genes containing highly diverged loci [i.e., genetic divergence (GD)] or showing expression divergence (ED) between populations were widely distributed in the genome and substantially enriched in functional categories related to stress responses, disease resistance, timing of flowering, cell cycle regulation, plant growth, and development. Nine genomic regions showing evidence of strong positive selection were overlapped with GD genes, which had significant differences between Oregon (a southernmost peripheral deme) and the other demes. However, we did not find evidence that genes under positive selection show an enrichment for ED. PD genes and genes under selection pertained to the same gene classes, such as SERINE/CYSTEINE PROTEASE, ABC TRANSPORTER, GLYCOSYLTRANSFERASEand other transferases. Our analysis also revealed that GD genes were polymorphic within the species (41.9 ± 3.66 biallelic variants per gene), as previously reported in herbaceous plants. By contrast, ED genes contained less genetic variants (10.73 ± 1.14) and were likely highly expressed. In addition, we found that trans- rather than cis-acting variants considerably contribute to the evolution of >90% PD genes. Overall, this study elucidates that cohorts of PD genes agree with the general attributes of known speciation genes and GD genes will provide substrates for positive selection to operate on.

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism’s potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea(Matt.) Liebl. and Quercus roburL.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species’ evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module “leaf size gradient” was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under “global regulation” in vegetative bud establishment and development. Lamina basal shape and intercalary veins’ number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits’ ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

We conducted an integrative system biology of metabolome and transcriptome profile analyses during pomegranate (Punica granatumL.) seed germination and utilized a weighted gene co-expression network analysis (WGCNA) to describe the functionality and complexity of the physiological and morphogenetic processes as well as gene expression and metabolic differences during seed germination stages. In total, 489 metabolites were detected, including 40 differentially accumulated metabolites. The transcriptomic analysis showed the expression of 6 984 genes changed significantly throughout the whole germination process. Using WGCNA, we identified modules related to the various seed germination stages and hub genes. In the initial imbibition stage (stage 1), the pivotal genes involved in RNA transduction and the glycolytic pathway were most active, while in the sprouting stage (stage 4), the pivotal genes were involved in multiple metabolic pathways. In terms of secondary metabolic pathways, we found flavonoid 4-reductase genes of anthocyanin biosynthesis pathway are most significantly affected during pomegranate seed germination, while the flavonol synthase gene was mainly involved in the regulation of isoflavonoid biosynthesis.

Marker-assisted selection (MAS) is implemented to improve Camellia oleiferayield and fruit attributes for meeting China’s increased demand for new varieties development. We conducted three-generational (G0, G1, and G2) hybridization (diallel mating) and selection experiment and used Sequence-Related Amplified Polymorphism (SRAP). SRAP markers to investigate their utility in a MAS framework. The utilized SRAP markers were instrumental in hybrid authenticity and the identification of matroclinal or patroclinal inheritance presence, thus guiding mating pair selection and direction (their role as male or females). Across the studied 3 generations, estimates of genetic diversity parameters showed steady increase with percentage increase of ((G0to G1and G1to G2) 9.25 and 9.05: observed number of alleles; 3.12 and 7.80: means effective number of alleles; 12.35 and 22.34: Nei‘s gene diversity; and 14.21 and 21.77: Shannon‘s index), indicating lack of diversity reduction associated with selection. Estimates of genetic distance and their correlation with heterosis were useful in guiding selection of mating pairs for achieving the desired yield and fruit attributes (fruit diameter, height, weight, and index, peel thickness, number of seeds per fruit, seed weight per fruit, and seed rate). Most yield and fruit attributes exhibited high broad-sense heritability with increasing trend over generation intervals, indicating the increased potential of hybrid breeding for this species.

Here, we perform cross-generational GS analysis on coastal Douglas-fir (Pseudotsuga menziesii), reflecting trans-generational selective breeding application. A total of 1321 trees, representing 37 full-sib F1families from 3 environments in British Columbia, Canada, were used as the training population for (1) EBVs (estimated breeding values) of juvenile height (HTJ) in the F1generation predicting genomic EBVs of HTJ of 136 individuals in the F2generation, (2) deregressed EBVs of F1HTJ predicting deregressed genomic EBVs of F2HTJ, (3) F1mature height (HT35) predicting HTJ EBVs in F2, and (4) deregressed F1HT35 predicting genomic deregressed HTJ EBVs in F2. Ridge regression best linear unbiased predictor (RR-BLUP), generalized ridge regression (GRR), and Bayes-B GS methods were used and compared to pedigree-based (ABLUP) predictions. GS accuracies for scenarios 1 (0.92, 0.91, and 0.91) and 3 (0.57, 0.56, and 0.58) were similar to their ABLUP counterparts (0.92 and 0.60, respectively) (using RR-BLUP, GRR, and Bayes-B). Results using deregressed values fell dramatically for both scenarios 2 and 4 which approached zero in many cases. Cross-generational GS validation of juvenile height in Douglas-fir produced predictive accuracies almost as high as that of ABLUP. Without capturing LD, GS cannot surpass the prediction of ABLUP. Here we tracked pedigree relatedness between training and validation sets. More markers or improved distribution of markers are required to capture LD in Douglas-fir. This is essential for accurate forward selection among siblings as markers that track pedigree are of little use for forward selection of individuals within controlled pollinated families.

Coastal Douglas-fir (Pseudotsuga menziesii(Mirb.) Franco) is one of the most commercially important softwood species in North America. In British Columbia, Canada, breeding has increased volume gains between 20 and 30%, while 97% of seedlings come from improved seed sources. Branching traits in particular, have a strong influence on strength and stiffness of Douglas-fir wood; however, they are rarely measured. Remotely Piloted Aerial Systems and Airborne Laser Scanning Systems (RPAS-LS) produce high-density three-dimensional point clouds that can be used for the creation of internal geometric features describing individual tree branching structures. We analyzed a Coastal Douglas-fir progeny test trial located in British Columbia, Canada, and developed a new method to estimate branch attributes from RPAS-LS data for inclusion as selection criteria in tree improvement programs. Branch length, angle, width, and volume were estimated for each tree. Narrow-sense heritability (the proportion of variation due to genetics) and genetic correlations were also estimated. The method extracted branch length with a correlation (r) of 0.93 compared to manual measurements. Using these branch attributes, results then show that branch angle had the highest heritability (0.277), while tree height and branch length had the highest genetic correlation (0.668). These findings are encouraging for forest managers as they indicate that branch level metrics should be considered when selecting trees in breeding programs.

The genus Rhododendron(Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Abstract: We present a probabilistic model to minimize the fingerprinting effort associated with the implementation of the “breeding without breeding” scheme under partial pedigree reconstruction. Our approach is directed at achieving a declared target population’s minimum effective population size (N e ), following the pedigree reconstruction and genotypic selection and is based on the graph theory algorithm. The primary advantage of the proposed method is to reduce the cost associated with fingerprinting before the implementation of the pedigree reconstruction for seed parent–offspring derived from breeding arboreta and production or natural populations. Stochastic simulation was conducted to test the method’s efficiency assuming a simple polygenic model and a single trait. Hypothetical population consisted of 30 parental trees that were paired at random (selfing excluded), resulting in 600 individuals (potential candidates for forwards selection). The male parentage was assumed initially unknown. The model was used to estimate the minimum genotyping sample size needed to reaching the prescribed N e . Results were compared with the known pedigree data. The model was successful in revealing the true relationship pattern over the whole range of N e . Two to three offspring entered genotyping to meet the N e  = 2 while 41 to 43 were required to satisfy the N e  = 14. Importantly, genetic gain was affected at the lower limits of the genotyping effort. Doubling the number of parents resulted in considerable reduction of the genotyping effort at higher N e values.

Ginger (Zingiber officinale) is one of the most valued spice plants worldwide; it is prized for its culinary and folk medicinal applications and is therefore of high economic and cultural importance. Here, we present a haplotype-resolved, chromosome-scale assembly for diploid ginger anchored to 11 pseudochromosome pairs with a total length of 3.1 Gb. Remarkable structural variation was identified between haplotypes, and two inversions larger than 15 Mb on chromosome 4 may be associated with ginger infertility. We performed a comprehensive, spatiotemporal, genome-wide analysis of allelic expression patterns, revealing that most alleles are coordinately expressed. The alleles that exhibited the largest differences in expression showed closer proximity to transposable elements, greater coding sequence divergence, more relaxed selection pressure, and more transcription factor binding site differences. We also predicted the transcription factors potentially regulating 6-gingerol biosynthesis. Our allele-aware assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in ginger.

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

AbstractTo confirm patterns of diversity and differentiation found with isozymes and mitochondrial DNA, we surveyed 10 populations of the white pine weevil, Pissodes strobi(Peck), for randomly amplified polymorphic DNA (RAPD) markers. Four weevil populations were sampled from Sitka spruce, Picea sitchensis(Bong.) Carr (Pinaceae), five from the "interior" spruce of British Columbia [admixtures of white spruce, Picea glauca(Moench) Voss, and Engelmann spruce, Picea engelmanni(Parry)], and one from Jack pine, Pinus banksianaLamb. (Pinaceae), in Ontario. In each population, 30–60 weevils were assayed with 10 RAPD primers, yielding 74 RAPD markers. Genetic analyses showed that populations from interior spruce and Jack pine formed a distinct complex; as well, Vancouver Island populations formed a distinct group within the Sitka populations. Levels of diversity, both in terms of polymorphic loci and expected heterozygosity, declined from east to west, supporting the contention that P. strobioriginated in eastern North America and migrated west, and suggesting that biocontrol methods may be more effective on populations from Sitka spruce, owing to their reduced diversity. These results parallel an earlier isozyme study but, in contrast, the diversity differences and population relationships are demonstrated to be statistically significant, owing to both the much larger number of loci sampled and the attachment of statistical confidence intervals to estimates of diversity and differentiation.

The joint levels of outcrossing and contamination in two young northern Swedish Scots pine (Pinus sylvestris L.) seed orchards were studied with the aid of allozyme markers. High levels of contamination (21-36%) were estimated and were attributed to the observed low within orchard pollen production. Outcrossing rate estimates were high (0.957-0.961), though they departed significantly from complete outcrossing (t= 1.0) in spite of the observed high contamination levels. These results confirmed reported relationships between outcrossing and contamination levels in open-pollinated seed orchards. It was concluded that a crop management option relying on supplemental mass-pollination would be effective in reducing both contamination and inbreeding levels.

Microsatellites are codominant and highly polymorphic genetic markers. They are widely used in studying parentage and reproductive success in many insect species. We developed six microsatellite markers in Pissodes storbi. The number of alleles per locus ranges from 5 to 16 alleles, with a mean number of 10. The high allelic diversity observed in this study showed that they could be useful in the study of parentage and reproductive success in P. strobi.

Among conifer species, the resolution of malate dehydrogenase (MDH) band patterns was greatest in Douglas-fir in which there was evidence of four loci, two of which regularly form a heterodimeric product. MDH in lodgepole pine, white spruce, and Sitka spruce also appear to be under control of four loci, while MDH in the true firs, white and grand fir, and western red-cedar appear to be controlled by only three loci.