Multiple continental radiations and correlates of diversification in Lupinus (Leguminosae): testing for key innovation with incomplete taxon sampling

Replicate radiations provide powerful comparative systems to address questions about the interplay between opportunity and innovation in driving episodes of diversification and the factors limiting their subsequent progression. However, such systems have been rarely documented at intercontinental scales. Here, we evaluate the hypothesis of multiple radiations in the genus Lupinus (Leguminosae), which exhibits some of the highest known rates of net diversification in plants. Given that incomplete taxon sampling, background extinction, and lineage-specific variation in diversification rates can confound macroevolutionary inferences regarding the timing and mechanisms of cladogenesis, we used Bayesian relaxed clock phylogenetic analyses as well as MEDUSA and BiSSE birth-death likelihood models of diversification, to evaluate the evolutionary patterns of lineage accumulation in Lupinus. We identified 3 significant shifts to increased rates of net diversification (r) relative to background levels in the genus (r = 0.18-0.48 lineages/myr). The primary shift occurred approximately 4.6 Ma (r = 0.48-1.76) in the montane regions of western North America, followed by a secondary shift approximately 2.7 Ma (r = 0.89-3.33) associated with range expansion and diversification of allopatrically distributed sister clades in the Mexican highlands and Andes. We also recovered evidence for a third independent shift approximately 6.5 Ma at the base of a lower elevation eastern South American grassland and campo rupestre clade (r = 0.36-1.33). Bayesian ancestral state reconstructions and BiSSE likelihood analyses of correlated diversification indicated that increased rates of speciation are strongly associated with the derived evolution of perennial life history and invasion of montane ecosystems. Although we currently lack hard evidence for "replicate adaptive radiations" in the sense of convergent morphological and ecological trajectories among species in different clades, these results are consistent with the hypothesis that iteroparity functioned as an adaptive key innovation, providing a mechanism for range expansion and rapid divergence in upper elevation regions across much of the New World.     

The impact of fossils and taxon sampling on ancient molecular dating analyses

The number and complexity of molecular dating studies has increased over the past decade. Along with a broadening acceptance of their utility has come significant controversy over the methods and models that are appropriate, as well as the accuracy of the estimates yielded by molecular clock analyses. Radically different age estimates have been published for the same divergences from analyses of different datasets with different fossil constraints obtained with different methods, and the underlying explanation for these differences is often unclear. Here we utilize two previously published datasets to examine the effect of fossil calibrations and taxon sampling on the age estimates for two deep eukaryote divergences in an attempt to discern the relative impact of these factors. Penalized likelihood, non-parametric rate smoothing, and Bayesian methods were utilized to generate age estimates for the origin of the Metazoa from a 7-gene dataset and for the divergence of Eukaryotes from a 129-gene dataset. From these analyses, it is clear that the fossil calibrations chosen and the method for applying constraints to these nodes have a large impact on age estimates, while the degree of taxon sampling within a dataset is less important in terms of the resulting age estimates. Concerns and recommendations for addressing these two factors when initiating a dating analysis are discussed.     

Into and out of the tropics: global diversification patterns in a hyperdiverse clade of ectomycorrhizal fungi

Ectomycorrhizal (ECM) fungi, symbiotic mutualists of many dominant tree and shrub species, exhibit a biogeographic pattern counter to the established latitudinal diversity gradient of most macroflora and fauna. However, an evolutionary basis for this pattern has not been explicitly tested in a diverse lineage. In this study, we reconstructed a mega‐phylogeny of a cosmopolitan and hyperdiverse genus of ECM fungi, Russula, sampling from annotated collections and utilizing publically available sequences deposited in GenBank. Metadata from molecular operational taxonomic unit cluster sets were examined to infer the distribution and plant association of the genus. This allowed us to test for differences in patterns of diversification between tropical and extratropical taxa, as well as how their associations with different plant lineages may be a driver of diversification. Results show that Russula is most species‐rich at temperate latitudes and ancestral state reconstruction shows that the genus initially diversified in temperate areas. Migration into and out of the tropics characterizes the early evolution of the genus, and these transitions have been frequent since this time. We propose the ‘generalized diversification rate’ hypothesis to explain the reversed latitudinal diversity gradient pattern in Russula as we detect a higher net diversification rate in extratropical lineages. Patterns of diversification with plant associates support host switching and host expansion as driving diversification, with a higher diversification rate in lineages associated with Pinaceae and frequent transitions to association with angiosperms.     

Phylogeography of Australian and New Zealand spray zone spiders (Anyphaenidae: Amaurobioides): Moa's Ark loses a few more passengers

Amaurobioides are restricted to the spray zone of southern continents, where they live in small, isolated populations and hunt from silk retreats built in rock crevices. A Star BEAST species tree based on ITS1 nuclear and ND1 mitochondrial genes did not support the hypothesis that this unusual niche linked the evolutionary history of these spiders to geological events reshaping Gondwana into present‐day Australia and New Zealand. Instead, it showed that Amaurobioides reached Australia approximately 4.5 Mya and dispersed twice to New Zealand. Approximately 2.37 Mya, spiders from Tasmania colonized the Deep South of the South Island and, approximately 0.38 Mya, those from South Australia colonized more northern regions. Thus, the present study further limits the scope of the Moa's Ark hypothesis of vicariant New Zealand biogeography.     

Temporal patterns of diversification in Andean Eois, a species-rich clade of moths (Lepidoptera, Geometridae)

The timing of the origin of present day Neotropical animal diversity is still a matter of debate. For a long time, a preponderance of glacial (i.e. Pleistocene) radiations has been proposed. However, recent data from molecular clock studies indicate a preglacial origin for most of the examined taxa. We performed a fossil-calibrated molecular dating analysis of the genus Eois, which is a major component of one of the world's most diverse assemblages of herbivorous insects. We found that diversification of Eois took place in the Miocene following a pattern best explained by density-dependent diversification. A strong slowdown of diversification towards the present was detected. Diversification of Eois does overlap with increased Andean uplift and diversification of the most commonly used host plant genus Piper. These findings match the patterns found for the majority of Neotropical tetrapods and for three other unrelated, ecologically different lepidopteran genera.     

The evolutionary diversification of parrots supports a taxon pulse model with multiple trans-oceanic dispersal events and local radiations

Vicariance is thought to have played a major role in the evolution of modern parrots. However, as the relationships especially of the African taxa remained mostly unresolved, it has been difficult to draw firm conclusions about the roles of dispersal and vicariance. Our analyses using the broadest taxon sampling of old world parrots ever based on 3219 bp of three nuclear genes revealed well-resolved and congruent phylogenetic hypotheses. Agapornis of Africa and Madagascar was found to be the sister group to Loriculus of Australasia and Indo-Malayasia and together they clustered with the Australasian Loriinae, Cyclopsittacini and Melopsittacus. Poicephalus and Psittacus from mainland Africa formed the sister group of the Neotropical Arini and Coracopsis from Madagascar and adjacent islands may be the closest relative of Psittrichas from New Guinea. These biogeographic relationships are best explained by independent colonization of the African continent via trans-oceanic dispersal from Australasia and Antarctica in the Paleogene following what may have been vicariance events in the late Cretaceous and/or early Paleogene. Our data support a taxon pulse model for the diversification of parrots whereby trans-oceanic dispersal played a more important role than previously thought and was the prerequisite for range expansion into new continents.     

基于Rag1基因序列变异的鲤形目鱼类系统发育重建:采样策略对生命之树的影响

The phylogenetic relationships of species are fundamental to any biological investigation, including all evolutionary studies. Accurate inferences of sister group relationships provide the researcher with an historical framework within which the attributes or geographic origin of species (or supraspecific groups) evolved. Taken out of this phylogenetic context, interpretations of evolutionary processes or origins, geographic distributions, or speciation rates and mechanisms, are subject to nothing less than a biological experiment without controls. Cypriniformes is the most diverse clade of freshwater fishes with estimates of diversity of nearly 3,500 species. These fishes display an amazing array of morphological, ecological, behavioral, and geographic diversity and offer a tremendous opportunity to enhance our understanding of the biotic and abiotic factors associated with diversification and adaptation to environments. Given the nearly global distribution of these fishes, they serve as an important model group for a plethora of biological investigations, including indicator species for future climatic changes. The occurrence of the zebrafish, Danio rerio, in this order makes this clade a critical component in understanding and predicting the relationship between mutagenesis and phenotypic expressions in vertebrates, including humans. With the tremendous diversity in Cypriniformes, our understanding of their phylogenetic relationships has not proceeded at an acceptable rate, despite a plethora of morphological and more recent molecular studies. Most studies are pre-Hennigian in origin or include relatively small numbers of taxa. Given that analyses of small numbers of taxa for molecular characters can be compromised by peculiarities of long-branch attraction and nodal-density effect, it is critical that significant progress in our understanding of the relationships of these important fishes occurs with increasing sampling of species to mitigate these potential problems. The recent Cypriniformes Tree of Life initiative is an effort to achieve this goal with morphological and molecular (mitochondrial and nuclear) data. In this early synthesis of our understanding of the phylogenetic relationships of these fishes, all types of data have contributed historically to improving our understanding, but not all analyses are complementary in taxon sampling, thus precluding direct understanding of the impact of taxon sampling on achieving accurate phylogenetic inferences. However, recent molecular studies do provide some insight and in some instances taxon sampling can be implicated as a variable that can influence sister group relationships. Other instances may also exist but without inclusion of more taxa for both mitochondrial and nuclear genes, one cannot distinguish between inferences being dictated by taxon sampling or the origins of the molecular data.     

Parentage analysis with incomplete sampling of candidate parents and offspring

Many breeding systems include 'multiple mating' in which males or females mate with multiple partners. We identify two forms of multiple mating: 'single-sex', where the next-generation individuals (NGIs) are the product of multiple mating by one sex; and 'two-sex', where the NGIs are the product of multiple mating by both sexes. For both mating systems we develop models that estimate the proportion of NGIs that is fathered (paternity) or mothered (maternity) by the putative parents. The models only require genetic data from the parent or parents in question and the sample of NGIs, as well as an estimate of population allele frequencies. The models provide unbiased estimates, can accommodate loci with many alleles and are robust to violations of their assumptions. They allow researchers to address intractable problems such as the parentage of seeds found on the ground, juvenile fish in a stream, and nestlings in a communal breeding bird. We demonstrate the models using genetic data from a nest of the bluegill sunfish Lepomis macrochirus, where the NGIs may be from multiple females that have spawned with multiple males from different life histories (cuckolder and parental).     

Ecological diversification and phylogeny of emydid turtles

Ecological diversification is a central topic in ecology and evolutionary biology. We undertook the first comprehensive species-level phylogenetic analysis of Emydidae (an ecologically diverse group of turtles), and used the resulting phylogeny to test four general hypotheses about ecological diversification. Phylogenetic analyses were based on data from morphology (237 parsimony-informative characters) and mitochondrial DNA sequences (547 parsimony-informative characters) and included 39 of the 40 currently recognized emydid species. Combined analyses of all data provide a well-supported hypothesis for intergeneric relationships, and support monophyly of the two subfamilies (Emydinae and Deirochelyinae) and most genera (with the notable exception of Clemmys and Trachemys ). Habitat and diet were mapped onto the combined-data tree to test fundamental hypotheses about ecological diversification. Using continuous coding of ecological characters showed that lineages changed in habitat before diet, ecological change was most frequently from generalist to specialist, and habitat and diet rarely changed on the same branch of the phylogeny. However, we also demonstrate that the results of ancestral trait reconstructions can be highly sensitive to character coding method (i.e. continuous vs. discrete). Finally, we propose a simple model to describe the pattern of ecological diversification in emydid turtles and other lineages, which may reconcile the (seemingly) conflicting conclusions of our study and two recent reviews of ecological diversification.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 79, 577–610.     

Molecular phylogeny and ecological diversification in a clade of New World songbirds (genus Vireo)

We constructed a molecular phylogeny for a clade of eye-ringed vireos ( Vireo flavifrons and the V . solitarius complex) to examine existing hypotheses of speciation and ecological diversification. Complete sequences of the mtDNA cytochrome b gene were obtained from 47 individuals of this group plus four vireonid outgroups. Mean levels of sequence divergence in the clade varied from 0.29% to 5.7%. Differences were greatest between V . flavifrons and four taxa of ' V . solitarius '. The latter separated into three taxonomic, geographical and ecological groups: V . plumbeus plumbeus, V . cassinii cassinii , and V . solitarius solitarius plus V . solitarius alticola . These differed by an average of 2.6–3.2%. Populations within each group revealed low levels of sequence variation ( x- = 0.20%) and little geographical structuring. The mtDNA data generally corroborate results from allozymes. V. plumbeus shows a loss of yellow-green carotenoid pigmentation from the ancestral condition. The occupancy of relatively dry habitats by this species and V . cassinii represents a derived ecological shift from more-humid environments occupied by other species of vireonids. Ecological divergence in this clade occurred in allopatry and is associated with generic-level stability in morphometrics and foraging styles. Migratory behaviour and seasonal habitat shifts apparently evolved multiple times in vireos breeding in temperate environments. Present geographical and ecological distributions, and low levels of intrataxon genetic divergence, are hypothesized to be the result of postglacial regionalization of climate–plant associations and rapid northward expansion of breeding ranges.     

Origin and rapid diversification of a tropical moss

Molecular sequences rarely evolve at a constant rate. Yet, even in instances where a clock can be assumed or approximated for a particular set of sequences, fossils or clear patterns of vicariance are rarely available to calibrate the clock. Thus, obtaining absolute timing for diversification of natural lineages can prove difficult. Unfortunately, without absolute time we cannot develop a complete understanding of important evolutionary processes, including adaptive radiations and key innovations. In the present study, the coding sequence of the nuclear gene, glyceraldehyde 3-phosphate dehydrogenase (gpd), extracted from the paleotropical moss, Mitthyridium, was found to exhibit clocklike behavior and used to reconstruct the history of 80 distinct molecular lineages that cover the full geographic range of Mitthyridium. Two separate clades endemic to two geographically distinct oceanic archipelagos were revealed by this phylogenetic analysis. This allowed the use of island age (as derived from potassium-argon dating) as a maximum age of origin of each monophyletic group, providing two independent time anchors for the clock found in gpd, the final piece needed to study absolute time. Based on results from both maximum age calibrations, which separately yielded highly consistent estimates, the ancestor of this moss group arose approximately 8 million years ago, and then diversified at the rapid rate of 0.56 +/- 0.004 new lineages per million years. Such a rate is on par with the highest diversification rates reported in the literature including rapidly radiating insular groups like the Hawaiian silversword alliance, a classic example of an adaptive radiation. Using independent sources of data, it was found that neither the age nor diversification estimates were affected by the use of molecular lineages rather than species as the operational taxonomic units. Identifying the cause for this rapid diversification requires further testing, but it appears to be related to a general shift in reproductive strategy from sexual to asexual, which may be a key innovation for this young group.     

Adaptive radiation and the evolution of nectarivory in a large songbird clade

The accumulation of exceptional ecological diversity within a lineage is a key feature of adaptive radiation resulting from diversification associated with the subdivision of previously underutilized resources. The invasion of unoccupied niche space is predicted to be a key determinant of adaptive diversification, and this process may be particularly important if the diversity of competing lineages within the area, in which the radiation unfolds, is already high. Here, we test whether the evolution of nectarivory resulted in significantly higher rates of morphological evolution, more extensive morphological disparity, and a heightened build‐up of sympatric species diversity in a large adaptive radiation of passerine birds (the honeyeaters, about 190 species) that have diversified extensively throughout continental and insular settings. We find that a large increase in rates of body size evolution and general expansion in morphological space followed an ancestral shift to nectarivory, enabling the build‐up of large numbers of co‐occurring species that vary greatly in size, compared to related and co‐distributed nonnectarivorous clades. These results strongly support the idea that evolutionary shifts into novel areas of niche space play a key role in promoting adaptive radiation in the presence of likely competing lineages.     

Non-monophyly of most supraspecific taxa of calcareous sponges (Porifera, Calcarea) revealed by increased taxon sampling and partitioned Bayesian analysis of ribosomal DNA

Calcareous sponges (Porifera, Calcarea) play an important role for our understanding of early metazoan evolution, since several molecular studies suggested their closer relationship to Eumetazoa than to the other two sponge 'classes,' Demospongiae and Hexactinellida. The division of Calcarea into the subtaxa Calcinea and Calcaronea is well established by now, but their internal relationships remain largely unresolved. Here, we estimate phylogenetic relationships within Calcarea in a Bayesian framework, using full-length 18S and partial 28S ribosomal DNA sequences. Both genes were analyzed separately and in combination and were further partitioned by stem and loop regions, the former being modelled to take non-independence of paired sites into account. By substantially increasing taxon sampling, we show that most of the traditionally recognized supraspecific taxa within Calcinea and Calcaronea are not monophyletic, challenging the existing classification system, while monophyly of Calcinea and Calcaronea is again highly supported.     

The biogeographic legacy of an imperilled taxon provides a foundation for assessing lineage diversification,demography and conservation genetics

Aim To test alternative biogeographic hypotheses related to the diversification of a montane mammal (Zapus hudsonius luteus) endemic to the American Southwest. Location South‐western United States. Methods We used statistical phylogeographic analyses of mitochondrial DNA (1512 bp; two genes) from 93 individuals from six geographic regions to test diversification hypotheses. Species distribution models of climate and fossil records were integrated to assess contemporary and historical distributions and barriers to gene flow. We calculated dates of divergence and examined historical demography using coalescent simulations. Results We documented monophyly of Z. h. luteus represented by 19 segregated haplotypes. Predicted current distribution generally coincided with known localities, while predicted paleodistributions suggested that this lineage was widespread throughout lower elevations of the American Southwest and on the Edwards Plateau (as documented by the fossil record). Population size did not change substantially during a westward shift in range that occurred in the last 100 k generations. Results supported fragmentation of a common ancestor during the Holocene as the most plausible explanation for genetic structure. Main conclusions Monophyletic Z. h. luteus reflects fragmentation of a common ancestor with recent (Holocene) upslope colonization of disjunct montane areas. We refute the hypotheses of in situ divergence or origins from a Colorado Piedmont ancestor. Instead, westward colonization from the Edwards Plateau during the Wisconsin followed by Holocene fragmentation, which serves as a generalized biogeographic hypothesis for species associated with mesic graminoid habitats in the American Southwest. Further exploration of these signatures using independent nuclear DNA is warranted. Key conservation implications are (1) Z. h. luteus is a monophyletic lineage on an independent evolutionary trajectory; (2) Z. h. luteus shared a recent common ancestor with Z. h. pallidus (not Z. h. preblei); (3) mtDNA does not reflect recent population declines; and (4) coalescent simulations and species distribution models reflect Holocene fragmentation.     

Early diversification dynamics in a highly successful insular plant taxon are consistent with the general dynamic model of oceanic island biogeography

The general dynamic model (GDM) of oceanic island biogeography views oceanic islands predominantly as sinks rather than sources of dispersing lineages. To test this, we conducted a biogeographic analysis of a highly successful insular plant taxon, Cyrtandra, and inferred the directionality of dispersal and founder events throughout the four biogeographical units of the Indo-Australian Archipelago (IAA), namely Sunda, Wallacea,  Philippines, and Sahul. Sunda was recovered as the major source area, followed by Wallacea, a system of oceanic islands. The relatively high number of events originating from Wallacea is attributed to its central location in the IAA and its complex geological history selecting for increased dispersibility. We also tested if diversification dynamics in Cyrtandra follow predictions of adaptive radiation, which is the dominant process as per the GDM. Diversification dynamics of dispersing lineages of Cyrtandra in the Southeast Asian grade showed early bursts followed by a plateau, which is consistent with adaptive radiation. We did not detect signals of diversity-dependent diversification, and this is attributed to Southeast Asian cyrtandras occupying various niche spaces, evident by their wide morphological range in habit and floral characters. The Pacific clade, which arrived at the immaturity phase of the Pacific Islands, showed diversification dynamics predicted by the island immaturity speciation pulse model (IISP), wherein rates increase exponentially, and their morphological range is controlled by the least action effect favoring woodiness and fleshy fruits. Our study provides a first step toward a framework for investigating diversification dynamics as predicted by the GDM in highly successful insular taxa.     

Geographical structure,narrow species ranges,and Cenozoic diversification in a pantropical clade of epiphyllous leafy liverworts

The evolutionary history and classification of epiphyllous cryptogams are still poorly known. Leptolejeunea is a largely epiphyllous pantropical liverwort genus with about 25 species characterized by deeply bilobed underleaves, elliptic to narrowly obovate leaf lobes, the presence of ocelli, and vegetative reproduction by cladia. Sequences of three chloroplast regions (rbcL, trnL‐F, psbA) and the nuclear ribosomal ITS region were obtained for 66 accessions of Leptolejeunea and six outgroup species to explore the phylogeny, divergence times, and ancestral areas of this genus. The phylogeny was estimated using maximum‐likelihood and Bayesian inference approaches, and divergence times were estimated with a Bayesian relaxed clock method. Leptolejeunea likely originated in Asia or the Neotropics within a time interval from the Early Eocene to the Late Cretaceous (67.9 Ma, 95% highest posterior density [HPD]: 47.9–93.7). Diversification of the crown group initiated in the Eocene or early Oligocene (38.4 Ma, 95% HPD: 27.2–52.6). Most species clades were established in the Miocene. Leptolejeunea epiphylla and L. schiffneri originated in Asia and colonized African islands during the Plio‐Pleistocene. Accessions of supposedly pantropical species are placed in different main clades. Several monophyletic morphospecies exhibit considerable sequence variation related to a geographical pattern. The clear geographic structure of the Leptolejeunea crown group points to evolutionary processes including rare long‐distance dispersal and subsequent speciation. Leptolejeunea may have benefitted from the large‐scale distribution of humid tropical angiosperm forests in the Eocene.     

Small-sample inference for incomplete longitudinal data with truncation and censoring in tumor xenograft models

In cancer drug development, demonstrating activity in xenograft models, where mice are grafted with human cancer cells, is an important step in bringing a promising compound to humans. A key outcome variable is the tumor volume measured in a given period of time for groups of mice given different doses of a single or combination anticancer regimen. However, a mouse may die before the end of a study or may be sacrificed when its tumor volume quadruples, and its tumor may be suppressed for some time and then grow back. Thus, incomplete repeated measurements arise. The incompleteness or missingness is also caused by drastic tumor shrinkage (<0.01 cm3) or random truncation. Because of the small sample sizes in these models, asymptotic inferences are usually not appropriate. We propose two parametric test procedures based on the EM algorithm and the Bayesian method to compare treatment effects among different groups while accounting for informative censoring. A real xenograft study on a new antitumor agent, temozolomide, combined with irinotecan is analyzed using the proposed methods.