A method for determining the position and size of optimal sequence regions for phylogenetic analysis

Best practice for wildlife gut microbiome research: A comprehensive review of methodology for 16S rRNA gene investigations

Extensive research in well-studied animal models underscores the importance of commensal gastrointestinal (gut) microbes to animal physiology. Gut microbes have been shown to impact dietary digestion, mediate infection, and even modify behavior and cognition. Given the large physiological and pathophysiological contribution microbes provide their host, it is reasonable to assume that the vertebrate gut microbiome may also impact the fitness, health and ecology of wildlife. In accordance with this expectation, an increasing number of investigations have considered the role of the gut microbiome in wildlife ecology, health, and conservation. To help promote the development of this nascent field, we need to dissolve the technical barriers prohibitive to performing wildlife microbiome research. The present review discusses the 16S rRNA gene microbiome research landscape, clarifying best practices in microbiome data generation and analysis, with particular emphasis on unique situations that arise during wildlife investigations. Special consideration is given to topics relevant for microbiome wildlife research from sample collection to molecular techniques for data generation, to data analysis strategies. Our hope is that this article not only calls for greater integration of microbiome analyses into wildlife ecology and health studies but provides researchers with the technical framework needed to successfully conduct such investigations.

DNA extraction from museum specimens of parasitic Hymenoptera

At the same time that molecular researchers are improving techniques to extract DNA from museum specimens, this increased demand for access to museum specimens has created tension between the need to preserve specimens for maintaining collections and morphological research and the desire to conduct molecular analyses. To address these concerns, we examined the suitability of non-invasive DNA extraction techniques on three species of parasitic Hymenoptera (Braconidae), and test the effects of body size (parasitoid species), age (time since collection), and DNA concentration from each extract on the probability of amplifying meaningful fragments of two commonly used genetic loci. We found that age was a significant factor for determining the probability of success for sequencing both 28S and COI fragments. While the size of the braconid parasitoids significantly affected the total amount of extracted DNA, neither size nor DNA concentration were significant factors for the amplification of either gene region. We also tested several primer combinations of various lengths, but were unable to amplify fragments longer than ∼150 base pairs. These short fragments of 28S and COI were however sufficient for species identification, and for the discovery of within species genetic variation.

Genome-scale evidence of the nematode-arthropod clade

The most extensive phylogenetic analysis carried out to date, including 11 complete genomes, is shown to support the Ecdysozoa hypothesis in the open-ended debate of the Coelomata-Ecdysozoa evolutionary problem.

Background

The issue of whether coelomates form a single clade, the Coelomata, or whether all animals that moult an exoskeleton (such as the coelomate arthropods and the pseudocoelomate nematodes) form a distinct clade, the Ecdysozoa, is the most puzzling issue in animal systematics and a major open-ended subject in evolutionary biology. Previous single-gene and genome-scale analyses designed to resolve the issue have produced contradictory results. Here we present the first genome-scale phylogenetic evidence that strongly supports the Ecdysozoa hypothesis.

Results

Through the most extensive phylogenetic analysis carried out to date, the complete genomes of 11 eukaryotic species have been analyzed in order to find homologous sequences derived from 18 human chromosomes. Phylogenetic analysis of datasets showing an increased adjustment to equal evolutionary rates between nematode and arthropod sequences produced a gradual change from support for Coelomata to support for Ecdysozoa. Transition between topologies occurred when fast-evolving sequences of Caenorhabditis elegans were removed. When chordate, nematode and arthropod sequences were constrained to fit equal evolutionary rates, the Ecdysozoa topology was statistically accepted whereas Coelomata was rejected.

Conclusions

The reliability of a monophyletic group clustering arthropods and nematodes was unequivocally accepted in datasets where traces of the long-branch attraction effect were removed. This is the first phylogenomic evidence to strongly support the 'moulting clade' hypothesis.

Amino acid vs. nucleotide characters: challenging preconceived notions

The 567-terminal analysis of atpB, rbcL, and 18S rDNA was used as an empirical example to test the use of amino acid vs. nucleotide characters for protein-coding genes at deeper taxonomic levels. Nucleotides for atpB and rbcL had 6.5 times the amount of possible synapomorphy as amino acids. Based on parsimony analyses with unordered character states, nucleotides outperformed amino acids for all three measures of phylogenetic signal used (resolution, branch support, and congruence with independent evidence). The nucleotide tree was much more resolved than the amino acid tree, for both large and small clades. Nearly twice the percentage of well-supported clades resolved in the 18S rDNA tree were resolved using nucleotides (91.8%) relative to amino acids (49.2%). The well-supported clades resolved by both character types were much better supported by nucleotides (98.7% vs. 83.8% average jackknife support). The faster evolving nucleotides with a smaller average character-state space outperformed the slower evolving amino acids with a larger average character-state space. Nucleotides outperformed amino acids even with 90% of the terminals deleted. The lack of resolution on the amino acid trees appears to be caused by a lack of congruence among the amino acids, not a lack of replacement substitutions.

Molecular epidemiology of serotype O foot-and-mouth disease virus with emphasis on West and South Africa

Genetic relationships of serotype O foot-and-mouth disease (FMD) viruses recovered from outbreaks of the disease in the West African countries of Niger, Burkina Faso and, Ghana (1988-1993) and those from South Africa (2000) were determined by partial VP1 gene characterization. A 581-bp fragment, corresponding to the C-terminus half of the ID (VP1 gene) region was amplified and sequenced. An homologous region of 495 nucleotides was ultimately used to determine genetic relationships of serotype O viruses from the Middle East, Europe, South America, North Africa, East Africa, southern Africa and Asia. Seven distinct type O genotypes were identified by phylogenetic reconstruction, consisting of viruses from the following geographical regions: Genotype A: Asia, the Middle East, and South Africa, Genotype B: East Africa, Genotype C: West and North Africa, Genotype D: Taiwan and Russia, Genotype E: Angola and Venezuela, Genotype F: Western Europe, and Genotype G: Europe and South America. The genotypes constitute three different evolutionary lineages (I-III), which correspond to three discrete continental regions, some of which display inter-continental distributions due to introductions. Results further indicate that the outbreaks in Burkina Faso (1992) and Ghana (1993) are part of the same epizootic and that the strain involved in a recent outbreak of the disease in South Africa is most closely related (97% sequence identity) to a 1997 Bangladesh strain.

Antigenic and genetic divergence of rabies viruses from bat species indigenous to Canada

Antigenic characterisation of over 350 chiropteran rabies viruses of the Americas, especially from species reported rabid in Canada, distinguished 13 viral types. In close accord with this classification, nucleotide sequencing of representative isolates, at both the N and G loci, identified four principal phylogenetic groups (I-IV), sub-groups of which circulated in particular bat species. Amongst the North American bat viruses, there was a notable division between group I specimens associated with colonial, non-migratory bats (Myotis sp. and Eptesicus fuscus) and those of group II harbored by solitary, migratory species (Lasiurus sp. and Lasionycteris noctivagans). Certain species of Myotis were clearly identified as rabies reservoirs, an observation often obscured previously by their frequent infection by viral variants of other chiroptera. An additional group (III) apparently circulates in E. fuscus, whilst viruses harbored by both insectivorous and haematophagus bats of Latin America clustered to a separate clade (group IV). Comparison of the predicted N and G proteins of these viruses with those of strains of terrestrial mammals indicated a similarity in structural organisation regardless of host species lifestyle. Finally, these sequences permitted examination of the evolutionary relationship of American bat rabies viruses within the Lyssavirus genus.

Phylogenetic information and experimental design in molecular systematics

Despite the widespread perception that evolutionary inference from molecular sequences is a statistical problem, there has been very little attention paid to questions of experimental design. Previous consideration of this topic has led to little more than an empirical folklore regarding the choice of suitable genes for analysis, and to dispute over the best choice of taxa for inclusion in data sets. I introduce what I believe are new methods that permit the quantification of phylogenetic information in a sequence alignment. The methods use likelihood calculations based on Markov-process models of nucleotide substitution allied with phylogenetic trees, and allow a general approach to optimal experimental design. Two examples are given, illustrating realistic problems in experimental design in molecular phylogenetics and suggesting more general conclusions about the choice of genomic regions, sequence lengths and taxa for evolutionary studies.

Primary genetic characterization of HIV-1 isolates from WHO-sponsored vaccine evaluation sites by the RNase-A mismatch method

We have analyzed by the RNase-A mismatch method 35 isolates from four WHO-sponsored vaccine evaluation sites as a secondary laboratory of the WHO Network for HIV Isolation and Characterization. The application of an estimator for the establishment of genetic distances based on the RNase-A digestion patterns in combination with the phylogenetic analysis has allowed us to construct a tree with five well defined groups of viruses. Because the clustering with known reference strains, samples from Brazil could be grouped as subtype B and the majority of those from Thailand were subtype E. Some of the samples from Uganda were classified as subtype D. Isolates from Rwanda and some from Uganda were identified as subtype A viruses. These results coincide with data obtained by heteroduplex mobility assay and nucleotide sequencing in env regions. The RNase-A mismatch method combined with phylogenetic analysis permitted the primary genetic classification of 33 of 35 samples from the WHO Network.