A phylogenetic and biogeographic study of Hamamelis (Hamamelidaceae), an eastern Asian and eastern North American disjunct genus

The chemical constituents from twigs of Hamamelis japonica and their antiviral activities

Three new monoterpenoid glycosides (1-3) and one new flavanol (4) along with 15 known compounds were isolated from the twig of Hamamelis japonica Sieb. et Zucc. The chemical constituent study of the twig of H. japonica has performed for the first time in the present investigation. Their structures were determined based on extensive spectroscopic methods including 1 D and 2 D NMR and CD spectra data. All isolated compounds were tested for their antiviral activities against HRV1B-, EV71-, PR8- and CVB3-infected Vero cells. Among the tested compounds, (-)-epigallocatechin 3-O-gallate exhibited the most consistent and effective antiviral activities against EV71 and PR8 infections.

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.-Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single-copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single-copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132-133 genes, including 85-87 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37-43 long repeats and 175-178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.

A New Pipeline for Removing Paralogs in Target Enrichment Data

Target enrichment (such as Hyb-Seq) is a well-established high throughput sequencing method that has been increasingly used for phylogenomic studies. Unfortunately, current widely used pipelines for analysis of target enrichment data do not have a vigorous procedure to remove paralogs in target enrichment data. In this study, we develop a pipeline we call Putative Paralogs Detection (PPD) to better address putative paralogs from enrichment data. The new pipeline is an add-on to the existing HybPiper pipeline, and the entire pipeline applies criteria in both sequence similarity and heterozygous sites at each locus in the identification of paralogs. Users may adjust the thresholds of sequence identity and heterozygous sites to identify and remove paralogs according to the level of phylogenetic divergence of their group of interest. The new pipeline also removes highly polymorphic sites attributed to errors in sequence assembly and gappy regions in the alignment. We demonstrated the value of the new pipeline using empirical data generated from Hyb-Seq and the Angiosperms353 kit for two woody genera Castanea (Fagaceae, Fagales) and Hamamelis (Hamamelidaceae, Saxifragales). Comparisons of data sets showed that the PPD identified many more putative paralogs than the popular method HybPiper. Comparisons of tree topologies and divergence times showed evident differences between data from HybPiper and data from our new PPD pipeline. We further evaluated the accuracy and error rates of PPD by BLAST mapping of putative paralogous and orthologous sequences to a reference genome sequence of Castanea mollissima. Compared to HybPiper alone, PPD identified substantially more paralogous gene sequences that mapped to multiple regions of the reference genome (31 genes for PPD compared with 4 genes for HybPiper alone). In conjunction with HybPiper, paralogous genes identified by both pipelines can be removed resulting in the construction of more robust orthologous gene data sets for phylogenomic and divergence time analyses. Our study demonstrates the value of Hyb-Seq with data derived from the Angiosperms353 probe set for elucidating species relationships within a genus, and argues for the importance of additional steps to filter paralogous genes and poorly aligned regions (e.g., as occur through assembly errors), such as our new PPD pipeline described in this study. [Angiosperms353; Castanea; divergence time; Hamamelis; Hyb-Seq, paralogs, phylogenomics.]

Genetic insights into the evolution of genera with the eastern Asia-eastern North America floristic disjunction: a transcriptomics analysis

PREMISE

Large disjunctions in species distributions provide excellent opportunities to study processes that shape biogeographic patterns. One such disjunction is the eastern Asia-eastern North America (EA-ENA) floristic disjunction. For many genera with this disjunction, species richness is greater in EA than in ENA; this pattern has been attributed, in part, to higher rates of molecular evolution and speciation in EA. Longer branch lengths have been found in some EA clades, relative to their ENA sister clades, suggesting that the EA lineages have evolved at a higher rate, possibly due to environmental heterogeneity, potentially contributing to the species richness anomaly.

METHODS

To evaluate whether rates of molecular evolution are elevated in EA relative to ENA, we used transcriptomes from species in 11 genera displaying this disjunction. Rates of molecular evolution were estimated for up to 385 orthologous nuclear loci per genus.

RESULTS

No statistically significant differences were identified in pairwise comparisons between EA and ENA sister species, suggesting equal rates of molecular evolution for both species; the data also suggest similar selection pressures in both regions. For larger genera, evidence likewise argues against more species-rich clades having higher molecular evolutionary rates, regardless of region. Our results suggest that genes across multiple gene ontology categories are evolving at similar rates under purifying selection in species in both regions.

CONCLUSIONS

Our data support the hypothesis that greater species richness in EA than ENA is due to factors other than an overall increase in rates of molecular evolution in EA.

Relationships between Tertiary relict and circumboreal woodland floras: a case study in Chimaphila (Ericaceae)

BACKGROUND AND AIMS

Tertiary relict and Arctic/circumboreal distributions are two major patterns of Northern Hemisphere intercontinental disjunctions with very different histories. Each has been well researched, but members of one biome have generally not been incorporated in the biogeographical analyses of the other, and links or transitions between these two biomes have rarely been addressed.

METHODS

Phylogenies of Chimaphila were generated based on cpDNA and nuclear ITS, using Bayesian and maximum likelihood methods. A time-calibrated phylogeny was generated using BEAST. Ancestral area reconstruction was inferred using both statistical dispersal-vicariance analysis and a dispersal-extinction-cladogenesis model.

KEY RESULTS

The Chimaphila crown group was estimated to have originated in the early Miocene. The lineages of C. umbellata diverged early, but its present circumboreal distribution was not achieved until around the middle Pliocene or later. Sister to this is a clade of four species with Tertiary relict distribution. Among these, two expansions occurred from North America to Asia, probably via the Bering Land Bridge, generating its current disjunctions.

CONCLUSIONS

Our data concur with a few other studies, indicating that the circumboreal woodland biome has an older origin than most true Arctic-alpine taxa, having gradually recruited taxa since the early Oligocene. For the origin of Asia-North America disjunctions in Chimaphila, an 'out-of-America' migration was supported. It is not clear in which direction Pyroloideae lineages moved between Tertiary relict disjunctions and Arctic/circumboreal distributions; each biome might have recruited species from the other.

Molecular Phylogeography and Ecological Niche Modeling of Sibbaldia procumbens s.l. (Rosaceae)

The phylogeographical analysis and ecological niche modeling (ENM) of the widely distributed Northern Hemisphere Sibbaldia procumbens s.l. can help evaluate how tectonic motion and climate change helped shape the current distribution patterns of this species. Three chloroplast regions (the atpI-atpH and trnL-trnF intergenic spacers and the trnL intron) were obtained from 332 (156 from present study and 176 from the previous study) individuals of S. procumbens s.l. An unrooted haplotype network was constructed using the software NETWORK, while BEAST was used to estimate the divergence times among haplotypes. ENM was performed by MAXENT to explore the historical dynamic distribution of S. procumbens s.l. The haplotype distribution demonstrates significant phylogeographical structure (N_ST > G_ST; P < 0.01). The best partitioning of genetic diversity by SAMOVA produced three groups, while the time to the most recent common ancestor of all haplotypes was estimated to originate during the Miocene, with most of the haplotype diversity having occurred during the Quaternary. The MAXENT analysis showed S. procumbens s.l. had a wider distribution range during the last glacial maximum and a narrower distribution range during the last interglacial, with predictions into the future showing the distribution range of S. procumbens s.l. shrinking.

Identification and development of microsatellite markers in Hamamelis mollis (Hamamelidaceae)

PREMISE OF THE STUDY

Hamamelis mollis (Hamamelidaceae) is a Tertiary relict species endemic to southern China. Polymorphic microsatellite markers were developed to reveal the genetic diversity of this species.

METHODS AND RESULTS

The genome of H. mollis was sequenced and de novo assembled into 642,351 contigs. A total of 72,097 paired primers were successfully designed from 80,282 simple sequence repeat (SSR) markers identified in 63,419 contigs. PCR amplification showed that 96 of the 136 synthesized primers could be successfully amplified, and 22 demonstrated polymorphism. The mean number of alleles, levels of observed heterozygosity, and levels of expected heterozygosity were 4.602 ± 0.140, 0.632 ± 0.020, and 0.696 ± 0.010, respectively. The majority of the 96 primer pairs could be amplified in at least one other Hamamelidaceae species, including Distylium myricoides (60), Loropetalum chinense (39), Exbucklandia populnea (24), and E. tonkinensis (24).

CONCLUSIONS

These microsatellite loci provide abundant genomic SSR markers to evaluate genetic diversity of this woody ornamental plant.

Evolution of biogeographic disjunction between eastern Asia and North America in Chamaecyparis: Insights from ecological niche models

The disjunct distribution of plants between eastern Asia (EA) and North America (NA) is one of the most well-known biogeographic patterns. However, the formation and historical process of this pattern have been long debated. Chamaecyparis is a good model to test previous hypotheses about the formation of this disjunct pattern as it contains six species disjunctly distributed in EA, western North America (WNA) and eastern North America (ENA). In this study, we applied ecological niche models to test the formation of the disjunct pattern of Chamaecyparis. The model calibrated with the EA species was able to predict the distribution of eastern NA species well, but not the western NA species. Furthermore, the eastern Asian species were shown to have higher niche overlap with the eastern North American species. The EA species were also shown to share more similar habitats with ENA species than with WNA species in the genus. Chamaecyparis species in WNA experienced a significant niche shift compared with congeneric species. Chamaecyparis had a low number of suitable regions in Europe and the middle and western NA during the Last Glacial Maximum (LGM) period, and became extinct in the former region whereas it retains residual distribution in the latter. The extirpations in western NA and Europe in response to the late Neogene and Quaternary climatic cooling and the more similar habitats between ENA and EA ultimately shaped the current intercontinental disjunct distribution of Chamaecyparis. Both current hypotheses may be also jointly applied to explain more eastern Asian and eastern North American disjunctions observed today.

Is Oligomeris (Resedaceae) indigenous to North America? Molecular evidence for a natural colonization from the Old World

Oligomeris linifolia constitutes one of the few examples of intercontinental disjunctions at the species level between the arid regions of the Old World and SW North America. The status of the American populations has been obscure, with some authors considering the populations to be introduced, whereas others believe them to be native. To clarify these conflicting opinions, we performed phylogeographic analyses using nuclear ribosomal ITS and plastid trnL-F and rps16 sequences to infer the origin of the disjunct American populations. Two independent molecular clock approaches based on ITS and cpDNA sequences (rbcL, matK, trnL-F) were used to estimate a divergence time of O. linifolia. Low levels of sequence divergence and estimates of relatively recent splits of Oligomeris lineages disagree with the vicariance hypotheses traditionally suggested to account for New-Old World disjunctions. In addition, significant genetic differentiation of American populations does not indicate a recent anthropogenic introduction. Morphological uniformity and the sharing of haplotypes between disjunct populations, together with the molecular clock results, suggest that a long-distance dispersal event from the Old Word to SW North America may have taken place during the Quaternary, in spite of limited dispersal mechanisms in Oligomeris.

Extreme long-distance dispersal of the lowland tropical rainforest tree Ceiba pentandra L. (Malvaceae) in Africa and the Neotropics

Many tropical tree species occupy continental expanses of rainforest and flank dispersal barriers such as oceans and mountains. The role of long-distance dispersal in establishing the range of such species is poorly understood. In this study, we test vicariance hypotheses for range disjunctions in the rainforest tree Ceiba pentandra, which is naturally widespread across equatorial Africa and the Neotropics. Approximate molecular clocks were applied to nuclear ribosomal [ITS (internal transcribed spacer)] and chloroplast (psbB-psbF) spacer DNA sampled from 12 Neotropical and five West African populations. The ITS (N=5) and psbB-psbF (N=2) haplotypes exhibited few nucleotide differences, and ITS and psbB-psbF haplotypes were shared by populations on both continents. The low levels of nucleotide divergence falsify vicariance explanations for transatlantic and cross-Andean range disjunctions. The study shows how extreme long-distance dispersal, via wind or marine currents, creates taxonomic similarities in the plant communities of Africa and the Neotropics.

Regional differences in rates of plant speciation and molecular evolution: a comparison between eastern Asia and eastern North America

The eastern Asian (EAS)-eastern North American (ENA) floristic disjunction is one of the best-known biogeographic patterns in the Northern Hemisphere. Recent paleontological and molecular analyses have illuminated the origins of the biogeographic pattern, but subsequent diversification and evolution of the disjunct floras in each of the two continents after isolation remains poorly understood. Although similar in climate and floristic composition, EAS has twice as many species as ENA in genera occurring in both regions. Explaining such differences in species diversity between regions with similar environmental conditions (diversity anomalies) is an important goal of the study of the global patterns of biodiversity. We used a phylogenetic approach to compare rates of net speciation and molecular evolution between the two regions. We first identified EAS-ENA disjunct sister clades from ten genera (Asarum, Buckleya, Carpinus, Carya, Cornus, Hamamelis, Illicium, Panax, Stewartia, and Styrax) that represent diverse angiosperm lineages using phylogenetic analyses of ITS (internal transcribed spacer of nuclear ribosomal DNA) sequence data. Species richness and substitution rate of ITS between sister clades were compared. The results revealed a pattern of greater species diversity in the EAS counterparts. A positive relationship between species diversity and ITS substitution rate was also documented. These results suggest greater net speciation and accelerated molecular evolution in EAS. The data support the idea that a regional difference in net speciation rate related to topographic heterogeneity contributes to the diversity anomaly between EAS and ENA. The close relationship between rates of ITS evolution and species richness further suggests that species production may be directly linked to rate of nucleotide substitution.

Anatomically preserved Liquidambar (Altingiaceae) from the middle Miocene of Yakima Canyon, Washington state, USA, and its biogeographic implications

Liquidambar changii Pigg, Ickert-Bond & Wen sp. nov. (Altingiaceae) is established for anatomically preserved, middle Miocene infructescences from Yakima Canyon, Washington, USA. Specimens are spherical, ∼2.5 cm in diameter, and have ∼25-30 tightly packed, bilocular fruits per head. Fruits are 3.4-4.7 mm wide × 2.6-3.5 mm long and wedge shaped, fused at the base, and free distally. Each locule contains 1-2 mature, elongate seeds proximally and 5-9 aborted seeds of more irregular shape distally. Mature seeds are 1.5 mm long × 1.2 mm wide, elongate, and triangular transversely, with a slight flange. Seeds have a seed coat for which three zones can be well defined, a uniseriate outer palisade layer, a middle region of isodiametric cells comprising most of the integument, and a uniseriate inner layer of tangentially elongate cells lining the embryo cavity. Liquidambar changii is most similar to the eastern Asian L. acalycina H.-T. Chang on features of infructescence, fruit, and seed morphology and quite unlike the North American L. styraciflua L. and other species. Such a close relationship between these two species supports a Beringian biogeographic track between eastern Asia and western North America during the Miocene. Previous phylogenetic and allozyme analysis of modern Liquidambar demonstrates a close relationship between North American-western Asian taxa and suggests a North Atlantic biogeographic track in the middle Miocene. Together, these biogeographic tracks underscore the complexity of the biogeographic history of the Altingiaceae in the Northern Hemisphere throughout the Neogene.

Phylogenetic relationships in Gleditsia (Leguminosae) based on ITS sequences

We used nucleotide sequences from the internal transcribed spacers and 5.8S gene of nuclear ribosomal DNA to test competing phylogenetic and biogeographic hypotheses in Gleditsia. Eleven of 13 Gleditsia species were sampled, along with two species of its sister genus, Gymnocladus. Analyses of ITS data and of a combined data set that included sequences of ITS and two chloroplast genes supported several conclusions that were interpreted in light of fossil data and current legume phylogeny. Gleditsia and Gymnocladus appear to have originated in eastern Asia during the Eocene. Eastern North American species of both genera most likely evolved from ancestors that migrated across the Bering land bridge, but the eastern Asian/eastern North American disjunction appears to be much older in Gymnocladus than in Gleditsia. Gleditsia amorphoides, from temperate South America, is sister to the rest of the genus, suggesting early long-distance dispersal from Asia. The remainder of Gleditsia is divided into three unresolved clades, possibly indicating a split early in the evolution of the genus. Two of those clades contain only Asian species, and one contains Asian and North American species. The North American species, Gleditsia triacanthos and Gleditsia aquatica, are polymorphic and paraphyletic with respect to their ITS and cpDNA sequences, which suggests recent diversification.

Phylogenetics and evolution of the eastern Asian-eastern North American disjunct aphid tribe, Hormaphidini (Hemiptera: Aphididae)

A conspicuous biogeographic pattern of the Northern Hemisphere is the disjunct occurrence of related taxa on different continents. Perhaps best studied in plants, this pattern includes disjunct distributions of genera in eastern Asia and eastern North America. Such continental disjunctions are thought to be the remnants of a mostly continuously distributed, mixed mesophytic forest dating to the Miocene, which subsequently became fragmented by geological and climatic changes. Some highly host-specific insects, namely aphids, live on descendants of the mixed mesophytic forest taxa and exhibit the same disjunct distributions as that of their host plants. We estimated the phylogeny of Hormaphidini aphids, which host-alternate between witch-hazel (Hamamelis; an eastern Asian-eastern North American disjunct genus) and birch (Betula). Based on partial nuclear elongation factor 1alpha and mitochondrial tRNA leucine/cytochrome oxidase II sequences, trees inferred from maximum-parsimony and maximum-likelihood showed strong support for two monophyletic genera (Hamamelistes and Hormaphis), each containing a clade of Japanese and a clade of North American species. The estimated divergence dates of Asian and North American clades in both genera was 20-30 million years ago, consistent with the idea that aphids may have experienced the same vicariance events as those of their host plants.

Many to flower, few to fruit: the reproductive biology of Hamamelis virginiana (Hamamelidaceae)

Hamamelis virginiana flowers from late September to late November. In 1977, we began studying the reproductive biology of this eastern North American arborescent shrub by examining floral phenology and rewards, pollen-ovule ratios, breeding system, pollination, pollinator and resource limitation, and seed dispersal. The homogamous, self-incompatible flowers emit a faint odor, bear nectar with sucrose ratios typical of bee- and fly-pollinated flowers, and produce abundant sticky pollen. Flowers were visited infrequently by insects representing six orders. Flies were the most common floral visitors, specifically members of the genus Bradysia, but small bees also carried high percentages of Hamamelis pollen. Despite high pollen/ovule ratios (11 445 grains/ovule), bees and flies are likely pollinators, as experiments indicate wind pollination is less likely. Pollen quantity and resource availability did not appear to limit reproductive output, but pollen quality did. Tests of >40 000 flowers showed natural fruit set to be <1%. The flowering time, breeding system, and clumped distribution of plants, likely due in part to limited seed dispersal, combine to yield this remarkably low fruit set. Because all other species of Hamamelis flower from late winter to early summer, it may be that H. virginiana evolved a fall flowering phenology to avoid competition for pollinators with the closely related H. vernalis.

Molecular phylogeny of the Magnoliaceae: the biogeography of tropical and temperate disjunctions

The boreotropical flora concept suggests that relictual tropical disjunctions between Asia and the Americas are a result of the expansion of the circumboreal tropical flora from the middle to the close of the Eocene. Subsequently, temperate species diverged at high latitudes and migrated to other continents. To test this concept, we conducted a molecular phylogenetic analysis (using cpDNA) of the Magnoliaceae, a former boreotropical element that currently contains both tropical and temperate disjuncts. Divergence times of the clades were estimated using sequences of matK and two intergenic regions consisting of psbA-trnH and atpB-rbcL. Results indicate the tropical American section Talauma branched first, followed by the tropical Asian clade and the West Indies clade. Within the remaining taxa, two temperate disjunctions were formed. Assuming the temperate disjunction of Magnolia acuminata and Asian relatives occurred 25 mya (late Oligocene; based on seed fossil records), section Talauma diverged 42 mya (mid-Eocene), and tropical Asian and the West Indies clades 36 mya (late Eocene). These events correlate with cooling temperatures during the middle to late Eocene and probably caused the tropical disjunctions.