Phylogenetic and evolutionary implications of nuclear ribosomal DNA variation in dwarf dandelions (Krigia,Lactuceae, Asteraceae)

Restriction site and length variations of nrDNA were examined for 51 populations of seven species ofKrigia. The nrDNA repeat ranged in size from 8.7 to 9.6 kilobase (kb). The transcribed region, including the two ITSs, was 5.35 kb long in all examinedKrigia populations. In contrast, the size of the nontranscribed IGS varied from 3.35 to 4.25 kb. Eight different types of length-variations were identified among the 51 populations, including distinct nrDNA lengths in the tetraploid and diploid populations of bothK. biflora andK. virginica. However, a few variations were detected among populations of the same species or within a cytotype. All populations ofKrigia sect.Cymbia share a 600 bp insertion in IGS near the 18 S gene, and this feature suggests monophyly of the section. AllKrigia spp. had a conjugated type of subrepeat composed of approximately 75 basepairs (bp) and 125 bp. Base modifications in the gene coding regions were highly conserved among species. Forty-five restriction sites from 15 enzymes were mapped, 24 of which were variable among populations. Only four of the variable sites occurred in the rRNA coding region while 20 variable sites were detected in the noncoding regions. Collectively, 25 enzymes generated about 66 restriction sites in each nrDNA; this amounts to about 4.3% of the nrDNA repeat. A total of 50 restriction sites was variable, 28 of which were phylogenetically informative. Phylogenetic analyses of site mutations indicated that two sections ofKrigia, sect.Cymbia and sect.Krigia, are monophyletic. In addition, relationships among several species were congruent with other sources of data, such as cpDNA restriction site variation and morphology. Both length and restriction site variation supported an allopolyploid origin of the hexaploidK. montana. The average sequence divergence value inKrigia nrDNA was 40 times greater than that of the chloroplast DNA. The rapid evolution of nrDNA sequences was primarily due to changes of the IGS sequences.     

Taxonomy ofHordeum caespitosum,H. jubatum,andH. lechleri (Poaceae: Triticeae)

Hordeum caespitosum Scribner,H. jubatum L., andH. lechleri (Steudel)Schenck are very similar in appearance and therefore until recently were mostly not recognized as separate entities. The first two are tetraploid and natives to North America, but the second occurs naturally in eastern Siberia and has been introduced in Europe and South America and may become a cosmopolitan weed. The third is hexaploid and South American. This study analyses their morphological diversity by means of selected multivariate techniques in order to determine if there is justification to recognize them as three separate morphological species. Logistic discrimination, although based on a reduced set of characters, yielded the highest percent of correct assignments. A linear discriminant function is provided and validated by 100 bootstrap repeats. Canonical discriminant analysis indicated three groups. It is subsequently concluded that the three are separate morphological species. Although a linear discriminant function is given, a traditional identification key is provided based on the palea length and triad (the group of three spikelets at each rachis node) length.     

Equivalence of regression curves sharing common parameters

In clinical trials, the comparison of two different populations is a common problem. Nonlinear (parametric) regression models are commonly used to describe the relationship between covariates, such as concentration or dose, and a response variable in the two groups. In some situations, it is reasonable to assume some model parameters to be the same, for instance, the placebo effect or the maximum treatment effect. In this paper, we develop a (parametric) bootstrap test to establish the similarity of two regression curves sharing some common parameters. We show by theoretical arguments and by means of a simulation study that the new test controls its significance level and achieves a reasonable power. Moreover, it is demonstrated that under the assumption of common parameters, a considerably more powerful test can be constructed compared with the test that does not use this assumption. Finally, we illustrate the potential applications of the new methodology by a clinical trial example.     

Hemocyanin gene family evolution in spiders (Araneae), with implications for phylogenetic relationships and divergence times in the infraorder Mygalomorphae

Hemocyanins are multimeric copper-containing hemolymph proteins involved in oxygen binding and transport in all major arthropod lineages. Most arachnids have seven primary subunits (encoded by paralogous genes a–g), which combine to form a 24-mer (4 × 6) quaternary structure. Within some spider lineages, however, hemocyanin evolution has been a dynamic process with extensive paralog duplication and loss. We have obtained hemocyanin gene sequences from numerous representatives of the spider infraorders Mygalomorphae and Araneomorphae in order to infer the evolution of the hemocyanin gene family and estimate spider relationships using these conserved loci. Our hemocyanin gene tree is largely consistent with the previous hypotheses of paralog relationships based on immunological studies, but reveals some discrepancies in which paralog types have been lost or duplicated in specific spider lineages. Analyses of concatenated hemocyanin sequences resolved deep nodes in the spider phylogeny and recovered a number of clades that are supported by other molecular studies, particularly for mygalomorph taxa. The concatenated data set is also used to estimate dates of higher-level spider divergences and suggests that the diversification of extant mygalomorphs preceded that of extant araneomorphs. Spiders are diverse in behavior and respiratory morphology, and our results are beneficial for comparative analyses of spider respiration. Lastly, the conserved hemocyanin sequences allow for the inference of spider relationships and ancient divergence dates.     

Dim Nighttime Illumination Interacts With Parametric Effects of Bright Light to Increase the Stability of Circadian Rhythm Bifurcation in Hamsters

The endogenous circadian pacemaker of mammals is synchronized to the environmental day by the ambient cycle of relative light and dark. The present studies assessed the actions of light in a novel circadian entrainment paradigm where activity rhythms are bifurcated following exposure to a 24-h light:dark:light:dark (LDLD) cycle. Bifurcated entrainment under LDLD reflects the temporal dissociation of component oscillators that comprise the circadian system and is facilitated when daily scotophases are dimly lit rather than completely dark. Although bifurcation can be stably maintained in LDLD, it is quickly reversed under constant conditions. Here the authors examine whether dim scotophase illumination acts to maintain bifurcated entrainment under LDLD through potential interactions with the parametric actions of bright light during the two daily photophases. In three experiments, wheel-running rhythms of Syrian hamsters were bifurcated under LDLD with dimly lit scotophases, and after several weeks, dim scotophase illumination was either retained or extinguished. Additionally, “full” and “skeleton” photophases were employed under LDLD cycles with dimly lit or completely dark scotophases to distinguish parametric from nonparametric effects of bright light. Rhythm bifurcation was more stable in full versus skeleton LDLD cycles. Dim light facilitated the maintenance of bifurcated entrainment under full LDLD cycles but did not prevent the loss of rhythm bifurcation in skeleton LDLD cycles. These studies indicate that parametric actions of bright light maintain the bifurcated entrainment state; that dim scotophase illumination increases the stability of the bifurcated state; and that dim light interacts with the parametric effects of bright light to increase the stability of rhythm bifurcation under full LDLD cycles. A further understanding of the novel actions of dim light may lead to new strategies for understanding, preventing, and treating chronobiological disturbances. (Author correspondence: jevans@msm.edu)     

Confidence interval of intrinsic optimum temperature estimated using thermodynamic SSI model

The intrinsic optimum temperature for the development of ectotherms is one of the most important factors not only for their physiological processes but also for ecological and evolutional processes. The Sharpe–Schoolfield–Ikemoto (SSI) model succeeded in defining the temperature that can thermodynamically meet the condition that at a particular temperature the probability of an active enzyme reaching its maximum activity is realized. Previously, an algorithm was developed by Ikemoto (Tropical malaria does not mean hot environments. Journal of Medical Entomology, 45, 963–969) to estimate model parameters, but that program was computationally very time consuming. Now, investigators can use the SSI model more easily because a full automatic computer program was designed by Shi et al. (A modified program for estimating the parameters of the SSI model. Environmental Entomology, 40, 462–469). However, the statistical significance of the point estimate of the intrinsic optimum temperature for each ectotherm has not yet been determined. Here, we provided a new method for calculating the confidence interval of the estimated intrinsic optimum temperature by modifying the approximate bootstrap confidence intervals method. For this purpose, it was necessary to develop a new program for a faster estimation of the parameters in the SSI model, which we have also done.     

Phylogeography of the dusky shrew, Sorex monticolus (Insectivora, Soricidae): insight into deep and shallow history in northwestern North America

Phylogenetic relationships among the dusky shrew (Sorexmonticolus) and eight related species (S. bairdi, S. bendirii, S. neomexicanus, S.ornatus, S. pacificus, S. palustris, S. sonomae and S.vagrans) were assessed using sequences from the mitochondrial cytochrome b gene (801 bp). Analyses using parsimony and maximum likelihood revealed significant molecular variation not reflected in previous morphological studies of these species. Conversely, three morphologically defined species (S.bairdi, S.neomexicanus and S.pacificus) were poorly differentiated. Sorexornatus and S.vagrans represented basal taxa for a more inclusive group that included: (i) a widespread Continental clade containing S.monticolus (Arizona to Alaska, including S. neomexicanus); (ii) a Coastal clade containing S.monticolus (Oregon to south-east Alaska, including S. bairdi and S. pacificus); (iii) the semiaquatic species (S. bendirii and S. palustris); and (iv) S.sonomae. Additional subdivision was observed within the Continental clade corresponding to populations from the northern and southern Rocky Mountains. Average uncorrected sequence divergence between the Coastal and Continental clades was 5.3% (range 4.5-6.2%), which exceeds many interspecific comparisons within this species complex and within the genus Sorex. Lack of resolution of internal nodes within topologies suggests a deep history of rapid diversification within this group. Late Pleistocene/Holocene glacial perturbations are reflected in the shallow phylogeographic structure within these clades in western North America. Our results suggest also that S. monticolus is not monophyletic under current taxonomic nomenclature. This perspective on phylogeographic history was developed within a growing comparative framework for other organisms in western North America.     

Bayes or bootstrap? A simulation study comparing the performance of Bayesian Markov chain Monte Carlo sampling and bootstrapping in assessing phylogenetic confidence

Bayesian Markov chain Monte Carlo sampling has become increasingly popular in phylogenetics as a method for both estimating the maximum likelihood topology and for assessing nodal confidence. Despite the growing use of posterior probabilities, the relationship between the Bayesian measure of confidence and the most commonly used confidence measure in phylogenetics, the nonparametric bootstrap proportion, is poorly understood. We used computer simulation to investigate the behavior of three phylogenetic confidence methods: Bayesian posterior probabilities calculated via Markov chain Monte Carlo sampling (BMCMC-PP), maximum likelihood bootstrap proportion (ML-BP), and maximum parsimony bootstrap proportion (MP-BP). We simulated the evolution of DNA sequence on 17-taxon topologies under 18 evolutionary scenarios and examined the performance of these methods in assigning confidence to correct monophyletic and incorrect monophyletic groups, and we examined the effects of increasing character number on support value. BMCMC-PP and ML-BP were often strongly correlated with one another but could provide substantially different estimates of support on short internodes. In contrast, BMCMC-PP correlated poorly with MP-BP across most of the simulation conditions that we examined. For a given threshold value, more correct monophyletic groups were supported by BMCMC-PP than by either ML-BP or MP-BP. When threshold values were chosen that fixed the rate of accepting incorrect monophyletic relationship as true at 5%, all three methods recovered most of the correct relationships on the simulated topologies, although BMCMC-PP and ML-BP performed better than MP-BP. BMCMC-PP was usually a less biased predictor of phylogenetic accuracy than either bootstrapping method. BMCMC-PP provided high support values for correct topological bipartitions with fewer characters than was needed for nonparametric bootstrap.