Inferring phylogenies from RAD sequence data

Reduced-representation genome sequencing represents a new source of data for systematics, and its potential utility in interspecific phylogeny reconstruction has not yet been explored. One approach that seems especially promising is the use of inexpensive short-read technologies (e.g., Illumina, SOLiD) to sequence restriction-site associated DNA (RAD)--the regions of the genome that flank the recognition sites of restriction enzymes. In this study, we simulated the collection of RAD sequences from sequenced genomes of different taxa (Drosophila, mammals, and yeasts) and developed a proof-of-concept workflow to test whether informative data could be extracted and used to accurately reconstruct "known" phylogenies of species within each group. The workflow consists of three basic steps: first, sequences are clustered by similarity to estimate orthology; second, clusters are filtered by taxonomic coverage; and third, they are aligned and concatenated for "total evidence" phylogenetic analysis. We evaluated the performance of clustering and filtering parameters by comparing the resulting topologies with well-supported reference trees and we were able to identify conditions under which the reference tree was inferred with high support. For Drosophila, whole genome alignments allowed us to directly evaluate which parameters most consistently recovered orthologous sequences. For the parameter ranges explored, we recovered the best results at the low ends of sequence similarity and taxonomic representation of loci; these generated the largest supermatrices with the highest proportion of missing data. Applications of the method to mammals and yeasts were less successful, which we suggest may be due partly to their much deeper evolutionary divergence times compared to Drosophila (crown ages of approximately 100 and 300 versus 60 Mya, respectively). RAD sequences thus appear to hold promise for reconstructing phylogenetic relationships in younger clades in which sufficient numbers of orthologous restriction sites are retained across species.     

A statistical test of phylogenies estimated from sequence data

A simple approach to testing the significance of the branching order, estimated from protein or DNA sequence data, of three taxa is proposed. The branching order is inferred by the transformed-distance method, under the assumption that one or two outgroups are available, and the branch lengths are estimated by the least-squares method. The inferred branching order is considered significant if the estimated internodal distance is significantly greater than zero. To test this, a formula for the variance of the internodal distance has been developed. The statistical test proposed has been checked by computer simulation. The same test also applies to the case of four taxa with no outgroup, if one considers an unrooted tree. Formulas for the variances of internodal distances have also been developed for the case of five taxa. Conditions are given under which it is more efficient to add the sequence of a fifth taxon than to do 25% more nucleotide sequencing in each of the original four. A method is presented for combining analyses of disparate data to get a single P value. Finally, the test, applied to the human-chimpanzee-gorilla problem, shows that the issue is not yet resolved.     

Generation and identification of DNA sequence flanking T-DNA integration site of Trichoderma atroviride mutants with high dichlorvos-degrading capacity

A protocol for efficient Agrobacterium tumefaciens-mediated transformation (ATMT) of biocontrol fungus Trichoderma atroviride strain T23 was developed to construct mutants with improved dichlorvos-degradation ability. A transformation frequency of 5 × 10⁻⁶ was achieved. Among 110 genetically stable T-DNA transformants of T. atroviride T23, two transformants, AMT-12 and AMT-28, confirmed by Southern blot analysis to have single-copy inserts of T-DNA, showed an increase in dichlorvos-degradation ability of more than 10% compared to that of the wild type, exhibited similar tolerance to the pesticide, but lower spore formation ability. Five transformants exhibited a reduction in degradation of more than 70%, exhibited wild-type spore formation, and tolerated up to 800 μg/mL of dichlorvos. The left-flanking sequence of the insertion site in AMT-12 was cloned as a 1845-bp fragment and shown to have 89% identity to the DNA from T. atroviride IMI 206040; however, the involvement of this DNA in dichlorvos degradation remains still to be determined. This study can promote both a more efficient isolation of DNA sequence flanking T-DNA integration site in T. atroviride mutants and a more rational utilization of these transformants in dichlorvos degradation.     

Phylogenetic relationships of Australian strongyloid nematodes inferred from ribosomal DNA sequence data

The sequence of the second internal transcribed spacer of the ribosomal DNA was determined for the following strongyloid nematodes: Cylicocyclus insignis, Chabertia ovina, Oesophagostomum venulosum, Cloacina communis, Cloacina hydriformis, Labiostrongylus labiostrongylus, Parazoniolaimus collaris, Macropostrongylus macropostrongylus, Macropostrongylus yorkei, Rugopharynx australis, Rugopharynx rosemariae, Macropostrongyloides baylisi, Oesophagostomoides longispicularis and Paramacropostrongylus toraliformis, and compared with published sequences for species of Strongylus and for Hypodontus macropi. The resultant phylogenetic trees supported current hypotheses based on morphological evidence for the separation of the families Strongylidae and Chabertiidae, but did not support the separation of the endemic Australian genera as a distinctive clade within the Chabertiidae. The implications of this finding for the phylogenetic origins of the Australian strongyloids are discussed.     

Restriction endonuclease analysis of ribosomal DNA from Saccharomyces cerevisiae

Summary Temperature-sensitive mutants of Escherichia coli that are unable to grow at high temperature can be obtained among those selected for resistance to streptovaricin or rifampicin at low temperature (Yura et al., 1970). One of these mutants (KY5323) that was supposed to carry a single mutation affecting both rifampicin resistance and temperature sensitivity was further investigated. Using purified RNA polymerase preparations obtained from the mutant and the wild type, it was found that the activity for DNA chain elongation is more sensitive to heat treatment than that for RNA chain initiation or DNA binding, and that the mutant enzyme is significantly more labile than the wild-type enzyme with respect to RNA chain elongation, when heat treatment is carried out at high salt concentration. These results, taken together with those of the enzyme reconstitution experiments, strongly suggest that the subunit of the polymerase is directly involved in both RNA chain initiation and elongation reactions. Enzyme reconstitution experiments using isolated subunits derived from the mutant and the wild-type polymerases demonstrate that the alteration of subunit is primarily responsible for both rifampicin resistance and thermolability of the mutant enzyme. In addition, the results suggested the apparent alteration of both and subunits in this mutant. Extensive transduction experiments provided genetic evidence that are consistent with the view that the strain KY5323 carries a second mutation affecting the subunit, beside the primary mutation affecting the subunit. The hypothetical subunit mutation seems to modify quantitatively the rifampicin resistance caused by the subunit mutation.     

Evolutionary relationships of trichostrongyloid nematodes (Strongylida) inferred from ribosomal DNA sequence data

The evolutionary relationships of 21 species of trichostrongyloid nematodes were determined by use of sequence data of the second internal transcribed spacer of the ribosomal DNA aligned according to secondary structure information. Irrespective of the method of analysis used, the topologies of the phylogenetic trees derived from the molecular data differed with respect to all four hypotheses proposed previously for the evolutionary relationships of the different subfamilies within the Trichostrongylidae based on morphological data. Thus, the molecular data set did not resolve the conflict between the four previous proposals for the subfamilial relationships. Nonetheless, all trees derived from the molecular data showed strong support for the exclusion of the genera Filarinema and Amidostomum from the clade containing the species within the family Trichostrongylidae. This represents a major difference from the most recent proposal of the systematics of the Trichostrongyloidea in which these two genera were included within the Trichostrongylidae. Therefore, the molecular data support an earlier systematic framework in which Filarinema and Amidostomum were considered to be sister groups of the Trichostrongyloidea.     

Two new Fusicoccum species from Acacia and Eucalyptus in Venezuela, based on morphology and DNA sequence data

Botryosphaeria spp. are common endophytes of woody plants, and they also include some serious pathogens of Eucalyptus and Acacia species. Numerous anamorphs have been associated with Botryosphaeria, of which the species Fusicoccum are amongst the most common. Here, we characterize two new Fusicoccum species, isolated from Eucalyptus and Acacia trees in Venezuela, based on morphological features in culture and comparisons of DNA sequence data. The two taxa named Fusicoccum andinum and F. stromaticum spp. nov, reside in two well-supported clades (BS values = 100 %) based on a combined data set of the ITS of the rDNA operon and translation elongation factor 1-α (EF1- α) gene sequences. The conidia of F. andinum are unusually large amongst Botryosphaeria anamorphs, and peripherally resemble those of B. mamane and B. melanops. F. stromaticum is characterized by large conidiomata in cultures, growth at 35 °C and slightly thickened conidial walls, which is different to most other Fusicoccum spp. No teleomorphs were observed for these fungi, but DNA sequence data show that they are anamorphs of Botryosphaeria.     

DNA sequence of Nodulin-45 from Lupinus angustifolius

A partial cDNA clone encoding Lupinus angustifolius Nodulin-45 was isolated by differential hybridisation. A genomic clone was also isolated, from which the DNA sequence was obtained for the 5′ end of the gene (including 1.2 kb of 5′ upstream region). The upstream region includes putative cis-elements, found upstream of other nodulin genes. Southern analysis indicates the presence of several Nodulin-45-like sequences in the lupin genome. The Nodulin-45 protein has a putative N-terminal endoplasmic reticulum-type signal sequence and also contains a large glycine-rich repeat sequence. The cDNA sequence is highly homologous to a Nodulin-45 cDNA sequence from Lupinus luteus (Szczyglowski et al., Plant Sci., 65 (1989) 87–95), although major sequence rearrangements are apparent between the L. luteus and L. angustifolius cDNAs.     

Congruency of phylogenies derived from different proteins

Summary This communication examines the question of phylogenetic congruency- i.e., whether or not the branching order of evolutionary trees is independent of the protein studied. It was found that trees constructed for birds on the basis of immunological comparison of their transferrins, albumins, and ovalbumins agree approximately with a published tree based on the amino acid sequences of their lysozymesc. This congruency is especially noteworthy with respect to the phylogenetic position of the chachalaca, a Mexican bird classified on morphological grounds in the family Cracidae of the order Galliformes. At the protein level, this species differs as much from non-cracid galliform birds as does the duck, which belongs to another order. Despite the organismal similarity between cracid and non-cracid galliform birds, the molecular relationship is remote. If this contrast between organismal and molecular results had been based on comparative studies with only lysozyme, one could have ascribed the contrast to the possibility that chachalaca lysozyme was paralogous, rather than orthologous, to the other bird lysozymesc. Examination of several proteins is thus desirable in cases of possible paralogy.This work was supported in part by grants GB-42028X from NSF and GM-21509 from NIH     

Trichoderma species form endophytic associations within Theobroma cacao trichomes

Trichoderma species are usually considered soil organisms that colonize plant roots, sometimes forming a symbiotic relationship. Recent studies demonstrate that Trichoderma species are also capable of colonizing the above ground tissues of Theobroma cacao (cacao) in what has been characterized as an endophytic relationship. Trichoderma species can be re-isolated from surface sterilized cacao stem tissue, including the bark and xylem, the apical meristem, and to a lesser degree from leaves. SEM analysis of cacao stems colonized by strains of four Trichoderma species (Trichoderma ovalisporum-DIS 70a, Trichoderma hamatum-DIS 219b, Trichoderma koningiopsis-DIS 172ai, or Trichoderma harzianum-DIS 219f) showed a preference for surface colonization of glandular trichomes versus non-glandular trichomes. The Trichoderma strains colonized the glandular trichome tips and formed swellings resembling appresoria. Hyphae were observed emerging from the glandular trichomes on surface sterilized stems from cacao seedlings that had been inoculated with each of the four Trichoderma strains. Fungal hyphae were observed under the microscope emerging from the trichomes as soon as 6 h after their isolation from surface sterilized cacao seedling stems. Hyphae were also observed, in some cases, emerging from stalk cells opposite the trichome head. Repeated single trichome/hyphae isolations verified that the emerging hyphae were the Trichoderma strains with which the cacao seedlings had been inoculated. Strains of four Trichoderma species were able to enter glandular trichomes during the colonization of cacao stems where they survived surface sterilization and could be re-isolated. The penetration of cacao trichomes may provide the entry point for Trichoderma species into the cacao stem allowing systemic colonization of this tissue.     

应用ITS序列分析探讨偏花报春的系统位置

对于偏花报春Primula secundiflora Franch. 的系统位置,主要有两种意见,一种认为偏花报春具有典型的钟状花冠,应置于钟花报春组Sect．Sikkimensis;而另一种意见则依据其他特征将其置于灯台报春组sect．Proliferae。通过对偏花报春、灯台报春组4种植物和钟花报春组3种植物核糖体DNA中的内转录间隔区(ITS)的序列测定及分析,并结合形态学及染色体特征的比较论证,认为偏花报春应置于灯台报春组。     

Nuclear ribosomal DNA sequence polymorphism and hybridization in checker mallows (Sidalcea, Malvaceae)

Checker mallows (Sidalcea, Malvaceae) constitute a western North American genus of annuals and perennials that have been regarded as taxonomically difficult because of complex patterns of morphological variation putatively stemming from hybridization and polyploidy. In recent molecular phylogenetic investigations extensive polymorphism was observed in the internal and external transcribed spacers (ITS and ETS) of 18S–26S nuclear ribosomal DNA for some Sidalcea samples. To resolve the evolutionary basis for this polymorphism and to readdress the evolutionary impact of hybridization in Sidalcea we cloned and sequenced the polymorphic DNAs and included the clones in phylogenetic analyses together with direct sequences of non-polymorphic samples. The positions of cloned spacer sequences in the phylogenetic trees suggest that S. reptans and two subspecies of S. malviflora may have been influenced by past hybridization with lineages of the “glaucescens” clade. Polymorphic sequence patterns in other taxa may be a result of extensive interbreeding within young clades, in keeping with the minimal sequence divergence, largely overlapping geographic distributions and morphology, and ploidy variation in these groups. Other possible explanations for polymorphic sequences in members of Sidalcea include slow concerted evolution relative to mutation rates, incomplete lineage sorting, and recent pseudogene formation.     

Bootstrapping phylogenies inferred from rearrangement data

ABSTRACT: Large-scale sequencing of genomes has enabled the inference of phylogenies based on the evolution of genomic architecture, under such events as rearrangements, duplications, and losses. Many evolutionary models and associated algorithms have been designed over the last few years and have found use in comparative genomics and phylogenetic inference. However, the assessment of phylogenies built from such data has not been properly addressed to date. The standard method used in sequence-based phylogenetic inference is the bootstrap, but it relies on a large number of homologous characters that can be resampled; yet in the case of rearrangements, the entire genome is a single character. Alternatives such as the jackknife suffer from the same problem, while likelihood tests cannot be applied in the absence of well established probabilistic models. We present a new approach to the assessment of distance-based phylogenetic inference from whole-genome data; our approach combines features of the jackknife and the bootstrap and remains nonparametric. For each feature of our method, we give an equivalent feature in the sequence-based framework; we also present the results of extensive experimental testing, in both sequence-based and genome-based frameworks. Through the feature-by-feature comparison and the experimental results, we show that our bootstrapping approach is on par with the classic phylogenetic bootstrap used in sequence-based reconstruction, and we establish the clear superiority of the classic bootstrap for sequence data and of our corresponding new approach for rearrangement data over proposed variants. Finally, we test our approach on a small dataset of mammalian genomes, verifying that the support values match current thinking about the respective branches. Our method is the first to provide a standard of assessment to match that of the classic phylogenetic bootstrap for aligned sequences. Its support values follow a similar scale and its receiver-operating characteristics are nearly identical, indicating that it provides similar levels of sensitivity and specificity. Thus our assessment method makes it possible to conduct phylogenetic analyses on whole genomes with the same degree of confidence as for analyses on aligned sequences. Extensions to search-based inference methods such as maximum parsimony and maximum likelihood are possible, but remain to be thoroughly tested.     

Phylogeny of Caragana (Fabaceae) based on DNA sequence data from rbcL, trnS–trnG, and ITS

Phylogenetic relationships of 48 species of Caragana (Fabaceae: tribe Hedysareae) and one representative each of Astragalus, Calophaca, Halimodendron, and Hedysarum are estimated from DNA sequences of the rbcL gene, trnS–trnG intron and spacer, and ITS region. At least one representative of all five sections and 12 series within Caragana are included. Analyses yielded strongly supported clades corresponding to sections Caragana, Bracteolatae, and Frutescentes. The species of section Jubatae are distributed among three strongly supported clades, i.e., one with the species of section Bracteolatae, another with two species of section Spinosae, and a third as sister to section Frutescentes. All but the last of these six clades are corroborated by at least one unambiguously traced morphological character. The placement of the other four species of section Spinosae are not well supported and lack unambiguous morphological synapomorphies, and the samples of Calophaca and Halimodendron nest within Caragana with weak support.     

Traditional infrageneric classification of Gymnopilus is not supported by ribosomal DNA sequence data

The traditional classification of Gymnopilus (Agaricales) recognizes two primary groups, Annulati and Gymnopilus, based on the presence or absence of a membranous partial veil. While our analyses of DNA sequence data from the nuclear ribosomal ITS1-5.8S-ITS2 (ITS) gene supports the monophyly of the genus, these traditional subgroups were not recovered. Five well-supported clades within the genus were identified through these analyses: 1) the spectabilis-imperialis group; 2) nevadensis-penetrans group; 3) a clade formed by G. underwoodii, G. validipes and G. cf. flavidellus; 4) aeruginosus-luteofolius group; and 5) lepidotus-subearlei group. Relationships among these subgroups were not resolved.     

Cyanide catabolizing enzymes in Trichoderma spp.

An investigation was made into the occurrence and distribution of the enzymes involved in HCN catabolism in different strains of the fungus Trichoderma. Three enzymes, cyanide hydratase, rhodanese and β-cyanoalanine synthase were studied. All the strains showed a high capacity to degrade cyanide via both the cyanide hydratase and rhodanese pathways. β-Cyanoalanine synthase, however, was not observed in any of the strains. The enzyme activities were found in varying levels in each of the Trichoderma strains. Experiments conducted with cyanide addition to the medium to assess whether the enzymes were induced in the presence of cyanide failed to show any statistically significant increase. This suggests a constitutive nature of both the enzymes in all the selected strains of Trichoderma used in this study.     

Biodegradation of cyanide by Trichoderma spp. and Fusarium spp.

The rate of degradation of cyanide by certain strains of the Trichoderma spp. was evaluated. For comparison two Fusarium spp., which had previously been demonstrated to degrade metallocyanides were also studied. Studies were carried out to assess the rate of degradation using cyanide as the sole source of carbon or in the presence of glucose. Biodegradation was observed in flask cultures using cyanide as the sole carbon source. Strong evidence of cyanide biodegradation and co-metabolism emerged from studies with flask cultures where glucose was provided as a co-substrate. The rate of degradation of 2000 ppm CN⁻ was enhanced almost three times in the presence of glucose. A concomitant increase in fungal biomass was also observed in all the strains over the experimental period. Growth yield calculations performed provided values that were comparable to those reported in literature for one-carbon substrates.