A new species of barred frog, Mixophyes (Anura: Myobatrachidae) from south-eastern Australia identified by molecular genetic analyses

Mixophyes are large ground-dwelling myobatrachid frogs from eastern Australia and New Guinea. Several of the species found in mid-eastern and south-eastern Australia are listed as threatened, due largely to declines presumably caused by the amphibian disease chytridiomycosis. Given the wide distribution of several of these species and that their distributions cross well-known biogeographic boundaries that often correspond to deep genetic breaks or species boundaries among closely related vertebrates, we undertook a molecular genetic assessment of population structure across the range of each species to determine the presence of undescribed species. Of the four species of Mixophyes subject to molecular population genetic analyses, one, the Stuttering Frog (Mixophyes balbus), showed a level of diversity consistent with the presence of two species. Morphometric, meristic and bioacoustic analyses corroborate these distinctions, and a new species is described for the populations south of the Macleay River valley in mid-eastern New South Wales to east Gippsland in Victoria. Applying the IUCN Red List threat criteria the new species meets the conservation status assessment criteria for Endangered 2B1a,b because its extent of occupancy and area of occupancy are below the threshold value and it has declined and disappeared from the southern two thirds of its distribution over the past 30 years.

Parallel decay of vision genes in subterranean water beetles

In the framework of neutral theory of molecular evolution, genes specific to the development and function of eyes in subterranean animals living in permanent darkness are expected to evolve by relaxed selection, ultimately becoming pseudogenes. However, definitive empirical evidence for the role of neutral processes in the loss of vision over evolutionary time remains controversial. In previous studies, we characterized an assemblage of independently-evolved water beetle (Dytiscidae) species from a subterranean archipelago in Western Australia, where parallel vision and eye loss have occurred. Using a combination of transcriptomics and exon capture, we present evidence of parallel coding sequence decay, resulting from the accumulation of frameshift mutations and premature stop codons, in eight phototransduction genes (arrestins, opsins, ninaC and transient receptor potential channel genes) in 32 subterranean species in contrast to surface species, where these genes have open reading frames. Our results provide strong evidence to support neutral evolutionary processes as a major contributing factor to the loss of phototransduction genes in subterranean animals, with the ultimate fate being the irreversible loss of a light detection system.

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.

Genetic Insights Into the Introduction History of Black Rats Into the Eastern Indian Ocean

Islands can be powerful demonstrations of how destructive invasive species can be on endemic faunas and insular ecologies. Oceanic islands in the eastern Indian Ocean have suffered dramatically from the impact of one of the world’s most destructive invasive species, the black rat, causing the loss of endemic terrestrial mammals and ongoing threats to ground-nesting birds. We use molecular genetic methods on both ancient and modern samples to establish the origins and minimum invasion frequencies of black rats on Christmas Island and the Cocos-Keeling Islands. We find that each island group had multiple incursions of black rats from diverse geographic and phylogenetic sources. Furthermore, contemporary black rat populations on these islands are highly admixed to the point of potentially obscuring their geographic sources. These hybridisation events between black rat taxa also pose potential dangers to human populations on the islands from novel disease risks. Threats of ongoing introductions from yet additional geographic sources is highlighted by genetic identifications of black rats found on ships, which provides insight into how recent ship-borne human smuggling activity to Christmas Island can negatively impact its endemic species.

Development and evaluation of a custom bait design based on 469 single-copy protein-coding genes for exon capture of isopods (Philosciidae: Haloniscus)

Transcriptome-based exon capture approaches, along with next-generation sequencing, are allowing for the rapid and cost-effective production of extensive and informative phylogenomic datasets from non-model organisms for phylogenetics and population genetics research. These approaches generally employ a reference genome to infer the intron-exon structure of targeted loci and preferentially select longer exons. However, in the absence of an existing and well-annotated genome, we applied this exon capture method directly, without initially identifying intron-exon boundaries for bait design, to a group of highly diverse Haloniscus (Philosciidae), paraplatyarthrid and armadillid isopods, and examined the performance of our methods and bait design for phylogenetic inference. Here, we identified an isopod-specific set of single-copy protein-coding loci, and a custom bait design to capture targeted regions from 469 genes, and analysed the resulting sequence data with a mapping approach and newly-created post-processing scripts. We effectively recovered a large and informative dataset comprising both short (<100 bp) and longer (>300 bp) exons, with high uniformity in sequencing depth. We were also able to successfully capture exon data from up to 16-year-old museum specimens along with more distantly related outgroup taxa, and efficiently pool multiple samples prior to capture. Our well-resolved phylogenies highlight the overall utility of this methodological approach and custom bait design, which offer enormous potential for application to future isopod, as well as broader crustacean, molecular studies.

The critical thermal maximum of diving beetles (Coleoptera: Dytiscidae): a comparison of subterranean and surface-dwelling species

Thermal tolerance limits in animals are often thought to be related to temperature and thermal variation in their environment. Recently, there has been a focus on studying upper thermal limits due to the likelihood for climate change to expose more animals to higher temperatures and potentially extinction. Organisms living in underground environments experience reduced temperatures and thermal variation in comparison to species living in surface habitats, but how these impact their thermal tolerance limits are unclear. In this study, we compare the thermal critical maximum (CTmax) of two subterranean diving beetles (Dytiscidae) to that of three related surface-dwelling species. Our results show that subterranean species have a lower CTmax (38.3-39.0°C) than surface species (42.0-44.5°C). The CTmax of subterranean species is ∼10°C higher than the highest temperature recorded within the aquifer. Groundwater temperature varied between 18.4°C and 28.8°C, and changes with time, depth and distance across the aquifer. Seasonal temperature fluctuations were 0.5°C at a single point, with the maximum heating rate being ∼1000x lower (0.008°C/hour) than that recorded in surface habitats (7.98°C/hour). For surface species, CTmax was 7-10°C higher than the maximum temperature in their habitats, with daily fluctuations from ∼1°C to 16°C and extremes of 6.9°C and 34.9°C. These findings suggest that subterranean dytiscid beetles are unlikely to reach their CTmax with a predicted warming of 1.3-5.1°C in the region by 2090. However, the impacts of long-term elevated temperatures on fitness, different life stages and other species in the beetle's trophic food web are unknown.

Taxonomic revision of south-eastern Australian giant burrowing frogs (Anura: Limnodynastidae: Heleioporus Gray)

The rarely encountered giant burrowing frog, Heleioporus australiacus, is distributed widely in a variety of sclerophyll forest habitats east of the Great Dividing Range in south-eastern Australia. Analyses of variation in nucleotide sequences of the mitochondrial ND4 gene and thousands of nuclear gene SNPs revealed the presence of two deeply divergent lineages. Multivariate morphological comparisons show the two lineages differ in body proportions with > 91% of individuals being correctly classified in DFA. The two lineages differ in the number and size of spots on the lateral surfaces and the degree by which the cloaca is surrounded by colour patches. The mating calls are significantly different in number of pulses in the note. The presence of a F2 hybrid in the area where the distribution of the two taxa come into closest proximity leads us to assign subspecies status to the lineages, as we have not been able to assess the extent of potential genetic introgression. In our sampling, the F2 hybrid sample sits within an otherwise unsampled gap of ~90km between the distributions of the two lineages. The nominate northern sub-species is restricted to the Sydney Basin bioregion, while the newly recognised southern subspecies occurs from south of the Kangaroo Valley in the mid-southern coast of New South Wales to near Walhalla in central Gippsland in Victoria. The habitat of the two subspecies is remarkably similar. Adults spend large portions of their lives on the forest floor where they forage and burrow in a variety of vegetation communities. The southern subspecies occurs most commonly in dry sclerophyll forests with an open understory in the south and in open forest and heath communities with a dense understory in the north of its distribution. The northern subspecies is also found in dry open forests and heaths in association with eroded sandstone landscapes in the Sydney Basin bioregion. Males of both taxa call from both constructed burrows and open positions on small streams, differing from the five Western Australian species of Heleioporus where males call only from constructed burrows. Using the IUCN Red List process, we found that the extent of occupancy and area of occupancy along with evidence of decline for both subspecies are consistent with the criteria for Endangered (A2(c)B2(a)(b)).

Unravelling the taxonomy and identification of a problematic group of benthic fishes from tropical rivers (Gobiidae: Glossogobius)

Flathead gobies (genus Glossogobius) include c. 40 small- to medium-sized benthic fishes found primarily in freshwater habitats across the Indo-Pacific, having biodiversity value as well as cultural and economic value as food fishes, especially in developing countries. To help resolve considerable confusion regarding the identification of some of the larger-growing Glossogobius species, a systematic framework was established using nuclear genetic markers, mitochondrial DNA barcoding and phenotypic evidence for a geographically widespread collection of individuals from the waterways of tropical northern Australia. Species boundaries and distribution patterns were discordant with those previously reported, most notably for the tank goby Glossogobius giuris, which included a cryptic species. Genetic divergence was matched with accompanying unique visual characters that aid field identification. Additional taxonomic complexity was also evident, by comparison with DNA barcodes from international locations, suggesting that the specific names applicable for two of the candidate species in Australia remain unresolved due to confusion surrounding type specimens. Although flathead gobies are assumed to be widespread and common, this study demonstrates that unrealised taxonomic and ecological complexity is evident, and this will influence assessments of tropical biodiversity and species conservation. This study supports the need for taxonomic studies of freshwater fishes to underpin management in areas subject to significant environmental change.

Surprising Pseudogobius: Molecular systematics of benthic gobies reveals new insights into estuarine biodiversity (Teleostei: Gobiiformes)

Snubnose gobies (genus Pseudogobius: Gobionellinae) are ubiquitous to, and important components of, estuarine ecosystems of the Indo-west Pacific. These small benthic fishes occur in freshwater, brackish and marine habitats such as mangroves, sheltered tide pools and lowland streams, and represent a model group for understanding the biodiversity and biogeography of estuarine fauna. To develop the species-level framework required for a concurrent morphological taxonomic appraisal, we undertook thorough sampling around the extensive Australian coastline, referenced to international locations, as part of a molecular systematic review using both nuclear and mitochondrial markers. The results indicate that while there are currently eight recognised species, the true diversity is close to double this, with a hotspot of endemism located in Australia. Complicated patterns were observed in southern Australia owing to two differing zones of introgression/admixture. Key drivers of diversity in the group appear to include plate tectonics, latitude, and historic barriers under glacial maxima, where an interplay between ready dispersal and habitat specialisation has led to regional panmixia but frequent geographic compartmentalisation within past and present landscapes. The findings have significant implications for biodiversity conservation, coastal and estuarine development, the basic foundations of field ecology, and for applied use such as in biomonitoring.

The tuatara genome reveals ancient features of amniote evolution

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

A draft genome assembly of the eastern banjo frog Limnodynastes dumerilii dumerilii (Anura: Limnodynastidae)

Amphibian genomes are usually challenging to assemble due to their large genome size and high repeat content. The Limnodynastidae is a family of frogs native to Australia, Tasmania and New Guinea. As an anuran lineage that successfully diversified on the Australian continent, it represents an important lineage in the amphibian tree of life but lacks reference genomes. Here we sequenced and annotated the genome of the eastern banjo frog Limnodynastes dumerilii dumerilii to fill this gap. The total length of the genome assembly is 2.38 Gb with a scaffold N50 of 285.9 kb. We identified 1.21 Gb of non-redundant sequences as repetitive elements and annotated 24,548 protein-coding genes in the assembly. BUSCO assessment indicated that more than 94% of the expected vertebrate genes were present in the genome assembly and the gene set. We anticipate that this annotated genome assembly will advance the future study of anuran phylogeny and amphibian genome evolution.

Horizontal transfer of BovB and L1 retrotransposons in eukaryotes

BACKGROUND

Transposable elements (TEs) are mobile DNA sequences, colloquially known as jumping genes because of their ability to replicate to new genomic locations. TEs can jump between organisms or species when given a vector of transfer, such as a tick or virus, in a process known as horizontal transfer. Here, we propose that LINE-1 (L1) and Bovine-B (BovB), the two most abundant TE families in mammals, were initially introduced as foreign DNA via ancient horizontal transfer events.

RESULTS

Using analyses of 759 plant, fungal and animal genomes, we identify multiple possible L1 horizontal transfer events in eukaryotic species, primarily involving Tx-like L1s in marine eukaryotes. We also extend the BovB paradigm by increasing the number of estimated transfer events compared to previous studies, finding new parasite vectors of transfer such as bed bug, leech and locust, and BovB occurrences in new lineages such as bat and frog. Given that these transposable elements have colonised more than half of the genome sequence in today's mammals, our results support a role for horizontal transfer in causing long-term genomic change in new host organisms.

CONCLUSIONS

We describe extensive horizontal transfer of BovB retrotransposons and provide the first evidence that L1 elements can also undergo horizontal transfer. With the advancement of genome sequencing technologies and bioinformatics tools, we anticipate our study to be a valuable resource for inferring horizontal transfer from large-scale genomic data.

Characterisation of major histocompatibility complex class I transcripts in an Australian dragon lizard

Characterisation of squamate major histocompatibility complex (MHC) genes has lagged behind other taxonomic groups. MHC genes encode cell-surface glycoproteins that present self- and pathogen-derived peptides to T cells and play a critical role in pathogen recognition. Here we characterise MHC class I transcripts for an agamid lizard (Ctenophorus decresii) and investigate the evolution of MHC class I in Iguanian lizards. An iterative assembly strategy was used to identify six full-length C. decresii MHC class I transcripts, which were validated as likely to encode classical class I MHC molecules. Evidence for exon shuffling recombination was uncovered for C. decresii transcripts and Bayesian phylogenetic analysis of Iguanian MHC class I sequences revealed a pattern expected under a birth-and-death mode of evolution. This work provides a stepping stone towards further research on the agamid MHC class I region.

DNA barcoding of microgastrine parasitoid wasps (Hymenoptera: Braconidae) using high-throughput methods more than doubles the number of species known for Australia

The Microgastrinae are a hugely diverse subfamily of endoparasitoid wasps of lepidopteran caterpillars. They are important in agriculture as biological control agents and play a significant ecological role in the regulation of caterpillar populations. Whilst the group has been the focus of intensive rearing and DNA barcoding studies in the Northern Hemisphere, the Australian fauna has received little attention. In total, 99 species have been described from or have been introduced into Australia, but the real species diversity for the region is clearly much larger than this. In this study, museum ethanol samples and recent field collections were mined for hundreds of specimens of microgastrine wasps, which were then barcoded for the COI region, ITS2 ribosomal spacer and the wingless nuclear genes, using a pooled sequencing approach on an Illumina Miseq system. Full COI sequences were obtained for 525 individuals which, when combined with 162 publicly available sequences, represented 417 haplotypes, and a total of 236 species were delimited using a consensus approach. By more than doubling the number of known microgastrine wasp species in Australia, our study highlights the value of DNA barcoding in the context of employing high-throughput sequencing methods of bulk ethanol museum collections for biodiversity assessment.

Superior ab initio identification, annotation and characterisation of TEs and segmental duplications from genome assemblies

Transposable Elements (TEs) are mobile DNA sequences that make up significant fractions of amniote genomes. However, they are difficult to detect and annotate ab initio because of their variable features, lengths and clade-specific variants. We have addressed this problem by refining and developing a Comprehensive ab initio Repeat Pipeline (CARP) to identify and cluster TEs and other repetitive sequences in genome assemblies. The pipeline begins with a pairwise alignment using krishna, a custom aligner. Single linkage clustering is then carried out to produce families of repetitive elements. Consensus sequences are then filtered for protein coding genes and then annotated using Repbase and a custom library of retrovirus and reverse transcriptase sequences. This process yields three types of family: fully annotated, partially annotated and unannotated. Fully annotated families reflect recently diverged/young known TEs present in Repbase. The remaining two types of families contain a mixture of novel TEs and segmental duplications. These can be resolved by aligning these consensus sequences back to the genome to assess copy number vs. length distribution. Our pipeline has three significant advantages compared to other methods for ab initio repeat identification: 1) we generate not only consensus sequences, but keep the genomic intervals for the original aligned sequences, allowing straightforward analysis of evolutionary dynamics, 2) consensus sequences represent low-divergence, recently/currently active TE families, 3) segmental duplications are annotated as a useful by-product. We have compared our ab initio repeat annotations for 7 genome assemblies to other methods and demonstrate that CARP compares favourably with RepeatModeler, the most widely used repeat annotation package.

LINEs between Species: Evolutionary Dynamics of LINE-1 Retrotransposons across the Eukaryotic Tree of Life

LINE-1 (L1) retrotransposons are dynamic elements. They have the potential to cause great genomic change because of their ability to ‘jump’ around the genome and amplify themselves, resulting in the duplication and rearrangement of regulatory DNA. Active L1, in particular, are often thought of as tightly constrained, homologous and ubiquitous elements with well-characterized domain organization. For the past 30 years, model organisms have been used to define L1s as 6–8 kb sequences containing a 5′-UTR, two open reading frames working harmoniously in cis, and a 3′-UTR with a polyA tail. In this study, we demonstrate the remarkable and overlooked diversity of L1s via a comprehensive phylogenetic analysis of elements from over 500 species from widely divergent branches of the tree of life. The rapid and recent growth of L1 elements in mammalian species is juxtaposed against the diverse lineages found in other metazoans and plants. In fact, some of these previously unexplored mammalian species (e.g. snub-nosed monkey, minke whale) exhibit L1 retrotranspositional ‘hyperactivity’ far surpassing that of human or mouse. In contrast, non-mammalian L1s have become so varied that the current classification system seems to inadequately capture their structural characteristics. Our findings illustrate how both long-term inherited evolutionary patterns and random bursts of activity in individual species can significantly alter genomes, highlighting the importance of L1 dynamics in eukaryotes.

Isolation and characterisation of 12 polymorphic microsatellite loci for the threatened mound-building malleefowl, Leipoa ocellata (Aves : Megapodiidae)

Using 454 pyrosequencing and genomic enrichment techniques we developed 12 polymorphic markers for the endangered megapode, the malleefowl (Leipoa ocellata). Pyrosequencing on a 454 instrument resulted in 65 536 reads, with 3469 containing microsatellite repeats. Of these, 232 contained unique flanking sequences and had more than 8 repeat motifs. We chose 13 loci based on reliability of amplification and, from these, 12 unlinked loci were selected for genotyping. In a single population (n = 19), the 12 markers were moderately polymorphic (number of alleles per locus range = 3–7) and showed moderate to high levels of heterozygosity (0.285–0.882). Nine microsatellite primer pairs developed from the brush turkey (Alectura lathami), the closest living relative of the malleefowl in the family, Megapodiidae, failed to reliably amplify malleefowl DNA.

A new species of Near-shore Marine Goby (Pisces: Gobiidae: Nesogobius) from Kangaroo Island, Australia

Nesogobius is one of two goby genera with all species wholly restricted to temperate Australian waters. Described here is a new member of the genus discovered during near-shore marine and estuarine fish sampling along the central southern Australian coastline. The tiger sandgoby Nesogobius tigrinus sp. nov. is distinguished from other congeners by a combination of colouration including four prominent vertical black bars on males; morphological characters involving body scales (large), head scales (naked), body depth (slender) and gill opening (wide); meristic counts including a lack of second dorsal and anal fin spines; and mitochondrial DNA sequence divergence. The species appears to be a narrow range endemic, restricted to specific sub-tidal habitat in the unique sheltered embayments of northeast Kangaroo Island. This study forms part of ongoing investigations to more fully describe the biodiversity and conservation requirements of the regional ichthyofauna.

Opsin transcripts of predatory diving beetles: a comparison of surface and subterranean photic niches

The regressive evolution of eyes has long intrigued biologists yet the genetic underpinnings remain opaque. A system of discrete aquifers in arid Australia provides a powerful comparative means to explore trait regression at the genomic level. Multiple surface ancestors from two tribes of diving beetles (Dytiscidae) repeatedly invaded these calcrete aquifers and convergently evolved eye-less phenotypes. We use this system to assess transcription of opsin photoreceptor genes among the transcriptomes of two surface and three subterranean dytiscid species and test whether these genes have evolved under neutral predictions. Transcripts for UV, long-wavelength and ciliary-type opsins were identified from the surface beetle transcriptomes. Two subterranean beetles showed parallel loss of all opsin transcription, as expected under 'neutral' regressive evolution. The third species Limbodessus palmulaoides retained transcription of a long-wavelength opsin (lwop) orthologue, albeit in an aphotic environment. Tests of selection on lwop indicated no significant differences between transcripts derived from surface and subterranean habitats, with strong evidence for purifying selection acting on L. palmulaoides lwop. Retention of sequence integrity and the lack of evidence for neutral evolution raise the question of whether we have identified a novel pleiotropic role for lwop, or an incipient phase of pseudogene development.

Opsin transcripts of predatory diving beetles: a comparison of surface and subterranean photic niches

The regressive evolution of eyes has long intrigued biologists yet the genetic underpinnings remain opaque. A system of discrete aquifers in arid Australia provides a powerful comparative means to explore trait regression at the genomic level. Multiple surface ancestors from two tribes of diving beetles (Dytiscidae) repeatedly invaded these calcrete aquifers and convergently evolved eye-less phenotypes. We use this system to assess transcription of opsin photoreceptor genes among the transcriptomes of two surface and three subterranean dytiscid species and test whether these genes have evolved under neutral predictions. Transcripts for UV, long-wavelength and ciliary-type opsins were identified from the surface beetle transcriptomes. Two subterranean beetles showed parallel loss of all opsin transcription, as expected under 'neutral' regressive evolution. The third species Limbodessus palmulaoides retained transcription of a long-wavelength opsin (lwop) orthologue, albeit in an aphotic environment. Tests of selection on lwop indicated no significant differences between transcripts derived from surface and subterranean habitats, with strong evidence for purifying selection acting on L. palmulaoides lwop. Retention of sequence integrity and the lack of evidence for neutral evolution raise the question of whether we have identified a novel pleiotropic role for lwop, or an incipient phase of pseudogene development.

Evaluating evolutionary history in the face of high gene tree discordance in Australian Gehyra (Reptilia: Gekkonidae)

Species tree methods have provided improvements for estimating species relationships and the timing of diversification in recent radiations by allowing for gene tree discordance. Although gene tree discordance is often observed, most discordance is attributed to incomplete lineage sorting rather than other biological phenomena, and the causes of discordance are rarely investigated. We use species trees from multi-locus data to estimate the species relationships, evolutionary history and timing of diversification among Australian Gehyra-a group renowned for taxonomic uncertainty and showing a large degree of gene tree discordance. We find support for a recent Asian origin and two major clades: a tropically adapted clade and an arid adapted clade, with some exceptions, but no support for allopatric speciation driven by chromosomal rearrangement in the group. Bayesian concordance analysis revealed high gene tree discordance and comparisons of Robinson-Foulds distances showed that discordance between gene trees was significantly higher than that generated by topological uncertainty within each gene. Analysis of gene tree discordance and incomplete taxon sampling revealed that gene tree discordance was high whether terminal taxon or gene sampling was maximized, indicating discordance is due to biological processes, which may be important in contributing to gene tree discordance in many recently diversified organisms.

Random non-coding fragments of lizard DNA: anonymous nuclear loci for the Australian skink, Tiliqua rugosa, and their utility in other Egernia-group species

We report the development of 48 anonymous nuclear loci from the Australian skink Tiliqua rugosa using 454 sequencing. These loci amplified across a Western Australian lineage (47 loci), a ‘northern’ lineage (48 loci) and a ‘southern’ lineage (46 loci). We further tested amplification for the related T. adelaidensis and Egernia stokesii where 37 and 34 loci amplified respectively. The loci showed variability within T. rugosa (22 polymorphic loci) and at least 27 loci also exhibited variation among the three species, highlighting the usefulness of these markers for phylogenetic, phylogeographic and population genetic analyses in T. rugosa and related species.

Multilocus phylogeny and recent rapid radiation of the viviparous sea snakes (Elapidae: Hydrophiinae)

The viviparous sea snakes (Hydrophiinae: Hydrophiini) comprise a young but morphologically and ecologically diverse clade distributed throughout the Indo-Pacific. Despite presenting a very promising model for marine diversification studies, many relationships among the 62 species and 16 genera in Hydrophiini remain unresolved. Here, we extend previous taxonomic and genomic sampling for Hydrophiini using three mitochondrial fragments and five nuclear loci for multiple individuals of 39 species in 15 genera. Our results highlight many of the impediments to inferring phylogenies in recent rapid radiations, including low variation at all five nuclear markers, and conflicting relationships supported by mitochondrial and nuclear trees. However, concatenated Bayesian and likelihood analyses, and a multilocus coalescent tree, recovered concordant support for primary clades and several previously unresolved inter-specific groupings. The Aipysurus group is monophyletic, with egg-eating specialists forming separate, early-diverging lineages. All three monotypic semi-aquatic genera (Ephalophis, Parahydrophis and Hydrelaps) are robustly placed as early diverging lineages along the branch leading to the Hydrophis group, with Ephalophis recovered as sister to Parahydrophis. The molecular phylogeny implies extensive evolutionary convergence in feeding adaptations within the Hydrophis group, especially the repeated evolution of small-headed (microcephalic) forms. Microcephalophis (Hydrophis) gracilis is robustly recovered as a relatively distant sister lineage to all other sampled Hydrophis group species, here termed the 'core Hydrophis group'. Within the 'core Hydrophis group', Hydrophis is recovered as broadly paraphyletic, with several other genera nested within it (Pelamis, Enhydrina, Astrotia, Thalassophina, Acalyptophis, Kerilia, Lapemis, Disteira). Instead of erecting multiple new genera, we recommend dismantling the latter (mostly monotypic) genera and recognising a single genus, Hydrophis Latreille 1802, for the core Hydrophis group. Estimated divergence times suggest that all Hydrophiini last shared a common ancestor ∼6million years ago, but that the majority of extant lineages diversified over the last ∼3.5million years. The core Hydrophis group is a young and rapidly speciating clade, with 26 sampled species and 9 genera and dated at only ∼1.5-3million years old.

Rise of the machines--recommendations for ecologists when using next generation sequencing for microsatellite development

Next generation sequencing is revolutionizing molecular ecology by simplifying the development of molecular genetic markers, including microsatellites. Here, we summarize the results of the large-scale development of microsatellites for 54 nonmodel species using next generation sequencing and show that there are clear differences amongst plants, invertebrates and vertebrates for the number and proportion of motif types recovered that are able to be utilized as markers. We highlight that the heterogeneity within each group is very large. Despite this variation, we provide an indication of what number of sequences and consequent proportion of a 454 run are required for the development of 40 designable, unique microsatellite loci for a typical molecular ecological study. Finally, to address the challenges of choosing loci from the vast array of microsatellite loci typically available from partial genome runs (average for this study, 2341 loci), we provide a microsatellite development flowchart as a procedural guide for application once the results of a partial genome run are obtained.

Positive selection in the N-terminal extramembrane domain of lung surfactant protein C (SP-C) in marine mammals

Maximum-likelihood models of codon and amino acid substitution were used to analyze the lung-specific surfactant protein C (SP-C) from terrestrial, semi-aquatic, and diving mammals to identify lineages and amino acid sites under positive selection. Site models used the nonsynonymous/synonymous rate ratio (omega) as an indicator of selection pressure. Mechanistic models used physicochemical distances between amino acid substitutions to specify nonsynonymous substitution rates. Site models strongly identified positive selection at different sites in the polar N-terminal extramembrane domain of SP-C in the three diving lineages: site 2 in the cetaceans (whales and dolphins), sites 7, 9, and 10 in the pinnipeds (seals and sea lions), and sites 2, 9, and 10 in the sirenians (dugongs and manatees). The only semi-aquatic contrast to indicate positive selection at site 10 was that including the polar bear, which had the largest body mass of the semi-aquatic species. Analysis of the biophysical properties that were influential in determining the amino acid substitutions showed that isoelectric point, chemical composition of the side chain, polarity, and hydrophobicity were the crucial determinants. Amino acid substitutions at these sites may lead to stronger binding of the N-terminal domain to the surfactant phospholipid film and to increased adsorption of the protein to the air-liquid interface. Both properties are advantageous for the repeated collapse and reinflation of the lung upon diving and resurfacing and may reflect adaptations to the high hydrostatic pressures experienced during diving.

Host-defence peptide profiles of the skin secretions of interspecific hybrid tree frogs and their parents, female Litoria splendida and male Litoria caerulea

Five healthy adult female first-generation hybrid tree frogs were produced by interspecific breeding of closely related tree frogs Litoria splendida and L. caerulea in a cage containing large numbers of males and females of both species. Phylogenetic analysis of mitochondrial DNA sequences established the female parent to be L. splendida. The peptide profile of the hybrid frogs included the neuropeptide caerulein, four antibiotics of the caerin 1 family and several neuronal nitric oxide synthase inhibitors of the caerin 1 and 2 classes of peptides. The skin secretions of the hybrids contained some peptides common to only one parent, some produced by both parental species, and four peptides expressed by the hybrids but not the parental species.

Comparison of ELISA and RT-PCR for the detection of Prunus necrotic ring spot virus and prune dwarf virus in almond (Prunus dulcis)

A technique based on the reverse transcriptase-polymerase chain reaction (RT-PCR) has been developed to detect the presence of Prunus necrotic ringspot virus (PNRSV) and prune dwarf virus (PDV) simultaneously in almond. This paper presents the results of a 3-year study comparing both enzyme-linked immunosorbent assay (ELISA) and RT-PCR for the detection of PNRSV and PDV using 175 almond leaf samples. Multiplex RT-PCR was found to be more sensitive than ELISA, especially when followed by nested PCR for the detection of PDV. The RT-PCR technique has the added advantage that plant material can be tested at any time throughout the growing season.

Detection of Prunus necrotic ringspot virus in almond: effect of sampling time on the efficiency of serological and biological indexing methodologies

The reliability of enzyme-linked immunosorbent assay (ELISA) and indexing methodologies for the detection of Prunus necrotic ringspot virus (PNRSV) in almond was assessed under local conditions. Thirteen field-grown almond trees were sampled fortnightly throughout the growing season. Petal and leaf homogenates were used for ELISA and to inoculate herbaceous indicator species, while buds collected from October onwards were budded to woody indicators. ELISA reliably detected PNRSV in petals and young leaves from bud-break until the cessation of stem elongation. While herbaceous indexing was not as reliable as ELISA, woody indicators could differentiate positive from negative samples reliably throughout the testing period. However, for mass screening of foundation plantings, nursery stock and elite germplasm, woody indexing is too costly and slow to give results. The use of ELISA can overcome these limitations but the timing of sample collection is critical. Petal or leaf tissue collected early in the season will yield the most reliable results.

Identification of Seed Gall Nematodes of Agronomic and Regulatory Concern with PCR-RFLP of ITS1

A molecular analysis of eight described species of seed gall nematode, along with six undescribed isolates from different hosts, has revealed a strong association between nucleotide sequence polymorphism and host status. Each anguinid nematode associated with a unique host produced a unique PCR-RFLP pattern for the ITS1 region. Anguina species that had been synonymized in the past, Anguina agrostis, A. funesta, and A. wevelli (Afrina wevelli), were readily discriminated. Two undescribed species from northern New South Wales and southeastern South Australia, reported to be vectors of Rathyaibacter toxicus in the disease called ''floodplain staggers,'' were differentiated by a single restriction enzyme, and both could be separated easily from A. funesta, the vector of R. toxicus in annual ryegrass toxicity. Other species differentiated in this study include A. agropyronifloris, A. graminis, A. microlaenae, A. pacificae, and undescribed species from host species Dactylis glomerata, Agrostis avenacea, Polypogon monospeliensis, Stipa sp., Astrebla pectinata, and Holcus lanatus. Phylogenetic analysis of the ITS1 region suggests that considerable anguinid genetic diversification has accompanied specialization on different host species.