Population divergence in the chloroplast genome of Eucalyptus nitens

The complete chloroplast genome sequencing analysis revealed an unusual IRs reduction in three species of subfamily Zygophylloideae

Tetraena mongolica, Zygophyllum xanthoxylon, and Z. fabago are three typical dryland plants with important ecological values in subfamily Zygophylloideae of Zygophyllaceae. Studies on the chloroplast genomes of them are favorable for understanding the diversity and phylogeny of Zygophyllaceae. Here, we sequenced and assembled the whole chloroplast genomes of T. mongolica, Z. xanthoxylon, and Z. fabago, and performed comparative genomic and phylogenetic analysis. The total size, structure, gene content and orders of these three chloroplast genomes were similar, and the three chloroplast genomes exhibited a typical quadripartite structure with a large single-copy region (LSC; 79,696–80,291 bp), a small single-copy region (SSC; 16,462–17,162 bp), and two inverted repeats (IRs; 4,288–4,413 bp). A total of 107 unique genes were identified from the three chloroplast genomes, including 70 protein-coding genes, 33 tRNAs, and 4 rRNAs. Compared with other angiosperms, the three chloroplast genomes were significantly reduced in overall length due to an unusual 16–24 kb shrinkage of IR regions and loss of the 11 genes which encoded subunits of NADH dehydrogenase. Genome-wide comparisons revealed similarities and variations between the three species and others. Phylogenetic analysis based on the three chloroplast genomes supported the opinion that Zygophyllaceae belonged to Zygophyllales in Fabids, and Z. xanthoxylon and Z. fabago belonged to Zygophyllum. The genome-wide comparisons revealed the similarity and variations between the chloroplast genomes of the three Zygophylloideae species and other plant species. This study provides a valuable molecular biology evidence for further studies of phylogenetic status of Zygophyllaceae.

Chloroplast genome analysis of box-ironbark Eucalyptus

Eucalyptus L'Hérit. (Myrtaceae) is a taxonomically complex and highly speciose genus that dominates much of Australia's woody vegetation. However, very little information is available about the molecular biology and chloroplast diversity of certain groups, such as Eucalyptus section Adnataria, which is found in many woodland habitats of eastern Australia. We report four new complete chloroplast genomes of Eucalyptus, including three genomes from species previously lacking any chloroplast reference sequences. Plastomes of E. albens, E. conica, E. crebra and E. melliodora assembled using a de novo approach were shown to be largely identical to each other, and similar in size and structure to previously published chloroplast genomes from Eucalyptus. A total of 132 genes (114 single-copy genes and 18 duplicated genes in the IR regions) were identified, and shown to be highly conserved in terms of gene order, content and organization. Slightly higher divergence in the intergenic spacers was identified through comparative genomic analyses. Chloroplast sequences of 35 additional individuals representing 12 species were assembled using a reference guided approach. Rates of nucleotide substitution varied among the protein coding genes, with 17 genes under possible positive selection, and 29 invariant genes. Phylogenetic analysis of either the whole reconstructed plastome sequences or the individual genes revealed extreme discordance with expected species boundaries or higher-level relationships. Plastome relationships were better predicted by geography than by nuclear DNA or taxonomic relationships, suggesting a substantial influence of gene flow over and above the effects of incomplete lineage sorting. These results provide resources for future research and valuable insights into the prevalence of interspecific gene flow among Eucalyptus species.

Frequent gene flow blurred taxonomic boundaries of sections in Lilium L. (Liliaceae)

Gene flow between species may last a long time in plants. Reticulation inevitably causes difficulties in phylogenetic reconstruction. In this study, we looked into the genetic divergence and phylogeny of 20 Lilium species based on multilocus analyses of 8 genes of chloroplast DNA (cpDNA), the internally transcribed nuclear ribosomal DNA (nrITS) spacer and 20 loci extracted from the expressed sequence tag (EST) libraries of L. longiflorum Thunb. and L. formosanum Wallace. The phylogeny based on the combined data of the maternally inherited cpDNA and nrITS was largely consistent with the taxonomy of Lilium sections. This phylogeny was deemed the hypothetical species tree and uncovered three groups, i.e., Cluster A consisting of 4 taxa from the sections Pseudolirium and Liriotypus, Cluster B consisting of the 4 taxa from the sections Leucolirion, Archelirion and Daurolirion, and Cluster C comprising 10 taxa mostly from the sections Martagon and Sinomartagon. In contrast, systematic inconsistency occurred across the EST loci, with up to 19 genes (95%) displaying tree topologies deviating from the hypothetical species tree. The phylogenetic incongruence was likely attributable to the frequent genetic exchanges between species/sections, as indicated by the high levels of genetic recombination and the IMa analyses with the EST loci. Nevertheless, multilocus analysis could provide complementary information among the loci on the species split and the extent of gene flow between the species. In conclusion, this study not only detected frequent gene flow among Lilium sections that resulted in phylogenetic incongruence but also reconstructed a hypothetical species tree that gave insights into the nature of the complex relationships among Lilium species.

Phylogeographic evidence for two mesic refugia in a biodiversity hotspot

Phylogeographic studies of flora in species-rich south-western Australia point to complex evolutionary histories, reflecting patterns of persistence and resilience to climatic changes during the Pleistocene. We asked whether coastal areas of the mid-west and south, as well as granite outcrops and inland ranges, have acted as major refugia within this region during Pleistocene climatic fluctuations by analysing phylogeographic patterns in the shrub Calothamnus quadrifidus R.Br. (Myrtaceae). We determined variation in chloroplast DNA data for 41 populations across the geographic range. Relationships and major clades were resolved using parsimony and Bayesian analyses. We tested for demographic and spatial expansion of the major clades and estimated clade divergence dates using an uncorrelated, lognormal relaxed clock based on two conservative chloroplast mutation rates. Two distinct phylogeographic clades were identified showing divergence during the Pleistocene, consistent with other phylogeographic studies of south-west Australian flora, emphasising the impact of climatic oscillations in driving divergence in this landscape. The southern clade was more diverse, having higher haplotype diversity and greater genetic structure, while the northern clade showed evidence of fluctuation in population size. Regions of high haplotype diversity with adjacent areas of low diversity observed in each clade indicated the locations of two coastal refugia: one on the south coast and another along the mid-west coast. This is the first evidence for major Pleistocene refugia using chloroplast genetic data in a common, widespread species from this region.

The major Australian cool temperate rainforest tree Nothofagus cunninghamii withstood Pleistocene glacial aridity within multiple regions: evidence from the chloroplast

Glacial aridity of the Pleistocene was inhospitable for the cool temperate rainforest tree Nothofagus cunninghamii over most of its current range in southeastern Australia, particularly in eastern Tasmania. A chloroplast DNA phylogeographic study was undertaken to investigate whether this species was likely to have survived in situ or conforms to a dispersal model of postglacial recovery. Twenty-three chloroplast haplotypes were identified by PCR-RFLP and direct sequencing of 2164 base pairs from 213 N. cunninghamii individuals collected in a range-wide survey. Fine-scale haplotype distribution was investigated using PCR-RFLP in eastern Tasmania. Deep chloroplast divergence occurred in N. cunninghamii. The single haplotype of the sister species, N. moorei, was nested among N. cunninghamii haplotypes. The distribution of N. cunninghamii haplotypes supports: multiple glacial refugia in coastal and inland western Tasmania, the centre of haplotype diversity; glacial survival in the central highlands of Victoria, corroborating pollen data; and the long-term occupation of eastern Tasmania because of the presence of a unique deeply diverged chloroplast lineage. Nothofagus cunninghamii withstood glacial aridity within multiple regions in apparently nonequable climates. This finding contributes to a growing understanding of how the resilience of temperate species during glacial periods has shaped modern biota.

Characterization of polymorphic microsatellite DNA markers in Banksia nivea, formerly Dryandra nivea

A genomic library was constructed and 10 novel polymorphic microsatellite markers were developed for Banksia nivea subsp. uliginosa, a woody shrub endemic to the southwest of Western Australia. Polymorphism was investigated for these markers in 30 individuals from one population located in restricted habitat at the base of the Whicher Range east of Busselton. Expected heterozygosity ranged from 0.124 to 0.898 and averaged 0.570. Observed heterozygosity ranged from 0.133 to 0.933 and averaged 0.544. The number of alleles per locus ranged from 2.0 to 16.0 and averaged 7.3. These markers also produced amplification products in three other Banksia species and will be used to assay genetic diversity in these species.

Identification of quantitative trait loci influencing foliar concentrations of terpenes and formylated phloroglucinol compounds in Eucalyptus nitens

Leaves of eucalypt species contain a variety of plant secondary metabolites, including terpenoids and formylated phloroglucinol compounds (FPCs). Both terpene and FPC concentrations are quantitative traits that can show large variation within a population and have been shown to be heritable. The molecular genetic basis of this variation is currently unknown. Progeny from a field trial of a three-generation mapping pedigree of Eucalyptus nitens were assayed for terpenes and FPCs. Quantitative trait loci (QTL) analyses were conducted using a map constructed from 296 markers to locate regions of the genome influencing foliar concentrations of these plant secondary compounds. A large number of significant QTL for 14 traits were located across nine linkage groups, with significant clustering of QTL on linkage groups 7, 8 and 9. As expected, QTL for biosynthetically related compounds commonly colocated, but QTL for unrelated monterpenes and FPCs also mapped closely together. Colocation of these QTL with mapped candidate genes from the various biosynthetic pathways, and subsequent use of these genes in association mapping, will assist in determining the causes of variation in plant secondary metabolites in eucalypts.

Genomic research in Eucalyptus

Eucalyptus L'Hérit. is a genus comprised of more than 700 species that is of vital importance ecologically to Australia and to the forestry industry world-wide, being grown in plantations for the production of solid wood products as well as pulp for paper. With the sequencing of the genomes of Arabidopsis thaliana and Oryza sativa and the recent completion of the first tree genome sequence, Populus trichocarpa, attention has turned to the current status of genomic research in Eucalyptus. For several eucalypt species, large segregating families have been established, high-resolution genetic maps constructed and large EST databases generated. Collaborative efforts have been initiated for the integration of diverse genomic projects and will provide the framework for future research including exploiting the sequence of the entire eucalypt genome which is currently being sequenced. This review summarises the current position of genomic research in Eucalyptus and discusses the direction of future research.

Phylogeography and divergence in the chloroplast genome of Western Australian Sandalwood (Santalum spicatum)

Western Australian sandalwood (Santalum spicatum) is widespread throughout Western Australia across the semiarid and arid regions. The diversity and phylogeographic patterns within the chloroplast genome of S. spicatum were investigated using restriction fragment length polymorphism analysis of 23 populations. The chloroplast diversity was structured into two main clades that were geographically separated, one centred in the southern (semiarid region) and the other in the northern (arid) region. Fragmentation due to climatic instability was identified as the most likely influence on the differentiation of the lineages. The lineage in the arid region showed a greater level of differentiation than that in the southern region, suggesting a higher level of gene flow or a more recent range expansion of sandalwood in the southern region. The phylogeographic pattern in the chloroplast genome is congruent with that detected in the nuclear genome, which identified different genetic influences between the regions and also suggested a more recent expansion of sandalwood in the southern region.

Long-term population isolation in the endangered tropical tree species Caesalpinia echinata Lam. revealed by chloroplast microsatellites

Habitat fragmentation represents the single most serious threat to the survival of tropical ecosystems. In formulating strategies to counteract the detrimental effects of fragmentation, knowledge of the levels and patterns of genetic diversity within and between natural populations is vital to the establishment of any conservation programme. We utilized polymorphic chloroplast microsatellite markers to analyse genetic diversity in populations of the endangered tropical tree Caesalpinia echinata Lam. representing the entire extant range of the species. Levels of within-population diversity were low, with only two of seven populations studied displaying any variation. The vast majority of the genetic variation was partitioned between geographical regions (36%) and between populations within regions (55%). These levels of genetic structuring, coupled with a calculated pollen-to-seed flow ratio of approximately 6.7:1, suggest that there has been little gene flow between the three major geographical regions over an extended period. Thus, the current tripartite distribution of the species is more consistent with the existence of separate glacial refugia, rather than reflecting any anthropogenic effects.

Regional genetic differentiation in Western Australian sandalwood (Santalum spicatum) as revealed by nuclear RFLP analysis

Western Australian sandalwood, Santalum spicatum, is widespread in the semi-arid and arid regions of Western Australia, and there is some morphological variation suggestive of two ecotypes. The level and structuring of genetic diversity within the species was investigated using anonymous nuclear RFLP loci. Santalum spicatum showed moderate levels of genetic diversity compared to other Australian tree species. The northern populations in the arid region showed greater levels of diversity and less population differentiation than the southern populations in the semi-arid region due to differences in the distribution of rare alleles. Equilibrium between drift and gene flow in the northern populations indicated that they have been established for a long period of time with stable conditions conducive to gene flow. In contrast, the southern populations showed a relationship between drift and gene flow indicative of a pattern of fragmentation and isolation where drift has greater effect than gene flow. The different patterns of diversity suggest that the ecotypes in the two regions have been subject to differences in the relative influences of drift and gene flow during their evolutionary history.

Genetic structure and reproduction dynamics of Salix reinii during primary succession on Mount Fuji, as revealed by nuclear and chloroplast microsatellite analysis

The early stage of volcanic desert succession is underway on the southeastern slope of Mount Fuji. We used markers of nuclear microsatellites (simple sequence repeats; SSR) and chloroplast microsatellites (cpSSR) to investigate the population genetic structure and reproduction dynamics of Salix reinii, one of the dominant pioneer shrubs in this area. The number of S. reinii genets in a patch and the area of the largest genet within the patch increased with patch area, suggesting that both clonal growth and seedling recruitment are involved in the reproduction dynamics of S. reinii. Five polymorphic cpSSR markers were developed for S. reinii by sequencing the noncoding regions between universal sequences in the chloroplast genome. Nineteen different cpSSR haplotypes were identified, indicating that S. reinii pioneer genets were created by the long-distance dispersal of seeds originating from different mother genets around the study site, where all vegetation was destroyed during the last eruption. Furthermore, the clustered distributions of different haplotypes within each patch or plot suggested that newly colonized genets tended to be generated from seeds dispersed near the initially established mother genets. These results revealed that the establishment of the S. reinii population on the southeastern slope of Mount Fuji involved two sequential modes of seed dispersal: long-distance dispersal followed by short-distance dispersal.

Genetic differentiation between mallee and tree forms in the Eucalyptus loxophleba complex

The Eucalyptus loxophleba complex comprises two tree taxa (E. loxophleba ssp. loxophleba and ssp. supralaevis), two mallee taxa (E. loxophleba ssp. lissophloia and ssp. gratiae), and three rare mallee species, of which one, E. blaxellii, is included in this study. The genetic diversity and phylogenetic relationships between the taxa, particularly in relation to habit, were assessed using anonymous nuclear RFLP loci. The level of diversity in the taxa was high and similar to that detected in other eucalypt species. The populations showed low differentiation at both the subspecies and the species levels. Phylogenetic relationships showed some genetic separation between the tree and mallee habit but no separation of the two taxa within the tree habit or within the mallee habit. The genetic analysis does not support the recognition of E. gratiae as a separate species. The geographically restricted E. blaxellii showed similar levels of diversity to populations of the other widespread taxa of E. loxophleba.

Chloroplast DNA haplotype variation and population differentiation in Sorbus aucuparia L. (Rosaceae: Maloideae)

Intra-specific chloroplast DNA (cpDNA) variation was studied in Sorbus aucuparia L., an entomophilous, mid-or early successional tree producing fleshy fruits. Eight PCR-amplified fragments of the chloroplast genome were screened for restriction fragment length polymorphisms, using one or two 4 bp-cutter restriction endonucleases. cpDNA variation was investigated on two geographical scales: (1) among four regions in France and Belgium; and (2) within the Belgian region. A total of 150 individuals from six populations were analysed. Fourteen polymorphisms were detected in six of the cpDNA fragments. All polymorphisms probably resulted from insertions or deletions, and allowed the identification of 12 haplotypes. The level of genetic differentiation computed on the basis of haplotype frequencies was similar on the two geographical scales considered (G(STc) = 0.286 among regions, G(STc) = 0.259 among populations within the Belgian region). These values are much lower than those obtained in nine previously studied temperate tree species, which are all wind-pollinated, late-successional species producing dry fruits. These results might primarily be accounted for by the contrasting life history traits of S. aucuparia. In order to obtain insights into the relative contribution of pollen and seeds to gene flow, G(STc) was also compared with previously obtained G(ST) estimates based on allozyme data.

Divergence in the chloroplast genome and nuclear rDNA of the rare western australian plant lambertia orbifolia gardner (Proteaceae)

The population genetic structure of the Australian plant Lambertia orbifolia was investigated for chloroplast DNA (cpDNA) and rDNA based on restriction fragment length polymorphism. Variation was assessed in 14-20 individuals from six populations with probes covering the majority of the chloroplast genome and the whole rRNA gene unit. For cpDNA, eight mutations were detected which were distributed over five haplotypes. Nucleotide diversity in the species was high and the majority of this diversity was distributed between populations with diversity within populations restricted to a single population. There was significant differentiation between the two regions in the species distribution with the Narrikup region being distinguished by a single haplotype that was characterized by six unique mutations. Variation in rDNA was detected with three gene length variants present in most individuals. However, the Narrikup region was characterized by homogenization of the gene unit to a single length variant in all individuals. The divergence of the Narrikup region suggests that the disjunction in the species distribution has been present for a long time and the two regions represent separate evolutionary lineages.

Molecular phylogeography, intraspecific variation and the conservation of tree species

Tree species are becoming the focus of increasing conservation concern, with some 9000 species now threatened globally. Studies of intraspecific variation can contribute to the development of conservation strategies, by identifying appropriate units for conservation. The recent application of molecular techniques to a variety of tree species has highlighted a far higher degree of population differentiation than indicated by previous isozyme analyses, a result consistent with theoretical predictions. Analysis of the geographic distribution of cpDNA lineages has also enabled current patterns of population differentiation to be related to postglacial migration routes from different forest refugia. Such results highlight the importance of refugial areas for conservation of intraspecific variation in tree species.

Abundant mitochondrial genome diversity, population differentiation and convergent evolution in pines

We examined mitochondrial DNA polymorphisms via the analysis of restriction fragment length polymorphisms in three closely related species of pines from western North America: knobcone (Pinus attenuata Lemm.), Monterey (P. radiata D. Don), and bishop (P. muricata D. Don). A total of 343 trees derived from 13 populations were analyzed using 13 homologous mitochondrial gene probes amplified from three species by polymerase chain reaction. Twenty-eight distinct mitochondrial DNA haplotypes were detected and no common haplotypes were found among the species. All three species showed limited variability within populations, but strong differentiation among populations. Based on haplotype frequencies, genetic diversity within populations (HS) averaged 0.22, and population differentiation (GST and theta) exceeded 0.78. Analysis of molecular variance also revealed that >90% of the variation resided among populations. For the purposes of genetic conservation and breeding programs, species and populations could be readily distinguished by unique haplotypes, often using the combination of only a few probes. Neighbor-joining phenograms, however, strongly disagreed with those based on allozymes, chloroplast DNA, and morphological traits. Thus, despite its diagnostic haplotypes, the genome appears to evolve via the rearrangement of multiple, convergent subgenomic domains.

Chloroplast DNA phylogeography of the argan tree of Morocco

Polymorphisms in the chloroplast genome of the argan tree (Sapotaceae), an endemic species of south-western Morocco, have been detected by restriction site studies of PCR-amplified fragments. A total of 12 chloroplast DNA (cpDNA) and two mitochondrial DNA (mtDNA) fragments were amplified and digested with a single restriction enzyme (HinfI). Polymorphisms were identified in six of the cpDNA fragments, whereas no mtDNA polymorphisms were detected in a survey of 95 individuals from 19 populations encompassing most of the natural range of the species. The cpDNA polymorphisms allowed the identification of 11 haplotypes. Two lineages, one in the south-east and the other in the north-west, divide the range of the argan tree into two distinct areas. The level of genetic differentiation measured at the haplotype level (GSTc = 0.60) (i.e. with unordered haplotypes) was smaller than when phylogenetic relationships were taken into account (NSTc = 0.71-0.74) (ordered haplotypes), indicating that population history must be considered in the study of the geographical distribution of cpDNA lineages in this species. If contrasted with the level of nuclear genetic differentiation measured in a previous study with isozymes (GSTn = 0.25), the results indicate a relatively high level of gene flow by seeds, or conversely a relatively low level of gene flow by pollen, as compared with other tree species. Goats and camels could have played an important role in disseminating the fruits of this tree.

A population genetic analysis of chloroplast DNA in Phacelia

Hierarchical sampling from populations, incipient and recognized varieties within Phacelia dubia and P. maculata has revealed high levels of intraspecific polymorphism in chloroplast DNA. Much of the variation is partitioned between populations as evidenced by population-specific variants at fixation in all three populations of P. dubia var. interior and in both populations of P. maculata. Nine of 16 populations were polymorphic for cpDNA haplotypes. A total of 16 haplotypes was found in a sample of 106 individuals; the most common occurred in eight of the 16 populations and in 31 per cent of the individuals in the entire sample. A phylogenetic analysis revealed four basic plastome types. The two major groups of plastomes were separated by four independent base-pair mutations which suggests an ancient split in the evolution of plastid genomes. Representatives from each major plastome division were found in each of five populations spanning two allopatric varieties of P. dubia. The geographical distribution of haplotypes and lack of evidence for recent admixture argue against migration as a source of the polymorphism. It is more likely that the current taxonomic varieties are descendants of a polymorphic common ancestor.

An integrated genetic linkage map for eucalypts using RFLP, RAPD and isozyme markers

An integrated genetic linkage map for E. nitens was constructed in an outbred three-generation pedigree. Analysis of 210 RFLP, 125 RAPD and 4 isozyme loci resulted in 330 markers linked in 12 linkage groups covering 1462 cM (n=11 in eucalypts). The 12th linkage group is comprised of only 5 markers and will probably coalesce with another linkage group when further linked loci are located. Co-dominant RFLP loci segregating in both parents were used to integrate linkages identified in the male and female parents. Differences in recombination frequencies in the two parents were observed for a number of pairs of loci, and duplication of sequences was identified both within and between linkage groups. The markers were distributed randomly across the genome except for the RFLPs in linkage group 10 and for some loci showing segregation distortion, which were clustered into three regions of the map. The use of a large number of co-dominant RFLP loci in this map enables it to be used in other pedigrees of E. nitens and forms a basis for the detection and location of QTL in E. nitens and other eucalypt species.

Detection and inheritance of RFLPs in Eucalyptus nitens

The level of polymorphism using genomic and cDNA probes with a number of restriction enzymes and the inheritance of the RFLP loci was investigated in E. nitens. The polymorphism detected with 366 genomic and cDNA probes and three to six restriction enzymes was analysed in three-generation outbred pedigrees. No difference in the level of polymorphism detected with genomic versus cDNA probes was observed. There was a difference in the efficiency of detection of polymorphism with six different restriction enzymes, with three of the enzymes (BglII, DraI and EcoRI) showing substantially more polymorphism than the others. There was no significant correlation between the size of the DNA fragments generated by the enzymes and the detection of polymorphism. Several cases of restriction-site mutations resulting in a polymorphism were observed. The inheritance of 69 loci was analysed in two pedigrees resulting from interpopulational crosses. The majority of the loci segregated according to expected ratios with distortion observed in only 3% of loci. Probes from the cDNA library detected a greater proportion of loci with more than two alleles than did probes from the genomic library. The high polymorphism, large number of alleles, and ease of interpretation of RFLPs in E. nitens means that they will be useful in a range of applications such as genetic linkage maps and paternity analysis.