Cryptic biodiversity in the freshwater fishes of the Kimberley endemism hotspot, northwestern Australia

Cryptic species, biogeography, and patterns of introgression in the fish genus Mogurnda (Eleotridae) from the Australian wet tropics: A purple patch for purple-spots

Accurately delimiting species is an essential first step towards understanding the true biodiversity of an ecosystem and any subsequent efforts to identify and protect taxa at risk of extinction. Current molecular evidence suggests that purple-spotted gudgeons (genus Mogurnda) harbour high levels of cryptic biodiversity across their broad distributional range. The present study uses a large single nucleotide polymorphism (SNP) dataset plus a companion allozyme dataset to clarify taxonomic uncertainty, patterns of introgression, and biogeographic relationships among Mogurnda populations within the Queensland Wet Tropics (QWT), a known biodiversity hotspot. Both datasets were strongly concordant in identifying a total of seven taxa split among distinct southern, northern, and lowlands groups. No two taxa were found in strict sympatry, but many appear to be parapatric and occur within the same drainage basin. Although clear evidence of introgression was only evident at six sites (∼4%), the genomic signature of modest historic admixture between proximally-distributed taxa was detected at multiple other sites. Nevertheless, all primary genetic and phylogenetic analyses strongly supported the integrity and diagnosability of these seven taxa. We therefore nominate these as novel candidate species for what appears to be yet another hyper-cryptic complex within the Australian freshwater ichthyological fauna. These results offer up intriguing ecological scenarios and conservation implications for multiple candidate species with narrow ranges in specialised habitat. We conclude by exploring the major biogeographic patterns displayed by QWT Mogurnda.

Cryptic biodiversity of freshwater fish species in Bangladesh

Unrecognized cryptic species impede conservation planning and biodiversity assessments. DNA barcoding has tremendously expanded the number of novel and cryptic species in biological science. Despite few sporadic studies, the exact number of freshwater species found in Bangladesh is not known. To assess this biodiversity, we sequenced the COI gene of 124 freshwater specimens, which were gathered from various localities around Bangladesh. Seven cryptic species hidden among the currently studied specimens were identified based on the findings of phylogenetic and species delimitation analyses. The preliminary assessment also encompassed a restricted morphological examination of these cryptic taxa. The appearance of cryptic species, some of them possibly endemic, has been hypothesized. This raises concerns regarding the true diversity and evolutionary history of freshwater species in Bangladesh, which are significantly underrepresented in the current systematic frameworks that do not account for DNA data. Our current study provides baseline data that might aid local ichthyologists in their quest to identify additional new species by combining several variables (morphology and ecology). Further research is warranted to protect the priceless freshwater species in Bangladesh.

Museum Skins Enable Identification of Introgression Associated with Cytonuclear Discordance

Increased sampling of genomes and populations across closely related species has revealed that levels of genetic exchange during and after speciation are higher than previously thought. One obvious manifestation of such exchange is strong cytonuclear discordance, where the divergence in mitochondrial DNA (mtDNA) differs from that for nuclear genes more (or less) than expected from differences between mtDNA and nuclear DNA (nDNA) in population size and mutation rate. Given genome-scale data sets and coalescent modeling, we can now confidently identify cases of strong discordance and test specifically for historical or recent introgression as the cause. Using population sampling, combining exon capture data from historical museum specimens and recently collected tissues we showcase how genomic tools can resolve complex evolutionary histories in the brachyotis group of rock-wallabies (Petrogale). In particular, applying population and phylogenomic approaches we can assess the role of demographic processes in driving complex evolutionary patterns and assess a role of ancient introgression and hybridization. We find that described species are well supported as monophyletic taxa for nDNA genes, but not for mtDNA, with cytonuclear discordance involving at least 4 operational taxonomic units across 4 species which diverged 183-278 kya. ABC modeling of nDNA gene trees supports introgression during or after speciation for some taxon pairs with cytonuclear discordance. Given substantial differences in body size between the species involved, this evidence for gene flow is surprising. Heterogenous patterns of introgression were identified but do not appear to be associated with chromosome differences between species. These and previous results suggest that dynamic past climates across the monsoonal tropics could have promoted reticulation among related species.

A revision of the gudgeon genus Hypseleotris (Gobiiformes: Gobioidei: Eleotridae) of northwest Australia, describing three new species and synonymizing the genus Kimberleyeleotris

Species within the northwest Australian clade of Hypseleotris (six species) and the genus Kimberleyeleotris (two species) are reviewed following the recording of new populations in the region and a molecular study of the group that identified three undescribed candidate species. Based on the analysis of extensive morphological and nuclear and mitochondrial molecular datasets, Kimberleyeleotris is here formally synonymised with Hypseleotris. Furthermore, three species from the Kimberley region, Western Australia, are described to science: Hypseleotris maranda sp. nov., Hypseleotris wunduwala sp. nov., and Hypseleotris garawudjirri sp. nov. The presence of, or number of scales across the head and body, the pattern of sensory papillae on the head, fin ray counts, dorsal and anal fin colouration (particularly in breeding males), and body depth, can be used to distinguish the members of the northwest Australia lineage. Furthermore, the newly described species were genetically separated from all northwest Australian congeners by K2P distances ranging from 7.8-11.3% based on the CO1 gene, and 7.7-16.3 % based on the entire mitochondrial genome. Two of the new species, H. maranda sp. nov. and H. wunduwala sp. nov., have extremely narrow ranges being found in single sub-catchments of the Roe and King Edward Rivers respectively. On the other hand, H. garawudjirri sp. nov. is moderately widespread, being found across the Charnley, Calder, and Sale rivers. While the conservation risk to H. maranda sp. nov. and H. wunduwala sp. nov. is inherently high due to their small range, there are currently no obvious local threatening processes to either of these species given their remote locations that are little impacted by human activities.

Plotting for change: an analytical framework to aid decisions on which lineages are candidate species in phylogenomic species discovery

A recent study argued that coalescent-based models of species delimitation mostly delineate population structure, not species, and called for the validation of candidate species using biological information additional to the genetic information, such as phenotypic or ecological data. Here, we introduce a framework to interrogate genomic datasets and coalescent-based species trees for the presence of candidate species in situations where additional biological data are unavailable, unobtainable or uninformative. For de novo genomic studies of species boundaries, we propose six steps: (1) visualize genetic affinities among individuals to identify both discrete and admixed genetic groups from first principles and to hold aside individuals involved in contemporary admixture for independent consideration; (2) apply phylogenetic techniques to identify lineages; (3) assess diagnosability of those lineages as potential candidate species; (4) interpret the diagnosable lineages in a geographical context (sympatry, parapatry, allopatry); (5) assess significance of difference or trends in the context of sampling intensity; and (6) adopt a holistic approach to available evidence to inform decisions on species status in the difficult cases of allopatry. We adopt this approach to distinguish candidate species from within-species lineages for a widespread species complex of Australian freshwater fishes (Retropinna spp.). Our framework addresses two cornerstone issues in systematics that are often not discussed explicitly in genomic species discovery: diagnosability and how to determine it, and what criteria should be used to decide whether diagnosable lineages are conspecific or represent different species.

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.

Gene Flow Increases Phylogenetic Structure and Inflates Cryptic Species Estimations: A Case Study on Widespread Philippine Puddle Frogs (Occidozyga laevis)

In cryptic amphibian complexes, there is a growing trend to equate high levels of genetic structure with hidden cryptic species diversity. Typically, phylogenetic structure and distance-based approaches are used to demonstrate the distinctness of clades and justify the recognition of new cryptic species. However, this approach does not account for gene flow, spatial, and environmental processes that can obfuscate phylogenetic inference and bias species delimitation. As a case study, we sequenced genome-wide exons and introns to evince the processes that underlie the diversification of Philippine Puddle Frogs-a group that is widespread, phenotypically conserved, and exhibits high levels of geographically-based genetic structure. We showed that widely adopted tree- and distance-based approaches inferred up to 20 species, compared to genomic analyses that inferred an optimal number of five distinct genetic groups. Using a suite of clustering, admixture, and phylogenetic network analyses, we demonstrate extensive admixture among the five groups and elucidate two specific ways in which gene flow can cause overestimations of species diversity: (1) admixed populations can be inferred as distinct lineages characterized by long branches in phylograms; and (2) admixed lineages can appear to be genetically divergent, even from their parental populations when simple measures of genetic distance are used. We demonstrate that the relationship between mitochondrial and genome-wide nuclear p-distances is decoupled in admixed clades, leading to erroneous estimates of genetic distances and, consequently, species diversity. Additionally, genetic distance was also biased by spatial and environmental processes. Overall, we showed that high levels of genetic diversity in Philippine Puddle Frogs predominantly comprise metapopulation lineages that arose through complex patterns of admixture, isolation-by-distance, and isolation-by-environment as opposed to species divergence. Our findings suggest that speciation may not be the major process underlying the high levels of hidden diversity observed in many taxonomic groups and that widely-adopted tree- and distance-based methods overestimate species diversity in the presence of gene flow.

A fossil-calibrated time-tree of all Australian freshwater fishes

Australian freshwater fishes are a relatively species-poor assemblage, mostly comprising groups derived from older repeated freshwater invasions by marine ancestors, plus a small number of Gondwanan lineages. These taxa are both highly endemic and highly threatened, but a comprehensive phylogeny for Australian freshwater fishes is lacking. This has hampered efforts to study their phylogenetic diversity, distribution of extinction risk, speciation rates, and rates of trait evolution. Here, we present a comprehensive dated phylogeny of 412 Australian fishes. We include all formally recognized freshwater species plus a number of genetically distinct subpopulations, species awaiting formal description, and predominantly brackish-water species that sometimes enter fresh water. The phylogeny was inferred using maximum-likelihood analysis of a multilocus data set comprising six mitochondrial and three nuclear genes from 326 taxa. We inferred the evolutionary timescale using penalized likelihood, then used a statistical approach to add 86 taxa for which no molecular data were available. The time-tree inferred in our study will provide a useful resource for macroecological studies of Australian freshwater fishes by enabling corrections for phylogenetic non-independence in evolutionary and ecological comparative analyses.

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.

A revision of the bioregionalisation of freshwater fish communities in the Australian Monsoonal Tropics

The Australian freshwater fish fauna is very unique, but poorly understood. In the Australian Monsoonal Tropics (AMT) biome of northern Australia, the number of described and candidate species has nearly doubled since the last attempt to analyse freshwater fish species composition patterns and determine a bioregionalisation scheme. Here, we utilise the most complete database of catchment-scale freshwater fish distributions from the AMT to date to: (a) reanalyze spatial patterns of species richness, endemism and turnover of freshwater fishes; (b) propose a biogeographic regionalisation based on species turnover; (c) assess the relationship between species turnover and patterns of environmental change and historic drainage connectivity; and (d) identify sampling gaps. Biogeographic provinces were identified using an agglomerative cluster analysis of a Simpson's beta (β sim) dissimilarity matrix. A generalised dissimilarity model incorporating eighteen environmental variables was used to investigate the environmental correlates of species turnover. Observed and estimated species richness and endemism were calculated and inventory completeness was estimated based on the ratio of observed to estimated species richness. Three major freshwater fish biogeographic provinces and 14 subprovinces are proposed. These differ substantially from the current bioregionalisation scheme. Species turnover was most strongly influenced by environmental variables that are interpreted to reflect changes in terrain (catchment relief and confinement), geology and climate (runoff perenniality, stream density), and biotic responses to climate (net primary productivity). Past connectivity between rivers during low sea-level events is also influential highlighting the importance of historical processes in explaining contemporary patterns of biodiversity in the AMT. The inclusion of 49 newly discovered species and candidate species only reinforced known focal points of species richness and endemism in the AMT. However, a number of key sampling gaps remain that need to be filled to fully characterise the proposed bioregionalisation.