A targeted capture approach to generating reference sequence databases for chloroplast gene regions.

Metabarcoding has improved the way we understand plants within our environment, from their ecology and conservation to invasive species management. The notion of identifying plant taxa within environmental samples relies on the ability to match unknown sequences to known reference libraries. Without comprehensive reference databases, species can go undetected or be incorrectly assigned, leading to false-positive and false-negative detections. To improve our ability to generate reference sequence databases, we developed a targeted capture approach using the OZBaits_CP V1.0 set, designed to capture chloroplast gene regions across the entirety of flowering plant diversity. We focused on generating a reference database for coastal temperate plant species given the lack of reference sequences for these taxa. Our approach was successful across all specimens with a target gene recovery rate of 92%, which was achieved in a single assay (i.e., samples were pooled), thus making this approach much faster and more efficient than standard barcoding. Further testing of this database highlighted 80% of all samples could be discriminated to family level across all gene regions with some genes achieving greater resolution than others-which was also dependent on the taxon of interest. Thus, we demonstrate the importance of generating reference sequences across multiple chloroplast gene regions as no single loci are sufficient to discriminate across all plant groups. The targeted capture approach outlined in this study provides a way forward to achieve this.

Horizontal Transposon Transfer and Its Implications for the Ancestral Ecology of Hydrophiine Snakes.

Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids ~30 Mya and have since rapidly diversified into six amphibious, ~60 marine and ~100 terrestrial species that fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route.

Bioindicator snake shows genomic signatures of natural and anthropogenic barriers to gene flow.

Urbanisation alters landscapes, introduces wildlife to novel stressors, and fragments habitats into remnant 'islands'. Within these islands, isolated wildlife populations can experience genetic drift and subsequently suffer from inbreeding depression and reduced adaptive potential. The Western tiger snake (Notechis scutatus occidentalis) is a predator of wetlands in the Swan Coastal Plain, a unique bioregion that has suffered substantial degradation through the development of the city of Perth, Western Australia. Within the urban matrix, tiger snakes now only persist in a handful of wetlands where they are known to bioaccumulate a suite of contaminants, and have recently been suggested as a relevant bioindicator of ecosystem health. Here, we used genome-wide single nucleotide polymorphism (SNP) data to explore the contemporary population genomics of seven tiger snake populations across the urban matrix. Specifically, we used population genomic structure and diversity, effective population sizes (Ne), and heterozygosity-fitness correlations to assess fitness of each population with respect to urbanisation. We found that population genomic structure was strongest across the northern and southern sides of a major river system, with the northern cluster of populations exhibiting lower heterozygosities than the southern cluster, likely due to a lack of historical gene flow. We also observed an increasing signal of inbreeding and genetic drift with increasing geographic isolation due to urbanisation. Effective population sizes (Ne) at most sites were small (< 100), with Ne appearing to reflect the area of available habitat rather than the degree of adjacent urbanisation. This suggests that ecosystem management and restoration may be the best method to buffer the further loss of genetic diversity in urban wetlands. If tiger snake populations continue to decline in urban areas, our results provide a baseline measure of genomic diversity, as well as highlighting which 'islands' of habitat are most in need of management and protection.

A new species of turtle-headed sea Snake (Emydocephalus: Elapidae) endemic to Western Australia.

We describe a new species of turtle-headed sea snake Emydocephalus orarius sp. nov. (Elapidae) from Western Australia's Coral Coast, Pilbara and Kimberley regions. Phylogenetic analysis of mitochondrial markers places the new species as the sister lineage to the two currently recognised species in Emydocephalus: E. annulatus from the Timor Sea reefs and Coral Sea, and E. ijimae from the Ryukyu Islands. Analysis of nuclear SNP data from the new species and E. annulatus from Australia and New Caledonia provides additional independent evidence of their evolutionary distinctiveness. The new taxon is usually morphologically diagnosable from its congeners using a combination of scalation and colour pattern characters, and appears to reach greater total lengths (>1 m in the new species versus typically ~80 cm in E. annulatus/E. ijimae). The new species is known largely from soft-bottomed trawl grounds, unlike E. annulatus and E.ijimae which usually inhabit coral reefs. The discovery of this new species brings the number of sea snake species endemic to Western Australia to six.

Molecular resolution to a morphological controversy: The case of North American fossil muskoxen Bootherium and Symbos.

The musk ox (Ovibos moschatus) is the only surviving member of a group of Pleistocene North American musk ox genera (Praeovibos, Ovibos, Bootherium, Euceratherium, and Soergelia) whose taxonomy is uncertain. The helmeted musk ox (Bootherium bombifrons) and the woodland musk ox (Symbos cavifrons) have been synonymised as male and female forms of a single Nearctic species found from Alaska, in the north, to Texas, in the south. However, this reclassification has not been tested using molecular data, despite the potential to use ancient DNA to examine these late Pleistocene taxa. In the present study, we sequenced mitochondrial genomes from seven subfossil musk ox specimens (originally identified as Bootherium and/or Symbos), allowing us to evaluate the identity of these muskoxen, explore their phylogeography, and estimate the timeline for their evolution. We also used nuclear genomic data to determine the sex of six of our seven samples. Ultimately, our molecular data support the synonymisation of the North American muskoxen Bootherium and Symbos.

Molecular phylogenetics supports the origin of an endemic Balearic shrew lineage (Nesiotites) coincident with the Messinian Salinity Crisis.

The red-toothed shrews (Soricinae) are the most widespread subfamily of shrews, distributed from northern South America to North America and Eurasia. Within this subfamily, the tribe Nectogalini includes the fossil species Nesiotites hidalgo recorded from the Late Pleistocene to Holocene of the Balearic Islands (Western Mediterranean). Although there is a consensus about the close relationship between the extinct red-toothed shrew genera Nesiotites and Asoriculus based on morphology, molecular data are necessary to further evaluate the phylogenetic relationships of the Balearic fossils. We obtained a near complete mitochondrial genome of N. hidalgo, allowing the first molecular phylogenetic analysis of this species. Analyses based on 15,167 bp of the mitochondrial genome placed N. hidalgo as close relative to the extant Himalayan shrew (Soriculus nigrescens), and a combined analysis using molecular and morphological data confirm that N. hidalgo and Asoriculus gibberodon are sister-taxa with S. nigrescens as the immediate outgroup. Molecular clock and divergence estimates suggest that the split between N. hidalgo and its closest living relative occurred around 6.44 Ma, which is in agreement with the previously proposed colonisation of the Balearic Islands from mainland Europe by nectogaline shrews during the Messinian Salinity Crisis (5.97-5.33 My ago). Our results highlight that it is possible to retrieve genetic data from extinct small mammals from marginal environments for DNA preservation. Additional finds from the fossil record of Soricinae from the Eurasian Late Miocene/Early Pliocene are needed to shed further light on the still confusing taxonomy and paleobiogeography of this clade.

Genetic diversity and drivers of dwarfism in extinct island emu populations.

Australia's iconic emu (Dromaius novaehollandiae novaehollandiae) is the only living representative of its genus, but fossil evidence and reports from early European explorers suggest that three island forms (at least two of which were dwarfs) became extinct during the nineteenth century. While one of these-the King Island emu-has been found to be conspecific with Australian mainland emus, little is known about how the other two forms-Kangaroo Island and Tasmanian emus-relate to the others, or even the size of Tasmanian emus. We present a comprehensive genetic and morphological analysis of Dromaius diversity, including data from one of the few definitively genuine Tasmanian emu specimens known. Our genetic analyses suggest that all the island populations represent sub-populations of mainland Dnovaehollandiae Further, the size of island emus and those on the mainland appears to scale linearly with island size but not time since isolation, suggesting that island size-and presumably concomitant limitations on resource availability-may be a more important driver of dwarfism in island emus, though its precise contribution to emu dwarfism remains to be confirmed.

Potential Causative Mutation for Melanism in Rats Identified in the Agouti Signaling Protein Gene (Asip) of the Rattus rattus Species Complex on Okinawa Island, Japan.

The occurrence of black fur, or melanism, in many mammalian species is known to be linked to DNA sequence variation in the agouti signaling protein (Asip) gene, which is a major determinant of eumelanin and pheomelanin pigments in coat color. We investigated 38 agouti (i.e., banded wildtype) and four melanistic Rattus rattus species complex (RrC) lineage II specimens from Okinawa Island, Ryukyu Islands, Japan, for genetic variation in three exons and associated flanking regions in the Asip gene. On Okinawa, a predicted loss-of-function mutation caused by a cysteine to serine amino acid change at p.124C>S (c.370T>A) in the highly conserved functional domain of Asip was found in melanistic rats, but was absent in agouti specimens, suggesting that the p.124C>S mutation is responsible for the observed melanism. Phylogeographic analysis found that Asip sequences from Okinawan RrC lineage II, including both agouti and melanistic specimens, differed from: 1) both agouti and melanistic RrC lineage I from Otaru, Hokkaido, Japan, and 2) agouti RrC lineages I and II from South Australia. This suggests the possibility of in-situ mutation of the Asip gene, either within the RrC lineage II population on Okinawa or in an unsampled RrC lineage II population with biogeographic links to Okinawa, although incomplete lineage sorting could not be ruled out.

East African origins for Madagascan chickens as indicated by mitochondrial DNA.

The colonization of Madagascar by Austronesian-speaking people during AD 50-500 represents the most westerly point of the greatest diaspora in prehistory. A range of economically important plants and animals may have accompanied the Austronesians. Domestic chickens (Gallus gallus) are found in Madagascar, but it is unclear how they arrived there. Did they accompany the initial Austronesian-speaking populations that reached Madagascar via the Indian Ocean or were they late arrivals with Arabian and African sea-farers? To address this question, we investigated the mitochondrial DNA control region diversity of modern chickens sampled from around the Indian Ocean rim (Southeast Asia, South Asia, the Arabian Peninsula, East Africa and Madagascar). In contrast to the linguistic and human genetic evidence indicating dual African and Southeast Asian ancestry of the Malagasy people, we find that chickens in Madagascar only share a common ancestor with East Africa, which together are genetically closer to South Asian chickens than to those in Southeast Asia. This suggests that the earliest expansion of Austronesian-speaking people across the Indian Ocean did not successfully introduce chickens to Madagascar. Our results further demonstrate the complexity of the translocation history of introduced domesticates in Madagascar.

Phylogeography of the antilopine wallaroo (Macropus antilopinus) across tropical northern Australia.

The distribution of antilopine wallaroo, Macropus antilopinus, is marked by a break in the species' range between Queensland and the Northern Territory, coinciding with the Carpentarian barrier. Previous work on M. antilopinus revealed limited genetic differentiation between the Northern Territory and Queensland M. antilopinus populations across this barrier. The study also identified a number of divergent lineages in the Northern Territory, but was unable to elucidate any geographic structure. Here, we re-examine these results to (1) determine phylogeographic patterns across the range of M. antilopinus and (2) infer the biogeographic barriers associated with these patterns. The tropical savannahs of northern Australia: from the Cape York Peninsula in the east, to the Kimberley in the west. We examined phylogeographic patterns in M. antilopinus using a larger number of samples and three mtDNA genes: NADH dehydrogenase subunit 2, cytochrome b, and the control region. Two datasets were generated and analyzed: (1) a subset of samples with all three mtDNA regions concatenated together and (2) all samples for just control region sequences that included samples from the previous study. Analysis included generating phylogenetic trees based on Bayesian analysis and intraspecific median-joining networks. The contemporary spatial structure of M. antilopinus mtDNA lineages revealed five shallow clades and a sixth, divergent lineage. The genetic differences that we found between Queensland and Northern Territory M. antilopinus samples confirmed the split in the geographic distribution of the species. We also found weak genetic differentiation between Northern Territory samples and those from the Kimberley region of Western Australia, possibly due to the Kimberley Plateau-Arnhem Land barrier. Within the Northern Territory, two clades appear to be parapatric in the west, while another two clades are broadly sympatric across the Northern Territory. MtDNA diversity of M. antilopinus revealed an unexpectedly complex evolutionary history involving multiple sympatric and parapatric mtDNA clades across northern Australia. These phylogeographic patterns highlight the importance of investigating genetic variation across distributions of species and integrating this information into biodiversity conservation.

Phylogeography of Eastern Grey Kangaroos, Macropus giganteus, Suggests a Mesic Refugium in Eastern Australia.

Phylogeographic studies around the world have identified refugia where fauna were able to persist during unsuitable climatic periods, particularly during times of glaciation. In Australia the effects of Pleistocene climate oscillations on rainforest taxa have been well studied but less is known about the effects on mesic-habitat fauna, such as the eastern grey kangaroo (Macropus giganteus). The eastern grey kangaroo is a large mammal that is common and widespread throughout eastern Australia, preferring dry mesic habitat, rather than rainforest. As pollen evidence suggests that the central-eastern part of Australia (southeast Queensland and northern New South Wales) experienced cycles of expansion in mesic habitat with contraction in rainforests, and vice versa during glacial and interglacial periods, respectively, we hypothesise that the distribution of the eastern grey kangaroo was affected by these climate oscillations and may have contracted to mesic habitat refugia. From 375 mitochondrial DNA control region sequences from across the distribution of eastern grey kangaroos we obtained 108 unique haplotypes. Phylogenetic analysis identified two clades in Queensland, one of which is newly identified and restricted to a small coastal region in southern Queensland north of Brisbane, known as the Sunshine Coast. The relatively limited geographic range of this genetically isolated clade suggests the possibility of a mesic habitat refugium forming during rainforest expansion during wetter climate cycles. Other potential, although less likely, reasons for the genetic isolation of the highly distinct clade include geographic barriers, separate northward expansions, and strong local adaptation.

Late pleistocene Australian marsupial DNA clarifies the affinities of extinct megafaunal kangaroos and wallabies.

Understanding the evolution of Australia's extinct marsupial megafauna has been hindered by a relatively incomplete fossil record and convergent or highly specialized morphology, which confound phylogenetic analyses. Further, the harsh Australian climate and early date of most megafaunal extinctions (39-52 ka) means that the vast majority of fossil remains are unsuitable for ancient DNA analyses. Here, we apply cross-species DNA capture to fossils from relatively high latitude, high altitude caves in Tasmania. Using low-stringency hybridization and high-throughput sequencing, we were able to retrieve mitochondrial sequences from two extinct megafaunal macropodid species. The two specimens, Simosthenurus occidentalis (giant short-faced kangaroo) and Protemnodon anak (giant wallaby), have been radiocarbon dated to 46-50 and 40-45 ka, respectively. This is significantly older than any Australian fossil that has previously yielded DNA sequence information. Processing the raw sequence data from these samples posed a bioinformatic challenge due to the poor preservation of DNA. We explored several approaches in order to maximize the signal-to-noise ratio in retained sequencing reads. Our findings demonstrate the critical importance of adopting stringent processing criteria when distant outgroups are used as references for mapping highly fragmented DNA. Based on the most stringent nucleotide data sets (879 bp for S. occidentalis and 2,383 bp for P. anak), total-evidence phylogenetic analyses confirm that macropodids consist of three primary lineages: Sthenurines such as Simosthenurus (extinct short-faced kangaroos), the macropodines (all other wallabies and kangaroos), and the enigmatic living banded hare-wallaby Lagostrophus fasciatus (Lagostrophinae). Protemnodon emerges as a close relative of Macropus (large living kangaroos), a position not supported by recent morphological phylogenetic analyses.

Using ancient DNA to study the origins and dispersal of ancestral Polynesian chickens across the Pacific.

The human colonization of Remote Oceania remains one of the great feats of exploration in history, proceeding east from Asia across the vast expanse of the Pacific Ocean. Human commensal and domesticated species were widely transported as part of this diaspora, possibly as far as South America. We sequenced mitochondrial control region DNA from 122 modern and 22 ancient chicken specimens from Polynesia and Island Southeast Asia and used these together with Bayesian modeling methods to examine the human dispersal of chickens across this area. We show that specific techniques are essential to remove contaminating modern DNA from experiments, which appear to have impacted previous studies of Pacific chickens. In contrast to previous reports, we find that all ancient specimens and a high proportion of the modern chickens possess a group of unique, closely related haplotypes found only in the Pacific. This group of haplotypes appears to represent the authentic founding mitochondrial DNA chicken lineages transported across the Pacific, and allows the early dispersal of chickens across Micronesia and Polynesia to be modeled. Importantly, chickens carrying this genetic signature persist on several Pacific islands at high frequencies, suggesting that the original Polynesian chicken lineages may still survive. No early South American chicken samples have been detected with the diagnostic Polynesian mtDNA haplotypes, arguing against reports that chickens provide evidence of Polynesian contact with pre-European South America. Two modern specimens from the Philippines carry haplotypes similar to the ancient Pacific samples, providing clues about a potential homeland for the Polynesian chicken.

Pleistocene Chinese cave hyenas and the recent Eurasian history of the spotted hyena, Crocuta crocuta.

The living hyena species (spotted, brown, striped and aardwolf) are remnants of a formerly diverse group of more than 80 fossil species, which peaked in diversity in the Late Miocene (about 7-8 Ma). The fossil history indicates an African origin, and morphological and ancient DNA data have confirmed that living spotted hyenas (Crocuta crocuta) of Africa were closely related to extinct Late Pleistocene cave hyenas from Europe and Asia. The current model used to explain the origins of Eurasian cave hyena populations invokes multiple migrations out of Africa between 3.5-0.35 Ma. We used mitochondrial DNA sequences from radiocarbon-dated Chinese Pleistocene hyena specimens to examine the origin of Asian populations, and temporally calibrate the evolutionary history of spotted hyenas. Our results support a far more recent evolutionary timescale (430-163 kya) and suggest that extinct and living spotted hyena populations originated from a widespread Eurasian population in the Late Pleistocene, which was only subsequently restricted to Africa. We developed statistical tests of the contrasting population models and their fit to the fossil record. Coalescent simulations and Bayes Factor analysis support the new radiocarbon-calibrated timescale and Eurasian origins model. The new Eurasian biogeographic scenario proposed for the hyena emphasizes the role of the vast steppe grasslands of Eurasia in contrast to models only involving Africa. The new methodology for combining genetic and geological data to test contrasting models of population history will be useful for a wide range of taxa where ancient and historic genetic data are available.