Monitoring genetic diversity in tropical trees with multilocus dominant markers

Since no universal codominant markers are currently available, dominant genetic markers, such as amplified fragment length polymorphism (AFLP), are valuable tools for assessing genetic diversity in tropical trees. However, the measurement of genetic diversity (H) with dominant markers depends on the frequency of null homozygotes (Q) and the fixation index (F) of populations. While Q can be estimated for AFLP loci, F is less accessible. Through a modelling approach, we show that the monolocus estimation of genetic diversity is strongly dependent on the value of F, but that the multilocus diversity estimate is surprisingly robust to variations in F. The robustness of the estimate is due to a mechanistic effect of compensation between negative and positive biases of H by different AFLP loci exhibiting contrasting frequency profiles of Q. The robustness was tested across contrasting theoretical frequency profiles of Q and verified for 10 neotropical species. Practical recommendations for the implementation of this analytical method are given for genetic surveys in tropical trees, where such markers are widely applied.     

Natural hybridization within seed sources of shortleaf pine (Pinus echinata Mill.) and loblolly pine (Pinus taeda L.)

Shortleaf and loblolly pine trees (n = 93 and 102, respectively) from 22 seed sources of the Southwide Southern Pine Seed Source Study plantings or equivalent origin were evaluated for amplified fragment length polymorphism (AFLP) variation. These sampled trees represent shortleaf pine and loblolly pine, as they existed across their native geographic ranges before intensive forest management. Using 17 primer pairs, a total of 96 AFLPs between shortleaf pine and loblolly pine were produced and scored on the sample trees and two control-pollinated F1 interspecies hybrids and their parents. In addition, the well known isocitrate dehydrogenase (IDH) isozyme marker was scored for all trees. IDH detected two putative hybrids among the loblolly pine samples and two among the shortleaf pine samples, while either 13 or 12 putative hybrids were detected using all AFLP markers and IDH and either NewHybrids or Structure software, respectively. Results of this study show that later generation hybrids can be reliably identified using AFLP markers and confirmed that IDH is not a definitive marker for detecting hybrids; that is, at least in some seed sources, the alternative species’ IDH allele resides in the source species. Based on all the markers, hybridization frequency varied geographically, ranging from 30% in an Arkansas seed source to 0% in several other seed sources. The hybridization level was higher in populations west of the Mississippi River than in populations east of the river; the shortleaf pine hybridization rates were 16.3% and 2.4% and the loblolly pine rates were 4.5% and 3.3%, west and east of the river, respectively. The results suggest that hybridization between these two species is significant but varies by seed source and species, and the potential for the unintended creation of hybrids should be considered in forest management decisions regarding both natural and artificial regeneration.     

Phylogeographic insights into an irruptive pest outbreak

Irruptive forest insect pests cause considerable ecological and economic damage, and their outbreaks have been increasing in frequency and severity. We use a phylogeographic approach to understand the location and progression of an outbreak by the MPB (Dendroctonus ponderosae Hopkins), an irruptive bark beetle that has caused unprecedented damage to lodgepole pine forests in western North America and is poised to expand its range across the boreal forest. We sampled MPB populations across British Columbia and Alberta and used phylogeographic methods to describe lineage diversification, characterize population structure, investigate expansion dynamics, and identify source populations of the outbreak. Using 1181 bp of mitochondrial DNA sequence from 267 individuals, we found high haplotype diversity, low nucleotide diversity, and limited lineage diversification. The overall pattern was consistent with isolation by distance at a continental scale, and with reduced diversity and population structure in the northerly, outbreak regions. Post-Pleistocene expansion was detected, however more recent expansion signals were not detected, potentially due to the size and rapid rate of range expansion. Based on the limited genetic structure, there were likely multiple source populations in southern British Columbia, although the magnitude of the demographic expansion and rate of spread have obscured the signature of these source populations. Our data highlight the need for caution in interpreting phylogeographic results for species with similar demographics.     

Wildfire provides refuge from local extinction but is an unlikely driver of outbreaks by mountain pine beetle

Bark beetle outbreaks and wildfire are important disturbances in conifer ecosystems, yet their interactions are not well understood. We evaluated whether fire injury increased susceptibility of lodgepole pines (Pinus contorta) to mountain pine beetle (Dendroctonus ponderosae Hopkins), how it influenced beetle reproductive success, and whether beetle population phase altered this interaction. Eight sites that experienced wildfire and eight unburned sites were examined in the Greater Yellowstone Ecosystem (USA). Half were in areas where D. ponderosae was undergoing outbreaks, and half were in areas with low populations. We examined 2056 trees one year after fire for burn injury and beetle attack. We quantified beetle reproductive success in a random sample of 106 trees, and measured gallery areas of D. ponderosae and competing subcortical herbivores in 79 additional trees. Baited flight traps sampled stand-level populations of subcortical herbivores and predators. Wildfire predisposed trees to D. ponderosae attack, but nonlinearly, with moderately injured trees being most preferred. This tree-level interaction was influenced by stand-level beetle population size, in that both healthy and fire-injured trees of all classes were attacked where populations were high, but no healthy trees, and only low and moderately injured trees were killed where populations were low. The number of adult brood produced per female was likewise curvilinear, being highest in moderately injured trees. This reflected an apparent trade-off, with high intraspecific competition arising from the large number of beetles needed to overcome defenses in healthy trees, vs. high interspecific competition and low substrate quality in more injured trees. These results suggest that fire-injured trees can provide a reservoir for D. ponderosae during periods when populations are too low to overcome defenses of healthy trees, and might otherwise face localized extinction. However, the likelihood of populations increasing from endemic to outbreak levels in response to increased susceptibility is offset by the opposing constraints of lower substrate quality and higher competitor load in severely injured hosts, and the relative scarcity of moderately injured trees. Wildfire may confer some reproductive increases to populations already outbreaking. We present a conceptual model of how these disturbances and inherent feedbacks interact to affect beetle population dynamics.     

Invasion of diverse habitats by few Japanese knotweed genotypes is correlated with epigenetic differentiation

The expansion of invasive species challenges our understanding of the process of adaptation. Given that the invasion process often entails population bottlenecks, it is surprising that many invasives appear to thrive even with low levels of sequence-based genetic variation. Using Amplified Fragment Length Polymorphism (AFLP) and methylation sensitive-AFLP (MS-AFLP) markers, we tested the hypothesis that differentiation of invasive Japanese knotweed in response to new habitats is more correlated with epigenetic variation than DNA sequence variation. We found that the relatively little genetic variation present was differentiated among species, with less differentiation among sites within species. In contrast, we found a great deal of epigenetic differentiation among sites within each species and evidence that some epigenetic loci may respond to local microhabitat conditions. Our findings indicate that epigenetic effects could contribute to phenotypic variation in genetically depauperate invasive populations. Deciphering whether differences in methylation patterns are the cause or effect of habitat differentiation will require manipulative studies.     

Microsatellite versus AFLP analyses of pre-management introgression levels in loblolly pine (<Emphasis Type="Italic">Pinus taeda</Emphasis> L.) and shortleaf pine (<Emphasis Type="Italic">P. echinata</Emphasis> Mill.)

Loblolly pine and shortleaf pine are known to naturally hybridize. In this study, we used 42 microsatellite markers and isocitrate dehydrogenase isozyme to create genetic profiles of 202 loblolly and shortleaf pine trees grown from seed collected in the 1950s for the Southwide Southern Pine Seed Source Study. Estimated Φ_PT was low in both loblolly (0.061) and shortleaf (0.080) pines, indicating that most of the diversity is accounted for within seed sources. However, both loblolly and shortleaf pines showed significant correlations between seed sources’ genetic and geographic distances, with R ² of 0.43 and 0.17, respectively. The hybridization rate was 4.0%, with more hybrids west of the Mississippi River (8.1%) than east of the river (2.1%). Additionally, about the same proportion of both species (4.5% of loblolly and 3.3% of shortleaf pine) were identified as hybrids. These results are consistent with prior studies on these two species but do contrast with the results from an amplified fragment length polymorphism (AFLP) analysis of the same samples. For example, the AFLP study concluded that 6.3% of the trees were hybrids, or 1.4 times higher than determined by this study. Of the 12 hybrids identified in the AFLP study, six were not identified as hybrids here, and of the eight hybrids identified here, only four were identified in the AFLP study. Although similar in overall results, we suggest that the microsatellite analysis is more convincing than the AFLP analysis because microsatellites provide more information per genetic locus than do AFLPs.     

Mountain pine beetle attack associated with low levels of 4-allylanisole in ponderosa pine

Mountain pine beetle (Dendroctonus ponderosae) is the most important insect pest in southern Rocky Mountain ponderosa pine (Pinus ponderosa) forests. Tree mortality is hastened by the various fungal pathogens that are symbiotic with the beetles. The phenylpropanoid 4-allylanisole is an antifungal and semiochemical for some pine beetle species. We analyzed 4-allylanisole and monoterpene profiles in the xylem oleoresin from a total of 107 trees at six sites from two chemotypes of ponderosa pine found in Colorado and New Mexico using gas chromatography-mass spectroscopy (GC-MS). Although monoterpene profiles were essentially the same in attacked and nonattacked trees, significantly lower levels of 4-allylanisole were found in attacked trees compared with trees that showed no evidence of attack for both chemotypes.     

Genetic diversity within and among populations of shortleaf pine (Pinus echinata Mill.) and loblolly pine (Pinus taeda L.)

Shortleaf pine (n = 93) and loblolly pine (n = 112) trees representing 22 seed sources or 16 physiographic populations were sampled from Southwide Southern Pine Seed Source Study plantings located in Oklahoma, Arkansas, and Mississippi. The sampled trees were grown from shortleaf pine and loblolly pine seeds formed in 1951 and 1952, prior to the start of intensive forest management across their native ranges. Amplification fragment length polymorphism (AFLP) markers were developed and used to study genetic diversity and its structure in these pine species. After screening 48 primer pairs, 17 and 21 pairs were selected that produced 794 and 647 AFLPs in shortleaf pine and loblolly pine, respectively. High-AFLP-based genetic diversity exists within shortleaf pine and loblolly pine, and most (84.73% in shortleaf pine; 87.69% in loblolly pine) of this diversity is maintained within physiographic populations. The high value of unbiased measures of genetic identity and low value of genetic distance for all pairwise comparisons indicates that the populations have similar genetic structures. For shortleaf pine, there was no significant correlation between geographic distance and genetic distance (r = 0.28), while for loblolly pine there was a weak but significant correlation (r = 0.51).     

Phenotypic plasticity and genetic isolation-by-distance in the freshwater mussel <Emphasis Type="Italic">Unio pictorum</Emphasis> (Mollusca: Unionoida)

Freshwater mussels (Unionoida) show high intraspecific morphological variability, and some shell morphological traits are believed to be associated with habitat conditions. It is not known whether and which of these ecophenotypic differences reflect underlying genetic differentiation or are the result of phenotypic plasticity. Using 103 amplified fragment length polymorphism (AFLP) markers, we studied population genetics of three paired Unio pictorum populations sampled from two different habitat types (marina and river) along the River Thames. We found genetic differences along the Thames which were consistent with a pattern of isolation by distance and probably reflect limited dispersal via host fish species upon which unionoid larvae are obligate parasites. No consistent genetic differences were found between the two different habitat types suggesting that morphological differences in the degree of shell elongation and the shape of dorso-posterior margin are caused by phenotypic plasticity. Our study provides the first good evidence for phenotypic plasticity of shell shape in a European unionoid and illustrates the need to include genetic data in order properly to interpret geographic patterns of morphological variation.     

ALLOZYME AND MORPHOLOGICAL DIFFERENTIATION OF MOUNTAIN PINE BEETLES DENDROCTONUS PONDEROSAE HOPKINS (COLEOPTERA: SCOLYTIDAE) ASSOCIATED WITH HOST TREE

The purpose of this study was to determine whether mountain pine beetles utilizing different host species were differentiated for either morphological or protein variation. Genetic differentiation among host species has been reported for the southern pine beetle, the Douglas-fir beetle, the jeffrey pine beetle, and the mountain pine beetle. However, in these studies, the host trees were sampled at separate sites, and hence geographic variation and variation due to host tree were confounded. The mountain pine beetle occasionally utilizes three host trees (ponderosa pine, lodgepole pine, and limber pine) at single sites in Colorado. Five polymorphic enzyme loci and six morphological characters were used to describe beetles resident in different hosts. Differentiation within a site among host trees was detected at two of five polymorphic proteins, and for both size and morphological shape. The magnitude of genetic differentiation among hosts within a site was approximately equivalent to the magnitude of differentiation among sites. These data suggest that the species of host tree may be an important biotic factor associated with the genetic structure of bark beetle communities. The results are discussed in terms of their potential role in the process of speciation by host race formation.     

The impact of phloem nutrients on overwintering mountain pine beetles and their fungal symbionts

In the low nutrient environment of conifer bark, subcortical beetles often carry symbiotic fungi that concentrate nutrients in host tissues. Although bark beetles are known to benefit from these symbioses, whether this is because they survive better in nutrient-rich phloem is unknown. After manipulating phloem nutrition by fertilizing lodgepole pine trees (Pinus contorta Douglas var. latifolia), we found bolts from fertilized trees to contain more living individuals, and especially more pupae and teneral adults than bolts from unfertilized trees at our southern site. At our northern site, we found that a larger proportion of mountain pine beetle (Dendroctonus ponderosae Hopkins) larvae built pupal chambers in bolts from fertilized trees than in bolts from unfertilized trees. The symbiotic fungi of the mountain pine beetle also responded to fertilization. Two mutualistic fungi of bark beetles, Grosmannia clavigera (Rob.-Jeffr. & R. W. Davidson) Zipfel, Z. W. de Beer, & M. J. Wingf. and Leptographium longiclavatum Lee, S., J. J. Kim, & C. Breuil, doubled the nitrogen concentrations near the point of infection in the phloem of fertilized trees. These fungi were less capable of concentrating nitrogen in unfertilized trees. Thus, the fungal symbionts of mountain pine beetle enhance phloem nutrition and likely mediate the beneficial effects of fertilization on the survival and development of mountain pine beetle larvae.     

Comparative population genetic study of two oligophagous insects associated with the same hosts

A parallel study of the genetic structure of two oligophagous species associated with the same hosts was conducted to determine the main factors shaping the distribution of genetic diversity. The bark beetle Tomicus piniperda and the pine processionary moth (PPM) Thaumetopoea pityocampa are both associated with the genus Pinus and belong to different guilds (xylophagous vs defoliating species). The PPM is an ectophagous species that feeds on the needles of living trees, whereas T. piniperda is endophagous and bores galleries in the inner bark of weakened trees. Both species were sampled in the main regions of France, and their genetic structure was assessed after genotyping with five microsatellite markers. Populations of the PPM were significantly structured. A pattern of isolation by distance was found when distances were calculated as bypassing the Massif Central, whereas no such pattern could be found with raw geographic distances. On the contrary, most populations of T. piniperda were not differentiated. No effect of host species could be detected in either of the two species. We conclude that the two taxa have contrasting effective dispersal rates per generation, and we hypothesize that this reflects the different selection pressures acting on individual fitness via different strategies of host use.     

Contrasting patterns of genetic diversity at three different genetic markers in a marine mammal metapopulation

Many studies use genetic markers to explore population structure and variability within species. However, only a minority use more than one type of marker and, despite increasing evidence of a link between heterozygosity and individual fitness, few ask whether diversity correlates with population trajectory. To address these issues, we analysed data from the Steller's sea lion, Eumetiopias jubatus , where three stocks are distributed over a vast geographical range and where both genetic samples and detailed demographic data have been collected from many diverse breeding colonies. To previously published mitochondrial DNA (mtDNA) and microsatellite data sets, we have added new data for amplified fragment length polymorphism (AFLP) markers, comprising 238 loci scored in 285 sea lions sampled from 23 natal rookeries. Genotypic diversity was low relative to most vertebrates, with only 37 loci (15.5%) being polymorphic. Moreover, contrasting geographical patterns of genetic diversity were found at the three markers, with Nei's gene diversity tending to be higher for AFLPs and microsatellites in rookeries of the western and Asian stocks, while the highest mtDNA values were found in the eastern stock. Overall, and despite strongly contrasting demographic histories, after applying phylogenetic correction we found little correlation between genetic diversity and either colony size or demography. In contrast, we were able to show a highly significant positive relationship between AFLP diversity and current population size across a range of pinniped species, even though equivalent analyses did not reveal significant trends for either microsatellites or mtDNA.     

Evaluation of AFLP for genetic mapping in Pinus radiata D. Don

Efficient construction of reasonable density genetic linkage maps is an essential component of QTL detection programmes. The AFLP technique has been used to produce genetic linkage maps in a range of species. We have developed protocols to generate reproducible AFLP profiles in Pinus radiata and have evaluated the inheritance and informativeness of AFLP markers in this important timber species. The large genome size of P. radiata necessitated increased levels of selection at both the pre-amplification and selective amplification steps of the AFLP protocol to generate reproducible AFLP profiles. Once optimised ca. 41.3 scorable AFLP bands were resolvable through denaturing gels, of which 48.4% were polymorphic in a screen of eight unrelated trees. This level of polymorphism is ca. three times higher than with RAPD markers. The total number of bands and the number of polymorphismic bands per PCR were ca. halved when AFLPs were electrophoresed on non-denaturing gels and stained with ethidium bromide. Using the protocols developed, AFLP is an efficient method for generating the DNA markers required for genetic linkage map construction in P. radiata. This revised version was published online in June 2006 with corrections to the Cover Date.     

Assessment of the validity of the sections in Musa (musaceae) using AFLP

Musa L. (Musaceae) is currently separated into five sections (Musa. Rhodochlamys, Callimusa, Australimusa and Ingentimusa) based on chromosome numbers and morphological characters. However, the validation of this classification system is questioned due to the common occurrence of hybridizations across sections and the system not accommodating anomalous species. This study employed amplified fragment length polymorphism (AFLP) in a phenetic examination of the relationships among four sections (material of sect. Ingentimusa was not available) to evaluate whether their genetic differences justify distinction into separate groups. Using eight primer combinations, a total of 276 bands was scored, of which 275 were polymorphic. Among the monomorphic bands, 11 unique markers were identified that revealed the distinct separation of the 11-chromosome species from the 10-chromosome species. AFLP results suggest that species of sect. Rhodochlamys should be combined into a single section with species of sect. Musa, and likewise for species of sect. Australimnusa to be merged with those of sect. Callimusa.     

Testing for host‐associated differentiation in two egg parasitoids of a forest herbivore

Several studies underline the importance of ecological barriers and differential selection in driving sympatric speciation. Host‐associated differentiation (HAD) has been proposed as one of the mechanisms leading to sympatric speciation. However, it is still unclear how common HAD is or which are the factors that could promote it. In particular, not much is known about HAD in predators and parasitoids of herbivorous insects. One of the characteristics postulated to pre‐dispose insects to HAD is parthenogenesis as it may favour adaptive responses to particular environments, amplifying selected gene complexes. In this study, we used amplified fragment length polymorphism (AFLP) markers to determine whether HAD is present in two parthenogenetic egg parasitoids attacking the same herbivore species – the pine processionary moth, Thaumetopoea pityocampa (Denis & Schiffermüller) (Lepidoptera: Notodontidae) – on two host Pinus species. A total of 100 loci for 59 individuals sampled in four populations of Baryscapus servadeii (Domenichini) (Hymenoptera: Eulophidae), a specialist parasitoid, and 106 loci for 117 individuals sampled in six populations of Ooencyrtus pityocampae Mercet (Hymenoptera: Encyrtidae), a generalist parasitoid, were analysed. Levels of genetic differentiation were also assessed with an outlier analysis, checking for alleles associated to host plants. No evidence of HAD was detected in any of the two parasitoid species. We hypothesize that both the lack of strict parthenogenetic reproduction and the ectophagous nature of the insect host could explain the absence of HAD. The genetic variation observed in the generalist parasitoid responded to a pattern of local adaptation, whereas no relationship with either host or geography was found in the specialist parasitoid.     

Genetic structure in the seabuckthorn carpenter moth (Holcocerus hippophaecolus) in China: the role of outbreak events, geographical and host factors

Understanding factors responsible for structuring genetic diversity is of fundamental importance in evolutionary biology. The seabuckthorn carpenter moth (Holcocerus hippophaecolus Hua) is a native species throughout the north of China and is considered the main threat to seabuckthorn, Hippophae rhamnoides L. We assessed the influence of outbreaks, environmental factors and host species in shaping the genetic variation and structure of H. hippophaecolus by using Amplified Fragment Length Polymorphism (AFLP) markers. We rejected the hypothesis that outbreak-associated genetic divergence exist, as evidenced by genetic clusters containing a combination of populations from historical outbreak areas, as well as non-outbreak areas. Although a small number of markers (4 of 933 loci) were identified as candidates under selection in response to population densities. H. hippophaecolus also did not follow an isolation-by-distance pattern. We rejected the hypothesis that outbreak and drought events were driving the genetic structure of H. hippophaecolus. Rather, the genetic structure appears to be influenced by various confounding bio-geographical factors. There were detectable genetic differences between H. hippophaecolus occupying different host trees from within the same geographic location. Host-associated genetic divergence should be confirmed by further investigation.     

Local genetic structure in a North American epiphytic lichen, Ramalina menziesii (Ramalinaceae)

Epiphytic lichens possess unique life history traits that can have conflicting effects on genetic structure: symbiotic mutualism between a fungus with its algal or cyanobacterial photobiont, association with a host plant, and ability to reproduce sexually and asexually. Our study species, Ramalina menziesii, has small ascospores that can facilitate long-distance gene movement, and it is capable of clonal reproduction. The goals of this study are to test whether different haplotypes were differentially distributed across host plant species, to look for evidence of asexual vs. sexual reproduction, and to assess the local genetic structure of the population. We sampled individuals from multiple trees of three oak species in four lichen subpopulations within a savanna ecosystem. Using DNA sequence data from four fungal nuclear loci, we found no tendency for host specialization. Alleles were randomly distributed across subpopulations. The frequency of multilocus genotypes was consistent with a randomly mating population. Sexual reproduction involving relichenization appeared to be the predominant mode of reproduction of R. menziesii at this study site. We found no significant local genetic structure suggesting widespread gene flow at the local scale. The genetic structure of this lichen is comparable to that of widely distributed epiphytic plants.     

Modeling cold tolerance in the mountain pine beetle, Dendroctonus ponderosae

Cold-induced mortality is a key factor driving mountain pine beetle, Dendroctonus ponderosae, population dynamics. In this species, the supercooling point (SCP) is representative of mortality induced by acute cold exposure. Mountain pine beetle SCP and associated cold-induced mortality fluctuate throughout a generation, with the highest SCPs prior to and following winter. Using observed SCPs of field-collected D. ponderosae larvae throughout the developmental season and associated phloem temperatures, we developed a mechanistic model that describes the SCP distribution of a population as a function of daily changes in the temperature-dependent processes leading to gain and loss of cold tolerance. It is based on the changing proportion of individuals in three states: (1) a non cold-hardened, feeding state, (2) an intermediate state in which insects have ceased feeding, voided their gut content and eliminated as many ice-nucleating agents as possible from the body, and (3) a fully cold-hardened state where insects have accumulated a maximum concentration of cryoprotectants (e.g. glycerol). Shifts in the proportion of individuals in each state occur in response to the driving variables influencing the opposite rates of gain and loss of cold hardening. The level of cold-induced mortality predicted by the model and its relation to extreme winter temperature is in good agreement with a range of field and laboratory observations. Our model predicts that cold tolerance of D. ponderosae varies within a season, among seasons, and among geographic locations depending on local climate. This variability is an emergent property of the model, and has important implications for understanding the insect's response to seasonal fluctuations in temperature, as well as population response to climate change. Because cold-induced mortality is but one of several major influences of climate on D. ponderosae population dynamics, we suggest that this model be integrated with others simulating the insect's biology.