Region-Specific Genetic Diversity of Black Rats (Rattus rattus Complex) in Southeast and East Asia Shaped by Rapid Population Expansion Events.

Among the six mitochondrial DNA lineages of the black rat (Rattus rattus Complex; RrC), lineages II and IV are widespread in Southeast and East Asia. This study explored their demographic history using 17 new sequences from the Miyako Islands in the Ryukyu archipelago, together with 178 publicly available cytochrome b sequences. We defined six and two haplotype groups showing rapid population expansion signals in Lineages II and IV, respectively. The six haplotype groups of Lineage II were represented by haplotypes from 1) Myanmar/Bangladesh/Northeast India, 2) Laos, 3) Thailand, 4) Indonesia/Philippines, 5) Vietnam/southern China, and 6) the Ryukyu archipelago. These expansion times were estimated using time-dependent evolutionary rates to be 115,300 years ago (ya), 128,500 ya, 9600 ya, 10,600 ya, 7200 ya, and 1400 ya, respectively, although all had large confidence intervals. The two groups of Lineage IV were recovered from the mainland and islands of Southeast Asia with predicted expansion times of 197,000 ya and 5800 ya, respectively. These results suggest that climatic fluctuations during the last 200,000 years of the Quaternary, affected the population dynamics in subtropical areas at different times. Furthermore, the results of the younger rapid expansion events of RrC suggest the possibility of agricultural advancement and dispersal of Neolithic farmers to different areas within the mainland and islands of Southeast Asia during the Holocene. A subset of rats from the Miyako Islands were found to have the same lineage IV haplotypes as those in Southeast Asia, suggesting a recent introduction of these new lineages.

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.

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.

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.

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.

Whole-Genome Sequencing Can Identify Clinically Relevant Variants from a Single Sub-Punch of a Dried Blood Spot Specimen.

The collection of dried blood spots (DBS) facilitates newborn screening for a variety of rare, but very serious conditions in healthcare systems around the world. Sub-punches of varying sizes (1.5-6 mm) can be taken from DBS specimens to use as inputs for a range of biochemical assays. Advances in DNA sequencing workflows allow whole-genome sequencing (WGS) libraries to be generated directly from inputs such as peripheral blood, saliva, and DBS. We compared WGS metrics obtained from libraries generated directly from DBS to those generated from DNA extracted from peripheral blood, the standard input for this type of assay. We explored the flexibility of DBS as an input for WGS by altering the punch number and size as inputs to the assay. We showed that WGS libraries can be successfully generated from a variety of DBS inputs, including a single 3 mm or 6 mm diameter punch, with equivalent data quality observed across a number of key metrics of importance in the detection of gene variants. We observed no difference in the performance of DBS and peripheral-blood-extracted DNA in the detection of likely pathogenic gene variants in samples taken from individuals with cystic fibrosis or phenylketonuria. WGS can be performed directly from DBS and is a powerful method for the rapid discovery of clinically relevant, disease-causing gene variants.

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.

An automated 13.5 hour system for scalable diagnosis and acute management guidance for genetic diseases.

While many genetic diseases have effective treatments, they frequently progress rapidly to severe morbidity or mortality if those treatments are not implemented immediately. Since front-line physicians frequently lack familiarity with these diseases, timely molecular diagnosis may not improve outcomes. Herein we describe Genome-to-Treatment, an automated, virtual system for genetic disease diagnosis and acute management guidance. Diagnosis is achieved in 13.5 h by expedited whole genome sequencing, with superior analytic performance for structural and copy number variants. An expert panel adjudicated the indications, contraindications, efficacy, and evidence-of-efficacy of 9911 drug, device, dietary, and surgical interventions for 563 severe, childhood, genetic diseases. The 421 (75%) diseases and 1527 (15%) effective interventions retained are integrated with 13 genetic disease information resources and appended to diagnostic reports ( https://gtrx.radygenomiclab.com ). This system provided correct diagnoses in four retrospectively and two prospectively tested infants. The Genome-to-Treatment system facilitates optimal outcomes in children with rapidly progressive genetic diseases.

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.

European and Asian contribution to the genetic diversity of mainland South American chickens.

Chickens (Gallus gallus domesticus) from the Americas have long been recognized as descendants of European chickens, transported by early Europeans since the fifteenth century. However, in recent years, a possible pre-Columbian introduction of chickens to South America by Polynesian seafarers has also been suggested. Here, we characterize the mitochondrial control region genetic diversity of modern chicken populations from South America and compare this to a worldwide dataset in order to investigate the potential maternal genetic origin of modern-day chicken populations in South America. The genetic analysis of newly generated chicken mitochondrial control region sequences from South America showed that the majority of chickens from the continent belong to mitochondrial haplogroup E. The rest belongs to haplogroups A, B and C, albeit at very low levels. Haplogroup D, a ubiquitous mitochondrial lineage in Island Southeast Asia and on Pacific Islands is not observed in continental South America. Modern-day mainland South American chickens are, therefore, closely allied with European and Asian chickens. Furthermore, we find high levels of genetic contributions from South Asian chickens to those in Europe and South America. Our findings demonstrate that modern-day genetic diversity of mainland South American chickens appear to have clear European and Asian contributions, and less so from Island Southeast Asia and the Pacific Islands. Furthermore, there is also some indication that South Asia has more genetic contribution to European chickens than any other Asian chicken populations.

A simple ATAC-seq protocol for population epigenetics.

We describe here a protocol for the generation of sequence-ready libraries for population epigenomics studies. The protocol is a streamlined version of the Assay for transposase accessible chromatin with high-throughput sequencing (ATAC-seq) that provides a positive display of accessible, presumably euchromatic regions. The protocol is straightforward and can be used with small individuals such as daphnia and schistosome worms, and probably many other biological samples of comparable size, and it requires little molecular biology handling expertise.

A perspective for resolving the systematics of Rattus, the vertebrates with the most influence on human welfare.

The murid rodent genus Rattus Fischer 1803 contains several species that are responsible for massive loss of crops and food, extinction of other species and the spread of zoonotic diseases to humans, as well as a laboratory species used to answer important questions in physiology, immunology, pharmacology, toxicology, nutrition, behaviour and learning. Despite the well-known significant impacts of Rattus, a definitive evolutionary based systematic framework for the genus is not yet available. The past 75 years have seen more dramatic changes in membership of Rattus than in almost any other genus of mammals. In fact, the Rattus genus has been a receptacle for any generalised Old World murine that lacked morphological specialisation and at one point, has included more than 560 species and/or subspecies, spread across Eurasia, Africa and the Australo-Papuan region. The dissolution of Rattus is ongoing as many of its constituent species and many genera of Rattini remain unsampled in any molecular study. To address this sampling limitation, we sequenced the mitochondrial cytochrome b (cytb) gene and examined phylogenetic relationships using both Bayesian and Maximum Likelihood algorithms for an expanded set of taxa within Rattus and among closely related genera. Here we place previously unsampled taxa in a phylogenetic context for the first time, including R. burrus, R. hoogerwerfi, R. lugens, and R. mindorensis within the Asian Rattus group, R. facetus within the Australo-Papuan Rattus radiation, and the undescribed 'Bisa Rat' described by Flannery as sister to the recently described genus Halmaheramys. We also present an exploratory foray into the wider topic of Rattus phylogenetics and propose that a reorganisation of the Rattus genus should require that it be a monophyletic group, include at least the type species R. norvegicus and R. rattus (plus their close allies); and exclude the Bandicota/Nesokia clade and other such specialised genera.

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.