An automated system for interrogating the evolution of microbial endosymbiosis

Inter-kingdom endosymbiotic interactions between bacteria and eukaryotic cells are critical to human health and disease. However, the molecular mechanisms that drive the emergence of endosymbiosis remain obscure. Here, we describe the development of a microfluidic system, named SEER (S̲ystem for the E̲volution of E̲ndosymbiotic R̲elationships), that automates the evolutionary selection of bacteria with enhanced intracellular survival and persistence within host cells, hallmarks of endosymbiosis. Using this system, we show that a laboratory strain of Escherichia coli that initially possessed limited abilities to survive within host cells, when subjected to SEER selection, rapidly evolved to display a 55-fold enhancement in intracellular survival. Notably, molecular dissection of the evolved strains revealed that a single-point mutation in a flexible loop of CpxR, a gene regulator that controls bacterial stress responses, substantially contributed to this intracellular survival. Taken together, these results establish SEER as the first microfluidic system for investigating the evolution of endosymbiosis, show the importance of CpxR in endosymbiosis, and set the stage for evolving bespoke inter-kingdom endosymbiotic systems with novel or emergent properties.

Identification of new protein-protein and protein-DNA interactions linked with wood formation in Populus trichocarpa

Cellular processes, such as signal transduction and cell wall deposition, are organized by macromolecule interactions. Experimentally determined protein-protein interactions (PPIs) and protein-DNA interactions (PDIs) relevant to woody plant development are sparse. To begin to develop a Populus trichocarpa Torr. & A. Gray wood interactome, we applied the yeast-two-hybrid (Y2H) assay in different ways to enable the discovery of novel PPIs and connected networks. We first cloned open reading frames (ORFs) for 361 genes markedly upregulated in secondary xylem compared with secondary phloem and performed a binary Y2H screen with these proteins. By screening a xylem cDNA library for interactors of a subset of these proteins and then recapitulating the process by using a subset of the interactors as baits, we ultimately identified 165 PPIs involving 162 different ORFs. Thirty-eight transcription factors (TFs) included in our collection of P. trichocarpa wood ORFs were used in a Y1H screen for binding to promoter regions of three genes involved in lignin biosynthesis resulting in 40 PDIs involving 20 different TFs. The network incorporating both the PPIs and PDIs included 14 connected subnetworks, with the largest having 132 members. Protein-protein interactions and PDIs validated previous reports and also identified new candidate wood formation proteins and modules through their interactions with proteins and promoters known to be involved in secondary cell wall synthesis. Selected examples are discussed including a PPI between Mps one binder (MOB1) and a mitogen-activated protein kinase kinase kinase kinase (M4K) that was further characterized by assays confirming the PPI as well as its effect on subcellular localization. Mapping of published transcriptomic data showing developmentally detailed expression patterns across a secondary stem onto the network supported that the PPIs and PDIs are relevant to wood formation, and also illustrated that wood-associated interactions involve gene products that are not upregulated in secondary xylem.

Abstract 066: Enhancer Repertoires That Define Renin Cell Identity

Control of the renin cell phenotype is crucial for the regulation of blood pressure and fluid- electrolyte homeostasis. Enhancers are cis -acting DNA sequences that harbor distinct chromatin features and regulate gene expression in an orientation-independent manner. Recently, clusters of enhancers or super-enhancers (SE) highly enriched with master transcription factors, possessing open chromatin configuration and in close proximity to cell-identity genes have been proposed. We tested the hypothesis that renin cells have unique repertoires of enhancers and super-enhancers, distinct from other cell types. Those regulatory clusters may in turn confer the identity of renin cells. To define the genome-wide enhancer landscape characteristic of renin cells, we studied As4.1 cells, kidney tumor cells that express renin constitutively, and native renin cells sorted from the kidneys of Ren1cKO-YFP + mice. In these mice, the renin promoter drives YFP expression thus marking the renin cells. We used genome-wide ChIP-Seq for Med1 (subunit 1 of the Mediator complex), H3K27Ac (active enhancers) and Pol II (to visualize putative genomic areas undergoing transcription). The ROSE algorithm we used to ascertain super-enhancers. Chromatin accessibility genome-wide was assessed using ATAC-Seq. The results were compared to twenty-one other cell types that do not express renin. In As4.1 cells, we identified 14,871 enhancers based on H3K27Ac. Of those, 888 were classified as super-enhancers. The Med1 signal in As4.1 cells showed a SE localized 5kb upstream the Ren1 gene, which was ranked at position 25 among other SEs. The H3K27Ac signal showed highest occupancy in the same region. ChIP-Seq for H3K27Ac in YFP + cells showed 211 SEs of 2,987 peaks. The SE for the renin gene possessed the highest signal and ranked number 1, indicating its importance in renin cells. One hundred and thirteen SEs were unique to renin cells, including the SE associated with the renin gene. ATAC-Seq signals overlapped with the renin SE and the classical enhancer indicating that the chromatin was accessible for transcription. In summary, renin-expressing cells possess distinct repertoires of unique enhancers and super-enhancers that acting in concert are likely to determine the renin phenotype.

Characterization and RNA-seq analysis of underperformer, an activation-tagged potato mutant

The potato cv. Bintje and a Bintje activation-tagged mutant, underperformer (up) were compared. Mutant up plants grown in vitro were dwarf, with abundant axillary shoot growth, greater tuber yield, altered tuber traits and early senescence compared to wild type. Under in vivo conditions, the dwarf and early senescence phenotypes of the mutant remained, but the up plants exhibited a lower tuber yield and fewer axillary shoots compared to wild type. Southern blot analyses indicated a single T-DNA insertion in the mutant, located on chromosome 10. Initial PCR-based gene expression studies indicated transcriptional activation/repression of several genes in the mutant flanking the insertion. The gene immediately flanking the right border of the T-DNA insertion, which encoded an uncharacterized Broad complex, Tramtrac, Bric-a-brac; also known as Pox virus and Zinc finger (BTB/POZ) domain-containing protein (StBTB/POZ1) containing an Armadillo repeat region, was up-regulated in the mutant. Global gene expression comparisons between Bintje and up using RNA-seq on leaves from 60 day-old plants revealed a dataset of over 1,600 differentially expressed genes. Gene expression analyses suggested a variety of biological processes and pathways were modified in the mutant, including carbohydrate and lipid metabolism, cell division and cell cycle activity, biotic and abiotic stress responses, and proteolysis.

Comparative phylogenomics and evolution of the Brucellae reveal a path to virulence

Brucella species include important zoonotic pathogens that have a substantial impact on both agriculture and human health throughout the world. Brucellae are thought of as "stealth pathogens" that escape recognition by the host innate immune response, modulate the acquired immune response, and evade intracellular destruction. We analyzed the genome sequences of members of the family Brucellaceae to assess its evolutionary history from likely free-living soil-based progenitors into highly successful intracellular pathogens. Phylogenetic analysis split the genus into two groups: recently identified and early-dividing "atypical" strains and a highly conserved "classical" core clade containing the major pathogenic species. Lateral gene transfer events brought unique genomic regions into Brucella that differentiated them from Ochrobactrum and allowed the stepwise acquisition of virulence factors that include a type IV secretion system, a perosamine-based O antigen, and systems for sequestering metal ions that are absent in progenitors. Subsequent radiation within the core Brucella resulted in lineages that appear to have evolved within their preferred mammalian hosts, restricting their virulence to become stealth pathogens capable of causing long-term chronic infections.

An active ac/ds transposon system for activation tagging in tomato cultivar m82 using clonal propagation

Tomato (Solanum lycopersicum) is a model organism for Solanaceae in both molecular and agronomic research. This project utilized Agrobacterium tumefaciens transformation and the transposon-tagging construct Activator (Ac)/Dissociator (Ds)-ATag-Bar_gosGFP to produce activation-tagged and knockout mutants in the processing tomato cultivar M82. The construct carried hygromycin resistance (hyg), green fluorescent protein (GFP), and the transposase (TPase) of maize (Zea mays) Activator major transcript X054214.1 on the stable Ac element, along with a 35S enhancer tetramer and glufosinate herbicide resistance (BAR) on the mobile Ds-ATag element. An in vitro propagation strategy was used to produce a population of 25 T0 plants from a single transformed plant regenerated in tissue culture. A T1 population of 11,000 selfed and cv M82 backcrossed progeny was produced from the functional T0 line. This population was screened using glufosinate herbicide, hygromycin leaf painting, and multiplex polymerase chain reaction (PCR). Insertion sites of transposed Ds-ATag elements were identified through thermal asymmetric interlaced PCR, and resulting product sequences were aligned to the recently published tomato genome. A population of 509 independent, Ds-only transposant lines spanning all 12 tomato chromosomes has been developed. Insertion site analysis demonstrated that more than 80% of these lines harbored Ds insertions conducive to activation tagging. The capacity of the Ds-ATag element to alter transcription was verified by quantitative real-time reverse transcription-PCR in two mutant lines. The transposon-tagged lines have been immortalized in seed stocks and can be accessed through an online database, providing a unique resource for tomato breeding and analysis of gene function in the background of a commercial tomato cultivar.

Transposon tagging in diploid strawberry

Fragaria vesca was transformed with a transposon tagging construct harbouring amino terminally deleted maize transposase and EGFP (Ac element), NPTII, CaMV 35S promoter (P35S) driving transposase and mannopine synthase promoter (Pmas) driving EGFP (Ds element). Of 180 primary transgenics, 48 were potential launch pads, 72 were multiple insertions or chimaeras, and 60 exhibited somatic transposition. T₁ progeny of 32 putative launch pads were screened by multiplex PCR for transposition. Evidence of germ-line transposition occurred in 13 putative launch pads; however, the transposition frequency was too low in three for efficient recovery of transposants. The transposition frequency in the remaining launch pads ranged from 16% to 40%. After self-pollination of the T₀ launch pads, putative transposants in the T₁ generation were identified by multiplex PCR. Sequencing of hiTAIL-PCR products derived from nested primers within the Ds end sequences (either P35S at the left border or the inverted repeat at the right border) of T₁ plants revealed transposition of the Ds element to distant sites in the strawberry genome. From more than 2400 T₁ plants screened, 103 unique transposants have been identified, among which 17 were somatic transpositions observed in the T₀ generation. Ds insertion sites were dispersed among various gene elements [exons (15%), introns (23%), promoters (30%), 3' UTRs (17%) as well as intergenically (15%)]. Three-primer (one on either side of the Ds insertion and one within the Ds T-DNA) PCR could be used to identify homozygous T₂ transposon-tagged plants. The mutant collection has been catalogued in an on-line database.

SNP discovery and genetic mapping of T-DNA insertional mutants in Fragaria vesca L

As part of a program to develop forward and reverse genetics platforms in the diploid strawberry [Fragaria vesca L.; (2n = 2x = 14)] we have generated insertional mutant lines by T-DNA mutagenesis using pCAMBIA vectors. To characterize the T-DNA insertion sites of a population of 108 unique single copy mutants, we utilized thermal asymmetric interlaced PCR (hiTAIL-PCR) to amplify the flanking region surrounding either the left or right border of the T-DNA. Bioinformatics analysis of flanking sequences revealed little preference for insertion site with regard to G/C content; left borders tended to retain more of the plasmid backbone than right borders. Primers were developed from F. vesca flanking sequences to attempt to amplify products from both parents of the reference F. vesca 815 x F. bucharica 601 mapping population. Polymorphism occurred as: presence/absence of an amplification product for 16 primer pairs and different size products for 12 primer pairs, For 46 mutants, where polymorphism was not found by PCR, the amplification products were sequenced to reveal SNP polymorphism. A cleaved amplified polymorphic sequence/derived cleaved amplified polymorphism sequence (CAPS/dCAPS) strategy was then applied to find restriction endonuclease recognition sites in one of the parental lines to map the SNP position of 74 of the T-DNA insertion lines. BLAST search of flanking regions against GenBank revealed that 46 of 108 flanking sequences were close to presumed strawberry genes related to annotated genes from other plants.

Genome sequences of three agrobacterium biovars help elucidate the evolution of multichromosome genomes in bacteria

The family Rhizobiaceae contains plant-associated bacteria with critical roles in ecology and agriculture. Within this family, many Rhizobium and Sinorhizobium strains are nitrogen-fixing plant mutualists, while many strains designated as Agrobacterium are plant pathogens. These contrasting lifestyles are primarily dependent on the transmissible plasmids each strain harbors. Members of the Rhizobiaceae also have diverse genome architectures that include single chromosomes, multiple chromosomes, and plasmids of various sizes. Agrobacterium strains have been divided into three biovars, based on physiological and biochemical properties. The genome of a biovar I strain, A. tumefaciens C58, has been previously sequenced. In this study, the genomes of the biovar II strain A. radiobacter K84, a commercially available biological control strain that inhibits certain pathogenic agrobacteria, and the biovar III strain A. vitis S4, a narrow-host-range strain that infects grapes and invokes a hypersensitive response on nonhost plants, were fully sequenced and annotated. Comparison with other sequenced members of the Alphaproteobacteria provides new data on the evolution of multipartite bacterial genomes. Primary chromosomes show extensive conservation of both gene content and order. In contrast, secondary chromosomes share smaller percentages of genes, and conserved gene order is restricted to short blocks. We propose that secondary chromosomes originated from an ancestral plasmid to which genes have been transferred from a progenitor primary chromosome. Similar patterns are observed in select Beta- and Gammaproteobacteria species. Together, these results define the evolution of chromosome architecture and gene content among the Rhizobiaceae and support a generalized mechanism for second-chromosome formation among bacteria.

A robust species tree for the alphaproteobacteria

The branching order and coherence of the alphaproteobacterial orders have not been well established, and not all studies have agreed that mitochondria arose from within the Rickettsiales. A species tree for 72 alphaproteobacteria was produced from a concatenation of alignments for 104 well-behaved protein families. Coherence was upheld for four of the five orders with current standing that were represented here by more than one species. However, the family Hyphomonadaceae was split from the other Rhodobacterales, forming an expanded group with Caulobacterales that also included Parvularcula. The three earliest-branching alphaproteobacterial orders were the Rickettsiales, followed by the Rhodospirillales and then the Sphingomonadales. The principal uncertainty is whether the expanded Caulobacterales group is more closely associated with the Rhodobacterales or the Rhizobiales. The mitochondrial branch was placed within the Rickettsiales as a sister to the combined Anaplasmataceae and Rickettsiaceae, all subtended by the Pelagibacter branch. Pelagibacter genes will serve as useful additions to the bacterial outgroup in future evolutionary studies of mitochondrial genes, including those that have transferred to the eukaryotic nucleus.

GeneTrees: a phylogenomics resource for prokaryotes

The GeneTrees phylogenomics system pursues comparative genomic analyses from the perspective of gene phylogenies for individual genes. The GeneTrees project has the goal of providing detailed evolutionary models for all protein-coding gene components of the fully sequenced genomes. Currently, a database of alignments and trees for all protein sequences for 325 fully sequenced and annotated prokaryote genomes is available. The prokaryote database contains 890 000 protein sequences organized into over 100 000 alignments, each described by a phylogenetic tree. An original homology group discovery tool assembles sets of related proteins from all versus all pairwise alignments. Multiple alignments for each homology group are stored and subjected to phylogenetic tree inference. A graphical web interface provides visual exploration of the GeneTrees database. Homology groups can be queried by sequence identifiers or annotation terms. Genomes can be browsed visually on a gene map of each chromosome or plasmid. Phylogenetic trees with support values are displayed in conjunction with the associated sequence alignment. A variety of classes of information can be selected to label the tree tips to aid in visual evaluation of annotation and gene function. This web interface is available at .

Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis

Draft genome sequences have been determined for the soybean pathogen Phytophthora sojae and the sudden oak death pathogen Phytophthora ramorum. Oömycetes such as these Phytophthora species share the kingdom Stramenopila with photosynthetic algae such as diatoms, and the presence of many Phytophthora genes of probable phototroph origin supports a photosynthetic ancestry for the stramenopiles. Comparison of the two species' genomes reveals a rapid expansion and diversification of many protein families associated with plant infection such as hydrolases, ABC transporters, protein toxins, proteinase inhibitors, and, in particular, a superfamily of 700 proteins with similarity to known oömycete avirulence genes.

The xylem and phloem transcriptomes from secondary tissues of the Arabidopsis root-hypocotyl

The growth of secondary xylem and phloem depends on the division of cells in the vascular cambium and results in an increase in the diameter of the root and stem. Very little is known about the genetic mechanisms that control cambial activity and the differentiation of secondary xylem and phloem cell types. To begin to identify new genes required for vascular cell differentiation and function, we performed genome-wide expression profiling of xylem and phloem-cambium isolated from the root-hypocotyl of Arabidopsis (Arabidopsis thaliana). Gene expression in the remaining nonvascular tissue was also profiled. From these transcript profiles, we assembled three sets of genes with expression significantly biased toward xylem, phloem-cambium, or nonvascular tissue. We also assembled three two-tissue sets of genes with expression significantly biased toward xylem/phloem-cambium, xylem/nonvascular, or phloem-cambium/nonvascular tissues. Localizations predicted by transcript profiles were supported by results from promoter-reporter and reverse transcription-polymerase chain reaction experiments with nine xylem- or phloem-cambium-biased genes. An analysis of the members of the phloem-cambium gene set suggested that some genes involved in regulating primary meristems are also regulators of the cambium. Secondary phloem was implicated in the synthesis of auxin, glucosinolates, cytokinin, and gibberellic acid. Transcript profiles also supported the importance of class III HD ZIP and KANADI transcription factors as regulators of radial patterning during secondary growth, and identified several members of the G2-like, NAC, AP2, MADS, and MYB transcription factor families that may play roles as regulators of xylem or phloem cell differentiation and activity.

PIML: the Pathogen Information Markup Language

MOTIVATION

A vast amount of information about human, animal and plant pathogens has been acquired, stored and displayed in varied formats through different resources, both electronically and otherwise. However, there is no community standard format for organizing this information or agreement on machine-readable format(s) for data exchange, thereby hampering interoperation efforts across information systems harboring such infectious disease data.

RESULTS

The Pathogen Information Markup Language (PIML) is a free, open, XML-based format for representing pathogen information. XSLT-based visual presentations of valid PIML documents were developed and can be accessed through the PathInfo website or as part of the interoperable web services federation known as ToolBus/PathPort. Currently, detailed PIML documents are available for 21 pathogens deemed of high priority with regard to public health and national biological defense. A dynamic query system allows simple queries as well as comparisons among these pathogens. Continuing efforts are being taken to include other groups' supporting PIML and to develop more PIML documents.

AVAILABILITY

All the PIML-related information is accessible from http://www.vbi.vt.edu/pathport/pathinfo/

The S8 serine, C1A cysteine and A1 aspartic protease families in Arabidopsis

The Arabidopsis thaliana genome has over 550 protease sequences representing all five catalytic types: serine, cysteine, aspartic acid, metallo and threonine (MEROPS peptidase database, http://merops.sanger.ac.uk/), which probably reflect a wide variety of as yet unidentified functions performed by plant proteases. Recent indications that the 26S proteasome, a T1 family-threonine protease, is a regulator of light and hormone responsive signal transduction highlight the potential of proteases to participate in many aspects of plant growth and development. Recent discoveries that proteases are required for stomatal distribution, embryo development and disease resistance point to wider roles for four additional multigene families that include some of the most frequently studied (yet poorly understood) plant proteases: the subtilisin-like, serine proteases (family S8), the papain-like, cysteine proteases (family C1A), the pepsin-like, aspartic proteases (family A1) and the plant matrixin, metalloproteases (family M10A). In this report, 54 subtilisin-like, 30 papain-like and 59 pepsin-like proteases from Arabidopsis, are compared with S8, C1A and A1 proteases known from other plant species at the functional, phylogenetic and gene structure levels. Examples of structural conservation between S8, C1A and A1 genes from rice, barley, tomato and soybean and those from Arabidopsis are noted, indicating that some common, essential plant protease roles were established before the divergence of monocots and eudicots. Numerous examples of tandem duplications of protease genes and evidence for a variety of restricted expression patterns suggest that a high degree of specialization exists among proteases within each family. We propose that comprehensive analysis of the functions of these genes in Arabidopsis will firmly establish serine, cysteine and aspartic proteases as regulators and effectors of a wide range of plant processes.

The Arabidopsis Information Resource (TAIR): a comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant

Arabidopsis thaliana, a small annual plant belonging to the mustard family, is the subject of study by an estimated 7000 researchers around the world. In addition to the large body of genetic, physiological and biochemical data gathered for this plant, it will be the first higher plant genome to be completely sequenced, with completion expected at the end of the year 2000. The sequencing effort has been coordinated by an international collaboration, the Arabidopsis Genome Initiative (AGI). The rationale for intensive investigation of Arabidopsis is that it is an excellent model for higher plants. In order to maximize use of the knowledge gained about this plant, there is a need for a comprehensive database and information retrieval and analysis system that will provide user-friendly access to Arabidopsis information. This paper describes the initial steps we have taken toward realizing these goals in a project called The Arabidopsis Information Resource (TAIR) (www.arabidopsis.org).

Plant disease resistance genes encode members of an ancient and diverse protein family within the nucleotide-binding superfamily

The nucleotide binding site (NBS) is a characteristic domain of many plant resistance gene products. An increasing number of NBS-encoding sequences are being identified through gene cloning, PCR amplification with degenerate primers, and genome sequencing projects. The NBS domain was analyzed from 14 known plant resistance genes and more than 400 homologs, representing 26 genera of monocotyledonous, dicotyle-donous and one coniferous species. Two distinct groups of diverse sequences were identified, indicating divergence during evolution and an ancient origin for these sequences. One group was comprised of sequences encoding an N-terminal domain with Toll/Interleukin-1 receptor homology (TIR), including the known resistance genes, N, M, L6, RPP1 and RPP5. Surprisingly, this group was entirely absent from monocot species in searches of both random genomic sequences and large collections of ESTs. A second group contained monocot and dicot sequences, including the known resistance genes, RPS2, RPM1, I2, Mi, Dm3, Pi-B, Xa1, RPP8, RPS5 and Prf. Amino acid signatures in the conserved motifs comprising the NBS domain clearly distinguished these two groups. The Arabidopsis genome is estimated to contain approximately 200 genes that encode related NBS motifs; TIR sequences were more abundant and outnumber non-TIR sequences threefold. The Arabidopsis NBS sequences currently in the databases are located in approximately 21 genomic clusters and 14 isolated loci. NBS-encoding sequences may be more prevalent in rice. The wide distribution of these sequences in the plant kingdom and their prevalence in the Arabidopsis and rice genomes indicate that they are ancient, diverse and common in plants. Sequence inferences suggest that these genes encode a novel class of nucleotide-binding proteins.

Plant disease resistance genes encode members of an ancient and diverse protein family within the nucleotide-binding superfamily

The nucleotide binding site (NBS) is a characteristic domain of many plant resistance gene products. An increasing number of NBS-encoding sequences are being identified through gene cloning, PCR amplification with degenerate primers, and genome sequencing projects. The NBS domain was analyzed from 14 known plant resistance genes and more than 400 homologs, representing 26 genera of monocotyledonous, dicotyle-donous and one coniferous species. Two distinct groups of diverse sequences were identified, indicating divergence during evolution and an ancient origin for these sequences. One group was comprised of sequences encoding an N-terminal domain with Toll/Interleukin-1 receptor homology (TIR), including the known resistance genes, N, M, L6, RPP1 and RPP5. Surprisingly, this group was entirely absent from monocot species in searches of both random genomic sequences and large collections of ESTs. A second group contained monocot and dicot sequences, including the known resistance genes, RPS2, RPM1, I2, Mi, Dm3, Pi-B, Xa1, RPP8, RPS5 and Prf. Amino acid signatures in the conserved motifs comprising the NBS domain clearly distinguished these two groups. The Arabidopsis genome is estimated to contain approximately 200 genes that encode related NBS motifs; TIR sequences were more abundant and outnumber non-TIR sequences threefold. The Arabidopsis NBS sequences currently in the databases are located in approximately 21 genomic clusters and 14 isolated loci. NBS-encoding sequences may be more prevalent in rice. The wide distribution of these sequences in the plant kingdom and their prevalence in the Arabidopsis and rice genomes indicate that they are ancient, diverse and common in plants. Sequence inferences suggest that these genes encode a novel class of nucleotide-binding proteins.

DNA hybridization evidence for the Australasian affinity of the American marsupial Dromiciops australis

DNA hybridization was used to compare representatives of the major groups of marsupials and a eutherian outgroup. Because of the large genetic distances separating marsupial families, trees were calculated from normalized percentages of hybridization; thermal-melting statistics, however, gave identical topologies for the well-supported clades. The most notable results were the association of the only extant microbiotheriid, Dromiciops australis, an American marsupial, with the Australasian Diprotodontia, and of both together with the Dasyuridae. Estimates of the rate of divergence among marsupial genomes suggest that the Dromiciops-Diprotodontia split occurred approximately 50 million years ago, well after the establishment of the major clades of marsupials but before deep oceanic barriers prohibited dispersal among Australia, Antarctica, and South America. Because Dromiciops is nested within an Australasian group, it seems likely that dispersal from Australia accounts for its present distribution.

DNA/DNA hybridization studies of the carnivorous marsupials. I: The intergeneric relationships of bandicoots (Marsupialia: Perameloidea)

A complete suite of comparisons among six bandicoot species and one outgroup marsupial was generated using the hydroxyapatite chromatography method of DNA/DNA hybridization; heterologous comparisons were also made with three other bandicoot taxa. Matrices of delta Tm's, delta modes, and delta T50Hs were generated and corrected for nonreciprocity, homoplasy, and, in the case of delta Tm's, normalized percent hybridization; these matrices were analyzed using the FITCH algorithm in Felsenstein's PHYLIP (version 3.1). Uncorrected and nonreciprocity-corrected matrices were also jackknifed and analyzed with FITCH to test for consistency. Finally, sample scores for delta Tm, delta mode, and delta T50H matrices were bootstrapped and then subjected to phylogenetic analysis. These manipulations were carried out, in part, to address criticisms of the statistics used to summarize DNA/DNA hybridization (especially T50H) and the method itself. However, with the exception of an unresolved trichotomy among the two Echymipera species and Peroryctes longicauda, all trees showed the same branchpoints. Except in the case of the tree generated from reciprocal-corrected delta Tm data, nodes were stable under jackknifing; and, again excepting the above-mentioned trichotomy, all nodes were supported by 95% or more of the bootstrapped trees. These results suggest that, despite arguments to the contrary, all three summary statistics can be valid for DNA/DNA hybridization data. Of taxonomic interest is the placement of Echymipera spp. and Peroryctes longicauda together and separate from the more distant Peroryctes raffrayanus; the genus Peroryctes is thus at least paraphyletic. The trees further grouped Echymipera-plus-Peroryctes as the sister group of Isoodon-plus-Perameles. Limited hybridizations with Macrotis lagotis suggest that its current position as representative of an entirely distinct family of perameloids is correct.