Phenotypic variation in seedlings of a “keystone” tree species (Quercus douglasii): the interactive effects of acorn source and competitive environment

Blue oak (Quercus douglasii) is a deciduous tree species endemic to California that currently exhibits poor seedling survival to sapling age classes. We used common garden techniques to examine how genetic variation at regional and local scales affected phenotypic expression in traits affecting oak seedling growth and survival. Between-population variation was examined for seedlings grown from acorns collected from a northern, mesic population and a southern, xeric population. Within-population variation was examined by comparing seedlings from different maternal families within the mesic population. Acorns were planted into neighborhoods of an annual dicot (Erodium botrys), an annual grass (Bromus diandrus), and a perennial bunchgrass (Nassella pulchra). By varying the species composition of herbaceous neighborhoods into which acorns were planted, the interactive effects of competition and acorn germplasm source on phenotypic expression could also be examined. Potential maternal effects, expressed as variation in acorn size, were assessed by weighing each acorn before planting. Probability of seedling emergence increased significantly with acorn size in the xeric population but not in the mesic population. Similarly, the effect of acorn size on seedling leaf area, stem weight, and root weight was also population-dependent. At a within-population level, acorn size effects on seedling traits varied significantly among maternal families. In addition to acorn size effects, rates of oak seedling emergence were also dependent on an interaction of population source and competitive environment. Interactions between maternal family and competitive environment in the expression of seedling leaf characters suggest the possibility of genetic variation for plasticity in traits such as specific leaf area. Using carbon isotope discrimination () as an index of relative water-use efficiency (WUE), higher water use efficiency was indicated for oak seedlings grown in the annual plant neighborhoods compared to seedlings grown in the bunchgrass neighborhood. This trend may represent an adaptive plastic response because, compared to the bunchgrass neighborhood, soil water depletion was more rapid within annual plant neighborhoods.     

Estimating heritability using family-pooled phenotypic and genotypic data: a simulation study applied to aquaculture

Estimating heritability based on individual phenotypic and genotypic measurements can be expensive and labour-intensive in commercial aquaculture breeding. Here, the feasibility of estimating heritability using within-family means of phenotypes and allelic frequencies was investigated. Different numbers of full-sib families and family sizes across ten generations with phenotypic and genotypic information on 10 K SNPs were analysed in ten replicates. Three scenarios, representing differing numbers of pools per family (one, two and five) were considered. The results showed that using one pool per family did not reliably estimate the heritability of family means. Using simulation parameters appropriate for aquaculture, at least 200 families of 60 progeny per family divided equally in two pools per family was required to estimate the heritability of family means effectively. Although application of five pools generated more within- and between- family relationships, it reduced the number of individuals per pool and increased within-family residual variation, hence, decreased the heritability of family means. Moreover, increasing the size of pools resulted in increasing the heritability of family means towards one. In addition, heritability of family mean estimates were higher than family heritabilities obtained from Falconer’s formula due to lower intraclass correlation estimate compared to the coefficient of relationship.Subject terms: Genome evolution     

Selection bias in quantitative trait loci mapping

A simulation study was performed to see whether selection affected quantitative trait loci (QTL) mapping. Populations under random selection, under selection among full-sib families, and under selection within a full-sib family were simulated each with heritability of 0.3, 0.5, and 0.7. They were analyzed with the marker spacing of 10 cM and 20 cM. The accuracy for QTL detection decreased for the populations under selection within full-sib family. Estimates of QTL effects and positions differed (P < .05) from their input values. The problems could be ignored when mapping a QTL for the populations under selection among full-sib families. A large heritability helped reduction of such problems. When the animals were selected within a full-sib family, the QTL was detected for the populations with heritability of 0.5 or larger using the marker spacing of 10 cM, and with heritability of 0.7 using the marker spacing of 20 cM. This study implied that when selection was introduced, the accuracy for QTL detection decreased and the estimates of QTL effects were biased. A caution was warranted on the decision of data (including selected animals to be genotyped) for QTL mapping.     

Inferring relatedness and heritability using molecular markers in radiata pine

A set of eight unlinked microsatellite markers was used to estimate relatedness among 355 individuals of a Pinus radiata breeding population. The average performance of open-pollinated progeny of each individual, for wood density, was considered to represent the phenotype of all 355 individuals. Marker-based estimates of relationship were compared with the pedigree-based coefficients of relationships. The phenotypic similarity among all pairs of individuals was regressed on marker-estimated relatedness to estimate the inheritance of wood density. The marker-based estimate of heritability was compared with that obtained using classical quantitative genetic methods. Overall, a low correlation (0.13) was observed between marker-based and pedigree-based estimates of relatedness. After discarding negative estimates of relatedness, the average coefficient of relationship among known groups of maternal half-sibs, full-sibs and unrelated individuals, increased from 0.24 to 0.29 (0.25 expected), from 0.43 to 0.48 (0.50 expected) and from –0.04 to 0.15 (0 expected), respectively. Marker-based and conventional estimates of heritability of wood density were 0.79 and 0.38, respectively. However, by using only marker loci with expected Hardy–Weinberg frequencies, marker-based estimate of heritability was 0.33, which is very similar to that obtained from conventional approaches. The use of molecular markers to understand quantitative genetic variation is discussed.     

A comparison of four experimental designs for the estimation of heritability

Summary The partial diallel cross, the complete diallel cross, and the designs known as North Carolina Experiments 1 and 2 are compared for their usefulness in estimating heritability. It is first shown that reliable values for the sampling mean and variance of heritability estimates are obtained from approximate expressions based on the moments of the chi-square distribution. These expressions are then applied to determine the optimum experimental designs for a range of situations.The main basis for discrimination is the amount of information per unit, defined as i = 1/(N var( ²)), where ² is the estimate of the heritability h ² and N is the number of units in the experiment, either individuals or families.The two parameters considered were the heritability of individuals and the heritability of full-sib families, and for each of these the partial diallel cross was the most preferred, followed in decreasing order of preference by design NC2, the complete diallel, and design NC1.It is first shown that there is no optimum number of parents for a partial diallel cross or male parents for designs NC1 and NC2. The number of crosses per parent for a partial diallel or dams per sire for designs NC2 and NC2 should generally be six or less. Any expansion should be in the direction of using more parents in the case of the partial diallel, or more male parents in the case of designs NC1 and NC2. For the two heritability parameters considered in this study it is inefficient to increase the number of replicates beyond two.     

The beta-binomial model for estimating heritabilities of binary traits

Calculations of individual narrow-sense heritability and family mean heritability of a binary trait in stochastically simulated sib trials in completely randomized block experiments showed that in some situations estimates of realized heritabilities obtained from the mixed linear threshold model could be improved by application of a proposed beta-binomial model. The proposed model adopts the beta-binomial as the conjugate-prior for the distribution of probabilities of observing the binary trait in a genetic entry. Estimation of the beta parameters allows an estimation of selection response and, by linkage to a threshold model for the individual observations, the desired heritabilities can be obtained. The average bias in the betabinomial estimates of heritability and family mean heritability was less than 2%. Improvements over existing procedures were especially manifest at heritabilities above 0.3 and at low overall probabilities of observing the trait (p < 0.30). The lowest root mean square errors were consistently obtained with the algorithm proposed by Harville and Mee (1984). The beta-binomial framework, although restricted to a single random additive genetic effect, further facilitates general analysis, estimation of selection response, and calculation of reliable family mean heritability. Intraclass correlations can be estimated directly from the beta-binomial parameters.     

Characterization of the contribution of shared environmental and genetic factors to metabolic syndrome methylation heritability and familial correlations

Background

Transgenerational epigenetic inheritance has been posited as a possible contributor to the observed heritability of metabolic syndrome (MetS). Yet the extent to which estimates of epigenetic inheritance for DNA methylation sites are inflated by environmental and genetic covariance within families is still unclear. We applied current methods to quantify the environmental and genetic contributors to the observed heritability and familial correlations of four previously associated MetS methylation sites at three genes (CPT1A, SOCS3 and ABCG1) using real data made available through the GAW20.

Results

Our findings support the role of both shared environment and genetic variation in explaining the heritability of MetS and the four MetS cytosine-phosphate-guanine (CpG) sites, although the resulting heritability estimates were indistinguishable from one another. Familial correlations by type of relative pair generally followed our expectation based on relatedness, but in the case of sister and parent pairs we observed nonsignificant trends toward greater correlation than expected, as would be consistent with the role of shared environmental factors in the inflation of our estimated correlations.

Conclusions

Our work provides an interesting and flexible statistical framework for testing models of epigenetic inheritance in the context of human family studies. Future work should endeavor to replicate our findings and advance these methods to more robustly describe epigenetic inheritance patterns in human populations.

     

Growth and cold hardiness of intervarietal hybrids of douglas-fir

Summary Potentials for improvement of the interior variety of Douglas-fir (Pseudotsuga menziesii var. glauca) by hybridization with the coastal variety (P.m. var. menziesii) were explored. The primary objective was to assess possibilities for increasing growth potentials of the interior variety while maintaining adaptation to relatively cold inland environments. Seventy full-sib hybrid families and their half-sib parental lines were grown in two contrasting environments common to the interior variety. Nine traits related to growth, phenology, and freezing tolerance in 4-year-old trees were compared.For traits related to growth (height and diameter) in the inland environment, hybrids equaled the growth of the coastal variety and exceeded the interior variety by 40 percent. For traits related to adaptation (bud burst, bud set, frost damage, tree form, and freezing tolerance), hybrids were intermediate but approached levels characteristic of the interior variety. Survival of hybrids equaled that of the interior variety and was superior to that of the coastal variety.Hybrid characters could not be predicted reliably from those of parental lines. Yet, quantitative genetic analyses suggest that expression of characters related to growth depends on nonadditive genetic effects, but expression of those related to adaptation is somewhat dependent on additive effects.Realization of the tremendous potential of hybridization for improvement of the interior variety will require at least one backcross generation or additional crosses utilizing introgressed populations.     

Segregation analysis of α-L-fucosidase activity

Summary Complex segregation analysis of plasma -L-fucosidase in 45 British families provides evidence for an additive major gene causing low activities of fucosidase. There was no significant evidence of polygenic heritability or common family environment.     

Multienvironment genomic variance decomposition analysis of open-pollinated Interior spruce (<Emphasis Type="Italic">Picea glauca</Emphasis> x <Emphasis Type="Italic">engelmannii</Emphasis>)

The advantages of open-pollinated (OP) family testing over controlled crossing (i.e., structured pedigree) are the potential to screen and rank a large number of parents and offspring with minimal cost and efforts; however, the method produces inflated genetic parameters as the actual sibling relatedness within OP families rarely meets the half-sib relatedness assumption. Here, we demonstrate the unsurpassed utility of OP testing after shifting the analytical mode from pedigree- (ABLUP) to genomic-based (GBLUP) relationship using phenotypic tree height (HT) and wood density (WD) and genotypic (30k SNPs) data for 1126 38-year-old Interior spruce (Picea glauca (Moench) Voss x P. engelmannii Parry ex Engelm.) trees, representing 25 OP families, growing on three sites in Interior British Columbia, Canada. The use of the genomic realized relationship permitted genetic variance decomposition to additive, dominance, and epistatic genetic variances, and their interactions with the environment, producing more accurate narrow-sense heritability and breeding value estimates as compared to the pedigree-based counterpart. The impact of retaining (random folding) vs. removing (family folding) genetic similarity between the training and validation populations on the predictive accuracy of genomic selection was illustrated and highlighted the former caveats and latter advantages. Moreover, GBLUP models allowed breeding value prediction for individuals from families that were not included in the developed models, which was not possible with the ABLUP. Response to selection differences between the ABLUP and GBLUP models indicated the presence of systematic genetic gain overestimation of 35 and 63% for HT and WD, respectively, mainly caused by the inflated estimates of additive genetic variance and individuals’ breeding values given by the ABLUP models. Extending the OP genomic-based models from single to multisite made the analysis applicable to existing OP testing programs.     

Within-population genetic diversity of Taraxacum officinale (Asteraceae): differential genotype response and effect on interspecific competition

Clonal offspring of five morphologically distinct individuals of Taraxacum officinale were planted in a greenhouse experiment with each of three competitors, Plantago major, Poa pratensis and Trifolium pratense. The competitors were chosen to represent a series of competitive environments experienced by a natural population of T. officinale through the year. Differences in size, morphology, and response to the competitive environments were found among clones and support classification of the five individuals as distinct genotypes. Both differential competitive responses (alteration in performance) and competitive effects (impediment by competitor performance) were exhibited among genotypes. The differential response by the T. officinale genotypes to the competitors indicates that the biotic environment may influence the genetic structure of a population. The biotic environment in this case is determined by the sequential appearance and dominance of competitors in a field rather than the spatial distribution of these competitors. Since competitors change in relative dominance across seasons, competition is likely to be a component of the genotype by season interaction that had been observed in the natural population. Thus, differential responses among genotypes to a temporally and spatially fluctuating biotic environment may contribute to the maintenance of within-population genetic diversity.     

Heritability and shared environment estimates for myopia and associated ocular biometric traits: the Genes in Myopia (GEM) family study

To examine the familial correlations, heritability (h ²) and common environmental components (c ²) of myopia and ocular biometric traits (all treated as continuous outcomes) in families collected through the Genes in Myopia (GEM) family study in Australia. A total of 132 pedigrees (723 participants) were recruited for this study. All individuals completed a risk factor questionnaire and underwent a detailed eye examination including spherical equivalent (SphE) and ocular biometric measurements of axial length (AL), anterior chamber depth (ACD) and corneal curvature (CC). Familial correlations were calculated and h ² and c² were estimated using a variance component model that assumes a multivariate t distribution within each pedigree. Two definitions of common environments (c ²) were considered: nuclear family (current) shared environment (Model 1) and sib-ship (childhood) shared environment (Model 2). Population ascertainment adjustment was performed using the Blue Mountains eye study dataset. The trends observed for familial correlations suggested that SphE is influenced by both environmental and genetic factors whereas AL, ACD and CC are predominantly genetically determined. This was largely confirmed by variance components modelling. Heritability estimates (adjusted for age, sex and years of education) from the best fitting ACE model (Model 2, childhood shared environment) were 0.50 ± 0.05 for SphE, 0.73 ± 0.04 for AL, 0.78 ± 0.04 for ACD and 0.16 ± 0.06 for CC. Childhood environmental effects were significant with c ² estimated to be 0.33 ± 0.04 for SphE, 0.06 ± 0.03 for AL, 0.22 ± 0.04 for ACD and 0.10 ± 0.05 for CC. Age was associated with SphE, total years of education was associated with AL and sex was associated with all traits studied. We used a novel and conservative approach to account for and estimate common environmental effects by specifying either nuclear family or sib-ship environment when estimating heritability estimates and showed that all traits examined (SphE, AL, ACD and CC) are heritable, thus reflecting a genetic component. These traits therefore all represent candidates for quantitative trait linkage analyses.     

Effects of propagule type on genetic parameters of wood density and growth in a loblolly pine progeny test at ages 10 and 11?years

Nine full-sib families of loblolly pine (Pinus taeda L.) were produced by a 3 × 3 factorial mating design. Rooted cuttings and seedlings of full-sib families were tested together in two field locations. Twelve-millimeter wood increment cores were collected from 10- and 11-year-old test trees. On each of the two sites, there were six blocks and a split-plot design, with propagule type as the whole plot and family as the sub-plot. In addition to the collection of wood samples, height and diameter of 1,600 trees were measured. No significant differences were found between cuttings and seedlings for wood density and growth traits. Significant family variation was found for growth and wood density. Genetic parameters estimated for wood density and growth traits using seedlings and rooted cuttings showed that individual-tree and family heritability estimates from rooted cuttings were similar to or higher than those from seedlings for all traits. Half-sib breeding values for parents were highly correlated based on seedling and rooted cutting estimates for height (0.95) and wood density (0.99) but not for diameter (0.56), which suggests that wood density and height breeding value estimates from rooted cuttings in clonal progeny tests can be estimated by traditional seedling tests, but not for tree diameter.     

Iterative kriging for removing spatial autocorrelation in analysis of forest genetic trials

Conventional analysis of spatially correlated data in inadequately blocked field genetic trials may give erroneous results that would seriously affect breeding decisions. Forest genetic trials are commonly very large and strongly heterogeneous, so adjustments for micro-environmental heterogeneity become indispensable. This study explores the use of geostatistics to account for the spatial autocorrelation in four Pinus pinaster Ait. progeny trials established on hilly and irregular terrains with a randomized complete block design and large blocks. Data of five different traits assessed at age 8 were adjusted using an iterative method based on semivariograms and kriging, and the effects on estimates of variance components, heritability, and family effects were evaluated in relation to conventional analysis. Almost all studied traits showed nonrandom spatial structures. Therefore, after the adjustments for spatial autocorrelation, the block and family × block variance components, which were extremely high in the conventional analysis, almost disappeared. The reduction of the interaction variance was recovered by the family variance component, resulting in higher heritability estimates. The removal of the spatial autocorrelation also affected the estimation of family effects, resulting in important changes in family ranks after the spatial adjustments. Comparison among families was also greatly improved due to higher accuracy of the family effect estimations. The analysis improvement was larger for growth traits, which showed the strongest spatial heterogeneity, but was also evident for other traits such as straightness or number of whorls. The present paper demonstrates how spatial autocorrelation can drastically affect the analysis of forest genetic trials with large blocks. The iterative kriging procedure presented in this paper is a promising tool to account for this spatial heterogeneity.     

Variation in the relative magnitude of intraspecific and interspecific competitive effects in novel versus familiar environments in two<Emphasis Type="Italic">Drosophila</Emphasis> species

Models of competitor coevolution, especially the genetic feedback hypothesis, suggest that a negative correlation between intraspecific and interspecific competitive effects may be important in sustaining competitor coexistence, and can give rise to oscillatory dynamics with repeated reversals of competitive superiority. I reanalyzed previously published census data from an experiment in which populationsof Drosophila melanogaster andD. simulans underwent competitive coevolution in one familiar and two novel environments, to specifically look for any evidence of a negative relationship between intraspecific and interspecific competitive effects on population growth rates, and for any indication of short period cycling in the relative magnitude of intraspecific and interspecific competitive effects. While there was considerable variation in the relative magnitude of intraspecific and interspecific competitive effects over generations, among both populations and environments, there was no clear evidence supporting the genetic feedback hypothesis. Intraspecific and interspecific competitive effects on population growth rates were strongly positively correlated in novel environments, and uncorrelated in the familiar environment. Data from the familiar environment indicated that indices of competition of populations of the initially superior competitor,D. melanogaster, might be showing some cyclic behaviour, but I argue that this is likely to be transient, and not suggestive of sustained oscillatory dynamics predicted by the genetic feedback model. I discuss the results in the context of the importance of the genetic architecture of intraspecific and interspecific competitive abilities in determining the coevolutionary trajectory of competitive interactions.     

Spacing and Competition Between Planted Grass Plugs and Preexisting Perennial Grasses in a Restoration Site in Oregon

Planting native species into restoration settings where other natives already occur is a common practice. However, the competitive consequences of such plantings are rarely studied. Planting density also affects restoration costs. Here we examined the effects of established individuals of Lemmon's needlegrass (Achnatherum lemmonii) on plugs of bluebunch wheatgrass (Pseudoroegneria spicata) and Idaho fescue (Festuca idahoensis) in a restoration site in Oregon. All three of these grasses are local native perennials. Plugs were planted at 6, 12, and 18 cm from established A. lemmonii bunchgrasses and also in plots without A. lemmonii neighbors. Plug survival was uniformly high, averaging more than 98%. Plugs planted at 6 cm from established grasses showed significantly lower growth and reproduction than plugs planted at 18 cm, which had similar values to plugs not planted in the vicinity of A. lemmonii. These results suggest that interplanting distances of as little as 18 cm were sufficient to greatly reduce competitive effects on newly planted plugs, at least in early establishment at this site.     

Contribution of photosynthetic rate to growth and reproduction in Amaranthus hybridus

While it is known that genetic variation for photosynthetic and growth traits exists in natural populations, the functional significance of this variation remains unclear, particularly for photosynthetic traits. To test the hypothesis that photosynthetic rate has direct effects on reproduction as well as contributing indirectly to reproduction through effects on growth, we compared wild-type Amaranthus hybridus families to those with a single gene mutation that confers a lower photosynthetic rate. Wild-type and photosynthetic-mutant families were grown in competitive and non-competitive environments and we compared size, biomass allocation, architecture, and reproduction at three developmental stages. To assess the contributions of individual growth traits to reproduction, we calculated covariances between standardized traits and relative fitness (selection differentials), and compared selection between the two biotypes. Finally, we used path analysis to calculate the indirect effects of photosynthetic rate on fitness through growth. The size, allocation, and architecture of photosynthetic mutants did not differ from those of the wild type in either the competitive or non-competitive environment, with the exception that they were taller by the last developmental stage. However, the reproductive biomass of the photosynthetic mutants was significantly reduced compared to the wild type. In the competitive environment, the wild type achieved greater fitness because, while similar in size to the mutants, at any given size it produced more reproductive biomass. This suggests that photosynthetic rate affected the linkage between plant size and reproduction and is evidence of an indirect contribution to fitness. In the non-competitive environment, there were fewer differences in selection differentials between the two plant genotypes, suggesting fewer indirect effects. Path analysis showed that variation in photosynthetic biotype had indirect effects on reproductive biomass, via growth traits, and that there were no direct effects. Photosynthetic rate appears to have fitness consequences primarily through multiple contributions to growth throughout development. Received: 27 March 1998 / Accepted: 28 August 1998     

Rapid evolution of a heteroplasmic repetitive sequence in the mitochondrial DNA control region of carnivores

We describe a repetitive DNA region at the 3 end of the mitochondrial DNA (mtDNA) control region and compare it in 21 carnivore species representing eight carnivore families. The sequence and organization of the repetitive motifs can differ extensively between arrays; however, all motifs appear to be derived from the core motif ACGT. Sequence data and Southern blot analysis demonstrate extensive heteroplasmy. The general form of the array is similar between heteroplasmic variants within an individual and between individuals within a species (varying primarily in the length of the array, though two clones from the northern elephant seal are exceptional). Within certain families, notably ursids, the array structure is also similar between species. Similarity between species was not apparent in other carnivore families, such as the mustelids, suggesting rapid changes in the organization and sequence of some arrays. The pattern of change seen within and between species suggests that a dominant mechanism involved in the evolution of these arrays is DNA slippage. A comparative analysis shows that the motifs that are being reiterated or deleted vary within and between arrays, suggesting a varying rate of DNA turnover. We discuss the evolutionary implications of the observed patterns of variation and extreme levels of heteroplasmy.By acceptance of this article, the publisher acknowledges the right of the US Government to retain non-exclusive, royalty-free license in and to any copyright covering the article. Correspondence to: A.R. Hoetzel     

Effects of nutrient amendments on fine root biomass in a primary successional forest in Hawai'i

Summary Thalli of Lobaria pulmonaria (L.) Hoffm., a nitrogen-fixing epiphyte common in mesic temperate forests, were collected in a Douglas-fir (Pseudotsuga menziesii Franco) forest near Corvallis, Oregon, and maintained for 20 to 40 days in controlled-environment chambers with atmospheric CO₂ concentrations of 374 and 700 ll^-1. Nitrogenase activity, which was assayed by the acetylene reduction method, was approximately doubled in the lichen maintained in elevated CO₂. Increases in nitrogen fixation by lichens may be an important part of the integrated ecosystem response to rising CO₂.     

The mating system in natural and shelterwood stands of Douglas-fir

Summary Mating systems in two pairs of old-growth uncut and adjacent shelterwood stands of Douglas-fir (Pseudotsuga menziesii var. menziesii) were compared by estimating the proportions of viable progenies due to outcrossing (t) with both single-locus and multilocus techniques. Single-locus population estimates (_s), ranging from 0.41 to 1.16, were significantly (P<0.05) heterogeneous among loci in three of four stands; mean single-locus estimates for shelterwoods were not significantly different from those for uncut stands. Multilocus population estimates (_m) ranged from 0.94 to 1.00; again, estimates for shelterwoods were not significantly different from those for uncut stands. Multilocus estimates were slightly higher than mean single-locus estimates for uncut stands but were nearly equivalent for shelterwoods, suggesting that related matings other than selfs may be associated with uncut stands, but not shelterwoods. Individual-tree outcrossing rates (_{m i}), estimated for six trees in each shelterwood, ranged from 0.90 to 1.10 and were significantly heterogeneous among trees at one of the two shelterwoods. Outcrossing was high (> 0.90) in both uncult and shelterwood stands, and no evidence indicated that low parent-tree density had affected stand outcrossing rates.FRL 1918, Forest Research Laboratory. Oregon State University, Corvallis, USA