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.

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.

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.

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.

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.

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.