Conservation genomics and systematics of a near-extinct island radiation.

Conservation benefits from incorporating genomics to explore the impacts of population declines, inbreeding, loss of genetic variation and hybridization. Here we use the near-extinct Mariana Islands reedwarbler radiation to showcase how ancient DNA approaches can allow insights into the population dynamics of extinct species and threatened populations for which historical museum specimens or material with low DNA yield (e.g., scats, feathers) are the only sources for DNA. Despite their having paraphyletic mitochondrial DNA (mtDNA), nuclear single nucleotide polymorphisms (SNPs) support the distinctiveness of critically endangered Acrocephalus hiwae and the other three species in the radiation that went extinct between the 1960s and 1990s. Two extinct species, A. yamashinae and A. luscinius, were deeply divergent from each other and from a third less differentiated lineage containing A. hiwae and extinct A. nijoi. Both mtDNA and SNPs suggest that the two isolated populations of A. hiwae from Saipan and Alamagan Islands are sufficiently distinct to warrant subspecies recognition and separate conservation management. We detected no significant differences in genetic diversity or inbreeding between Saipan and Alamagan, nor strong signatures of geographical structuring within either island. However, the implications of possible signatures of inbreeding in both Saipan and Alamagan, and long-term population declines in A. hiwae that pre-date modern anthropogenic threats require further study with denser population sampling. Our study highlights the value that conservation genomics studies of island radiations have as windows onto the possible future for the world's biota as climate change and habitat destruction increasingly fragment their ranges and contribute to rapid declines in population abundances.

Comparative mitochondrial and chloroplast genomics of a genetically distinct form of Sargassum contributing to recent "Golden Tides" in the Western Atlantic.

Over the past 5 years, massive accumulations of holopelagic species of the brown macroalga Sargassum in coastal areas of the Caribbean have created "golden tides" that threaten local biodiversity and trigger economic losses associated with beach deterioration and impact on fisheries and tourism. In 2015, the first report identifying the cause of these extreme events implicated a rare form of the holopelagic species Sargassum natans (form VIII). However, since the first mention of S. natans VIII in the 1930s, based solely on morphological characters, no molecular data have confirmed this identification. We generated full-length mitogenomes and partial chloroplast genomes of all representative holopelagic Sargassum species, S. fluitans III and S. natans I alongside the putatively rare S. natans VIII, to demonstrate small but consistent differences between S. natans I and VIII (7 bp differences out of the 34,727). Our comparative analyses also revealed that both S. natans I and S. natans VIII share a very close phylogenetic relationship with S. fluitans III (94- and 96-bp differences of 34,727). We designed novel primers that amplified regions of the cox2 and cox3 marker genes with consistent polymorphic sites that enabled differentiation between the two S. natans forms (I and VIII) from each other and both from S. fluitans III in over 150 Sargassum samples including those from the 2014 golden tide event. Despite remarkable gene synteny and sequence conservation, the three Sargassum forms differ in morphology, ecology, and distribution patterns, warranting more extensive interrogation of holopelagic Sargassum genomes as a whole.

Initial Evidence for Adaptive Selection on the NADH Subunit Two of Freshwater Dolphins by Analyses of Mitochondrial Genomes.

A small number of cetaceans have adapted to an entirely freshwater environment, having colonized rivers in Asia and South America from an ancestral origin in the marine environment. This includes the 'river dolphins', early divergence from the odontocete lineage, and two species of true dolphins (Family Delphinidae). Successful adaptation to the freshwater environment may have required increased demands in energy involved in processes such as the mitochondrial osmotic balance. For this reason, riverine odontocetes provide a compelling natural experiment in adaptation of mammals from marine to freshwater habitats. Here we present initial evidence of positive selection in the NADH dehydrogenase subunit 2 of riverine odontocetes by analyses of full mitochondrial genomes, using tests of selection and protein structure modeling. The codon model with highest statistical support corresponds to three discrete categories for amino acid sites, those under positive, neutral, and purifying selection. With this model we found positive selection at site 297 of the NADH dehydrogenase subunit 2 (dN/dS>1.0,) leading to a substitution of an Ala or Val from the ancestral state of Thr. A phylogenetic reconstruction of 27 cetacean mitogenomes showed that an Ala substitution has evolved at least four times in cetaceans, once or more in the three 'river dolphins' (Families Pontoporidae, Lipotidae and Inidae), once in the riverine Sotalia fluviatilis (but not in its marine sister taxa), once in the riverine Orcaella brevirostris from the Mekong River (but not in its marine sister taxa) and once in two other related marine dolphins. We located the position of this amino acid substitution in an alpha-helix channel in the trans-membrane domain in both the E. coli structure and Sotalia fluviatilis model. In E. coli this position is located in a helix implicated in a proton translocation channel of respiratory complex 1 and may have a similar role in the NADH dehydrogenases of cetaceans.

Low diversity in the mitogenome of sperm whales revealed by next-generation sequencing.

Large population sizes and global distributions generally associate with high mitochondrial DNA control region (CR) diversity. The sperm whale (Physeter macrocephalus) is an exception, showing low CR diversity relative to other cetaceans; however, diversity levels throughout the remainder of the sperm whale mitogenome are unknown. We sequenced 20 mitogenomes from 17 sperm whales representative of worldwide diversity using Next Generation Sequencing (NGS) technologies (Illumina GAIIx, Roche 454 GS Junior). Resequencing of three individuals with both NGS platforms and partial Sanger sequencing showed low discrepancy rates (454-Illumina: 0.0071%; Sanger-Illumina: 0.0034%; and Sanger-454: 0.0023%) confirming suitability of both NGS platforms for investigating low mitogenomic diversity. Using the 17 sperm whale mitogenomes in a phylogenetic reconstruction with 41 other species, including 11 new dolphin mitogenomes, we tested two hypotheses for the low CR diversity. First, the hypothesis that CR-specific constraints have reduced diversity solely in the CR was rejected as diversity was low throughout the mitogenome, not just in the CR (overall diversity π = 0.096%; protein-coding 3rd codon = 0.22%; CR = 0.35%), and CR phylogenetic signal was congruent with protein-coding regions. Second, the hypothesis that slow substitution rates reduced diversity throughout the sperm whale mitogenome was rejected as sperm whales had significantly higher rates of CR evolution and no evidence of slow coding region evolution relative to other cetaceans. The estimated time to most recent common ancestor for sperm whale mitogenomes was 72,800 to 137,400 years ago (95% highest probability density interval), consistent with previous hypotheses of a bottleneck or selective sweep as likely causes of low mitogenome diversity.

A spatially-explicit model of acarological risk of exposure to Borrelia burgdorferi-infected Ixodes pacificus nymphs in northwestern California based on woodland type, temperature, and water vapor.

In the far-western United States, the nymphal stage of the western black-legged tick, Ixodes pacificus, has been implicated as the primary vector to humans of Borrelia burgdorferi sensu stricto (hereinafter referred to as B. burgdorferi), the causative agent of Lyme borreliosis in North America. In the present study, we sought to determine if infection prevalence with B. burgdorferi in I. pacificus nymphs and the density of infected nymphs differ between dense-woodland types within Mendocino County, California, and to develop and evaluate a spatially-explicit model for density of infected nymphs in dense woodlands within this high-incidence area for Lyme borreliosis. In total, 4.9% (264) of 5431 I. pacificus nymphs tested for the presence of B. burgdorferi were infected. Among the 78 sampling sites, infection prevalence ranged from 0 to 22% and density of infected nymphs from 0 to 2.04 per 100 m(2). Infection prevalence was highest in woodlands dominated by hardwoods (6.2%) and lowest for redwood (1.9%) and coastal pine (0%). Density of infected nymphs also was higher in hardwood-dominated woodlands than in conifer-dominated ones that included redwood or pine. Our spatial risk model, which yielded an overall accuracy of 85%, indicated that warmer areas with less variation between maximum and minimum monthly water vapor in the air were more likely to include woodlands with elevated acarological risk of exposure to infected nymphs. We found that 37% of dense woodlands in the county were predicted to pose an elevated risk of exposure to infected nymphs, and that 94% of the dense-woodland areas that were predicted to harbor elevated densities of infected nymphs were located on privately-owned land.

Multiple diverse circoviruses infect farm animals and are commonly found in human and chimpanzee feces.

Circoviruses are known to infect birds and pigs and can cause a wide range of severe symptoms with significant economic impact. Using viral metagenomics, we identified circovirus-like DNA sequences and characterized 15 circular viral DNA genomes in stool samples from humans in Pakistan, Nigeria, Tunisia, and the United States and from wild chimpanzees. Distinct genomic features and phylogenetic analysis indicate that some viral genomes were part of a previously unrecognized genus in the Circoviridae family we tentatively named "Cyclovirus" whose genetic diversity is comparable to that of all the known species in the Circovirus genus. Circoviridae detection in the stools of U.S. adults was limited to porcine circoviruses which were also found in most U.S. pork products. To determine whether the divergent cycloviruses found in non-U.S. human stools were of dietary origin, we genetically compared them to the cycloviruses in muscle tissue samples of commonly eaten farm animals in Pakistan and Nigeria. Limited genetic overlap between cycloviruses in human stool samples and local cow, goat, sheep, camel, and chicken meat samples indicated that the majority of the 25 Cyclovirus species identified might be human viruses. We show that the genetic diversity of small circular DNA viral genomes in various mammals, including humans, is significantly larger than previously recognized, and frequent exposure through meat consumption and contact with animal or human feces provides ample opportunities for cyclovirus transmission. Determining the role of cycloviruses, found in 7 to 17% of non-U.S. human stools and 3 to 55% of non-U.S. meat samples tested, in both human and animal diseases is now facilitated by knowledge of their genomes.

Metagenomic analyses of viruses in stool samples from children with acute flaccid paralysis.

We analyzed viral nucleic acids in stool samples collected from 35 South Asian children with nonpolio acute flaccid paralysis (AFP). Sequence-independent reverse transcription and PCR amplification of capsid-protected, nuclease-resistant viral nucleic acids were followed by DNA sequencing and sequence similarity searches. Limited Sanger sequencing (35 to 240 subclones per sample) identified an average of 1.4 distinct eukaryotic viruses per sample, while pyrosequencing yielded 2.6 viruses per sample. In addition to bacteriophage and plant viruses, we detected known enteric viruses, including rotavirus, adenovirus, picobirnavirus, and human enterovirus species A (HEV-A) to HEV-C, as well as numerous other members of the Picornaviridae family, including parechovirus, Aichi virus, rhinovirus, and human cardiovirus. The viruses with the most divergent sequences relative to those of previously reported viruses included members of a novel Picornaviridae genus and four new viral species (members of the Dicistroviridae, Nodaviridae, and Circoviridae families and the Bocavirus genus). Samples from six healthy contacts of AFP patients were similarly analyzed and also contained numerous viruses, particularly HEV-C, including a potentially novel Enterovirus genotype. Determining the prevalences and pathogenicities of the novel genotypes, species, genera, and potential new viral families identified in this study in different demographic groups will require further studies with different demographic and patient groups, now facilitated by knowledge of these viral genomes.

Phylogeny of swallows (Aves: Hirundinidae) estimated from nuclear and mitochondrial DNA sequences.

The phylogeny of swallows was reconstructed by comparing segments of three genes, nuclear beta-fibrinogen intron 7 (betafib7), mitochondrial cytochrome b (cytb), and mitochondrial ND2, in a variety of combinations using maximum likelihood and Bayesian methods. betafib7 was sequenced for 47 species, cytb for 74 species, and ND2 for 61 species to yield comparisons among 75 of the 84 currently recognized swallow species. The family Hirundinidae was confirmed to consist of two clades, Pseudochelidoninae (river martins) and Hirundininae (typical swallows). The Hirundininae is further divided into mud nesters (Hirundo sensu lato), core martins (Phedina, Riparia, and New World endemic genera), and basal relicts (Psalidoprocne, Cheramoeca, and Pseudhirundo). We did not resolve the hierarchy among these three hirundinine groups, but discovered many relationships within them. Mud-nesting genera have the following relationships: (Hirundo sensu stricto, Ptyonoprogne), (Delichon, (Petrochelidon, Cecropis)). Core martins have the following topology: (Phedina, Riparia cincta), (Riparia sensu stricto, Tachycineta, ((Stelgidopteryx, Progne), (Neotropical endemic genera))). Interspecific relationships among the Neotropical endemics were resolved completely; Atticora and Notiochelidon are paraphyletic, and all Neotropical endemics probably should be lumped into one or two genera. The final group of hirundinines, the basal relicts, consists of a sister pair, the Australian Cheramoeca and African Pseudhirundo. The African saw-wings (Psalidoprocne) are their likely sister group.

Phylogeny of the tree swallow genus, Tachycineta (Aves: Hirundinidae), by Bayesian analysis of mitochondrial DNA sequences.

To set the stage for historical analyses of the ecology and behavior of tree swallows and their allies (genus Tachycineta), we reconstructed the phylogeny of the nine Tachycineta species by comparing DNA sequences of six mitochondrial genes: Cytochrome b (990 base pairs), the second subunit of nicotinamide adenine dinucleotide dehydrogenase (839 base pairs), cytochrome oxidase II (85 base pairs), ATPase 8 (158 base pairs), tRNA-lysine (73 base pairs), and tRNA-methionine (25 base pairs). The phylogeny consisted of two main clades: South and Central American species ((T. stolzmanni, T. albilinea, T. albiventris), (T. leucorrhoa, T. meyeni)), and North American and Caribbean species (T. bicolor, (T. thalassina, T. euchrysea, T. cyaneoviridis)). The genetic distances among the species suggested that Tachycineta is a relatively old group compared to other New World swallow genera. One interesting biogeographic discovery was the close relationship between Caribbean and western North American taxa. This historical connection occurs in other groups of swallows and swifts as well. To reconstruct the phylogeny, we employed Bayesian as well as traditional maximum-likelihood methods. The Bayesian approach provided probability values for trees produced from the different genes and gene combinations, as well as probabilities of branches within those trees. We compared Bayesian and maximum-likelihood bootstrap branch support and found that all branches with Bayesian probabilities > or = 95% received bootstrap support >70%.

RECOGNIZING AND TESTING HOMOLOGY OF COURTSHIP DISPLAYS IN STORKS (AVES: CICONIIFORMES: CICONIIDAE).

Ethological studies in the 1940s and 1950s, most notably those of Lorenz and Tinbergen, emphasized a historical perspective. By the 1970s, the notion that behavioral traits are too plastic to retain historical information became prevalent, and evolutionary approaches in behavioral studies were largely abandoned. However, several recent studies have demonstrated that behavioral characters are remarkably consistent with phylogenies obtained from other data and not particularly prone to homoplasy. In this study, I coded descriptions of courtship display behaviors in stork species (Aves: Ciconiiformes: Ciconiidae) as a matrix of discrete characters. I mapped each behavioral character onto a phylogeny based on DNA-DNA hybridization distances to test the homology of individual characters. Generally, displays occurring early in courtship were congruent with phylogenetic relationships and showed little homoplasy, while displays occurring late in courtship were more homoplastic. I also performed a phylogenetic analysis of the behavioral data matrix using maximum parsimony. The strict consensus of the 24 most-parsimonious trees was congruent with the DNA-DNA hybridization tree in all nodes having greater than 70% bootstrap support.