Phylogenetic relationships in Taxodiaceae and Cupressaceae sensu stricto based on matK gene, chlL gene, trnL-trnF IGS region, and trnL intron sequences

Squamabietenols A-F, undescribed abietane-O-abietane dimeric diterpenoids from the ornamental conifer Juniperus squamata and their ATP-citrate lyase inhibitory activities

Six undescribed naturally occurring abietane-O-abietane dimers (squamabietenols A-F) together with one 3,4-seco-totarane-type, a pimarane-type, and 17 related known mono-/dimeric diterpenoids were isolated and characterized from the needles and twigs of the ornamental conifer Juniperus squamata. The undescribed structures and their absolute configurations were established by extensive spectroscopic methods, GIAO NMR calculations with DP4+ probability analyses, and ECD calculations. Squamabietenols A and B showed significant inhibitory effects against ATP-citrate lyase (ACL, a novel drug target for hyperlipidemia and other metabolic disorders), with IC₅₀ values of 8.82 and 4.49 μM, respectively. A molecular docking study corroborated the findings by highlighting the interactions between the bioactive compounds and the ACL enzyme (binding affinities: -7.1 to -9.0 kcal/mol). The unique abietane-O-abietane dimeric diterpenoids are quite rare in the vegetable kingdom, and they are of chemotaxonomic significance for the Cupressaceae family.

A DNA barcode reference library for the native woody seed plants of Japan

DNA barcode databases are increasingly available for a range of organisms, facilitating the wide application of DNA barcode-based studies. Here we announce the development of a comprehensive DNA barcode reference library of Japanese native woody seed plants representing 43 orders, 99 families, 303 genera and 834 species, and comprising 77.3% of the genera and 72.2% of the species of native woody seed plants in Japan. A total of 6216 plant specimens were collected from 223 sites across the subtropical, temperate, boreal and alpine biomes in Japan with most species represented by multiple accessions. This reference library utilized three chloroplast DNA regions (rbcL, trnH-psbA and matK) and consists of 14,403 barcode sequences. Individual regions varied in their identification rates, with species-level and genus-level rates for rbcL, trnH-psbA and matK based on blast being 57.4%/96.2%, 78.5%/99.1% and 67.8%/98.1%, respectively. Identification rates were higher using region combinations, with total species-level rates for two region combinations (rbcL & trnH-psbA, rbcL & matK and trnH-psbA & matK) ranging between 90.6% and 95.8%, and for all three regions being equal to 98.6%. Genus-level identification rates were even higher, ranging between 99.7% and 100% for two region combinations and being 100% for the three regions. These results indicate that this DNA barcode reference library is an effective resource for investigations of native woody seed plants in Japan using DNA barcodes and provides a useful template for the development of libraries for other components of the Japanese flora.

Phylogeny and evolution of Cupressaceae: Updates on intergeneric relationships and new insights on ancient intergeneric hybridization

After the merger of the former Taxodiaceae and Cupressaceae s.s., currently the conifer family Cupressaceae (sensu lato) comprises seven subfamilies and 32 genera, most of which are important components of temperate and mountainous forests. With the exception of a recently published genus-level phylogeny of gymnosperms inferred from sequence analysis of 790 orthologs, previous phylogenetic studies of Cupressaceae were based mainly on morphological characters or a few molecular markers, and did not completely resolve the intergeneric relationships. In this study, we reconstructed a robust and well-resolved phylogeny of Cupressaceae represented by all 32 genera, using 1944 genes (Orthogroups) generated from transcriptome sequencing. Reticulate evolution analyses detected a possible ancient hybridization that occurred between ancestors of two subclades of Cupressoideae, including Microbiota-Platycladus-Tetraclinis (MPT) and Juniperus-Cupressus-Hesperocyparis-Callitropsis-Xanthocyparis (JCHCX), although both concatenation and coalescent trees are highly supported. Moreover, divergence time estimation and ancestral area reconstruction indicate that Cupressaceae very likely originated in Asia in the Triassic, and geographic isolation caused by continental separation drove the vicariant evolution of the two subfamilies Cupressoideae and Callitroideae in the northern and southern hemispheres, respectively. Evolutionary analyses of some morphological characters suggest that helically arranged linear-acicular leaves and imbricate bract-scale complexes represent ancestral states, and the shift from linear-acicular leaves to scale-like leaves was associated with the shift from helical to decussate arrangement. Our study sheds new light on phylogeny and evolutionary history of Cupressaceae, and strongly suggests that both dichotomous phylogenetic and reticulate evolution analyses be conducted in phylogenomic studies.

A protocol for Agrobacterium-mediated transformation of Japanese cedar, Sugi (Cryptomeria japonica D. Don) using embryogenic tissue explants

Sugi (Cryptomeria japonica D. Don) is the most important afforestation coniferous tree in Japan. Coniferous trees normally have a long juvenile period and require a long cultivation time for breeding. Through a traditional breeding project that began in the 1950s, first generation plus trees with excellent traits were selected primarily from artificial forests and used as seedlings. Recently, the second generation plus trees obtained by crossing between plus trees have been selected. In light of this situation, the improvement of Sugi by a transgenic approach is effective in terms of shortening the breeding period compared with conventional crossing-dependent approaches. There are three key points to an efficient Agrobacterium-mediated transformation system: (1) establishment of explants with high regeneration ability, (2) optimal co-cultivation conditions for explants and Agrobacterium, and (3) efficient elimination of Agrobacterium. Here we describe a protocol for Agrobacterium-mediated transformation of Sugi that meets the above criteria using embryogenic tissues as explants isolated from immature seeds obtained by crossing.

High-density genetic map construction and quantitative trait loci identification for growth traits in (Taxodium distichum var. distichum × T. mucronatum) × T. mucronatum

BACKGROUND

'Zhongshanshan' is the general designation for the superior interspecific hybrid clones of Taxodium species, which is widely grown for economic and ecological purposes in southern China. Growth is the priority objective in 'Zhongshanshan' tree improvement. A high-density linkage map is vital to efficiently identify key quantitative trait loci (QTLs) that affect growth.

RESULTS

In total, 403.16 Gb of data, containing 2016,336 paired-end reads, was obtained after preprocessing. The average sequencing depth was 28.49 in T. distichum var. distichum, 25.18 in T. mucronatum, and 11.12 in each progeny. In total, 524,662 high-quality SLAFs were detected, of which 249,619 were polymorphic, and 6166 of the polymorphic markers met the requirements for use in constructing a genetic map. The final map harbored 6156 SLAF markers on 11 linkage groups, and was 1137.86 cM in length, with an average distance of 0.18 cM between adjacent markers. Separate QTL analyses of traits in different years by CIM detected 7 QTLs. While combining multiple-year data, 13 QTLs were detected by ICIM. 5 QTLs were repeatedly detected by the two methods, and among them, 3 significant QTLs (q6-2, q4-2 and q2-1) were detected in at least two traits. Bioinformatic analysis discoveried a gene annotated as a leucine-rich repeat receptor-like kinase gene within q4-2.

CONCLUSIONS

This map is the most saturated one constructed in a Taxodiaceae species to date, and would provide useful information for future comparative mapping, genome assembly, and marker-assisted selection.

An overview of extant conifer evolution from the perspective of the fossil record

PREMISE OF THE STUDY

Conifers are an important living seed plant lineage with an extensive fossil record spanning more than 300 million years. The group therefore provides an excellent opportunity to explore congruence and conflict between dated molecular phylogenies and the fossil record.

METHODS

We surveyed the current state of knowledge in conifer phylogenetics to present a new time-calibrated molecular tree that samples ~90% of extant species diversity. We compared phylogenetic relationships and estimated divergence ages in this new phylogeny with the paleobotanical record, focusing on clades that are species-rich and well known from fossils.

KEY RESULTS

Molecular topologies and estimated divergence ages largely agree with the fossil record in Cupressaceae, conflict with it in Araucariaceae, and are ambiguous in Pinaceae and Podocarpaceae. Molecular phylogenies provide insights into some fundamental questions in conifer evolution, such as the origin of their seed cones, but using them to reconstruct the evolutionary history of specific traits can be challenging.

CONCLUSIONS

Molecular phylogenies are useful for answering deep questions in conifer evolution if they depend on understanding relationships among extant lineages. Because of extinction, however, molecular datasets poorly sample diversity from periods much earlier than the Late Cretaceous. This fundamentally limits their utility for understanding deep patterns of character evolution and resolving the overall pattern of conifer phylogeny.

The Complete Chloroplast Genome Sequence of a Relict Conifer Glyptostrobus pensilis: Comparative Analysis and Insights into Dynamics of Chloroplast Genome Rearrangement in Cupressophytes and Pinaceae

Glyptostrobus pensilis, belonging to the monotypic genus Glyptostrobus (Family: Cupressaceae), is an ancient conifer that is naturally distributed in low-lying wet areas. Here, we report the complete chloroplast (cp) genome sequence (132,239 bp) of G. pensilis. The G. pensilis cp genome is similar in gene content, organization and genome structure to the sequenced cp genomes from other cupressophytes, especially with respect to the loss of the inverted repeat region A (IRA). Through phylogenetic analysis, we demonstrated that the genus Glyptostrobus is closely related to the genus Cryptomeria, supporting previous findings based on physiological characteristics. Since IRs play an important role in stabilize cp genome and conifer cp genomes lost different IR regions after splitting in two clades (cupressophytes and Pinaceae), we performed cp genome rearrangement analysis and found more extensive cp genome rearrangements among the species of cupressophytes relative to Pinaceae. Additional repeat analysis indicated that cupressophytes cp genomes contained less potential functional repeats, especially in Cupressaceae, compared with Pinaceae. These results suggested that dynamics of cp genome rearrangement in conifers differed since the two clades, Pinaceae and cupressophytes, lost IR copies independently and developed different repeats to complement the residual IRs. In addition, we identified 170 perfect simple sequence repeats that will be useful in future research focusing on the evolution of genetic diversity and conservation of genetic variation for this endangered species in the wild.

Whole genome duplication in coast redwood (Sequoia sempervirens) and its implications for explaining the rarity of polyploidy in conifers

Polyploidy is common and an important evolutionary factor in most land plant lineages, but it is rare in gymnosperms. Coast redwood (Sequoia sempervirens) is one of just two polyploid conifer species and the only hexaploid. Evidence from fossil guard cell size suggests that polyploidy in Sequoia dates to the Eocene. Numerous hypotheses about the mechanism of polyploidy and parental genome donors have been proposed, based primarily on morphological and cytological data, but it remains unclear how Sequoia became polyploid and why this lineage overcame an apparent gymnosperm barrier to whole-genome duplication (WGD). We sequenced transcriptomes and used phylogenetic inference, Bayesian concordance analysis and paralog age distributions to resolve relationships among gene copies in hexaploid coast redwood and close relatives. Our data show that hexaploidy in coast redwood is best explained by autopolyploidy or, if there was allopolyploidy, it happened within the Californian redwood clade. We found that duplicate genes have more similar sequences than expected, given the age of the inferred polyploidization. Conflict between molecular and fossil estimates of WGD can be explained if diploidization occurred very slowly following polyploidization. We extrapolate from this to suggest that the rarity of polyploidy in gymnosperms may be due to slow diploidization in this clade.

Phylogeny and divergence times of gymnosperms inferred from single-copy nuclear genes

Phylogenetic reconstruction is fundamental to study evolutionary biology and historical biogeography. However, there was not a molecular phylogeny of gymnosperms represented by extensive sampling at the genus level, and most published phylogenies of this group were constructed based on cytoplasmic DNA markers and/or the multi-copy nuclear ribosomal DNA. In this study, we use LFY and NLY, two single-copy nuclear genes that originated from an ancient gene duplication in the ancestor of seed plants, to reconstruct the phylogeny and estimate divergence times of gymnosperms based on a complete sampling of extant genera. The results indicate that the combined LFY and NLY coding sequences can resolve interfamilial relationships of gymnosperms and intergeneric relationships of most families. Moreover, the addition of intron sequences can improve the resolution in Podocarpaceae but not in cycads, although divergence times of the cycad genera are similar to or longer than those of the Podocarpaceae genera. Our study strongly supports cycads as the basal-most lineage of gymnosperms rather than sister to Ginkgoaceae, and a sister relationship between Podocarpaceae and Araucariaceae and between Cephalotaxaceae-Taxaceae and Cupressaceae. In addition, intergeneric relationships of some families that were controversial, and the relationships between Taxaceae and Cephalotaxaceae and between conifers and Gnetales are discussed based on the nuclear gene evidence. The molecular dating analysis suggests that drastic extinctions occurred in the early evolution of gymnosperms, and extant coniferous genera in the Northern Hemisphere are older than those in the Southern Hemisphere on average. This study provides an evolutionary framework for future studies on gymnosperms.

Climate, not Aboriginal landscape burning, controlled the historical demography and distribution of fire-sensitive conifer populations across Australia

Climate and fire are the key environmental factors that shape the distribution and demography of plant populations in Australia. Because of limited palaeoecological records in this arid continent, however, it is unclear as to which factor impacted vegetation more strongly, and what were the roles of fire regime changes owing to human activity and megafaunal extinction (since ca 50 kya). To address these questions, we analysed historical genetic, demographic and distributional changes in a widespread conifer species complex that paradoxically grows in fire-prone regions, yet is very sensitive to fire. Genetic demographic analysis showed that the arid populations experienced strong bottlenecks, consistent with range contractions during the Last Glacial Maximum (ca 20 kya) predicted by species distribution models. In southern temperate regions, the population sizes were estimated to have been mostly stable, followed by some expansion coinciding with climate amelioration at the end of the last glacial period. By contrast, in the flammable tropical savannahs, where fire risk is the highest, demographic analysis failed to detect significant population bottlenecks. Collectively, these results suggest that the impact of climate change overwhelmed any modifications to fire regimes by Aboriginal landscape burning and megafaunal extinction, a finding that probably also applies to other fire-prone vegetation across Australia.

Extended linkage disequilibrium in noncoding regions in a conifer, Cryptomeria japonica

We measured linkage disequilibrium in mostly noncoding regions of Cryptomeria japonica, a conifer belonging to Cupressaceae. Linkage disequilibrium was extensive and did not decay even at a distance of 100 kb. The average estimate of the population recombination rate per base pair was 1.55 × 10(-5) and was <1/70 of that in the coding regions. We discuss the impact of low recombination rates in a large part of the genome on association studies.

Multilocus patterns of nucleotide polymorphism and demographic change in Taxodium distichum (Cupressaceae) in the lower Mississippi River alluvial valley

PREMISE OF THE STUDY

Studies of the geographic patterns of genetic variation can give important insights into the past population structure of species. Our study species, Taxodium distichum L. (bald-cypress), prefers riparian and wetland habitats and is widely distributed in southeastern North America and Mexico. We compared the genetic variation of T. distichum with that of its close relative, Cryptomeria japonica, which is endemic to Japan. •

METHODS

Nucleotide polymorphisms of T. distichum in the lower Mississippi River alluvial valley, USA, were examined at 10 nuclear loci. •

KEY RESULTS

The average nucleotide diversity at silent sites, π(sil), across the 10 loci in T. distichum was higher than that of C. japonica (π(sil) = 0.00732 and 0.00322, respectively). In T. distichum, Tajima's D values were each negative at 9 out of 10 loci, which suggests a recent population expansion. Maximum-likelihood and Bayesian estimations of the exponential population growth rate (g) of T. distichum populations indicated that this species had expanded approximately at the rate of 1.7-1.0 × 10(-6) per year in the past. •

CONCLUSIONS

Taxodium distichum had significantly higher nucleotide variation than C. japonica, and its patterns of polymorphism contrasted strikingly with those of the latter, which previously has been inferred to have experienced a reduction in population size.

Diversification and biogeography of Juniperus (Cupressaceae): variable diversification rates and multiple intercontinental dispersals

• A central aim of biogeography is to understand when and how modern patterns of species diversity and distribution developed. Many plant groups have disjunct distributions within the Northern Hemisphere, but among these very few have been studied that prefer warm semi-arid habitats. • Here we examine the biogeography and diversification history of Juniperus, which occurs in semi-arid habitats through much of the Northern Hemisphere. A phylogeny was generated based on > 10,000 bp of cpDNA for 51 Juniperus species plus many outgroups. Phylogenies based on fewer species were also constructed based on nuclear internal transcribed spacer (nrITS) and combined nrITS/cpDNA data sets to check for congruence. Divergence time-scales and ancestral distributions were further inferred. • Both long dispersal and migration across land bridges probably contributed to the modern range of Juniperus, while long-term climatic changes and the uplift of the Qinghai-Tibetan plateau probably drove its diversification. Diversification apparently slowed down during climate-stable period of the Oligocene, and then speeded up from the Miocene onwards. • Juniperus probably originated in Eurasia, and was a part of the south Eurasian Tethyan vegetation of the Eocene to Oligocene. It reached America once at this time, once in the Miocene and once more recently.

Genetic diversity and structure of natural fragmented Chamaecyparis obtusa populations as revealed by microsatellite markers

The genetic diversity and population structure of hinoki (Chamaecyparis obtusa) in Japan were investigated by examining the distribution of alleles at 13 microsatellite loci in 25 natural populations from Iwaki in northern Japan to Yakushima Island in southern Japan. On average, 26.9 alleles per locus were identified across all populations and 4.0% of the genetic variation was retained among populations (G(ST) = 0.040). According to linkage disequilibrium analysis, estimates of effective population size and detected evidence of bottleneck events, the genetic diversity of some populations may have declined as a result of fragmentation and/or over-exploitation. The central populations located in the Chubu district appear to have relatively large effective population sizes, while marginal populations, such as the Yakushima, Kobayashi and Iwaki populations, have smaller effective population sizes and are isolated from the other populations. Microsatellite analysis revealed the genetic uniqueness of the Yakushima population. Although genetic differentiation between populations was low, we detected a gradual cline in the genetic structure and found that locus Cos2619 may be non-neutral with respect to natural selection.

A frameshift mutation of the chloroplast matK coding region is associated with chlorophyll deficiency in the Cryptomeria japonica virescent mutant Wogon-Sugi

Wogon-Sugi has been reported as a cytoplasmically inherited virescent mutant selected from a horticultural variety of Cryptomeria japonica. Although previous studies of plastid structure and inheritance indicated that at least some mutations are encoded by the chloroplast genome, the causative gene responsible for the primary chlorophyll deficiency in Wogon-Sugi, has not been identified. In this study, we identified this gene by genomic sequencing of chloroplast DNA and genetic analysis. Chloroplast DNA sequencing of 16 wild-type and 16 Wogon-Sugi plants showed a 19-bp insertional sequence in the matK coding region in the Wogon-Sugi. This insertion disrupted the matK reading frame. Although an indel mutation in the ycf1 and ycf2 coding region was detected in Wogon-Sugi, sequence variations similar to that of Wogon-Sugi were also detected in several wild-type lines, and they maintained the reading frame. Genetic analysis of the 19 bp insertional mutation in the matK coding region showed that it was found only in the chlorophyll-deficient sector of 125 full-sibling seedlings. Therefore, the 19-bp insertion in the matK coding region is the most likely candidate at present for a mutation underlying the Wogon-Sugi phenotype.

Adaptive and slightly deleterious evolution in a conifer, Cryptomeria japonica

In order to evaluate effects of the population structure and natural selection on organisms having long generation times, we surveyed DNA polymorphisms at five loci encoding 9-cis-epoxycarotenoid dioxygenase (NCED), ammonium transporter, calmodulin, aquaporin, and the second major allergen with polymethylgalacturonase enzyme activity in the pollen (Cryj2) in a conifer, Cryptomeria japonica. The average nucleotide diversity at silent sites across 12 loci including the previously analyzed seven loci was 0.0044. The population recombination rate (4Nr, where N and r are the effective population size and recombination rate per base per generation, respectively) was estimated as 0.00046 and a slow reduction in the population size was indicated, according to the maximum likelihood method implemented in LAMARC. At NCED, the McDonald-Kreitman (MK) test revealed an excess of replacement polymorphisms, suggesting contributions of slightly deleterious mutations. In contrast, the MK test revealed an excess of replacement divergence at Cryj2 and a maximum likelihood approach using the PAML package revealed that certain amino acid sites had a nonsynonymous/synonymous substitution rate ratio (omega) > 4.0, indicating adaptive evolution at this locus. The overall analysis of the 12 loci suggested that adaptive, neutral, and slightly deleterious evolution played important roles in the evolution of C. japonica.

Molecular phylogeny and biogeography of Picea (Pinaceae): Implications for phylogeographical studies using cytoplasmic haplotypes

The center of diversity is not necessarily the place of origin, as has been established by many plant molecular phylogenies. Picea is a complicated but very important genus in coniferous forests of the Northern Hemisphere, with a high species diversity in Asia. Its phylogeny and biogeography were investigated here using sequence analysis of the paternally inherited chloroplast trnC-trnD and trnT-trnF regions and the maternally inherited mitochondrial nad5 intron 1. We found that the North American P. breweriana and P. sitchensis were basal to the other spruces that were further divided into three clades in the cpDNA phylogeny, and that the New World species harbored four of five mitotypes detected, including two ancestral ones and three endemics. These results, combined with biogeographic analyses using DIVA and MacClade and fossil evidence, suggest that Picea originated in North America, and that its present distribution could stem from two times of dispersal from North America to Asia by the Beringian land bridge, and then from Asia to Europe. Most of the northeastern Asian species and the European P. abies could arise from a recent radiation given the very low interspecific genetic differentiation and pure mitotype of them. Considering that the ancestral mtDNA polymorphism can be preserved in many descendant species, even distantly related ones, we suggest that more species, at least the closely related ones, should be sampled in the phylogeographical study using cytoplasmic haplotypes if possible. In addition, we also discussed the evolution and phylogenetic utility of morphological characters in Picea.

Relaxation of Functional Constraint on Light-Independent Protochlorophyllide Oxidoreductase in Thuja

The light-independent protochlorophyllide oxidoreductase (DPOR) plays a key role in the ability of nonflowering plants and algae to synthesize chlorophyll in darkness. This enzyme consists of three subunits encoded by the chlB, chlL, and chlN genes in the plastid genome. Previously, we found a high nonsynonymous substitution rate (dN) of the chlL gene in the lineage of Thuja standishii, a conifer belonging to the Cupressaceae. Here we revealed that the acceleration of dN in the chlL occurred as well in other species of Thuja, Thuja occidentalis and Thuja plicata. In addition, dark-grown seedlings of T. occidentalis were found to exhibit a pale yellowish color, and their chlorophyll concentration was much lower than that of other species of Cupressaceae. The results suggested that the species of Thuja have lost the ability to synthesize chlorophyll in darkness, and the functional constraint on the DPOR would thus be expected to be relaxed in this genus. Therefore, we expected to find that the evolutionary rates of all subunits of DPOR would in this case be accelerated. Sequence analyses of the chlN and chlB (encoding the other subunits of DPOR) in 18 species of Cupressaceae revealed that the dN of the chlN gene was accelerated in Thuja as was the dN of the chlL gene, but the dN of the chlB gene did not appear to differ significantly among the species of Cupressaceae. Sequencing of reverse transcription-polymerase chain reaction (RT-PCR) products of these genes showed that RNA editing was rare and unlikely to have contributed to the acceleration. Moreover, the RT-PCR analysis indicated that all chl genes were still transcriptionally active in T. occidentalis. Based on these results, it appears that species of Thuja still bear the DPOR protein, although the enzyme has lost its activity because of nonsynonymous mutations of some of the chl genes. The lack of acceleration of the dN of the chlB gene might be accounted for by various unknown functions of its gene product.

Contrasting patterns of DNA variation in natural populations of two related conifers, Cryptomeria japonica and Taxodium distichum (Cupressaceae sensu lato)

We investigated DNA variation within and between two closely related conifers, Cryptomeria japonica and Taxodium distichum, at nuclear loci encoding ferredoxin, glutamyl-tRNA reductase, lycopene beta cyclase, and phosphoribosylanthranilate transferase. Average nucleotide diversity at silent sites was estimated to be 0.0035 (SE 0.0012) in C. japonica and 0.0058 (SE 0.0006) in T. distichum. One population in C. japonica was differentiated from the others but generally there was not much differentiation among populations or varieties within the two species. However, the two species seemed to differ in frequency spectra of DNA polymorphisms. Excesses of intermediate-frequency variants were found in C. japonica, whereas excesses of both rare and high-frequency variants were found in T. distichum, which suggested different histories of population structures in the two species. Deviations from the standard neutral expectations in DNA polymorphisms were found by applications of neutrality tests. The results show that actions of selection to respective loci seem to differ between the two species, indicating differences of interaction among evolutionary factors.

Comparative analysis of expressed sequence tags of conifers and angiosperms reveals sequences specifically conserved in conifers

To identify and characterize lineage-specific genes of conifers, two sets of ESTs (with 12791 and 5902 ESTs, representing 5373 and 3018 gene transcripts, respectively) were generated from the Cupressaceae species Cryptomeria japonica and Chamaecyparis obtusa. These transcripts were compared with non-redundant sets of genes generated from Pinaceae species, other gymnosperms and angiosperms. About 6% of tentative unique genes (Unigenes) of C. japonica and C. obtusa had homologs in other conifers but not angiosperms, and about 70% had apparent homologs in angiosperms. The calculated GC contents of orthologous genes showed that GC contents of coniferous genes are likely to be lower than those of angiosperms. Comparisons of the numbers of homologous genes in each species suggest that copy numbers of genes may be correlated between diverse seed plants. This correlation suggests that the multiplicity of such genes may have arisen before the divergence of gymnosperms and angiosperms.

The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis

Chloroplast DNA sequences are a primary source of data for plant molecular systematic studies. A few key papers have provided the molecular systematics community with universal primer pairs for noncoding regions that have dominated the field, namely trnL-trnF and trnK/matK. These two regions have provided adequate information to resolve species relationships in some taxa, but often provide little resolution at low taxonomic levels. To obtain better phylogenetic resolution, sequence data from these regions are often coupled with other sequence data. Choosing an appropriate cpDNA region for phylogenetic investigation is difficult because of the scarcity of information about the tempo of evolutionary rates among different noncoding cpDNA regions. The focus of this investigation was to determine whether there is any predictable rate heterogeneity among 21 noncoding cpDNA regions identified as phylogenetically useful at low levels. To test for rate heterogeneity among the different cpDNA regions, we used three species from each of 10 groups representing eight major phylogenetic lineages of phanerogams. The results of this study clearly show that a survey using as few as three representative taxa can be predictive of the amount of phylogenetic information offered by a cpDNA region and that rate heterogeneity exists among noncoding cpDNA regions.

Evaluation of cleaved amplified polymorphic sequence markers for Chamaecyparis obtusa based on expressed sequence tag information from Cryptomeria japonica

We have developed and evaluated sequence-tagged site (STS) primers based on expressed sequence-tag information derived from sugi (Cryptomeria japonica) for use in hinoki (Chamaecyparis obtusa), a species that belongs to a different family (although it appears to be fairly closely related to sugi). Of the 417 C. japonica STS primer pairs we screened, 120 (approximately 30%) were transferable and provided specific PCR amplification products from 16 C. obtusa plus trees. We used haploid megagametophytes to investigate the homology of 80 STS fragments between C. obtusa and C. japonica and to identify orthologous loci. Nearly 90% of the fragments showed high (>70%) degrees of similarity between the species, and 35 STSs indicated homology to entries with the same putative function in a public DNA database. Of the 120 STS fragments amplified, 72 showed restriction fragment length polymorphisms; in addition, the CC2430 primers detected amplicon length polymorphism. We assessed the inheritance pattern of 27 cleaved amplified polymorphic sequence markers, using 20 individuals from the segregation population. All the markers analyzed were consistent with the marker inheritance patterns obtained from the screening panel, and no markers (except CC2716) showed significant (P<0.01) deviation from the expected segregation ratio. In total, 136 polymorphic markers were developed using C. japonica-based STS primers without any sequence modification. In addition, the applicability of STS-based markers developed in one species to other species was found to closely reflect the evolutionary distance between the species, which is roughly concordant with the difference between their rbcL sequences. We plan to use these markers for genetic studies in C. obtusa. Most of the markers should also provide reliable anchor loci for comparative mapping studies of the C. obtusa and C. japonica genomes.

The circumscription and phylogenetic relationships of Callitropsis and the newly described genus Xanthocyparis (Cupressaceae)

A new species of conifer was recently discovered in northern Vietnam. In a preliminary phylogenetic analysis of morphological data a possible sister species, Chamaecyparis nootkatensis (D. Don) Spach, was identified; however, because of the presumed phylogenetic remoteness of these two species to the remainder of the Cupressaceae, a new genus-Xanthocyparis-was described to accommodate both species. Here an analysis of ITS (nrDNA), matK, and rbcL sequence data in combination with 58 informative morphological characters was aimed at testing the monophyly of the remainder of Chamaecyparis and evaluating the placement and monophyly of Xanthocyparis. Chamaecyparis, minus C. nootkatensis, was resolved as a monophyletic group, remote from Cupressus and Xanthocyparis. Cupressus, Juniperus, and Xanthocyparis formed a very highly supported monophyletic group. However, Cupressus was not monophyletic. Instead the Old World species sampled were resolved sister to a clade containing a monophyletic Juniperus, a monophyletic Xanthocyparis, and a clade of New World Cupressus species. If both species of Xanthocyparis are to be treated as members of the same genus, then due to the principal of priority they will have to be recognized in the genus Callitropsis. Research is continuing to resolve the status of New World and Old World Cupressus.

Addressing the "hardest puzzle in American pomology:" Phylogeny of Prunus sect. Prunocerasus (Rosaceae) based on seven noncoding chloroplast DNA regions

Prunus subg. Prunus sect. Prunocerasus (Rosaceae) is a North American taxon with 17 commonly recognized taxa. To test the hypothesis of monophyly for the section we sequenced the trnG and rpL16 introns and the trnH-psbA and trnS-trnG intergenic spacers for at least two representatives of each of the five subgenera in Prunus. Additionally we sampled heavily among Prunus subg. Prunus sections Prunus and Armeniaca and Prunus subg. Amygdalus because these groups are putatively most closely related to Prunocerasus. Once monophyly of sect. Prunocerasus was shown we added the sequences of trnL and rpS16 introns and the trnL-trnF spacer in an attempt to increase resolution within the section. The species of sect. Prunocerasus showed an initial split with P. subcordata, the only species from western North America, sister to the rest of the group. The remaining species fell into three primary clades. Within each of the three primary clades there was little phylogenetic resolution. Lastly, we present evidence that P. texana, previously classified in subg. Amygdalus, may be a plum or at least contain a Prunocerasus chloroplast. This is the first phylogenetic hypothesis presented for sect. Prunocerasus, and the clades recovered contrast sharply with previously defined groups based on morphological characters.

Distribution of calcium oxalate crystals in the secondary phloem of conifers: a constitutive defense mechanism?

•   Calcium oxalate (CaOx) crystals were studied in the secondary phloem of 46 conifer species from six families to characterize their distribution relative to defense. •   Scanning electron microscopy, and polarized light microscopy of stem tissue allowed analysis of the quantity, location, morphology and size of the crystals. •   Vastly different patterns of CaOx crystal deposition were observed in stems of Pinaceae and nonPinaceae lineages. The CaOx crystals were present in all species but the highest density occurred along the nonPinaceae lineage. In Pinaceae, all species accumulated crystals intracellularly in crystalliferous parenchyma, whereas all nonPinaceae species had only extracellular crystals. A possible relationship between the number of aggressive bark beetles species and the amount of CaOx accumulation was noted where increased crystal accumulation appears to be antagonistic to beetle attack. Mapping along with this trend was the presence of phloem fibers. •   We conclude that in conifer stems the patterns and frequency of CaOx crystals function as a constitutive defense and in combination with fiber rows, provides an effective barrier against small bark-boring insects.

DNA variation in a conifer, Cryptomeria japonica (Cupressaceae sensu lato)

We investigated the nucleotide variation of a conifer, Cryptomeria japonica, and the divergence between this species and its closest relative, Taxodium distichum, at seven nuclear loci (Acl5, Chi1, Ferr, GapC, HemA, Lcyb, and Pat). Samples of C. japonica were collected from three areas, Kantou-Toukai, Hokuriku, and Iwate. No apparent geographic differentiation was found among these samples. However, the frequency spectrum of the nucleotide polymorphism revealed excesses of intermediate-frequency variants, which suggests that the population was not panmictic and a constant size in the past. The average nucleotide diversity, pi, for silent sites was 0.00383. However, values of pi for silent sites vary among loci. Comparisons of polymorphism to divergence among loci (the HKA test) showed that the polymorphism at the Acl5 locus was significantly lower. We also observed a nearly significant excess of replacement polymorphisms at the Lcyb locus. These results suggested possibilities of natural selection acting at some of the loci. Intragenic recombination was detected only once at the Chi1 locus and was not detected at the other loci. The low level of population recombination rate, 4Nr, seemed to be due to both low level of recombination, r, and small population size, N.

Isolation, characterization, and cross-species utility of microsatellites in yellow cedar (Chamaecyparis nootkatensis)

Chamaecyparis nootkatensis is an ecologically and economically important conifer of the north Pacific coastal forests. To aid in studies of clonal structure and genetic differentiation of this and related species, we isolated and characterized microsatellites from C. nootkatensis. A microsatellite-enriched library yielded 75 repeat-containing sequences for which primer pairs were designed. Only five showed reliable amplification and polymorphism, with an average of 13.7 alleles/locus and a mean expected heterozygosity of 0.592. In progeny tests with four families, few null alleles were directly detected and loci segregated according to Mendelian expectations. However, in one primer pair, high heterozygote deficiency was observed, suggesting the presence of a null allele. The ability of primer pairs to cross amplify was tested on 18 species of the Cupressaceae sensu lato; three primer pairs yielded polymorphic loci in Cupressus and Juniperus species, but not in other Chamaecyparis species. This also supports recent findings of a closer affinity of C. nootkatensis with Cupressus over other Chamaecyparis species.

Molecular evolution of nuclear genes in Cupressacea, a group of conifer trees

We surveyed the molecular evolutionary characteristics of 11 nuclear genes from 10 conifer trees belonging to the Taxodioideae, the Cupressoideae, and the Sequoioideae. Comparisons of substitution rates among the lineages indicated that the synonymous substitution rates of the Cupressoideae lineage were higher than those of the Taxodioideae. This result parallels the pattern previously found in plastid genes. Likelihood-ratio tests showed that the nonsynonymous-synonymous rate ratio did not change significantly among lineages. In addition, after adjustments for lineage effects, the dispersion indices of synonymous and nonsynonymous substitutions were considerably reduced, and the latter was close to 1. These results indicated that the acceleration of evolutionary rates in the Cupressoideae lineage occurred in both the nuclear and plastid genomes, and that generally, this lineage effect affected synonymous and nonsynonymous substitutions similarly. We also investigated the relationship of synonymous substitution rates with the nonsynonymous substitution rate, base composition, and codon bias in each lineage. Synonymous substitution rates were positively correlated with nonsynonymous substitution rates and GC content at third codon positions, but synonymous substitution rates were not correlated with codon bias. Finally, we tested the possibility of positive selection at the protein level, using maximum likelihood models, assuming heterogeneous nonsynonymous-synonymous rate ratios among codon (amino acid) sites. Although we did not detect strong evidence of positively selected codon sites, the analysis suggested that significant variation in nonsynonymous-synonymous rate ratio exists among the sites. The most likely sites for action of positive selection were found in the ferredoxin gene, which is an important component of the apparatus for photosynthesis.

Nucleotide sequence variation of the chloroplast trnK/matK region in two wild Fagopyrum(Polygonaceae) species, F. leptopodum and F. statice

Nucleotide sequence polymorphisms of the intron of the chloroplast trnK (UUU) gene, including a matK gene, were investigated within two wild Fagopyrum species, F. leptopodum and F. statice, to assess the degree and pattern of the inter- and intraspecific differences in coding and noncoding chloroplast DNA regions in higher plants. Ten and five accessions were used for F. leptopodum and F. statice, respectively. The length of the trnK intron region in these species ranged from 2494 to 2506 bp. In the trnK intron, the net nucleotide substitution number per site (Da) between the two species was 0.00109, lower than the nucleotide diversity (pi), 0.00195 for F. leptopodum and 0.00144 for F. statice, suggesting a low level of interspecific divergence. This result seems to be due to the phylogenetic pattern that both species are interspersed with each other, which was revealed by the phylogenetic analyses based on the nucleotide substitutions and indels. In the matK gene region (1524 bp), seven and two nucleotide substitutions were found within F. leptopodum and F. statice, respectively. All of the nine nucleotide substitutions (eight of which were nonsynonymous) within and between F. leptopodum and F. statice were clustered in the 5' part of the matK gene region, and no variation was found in the 3' part. This suggests that most of the 3' part is occupied by the conserved domains that are important for the binding activity of the gene product to the precursor mRNA, and therefore implies that the 3' part is more functionally constrained than the 5' part.