Barrier to gene flow between two ecologically divergent Populus species, P. alba (white poplar) and P. tremula (European aspen): the role of ecology and life history in gene introgression

The renewed interest in the use of hybrid zones for studying speciation calls for the identification and study of hybrid zones across a wide range of organisms, especially in long-lived taxa for which it is often difficult to generate interpopulation variation through controlled crosses. Here, we report on the extent and direction of introgression between two members of the "model tree" genus Populus: Populus alba (white poplar) and Populus tremula (European aspen), across a large zone of sympatry located in the Danube valley. We genotyped 93 hybrid morphotypes and samples from four parental reference populations from within and outside the zone of sympatry for a genome-wide set of 20 nuclear microsatellites and eight plastid DNA restriction site polymorphisms. Our results indicate that introgression occurs preferentially from P. tremula to P. alba via P. tremula pollen. This unidirectional pattern is facilitated by high levels of pollen vs. seed dispersal in P. tremula (pollen/seed flow = 23.9) and by great ecological opportunity in the lowland floodplain forest in proximity to P. alba seed parents, which maintains gene flow in the direction of P. alba despite smaller effective population sizes (N(e)) in this species (P. alba N(e)c. 500-550; P. tremula N(e)c. 550-700). Our results indicate that hybrid zones will be valuable tools for studying the genetic architecture of the barrier to gene flow between these two ecologically divergent Populus species.     

Clonality and spatial genetic structure in Populus x canescens and its sympatric backcross parent P. alba in a Central European hybrid zone.

Spatial genetic structure (SGS) holds the key to understanding the role of clonality in hybrid persistence, but multilocus SGS in hybrid zones has rarely been quantified. Here, the aim was to fill this gap for natural hybrids between two diploid, ecologically divergent European tree species with mixed sexual/asexual reproduction, Populus alba and P. tremula. Nuclear microsatellites were used to quantify clonality, SGS, and historical gene dispersal distances in up to 407 trees from an extensive Central European hybrid zone including three subpopulation replicates. The focus was on P. x canescens and its backcross parent P. alba, as these two genotypic classes co-occur and interact directly. Sexual recombination in both taxa was more prominent than previously thought, but P. x canescens hybrids tended to build larger clones extending over larger areas than P. alba. The 3.4 times stronger SGS in the P. x canescens genet population was best explained by a combination of interspecific gene flow, assortative mating, and increased clonality in hybrids. Clonality potentially contributes to the maintenance of hybrid zones of P. alba and P. tremula in time and space. Both clonality and SGS need to be taken into account explicitly when designing population genomics studies of locus-specific effects in hybrid zones.     

Genetic analysis of post-mating reproductive barriers in hybridizing European Populus species

Molecular genetic analyses of experimental crosses provide important information on the strength and nature of post-mating barriers to gene exchange between divergent populations, which are topics of great interest to evolutionary geneticists and breeders. Although not a trivial task in long-lived organisms such as trees, experimental interspecific recombinants can sometimes be created through controlled crosses involving natural F(1)'s. Here, we used this approach to understand the genetics of post-mating isolation and barriers to introgression in Populus alba and Populus tremula, two ecologically divergent, hybridizing forest trees. We studied 86 interspecific backcross (BC(1)) progeny and >350 individuals from natural populations of these species for up to 98 nuclear genetic markers, including microsatellites, indels and single nucleotide polymorphisms, and inferred the origin of the cytoplasm of the cross with plastid DNA. Genetic analysis of the BC(1) revealed extensive segregation distortions on six chromosomes, and >90% of these (12 out of 13) favored P. tremula donor alleles in the heterospecific genomic background. Since selection was documented during early diploid stages of the progeny, this surprising result was attributed to epistasis, cyto-nuclear coadaptation, heterozygote advantage at nuclear loci experiencing introgression or a combination of these. Our results indicate that gene flow across 'porous' species barriers affects these poplars and aspens beyond neutral, Mendelian expectations and suggests the mechanisms responsible. Contrary to expectations, the Populus sex determination region is not protected from introgression. Understanding the population dynamics of the Populus sex determination region will require tests based on natural interspecific hybrid zones.     

A set of novel DNA polymorphisms within candidate genes potentially involved in ecological divergence between Populus alba and P. tremula, two hybridizing European forest trees

We have identified 53 DNA (single nucleotide, microsatellite, and insertion-deletion) polymorphisms within 12 candidate genes potentially involved in ecological differences between Populus alba and P. tremula, two hybridizing European forest trees. The genes represent candidates for functional roles associated with abiotic or biotic stress response, cross-talk, phenology, and leaf development. Distributions within sequences, intraspecific levels of diversity, and genetic divergence (F(ST) ) between species are reported for each polymorphism, as are haplotype frequencies for each gene. The markers will be used for population genomic studies of the barrier to gene flow between these two ecologically divergent forest trees.     

Evaluation of qPCR reference genes in two genotypes of Populus for use in photoperiod and low-temperature studies

ABSTRACT: BACKGROUND: Quantitative PCR (qPCR) is a widely used technique for gene expression analysis. A common normalization method for accurate qPCR data analysis involves stable reference genes to determine relative gene expression. Despite extensive research in the forest tree species Populus, there is not a resource for reference genes that meet the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) standards for qPCR techniques and analysis. Since Populus is a woody perennial species, studies of seasonal changes in gene expression are important towards advancing knowledge of this important developmental and physiological trait. The objective of this study was to evaluate reference gene expression stability in various tissues and growth conditions in two important Populus genotypes (P. trichocarpa "Nisqually 1" and P. tremula X P. alba 717 1-B4) following MIQE guidelines. RESULTS: We evaluated gene expression stability in shoot tips, young leaves, mature leaves and bark tissues from P. trichocarpa and P. tremula. x P. alba grown under long-day (LD), short-day (SD) or SD plus low-temperatures conditions. Gene expression data were analyzed for stable reference genes among 18S rRNA, ACT2, CDC2, CYC063, TIP4-like, UBQ7, PT1 and ANT using two software packages, geNormPLUS and BestKeeper. GeNormPLUS ranked TIP4-like and PT1 among the most stable genes in most genotype/tissue combinations while BestKeeper ranked CDC2 and ACT2 among the most stable genes. CONCLUSIONS: This is the first comprehensive evaluation of reference genes in two important Populus genotypes and the only study in Populus that meets MIQE standards. Both analysis programs identified stable reference genes in both genotypes and all tissues grown under different photoperiods. This set of reference genes was found to be suitable for either genotype considered here and may potentially be suitable for other Populus species and genotypes. These results provide a valuable resource for the Populus research community.     

Matrix attachment regions (MARs) enhance transformation frequency and transgene expression in poplar

We tested the value of a matrix attachment region (MAR) fragment derived from a tobacco gene for increasing the frequency of Agrobacterium-mediated transformation. A binary vector that carried a GUS reporter gene containing an intron and an nptII gene was modified to contain flanking MAR elements within the T-DNA borders. Vectors containing or lacking MARs were then used to transform tobacco, a readily transformabl e poplar clone (Populus tremula × P. alba), and a recalcitrant poplar clone (Populus trichocarpa × P. deltoides). MARs increased GUS gene expression approximately 10-fold in the two hybrid poplar clones and twofold in tobacco one month after cocultivation with Agrobacterium; MARs also increased the frequency of kanamycin-resistant poplar shoots recovered     

Comparative Nucleotide Diversity Across North American and European Populus Species

Nucleotide polymorphisms in two North American balsam poplars (Populus trichocarpa Torr. & Gray and P. balsamifera L.; section Tacamahaca), and one Eurasian aspen (P. tremula L.; section Populus) were compared using nine loci involved in defense, stress response, photoperiodism, freezing tolerance, and housekeeping. Nucleotide diversity varied among species and was highest for P. tremula (θ(w) = 0.005, π(T) = 0.007) as compared to P. balsamifera (θ(w) = 0.004, π(T) = 0.005) or P. trichocarpa (θ(w) = 0.002, π(T) = 0.003). Across species, the defense and the stress response loci accounted for the majority of the observed level of nucleotide diversity. In general, the studied loci did not deviate from neutral expectation either at the individual locus (non-significant normalized Fay and Wu's H) or at the multi-locus level (non-significant HKA test). Using molecular clock analysis, section Tacamahaca probably shared a common ancestor with section Populus approximately 4.5 million year ago. Divergence between the two closely related balsam poplars was about 0.8 million years ago, a pattern consistent with an isolation-with-migration (IM) model. As expected, P. tremula showed a five-fold higher substitution rate (2 × 10(-8) substitution/site/year) compared to the North American species (0.4 × 10(-8) substitution/site/year), probably reflecting its complex demographic history. Linkage disequilibrium (LD) varied among species with a more rapid decay in the North American species (<400 bp) in comparison to P. tremula (?400 bp). The similarities in nucleotide diversity pattern and LD decay of the two balsam poplar species likely reflects the recent time of their divergence.     

Comparative genome mapping among Populus adenopoda, P. alba, P. deltoides, P. euramericana and P. trichocarpa

Among the genus Populus, the sections Populus (white poplar), Aigeiros Duby (black poplar) and Tacamahaca Spach contain many tree species of economical and ecological important properties. Two parental maps for the inter-specific hybrid population of Populus adenopoda × P. alba (two species of Populus section) were constructed based on SSR and SRAP markers by means of a two-way pseudo-test cross mapping strategy. The same set of SSR markers developed from the P. trichocarpa (belonging to Tacamahaca section) genome which were used to construct the maps of P. deltoides and P. euramericana (two species of Aigeiros section) was chosen to analyze the genotype of the experimental population of P. adenopoda × P. alba. Using the mapped SSR markers as allelic bridges, the alignment of the white and black poplar maps to each other and to the P. trichocarpa physical map was conducted. The alignment showed high degree of marker synteny and colinearity and the closer relationship between Aigeiros and Tacamahaca sections than that of Populus and Tacamahaca. Moreover, there was evidence for the chromosomal duplication and inter-chromosomal reorganization involving some poplar linkage groups, suggesting a complicated course of fission or fusion in one of the lineages. A poplar consensus map based on the comparisons could be constructed will be useful in practical applications including marker assisted selection.     

An excess of nonsynonymous polymorphism and extensive haplotype structure at the PtABI1B locus in European aspen (Populus tremula): a case of balancing selection in an obligately outcrossing plant?

Here, we describe an unusually pronounced haplotype structure at the PtABI1B locus in the obligately outcrossing tree Populus tremula. Both nucleotide diversity and divergence at PtABI1B was low compared to other P. tremula genes suggesting that the gene is located in a region with a low mutation rate. Despite this, PtABI1B shows a very marked excess of nonsynonymous polymorphisms across the entire coding region and linkage disequilibrium (LD) extending across the entire PtABI1B region of approximately 2.6 kb. Such extensive LD is normally not seen in P. tremula. The extensive LD at PtABI1B is caused by the presence of two distinct haplotypes. The haplotype structure is not caused by a lack of recombination in the region, because evidence of recombination can be detected. In addition, several statistical tests strongly reject neutrality for the PtABI1B region, suggesting that the unusual haplotype structure could be actively maintained by balancing selection.     

rDNA cloning and rapid hybrid identification inPopulus spp. (Salicaceae)

Ribosomal DNA genes fromP. deltoides have been cloned and specific sequences of the 25 S and 18 S rDNA region, labelled by digoxigenin, have been used to determine the rDNA structure ofPopulus tremula, P. fremontii, P. maximowiczii, P. yunnanensis, P. nigra, P. wislizenii, P. alba. The restriction maps of the coding region appeared to be similar among the examined species and with those ofP. deltoides andP. trichocarpa, reported in a previous paper. Inter- and intraspecific variation in rDNA repeat unit length have been revealed after EcoRI digestions. SstI and XbaI restriction sites have been found at different positions in the IGS of some species. The polymorphic fragments generated by SstI digestion allowed the identification of the hybrid origin of some genotypes. The number of rDNA genes in the genome ofP. deltoides has been estimated to be about 2 000 copies. Finally, the usefulness of these studies inPopulus spp. taxonomy and forestry genetics is discussed.Ribosomal RNA gene structure in somePopulus spp. (Salicaceae) and their hybrids 2.     

Regulation of glutathione synthesis in leaves of transgenic poplar (Populus tremula X P. alba) overexpressing glutathione synthetase

The poplar hybrid Populus tremula X P. alba was transformed with the Escherichia coli gene for glutathione synthetase ( gsh II ) targetted to the cytosol. Leaves of five lines of transgenic plants exhibited glutathione synthetase activities 15- to 60-fold higher than those of wild-type plants. Total glutathione levels and GSH/GSSG ratios were similar in transgenic and wild-type plants. Precursor feeding experiments with cysteine and γ-glutamylcysteine suggest that glutathione synthesis in the cytoplasm is controlled by a multistep procedure that includes (i) the availability of cysteine, (ii) the availability of γ-glutamylcysteine, and (iii) regulation of the activities of both γ-glutamylcysteine synthetase and glutathione synthetase. However step (ii) may set an upper limit for the cellular glutathione content.     

Salicylate and catechol levels are maintained in nahG transgenic poplar

Metabolic profiling was used to investigate the molecular phenotypes of a transgenic Populus tremula x P. alba hybrid expressing the nahG transgene, a bacterial gene encoding salicylate hydroxylase that converts salicylic acid to catechol. Despite the efficacy of this transgenic approach to reduce salicylic acid levels in other model systems and thereby elucidate roles for salicylic acid in plant signaling, transgenic poplars had similar foliar levels of salicylic acid and catechol to that of non-transformed controls and exhibited no morphological phenotypes. To gain a deeper understanding of the basis for these observations, we analyzed metabolic profiles of leaves as influenced by transgene expression. Expression of nahG decreased quinic acid conjugates and increased catechol glucoside, while exerting little effect on levels of salicylic acid and catechol, the substrate and product, respectively, of the nahG enzyme. This suggests a biological role of elevated constitutive salicylic acid levels in Populus, in contrast to other plant systems in which nahG dramatically reduces salicylic acid levels.     

Isolation and characterization of microsatellites in trembling aspen (Populus tremuloides)

 We have identified, isolated, and characterized microsatellite/simple sequence repeat (SSR) loci in trembling aspen (Populus tremuloides) by screening partial genomic libraries. We have also examined the compatibility and use of the P. tremuloides SSR primers to resolve microsatellites in other Populus species. Fourteen microsatellites were identified from 1600 clones screened. The TC/AG microsatellites were the most abundant. A total of 29 alleles were detected in 36 P. tremuloides individuals at the four SSR loci (two each of di- and tri-nucleotide repeats) characterized. The number of alleles at the SSR loci ranged from 5 to 11, with an average of 7.25 alleles per locus, and the observed heterozygosity ranged from 0.19 to 0.82, with a mean of 0.46 per locus. Although the highest polymorphism was observed for a dinucleotide SSR locus, the trinucleotide SSR loci showed substantial polymorphism. There were 34 unique multilocus genotypes among the 36 P. tremuloides individuals examined, and 89% of the individuals had unique multilocus genotypes. Two pairs of SSR primers were successful in PCR, amplifying genomic DNA and resolving microsatellites of comparable size from Populus deltoides, P. nigra, P.×canadensis, and P. maximowiczii. The microsatellite DNA markers developed could be used for clonal fingerprinting, certification of controlled crosses, genome mapping, marker-assisted early selection, genetic diversity assessments, and conservation and sustainable management of poplar genetic resources. Received: 14 November 1997 / Accepted: 17 November 1997     

RT-PCR analysis and stress response capacity of transgenic gshI-poplar clones (Populus x canescens) in response to paraquat exposure

Stress response capacity (Fv/Fm at 690 nm and F690/F735 at Fmax) of untransformed hybrid poplar, Populus x canescens (P tremula x P alba), and two transgenic lines overexpressing gamma-ECS (gamma-glutamylcysteine synthetase) either in the cytosol (cyt-ECS) or in the chloroplast (chl-ECS) was studied in response to the herbicide paraquat (4.0 x 10(-9) to 4.0 x 10(-6) M) for 21 days. Significant differences at sublethal (4.0 x 10(-7) M) and bleaching (4.0 x 10(-6) M) concentrations of paraquat were observed with about a two-fold and eight-fold decrease in the photosynthetic activity (Fv/Fm at 690 nm and F690/F735 at Fmax), respectively. None of the gshI transgenic lines (cyt-ECS, chl-ECS) with elevated GSH content exhibited significant tolerance to paraquat. Semiquantitative RT-PCR of the cyt-ECS clone was used for gene expression analysis of the nuclear encoded rbcS gene and the stress responsive gst gene. Expression of the constitutively expressed 26SrRNA ribosomal gene was probed as a control for all RT-PCR reactions. The relative intensities of gene expressions normalized to the level of 26SrRNA intensity showed a 50% decrease in the nuclear encoded rbcS expression and a 120% increase in the stress responsive gst gene expression of the paraquat treated (4.0 x 10(-7) M) samples of the transgenic poplar line (cyt-ECS).     

Study on Origin of Populus tomentosa Carr.

This study analyzed five varieties of P. tomentosa and their putative parents, P. alba, P. davidiana, P. adenopoda and P. tremula by RAPD and artificial cross. Of 78 bands revealed by 13 arbitrary primers selected, 58 (74.36%) bands were polymorphic among P. tomentosa and four putative parents. The RAPD statistical results showed the relationship between P. tomentosa and P. alba or P. adenopoda was closer than that between P. tomentosa and P. davidiana or P. tremula . There were higher similarity and lower genetic distance between P. tomentosa and P. alba or P. adenopoda than that between P. tomentosa and the other two putative parents. The results from UPGMA, Fitch margoliash phylogenetic dendrogram were similar. In addition, artificial hybrids including reciprocal cross between P. alba and P. adenopoda were used to compare with five varieties of P. tomentosa . And those hybrids between P. alba as female and P. adenopoda as male were closer to P. tomentosa than those from reciprocal cross. The results demonstrated that P. tomentosa was a natural hybrid between P. alba as female parent and P. adenopoda as male parent.     

用RAPD探讨毛白杨起源

利用RAPD分子标记对山杨、欧洲山杨、银白杨、响叶杨及毛白杨的5个不同类型,即抱头毛白杨、截叶毛白杨、易县毛白杨、陕西毛白杨和南京栽培的毛白杨进行了分析,并进行了银白杨×响叶杨正反交实验,研究杂种后代与毛白杨的关系。RAPD分析说明,毛白杨的亲缘关系与银白杨和响叶杨较近,而与山杨和欧洲山杨的亲缘关系则较远。毛白杨在聚类图上首先是与银白杨和响叶杨聚为一类,再与山杨和欧洲山杨聚在一起。通过比较毛白杨与银白杨和响叶杨的正反交子代,发现毛白杨更接近银白杨×响叶杨的子代,而与反交子代有所区别。由此可见,毛白杨为天然杂种是没有疑义的,而且很可能是银白杨×响叶杨的杂种。本研究筛选出的13个引物,产生了78条谱带,覆盖了基因组的一定区域,其可靠性是较高的。     

Genetic structure of hybrid zones between Pinus pumila and P. parviflora var. pentaphylla (Pinaceae) revealed by molecular hybrid index analysis

Pinus species have three differently inherited genomes: paternal chloroplast, maternal mitochondrial, and biparental nuclear. Our previous study on the hybrid zones between alpine Pinus pumila and montane to subalpine P. parviflora var. pentaphylla demonstrated contrasting patterns of introgression of two cytoplasmic genomes, i.e., the paternal cpDNA flowed from P. parviflora var. pentaphylla to P. pumila, and the maternal mtDNA flowed in the reverse direction. In the present study, we developed codominant nuclear DNA markers diagnostic or mostly diagnostic for each parental species by single-strand conformation polymorphism (SSCP) of polymerase chain reaction (PCR) products, using expressed sequence tag (EST) primers of Pinus taeda. To describe the introgressive patterns of the nuclear genes, the molecular hybrid index (MHI) showing the overall proportion of alleles inferred to be derived from P. pumila was determined for each plant collected in hybrid zones on Mt. Asahidake and Mt. Higashiazuma, Japan. At Mt. Asahidake, the MHI values changed clinally according to the altitudes at which the plants were collected. However, at Mt. Higashiazuma, there was a gap in the MHI values between the plants above and below the Abies and Tsuga forest zone (alt. 1800-1900 m). This suggested that the zone plays a role in creating an effective barrier to gene flow in the hybrid zone.     

First isolation of an isoprene synthase gene from poplar and successful expression of the gene in Escherichia coli

For the first time, the complete functional gene for isoprene synthase has been isolated from poplar (Populus alba x Populus tremula). The gene was quite similar to known limonene and other monoterpene synthases, but was found to specifically catalyze the formation of isoprene from the precursor dimethylallyl diphosphate with only a marginal activity for the formation of the monoterpene limonene from geranyl diphosphate as compared with limonene synthases. Omitting the part of the gene that putatively encoded the signal peptide necessary for transport into the chloroplast led to an enhanced rate of isoprene formation by the recombinant protein.