Past climate cooling and orogenesis of the Hengduan Mountains have influenced the evolution of Impatiens sect. Impatiens (Balsaminaceae) in the Northern Hemisphere.

BACKGROUND: Impatiens sect. Impatiens is distributed across the Northern Hemisphere and has diversified considerably, particularly within the Hengduan Mountains (HDM) in southwest China. Yet, the infra-sectional phylogenetic relationships are not well resolved, largely due to limited taxon sampling and an insufficient number of molecular markers. The evolutionary history of its diversification is also poorly understood. In this study, plastome data and the most complete sampling to date were used to reconstruct a robust phylogenetic framework for this section. The phylogeny was then used to investigate its biogeographical history and diversification patterns, specifically with the aim of understanding the role played by the HDM and past climatic changes in its diversification. RESULTS: A stable phylogeny was reconstructed that strongly supported both the monophyly of the section and its division into seven major clades (Clades I-VII). Molecular dating and ancestral area reconstruction suggest that sect. Impatiens originated in the HDM and Southeast China around 11.76 Ma, after which different lineages dispersed to Northwest China, temperate Eurasia, and North America, mainly during the Pliocene and Pleistocene. An intercontinental dispersal event from East Asia to western North America may have occurred via the Bering Land Bridge or Aleutian Islands. The diversification rate was high during its early history, especially with the HDM, but gradually decreased over time both within and outside the HDM. Multiple linear regression analysis showed that the distribution pattern of species richness was strongly associated with elevation range, elevation, and mean annual temperature. Finally, ancestral niche analysis indicated that sect. Impatiens originated in a relatively cool, middle-elevation area. CONCLUSIONS: We inferred the evolutionary history of sect. Impatiens based on a solid phylogenetic framework. The HDM was the primary source or pump of its diversity in the Northern Hemisphere. Orogeny and climate change may have also shaped its diversification rates, as a steady decrease in the diversification rate coincided with the uplift of the HDM and climate cooling. These findings provide insights into the distribution pattern of sect. Impatiens and other plants in the Northern Hemisphere.

De-novo whole genome assembly of the orange jewelweed, Impatiens capensis Meerb. (Balsaminaceae) using nanopore long-read sequencing.

The plant family Balsaminaceae comprises only two genera, and they are a study in contrasts. While Impatiens is an impressively prolific genus, with over 1,000 species and more being discovered each year, its sister genus, Hydrocera, has one solitary species, H. triflora. The two genera also differ in geographic distribution and habitat type (Impatiens species are widely distributed in much of the Old World and N. America, while H. triflora is confined to wetlands specific to S. India, Sri Lanka, and SE Asia). Other contrasting features include plant habit, habitat, floral architecture, mode of seed dispersal, and a host of other traits. The family Balsaminaceae is therefore an excellent model for studying speciation and character evolution as well as understanding the proximal and evolutionary forces that have driven the two genera to adopt such contrasting evolutionary paths. Various species of the Impatiens genus are also commercially important in the ornamental flower industry and as sources of phytochemicals that are of medicinal and other commercial value. As a preliminary step towards studying the genomic basis of the contrasting features of the two genera, we have sequenced and assembled, de novo, the genome of an iconic Impatiens species from N. America, namely I. capensis, and report our findings here.

PhyloM: A Computer Program for Phylogenetic Inference from Measurement or Binary Data, with Bootstrapping.

Quantitative and binary results are ubiquitous in biology. Inasmuch as an underlying genetic basis for the observed variation in these observations can be assumed, it is pertinent to infer the evolutionary relationships among the entities being measured. I present a computer program, PhyloM, that takes measurement data or binary data as input, using which, it directly generates a pairwise distance matrix that can then be subjected to the popular neighbor-joining (NJ) algorithm to produce a phylogenetic tree. PhyloM also has the option of nonparametric bootstrapping for testing the level of support for the inferred phylogeny. Finally, PhyloM also allows the user to root the tree on any desired branch. PhyloM was tested on Biolog Gen III growth data from isolates within the genus Chromobacterium and the closely related Aquitalea sp. This allowed a comparison with the genotypic tree inferred from whole-genome sequences for the same set of isolates. From this comparison, it was possible to infer parallel evolution. PhyloM is a stand-alone and easy-to-use computer program with a user-friendly graphical user interface that computes pairwise distances from measurement or binary data, which can then be used to infer phylogeny using NJ using a utility in the same program. Alternatively, the distance matrix can be downloaded for use in another program for phylogenetic inference or other purposes. It does not require any software to be installed or computer code written and is open source. The executable and computer code are available on GitHub.

Insights into plant biodiversity conservation in large river valleys in China: A spatial analysis of species and phylogenetic diversity.

Large river valleys (LRVs) are heterogeneous in habitat and rich in biodiversity, but they are largely overlooked in policies that prioritize conservation. Here, we aimed to identify plant diversity hotspots along LRVs based on species richness and spatial phylogenetics, evaluate current conservation effectiveness, determine gaps in the conservation networks, and offer suggestions for prioritizing conservation. We divided the study region into 50 km × 50 km grid cells and determined the distribution patterns of seed plants by studying 124,927 occurrence points belonging to 14,481 species, using different algorithms. We generated phylogenies for the plants using the "V. PhyloMaker" R package, determined spatial phylogenetics, and conducted correlation analyses between different distribution patterns and spatial phylogenetics. We evaluated the effectiveness of current conservation practices and discovered gaps of hotspots within the conservation networks. In the process, we identified 36 grid cells as hotspots (covering 10% of the total area) that contained 83.4% of the species. Fifty-eight percent of the hotspot area falls under the protection of national nature reserves (NNRs) and 83% falls under national and provincial nature reserves (NRs), with 42% of the area identified as conservation gaps of NNRs and 17% of the area as gaps of NRs. The hotspots contained high proportions of endemic and threatened species, as did conservation gaps. Therefore, it is necessary to optimize the layout of current conservation networks, establish micro-nature reserves, conduct targeted conservation priority planning focused on specific plant groups, and promote conservation awareness. Our results show that the conservation of three hotspots in Southwest China, in particular, is likely to positively affect the protection of biodiversity in the LRVs, especially with the participation of the neighboring countries, India, Myanmar, and Laos.

Impatiensyunlingensis (Balsaminaceae), a new species from Yunnan, China.

Impatiensyunlingensis S.X. Yu, Chang Y. Xia & J.H. Yu (Balsaminaceae), a species new to science discovered in Yunnan, China, is described and illustrated here, along with its phylogenetic position among other Impatiens species. Morphological, micro-morphological and molecular evidence is presented as an attestation of its novelty. Impatiensyunlingensis is similar to I.delavayi in having coarsely crenate leave margins, bracts in the upper part, ca. 4/5 length of the pedicels, saccate lower sepal with shallowly bifid spur, linear capsules, and elliptic-oblong, tuberculate seeds, but differs from I.delavayi with lateral sepals 4 (vs. 2), lateral united petal basal lobes subtriangular (vs. dolabriform), and seeds' surface equipped with tubercular ornamentation mostly covered with grain shaped appendages (vs. glabrous and without grain shaped appendages on the top).

Computational tools for fitting the Hill equation to dose-response curves.

INTRODUCTION: Many biological response curves commonly assume a sigmoidal shape that can be approximated well by means of the 4-parameter nonlinear logistic equation, also called the Hill equation. However, estimation of the Hill equation parameters requires access to commercial software or the ability to write computer code. Here we present two user-friendly and freely available computer programs to fit the Hill equation - a Solver-based Microsoft Excel template and a stand-alone GUI-based "point and click" program, called HEPB. METHODS: Both computer programs use the iterative method to estimate two of the Hill equation parameters (EC50 and the Hill slope), while constraining the values of the other two parameters (the minimum and maximum asymptotes of the response variable) to fit the Hill equation to the data. In addition, HEPB draws the prediction band at a user-defined confidence level, and determines the EC50 value for each of the limits of this band to give boundary values that help objectively delineate sensitive, normal and resistant responses to the drug being tested. RESULTS: Both programs were tested by analyzing twelve datasets that varied widely in data values, sample size and slope, and were found to yield estimates of the Hill equation parameters that were essentially identical to those provided by commercial software such as GraphPad Prism and nls, the statistical package in the programming language R. DISCUSSION: The Excel template provides a means to estimate the parameters of the Hill equation and plot the regression line in a familiar Microsoft Office environment. HEPB, in addition to providing the above results, also computes the prediction band for the data at a user-defined level of confidence, and determines objective cut-off values to distinguish among response types (sensitive, normal and resistant). Both programs are found to yield estimated values that are essentially the same as those from standard software such as GraphPad Prism and the R-based nls. Furthermore, HEPB also has the option to simulate 500 response values based on the range of values of the dose variable in the original data and the fit of the Hill equation to that data.

An airtight approach to the inebriometer: from construction to application with volatile anesthetics.

The mechanism by which inhaled anesthetics work is not fully understood, although they have been extensively used. Much research has been done showing the likelihood that there is more than one pathway or mechanism of action. A long-term goal of our laboratory is to decipher these mechanisms using Drosophila melanogaster, an excellent model organism for this purpose. In order to do this, we have modified and constructed an apparatus called the inebriometer to quantitatively analyze the response of flies to inhaled anesthetics. While the inebriometer is not new to the fly community, our updated design provides a relatively low-labor, high-throughput means for performing screens in search of genes involved in the anesthetic mechanism. Here we describe our construction of an airtight inebriometer that we have designed for this purpose. We also provide data that validates this apparatus and establishes a procedure for its use.

Chromobacterium vaccinii sp. nov., isolated from native and cultivated cranberry (Vaccinium macrocarpon Ait.) bogs and irrigation ponds.

A large number of Gram-negative, motile, mesophilic, violacein-producing bacteria were isolated from the soils and roots of Vaccinium macrocarpon Ait. and Kalmia angustifolia L. plants and from irrigation ponds associated with wild and cultivated cranberry bogs in Massachusetts, USA. Phylogenetic analyses of 16S rRNA gene sequences placed these isolates in a clade with Chromobacterium species, but the specialized environment from which they were isolated, their low genomic DNA relatedness with Chromobacterium violaceum ATCC 12472(T) and C. subtsugae PRAA4-1(T), significant differences in fatty acid composition and colony morphology indicate that the cranberry and Kalmia isolates comprise a separate species of Chromobacterium, for which the name Chromobacterium vaccinii sp. nov. is proposed. Strain MWU205(T) ( = ATCC BAA-2314(T)  = DSM 25150(T)) is proposed as the type strain for the novel species. Phenotypic analysis of 26 independent isolates of C. vaccinii sp. nov. indicates that, despite close geographical and biological proximity, there is considerable metabolic diversity among individuals within the population.

Tubulin evolution in insects: gene duplication and subfunctionalization provide specialized isoforms in a functionally constrained gene family.

BACKGROUND: The completion of 19 insect genome sequencing projects spanning six insect orders provides the opportunity to investigate the evolution of important gene families, here tubulins. Tubulins are a family of eukaryotic structural genes that form microtubules, fundamental components of the cytoskeleton that mediate cell division, shape, motility, and intracellular trafficking. Previous in vivo studies in Drosophila find a stringent relationship between tubulin structure and function; small, biochemically similar changes in the major alpha 1 or testis-specific beta 2 tubulin protein render each unable to generate a motile spermtail axoneme. This has evolutionary implications, not a single non-synonymous substitution is found in beta 2 among 17 species of Drosophila and Hirtodrosophila flies spanning 60 Myr of evolution. This raises an important question, How do tubulins evolve while maintaining their function? To answer, we use molecular evolutionary analyses to characterize the evolution of insect tubulins. RESULTS: Sixty-six alpha tubulins and eighty-six beta tubulin gene copies were retrieved and subjected to molecular evolutionary analyses. Four ancient clades of alpha and beta tubulins are found in insects, a major isoform clade (alpha 1, beta 1) and three minor, tissue-specific clades (alpha 2-4, beta 2-4). Based on a Homarus americanus (lobster) outgroup, these were generated through gene duplication events on major beta and alpha tubulin ancestors, followed by subfunctionalization in expression domain. Strong purifying selection acts on all tubulins, yet maximum pairwise amino acid distances between tubulin paralogs are large (0.464 substitutions/site beta tubulins, 0.707 alpha tubulins). Conversely orthologs, with the exception of reproductive tissue isoforms, show little sequence variation except in the last 15 carboxy terminus tail (CTT) residues, which serve as sites for post-translational modifications (PTMs) and interactions with microtubule-associated proteins. CTT residues overwhelming comprise the co-evolving residues between Drosophila alpha 2 and beta 3 tubulin proteins, indicating CTT specializations can be mediated at the level of the tubulin dimer. Gene duplications post-dating separation of the insect orders are unevenly distributed, most often appearing in major alpha 1 and minor beta 2 clades. More than 40 introns are found in tubulins. Their distribution among tubulins reveals that insertion and deletion events are common, surprising given their potential for disrupting tubulin coding sequence. Compensatory evolution is found in Drosophila beta 2 tubulin cis-regulation, and reveals selective pressures acting to maintain testis expression without the use of previously identified testis cis-regulatory elements. CONCLUSION: Tubulins have stringent structure/function relationships, indicated by strong purifying selection, the loss of many gene duplication products, alpha-beta co-evolution in the tubulin dimer, and compensatory evolution in beta 2 tubulin cis-regulation. They evolve through gene duplication, subfunctionalization in expression domain and divergence of duplication products, largely in CTT residues that mediate interactions with other proteins. This has resulted in the tissue-specific minor insect isoforms, and in particular the highly diverse alpha3, alpha4, and beta2 reproductive tissue-specific tubulin isoforms, illustrating that even a highly conserved protein family can participate in the adaptive process and respond to sexual selection.

Phylogenetic inference under varying proportions of indel-induced alignment gaps.

BACKGROUND: The effect of alignment gaps on phylogenetic accuracy has been the subject of numerous studies. In this study, we investigated the relationship between the total number of gapped sites and phylogenetic accuracy, when the gaps were introduced (by means of computer simulation) to reflect indel (insertion/deletion) events during the evolution of DNA sequences. The resulting (true) alignments were subjected to commonly used gap treatment and phylogenetic inference methods. RESULTS: (1) In general, there was a strong--almost deterministic--relationship between the amount of gap in the data and the level of phylogenetic accuracy when the alignments were very "gappy", (2) gaps resulting from deletions (as opposed to insertions) contributed more to the inaccuracy of phylogenetic inference, (3) the probabilistic methods (Bayesian, PhyML & "MLepsilon, " a method implemented in DNAML in PHYLIP) performed better at most levels of gap percentage when compared to parsimony (MP) and distance (NJ) methods, with Bayesian analysis being clearly the best, (4) methods that treat gapped sites as missing data yielded less accurate trees when compared to those that attribute phylogenetic signal to the gapped sites (by coding them as binary character data--presence/absence, or as in the MLepsilon method), and (5) in general, the accuracy of phylogenetic inference depended upon the amount of available data when the gaps resulted from mainly deletion events, and the amount of missing data when insertion events were equally likely to have caused the alignment gaps. CONCLUSION: When gaps in an alignment are a consequence of indel events in the evolution of the sequences, the accuracy of phylogenetic analysis is likely to improve if: (1) alignment gaps are categorized as arising from insertion events or deletion events and then treated separately in the analysis, (2) the evolutionary signal provided by indels is harnessed in the phylogenetic analysis, and (3) methods that utilize the phylogenetic signal in indels are developed for distance methods too. When the true homology is known and the amount of gaps is 20 percent of the alignment length or less, the methods used in this study are likely to yield trees with 90-100 percent accuracy.

Maximum likelihood outperforms maximum parsimony even when evolutionary rates are heterotachous.

Heterotachy occurs when the relative evolutionary rates among sites are not the same across lineages. Sequence alignments are likely to exhibit heterotachy with varying severity because the intensity of purifying selection and adaptive forces at a given amino acid or DNA sequence position is unlikely to be the same in different species. In a recent study, the influence of heterotachy on the performance of different phylogenetic methods was examined using computer simulation for a four-species phylogeny. Maximum parsimony (MP) was reported to generally outperform maximum likelihood (ML). However, our comparisons of MP and ML methods using the methods and evaluation criteria employed in that study, but considering the possible range of proportions of sites involved in heterotachy, contradict their findings and indicate that, in fact, ML is significantly superior to MP even under heterotachy.

Inferring species phylogenies from multiple genes: concatenated sequence tree versus consensus gene tree.

Phylogenetic trees from multiple genes can be obtained in two fundamentally different ways. In one, gene sequences are concatenated into a super-gene alignment, which is then analyzed to generate the species tree. In the other, phylogenies are inferred separately from each gene, and a consensus of these gene phylogenies is used to represent the species tree. Here, we have compared these two approaches by means of computer simulation, using 448 parameter sets, including evolutionary rate, sequence length, base composition, and transition/transversion rate bias. In these simulations, we emphasized a worst-case scenario analysis in which 100 replicate datasets for each evolutionary parameter set (gene) were generated, and the replicate dataset that produced a tree topology showing the largest number of phylogenetic errors was selected to represent that parameter set. Both randomly selected and worst-case replicates were utilized to compare the consensus and concatenation approaches primarily using the neighbor-joining (NJ) method. We find that the concatenation approach yields more accurate trees, even when the sequences concatenated have evolved with very different substitution patterns and no attempts are made to accommodate these differences while inferring phylogenies. These results appear to hold true for parsimony and likelihood methods as well. The concatenation approach shows >95% accuracy with only 10 genes. However, this gain in accuracy is sometimes accompanied by reinforcement of certain systematic biases, resulting in spuriously high bootstrap support for incorrect partitions, whether we employ site, gene, or a combined bootstrap resampling approach. Therefore, it will be prudent to report the number of individual genes supporting an inferred clade in the concatenated sequence tree, in addition to the bootstrap support.