Search

Abstract What little we know about how microbiomes change over the course of host dispersal has been gleaned from simulations or snapshot sampling of microbiomes of hosts undertaking regular, cyclical migrations. These studies suggest that major changes in both microbiome richness and turnover occur in response to long-distance movements, but we do not yet know how rare or sporadic dispersal events for non-migratory organisms might affect the microbiomes of their hosts. Here we directly examine the microbiomes of rafting seaweed, leveraging host genomic analyses, amplicon sequencing, and oceanographic modelling to study the impacts of ecological dispersal of hosts on their microbiomes. We find that once dislodged from coastal shores and adrift, kelp-associated microbial communities change profoundly—the core microbes found on attached kelp give way to a few abundant taxa and many rare taxa. Changes in microbial species richness and composition are strongly linked to variability in sea surface temperature rather than length of time spent rafting. These changes are associated with increased contributions of neutral processes shaping community assembly. These findings highlight the role of environmental predictability in triggering major community successional changes and challenge the importance of host selection in determining the microbiome.

Abstract Background Invasive management strategies range from preventing new invasive species incursions to eliminating established populations, with all requiring effective monitoring to guide action. The use of DNA sampled from the environment (eDNA) is one such tool that provides the ability to surveille and monitor target invasive species through passive sampling. Technology being developed to eliminate invasive species includes genetic biocontrol in the form of gene drive. This approach would drive a trait through a population and could be used to eliminate or modify a target population. Once a gene drive organism is released into a population then monitoring changes in density of the target species and the spread of the drive in the population would be critical. Results In this paper, we use invasive Mus musculus as a model for development of an eDNA assay that detects wild-type M. musculus and gene drive M. musculus. We demonstrate successful development of an assay where environmental samples could be used to detect wild-type invasive M. musculus and the relative density of wild-type to gene drive M. musculus. Conclusions The development of a method that detects both wild-type M. musculus and a gene drive M. musculus (t CRISPR) from environmental samples expands the utility of environmental DNA. This method provides a tool that can immediately be deployed for invasive wild M. musculus management across the world. This is a proof-of-concept that a genetic biocontrol construct could be monitored using environmental samples.

Squamate reptiles are amongst the most successful terrestrial vertebrate lineages, with over 10,000 species across a broad range of ecosystems. Despite their success, squamates are also amongst the least studied lineages immunologically. Recently, a universal lack of γδ T cells in squamates due to deletions of the genes encoding the T cell receptor (TCR) γ and δ chains was discovered. Here, we begin to address how the loss of γδ T cells may have impacted the evolution of the squamate immune system. Using the skink Tiliqua rugosa, we found that squamates have not significantly increased the complexity of conventional T cell receptor beta (TCRβ or TRB) chain V regions compared to that of the nearest living squamate relative, the tuatara, Sphenodon punctatus or other amniotes. Our analyses include a putative new TCR locus. This novel locus contains V, D, and J gene segments that undergo V(D)J recombination, albeit with a limited number of gene segments in most squamate species. Based on conserved residues, the predicted protein chain would be expected to form a heterodimer with TCRα. This new TCR locus appears to be derived from an ancient duplication of the TRB locus and is homologous to the recently described T cell receptor epsilon (TRE). TRE is absent from the genomes of the tuatara and all Archosaurs examined and appears squamate specific.

ABSTRACT Environmental DNA (eDNA) offers an efficient approach to biodiversity monitoring and biosecurity surveillance in coastal ecosystems. Understanding eDNA dynamics in tidal areas is critical for guiding sampling strategy design and interpreting molecular detection results. In a field study conducted at Ōpua, Bay of Islands, Aotearoa, New Zealand, an international marina with a dynamic coastal environment, we investigated eDNA detection patterns across a local spatial and temporal gradient. Hourly samples were collected from six sampling stations around the marina, from artificial shore structures and a channel, over a 12‐h tidal cycle. Detected eukaryotic communities were assessed using eDNA metabarcoding analysis of the small‐subunit ribosomal RNA (18S rRNA) and mitochondrial Cytochrome C Oxidase subunit I (COI) genes. This core community was screened for marine nonindigenous species (NIS) using an online metabarcoding data screening app (ExPAT). Community composition varied significantly based on sampling location and time, with most of the core community captured around mid‐tide. The NIS community exhibited inconsistent detection for some species (e.g., the amphipod Monocorophium acherusicum) across space and time, whereas other NIS (e.g., the mussel Arcuatula senhousia) were consistently detected and behaved similar to the core community. Overall, species‐specific eDNA detection patterns may be linked to the biology and movement of the organism. The results indicated slightly higher NIS detection from artificial structures at the shore stations and within 2 h of low tide. For marine biosecurity applications, users may want to prioritize sampling near‐shore during low tide, while samples collected during mid‐tide may provide more comprehensive insights into the core community for biomonitoring. Further field studies across multiple tidal cycles and marinas may guide the standardization of molecular biomonitoring and biosurveillance sampling strategies and support their broader integration into marine biosecurity programs.

Abstract Combining genome assembly with population and functional genomics can provide valuable insights to development and evolution, as well as tools for species management. Here, we present a chromosome-level genome assembly of the common brushtail possum (Trichosurus vulpecula), a model marsupial threatened in parts of their native range in Australia, but also a major introduced pest in New Zealand. Functional genomics reveals post-natal activation of chemosensory and metabolic genes, reflecting unique adaptations to altricial birth and delayed weaning, a hallmark of marsupial development. Nuclear and mitochondrial analyses trace New Zealand possums to distinct Australian subspecies, which have subsequently hybridised. This admixture allowed phasing of parental alleles genome-wide, ultimately revealing at least four genes with imprinted, parent-specific expression not yet detected in other species (MLH1, EPM2AIP1, UBP1 and GPX7). We find that reprogramming of possum germline imprints, and the wider epigenome, is similar to eutherian mammals except onset occurs after birth. Together, this work is useful for genetic-based control and conservation of possums, and contributes to understanding of the evolution of novel mammalian epigenetic traits.

Abstract Analysis of environmental DNA (eDNA) has gained widespread usage for taxonomically based biodiversity assessment. While interest in applying noninvasive eDNA monitoring for population genetic assessments has grown, its usage in this sphere remains limited. One barrier to uptake is that the effectiveness of eDNA detection below the species level remains to be determined for multiple species and environments. Here, we test the utility of this emergent technology to obtain within‐species haplotypic variation of New Zealand (NZ) blackfoot pāua (Haliotis iris). We compare mitochondrial haplotype diversity recovered from marine eDNA samples against traditional tissue samples of blackfoot pāua collected at the same NZ coastal site. Targeting the ATP8‐ATP6 region, we recovered four mitochondrial haplotypes from eDNA versus six haplotypes from tissue samples. Three common haplotypes were recovered with both eDNA and tissue samples, while only one out of three rare haplotypes – represented in tissue samples by one individual each – was recovered with our eDNA methods. We demonstrate that eDNA monitoring is an effective tool for recovering common genetic diversity from pāua, although rare (

Abstract Environmental DNA and RNA (eDNA/eRNA) can serve as molecular tools for biodiversity monitoring and biosecurity surveillance. However, uncertainties still exist regarding the persistence and dynamics of marine nucleic acids in the environment and the effects of post‐sampling storage on species detectability. To bridge these gaps, an experiment was conducted in an Auckland marina, a known New Zealand entry point for marine non‐indigenous species (NIS). We targeted a prominent invader, the Mediterranean fanworm Sabella spallanzanii, for eDNA/eRNA‐based detection. Permeable dialysis bags filled with seawater collected near an S. spallanzanii colony were deployed in the marina to simulate environmental conditions, with a subset of bags stored on ice to mimic field storage conditions. Sabella spallanzanii eDNA/eRNA signal was quantified using droplet digital PCR on samples collected over 24 h of dialysis bag deployment. Results challenged traditional first‐order decay models, showing inconsistent eDNA/eRNA signal patterns and no significant concentration changes between 0 and 24 h. Consequently, total eDNA fragmentation was assessed using the Agilent 2100 Bioanalyzer® electrophoresis system, which revealed a rise in the number of total eDNA fragments and an unexpected increase in the median fragment size of total eDNA under field conditions over time, likely originating from the ambient microbiome. Additional analysis using long‐read sequencing (Oxford Nanopore Technologies, UK) revealed an increase in microbial eDNA reads within the in‐field samples, suggesting potential microbial growth within the dialysis bags. In contrast, the ice‐stored samples exhibited no significant changes in the number of reads assigned to microbial taxa, implying limited microbial growth in cold storage. These findings provide insights into total eDNA dynamics and its potential impact on targeted eDNA concentration and detection. Further comprehensive research on eDNA/eRNA dynamics, particularly focused on eRNA, is essential, as this understanding is crucial for refining survey interpretation and sampling design for effective environmental management.

Molecular biomonitoring programs increasingly use environmental DNA (eDNA) for detecting targeted species such as marine non-indigenous species (NIS) or endangered species. However, the current molecular detection workflow is cumbersome and time-demanding, and thereby can hinder management efforts and restrict the “opportunity window” for rapid management responses. Here, we describe a direct droplet digital PCR (direct-ddPCR) approach to detect species-specific free-floating extra-cellular eDNA (free-eDNA) signals, i.e., detection of species-specific eDNA without the need for filtration or DNA extraction, with seawater samples. This first proof-of-concept aquarium study was conducted with three distinct marine species: the Mediterranean fanworm Sabella spallanzanii, the ascidian clubbed tunicate Styela clava, and the brown bryozoan Bugula neritina to evaluate the detectability of free-eDNA in seawater. The detectability of targeted free-eDNA was assessed by directly analysing aquarium marine water samples using an optimized species-specific ddPCR assay. The results demonstrated the consistent detection of S. spallanzanii and B. neritina free-eDNA when these organisms were present in high abundance. Once organisms were removed, the free-eDNA signal exponentially declined, noting that free-eDNA persisted between 24–72 h. Results indicate that organism biomass, specimen characteristics (e.g., stress and viability), and species-specific biological differences may influence free-eDNA detectability. This study represents the first step in assessing the feasibility of direct-ddPCR technology for the detection of marine species. Our results provide information that could aid in the development of new technology, such as a field development of ddPCR systems, which could allow for automated continuous monitoring of targeted marine species, enabling point-of-need detection and rapid management responses.

Abstract Environmental DNA (eDNA) metabarcoding has shown great promise as an effective, non‐invasive monitoring method for marine biomes. However, long filtration times and the need for state‐of‐the‐art laboratories are restricting sample replication and in situ species detections. Methodological innovations, such as passive filtration and self‐contained DNA extraction protocols, have the potential to alleviate these issues. We explored the implementation of passive sampling and a self‐contained DNA extraction protocol by comparing fish diversity obtained from active filtration (1 L; 0.45 μm cellulose nitrate [CN] filters) to five passive substrates, including 0.45 μm CN filters, 5 μm nylon filters, 0.45 μm positively charged nylon filters, artificial sponges, and fishing net. Fish diversity was then compared between the PDQeX Nucleic Acid Extractor and the conventional Qiagen DNeasy Blood & Tissue protocol. Experiments were conducted in both a controlled mesocosm and in situ at the Portobello Marine Laboratory, New Zealand. No significant differences in fish diversity were observed among active filtration and more porous passive materials (artificial sponges and fishing net) for both the mesocosm and harbor waters. For the in situ comparison, all passive filter membranes detected a significantly lower number of fish species, resulting from partial sample drop‐out. While no significant differences in fish eDNA signal diversity were observed between either DNA extraction methods in the mesocosm, the PDQeX system was less effective at detecting fish for the in situ comparison. Our results demonstrate that a passive sampling approach using porous substrates can be effectively implemented to capture eDNA from seawater, eliminating vacuum filtration processing. The large variation in efficiency observed among the five substrate types, however, warrants further optimization of the passive sampling approach for routine eDNA applications. The PDQeX system can extract high‐abundance DNA in a mesocosm with further optimization to detect low‐abundance eDNA from the marine environment.

Abstract Background In a time of rapid environmental change, understanding how the challenges experienced by one generation can influence the fitness of future generations is critically needed. Using tolerance assays and transcriptomic and methylome approaches, we use zebrafish as a model to investigate cross-generational acclimation to hypoxia. Results We show that short-term paternal exposure to hypoxia endows offspring with greater tolerance to acute hypoxia. We detected two hemoglobin genes that are significantly upregulated by more than 6-fold in the offspring of hypoxia exposed males. Moreover, the offspring which maintained equilibrium the longest showed greatest upregulation in hemoglobin expression. We did not detect differential methylation at any of the differentially expressed genes, suggesting that other epigenetic mechanisms are responsible for alterations in gene expression. Conclusions Overall, our findings suggest that an epigenetic memory of past hypoxia exposure is maintained and that this environmentally induced information is transferred to subsequent generations, pre-acclimating progeny to cope with hypoxic conditions.

Active environmental DNA (eDNA) surveillance through species-specific amplification has shown increased sensitivity in the detection of non-indigenous species (NIS) compared to traditional approaches. When many NIS are of interest, however, active surveillance decreases in cost- and time-efficiency. Passive surveillance through eDNA metabarcoding takes advantage of the complex DNA signal in environmental samples and facilitates the simultaneous detection of multiple species. While passive eDNA surveillance has previously detected NIS, comparative studies are essential to determine the ability of eDNA metabarcoding to accurately describe the range of invasion for multiple NIS versus alternative approaches. Here, we surveyed twelve sites, covering nine rivers across Belarus for NIS with three different techniques, i.e. an ichthyological, hydrobiological and eDNA survey, whereby DNA was extracted from 500 ml surface water samples and amplified with two 16S rDNA primer assays targeting the fish and macroinvertebrate biodiversity. Nine non-indigenous fish and ten non-indigenous benthic macroinvertebrates were detected by traditional surveys, while seven NIS eDNA signals were picked up, including four fish, one aquatic and two benthic macroinvertebrates. Passive eDNA surveillance extended the range of invasion further north for two invasive fish and identified a new NIS for Belarus, the freshwater jellyfish Craspedacusta sowerbii. False-negative detections for the eDNA survey might be attributed to: (i) preferential amplification of aquatic over benthic macroinvertebrates from surface water samples and (ii) an incomplete reference database. The evidence provided in this study recommends the implementation of both molecular-based and traditional approaches to maximise the probability of early detection of non-native organisms.

Summary: Programmed DNA loss is a gene silencing mechanism that is employed by several vertebrate and nonvertebrate lineages, including all living jawless vertebrates and songbirds. Reconstructing the evolution of somatically eliminated (germline-specific) sequences in these species has proven challenging due to a high content of repeats and gene duplications in eliminated sequences and a corresponding lack of highly accurate and contiguous assemblies for these regions. Here, we present an improved assembly of the sea lamprey (Petromyzon marinus) genome that was generated using recently standardized methods that increase the contiguity and accuracy of vertebrate genome assemblies. This assembly resolves highly contiguous, somatically retained chromosomes and at least one germline-specific chromosome, permitting new analyses that reconstruct the timing, mode, and repercussions of recruitment of genes to the germline-specific fraction. These analyses reveal major roles of interchromosomal segmental duplication, intrachromosomal duplication, and positive selection for germline functions in the long-term evolution of germline-specific chromosomes.

Robert Macey et al. use deep sequencing to reveal two co-existing, deeply divergent, fully functional mitochondrial genomes in a single Tuatara from Lady Alice Island, New Zealand. This finding will stimulate further investigation, having profound implications for our understanding of mitochondrial evolution and function.

Summary: The kākāpō is a flightless parrot endemic to New Zealand. Once common in the archipelago, only 201 individuals remain today, most of them descending from an isolated island population. We report the first genome-wide analyses of the species, including a high-quality genome assembly for kākāpō, one of the first chromosome-level reference genomes sequenced by the Vertebrate Genomes Project (VGP). We also sequenced and analyzed 35 modern genomes from the sole surviving island population and 14 genomes from the extinct mainland population. While theory suggests that such a small population is likely to have accumulated deleterious mutations through genetic drift, our analyses on the impact of the long-term small population size in kākāpō indicate that present-day island kākāpō have a reduced number of harmful mutations compared to mainland individuals. We hypothesize that this reduced mutational load is due to the island population having been subjected to a combination of genetic drift and purging of deleterious mutations, through increased inbreeding and purifying selection, since its isolation from the mainland ∼10,000 years ago. Our results provide evidence that small populations can survive even when isolated for hundreds of generations. This work provides key insights into kākāpō breeding and recovery and more generally into the application of genetic tools in conservation efforts for endangered species.

Abstract Background The utility of environmental DNA (eDNA) metabarcoding surveys to accurately detect species depends on the degree of DNA dispersal. Multiple marine studies have observed only minimal eDNA transport by horizontal water movement across small spatial scales, leading to the conclusion that spatially specific eDNA signals accurately resemble in‐field species assemblages along a horizontal axis. Marine communities, however, are also structured vertically according to depth. In marine environments displaying permanent water stratification, vertical zonation patterns may be more apparent and present on smaller spatial scales (i.e., meters) than horizontal community structuring. The scale at which eDNA signals differ along a vertical transect and the accuracy of eDNA metabarcoding in revealing naturally stratified communities have yet to be assessed. Methods and results In this study, we determined the ability of eDNA metabarcoding surveys to distinguish vertically localized community assemblages. To test this, we sampled three vertical transects along a steep rock wall at three depths (0 m, 4 m, 15 m), covering two distinct communities that were separated by near‐permanent water column stratification in the form of a strong halocline at ~3 m. Using three metabarcoding assays, our eDNA metabarcoding survey detected 54 taxa, across 46 families and 7 phyla, including 19 fish, 15 crustacean, and 8 echinoderm species. Ordination and cluster analyses show distinct eDNA signals across the halocline for all three replicate transects, suggesting that vertical dispersal of eDNA between communities was limited. Furthermore, eDNA signals of individual taxa were only retrieved within their observed vertical distribution, providing biological validation for the obtained results. Our results demonstrate, for the first time, the need to take into consideration oceanographic (e.g. water column stratification) and biological processes (e.g. vertical community structuring) when designing sampling strategies for marine eDNA metabarcoding surveys.

Germline cells transfer genetic information to offspring, and in zebrafish, drive sex determination. Here the authors report that, unlike mammals, the germline of zebrafish does not undergo genome-wide DNA methylation erasure, while amplifying and demethylating sex-linked rDNA during feminisation.

Abstract DNA extraction from environmental samples (environmental DNA; eDNA) for metabarcoding‐based biodiversity studies is gaining popularity as a noninvasive, time‐efficient, and cost‐effective monitoring tool. The potential benefits are promising for marine conservation, as the marine biome is frequently under‐surveyed due to its inaccessibility and the consequent high costs involved. With increasing numbers of eDNA‐related publications have come a wide array of capture and extraction methods. Without visual species confirmation, inconsistent use of laboratory protocols hinders comparability between studies because the efficiency of target DNA isolation may vary. We determined an optimal protocol (capture and extraction) for marine eDNA research based on total DNA yield measurements by comparing commonly employed methods of seawater filtering and DNA isolation. We compared metabarcoding results of both targeted (small taxonomic group with species‐level assignment) and universal (broad taxonomic group with genus/family‐level assignment) approaches obtained from replicates treated with the optimal and a low‐performance capture and extraction protocol to determine the impact of protocol choice and DNA yield on biodiversity detection. Filtration through cellulose‐nitrate membranes and extraction with Qiagen's DNeasy Blood & Tissue Kit outperformed other combinations of capture and extraction methods, showing a ninefold improvement in DNA yield over the poorest performing methods. Use of optimized protocols resulted in a significant increase in OTU and species richness for targeted metabarcoding assays. However, changing protocols made little difference to the OTU and taxon richness obtained using universal metabarcoding assays. Our results demonstrate an increased risk of false‐negative species detection for targeted eDNA approaches when protocols with poor DNA isolation efficacy are employed. Appropriate optimization is therefore essential for eDNA monitoring to remain a powerful, efficient, and relatively cheap method for biodiversity assessments. For seawater, we advocate filtration through cellulose‐nitrate membranes and extraction with Qiagen's DNeasy Blood & Tissue Kit or phenol‐chloroform‐isoamyl for successful implementation of eDNA multi‐marker metabarcoding surveys.

The Illumina Infinium® MethylationEPIC BeadChip system (hereafter EPIC array) is considered to be the current gold standard detection method for assessing DNA methylation at the genome-wide level. EPIC arrays are often used for hypothesis generation or pilot studies, the natural conclusion to which is to validate methylation candidates and expand these in a larger cohort, in a targeted manner. As such, an accurate smaller-scale, targeted technique, that generates data at the individual CpG level that is equivalent to the EPIC array, is needed. Here, we tested an alternative DNA methylation detection technique, known as bisulfite-based amplicon sequencing (BSAS), to determine its ability to validate CpG sites detected in EPIC array studies. BSAS was able to detect differential DNA methylation at CpG sites to a degree which correlates highly with the EPIC array system at some loci. However, BSAS correlated less well with EPIC array data in some instances, and most notably, when the magnitude of change via EPIC array was greater than 5%. Therefore, our data suggests that BSAS can be used to validate EPIC array data, but each locus must be compared on an individual basis, before being taken forward into large scale screening. Further, BSAS does offer advantages compared to the probe-based EPIC array; BSAS amplifies a region of the genome (∼500 bp) around a CpG of interest, allowing analyses of other CpGs in the region that may not be present on the EPIC array, aiding discovery of novel CpG sites and differentially methylated regions of interest. We conclude that BSAS offers a valid investigative tool for specific regions of the genome that are currently not contained on the array system.

Many teleost fishes undergo natural sex change, and elucidating the physiological and molecular controls of this process offers unique opportunities not only to develop methods of controlling sex in aquaculture settings, but to better understand vertebrate sexual development more broadly. Induction of sex change in some sequentially hermaphroditic or gonochoristic fish can be achieved in vivo through social manipulation, inhibition of aromatase activity, or steroid treatment. However, the induction of sex change in vitro has been largely unexplored. In this study, we established an in vitro culture system for ovarian explants in serum-free medium for a model sequential hermaphrodite, the New Zealand spotty wrasse (Notolabrus celidotus). This culture technique enabled evaluating the effect of various treatments with 17β-estradiol (E2), 11-ketotestosterone (11KT) or cortisol (CORT) on spotty wrasse ovarian architecture for 21 days. A quantitative approach to measuring the degree of ovarian atresia within histological images was also developed, using pixel-based machine learning software. Ovarian atresia likely due to culture was observed across all treatments including no-hormone controls, but was minimised with treatment of at least 10 ng/mL E2. Neither 11KT nor CORT administration induced proliferation of spermatogonia (i.e., sex change) in the cultured ovaries indicating culture beyond 21 days may be needed to induce sex change in vitro. The in vitro gonadal culture and analysis systems established here enable future studies investigating the paracrine role of sex steroids, glucocorticoids and a variety of other factors during gonadal sex change in fish.

Mollusc shells are an abundant resource that have been long used to predict the structures of ancient ecological communities, examine evolutionary processes, reconstruct paleoenvironmental conditions, track and predict responses to climatic change, and explore the movement of hominids across the globe. Despite the ubiquity of mollusc shell in many environments, it remains relatively unexplored as a substrate for molecular genetic analysis. Here we undertook a series of experiments using the New Zealand endemic greenshell mussel, Perna canaliculus, to explore the utility of fresh, aged, beach-cast and cooked mollusc shell for molecular genetic analyses. We find that reasonable quantities of DNA (0.002–21.48 ng/mg shell) can be derived from aged, beach-cast and cooked mussel shell and that this can routinely provide enough material to undertake PCR analyses of mitochondrial and nuclear gene fragments. Mitochondrial PCR amplification had an average success rate of 96.5% from shell tissue extracted thirteen months after the animal’s death. A success rate of 93.75% was obtained for cooked shells. Amplification of nuclear DNA (chitin synthase gene) was less successful (80% success from fresh shells, decreasing to 10% with time, and 75% from cooked shells). Our results demonstrate the promise of mollusc shell as a substrate for genetic analyses targeting both mitochondrial and nuclear genes.

Abstract Background The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events. Results Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes. Conclusions The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda.

Abstract Background Sex hormones play important roles in teleost ovarian and testicular development. In zebrafish, ovarian differentiation appears to be dictated by an oocyte-derived signal via Cyp19a1a aromatase-mediated estrogen production. Androgens and aromatase inhibitors can induce female-to-male sex reversal, however, the mechanisms underlying gonadal masculinisation are poorly understood. We used histological analyses together with RNA sequencing to characterise zebrafish gonadal transcriptomes and investigate the effects of 17α-methyltestosterone on gonadal differentiation. Results At a morphological level, 17α-methyltestosterone (MT) masculinised gonads and accelerated spermatogenesis, and these changes were paralleled in masculinisation and de-feminisation of gonadal transcriptomes. MT treatment upregulated expression of genes involved in male sex determination and differentiation (amh, dmrt1, gsdf and wt1a) and those involved in 11-oxygenated androgen production (cyp11c1 and hsd11b2). It also repressed expression of ovarian development and folliculogenesis genes (bmp15, gdf9, figla, zp2.1 and zp3b). Furthermore, MT treatment altered epigenetic modification of histones in zebrafish gonads. Contrary to expectations, higher levels of cyp19a1a or foxl2 expression in control ovaries compared to MT-treated testes and control testes were not statistically significant during early gonad development (40 dpf). Conclusion Our study suggests that both androgen production and aromatase inhibition are important for androgen-induced gonadal masculinisation and natural testicular differentiation in zebrafish.

Fishes exhibit remarkably diverse, and plastic, patterns of sexual development, most striking of which is sequential hermaphroditism, where individuals readily reverse sex in adulthood. How this stunning example of phenotypic plasticity is controlled at a genetic level remains poorly understood. Several genes have been implicated in regulating sex change, yet the degree to which a conserved genetic machinery orchestrates this process has not yet been addressed. Using captive and in-the-field social manipulations to initiate sex change, combined with a comparative qPCR approach, we compared expression patterns of four candidate regulatory genes among three species of wrasses (Labridae)—a large and diverse teleost family where female-to-male sex change is pervasive, socially-cued, and likely ancestral. Expression in brain and gonadal tissues were compared among the iconic tropical bluehead wrasse (Thalassoma bifasciatum) and the temperate spotty (Notolabrus celidotus) and kyusen (Parajulus poecilepterus) wrasses. In all three species, gonadal sex change was preceded by downregulation of cyp19a1a (encoding gonadal aromatase that converts androgens to oestrogens) and accompanied by upregulation of amh (encoding anti-müllerian hormone that primarily regulates male germ cell development), and these genes may act concurrently to orchestrate ovary-testis transformation. In the brain, our data argue against a role for brain aromatase (cyp19a1b) in initiating behavioural sex change, as its expression trailed behavioural changes. However, we find that isotocin (it, that regulates teleost socio-sexual behaviours) expression correlated with dominant male-specific behaviours in the bluehead wrasse, suggesting it upregulation mediates the rapid behavioural sex change characteristic of blueheads and other tropical wrasses. However, it expression was not sex-biased in temperate spotty and kyusen wrasses, where sex change is more protracted and social groups may be less tightly-structured. Together, these findings suggest that while key components of the molecular machinery controlling gonadal sex change are phylogenetically conserved among wrasses, neural pathways governing behavioural sex change may be more variable.

We report the complete mitogenome of Sitona obsoletus, an agricultural pest in New Zealand and some European countries. Like other Sitona mitogenomes, the 6 tRNA gene box is ordered RNSAEF, supporting the hypothesis that this signature is common to, and potentially diagnostic, of this genus. The Trojan Female Technique (TFT) is a genetic pest control strategy that exploits mitochondrial DNA alleles that affect male, but not female fertility and fitness. The complete mitogenome is an essential first step in exploring the utility of TFT for the control of S. obsoletus.

Microsatellites are the genetic markers of choice for many population genetic studies, but must be isolated de novo using recombinant approaches where prior genetic data are lacking. Here we utilized high-throughput genomic sequencing technology to produce millions of base pairs of short fragment reads, which were screened with bioinformatics toolsets to identify primers that amplify polymorphic microsatellite loci. Using this approach we isolated 13 polymorphic microsatellites for the blue duck (Hymenolaimus malacorhynchos), a species for which limited genetic data were available. Our genomic approach eliminates recombinant genetic steps, significantly reducing the time and cost requirements of marker development compared with traditional approaches. While this application of genomic sequencing may seem obvious to many, this study is, to the best of our knowledge, the first attempt to describe the use of genomic sequencing for the development of microsatellite markers in a non-model organism or indeed any organism.

TaqMan-nuclease assays are widely used for the qualitative detection of single nucleotide polymorphisms (SNPs) and the determination of biallelic states in pooled or heterozygous DNA samples. These assays are highly specific, reproducible, and suitable for high-throughput approaches. A crucial limitation of this method, and others, is the detection of minor allele frequencies with detection limits of generally 3% to 9% for minor allele contributions. Here we describe the combination of customized TaqMan-nuclease assay and allele-specific restriction to increase the sensitivity of this method, allowing the qualitative detection of allele contributions as low as 0.05%.

Microsatellites are currently one of the most commonly used genetic markers. The application of bioinformatic tools has become common practice in the study of these short tandem repeats (STR). However, in silico studies can suffer from study bias. Using a meta-analysis on microsatellite distribution in yeast we show that estimates of numbers of repeats reported by different studies can differ in the order of several magnitudes, even within a single genome. These differences arise because varying definitions of microsatellites, spanning repeat size, array length and array composition, are used in different search paradigms, with minimum array length being the main influencing factor. Structural differences in the implemented search algorithm additionally contribute to variation in the number of repeats detected. We suggest that for future studies a consistent approach to STR searches is adopted in order to improve the power of intra- and interspecific comparisons