A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

Genome sequences of the human body louse and its primary endosymbiont provide insights into the permanent parasitic lifestyle

As an obligatory parasite of humans, the body louse (Pediculus humanus humanus) is an important vector for human diseases, including epidemic typhus, relapsing fever, and trench fever. Here, we present genome sequences of the body louse and its primary bacterial endosymbiont Candidatus Riesia pediculicola. The body louse has the smallest known insect genome, spanning 108 Mb. Despite its status as an obligate parasite, it retains a remarkably complete basal insect repertoire of 10,773 protein-coding genes and 57 microRNAs. Representing hemimetabolous insects, the genome of the body louse thus provides a reference for studies of holometabolous insects. Compared with other insect genomes, the body louse genome contains significantly fewer genes associated with environmental sensing and response, including odorant and gustatory receptors and detoxifying enzymes. The unique architecture of the 18 minicircular mitochondrial chromosomes of the body louse may be linked to the loss of the gene encoding the mitochondrial single-stranded DNA binding protein. The genome of the obligatory louse endosymbiont Candidatus Riesia pediculicola encodes less than 600 genes on a short, linear chromosome and a circular plasmid. The plasmid harbors a unique arrangement of genes required for the synthesis of pantothenate, an essential vitamin deficient in the louse diet. The human body louse, its primary endosymbiont, and the bacterial pathogens that it vectors all possess genomes reduced in size compared with their free-living close relatives. Thus, the body louse genome project offers unique information and tools to use in advancing understanding of coevolution among vectors, symbionts, and pathogens.

Evolutionary history of Ramphastos toucans: Molecular phylogenetics, temporal diversification, and biogeography

The toucan genus Ramphastos (Piciformes: Ramphastidae) has been a model in the formulation of Neotropical paleobiogeographic hypotheses. Weckstein (2005) reported on the phylogenetic history of this genus based on three mitochondrial genes, but some relationships were weakly supported and one of the subspecies of R. vitellinus (citreolaemus) was unsampled. This study expands on Weckstein (2005) by adding more DNA sequence data (including a nuclear marker) and more samples, including R v. citreolaemus. Maximum parsimony, maximum likelihood, and Bayesian methods recovered similar trees, with nodes showing high support. A monophyletic R. vitellinus complex was strongly supported as the sister-group to R. brevis. The results also confirmed that the southeastern and northern populations of R. vitellinus ariel are paraphyletic. X v. citreolaemus is sister to the Amazonian subspecies of the vitellinus complex. Using three protein-coding genes (COI, cytochrome-b and ND2) and interval-calibrated nodes under a Bayesian relaxed-clock framework, we infer that ramphastid genera originated in the middle Miocene to early Pliocene, Ramphastos species originated between late Miocene and early Pleistocene, and intra-specific divergences took place throughout the Pleistocene. Parsimony-based reconstruction of ancestral areas indicated that evolution of the four trans-Andean Ramphastos taxa (R. v. citreolaemus, R. a. swainsonii, R. brevis and R. sulfuratus) was associated with four independent dispersals from the cis-Andean region. The last pulse of Andean uplift may have been important for the evolution of R. sulfuratus, whereas the origin of the other trans-Andean Ramphastos taxa is consistent with vicariance due to drying events in the lowland forests north of the Andes. Estimated rates of molecular evolution were higher than the ""standard"" bird rate of 2% substitutions/site/million years for two of the three genes analyzed (cytochrome-b and ND2). (C) 2009 Elsevier Inc. All rights reserved.

Gene expression profiling reveals molecular marker candidates of laryngeal squamous cell carcinoma

Laryngeal squamous cell carcinoma is very common in head and neck cancer, with high mortality rates and poor prognosis. In this study, we compared expression profiles of clinical samples from 13 larynx tumors and 10 non-neoplastic larynx tissues using a custom-built cDNA microarray containing 331 probes for 284 genes previously identified by informatics analysis of EST databases as markers of head and neck tumors. Thirty-five genes showed statistically significant differences (SNR >= 11.01, p <= 0.001) in the expression between tumor and non-tumor larynx tissue samples. Functional annotation indicated that these genes are involved in cellular processes relevant to the cancer phenotype, such as apoptosis, cell cycle, DNA repair, proteolysis, protease inhibition, signal transduction and transcriptional regulation. Six of the identified transcripts map to intronic regions of protein-coding genes and may comprise non-annotated exons or as yet uncharacterized long ncRNAs with a regulatory role in the gene expression program of larynx tissue. The differential expression of 10 of these genes (ADCY6, AES, AL2SCR3, CRR9, CSTB, DUSP1, MAP3K5, PLAT, UBL1 and ZNF706) was independently confirmed by quantitative real-time RT-PCR. Among these, the CSTB gene product has cysteine protease inhibitor activity that has been associated with an antimetastatic function. Interestingly, CSTB showed a low expression in the tumor samples analyzed (p<0.0001). The set of genes identified here contribute to a better understanding of the molecular basis of larynx cancer, and provide candidate markers for improving diagnosis, prognosis and treatment of this carcinoma.

Mitochondrial DNA sequence and gene organization in the Australian Blacklip Abalone Haliotis rubra (Leach)

The complete mitochondrial DNA of the blacklip abalone Haliotis rubra (Gastropoda: Mollusca) was cloned and 16,907 base pairs were sequenced. The sequence represents an estimated 99.85% of the mitochondrial genome, and contains 2 ribosomal RNA, 22 transfer RNA, and 13 protein-coding genes found in other metazoan mtDNA. An AT tandem repeat and a possible C-rich domain within the putative control region could not be fully sequenced. The H. rubra mtDNA gene order is novel for mollusks, separated from the black chiton Katharina tunicata by the individual translocations of 3 tRNAs. Compared with other mtDNA regions, sequences from the ATP8, NAD2, NAD4L, NAD6, and 12S rRNA genes, as well as the control region, are the most variable among representatives from Mollusca, Arthropoda, and Rhynchonelliformea, with similar mtDNA arrangements to H. rubra. These sequences are being evaluated as genetic markers within commercially important Haliotis species, and some applications and considerations for their use are discussed.

Phylogeny of the Australian freshwater crayfish Cherax destructor-complex (Decapoda:Parastacidae) inferred from four mitochondrial gene regions

The phylogenetic relationships among 32 individuals of Australian freshwater crayfish belonging to the Cherax destructor-complex were investigated using a dataset comprising sequences from four mitochondrial gene regions: the large subunit rRNA (16S rRNA), cytochrome oxidase I (COI), adenosine triphosphatase 6 (ATPase 6), and cytochrome oxidase III (COIII). A total of 1602 bp was obtained, and a combined analysis of the data produced a tree with strong support (bootstrap values 94–100%) for three divergent lineages, verifying the phylogenetic hypotheses of relationships within the C. destructor species-complex suggested in previous studies. Overall, sequences from the 16S rRNA gene showed the least variation compared to those generated from protein coding genes, which presented considerably greater levels of divergence. The level of divergence within C. destructor was found to be greater than that observed in other species of freshwater crayfish, but interspecific variation among species examined in the present study was similar to that reported previously.

Complete genome sequence of the Medicago microsymbiont Ensifer (Sinorhizobium) medicae strain WSM419

Ensifer (Sinorhizobium) medicae is an effective nitrogen fixing microsymbiont of a diverse range of annual Medicago (medic) species. Strain WSM419 is an aerobic, motile, non-spore forming, Gram-negative rod isolated from a M. murex root nodule collected in Sardinia, Italy in 1981. WSM419 was manufactured commercially in Australia as an inoculant for annual medics during 1985 to 1993 due to its nitrogen fixation, saprophytic competence and acid tolerance properties. Here we describe the basic features of this organism, together with the complete genome sequence, and annotation. This is the first report of a complete genome sequence for a microsymbiont of the group of annual medic species adapted to acid soils. We reveal that its genome size is 6,817,576 bp encoding 6,518 protein-coding genes and 81 RNA only encoding genes. The genome contains a chromosome of size 3,781,904 bp and 3 plasmids of size 1,570,951, 1,245,408 and 219,313 bp. The smallest plasmid is a feature unique to this medic microsymbiont.

Complete genome sequence of Rhizobium leguminosarum bv trifolii strain WSM2304, an effective microsymbiont of the South American clover Trifolium polymorphum

Rhizobium leguminosarum bv. trifolii is the effective nitrogen fixing microsymbiont of a diverse range of annual and perennial Trifolium (clover) species. Strain WSM2304 is an aerobic, motile, non-spore forming, Gram-negative rod isolated from Trifolium polymorphum in Uruguay in 1998. This microsymbiont predominated in the perennial grasslands of Glencoe Research Station, in Uruguay, to competitively nodulate its host, and fix atmospheric nitrogen. Here we describe the basic features of WSM2304, together with the complete genome sequence, and annotation. This is the first completed genome sequence for a nitrogen fixing microsymbiont of a clover species from the American centre of origin. We reveal that its genome size is 6,872,702 bp encoding 6,643 protein-coding genes and 62 RNA only encoding genes. This multipartite genome was found to contain 5 distinct replicons; a chromosome of size 4,537,948 bp and four circular plasmids of size 4,537,948, 1,266,105, 501,946, 308,747 and 257,956 bp.

Complete genome sequence of rhizobium leguminosarum bv. trifolii strain WSM1325, an effective microsymbiont of annual Mediterranean clovers

Rhizobium leguminosarum bv trifolii is a soil-inhabiting bacterium that that has the capacity to be an effective nitrogen fixing microsymbiont of a diverse range of annual Trifolium (clover) species. Strain WSM1325 is an aerobic, motile, non-spore forming, Gram-negative rod isolated from root nodules collected in 1993 from the Greek Island of Serifos. WSM1325 is manufactured commercially in Australia as an inoculant for a broad range of annual clovers of Mediterranean origin due to its superior attributes of saprophytic competence, nitrogen fixation and acid-tolerance. Here we describe the basic features of this organism, together with the complete genome sequence, and annotation. This is the first completed genome sequence for a microsymbiont of annual clovers. We reveal that its genome size is 7,418,122 bp encoding 7,232 protein-coding genes and 61 RNA-only encoding genes. This multipartite genome contains 6 distinct replicons; a chromosome of size 4,767,043 bp and 5 plasmids of size 828,924, 660,973, 516,088, 350,312 and 294,782 bp.

Genome sequence of Ensifer sp. TW10, a Tephrosia wallichii (Biyani) microsymbiont native to the Indian Thar Desert

Ensifer sp. TW10 is novel N2-fixingbacterium isolated from a root nodule of the perennial legume Tephrosia wallichii Graham (known locally as Biyani) found in the Great Indian (or Thar) desert, a large arid regionin the northwestern part of the Indian subcontinent. Strain TW10 is a Gram-negative, rod shaped,aerobic, motile, non-spore forming, species of root nodule bacteria (RNB) that promiscuously nodulates legumes in Thar Desert alkaline soil. It is fast growing, acid-producing, and tolerates up to 2% NaCl and capable of growth at 40C. In this report we describe for the first time the primary features of this Thar Desert soil saprophyte together with genome sequence informationand annotation. The 6,802,256 bp genome has a GC content of 62% and is arranged into 57 scaffolds containing 6,470 protein-coding genes, 73 RNA genes and asingle rRNA operon. This genome is one of 100 RNB genomes sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.

Genome sequence of Ensifer medicae strain WSM1369, an effective microsymbiont of the annual legume Medicago sphaerocarpos

Ensifer medicae WSM1369 is an aerobic, motile, Gram-negative, non-spore-forming rod that can exist as a soil saprophyte or as a legume microsymbiont of Medicago. WSM1369 was isolated in 1993 from a nodule recovered from the roots of Medicago sphaerocarpos growing at San Pietro di Rudas, near Aggius in Sardinia (Italy). WSM1369 is an effective microsymbiont of the annual forage legumes M. polymorpha and M. sphaerocarpos. Here we describe the features of E. medicae WSM1369, together with genome sequence information and its annotation. The 6,402,557 bp standard draft genome is arranged into 307 scaffolds of 307 contigs containing 6,656 protein-coding genes and 79 RNA-only encoding genes. This rhizobial genome is one of 100 sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.

Genome sequence of Ensifer medicae strain WSM1115, an acid-tolerant Medicago-nodulating microsymbiont from Samothraki, Greece

Ensifer medicae (syn. Sinorhizobium medicae) strain WSM1115 forms effective nitrogen fixing symbioses with a range of annual Medicago species and is used in commercial inoculants in Australia. WSM1115 is an aerobic, motile, Gram-negative, non-spore-forming rod. It was isolated from a nodule recovered from the root of burr medic (Medicago polymorpha) collected on the Greek Island of Samothraki. WSM1115 has a broad host range for nodulation and N2 fixation capacity within the genus Medicago, although this does not extend to all medic species. WSM1115 is considered saprophytically competent in moderately acid soils (pHCaCl2 5.0) however, has failed to persist at field sites where soil salinity exceeded 10 ECe (dS/m). Here we describe the features of E. medicae strain WSM1115, together with genome sequence information and its annotation. The 6,861,065 bp high-quality-draft genome is arranged into 7 scaffolds of 28 contigs, contains 6,789 protein-coding genes and 83 RNA-only encoding genes, and is one of 100 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.

Genome sequence of Ensifer arboris strain LMG 14919T, a microsymbiont of the legume Prosopis chilensis growing in Kosti, Sudan

Ensifer arboris LMG 14919T is an aerobic, motile, Gram-negative, non-spore-forming rod that can exist as a soil saprophyte or as a legume microsymbiont of several species of legume trees. LMG 14919T was isolated in 1987 from a nodule recovered from the roots of the tree Prosopis chilensis growing in Kosti, Sudan. LMG 14919T is highly effective at fixing nitrogen with P. chilensis (Chilean mesquite) and Acacia senegal (gum Arabic tree or gum acacia). LMG 14919T does not nodulate the tree Leucena leucocephala, nor the herbaceous species Macroptilium atropurpureum, Trifolium pratense, Medicago sativa, Lotus corniculatus and Galega orientalis. Here we describe the features of E. arboris LMG 14919T, together with genome sequence information and its annotation. The 6,850,303 bp high-quality-draft genome is arranged into 7 scaffolds of 12 contigs containing 6,461 protein-coding genes and 84 RNA-only encoding genes, and is one of 100 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.

Genome sequence of Ensifer meliloti strain WSM1022, a highly effective microsymbiont of the model legume Medicago truncatula A17

Ensifer meliloti WSM1022 is an aerobic, motile, Gram-negative, non-spore-forming rod that can exist as a soil saprophyte or as a legume microsymbiont of Medicago. WSM1022 was isolated in 1987 from a nodule recovered from the roots of the annual Medicago orbicularis growing on the Cyclades Island of Naxos in Greece. WSM1022 is highly effective at fixing nitrogen with M. truncatula and other annual species such as M. tornata and M. littoralis and is also highly effective with the perennial M. sativa (alfalfa or lucerne). In common with other characterized E. meliloti strains, WSM1022 will nodulate but fixes poorly with M. polymorpha and M. sphaerocarpos and does not nodulate M. murex. Here we describe the features of E. meliloti WSM1022, together with genome sequence information and its annotation. The 6,649,661 bp high-quality-draft genome is arranged into 121 scaffolds of 125 contigs containing 6,323 protein-coding genes and 75 RNA-only encoding genes, and is one of 100 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.

Genome sequence of Burkholderia mimosarum strain LMG 23256T, a mimosa pigra microsymbiont from Anso, Taiwan

Burkholderia mimosarum strain LMG 23256T is an aerobic, motile, Gram-negative, non-spore-forming rod that can exist as a soil saprophyte or as a legume microsymbiont of Mimosa pigra (giant sensitive plant). LMG 23256T was isolated from a nodule recovered from the roots of the M. pigra growing in Anso, Taiwan. LMG 23256T is highly effective at fixing nitrogen with M. pigra. Here we describe the features of B. mimosarum strain LMG 23256T, together with genome sequence information and its annotation. The 8,410,967 bp high-quality-draft genome is arranged into 268 scaffolds of 270 contigs containing 7,800 protein-coding genes and 85 RNA-only encoding genes, and is one of 100 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.