A return-on-investment approach for prioritization of rigorous taxonomic research needed to inform responses to the biodiversity crisis.

Global biodiversity loss is a profound consequence of human activity. Disturbingly, biodiversity loss is greater than realized because of the unknown number of undocumented species. Conservation fundamentally relies on taxonomic recognition of species, but only a fraction of biodiversity is described. Here, we provide a new quantitative approach for prioritizing rigorous taxonomic research for conservation. We implement this approach in a highly diverse vertebrate group-Australian lizards and snakes. Of 870 species assessed, we identified 282 (32.4%) with taxonomic uncertainty, of which 17.6% likely comprise undescribed species of conservation concern. We identify 24 species in need of immediate taxonomic attention to facilitate conservation. Using a broadly applicable return-on-investment framework, we demonstrate the importance of prioritizing the fundamental work of identifying species before they are lost.

Comparison of Reptilian Genomes Reveals Deletions Associated with the Natural Loss of γδ T Cells in Squamates.

T lymphocytes or T cells are key components of the vertebrate response to pathogens and cancer. There are two T cell classes based on their TCRs, αß T cells and γδ T cells, and each plays a critical role in immune responses. The squamate reptiles may be unique among the vertebrate lineages by lacking an entire class of T cells, the γδ T cells. In this study, we investigated the basis of the loss of the γδ T cells in squamates. The genome and transcriptome of a sleepy lizard, the skink Tiliqua rugosa, were compared with those of tuatara, Sphenodon punctatus, the last living member of the Rhynchocephalian reptiles. We demonstrate that the lack of TCRγ and TCRδ transcripts in the skink are due to large deletions in the T. rugosa genome. We also show that tuataras are on a growing list of species, including sharks, frogs, birds, alligators, and platypus, that can use an atypical TCRδ that appears to be a chimera of a TCR chain with an Ab-like Ag-binding domain. Tuatara represents the nearest living relative to squamates that retain γδ T cells. The loss of γδTCR in the skink is due to genomic deletions that appear to be conserved in other squamates. The genes encoding the αßTCR chains in the skink do not appear to have increased in complexity to compensate for the loss of γδ T cells.

A giant armoured skink from Australia expands lizard morphospace and the scope of the Pleistocene extinctions.

There are more species of lizards and snakes (squamates) alive today than any other order of land vertebrates, yet their fossil record has been poorly documented compared with other groups. Here, we describe a gigantic Pleistocene skink from Australia based on extensive material that includes much of the skull and postcranial skeleton, and spans ontogenetic stages from neonate to adult. Tiliqua frangens substantially expands the known ecomorphological diversity of squamates. At approximately 2.4 kg, it was more than double the mass of any living skink, with an exceptionally broad, deep skull, squat limbs and heavy, ornamented body armour. It probably filled the armoured herbivore niche that land tortoises (testudinids), absent from Australia, occupy on other continents. Tiliqua frangens and other giant Plio-Pleistocene skinks suggest that small-bodied groups that dominate vertebrate biodiversity might have lost their largest and often most morphologically extreme representatives in the Late Pleistocene, expanding the scope of these extinctions.

Diversity and phylogenetic relationships of haemosporidian and hemogregarine parasites in Australian lizards.

Apicomplexa is a large monophyletic phylum of unicellular, parasitic organisms. Reptiles are hosts to both haemosporidian (Haemosporida) and hemogregarine (Eucoccidiorida) apicomplexan blood parasites. Within reptiles our understanding of their diversity remains limited, with a paucity of information from Australia, despite a high diversity of squamates (snakes and lizards). We provide a preliminary assessment of haemosporidian and hemogregarine diversity occurring in lizards across northern tropical Australia, building on existing data with results from a microscopy and genetic assessment. We screened total of 233 blood slides using microscopy and detected hemogregarines in 25 geckos, 2 skinks and 1 agamid, while haemosporidians were detected in 13 geckos. DNA sequencing of 28 samples of the hemogregarine 18S rRNA (∼900 bp) nuclear gene revealed five lineages of Australian lizard hemogregarines within heteroxenous adeleids. We sequenced 10 samples of Haemosporida mtDNA (cytb & coI: ∼1313 bp) and phylogenetic analysis with 30 previously published sequences revealed that the Australian Haemosporida grouped within the Haemoproteidae but were not supported as a monophyletic clade. Our results demonstrate that there is significant undocumented evolutionary diversity in Australian lizard haemosporidian and hemogregarine parasites, with preliminary evidence of significantly higher infection rates in geckos.

Personality, spatiotemporal ecological variation and resident/explorer movement syndromes in the sleepy lizard.

Individual variation in movement is profoundly important for fitness and offers key insights into the spatial and temporal dynamics of populations and communities. Nonetheless, individual variation in fine-scale movement behaviours is rarely examined even though animal tracking devices offer the long-term, high-resolution, repeatable data in natural conditions that are ideal for studying this variation. Furthermore, of the few studies that consider individual variation in movement, even fewer also consider the internal traits and environmental factors that drive movement behaviour which are necessary for contextualising individual differences in movement patterns. In this study, we GPS tracked a free-ranging population of sleepy lizards Tiliqua rugosa, each Austral spring over 5 years to examine consistent among-individual variation in movement patterns, as well as how these differences were mediated by key internal and ecological factors. We found that individuals consistently differed in a suite of weekly movement traits, and that these traits strongly covaried among-individuals, forming movement syndromes. Lizards fell on a primary movement continuum, from 'residents' that spent extended periods of time residing within smaller core areas of their home range, to 'explorers' that moved greater distances and explored vaster areas of the environment. Importantly, we also found that these consistent differences in lizard movement were related to two ecologically important animal personality traits (boldness and aggression), their sex, key features of the environment (including food availability, and a key water resource), habitat type and seasonal variation (cool/moist vs. hot/drier) in environmental conditions. Broadly, these movement specialisations likely reflect variation in life-history tactics including foraging and mating tactics that ultimately underlie key differences in space use. Such information can be used to connect phenotypic population structure to key ecological and evolutionary processes, for example social networks and disease-transmission pathways, further highlighting the value of examining individual variation in movement behaviour.

Viral metagenomic analysis of fecal samples reveals an enteric virome signature in irritable bowel syndrome.

BACKGROUND: Changes in the enteric microbiota have been suggested to contribute to gastrointestinal diseases, including irritable bowel syndrome. Most of the published work is on bacterial dysbiosis with meager data on the role of the virome in irritable bowel syndrome and other gastrointestinal diseases. In the current study, we therefore aimed to investigate the viral community composition of the gut and test for potential dysbiosis linked to irritable bowel syndrome. RESULTS: A metagenomics analysis on fecal samples of 50 individuals - 30 of whom met the Rome IV criteria for IBS and 20 healthy controls- was conducted. There was a noticeable alteration in viral taxa observed in association with irritable bowel syndrome when compared to healthy individuals - where some eukaryotic viral taxa noticeably prevail over others. We observed a significant decrease in the diversity and abundance of enteric virome particularly in eukaryotic viruses of Megavirales in patients with irritable bowel syndrome. CONCLUSIONS: These findings shed light on a new hypothesis that the alteration of the viral taxa contributes to the pathogenesis of irritable bowel syndrome and related symptoms, and therefore, pave the way for developing a new diagnostic biomarker or anti-viral drugs for the treatment of irritable bowel syndrome.

Temporal variation in thermal plasticity in a free-ranging subalpine lizard.

Thermally variable environments are particularly challenging for ectotherms as physiological functions are thermo-dependent. As a consequence, ectotherms in highly seasonal environments are predicted to have greater thermal plasticity. However, much of our understanding of thermal plasticity comes from controlled experiments in a laboratory setting. Relatively fewer studies investigate thermal plasticity in free-ranging animals living in their natural environment. We investigated the presence of thermal plasticity within a single activity season in adult males of a natural high elevation population of White's skink (Liopholis whitii) in south-eastern Australia. This species lives in a permanent home site (rock crevice and/or burrow), facilitating the repeated capture of the same individuals across the activity season. We monitored the thermal variation across the field site and over the activity season, and tested thermal tolerances and performance of male L. whitii on three occasions across their activity season. Maximum and average temperatures varied across the field site, and temperatures gradually increased across the study period. Evidence of temporal plasticity was identified in the critical thermal minimum and thermal tolerance breadth, but not in the critical thermal maximum. Thermal performance was also found to be plastic, but no temporal patterns were evident. Our temporal plasticity results are consistent which much of the previous literature, but this is one of the first studies to identify these patterns in a free-ranging population. In addition, our results indicate that performance may be more plastic than previous literature suggests. Overall, our study highlights the need to pair laboratory and field studies in order to understand thermal plasticity in an ecologically relevant context.

Plio-Pleistocene diversification and biogeographic barriers in southern Australia reflected in the phylogeography of a widespread and common lizard species.

Palaeoclimatic events and biogeographical processes since the mid-Tertiary have played an important role in shaping the evolution and distribution of Australian fauna. However, their impacts on fauna in southern and arid zone regions of Australia are not well understood. Here we investigate the phylogeography of an Australian scincid lizard, Tiliqua rugosa, across southern Australia using mitochondrial DNA (mtDNA) and 11 nuclear DNA markers (nuDNA), including nine anonymous nuclear loci. Phylogenetic analyses revealed three major mtDNA lineages within T. rugosa, geographically localised north and south of the Murray River in southern Australia, and west of the Nullarbor Plain. Molecular variance and population analyses of both mtDNA and nuDNA haplotypes revealed significant variation among the three populations, although potential introgression of nuDNA markers was also detected for the Northern and Southern population. Coalescent times for major mtDNA lineages coincide with an aridification phase, which commenced after the early Pliocene and increased in intensity during the Late Pliocene-Pleistocene. Species distribution modelling and a phylogeographic diffusion model suggest that the range of T. rugosa may have contracted during the Last Glacial Maximum and the locations of optimal habitat appear to coincide with the geographic origin of several distinct mtDNA lineages. Overall, our analyses suggest that Plio-Pleistocene climatic changes and biogeographic barriers associated with the Nullarbor Plain and Murray River have played a key role in shaping the present-day distribution of genetic diversity in T. rugosa and many additional ground-dwelling animals distributed across southern Australia.

Low food availability during gestation enhances offspring post-natal growth, but reduces survival, in a viviparous lizard.

The environment experienced by a mother can have profound effects on the fitness of her offspring (i.e., maternal effects). Maternal effects can be adaptive when the developmental environments experienced by offspring promote phenotypes that provide fitness benefits either via matching offspring phenotype to the post-developmental environment (also known as anticipatory maternal effects) or through direct effects on offspring growth and survival. We tested these hypotheses in a viviparous lizard using a factorial experimental design in which mothers received either high or low amounts of food during gestation, and resultant offspring were raised on either high or low amounts of food post-birth. We found no effect of food availability during gestation on reproductive traits of mothers or offspring traits at birth. However, offspring from mothers who received low food during gestation exhibited a greater increase in condition in the post-birth period, suggesting some form of priming of offspring by mothers to cope with an anticipated poor environment after birth. Offspring that received low food during gestation were also more likely to die, suggesting a trade-off for this accelerated growth. There were also significant effects of post-birth food availability on offspring snout-vent length and body condition growth, with offspring with high food availability post birth doing better. However, the effects of the pre- and post-natal resource evnironment on offspring growth were independent on one another, therefore, providing no support for the presence of anticipatory maternal effects in the traditional sense.

Selection outweighs drift at a fine scale: Lack of MHC differentiation within a family living lizard across geographically close but disconnected rocky outcrops.

The highly polymorphic genes of the major histocompatibility complex (MHC) are involved in disease resistance, mate choice and kin recognition. Therefore, they are widely used markers for investigating adaptive variation. Although selection is the key driver, gene flow and genetic drift also influence adaptive genetic variation, sometimes in opposing ways and with consequences for adaptive potential. To further understand the processes that generate MHC variation, it is helpful to compare variation at the MHC with that at neutral genetic loci. Differences in MHC and neutral genetic variation are useful for inferring the relative influence of selection, gene flow and drift on MHC variation. To date, such investigations have usually been undertaken at a broad spatial scale. Yet, evolutionary and ecological processes can occur at a fine spatial scale, particularly in small or fragmented populations. We investigated spatial patterns of MHC variation among three geographically close, naturally discrete, sampling sites of Egernia stokesii, an Australian lizard. The MHC of E. stokesii has recently been characterized, and there is evidence for historical selection on the MHC. We found E. stokesii MHC weakly differentiated among sites compared to microsatellites, suggesting selection, acting similarly at each site, has outweighed any effects of low gene flow or of genetic drift on E. stokesii MHC variation. Our findings demonstrate the strength of selection in shaping patterns of MHC variation or consistency at a fine spatial scale.

Avian predation intensity as a driver of clinal variation in colour morph frequency.

Phenotypic variation provides the framework for natural selection to work upon, enabling adaptive evolution. One of the most discernible manifestations of phenotypic variability is colour variation. When this variation is discrete, genetically based colour pattern morphs occur simultaneously within a population. Why and how colour polymorphisms are maintained is an evolutionary puzzle. Several evolutionary drivers have been hypothesized as influencing clinal patterns of morph frequency, with spatial variation in climate and predation being considered especially important. Despite this, no study has examined both of their roles simultaneously. The aims of this study were to: (a) examine the covariation of physiology, environmental variables and colouration at a local scale; and (b) determine if these factors and their interplay explain broad clinal variation in morph frequency. We used the lizard Liopholis whitii as a model system, as this species displays a discrete, heritable polymorphism for colour pattern (plain-backed, patterned morphs) whose morph frequency varies latitudinally. We measured reflectance, field activity temperatures and microhabitat structure to test for differences in crypsis, thermal biology and microhabitat selection of patterned and plain-backed morphs within a single population where colour morphs occur sympatrically. We then used data from the literature to perform a broad-scale analysis to identify whether these factors also explained the latitudinal variation of morph frequency in this species. At the local scale, plain-backed morphs were found to be less cryptic than patterned morphs while no other differences were detected in terms of thermal biology, dorsal reflectance and microhabitat use. At a broader scale, predation was the most influential factor mediating morph frequency across latitudes. However, the observed pattern of morph frequency is opposite to what the modelling results suggest in that the incidence of the least cryptic morph is highest where predation pressure is most severe. Clinal variation in the level of background matching between morphs or the potential reproductive advantage by the plain-backed morph may, instead, be driving the observed morph frequency. Together, these results provide key insights into the evolution of local adaptation as well as the ecological forces involved in driving the dynamics of colour polymorphism.

Egernia stokesii (gidgee skink) MHC I positively selected sites lack concordance with HLA peptide binding regions.

Genes of the major histocompatibility complex (MHC) play an important role in vertebrate disease resistance, kin recognition and mate choice. Mammalian MHC is the most widely characterised of all vertebrates, and attention is often given to the peptide binding regions of the MHC because they are presumed to be under stronger selection than non-peptide binding regions. For vertebrates where the MHC is less well understood, researchers commonly use the amino acid positions of the peptide binding regions of the human leukocyte antigen (HLA) to infer the peptide binding regions within the MHC sequences of their taxon of interest. However, positively selected sites within MHC have been reported to lack correspondence with the HLA in fish, frogs, birds and reptiles including squamates. Despite squamate diversity, the MHC has been characterised in few snakes and lizards. The Egernia group of scincid lizards is appropriate for investigating mechanisms generating MHC variation, as their inclusion will add a new lineage (i.e. Scincidae) to studies of selection on the MHC. We aimed to identify positively selected sites within the MHC of Egernia stokesii and then determine if these sites corresponded with the peptide binding regions of the HLA. Six positively selected sites were identified within E. stokesii MHC I, only two were homologous with the HLA. E. stokesii positively selected sites corresponded more closely to non-lizard than other lizard taxa. The characterisation of the MHC of more intermediate taxa within the squamate order is necessary to understand the evolution of the MHC across all vertebrates.

Genes and Group Membership Predict Gidgee Skink (Egernia stokesii) Reproductive Pairs.

Due to their role in mate choice, disease resistance and kin recognition, genes of the major histocompatibility complex (MHC) are good candidates for investigating genetic-based mate choice. MHC-based mate choice is context dependent and influenced by many factors including social structure. Social structure diversity makes the Egernia group of lizards suitable for comparative studies of MHC-based mate choice. We investigated mate choice in the gidgee skink (Egernia stokesii), a lizard that exhibits high levels of social group and spatial stability. Group membership was incorporated into tests of the good genes as heterozygosity and compatible genes hypotheses for adaptive (MHC) and neutral (microsatellite) genetic diversity (n = 47 individuals genotyped). Females were more likely to pair with a male with higher MHC diversity and with whom they had a lower degree of microsatellite relatedness. Males were more likely to pair with a female with higher microsatellite heterozygosity and with whom they shared a lower proportion of MHC alleles. Lizards were more likely to mate with an individual from within, rather than outside, their social group, which confirmed earlier findings for this species and indicated mate choice had already largely occurred prior to either social group formation or acceptance of an individual into an existing group. Thus, a combination of genes and group membership, rather than group membership alone, predicted mate choice in this species. This work will contribute to an enhanced understanding of squamate group formation and a deeper understanding of the evolution of sociality within all vertebrates.

Dispersal responses override density effects on genetic diversity during post-disturbance succession.

Dispersal fundamentally influences spatial population dynamics but little is known about dispersal variation in landscapes where spatial heterogeneity is generated predominantly by disturbance and succession. We tested the hypothesis that habitat succession following fire inhibits dispersal, leading to declines over time in genetic diversity in the early successional gecko Nephrurus stellatus We combined a landscape genetics field study with a spatially explicit simulation experiment to determine whether successional patterns in genetic diversity were driven by habitat-mediated dispersal or demographic effects (declines in population density leading to genetic drift). Initial increases in genetic structure following fire were likely driven by direct mortality and rapid population expansion. Subsequent habitat succession increased resistance to gene flow and decreased dispersal and genetic diversity inN. stellatus Simulated changes in population density alone did not reproduce these results. Habitat-mediated reductions in dispersal, combined with changes in population density, were essential to drive the field-observed patterns. Our study provides a framework for combining demographic, movement and genetic data with simulations to discover the relative influence of demography and dispersal on patterns of landscape genetic structure. Our results suggest that succession can inhibit connectivity among individuals, opening new avenues for understanding how disturbance regimes influence spatial population dynamics.

Widespread horizontal transfer of retrotransposons.

In higher organisms such as vertebrates, it is generally believed that lateral transfer of genetic information does not readily occur, with the exception of retroviral infection. However, horizontal transfer (HT) of protein coding repetitive elements is the simplest way to explain the patchy distribution of BovB, a long interspersed element (LINE) about 3.2 kb long, that has been found in ruminants, marsupials, squamates, monotremes, and African mammals. BovB sequences are a major component of some of these genomes. Here we show that HT of BovB is significantly more widespread than believed, and we demonstrate the existence of two plausible arthropod vectors, specifically reptile ticks. A phylogenetic tree built from BovB sequences from species in all of these groups does not conform to expected evolutionary relationships of the species, and our analysis indicates that at least nine HT events are required to explain the observed topology. Our results provide compelling evidence for HT of genetic material that has transformed vertebrate genomes.

Life history influences how fire affects genetic diversity in two lizard species.

'Fire mosaics' are often maintained in landscapes to promote successional diversity in vegetation with little understanding of how this will affect ecological processes in animal populations such as dispersal, social organization and re-establishment. To investigate these processes, we conducted a replicated, spatiotemporal landscape genetics study of two Australian woodland lizard species [Amphibolurus norrisi (Agamidae) and Ctenotus atlas (Scincidae)]. Agamids have a more complex social and territory structure than skinks, so fire might have a greater impact on their population structure and thus genetic diversity. Genetic diversity increased with time since fire in C. atlas and decreased with time since fire in A. norrisi. For C. atlas, this might reflect its increasing population size after fire, but we could not detect increased gene flow that would reduce the loss of genetic diversity through genetic drift. Using landscape resistance analyses, we found no evidence that postfire habitat succession or topography affected gene flow in either species and we were unable to distinguish between survival and immigration as modes of postfire re-establishment. In A. norrisi, we detected female-biased dispersal, likely reflecting its territorial social structure and polygynous mating system. The increased genetic diversity in A. norrisi in recently burnt habitat might reflect a temporary disruption of its territoriality and increased male dispersal, a hypothesis that was supported with a simulation experiment. Our results suggest that the effects of disturbance on genetic diversity will be stronger for species with territorial social organization.

Characterisation of novel and rare Y-chromosome short tandem repeat alleles in self-declared South Australian Aboriginal database.

Y-chromosome short tandem repeats (Y-STRs) are used in forensic science laboratories all over the world, as their application is wide and often vital in solving casework. Analysis of an in-house database of South Australian self-declared Aboriginal males held by Forensic Science South Australia (FSSA) using the Applied Biosystem's AmpFℓSTR® Yfiler™ PCR Amplification Kit revealed 43 variant Y-STR alleles at 6 of the 17 loci. All variant alleles were sequenced to determine the exact repeat structure for each. As a high level of admixture has previously been found within the SA Aboriginal database, samples were haplogrouped using Y-SNPs to determine their likely geographical origin. Although a number of variant alleles were associated with non-Aboriginal Y-haplogroups, a high frequency was observed within the Australian K-M9 lineage. Detailed knowledge of these variant alleles may have further application in the development of new DNA markers for identification purposes, and in population and evolutionary studies of Australian Aborigines.

Visualisation and detection of latent DNA deposited by pangolin scales onto plastic packaging materials.

We report on the detection and visualisation of latent DNA from pangolin scales deposited onto a plastic packaging material through the use of a nucleic acid staining dye. This latent DNA deposited by pangolin scales was subsequently isolated and analysed using DNA barcoding method. Pangolins are the most illegally traded mammalian species due to the demand for their scales and meat. The demand for their scales were mostly fuelled by its use in traditional medicines. The scales are usually packed into bags and transported globally via sea routes. This is the first report detailing the detection of trace latent DNA from processed wildlife products, on surfaces of bags that they were packaged in. Prior to this report, it was not known if the dried pangolin scales contained transferable quantities of biological material for DNA analyses. To address this, scales were removed from a roadkill Sunda pangolin (Manis javanica), processed by drying and packaged into one of five plastic bags. The presence of pangolin latent DNA was detected on the surface of the plastic bags and visualised using Diamond™ nucleic acid dye. Swabs were then used to recover the stained biological material from various locations in the five bags. The DNA was isolated and quantified using a newly designed quantitative PCR (qPCR) specific to M. javanica to amplify a fragment of the mitochondrial DNA cytochrome b gene. There was a positive correlation between the number of stained particles and DNA quantity, and a greater number of stained particles were found at the bottom of the bag than were found at the top. Conventional PCR targeting part of the cyt b gene amplified a product from all 30 samples taken from the bags and in all cases, sequence data generated matched that of the Sunda pangolin, as expected. All negative controls yielded no results. The method described here is the very first use of a nucleic acid staining dye to detect latent DNA from a mammalian species, other than humans, and highlights the opportunity for further use of Diamond™ nucleic acid dye in wildlife forensic science. It is anticipated that this method will be invaluable in retrieving latent DNA deposited by wildlife products from the environment in which they were contained, to determine the presence of these illegal wildlife products even when previously hidden, inaccessible, or no longer present physically. Further research is required to understand if the use on non-human mammalian wildlife species is feasible.

Anonymous nuclear loci in non-model organisms: making the most of high-throughput genome surveys.

MOTIVATION: When working with non-model organisms, few if any species-specific markers are available for phylogenetic, phylogeographic and population studies. Therefore, researchers often try to adapt markers developed in distantly related taxa, resulting in poor amplification and ascertainment bias in their target taxa. Markers can be developed de novo and anonymous nuclear loci (ANL) are proving to be a boon for researchers seeking large numbers of fast-evolving, independent loci. However, the development of ANL can be laboratory intensive and expensive. A workflow is described to identify suitable low-copy anonymous loci from high-throughput shotgun sequences, dramatically reducing the cost and time required to develop these markers and produce robust multilocus datasets. RESULTS: By successively removing repetitive and evolutionary conserved sequences from low coverage shotgun libraries, we were able to isolate thousands of potential ANL. Empirical testing of loci developed from two reptile taxa confirmed that our methodology yields markers with comparable amplification rates and nucleotide diversities to ANLs developed using other methodologies. Our approach capitalizes on next-generation sequencing technologies to enable the development of phylogenetic, phylogeographic and population markers for taxa lacking suitable genomic resources.