Tracking the geographical spread of avian influenza (H5N1) with multiple phylogenetic trees

Fundamentals of genomic epidemiology, lessons learned from the coronavirus disease 2019 (COVID-19) pandemic, and new directions

The coronavirus disease 2019 (COVID-19) pandemic was one of the significant causes of death worldwide in 2020. The disease is caused by severe acute coronavirus syndrome (SARS) coronavirus 2 (SARS-CoV-2), an RNA virus of the subfamily Orthocoronavirinae related to 2 other clinically relevant coronaviruses, SARS-CoV and MERS-CoV. Like other coronaviruses and several other viruses, SARS-CoV-2 originated in bats. However, unlike other coronaviruses, SARS-CoV-2 resulted in a devastating pandemic. The SARS-CoV-2 pandemic rages on due to viral evolution that leads to more transmissible and immune evasive variants. Technology such as genomic sequencing has driven the shift from syndromic to molecular epidemiology and promises better understanding of variants. The COVID-19 pandemic has exposed critical impediments that must be addressed to develop the science of pandemics. Much of the progress is being applied in the developed world. However, barriers to the use of molecular epidemiology in low- and middle-income countries (LMICs) remain, including lack of logistics for equipment and reagents and lack of training in analysis. We review the molecular epidemiology literature to understand its origins from the SARS epidemic (2002-2003) through influenza events and the current COVID-19 pandemic. We advocate for improved genomic surveillance of SARS-CoV and understanding the pathogen diversity in potential zoonotic hosts. This work will require training in phylogenetic and high-performance computing to improve analyses of the origin and spread of pathogens. The overarching goals are to understand and abate zoonosis risk through interdisciplinary collaboration and lowering logistical barriers.

Phylogenetic Concepts and Tools Applied to Epidemiologic Investigations of Infectious Diseases

In this review, which is a part of the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology of Infectious Diseases, I present an overview of the principles used to classify organisms in the field of phylogenetics, highlight the methods used to infer the interrelationships of organisms, and summarize how these concepts are applied to molecular epidemiologic analyses. I present steps in analyses that come downstream of the assembly of a set of genomes or genes and the production of a multiple-sequence alignment or other matrices of putative orthologs for comparison. I focus on the history of the problem of phylogenetic reconstruction and debates within the field about the most appropriate methods. I illustrate methods that bridge the gap between molecular epidemiology and traditional epidemiology, including phylogenetic character evolution and geographic visualization. Finally, I provide practical advice on how to conduct an example analysis in the appendix. *This article is part of a curated collection.

A high-precision rule-based extraction system for expanding geospatial metadata in GenBank records

OBJECTIVE

The metadata reflecting the location of the infected host (LOIH) of virus sequences in GenBank often lacks specificity. This work seeks to enhance this metadata by extracting more specific geographic information from related full-text articles and mapping them to their latitude/longitudes using knowledge derived from external geographical databases.

MATERIALS AND METHODS

We developed a rule-based information extraction framework for linking GenBank records to the latitude/longitudes of the LOIH. Our system first extracts existing geospatial metadata from GenBank records and attempts to improve it by seeking additional, relevant geographic information from text and tables in related full-text PubMed Central articles. The final extracted locations of the records, based on data assimilated from these sources, are then disambiguated and mapped to their respective geo-coordinates. We evaluated our approach on a manually annotated dataset comprising of 5728 GenBank records for the influenza A virus.

RESULTS

We found the precision, recall, and f-measure of our system for linking GenBank records to the latitude/longitudes of their LOIH to be 0.832, 0.967, and 0.894, respectively.

DISCUSSION

Our system had a high level of accuracy for linking GenBank records to the geo-coordinates of the LOIH. However, it can be further improved by expanding our database of geospatial data, incorporating spell correction, and enhancing the rules used for extraction.

CONCLUSION

Our system performs reasonably well for linking GenBank records for the influenza A virus to the geo-coordinates of their LOIH based on record metadata and information extracted from related full-text articles.

Anchored enrichment dataset for true flies (order Diptera) reveals insights into the phylogeny of flower flies (family Syrphidae)

BACKGROUND

Anchored hybrid enrichment is a form of next-generation sequencing that uses oligonucleotide probes to target conserved regions of the genome flanked by less conserved regions in order to acquire data useful for phylogenetic inference from a broad range of taxa. Once a probe kit is developed, anchored hybrid enrichment is superior to traditional PCR-based Sanger sequencing in terms of both the amount of genomic data that can be recovered and effective cost. Due to their incredibly diverse nature, importance as pollinators, and historical instability with regard to subfamilial and tribal classification, Syrphidae (flower flies or hoverflies) are an ideal candidate for anchored hybrid enrichment-based phylogenetics, especially since recent molecular phylogenies of the syrphids using only a few markers have resulted in highly unresolved topologies. Over 6200 syrphids are currently known and uncovering their phylogeny will help us to understand how these species have diversified, providing insight into an array of ecological processes, from the development of adult mimicry, the origin of adult migration, to pollination patterns and the evolution of larval resource utilization.

RESULTS

We present the first use of anchored hybrid enrichment in insect phylogenetics on a dataset containing 30 flower fly species from across all four subfamilies and 11 tribes out of 15. To produce a phylogenetic hypothesis, 559 loci were sampled to produce a final dataset containing 217,702 sites. We recovered a well resolved topology with bootstrap support values that were almost universally >95 %. The subfamily Eristalinae is recovered as paraphyletic, with the strongest support for this hypothesis to date. The ant predators in the Microdontinae are sister to all other syrphids. Syrphinae and Pipizinae are monophyletic and sister to each other. Larval predation on soft-bodied hemipterans evolved only once in this family.

CONCLUSIONS

Anchored hybrid enrichment was successful in producing a robustly supported phylogenetic hypothesis for the syrphids. Subfamilial reconstruction is concordant with recent phylogenetic hypotheses, but with much higher support values. With the newly designed probe kit this analysis could be rapidly expanded with further sampling, opening the door to more comprehensive analyses targeting problem areas in syrphid phylogenetics and ecology.

Closely-related taxa influence woody species discrimination via DNA barcoding: evidence from global forest dynamics plots

To determine how well DNA barcodes from the chloroplast region perform in forest dynamics plots (FDPs) from global CTFS-ForestGEO network, we analyzed DNA barcoding sequences of 1277 plant species from a wide phylogenetic range (3 FDPs in tropics, 5 in subtropics and 5 in temperate zone) and compared the rates of species discrimination (RSD). We quantified RSD by two DNA barcode combinations (rbcL + matK and rbcL + matK + trnH-psbA) using a monophyly-based method (GARLI). We defined two indexes of closely-related taxa (Gm/Gt and S/G ratios) and correlated these ratios with RSD. The combination of rbcL + matK averagely discriminated 88.65%, 83.84% and 72.51% at the local, regional and global scales, respectively. An additional locus trnH-psbA increased RSD by 2.87%, 1.49% and 3.58% correspondingly. RSD varied along a latitudinal gradient and were negatively correlated with ratios of closely-related taxa. Successes of species discrimination generally depend on scales in global FDPs. We suggested that the combination of rbcL + matK + trnH-psbA is currently applicable for DNA barcoding-based phylogenetic studies on forest communities.

Predicting the emergence of H3N2 influenza viruses reveals contrasted modes of evolution of HA and NA antigens

Vaccine design for rapidly changing viruses is based on empirical surveillance of strains circulating in a given season to assess those that will most likely spread during the next season. The choice of which strains to include in the vaccine is critical, as an erroneous decision can lead to a nonimmunized human population that will then be at risk in the face of an epidemic or, worse, a pandemic. Here, we present the first steps toward a very general phylogenetic approach to predict the emergence of novel viruses. Our genomic model builds upon natural features of viral evolution such as selection and recombination / reassortment, and incorporates episodic bursts of evolution and or of recombination. As a proof-of-concept, we assess the performance of this model in a retrospective study, focusing: (i) on the emergence of an unexpected H3N2 influenza strain in 2007, and (ii) on a longitudinal design. Based on the analysis of hemagglutinin (HA) and neuraminidase (NA) genes, our results show a lack of predictive power in both experimental designs, but shed light on the mode of evolution of these two antigens: (i) supporting the lack of significance of recombination in the evolution of this influenza virus, and (ii) showing that HA evolves episodically while NA changes gradually.

Insect mitochondrial genomics: implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

Transposable element evolution in Heliconius suggests genome diversity within Lepidoptera

BACKGROUND

Transposable elements (TEs) have the potential to impact genome structure, function and evolution in profound ways. In order to understand the contribution of transposable elements (TEs) to Heliconius melpomene, we queried the H. melpomene draft sequence to identify repetitive sequences.

RESULTS

We determined that TEs comprise ~25% of the genome. The predominant class of TEs (~12% of the genome) was the non-long terminal repeat (non-LTR) retrotransposons, including a novel SINE family. However, this was only slightly higher than content derived from DNA transposons, which are diverse, with several families having mobilized in the recent past. Compared to the only other well-studied lepidopteran genome, Bombyx mori, H. melpomene exhibits a higher DNA transposon content and a distinct repertoire of retrotransposons. We also found that H. melpomene exhibits a high rate of TE turnover with few older elements accumulating in the genome.

CONCLUSIONS

Our analysis represents the first complete, de novo characterization of TE content in a butterfly genome and suggests that, while TEs are able to invade and multiply, TEs have an overall deleterious effect and/or that maintaining a small genome is advantageous. Our results also hint that analysis of additional lepidopteran genomes will reveal substantial TE diversity within the group.

Sperm Cells of a Primitive Strepsipteran

The unusual life style of Strepsiptera has presented a long-standing puzzle in establishing its affinity to other insects. Although Strepsiptera share few structural similarities with other insect orders, all members of this order share a parasitic life style with members of two distinctive families in the Coleoptera-the order now considered the most closely related to Strepsiptera based on recent genomic evidence. Among the structural features of several strepsipteran families and other insect families that have been surveyed are the organization of testes and ultrastructure of sperm cells. For comparison with existing information on insect sperm structure, this manuscript presents a description of testes and sperm of a representative of the most primitive extant strepsipteran family Mengenillidae, Eoxenos laboulbenei. We compare sperm structure of E. laboulbenei from this family with that of the three other families of Strepsiptera in the other strepsipteran suborder Stylopidia that have been studied as well as with members of the beetle families Meloidae and Rhipiphoridae that share similar life histories with Strepsiptera. Meloids, Rhipiphorids and Strepsipterans all begin larval life as active and viviparous first instar larvae. This study examines global features of these insects' sperm cells along with specific ultrastructural features of their organelles.

No evidence for directional evolution of body mass in herbivorous theropod dinosaurs

The correlation between large body size and digestive efficiency has been hypothesized to have driven trends of increasing mass in herbivorous clades by means of directional selection. Yet, to date, few studies have investigated this relationship from a phylogenetic perspective, and none, to our knowledge, with regard to trophic shifts. Here, we reconstruct body mass in the three major subclades of non-avian theropod dinosaurs whose ecomorphology is correlated with extrinsic evidence of at least facultative herbivory in the fossil record--all of which also achieve relative gigantism (more than 3000 kg). Ordinary least-squares regressions on natural log-transformed mean mass recover significant correlations between increasing mass and geological time. However, tests for directional evolution in body mass find no support for a phylogenetic trend, instead favouring passive models of trait evolution. Cross-correlation of sympatric taxa from five localities in Asia reveals that environmental influences such as differential habitat sampling and/or taphonomic filtering affect the preserved record of dinosaurian body mass in the Cretaceous. Our results are congruent with studies documenting that behavioural and/or ecological factors may mitigate the benefit of increasing mass in extant taxa, and suggest that the hypothesis can be extrapolated to herbivorous lineages across geological time scales.

Identification of cortex herbs using the DNA barcode nrITS2

To discriminate the cortex herbs of Chinese Pharmacopoeia, the second internal transcribed spacer (ITS2) of ribosomal DNA of 51 samples belonging to 19 kinds of cortex herbs were analyzed in this paper. Sequence assembly was performed using the CodonCode Aligner. Phylogenetic study was performed using the molecular evolutionary genetics analysis (MEGA) 4.0 software in accordance with the Kimura 2-Parameter (K2P) model. Phylogenetic tree was constructed using the neighbor-joining (NJ) method. We found that the ITS2 sequences of the studied samples ranged from 207 to 256 bp in length and easy to be amplified. The intraspecific genetic distance of the cortex herbs was between 0 and 0.073, which was lower than their interspecific genetic distance (mean, 0.868; minimum, 0.101). In the cluster dendrogram, all the studied medicinal materials with several samples, were monophyletic. It is concluded that ITS2 barcode is suitable for the identification of cortex herbs of Chinese Pharmacopoeia, and it will play an important role in the field of identification of traditional Chinese medicine (TCM).

Alignment free characterization of the influenza-A hemagglutinin genes by the ISSCOR method

Analyses and visualizations by the ISSCOR method of the influenza virus hemagglutinin genes of three different A-subtypes revealed some rather striking temporal (for A/H3N3), and spatial relationships (for A/H5N1) between groups of individual gene subsets. The application to the A/H1N1 set revealed also relationships between the seasonal H1, and the swine-like novel 2009 H1v variants in a quick and unambiguous manner. Based on these examples we consider the application of the ISSCOR method for analysis of large sets of homologous genes as a worthwhile addition to a toolbox of genomics-it allows a rapid diagnostics of trends, and possibly can even aid an early warning of newly emerging epidemiological threats.

Eco-virological approach for assessing the role of wild birds in the spread of avian influenza H5N1 along the Central Asian Flyway

A unique pattern of highly pathogenic avian influenza (HPAI) H5N1 outbreaks has emerged along the Central Asia Flyway, where infection of wild birds has been reported with steady frequency since 2005. We assessed the potential for two hosts of HPAI H5N1, the bar-headed goose (Anser indicus) and ruddy shelduck (Tadorna tadorna), to act as agents for virus dispersal along this 'thoroughfare'. We used an eco-virological approach to compare the migration of 141 birds marked with GPS satellite transmitters during 2005-2010 with: 1) the spatio-temporal patterns of poultry and wild bird outbreaks of HPAI H5N1, and 2) the trajectory of the virus in the outbreak region based on phylogeographic mapping. We found that biweekly utilization distributions (UDs) for 19.2% of bar-headed geese and 46.2% of ruddy shelduck were significantly associated with outbreaks. Ruddy shelduck showed highest correlation with poultry outbreaks owing to their wintering distribution in South Asia, where there is considerable opportunity for HPAI H5N1 spillover from poultry. Both species showed correlation with wild bird outbreaks during the spring migration, suggesting they may be involved in the northward movement of the virus. However, phylogeographic mapping of HPAI H5N1 clades 2.2 and 2.3 did not support dissemination of the virus in a northern direction along the migration corridor. In particular, two subclades (2.2.1 and 2.3.2) moved in a strictly southern direction in contrast to our spatio-temporal analysis of bird migration. Our attempt to reconcile the disciplines of wild bird ecology and HPAI H5N1 virology highlights prospects offered by both approaches as well as their limitations.

Advances in insect phylogeny at the dawn of the postgenomic era

Most species on Earth are insects and thus, understanding their evolutionary relationships is key to understanding the evolution of life. Insect relationships are increasingly well supported, due largely to technological advances in molecular sequencing and phylogenetic computational analysis. In this postgenomic era, insect systematics will be furthered best by integrative methods aimed at hypothesis corroboration from molecular, morphological, and paleontological evidence. This review of the current consensus of insect relationships provides a foundation for comparative study and offers a framework to evaluate incoming genomic evidence. Notable recent phylogenetic successes include the resolution of Holometabola, including the identification of the enigmatic Strepsiptera as a beetle relative and the early divergence of Hymenoptera; the recognition of hexapods as a crustacean lineage within Pancrustacea; and the elucidation of Dictyoptera orders, with termites placed as social cockroaches. Regions of the tree that require further investigation include the earliest winged insects (Palaeoptera) and Polyneoptera (orthopteroid lineages).

Philosophy and evolution: minding the gap between evolutionary patterns and tree-like patterns

Ever since Darwin, the familiar genealogical pattern known as the Tree of Life (TOL) has been prominent in evolutionary thinking and has dominated not only systematics, but also the analysis of the units of evolution. However, recent findings indicate that the evolution of DNA, especially in prokaryotes and such DNA vehicles as viruses and plasmids, does not follow a unique tree-like pattern. Because evolutionary patterns track a greater range of processes than those captured in genealogies, genealogical patterns are in fact only a subset of a broader set of evolutionary patterns. This fact suggests that evolutionists who focus exclusively on genealogical patterns are blocked from providing a significant range of genuine evolutionary explanations. Consequently, we highlight challenges to tree-based approaches, and point the way toward more appropriate methods to study evolution (although we do not present them in technical detail). We argue that there is significant benefit in adopting wider range of models, evolutionary representations, and evolutionary explanations, based on an analysis of the full range of evolutionary processes. We introduce an ecosystem orientation into evolutionary thinking that highlights the importance of "type 1 coalitions" (functionally related units with genetic exchanges, aka "friends with genetic benefits"), "type 2 coalitions" (functionally related units without genetic exchanges), "communal interactions," and "emergent evolutionary properties." On this basis, we seek to promote the study of (especially prokaryotic) evolution with dynamic evolutionary networks, which are less constrained than the TOL, and to provide new ways to analyze an expanded range of evolutionary units (genetic modules, recombined genes, plasmids, phages and prokaryotic genomes, pangenomes, microbial communities) and evolutionary processes. Finally, we discuss some of the conceptual and practical questions raised by such network-based representation.

What is the phylogenetic signal limit from mitogenomes? The reconciliation between mitochondrial and nuclear data in the Insecta class phylogeny

BACKGROUND

Efforts to solve higher-level evolutionary relationships within the class Insecta by using mitochondrial genomic data are hindered due to fast sequence evolution of several groups, most notably Hymenoptera, Strepsiptera, Phthiraptera, Hemiptera and Thysanoptera. Accelerated rates of substitution on their sequences have been shown to have negative consequences in phylogenetic inference. In this study, we tested several methodological approaches to recover phylogenetic signal from whole mitochondrial genomes. As a model, we used two classical problems in insect phylogenetics: The relationships within Paraneoptera and within Holometabola. Moreover, we assessed the mitochondrial phylogenetic signal limits in the deeper Eumetabola dataset, and we studied the contribution of individual genes.

RESULTS

Long-branch attraction (LBA) artefacts were detected in all the datasets. Methods using Bayesian inference outperformed maximum likelihood approaches, and LBA was avoided in Paraneoptera and Holometabola when using protein sequences and the site-heterogeneous mixture model CAT. The better performance of this method was evidenced by resulting topologies matching generally accepted hypotheses based on nuclear and/or morphological data, and was confirmed by cross-validation and simulation analyses. Using the CAT model, the order Strepsiptera was recovered as sister to Coleoptera for the first time using mitochondrial sequences, in agreement with recent results based on large nuclear and morphological datasets. Also the Hymenoptera-Mecopterida association was obtained, leaving Coleoptera and Strepsiptera as the basal groups of the holometabolan insects, which coincides with one of the two main competing hypotheses. For the Paraneroptera, the currently accepted non-monophyly of Homoptera was documented as a phylogenetic novelty for mitochondrial data. However, results were not satisfactory when exploring the entire Eumetabola, revealing the limits of the phylogenetic signal that can be extracted from Insecta mitogenomes. Based on the combined use of the five best topology-performing genes we obtained comparable results to whole mitogenomes, highlighting the important role of data quality.

CONCLUSION

We show for the first time that mitogenomic data agrees with nuclear and morphological data for several of the most controversial insect evolutionary relationships, adding a new independent source of evidence to study relationships among insect orders. We propose that deeper divergences cannot be inferred with the current available methods due to sequence saturation and compositional bias inconsistencies. Our exploratory analysis indicates that the CAT model is the best dealing with LBA and it could be useful for other groups and datasets with similar phylogenetic difficulties.

At last, a Pennsylvanian stem-stonefly (Plecoptera) discovered

BACKGROUND

Stem-relatives of many winged insect orders have been identified among Pennsylvanian fossils (Carboniferous Period). Owing to their presumed 'basal' position in insect phylogeny, stoneflies were expected to occur at this period. However, no relative has ever been designated convincingly.

RESULTS

In this paper, we report specimens belonging to a new fossil insect species collected from the Tupo Formation (Pennsylvanian; China). The wing venation of Gulou carpenteri gen. et sp. nov. exhibits character states diagnostic of the order Plecoptera, but lack character states shared by unequivocal representatives of the order. Derived from this identification, the delimitation of the fossil species is ascertained based on comparison of several extant stonefly species. This comparative analysis allowed a trait present in G. carpenteri gen. et sp. nov., but rarely occurring in extant species, to be documented and highlighted as atavistic. Affinities of taxa formerly proposed as putative stem-stoneflies are reconsidered in the light of the new discovery.

CONCLUSIONS

Gulou carpenteri gen. et sp. nov. is considered the only genuine Plecoptera reported from the Pennsylvanian. Continuing efforts on the systematics of Pennsylvanian winged insects indicate a fauna more diverse than previously appreciated. It suggests that insects already had a long, yet undocumented, history by this time.

Telling the whole story in a 10,000-genome world

BACKGROUND

Genome sequencing has revolutionized our view of the relationships among genomes, particularly in revealing the confounding effects of lateral genetic transfer (LGT). Phylogenomic techniques have been used to construct purported trees of microbial life. Although such trees are easily interpreted and allow the use of a subset of genomes as "proxies" for the full set, LGT and other phenomena impact the positioning of different groups in genome trees, confounding and potentially invalidating attempts to construct a phylogeny-based taxonomy of microorganisms. Network and graph approaches can reveal complex sets of relationships, but applying these techniques to large data sets is a significant challenge. Notwithstanding the question of what exactly it might represent, generating and interpreting a Tree or Network of All Genomes will only be feasible if current algorithms can be improved upon.

RESULTS

Complex relationships among even the most-similar genomes demonstrate that proxy-based approaches to simplifying large sets of genomes are not alone sufficient to solve the analysis problem. A phylogenomic analysis of 1173 sequenced bacterial and archaeal genomes generated phylogenetic trees for 159,905 distinct homologous gene sets. The relationships inferred from this set can be heavily dependent on the inclusion of other taxa: for example, phyla such as Spirochaetes, Proteobacteria and Firmicutes are recovered as cohesive groups or split depending on the presence of other specific lineages. Furthermore, named groups such as Acidithiobacillus, Coprothermobacter and Brachyspira show a multitude of affiliations that are more consistent with their ecology than with small subunit ribosomal DNA-based taxonomy. Network and graph representations can illustrate the multitude of conflicting affinities, but all methods impose constraints on the input data and create challenges of construction and interpretation.

CONCLUSIONS

These complex relationships highlight the need for an inclusive approach to genomic data, and current methods with minor alterations will likely scale to allow the analysis of data sets with 10,000 or more genomes. The main challenges lie in the visualization and interpretation of genomic relationships, and the redefinition of microbial taxonomy when subsets of genomic data are so evidently in conflict with one another, and with the "canonical" molecular taxonomy.

The Supramap project: linking pathogen genomes with geography to fight emergent infectious diseases

Novel pathogens have the potential to become critical issues of national security, public health and economic welfare. As demonstrated by the response to Severe Acute Respiratory Syndrome (SARS) and influenza, genomic sequencing has become an important method for diagnosing agents of infectious disease. Despite the value of genomic sequences in characterizing novel pathogens, raw data on their own do not provide the information needed by public health officials and researchers. One must integrate knowledge of the genomes of pathogens with host biology and geography to understand the etiology of epidemics. To these ends, we have created an application called Supramap (http://supramap.osu.edu) to put information on the spread of pathogens and key mutations across time, space and various hosts into a geographic information system (GIS). To build this application, we created a web service for integrated sequence alignment and phylogenetic analysis as well as methods to describe the tree, mutations, and host shifts in Keyhole Markup Language (KML). We apply the application to 239 sequences of the polymerase basic 2 (PB2) gene of recent isolates of avian influenza (H5N1). We map a mutation, glutamic acid to lysine at position 627 in the PB2 protein (E627K), in H5N1 influenza that allows for increased replication of the virus in mammals. We use a statistical test to support the hypothesis of a correlation of E627K mutations with avian-mammalian host shifts but reject the hypothesis that lineages with E627K are moving westward. Data, instructions for use, and visualizations are included as supplemental materials at: http://supramap.osu.edu/sm/supramap/publications. © The Willi Hennig Society 2010.

Drippy Pod of White Lupine: A New Bacterial Disease Caused by a Pathovar of Brenneria quercina

Drippy pod is a unique bacterial disease of Mediterranean white lupine (Lupinus albus) that first appeared in commercial fields in Eastern Washington State in the mid-1980s. The disease is most noticeable in the field as water-soaked lesions on lupine pods that produce an abundance of whitish-colored ooze with a sticky and foamy consistency. As the disease progresses, yellowing of lupine plants occurs with ooze characteristically dripping down the infected pods and stems and solidifying. A gram-negative rod-shaped bacterium with facultative anaerobic growth was repeatedly isolated from infected lupine tissues, and subsequently confirmed by Koch's postulates to infect lupines. Physiological and biochemical tests, including the API 20E and 50CHE strip assays, showed a highly uniform phenotype for the lupine strains that was distinctive for the genus Brenneria and most closely resembled the oak pathogen Brenneria quercina. Furthermore, sequence analyses of the 16S rDNA gene and the 16S-23S intergenic region of lupine strains revealed the highest similarity (>97%) to the corresponding regions of B. quercina and less similarity to the next closest species, B. salicis. Fatty acid profiling demonstrated that lupine strains were qualitatively similar in composition to Brenneria spp., and supported placement of the drippy pod bacterium in the species B. quercina. Oak strains of B. quercina, however, did not incite drippy pod disease on lupine. Consequently, the lupine strains that cause bacterial drippy pod disease were classified as B. quercina pv. lupinicola pv. nov.

Some considerations for analyzing biodiversity using integrative metagenomics and gene networks

BACKGROUND

Improving knowledge of biodiversity will benefit conservation biology, enhance bioremediation studies, and could lead to new medical treatments. However there is no standard approach to estimate and to compare the diversity of different environments, or to study its past, and possibly, future evolution.

PRESENTATION OF THE HYPOTHESIS

We argue that there are two conditions for significant progress in the identification and quantification of biodiversity. First, integrative metagenomic studies - aiming at the simultaneous examination (or even better at the integration) of observations about the elements, functions and evolutionary processes captured by the massive sequencing of multiple markers - should be preferred over DNA barcoding projects and over metagenomic projects based on a single marker. Second, such metagenomic data should be studied with novel inclusive network-based approaches, designed to draw inferences both on the many units and on the many processes present in the environments.

TESTING THE HYPOTHESIS

We reached these conclusions through a comparison of the theoretical foundations of two molecular approaches seeking to assess biodiversity: metagenomics (mostly used on prokaryotes and protists) and DNA barcoding (mostly used on multicellular eukaryotes), and by pragmatic considerations of the issues caused by the 'species problem' in biodiversity studies.

IMPLICATIONS OF THE HYPOTHESIS

Evolutionary gene networks reduce the risk of producing biodiversity estimates with limited explanatory power, biased either by unequal rates of LGT, or difficult to interpret due to (practical) problems caused by type I and type II grey zones. Moreover, these networks would easily accommodate additional (meta)transcriptomic and (meta)proteomic data.

Genome informatics of influenza A: from data sharing to shared analytical capabilities

Emerging infectious diseases are critical issues of public health and the economic and social stability of nations. As demonstrated by the international response to the severe acute respiratory syndrome (SARS) and influenza A, rapid genomic sequencing is a crucial tool to understand diseases that occur at the interface of human and animal populations. However, our ability to make sense of sequence data lags behind our ability to acquire the data. The potential of sequence data on pathogens is not fully realized until raw data are translated into public health intelligence. Sequencing technologies have become highly mechanized. If the political will for data sharing remains strong, the frontier for progress in emerging infectious diseases will be in analysis of sequence data and translation of results into better public health science and policy. For example, applying analytical tools such as Supramap (http://supramap.osu.edu) to genomic data for pathogens, public health scientists can track specific mutations in pathogens that confer the ability to infect humans or resist drugs. The results produced by the Supramap application are compelling visualizations of pathogen lineages and features mapped into geographic information systems that can be used to test hypotheses and to follow the spread of diseases across geography and hosts and communicate the results to a wide audience.