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.     

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.     

Influence of mutational and sampling factors on the estimation of demographic parameters in a "continuous" population under isolation by distance

In numerous species, individual dispersal is restricted in space so that "continuous" populations evolve under isolation by distance. A method based on individual genotypes assuming a lattice population model was recently developed to estimate the product Dsigma2, where D is the population density and sigma2 is the average squared parent-offspring distance. We evaluated the influence on this method of both mutation rate and mutation model, with a particular reference to microsatellite markers, as well as that of the spatial scale of sampling. Moreover, we developed and tested a nonparametric bootstrap procedure allowing the construction of confidence intervals for the estimation of Dsigma2. These two objectives prompted us to develop a computer simulation algorithm based on the coalescent theory giving individual genotypes for a continuous population under isolation by distance. Our results show that the characteristics of mutational processes at microsatellite loci, namely the allele size homoplasy generated by stepwise mutations, constraints on allele size, and change of slippage rate with repeat number, have little influence on the estimation of Dsigma2. In contrast, a high genetic diversity (approximately 0.7-0.8), as is commonly observed for microsatellite markers, substantially increases the precision of the estimation. However, very high levels of genetic diversity (>0.85) were found to bias the estimation. We also show that statistics taking into account allele size differences give unreliable estimations (i.e., high variance of Dsigma2 estimation) even under a strict stepwise mutation model. Finally, although we show that this method is reasonably robust with respect to the sampling scale, sampling individuals at a local geographical scale gives more precise estimations of Dsigma2.     

Simplicial edge representation of protein structures and alpha contact potential with confidence measure

Protein representation and potential function are two important ingredients for studying protein folding, equilibrium thermodynamics, and sequence design. We introduce a novel geometric representation of protein contact interactions using the edge simplices from the alpha shape of the protein structure. This representation can eliminate implausible neighbors that are not in physical contact, and can avoid spurious contact between two residues when a third residue is between them. We developed statistical alpha contact potential using an odds-ratio model. A studentized bootstrap method was then introduced to assess the 95% confidence intervals for each of the 210 propensity parameters. We found, with confidence, that there is significant long-range propensity (>30 residues apart) for hydrophobic interactions. We tested alpha contact potential for native structure discrimination using several sets of decoy structures, and found that it often performs comparably with atom-based potentials requiring many more parameters. We also show that accurate geometric representation is important, and that alpha contact potential has better performance than potential defined by cutoff distance between geometric centers of side chains. Hierarchical clustering of alpha contact potentials reveals natural grouping of residues. To explore the relationship between shape and physicochemical representations, we tested the minimum alphabet size necessary for native structure discrimination. We found that there is no significant difference in performance of discrimination when alphabet size varies from 7 to 20, if geometry is represented accurately by alpha simplicial edges. This result suggests that the geometry of packing plays an important role, but the specific residue types are often interchangeable.