Phylogenomic analysis on the exceptionally diverse fish clade Gobioidei (Actinopterygii: Gobiiformes) and data-filtering based on molecular clocklikeness

The complete mitochondrial genome of a Biwa goby, Gymnogobius isaza (Tanaka, 1916)

We determined the complete mitochondrial DNA sequence of a Biwa goby, Gymnogobius isaza (Tanaka, 1916) using next-generation sequencing methods. The composition of its mitogenome is the same as that observed in most other vertebrates, comprising of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and two control regions. Our molecular phylogenetic analysis confirmed the close phylogenetic relationship between G. isaza and G. petschiliensis. This mitogenome information will be useful for distribution surveys using environmental DNA and the development of conservation strategies for this species.

Diversification of the sleepers (Gobiiformes: Gobioidei: Eleotridae) and evolution of the root gobioid families

Eleotridae (sleepers) and five smaller families are the earliest diverging lineages within Gobioidei. Most inhabit freshwaters in and around the Indo-Pacific, but Eleotridae also includes species that have invaded the Neotropics as well as several inland radiations in the freshwaters of Australia, New Zealand, and New Guinea. Previous efforts to infer phylogeny of these families have been based on sets of mitochondrial or nuclear loci and have yielded uncertain resolution of clades within Eleotridae. We expand the taxon sampling of previous studies and use genomic data from nuclear ultraconserved elements (UCEs) to infer phylogeny, then calibrate the hypothesis with recently discovered fossils. Our hypothesis clarifies ambiguously resolved relationships, provides a timescale for divergences, and indicates the core crown Eleotridae diverged over a short period 24.3-26.3 Ma in the late Oligocene. Within Eleotridae, we evaluate diversification dynamics with BAMM and find evidence for an overall slowdown in diversification over the past 35 Ma, but with a sharp increase 3.5 Ma in the genus Mogurnda, a clade of brightly colored species found in the freshwaters of Australia and New Guinea.

Amblyopinae Mitogenomes Provide Novel Insights into the Paraphyletic Origin of Their Adaptation to Mudflat Habitats

The water-to-land transition is one of the most important events in evolutionary history of vertebrates. However, the genetic basis underlying many of the adaptations during this transition remains unclear. Mud-dwelling gobies in the subfamily Amblyopinae are one of the teleosts lineages that show terrestriality and provide a useful system for clarifying the genetic changes underlying adaptations to terrestrial life. Here, we sequenced the mitogenome of six species in the subfamily Amblyopinae. Our results revealed a paraphyletic origin of Amblyopinae with respect to Oxudercinae, which are the most terrestrial fishes and lead an amphibious life in mudflats. This partly explains the terrestriality of Amblyopinae. We also detected unique tandemly repeated sequences in the mitochondrial control region in Amblyopinae, as well as in Oxudercinae, which mitigate oxidative DNA damage stemming from terrestrial environmental stress. Several genes, such as ND2, ND4, ND6 and COIII, have experienced positive selection, suggesting their important roles in enhancing the efficiency of ATP production to cope with the increased energy requirements for life in terrestrial environments. These results strongly suggest that the adaptive evolution of mitochondrial genes has played a key role in terrestrial adaptions in Amblyopinae, as well as in Oxudercinae, and provide new insights into the molecular mechanisms underlying the water-to-land transition in vertebrates.

Molecular phylogenetics of the Clupeiformes based on exon-capture data and a new classification of the order

The Clupeiformes, including among others herrings, anchovies, shads and menhadens are ecologically and commercially important, yet their phylogenetic relationships are still controversial. Previous classification of Clupeiformes were based on morphological characters or lack of synapomorphic characters. More recent studies based on molecular data as well as new morphological evidence are keeping challenging their phylogenetic relations and there is still no consensus on many interrelationships within the Clupeiformes. In this study, we collected nuclear sequence data from 4,434 single-copy protein coding loci using a gene-capture method. We obtained a robust phylogeny based on 1,165 filtered loci with less than 30 % missing data. Our major findings include: 1) reconfirmation of monophyly of the Clupeiformes, that is, Denticipitidae is sister to all other clupeiforms; 2) the polyphyletic nature of dussumieriids and early branching of Spratelloididae from all other clupeoids were confirmed using datasets curated for less missing data and more balanced base composition in the respective taxa. The next branching clade is the monophyletic Engraulidae. Pristigasteridae also is monophyletic, but it was nested in the previously defined "Clupeidae". Within Pristigasteridae there is no support for monophyletic Pelloninae. Chirocentrus is close to Dussumieria and not to engraulids. The miniaturized Sundasalanx is placed close to the ehiravine Clupeonella, however, with a relatively deep split. The genus Clupea, is not part of the diverse "Clupeidae", but part of a clade containing additionally Sprattus and Etrumeus. Within the crown group clades, Alosidae and Dorosomatidae are retrieved as sister clades. Based on new fossil calibration points, we found that major lineages of the clupeiforms diverged in the late Cretaceous and early Paleogene. The extinction event at the end of the Cretaceous may have created ecological niches, which could have fueled the diversification of clupeiform fishes. Based on the strong evidence of the present study, we propose an updated classification of Clupeiformes consisting of ten families: Denticipitidae; Spratelloididae; Engraulidae (Engraulinae + Coiliinae); Clupeidae; Chirocentridae; Dussumieriidae; Pristigasteridae; Ehiravidae; Alosidae, Dorosomatidae.

Rapid evolution fuels transcriptional plasticity to ocean acidification

Ocean acidification (OA) is postulated to affect the physiology, behavior, and life-history of marine species, but potential for acclimation or adaptation to elevated pCO2 in wild populations remains largely untested. We measured brain transcriptomes of six coral reef fish species at a natural volcanic CO2  seep and an adjacent control reef in Papua New Guinea. We show that elevated pCO2 induced common molecular responses related to circadian rhythm and immune system but different magnitudes of molecular response across the six species. Notably, elevated transcriptional plasticity was associated with core circadian genes affecting the regulation of intracellular pH and neural activity in Acanthochromis polyacanthus. Gene expression patterns were reversible in this species as evidenced upon reduction of CO2 following a natural storm-event. Compared with other species, Ac. polyacanthus has a more rapid evolutionary rate and more positively selected genes in key functions under the influence of elevated CO2 , thus fueling increased transcriptional plasticity. Our study reveals the basis to variable gene expression changes across species, with some species possessing evolved molecular toolkits to cope with future OA.

Phylogenetics And The Evolution Of Terrestriality In Mudskippers (Gobiidae: Oxudercinae)

The initial vertebrate conquest of land by stegocephalians (Sarcopterygia) allowed access to new resources and exploitation of untapped niches precipitating a major phylogenetic diversification. However, a paucity of fossils has left considerable uncertainties about phylogenetic relationships and the eco-morphological stages in this key transition in Earth history. Among extant actinopterygians, three genera of mudskippers (Gobiidae: Oxudercinae), Boleophthalmus, Periophthalmus and Periophthalmodon are the most terrestrialized, with vertebral, appendicular, locomotory, respiratory, and epithelial specializations enabling overland excursions up to 14 hours. Unlike early stegocephalians, the ecologies and morphologies of the 45 species of oxudercines are well known, making them viable analogs for the initial vertebrate conquest of land. Nevertheless, they have received little phylogenetic attention. We compiled the largest molecular dataset to date, with 29 oxudercine species, and 5 nuclear and mitochondrial loci. Phylogenetic and comparative analyses revealed strong support for two independent terrestrial transitions, and a complex suit of ecomorphological forms in estuarine environments. Furthermore, neither Oxudercinae nor their presumed sister-group the eel gobies (Amblyopinae, a group of elongated gobies) were monophyletic with respect to each other, requiring a merging of these two subfamilies and revealing an expansion of phenotypic variation within the "mudskipper" clade. We did not find support for the expected linear model of ecomorphological and locomotory transition from fully aquatic, to mudswimming, to pectoral-aided mudswimming, to lobe-finned terrestrial locomotion proposed by earlier morphological studies. This high degree of convergent or parallel transitions to terrestriality, and apparent divergent directions of estuarine adaptation, promises even greater potential for this clade to illuminate the conquest of land. Future work should focus on these less-studied species with "transitional" and other mud-habitat specializations to fully resolve the dynamics of this diversification.

Gonad morphology of Rhyacichthys aspro (Valenciennes, 1837), and the diagnostic reproductive morphology of gobioid fishes

Rhyacichthys aspro is a “basal” taxon in the Suborder Gobioidei of the teleost order Gobiiformes. We provide detailed descriptions of the reproductive morphology of adult males and females to assess the diagnostic reproductive morphological characters of this speciose clade of bony fishes. Female R. aspro are asynchronous spawners: they are able to spawn more than once in a breeding season. Oocytes are inferred to have short attachment filaments. A conspicuous feature of the external anatomy of the reproductive system (RSy) of female R. aspro is an ornate fimbriate pad upon which the urogenital papilla rests. The male reproductive system is characterized by an intralobar collection system in both the testicular and secretory lobes, termed the “sperm‐collecting canal” and “milt‐collecting canal,” respectively. These may provide additional storage for sperm and milt. The spermatogenic lobe, or testis, is that portion of the male gobioid RSy comprising seminiferous lobules and separate from other RSy components. The secretory lobe is that portion of the male gobioid reproductive system that consists of secretory lobules and is separated from other components of the male RSy. The secretory lobe has also been called, in English, the sperm‐duct gland, accessory gonadal structure, or seminal vesicle, and is endorsed as a synapomorphy of gobioid fishes.

Correlation between acoustic divergence and phylogenetic distance in soniferous European gobiids (Gobiidae; Gobius lineage)

In fish, species identity can be encoded by sounds, which have been thoroughly investigated in European gobiids (Gobiidae, Gobius lineage). Recent evolutionary studies suggest that deterministic and/or stochastic forces could generate acoustic differences among related animal species, though this has not been investigated in any teleost group to date. In the present comparative study, we analysed the sounds from nine soniferous gobiids and quantitatively assessed their acoustic variability. Our interspecific acoustic study, incorporating for the first time the representative acoustic signals from the majority of soniferous gobiids, suggested that their sounds are truly species-specific (92% of sounds correctly classified into exact species) and each taxon possesses a unique set of spectro-temporal variables. In addition, we reconstructed phylogenetic relationships from a concatenated molecular dataset consisting of multiple molecular markers to track the evolution of acoustic signals in soniferous gobiids. The results of this study indicated that the genus Padogobius is polyphyletic, since P. nigricans was nested within the Ponto-Caspian clade, while the congeneric P. bonelli turned out to be a sister taxon to the remaining investigated soniferous species. Lastly, by extracting the acoustic and genetic distance matrices, sound variability and genetic distance were correlated for the first time to assess whether sound evolution follows a similar phylogenetic pattern. The positive correlation between the sound variability and genetic distance obtained here emphasizes that certain acoustic features from representative sounds could carry the phylogenetic signal in soniferous gobiids. Our study was the first attempt to evaluate the mutual relationship between acoustic variation and genetic divergence in any teleost fish.

Exon-capture data and locus screening provide new insights into the phylogeny of flatfishes (Pleuronectoidei)

There is an extensive collection of literature on the taxonomy and phylogenetics of flatfishes (Pleuronectiformes) that extends over two centuries, but consensus on many of their evolutionary relationships remains elusive. Phylogenetic uncertainty stems from highly divergent results derived from morphological and genetic characters, and between various molecular datasets. Deciphering relationships is complicated by rapid diversification early in the Pleuronectiformes tree and an abundance of studies that incompletely and inconsistently sample taxa and genetic markers. We present phylogenies based on a genome-wide dataset (4,434 nuclear markers via exon-capture) and wide taxon sampling (86 species spanning 12 of 16 families) of the largest flatfish suborder (Pleuronectoidei). Nine different subsets of the data and two tree construction approaches (eighteen phylogenies in total) are remarkably consistent with other recent molecular phylogenies, and show strong support for the monophyly of all families included except Pleuronectidae. Analyses resolved a novel phylogenetic hypothesis for the family Rhombosoleidae as being within the Pleuronectoidea rather than the Soleoidea, and failed to support the subfamily Hippoglossinae as a monophyletic group. Our results were corroborated with evidence from previous phylogenetic studies to outline regions of persistent phylogenetic uncertainty and identify groups in need of further phylogenetic inference.

Excluding loci with substitution saturation improves inferences from phylogenomic data

The historical signal in nucleotide sequences becomes eroded over time by substitutions occurring repeatedly at the same sites. This phenomenon, known as substitution saturation, is recognized as one of the primary obstacles to deep-time phylogenetic inference using genome-scale data sets. We present a new test of substitution saturation and demonstrate its performance in simulated and empirical data. For some of the 36 empirical phylogenomic data sets that we examined, we detect substitution saturation in around 50% of loci. We found that saturation tends to be flagged as problematic in loci with highly discordant phylogenetic signals across sites. Within each data set, the loci with smaller numbers of informative sites are more likely to be flagged as containing problematic levels of saturation. The entropy saturation test proposed here is sensitive to high evolutionary rates relative to the evolutionary timeframe, while also being sensitive to several factors known to mislead phylogenetic inference, including short internal branches relative to external branches, short nucleotide sequences, and tree imbalance. Our study demonstrates that excluding loci with substitution saturation can be an effective means of mitigating the negative impact of multiple substitutions on phylogenetic inferences. [Phylogenetic model performance; phylogenomics; substitution model; substitution saturation; test statistics.]

Phylogenomics reveals accelerated late Cretaceous diversification of bee flies (Diptera: Bombyliidae)

Bombyliidae is a very species-rich and widespread family of parasitoid flies with more than 250 genera classified into 17 extant subfamilies. However, little is known about their evolutionary history or how their present-day diversity was shaped. Transcriptomes of 15 species and anchored hybrid enrichment (AHE) sequence captures of 86 species, representing 94 bee fly species and 14 subfamilies, were used to reconstruct the phylogeny of Bombyliidae. We integrated data from transcriptomes across each of the main lineages in our AHE tree to build a data set with more genes (550 loci versus 216 loci) and higher support levels. Our overall results show strong congruence with the current classification of the family, with 11 out of 14 included subfamilies recovered as monophyletic. Heterotropinae and Mythicomyiinae are successive sister groups to the remainder of the family. We examined the evolution of key morphological characters through our phylogenetic hypotheses and show that neither the "sand chamber subfamilies" nor the "Tomophthalmae" are monophyletic in our phylogenomic analyses. Based on our results, we reinstate two tribes at the subfamily level (Phthiriinae stat. rev. and Ecliminae stat. rev.) and we include the genus Sericosoma Macquart (previously incertae sedis) in the subfamily Oniromyiinae, bringing the total number of bee fly subfamilies to 19. Our dating analyses indicate a Jurassic origin of the family (165-194 Ma), with the sand chamber evolving early in bee fly evolution, in the late Jurassic or mid-Cretaceous (100-165 Ma). We hypothesize that the angiosperm radiation and the hothouse climate established during the late Cretaceous accelerated the diversification of bee flies, by providing an expanded range of resources for the parasitoid larvae and nectarivorous adults.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Comparing Ultraconserved Elements and Exons for Phylogenomic Analyses of Middle American Cichlids: When Data Agree to Disagree

Choosing among types of genomic markers to be used in a phylogenomic study can have a major influence on the cost, design, and results of a study. Yet few attempts have been made to compare categories of next-generation sequence markers limiting our ability to compare the suitability of these different genomic fragment types. Here, we explore properties of different genomic markers to find if they vary in the accuracy of component phylogenetic trees and to clarify the causes of conflict obtained from different data sets or inference methods. As a test case, we explore the causes of discordance between phylogenetic hypotheses obtained using a novel data set of ultraconserved elements (UCEs) and a recently published exon data set of the cichlid tribe Heroini. Resolving relationships among heroine cichlids has historically been difficult, and the processes of colonization and diversification in Middle America and the Greater Antilles are not yet well understood. Despite differences in informativeness and levels of gene tree discordance between UCEs and exons, the resulting phylogenomic hypotheses generally agree on most relationships. The independent data sets disagreed in areas with low phylogenetic signal that were overwhelmed by incomplete lineage sorting and nonphylogenetic signals. For UCEs, high levels of incomplete lineage sorting were found to be the major cause of gene tree discordance, whereas, for exons, nonphylogenetic signal is most likely caused by a reduced number of highly informative loci. This paucity of informative loci in exons might be due to heterogeneous substitution rates that are problematic to model (i.e., computationally restrictive) resulting in systematic errors that UCEs (being less informative individually but more uniform) are less prone to. These results generally demonstrate the robustness of phylogenomic methods to accommodate genomic markers with different biological and phylogenetic properties. However, we identify common and unique pitfalls of different categories of genomic fragments when inferring enigmatic phylogenetic relationships.

Evolutionary Rate Variation Among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference

Phylogenetic analyses of genomic data provide a powerful means of reconstructing the evolutionary relationships among organisms, yet such analyses are often hindered by conflicting phylogenetic signals among loci. Identifying the signals that are most influential to species-tree estimation can help to inform the choice of data for phylogenomic analysis. We investigated this in an analysis of 30 phylogenomic data sets. For each data set, we examined the association between several branch-length characteristics of gene trees and the distance between these gene trees and the corresponding species trees. We found that the distance of each gene tree to the species tree inferred from the full data set was positively associated with variation in root-to-tip distances and negatively associated with mean branch support. However, no such associations were found for gene-tree length, a measure of the overall substitution rate at each locus. We further explored the usefulness of the best-performing branch-based characteristics for selecting loci for phylogenomic analyses. We found that loci that yield gene trees with high variation in root-to-tip distances have a disproportionately distant signal of tree topology compared with the complete data sets. These results suggest that rate variation across lineages should be taken into consideration when exploring and even selecting loci for phylogenomic analysis.[Branch support; data filtering; nucleotide substitution model; phylogenomics; substitution rate; summary coalescent methods.]

Exon probe sets and bioinformatics pipelines for all levels of fish phylogenomics

Exon markers have a long history of use in phylogenetics of ray-finned fishes, the most diverse clade of vertebrates with more than 35,000 species. As the number of published genomes increases, it has become easier to test exons and other genetic markers for signals of ancient duplication events and filter out paralogues that can mislead phylogenetic analysis. We present seven new probe sets for current target-capture phylogenomic protocols that capture 1,104 exons explicitly filtered for paralogues using gene trees. These seven probe sets span the diversity of teleost fishes, including four sets that target five hyperdiverse percomorph clades which together comprise ca. 17,000 species (Carangaria, Ovalentaria, Eupercaria, and Syngnatharia + Pelagiaria combined). We additionally included probes to capture legacy nuclear exons and mitochondrial markers that have been commonly used in fish phylogenetics (despite some exons being flagged for paralogues) to facilitate integration of old and new molecular phylogenetic matrices. We tested these probes experimentally for 56 fish species (eight species per probe set) and merged new exon-capture sequence data into an existing data matrix of 1,104 exons and 300 ray-finned fish species. We provide an optimized bioinformatics pipeline to assemble exon capture data from raw reads to alignments for downstream analysis. We show that legacy loci with known paralogues are at risk of assembling duplicated sequences with target-capture, but we also assembled many useful orthologous sequences that can be integrated with many PCR-generated matrices. These probe sets are a valuable resource for advancing fish phylogenomics because targeted exons can easily be extracted from increasingly available whole genome and transcriptome data sets, and also may be integrated with existing PCR-based exon and mitochondrial data.

Migratory flexibility in native Hawai'ian amphidromous fishes

We assessed the prevalence of life history variation across four of the five native amphidromous Hawai'ian gobioids to determine whether some or all exhibit evidence of partial migration. Analysis of otolith Sr.: Ca concentrations affirmed that all are amphidromous and revealed evidence of partial migration in three of the four species. We found that 25% of Lentipes concolor (n = 8), 40% of Eleotris sandwicensis (n = 20) and 29% of Stenogobius hawaiiensis (n = 24) did not exhibit a migratory life-history. In contrast, all individuals of Sicyopterus stimpsoni (n = 55) included in the study went to sea as larvae. Lentipes concolor exhibited the shortest mean larval duration (LD) at 87 days, successively followed by E. sandwicensis (mean LD = 102 days), S. hawaiiensis (mean LD = 114 days) and S. stimpsoni (mean LD = 120 days). These findings offer a fresh perspective on migratory life histories that can help improve efforts to conserve and protect all of these and other at-risk amphidromous species that are subject to escalating anthropogenic pressures in both freshwater and marine environments.

Assexon: Assembling Exon Using Gene Capture Data

Exon capture across species has been one of the most broadly applied approaches to acquire multi-locus data in phylogenomic studies of non-model organisms. Methods for assembling loci from short-read sequences (eg, Illumina platforms) that rely on mapping reads to a reference genome may not be suitable for studies comprising species across a wide phylogenetic spectrum; thus, de novo assembling methods are more generally applied. Current approaches for assembling targeted exons from short reads are not particularly optimized as they cannot (1) assemble loci with low read depth, (2) handle large files efficiently, and (3) reliably address issues with paralogs. Thus, we present Assexon: a streamlined pipeline that de novo assembles targeted exons and their flanking sequences from raw reads. We tested our method using reads from Lepisosteus osseus (4.37 Gb) and Boleophthalmus pectinirostris (2.43 Gb), which are captured using baits that were designed based on genome sequence of Lepisosteus oculatus and Oreochromis niloticus, respectively. We compared performance of Assexon to PHYLUCE and HybPiper, which are commonly used pipelines to assemble ultra-conserved element (UCE) and Hyb-seq data. A custom exon capture analysis pipeline (CP) developed by Yuan et al was compared as well. Assexon accurately assembled more than 3400 to 3800 (20%-28%) loci than PHYLUCE and more than 1900 to 2300 (8%-14%) loci than HybPiper across different levels of phylogenetic divergence. Assexon ran at least twice as fast as PHYLUCE and HybPiper. Number of loci assembled using CP was comparable with Assexon in both tests, while Assexon ran at least 7 times faster than CP. In addition, some steps of CP require the user’s interaction and are not fully automated, and this user time was not counted in our calculation. Both Assexon and CP retrieved no paralogs in the testing runs, but PHYLUCE and Hybpiper did. In conclusion, Assexon is a tool for accurate and efficient assembling of large read sets from exon capture experiments. Furthermore, Assexon includes scripts to filter poorly aligned coding regions and flanking regions, calculate summary statistics of loci, and select loci with reliable phylogenetic signal. Assexon is available at https://github.com/yhadevol/Assexon.

The complete mitochondrial genome characterization of Paratrypauchen microcephalus (Gobiiformes: Oxudercidae) and phylogenetic consideration

The complete mitochondrial genome of Paratrypauchen microcephalus (Gobiiformes: Oxudercidae) was completely sequenced by high throughput sequencing method. The complete mitochondrial genome was 16,552 bp in length, consisted of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, a putative control region (CR), and an origin of replication on the light-strand (OL). The base composition values for the mitochondrial genome were 29.1%, 27.6%, 15.7%, and 27.6% for A, C, G, and T, respectively. The gene arrangement is identical to those in typical fishes. Phylogenetic tree based on 13 protein-coding genes shows that P. microcephalus has a close phylogenetic relationship with genus Trypauchen and belongs to Oxudercidae.

Multiple molecular markers reinforce the systematic framework of unique Australian cave fishes (Milyeringa : Gobioidei)

Australia was once thought to be a biodiversity desert when considering the subterranean world; however, recent work has revealed a fascinating collection of cave creatures, many with surprising biogeographic histories. This has especially been so in the karstic regions of north-western Australia (Cape Range peninsula, Barrow Island, Pilbara), which is home not only to a diverse collection of subterranean invertebrates, but also to the continent’s only known underworld-adapted vertebrates, which includes the cave fish in the genus Milyeringa. These cave gudgeons have recently been in a state of taxonomic flux, with species being both split and lumped, but this was done with limited data (incomplete geographic sampling and no nuclear DNA sequence data). Therefore, we have revisited the systematic status of Milyeringa in a total-evidence molecular approach by integrating all existing data (mitochondrial, allozymes) with new DNA sequences from nuclear and mitochondrial loci and new multilocus allozyme data. Our conclusion, that there are two species, matches the most recent taxonomic treatment, with Milyeringa veritas present on both the eastern and western sides of the Cape Range peninsula, and Milyeringa justitia on Barrow Island. This has implications for future research in the linked fields of biogeography and conservation.