Whole-genome shotgun sequence assembly enables rapid gene characterization in the tropical fish barramundi, Lates calcarifer

[Extract] Source description: Barramundi, Lates calcarifer, is a commercially important fish farmed throughout Australia and South-East Asia. Despite an increasing availability of genetic resources for the species (e.g. microsatellite and SNP markers, linkage and BAC-based maps and transcriptomic assemblies), the complete characterization of genes is still reliant on laborious molecular methods (e.g. genome walking/RACE PCR cloning/Sanger sequencing)

Populations genetically rifting within a complex geological system: the case of strong structure and low genetic diversity in the migratory freshwater catfish, Bagrus docmak, in East Africa

The complex geological history of East Africa has been a driving factor in the rapid evolution of teleost biodiversity. While there is some understanding of how macroevolutionary drivers have shaped teleost speciation in East Africa, there is a paucity of research into how the same biogeographical factors have affected microevolutionary processes within lakes and rivers. To address this deficiency, population genetic diversity, demography, and structure were investigated in a widely distributed and migratory (potamodromous) African teleost species, Ssemutundu (Bagrus docmak). Samples were acquired from five geographical locations in East Africa within two major drainage basins; the Albertine Rift and Lake Victoria Basin. Individuals (N = 175) were genotyped at 12 microsatellite loci and 93 individuals sequenced at the mitochondrial DNA control region. Results suggested populations from Lakes Edward and Victoria had undergone a severe historic bottleneck resulting in very low nucleotide diversity (π = 0.004 and 0.006, respectively) and negatively significant Fu values (−3.769 and −5.049; p < .05). Heterozygosity deficiencies and restricted effective population size (NeLD) suggested contemporary exposure of these populations to stress, consistent with reports of the species decline in the East African Region. High genetic structuring between drainages was detected at both historical (ɸST = 0.62 for mtDNA; p < .001) and contemporary (microsatellite FST = 0.460; p < .001) levels. Patterns of low genetic diversity and strong population structure revealed are consistent with speciation patterns that have been linked to the complex biogeography of East Africa, suggesting that these biogeographical features have operated as both macro- and micro-evolutionary forces in the formation of the East African teleost fauna

Microbiome diversity and composition varies across body areas in a freshwater turtle

There is increasing recognition that microbiomes are important for host health and ecology, and understanding host microbiomes is important for planning appropriate conservation strategies. However, microbiome data are lacking for many taxa, including turtles. To further our understanding of the interactions between aquatic microbiomes and their hosts, we used next generation sequencing technology to examine the microbiomes of the Krefft's river turtle (Emydura macquarii krefftii). We examined the microbiomes of the buccal (oral) cavity, skin on the head, parts of the shell with macroalgae and parts of the shell without macroalgae. Bacteria in the phyla Proteobacteria and Bacteroidetes were the most common in most samples (particularly buccal samples), but Cyanobacteria, Deinococcus-thermus and Chloroflexi were also common (particularly in external microbiomes). We found significant differences in community composition among each body area, as well as significant differences among individuals. The buccal cavity had lower bacterial richness and evenness than any of the external microbiomes, and it had many amplicon sequence variants (ASVs) with a low relative abundance compared to other body areas. Nevertheless, the buccal cavity also had the most unique ASVs. Parts of the shell with and without algae also had different microbiomes, with particularly obvious differences in the relative abundances of the families Methylomonaceae, Saprospiraceae and Nostocaceae. This study provides novel, baseline information about the external microbiomes of turtles and is a first step in understanding their ecological roles

Conservation genomics reveals fine-scale population structuring and recent declines in the Critically Endangered Australian Kuranda Treefrog

The Kuranda Treefrog occurs in tropical north-east Australia and is listed as Critically Endangered due to its small distribution and population size, with observed declines due to drought and human-associated impacts to habitat. Field surveys identified marked population declines in the mid-2000s, culminating in very low abundance at most sites in 2005 and 2006, followed by limited recovery. Here, samples from before (2001–2004) and after (2007–2009) this decline were analysed using 7132 neutral genome-wide SNPs to assess genetic connectivity among breeding sites, genetic erosion, and effective population size. We found a high level of genetic connectivity among breeding sites, but also structuring between the population at the eastern end of the distribution (Jumrum Creek) versus all other sites. Despite finding no detectable sign of genetic erosion between the two times periods, we observed a marked decrease in effective population size (Ne), from 1720 individuals pre-decline to 818 post-decline. This mirrors the decline detected in the field census data, but the magnitude of the decline suggested by the genetic data is greater. We conclude that the current effective population size for the Kuranda Treefrog remains around 800 adults, split equally between Jumrum Creek and all other sites combined. The Jumrum Creek habitat requires formal protection. Connectivity among all other sites must be maintained and improved through continued replanting of rainforest, and it is imperative that impacts to stream flow and water quality are carefully managed to maintain or increase population sizes and prevent genetic erosion

Gene flow and genetic structure in Nile perch, Lates niloticus, from African freshwater rivers and lakes

Background Geological evolution of the African continent has been subject to complex processes including uplift, volcanism, desert formation and tectonic rifting. This complex geology has created substantial biogeographical barriers, and coupled with anthropogenic introductions of freshwater fishes, has influenced the genetic diversity, connectivity and sub-structuring of the teleost fauna. Nile perch, Lates niloticus, is an iconic fish in Africa and is of high commercial importance, both in the species' native range and where it has been translocated. However, the species is in decline and there is a need to understand its population genetic structure to facilitate sustainable management of the fishery and aquaculture development. Methodology Nile perch tissue samples were acquired from two West and four East (Lakes; Albert, Kyoga, Victoria and Turkana) African locations. Nineteen polymorphic microsatellite loci were used to study the genetic variation among populations across regions (West and East Africa), as well as between native and introduced environments within East Africa. Principal findings and their significance Results revealed strong and significant genetic structuring among populations across the sampled distribution (divergence across regions, FCT = 0.26, P = 0.000). STRUCTURE analysis at a broad scale revealed K = 2 clusters, the West African individuals were assigned to one cluster, while all individuals from the East African region, regardless of whether native or introduced, were assigned to another cluster. The distinct genetic clusters identified in the current study between the West and East African Nile perch, appear to have been maintained by presence of biogeographic barriers and restricted gene flow between the two regions. Therefore, any translocations of Nile perch should be carefully considered across the regions of West and East Africa. Further analysis at a regional scale revealed further structuring of up to K = 3 genetic clusters in East African Nile perch. Significantly (P < 0.05) lower genetic diversity based on analysis of allelic richness (AR) was obtained for the two translocated populations of Lake Kyoga (AR = 3.61) and Lake Victoria (AR = 3.52), compared to Nile perch populations from their putative origins of Lakes Albert (AR = 4.12) and Turkana (AR = 4.43). The lower genetic diversity in the translocated populations may be an indication of previous bottlenecks and may also indicate a difficulty for these populations to persist and adapt to climatic changes and anthropogenic pressures that are currently present in the East African region

Genome-wide comparisons reveal evidence for a species complex in the black-lip pearl oyster Pinctada margaritifera (Bivalvia: Pteriidae)

Evolutionary relationships in the black-lip pearl oyster Pinctada margaritifera which is highly valued for pearl production remain poorly understood. This species possesses an 18,000 km Indo-Pacific natural distribution, and its current description includes six subspecies defined exclusively on morphological characters. To evaluate its taxonomic identity using molecular data, 14 populations in both the Indian and Pacific Oceans (n = 69), and the congeneric taxa P. maxima and P. mazatlanica (n = 29 and n = 10, respectively) were sampled. Phylogenomic reconstruction was carried out using both 8,308 genome-wide SNPs and 10,000 dominant loci (DArTseq PAVs). Reconstructions using neighbour-joining (Nei’s 1972 distances), maximum likelihood and Bayesian approaches all indicate that the taxonomy of P. margaritifera is quite complex, with distinct evolutionary significant units (ESUs) identified within Tanzanian and Iranian populations. Contrastingly, phylogenies generated for Pacific Ocean oysters resolved a large monophyletic clade, suggesting little support for two current morphological subspecies classifications. Furthermore, P. mazatlanica formed a basal clade closest to French Polynesian P. margaritifera, suggesting it may be conspecific. Collectively, these findings provide evidence that P. margaritifera comprises a species complex, perhaps as a result of population fragmentation and increased divergence at range limits

Comparing genomic signatures of selection between the Abbassa Strain and eight wild populations of Nile tilapia (Oreochromis niloticus) in Egypt

Domestication to captive rearing conditions, along with targeted selective breeding have genetic consequences that vary from those in wild environments. Nile tilapia (Oreochromis niloticus) is one of the most translocated and farmed aquaculture species globally, farmed throughout Asia, North and South America, and its African native range. In Egypt, a breeding program established the Abbassa Strain of Nile tilapia (AS) in 2002 based on local broodstock sourced from the Nile River. The AS has been intensively selected for growth and has gone through genetic bottlenecks which have likely shifted levels and composition of genetic diversity within the strain. Consequently, there are questions on the possible genetic impact AS escapees may have on endemic populations of Nile tilapia. However, to date there have been no genetic studies comparing genetic changes in the domesticated AS to local wild populations. This study used 9,827 genome-wide SNPs to investigate population genetic structure and signatures of selection in the AS (generations 9–11) and eight wild Nile tilapia populations from Egypt. SNP analyses identified two major genetic clusters (captive and wild populations), with wild populations showing evidence of isolation-by-distance among the Nile Delta and upstream riverine populations. Between genetic clusters, approximately 6.9% of SNPs were identified as outliers with outliers identified on all 22 O. niloticus chromosomes. A lack of localized outlier clustering on the genome suggests that no genes of major effect were presently detected. The AS has retained high levels of genetic diversity (Ho_All = 0.21 ± 0.01; He_All = 0.23 ± 0.01) when compared to wild populations (Ho_All = 0.18 ± 0.01; He_All = 0.17 ± 0.01) after 11 years of domestication and selective breeding. Additionally, 565 SNPs were unique within the AS line. While these private SNPs may be due to domestication signals or founder effects, it is suspected that introgression with blue tilapia (Oreochromis aureus) has occurred. This study highlights the importance of understanding the effects of domestication in addition to wild population structure to inform future management and dissemination decisions. Furthermore, by conducting a baseline genetic study of wild populations prior to the dissemination of a domestic line, the effects of aquaculture on these populations can be monitored over time

The El Niño Southern Oscillation drives multidirectional inter-reef larval connectivity in the Great Barrier Reef

The El Niño Southern Oscillation (ENSO) is the strongest source of interannual global climate variability, and extreme ENSO events are projected to increase in frequency under climate change. Interannual variability in the Coral Sea circulation has been associated with ENSO, although uncertainty remains regarding ENSO's influence on hydrodynamics and larval dispersal in the adjacent Great Barrier Reef (GBR). We investigated larval connectivity during ENSO events from 2010 to 2017 throughout the GBR, based on biophysical modelling of a widespread predatory reef fish, Lutjanus carponotatus. Our results indicate a well-connected system over the study period with high interannual variability in inter-reef connectivity associated with ENSO. Larval connectivity patterns were highly correlated to variations in the Southern Oscillation Index (SOI). During El Niño conditions and periods of weak SOI, larval dispersal patterns were predominantly poleward in the central and southern regions, reversing to a predominant equatorward flow during very strong SOI and extreme La Niña conditions. These ENSO-linked connectivity patterns were associated with positive connectivity anomalies among reefs. Our findings identify ENSO as an important source of variation in larval dispersal and connectivity patterns in the GBR, which can influence the stability of population dynamics and patterns of biodiversity in the region

microDecon: a highly accurate read‐subtraction tool for the post‐sequencing removal of contamination in metabarcoding studies

Contamination is a ubiquitous problem in microbiome research and can skew results, especially when small amounts of target DNA are available. Nevertheless, no clear solution has emerged for removing microbial contamination. To address this problem, we developed the R package microDecon (https://github.com/donaldtmcknight/microDecon), which uses the proportions of contaminant operational taxonomic units (OTUs) or amplicon sequence variants (ASVs) in blank samples to systematically identify and remove contaminant reads from metabarcoding data sets. We rigorously tested microDecon using a series of computer simulations and a sequencing experiment. We also compared it to the common practice of simply removing all contaminant OTUs/ASVs and other methods for removing contamination. Both the computer simulations and our sequencing data confirmed the utility of microDecon. In our largest simulation (100,000 samples), using microDecon improved the results in 98.1% of samples. Additionally, in the sequencing data and in simulations involving groups, it enabled accurate clustering of groups as well as the detection of previously obscured patterns. It also produced more accurate results than the existing methods for identifying and removing contamination. These results demonstrate that microDecon effectively removes contamination across a broad range of situations. It should, therefore, be widely applicable to microbiome studies, as well as to metabarcoding studies in general

SNP data reveals the complex and diverse evolutionary history of the blue-ringed octopus genus (Octopodidae: Hapalochlaena) in the Asia-Pacific

The blue-ringed octopus species complex (Hapalochlaena spp.), known to occur from Southern Australia to Japan, currently contains four formally described species (Hapalochlaena maculosa, Hapalochlaena fasciata, Hapalochlaena lunulata and Hapalochlaena nierstraszi). These species are distinguished based on morphological characters (iridescent blue rings and/or lines) along with reproductive strategies. However, the observation of greater morphological diversity than previously captured by the current taxonomic framework indicates that a revision is required. To examine species boundaries within the genus we used mitochondrial (12S rRNA, 16S rRNA, cytochrome c oxidase subunit 1 [COI], cytochrome c oxidase subunit 3 [COIII] and cytochrome b [Cytb]) and genome-wide SNP data (DaRT seq) from specimens collected across its geographic range including variations in depth from 3 m to >100 m. This investigation indicates substantially greater species diversity present within the genus Hapalochlaena than is currently described. We identified 10,346 SNPs across all locations, which when analysed support a minimum of 11 distinct clades. Bayesian phylogenetic analysis of the mitochondrial COI gene on a more limited sample set dates the diversification of the genus to ∼30 mya and corroborates eight of the lineages indicated by the SNP analyses. Furthermore, we demonstrate that the diagnostic lined patterning of H. fasciata found in North Pacific waters and NSW, Australia is polyphyletic and therefore likely the result of convergent evolution. Several “deep water” (>100 m) lineages were also identified in this study with genetic convergence likely to be driven by external selective pressures. Examination of morphological traits, currently being undertaken in a parallel morphological study, is required to describe additional species within the complex