A collaborative backbone resource for comparative studies of subterranean evolution: The World Asellidae database.

Transition to novel environments, such as groundwater colonization by surface organisms, provides an excellent research ground to study phenotypic evolution. However, interspecific comparative studies on evolution to groundwater life are few because of the challenge in assembling large ecological and molecular resources for species-rich taxa comprised of surface and subterranean species. Here, we make available to the scientific community an operational set of working tools and resources for the Asellidae, a family of freshwater isopods containing hundreds of surface and subterranean species. First, we release the World Asellidae database (WAD) and its web application, a sustainable and FAIR solution to producing and sharing data and biological material. WAD provides access to thousands of species occurrences, specimens, DNA extracts and DNA sequences with rich metadata ensuring full scientific traceability. Second, we perform a large-scale dated phylogenetic reconstruction of Asellidae to support phylogenetic comparative analyses. Of 424 terminal branches, we identify 34 pairs of surface and subterranean species representing independent replicates of the transition from surface water to groundwater. Third, we exemplify the usefulness of WAD for documenting phenotypic shifts associated with colonization of subterranean habitats. We provide the first phylogenetically controlled evidence that body size of males decreases relative to that of females upon groundwater colonization, suggesting competition for rare receptive females selects for smaller, more agile males in groundwater. By making these tools and resources widely accessible, we open up new opportunities for exploring how phenotypic traits evolve in response to changes in selective pressures and trade-offs during groundwater colonization.

Cryptic diversity and phylogeographic patterns of Deto echinata (Isopoda: Detonidae) in southern Africa.

Recent phylogeographic studies of poorly-dispersing coastal invertebrates in highly biodiverse regions have led to the discovery of high levels of cryptic diversity and complex phylogeographic patterns that suggest isolation, geological, and ecological processes have shaped their biodiversity. Studies of southern African coastal invertebrates have uncovered cryptic diversity for various taxa and phylogeographic patterns that, although sharing some similarities across taxa, do differ. These findings underscore the need for additional studies to better understand the biodiversity levels, distributional patterns, and processes responsible for producing coastal biodiversity in that region. The coastal isopod Deto echinata is of particular interest, as its complex taxonomic history, poor dispersal capabilities, and broad geographic distribution suggest the potential for cryptic diversity. We use mitochondrial and nuclear sequences to characterize D. echinata individuals from localities ranging from northern Namibia to Glentana, about 2,500 km along the coastline on the south coast of South Africa. These are used to assess whether D. echinata harbors cryptic genetic diversity and whether phylogeographic distributional patterns correlate with those previously documented for other coastal isopods in the region. Analysis of mitochondrial and nuclear sequences revealed two deeply-divergent lineages that exhibit a distributional break in the Cape Peninsula region. These findings suggest D. echinata is a cryptic species complex in need of taxonomic revision and highlight the need for further taxonomic and phylogeographic studies of similarly poorly-dispersing coastal invertebrates in southern Africa.

The circadian clock gene bmal1 is necessary for co-ordinated circatidal rhythms in the marine isopod Eurydice pulchra (Leach).

Circadian clocks in terrestrial animals are encoded by molecular feedback loops involving the negative regulators PERIOD, TIMELESS or CRYPTOCHROME2 and positive transcription factors CLOCK and BMAL1/CYCLE. The molecular basis of circatidal (~12.4 hour) or other lunar-mediated cycles (~15 day, ~29 day), widely expressed in coastal organisms, is unknown. Disrupting circadian clockworks does not appear to affect lunar-based rhythms in several organisms that inhabit the shoreline suggesting a molecular independence of the two cycles. Nevertheless, pharmacological inhibition of casein kinase 1 (CK1) that targets PERIOD stability in mammals and flies, affects both circadian and circatidal phenotypes in Eurydice pulchra (Ep), the speckled sea-louse. Here we show that these drug inhibitors of CK1 also affect the phosphorylation of EpCLK and EpBMAL1 and disrupt EpCLK-BMAL1-mediated transcription in Drosophila S2 cells, revealing a potential link between these two positive circadian regulators and circatidal behaviour. We therefore performed dsRNAi knockdown of Epbmal1 as well as the major negative regulator in Eurydice, Epcry2 in animals taken from the wild. Epcry2 and Epbmal1 knockdown disrupted Eurydice's circadian phenotypes of chromatophore dispersion, tim mRNA cycling and the circadian modulation of circatidal swimming, as expected. However, circatidal behaviour was particularly sensitive to Epbmal1 knockdown with consistent effects on the power, amplitude and rhythmicity of the circatidal swimming cycle. Thus, three Eurydice negative circadian regulators, EpCRY2, in addition to EpPER and EpTIM (from a previous study), do not appear to be required for the expression of robust circatidal behaviour, in contrast to the positive regulator EpBMAL1. We suggest a neurogenetic model whereby the positive circadian regulators EpBMAL1-CLK are shared between circadian and circatidal mechanisms in Eurydice but circatidal rhythms require a novel, as yet unknown negative regulator.

Chemosensory proteins as putative semiochemical carriers in the desert isopod Hemilepistus reaumurii.

The order Isopoda contains both aquatic and terrestrial species, among which Hemilepistus reaumurii, which lives in arid environments and is the most adapted to terrestrial life. Olfaction has been deeply investigated in insects while it has received very limited attention in other arthropods, particularly in terrestrial crustaceans. In insects, soluble proteins belonging to two main families, Odorant Binding Proteins (OBPs) and Chemosensory Proteins (CSPs), are contained in the olfactory sensillar lymph and are suggested to act as carriers of hydrophobic semiochemicals to or from membrane-bound olfactory receptors. Other protein families, namely Nieman-Pick type 2 (NPC2) and Lipocalins (LCNs) have been also reported as putative odorant carriers in insects and other arthropod clades. In this study, we have sequenced and analysed the transcriptomes of antennae and of the first pair of legs of H. reaumurii focusing on soluble olfactory proteins. Interestingly, we have found 13 genes encoding CSPs, whose sequences differ from those of the other arthropod clades, including non-isopod crustaceans, for the presence of two additional cysteine residues, besides the four conserved ones. Binding assays on two of these proteins showed strong affinities for fatty acids and long-chain unsaturated esters and aldehydes, putative semiochemicals for this species.

Homeotic transformation in a terrestrial isopod: insights into the appendage identity in crustaceans.

In many crustacean species, an individual possesses both uniramous and biramous appendages that enable us to compare the two types on the same genetic background. Therefore, among the diverse morphologies of arthropod appendages, crustacean biramous appendages provide interesting subjects for studying the developmental mechanisms underlying appendage modifications. In this study, we report a malformed specimen of the terrestrial isopod Porcellio scaber, in which one of the pleopods was transformed into a different structure. Morphological observations of exoskeletons and musculatures by confocal scanning laser microscopy revealed that the transformed appendage was three-segmented, with at least the apical two segments having pereopod-like musculoskeletal structures. The apical segment of the transformed appendage lacked muscles, and the following segment had a pair of muscle bundles. These findings together with those of some previous studies of gene expression patterns in this species suggest that this anomaly could be caused by homeotic transformation of a flap-like pleopod into a three-segmented pereopod tip, which may be a homologous structure of the pleopod.

History of the terrestrial isopod genus Ligidium in Japan based on phylogeographic analysis.

BACKGROUND: Phylogeographical approaches explain the genetic diversity of local organisms in the context of their geological and geographic environments. Thus, genetic diversity can be a proxy for geological history. Here we propose a genus of woodland isopod, Ligidium, as a marker of geological history in relation to orogeny and the Quaternary glacial cycle. RESULTS: Mitochondrial analysis of 721 individuals from 97 sites across Japan revealed phylogenetic divergence between the northeastern and southwestern Japan arcs. It also showed repeated population expansions in northeastern Japan in response to Quaternary glacial and interglacial cycles. Genome-wide analysis of 83 selected individuals revealed multiple genetic nuclear clusters. The genomic groupings were consistent with the local geographic distribution, indicating that the Ligidium phylogeny reflects its regional history. CONCLUSION: Ligidium DNA sequence analysis can provide insight into the geological, geographical, and paleoenvironmental history of the studied region.

Redescription and molecular characterization of Mothocya parvostis Bruce, 1986 (Crustacea: Isopoda: Cymothoidae) parasitic on Japanese halfbeak, Hyporhamphus sajori (Temminck & Schlegel, 1846) (Hemiramphidae) with Mothocya sajori Bruce, 1986 placed into synonymy.

Two species of Mothocya have previously been recorded from Hyporhamphus sajori: M. parvostis Bruce, 1986 and M. sajori Bruce, 1986. Mothocya parvostis is re-described based on the ovigerous female type and additional materials collected from the host from in and around the type locality. Morphological re-examination of fresh specimens and the type materials together with genetic data show that the M. sajori and M. parvostis are the same species, differing primarily in size, therefore we have placed Mothocya sajori Bruce, 1986 into a junior synonym of Mothocya parvostis Bruce, 1986. Mothocya parvostis is characterized by the following combinations of characters: 1) body slightly to moderately twisted to one side; 2) pereonite 7 posterior margin moderately to deeply recessed; 3) uropodal rami extending to pleotelson posterior margin; and 4) uropod rami bluntly rounded, exopod 1.5 times as long as peduncle. The differences of four morphological features for M. parvostis and M. sajori was quantified. Furthermore, a total of 635 isopods infesting H. sajori were collected from all over Japan to conduct quantitative morphological and molecular sequence analyses (mitochondrial cytochrome c oxidase subunit I and 16S rRNA). Although the four quantitative features did not overlap between the two species in type specimens, all quantitative morphological values of newly collected specimens in this study did not display a bimodal distribution. In addition, our molecular analyses found only a single clade for our newly collected specimens in neighbor-joining tree.

Combining morphological and mitochondrial DNA data to describe a new species of Austroniscus Vanhöffen, 1914 (Isopoda, Janiroidea, Nannoniscidae) linking abyssal and hadal depths of the Puerto Rico Trench.

Hadal trenches are perceived as a unique deep-sea ecosystem with fundamentally different communities compared to the nearby abyss. So far, however, scarce information exists about how populations are genetically linked within a trench and about mechanisms for species divergence. The present study presents the morphological and molecular-genetic characterization and description of a new nannoniscid species within the genus Austroniscus Vanhöffen, 1914 obtained from abyssal and hadal depths of the Puerto Rico Trench, NW Atlantic. Samples were collected as part of the Vema-TRANSIT expedition onboard RV Sonne in January 2015. Because of the large depth differences between sampling locations (4,552-8,338 m), we expected to find different species within the genus inhabiting abyssal and hadal sites. Initial morphological examination using traditional light microscopy and Confocal Laser Scanning Microscopy was paired with subsequent molecular analysis based on mtDNA (COI and 16S). Contrary to our assumptions, combined morphological and molecular species delimitation analyses (sGMYC, mPTP, ABGD) revealed the presence of only one species spanning the abyssal and hadal seafloor of the Puerto Rico Trench. In addition, comparison with type material could show that this species belongs to a new species, Austroniscus brandtae n. sp., which is described herein. Incongruence between some species delimitation methods suggesting the presence of multiple species is interpreted as strong genetic population structuring within the trench, which is also supported by the analysis of the haplotype networks. The geographic and bathymetric distribution of Austroniscus species is discussed. The species described herein represents the first in the genus Austroniscus from the Atlantic Ocean and the deepest record of the genus to date, and hence significantly expanding previously known limits of its geographic and bathymetric range.

Heterogeneous distribution of sex ratio distorters in natural populations of the isopod Armadillidium vulgare.

In the isopod Armadillidium vulgare, many females produce progenies with female-biased sex ratios, owing to two feminizing sex ratio distorters (SRD): Wolbachia endosymbionts and the f element. We investigated the distribution and population dynamics of these SRD and mitochondrial DNA variation in 16 populations from Europe and Japan. Confirming and extending results from the 1990s, we found that the SRD are present at variable frequencies in populations and that the f element is overall more frequent than Wolbachia. The two SRD never co-occur at high frequency in any population, suggesting an apparent mutual exclusion. We also detected Wolbachia or the f element in some males, which probably reflects insufficient titer to induce feminization or presence of masculinizing alleles. Our results are consistent with a single integration event of a Wolbachia genome in the A. vulgare genome at the origin of the f element, which contradicts an earlier hypothesis of frequent losses and gains. We identified strong linkage between Wolbachia strains and mitochondrial haplotypes, but no association between the f element and mitochondrial background. Our results open new perspectives on SRD evolutionary dynamics in A. vulgare, the evolution of genetic conflicts and their impact on the variability of sex determination systems.

Integrative taxonomy of the subfamily Orbioninae (Crustacea: Isopoda) based on mitochondrial and nuclear data with evidence that supports Epicaridea as a suborder.

Epicaridea is a group of isopods with high morphological diversity, reduction and loss of characters, and strong sexual dimorphism due to their parasitic lifestyles but their systematics is not well understood. Despite the use of nuclear and mitochondrial genes to test the phylogeny of many invertebrate groups, few molecular data from epicarideans are known, especially from the subfamily Orbioninae. Species in this group are obligate penaeoid shrimp parasites and the lack molecular data has hampered studies on the phylogeny of Orbioninae. To rectify this, mitochondrial and nuclear genes of 9 orbionine species are sequenced here. Compared to the isopod ground pattern, the sequences of orbionines seem to be more plastic near the control region and major translocations are located between rrns and cob. A phylogenetic analysis based on three data sets showed strong support for a monophyletic Orbioninae and that Epicaridea should be accepted at the rank of a suborder within Isopoda. The monophyly of Parapenaeon and Orbione is in doubt based on morphological and molecular data. The genus Parapenaeon is revised and a new genus Aparapenaeon is erected for Parapenaeon japonica and three closely related species.

Livoneca redmanii Leach, 1818 (Cymothoidae) a parasitic isopod infesting the gills of the European seabass, Dicentrarchus labrax (Linnaeus, 1758): morphological and molecular characterization study.

BACKGROUND: Prevalence, morphology, and molecular characteristics of isopodiosis in the European seabass, Dicentrarchus labrax, in Egypt were assessed using light and electron microscopy and polymerase chain reaction targeting the mitochondrial COI (cytochrome oxidase c) gene. RESULTS: Adult parasites were found mainly in the branchial cavity between gill arches and to a lesser extent in the buccal cavity. They were morphologically identified as the Cymothoidae Livoneca redmanii Leach, 1818 (Crustacea: Isopoda). Obviously, a 23% prevalence rate of isopods was reported in D. labrax from Egyptian Mediterranean waters. Destructive and degenerative necrotic alterations with complete sloughing of gill lamellae of the infested fish were observed. DNA sequencing of the mitochondrial COI gene confirmed the identification of the parasite which was deposited in the GenBank under accession numbers MW600099, MZ208984, and MZ208985. Furthermore, the phylogenetic analysis demonstrated that parasites emerged from a monophyletic clade closely affiliated with L. redmanii and were clearly distinguished from other isopod genospecies. CONCLUSION: The present investigation addresses L. redmanii infestation in D. labrax in Egypt and affirmed morphological properties via the scanning electron microscopy (SEM) and molecular characteristics of this isopod species. The drastic effects of this parasite on the infected fish were proven both clinically and histopathologically.

Diversification within an oceanic Mediterranean island: Insights from a terrestrial isopod.

Understanding intra-island patterns of evolutionary divergence, including cases of cryptic diversity, is a crucial step towards deciphering speciation processes. Cyprus is an oceanic island isolated for at least 5.3 Mya from surrounding continental regions, while it remains unclear whether it was ever connected to the mainland, even during the Messinian Salinity Crisis. The terrestrial isopod species Armadillo officinalis, that is widespread across the Mediterranean, offers the opportunity to explore intra-island divergence patterns that might exhibit geographical structure related also to the region's known paleogeography. Genome-wide ddRADseq, as well as Sanger sequencing for four mitochondrial and three nuclear loci data were generated for this purpose. In total, 71 populations from Cyprus, neighbouring continental sites, i.e., Israel, Lebanon and Turkey, and other Mediterranean regions, i.e. Greece, Italy, and Tunisia, were included in the analysis. Phylogenetic reconstructions and population structure analyses support the existence of at least six genetically discrete groups across the study area. Five of these distinct genetic clades occur on Cyprus, four of which are endemic to the island and one is widely distributed along the circum-Mediterranean countries. The sixth clade is distributed in Israel. The closest evolutionary relationship of endemic Cypriot populations is with those from Israel, while the evolutionary clade that is present in countries all around the Mediterranean is very shallow. Cladochronological analyses date the origin of the species on the island at ∼6 Mya. Estimated f4 and D statistics as well as FST values indicate the genetic isolation between the populations sampled from Cyprus and surrounding continental areas, while there is evident gene flow among populations within the island. Species delimitation and population genetic metrics support the existence of three distinct taxonomic units across the study area, two of which occur on the island and correspond to the endemic clade and the widespread circum-Mediterranean one, respectively, while the third corresponds to Israel's clade. The islands' paleogeographic history and recent human activities seem to have shaped current patterns of genetic diversity in this group of species.

Microbial Multitrophic Communities Drive the Variation of Antibiotic Resistome in the Gut of Soil Woodlice (Crustacea: Isopoda).

Multitrophic communities inhabit in soil faunal gut, including bacteria, fungi, and protists, which have been considered a hidden reservoir for antibiotic resistance genes (ARGs). However, there is a dearth of research focusing on the relationships between ARGs and multitrophic communities in the gut of soil faunas. Here, we studied the contribution of multitrophic communities to variations of ARGs in the soil woodlouse gut. The results revealed diverse and abundant ARGs in the woodlouse gut. Network analysis further exhibited strong connections between key ecological module members and ARGs, suggesting that multitrophic communities in the keystone ecological cluster may play a pivotal role in the variation of ARGs in the woodlouse gut. Moreover, long-term application of sewage sludge significantly altered the woodlice gut resistome and interkingdom communities. The variation portioning analysis indicated that the fungal community has a greater contribution to variations of ARGs than bacterial and protistan communities in the woodlice gut after long-term application of sewage sludge. Together, our results showed that changes in gut microbiota associated with agricultural practices (e.g., sewage sludge application) can largely alter the gut interkingdom network in ecologically relevant soil animals, with implications for antibiotic resistance, which advances our understanding of the microecological drivers of ARGs in terrestrial ecosystem.

Genome of a giant isopod, Bathynomus jamesi, provides insights into body size evolution and adaptation to deep-sea environment.

BACKGROUND: The deep-sea may be regarded as a hostile living environment, due to low temperature, high hydrostatic pressure, and limited food and light. Isopods, a species-rich group of crustaceans, are widely distributed across different environments including the deep sea and as such are a useful model for studying adaptation, migration, and speciation. Similar to other deep-sea organisms, giant isopods have larger body size than their shallow water relatives and have large stomachs and fat bodies presumably to store organic reserves. In order to shed light on the genetic basis of these large crustaceans adapting to the oligotrophic environment of deep-sea, the high-quality genome of a deep-sea giant isopod Bathynomus jamesi was sequenced and assembled. RESULTS: B. jamesi has a large genome of 5.89 Gb, representing the largest sequenced crustacean genome to date. Its large genome size is mainly attributable to the remarkable proliferation of transposable elements (84%), which may enable high genome plasticity for adaptive evolution. Unlike its relatives with small body size, B. jamesi has expanded gene families related to pathways of thyroid and insulin hormone signaling that potentially contribute to its large body size. Transcriptomic analysis showed that some expanded gene families related to glycolysis and vesicular transport were specifically expressed in its digestive organs. In addition, comparative genomics and gene expression analyses in six tissues suggested that B. jamesi has inefficient lipid degradation, low basal metabolic rate, and bulk food storage, suggesting giant isopods adopt a more efficient mechanism of nutrient absorption, storage, and utilization to provide sustained energy supply for their large body size. CONCLUSIONS: Taken together, the giant isopod genome may provide a valuable resource for understanding body size evolution and adaptation mechanisms of macrobenthic organisms to deep-sea environments.

Phylogenetic position of Bopyroides hippolytes, with comments on the rearrangement of the mitochondrial genome in isopods (Isopoda: Epicaridea: Bopyridae).

BACKGROUND: Classification of parasitic bopyrids has traditionally been based on morphological characteristics, but phylogenetic relationships have remained elusive due to limited information provided by morphological data and tendency for loss of morphological features as a result of parasitic lifestyle. Subfamily Argeiinae was separated from Bopyrinae based on morphological evidence, although the assignment of all genera has not been phylogenetically evaluated. Bopyroides hippolytes has been traditionally classified in Bopyrinae, but divergent morphological characters make this assignment questionable. To investigate the relationship of bopyrines, we sequenced the complete mitochondrial genome of B. hippolytes and four mitochondrial genes of two other Bopyrinae. RESULTS: The phylogenetic trees based on separate and combined cox1and 18S sequence data recovered Bopyridae as robustly monophyletic, but Bopyrinae as polyphyletic. Bopyroides hippolytes was a close sister to Argeia pugettensis, type species to Argeiinae. Mitochondrial phylogenomics also suggested that B. hippolytes was close to Argeiinae. We also found a novel gene order in B. hippolytes compared to other isopods. CONCLUSIONS: Bopyroides hippolytes should be excluded from the Bopyrinae and has a close affinity with Argeia pugettensis based on molecular and morphological data. The conserved syntenic blocks of mitochondrial gene order have distinctive characteristics at a subordinal level and may be helpful for understanding the higher taxonomic level relationships of Isopoda.

Molecular features of Probopyrus sp. (Isopoda: Bopyridae) from Brazilian Amazonia and the parasitism of inland populations of Macrobrachium amazonicum (Decapoda: Palaemonidae).

Bopyrid isopods of the genus Probopyrus are well-known parasites of freshwater prawns of the genus Macrobrachium. The parasitism of coastal populations of Macrobrachium amazonicum by Probopyrus bithynis, for example, has been documented since the late 1980s. Despite this, molecular data on different populations are not available for any Probopyrus species. The present study is the first to describe Probopyrus populations from distinct regions of the Amazon basin based on sequences of two genes, the mitochondrial cytochrome oxidase C subunit I (COI) and the nuclear 18S ribosomal DNA (18S rDNA) gene. The analyses indicated the presence of two Probopyrus species, each parasitizing either the coastal or the inland populations of M. amazonicum. The results indicated the potential use of the COI barcode for the identification of Probopyrus species. We discuss the potential implications of the findings for the taxonomy of Probopyrus bithynis and other species of the genus Probopyrus.

Recent speciation and hybridization in Icelandic deep-sea isopods: An integrative approach using genomics and proteomics.

The crustacean marine isopod species Haploniscus bicuspis (Sars, 1877) shows circum-Icelandic distribution in a wide range of environmental conditions and along well-known geographic barriers, such as the Greenland-Iceland-Faroe (GIF) Ridge. We wanted to explore population genetics, phylogeography and cryptic speciation as well as investigate whether previously described, but unaccepted subspecies have any merit. Using the same set of specimens, we combined mitochondrial COI sequences, thousands of nuclear loci (ddRAD), and proteomic profiles, plus selected morphological characters using confocal laser scanning microscopy (CLSM). Five divergent genetic lineages were identified by COI and ddRAD, two south and three north of the GIF Ridge. Assignment of populations to the three northern lineages varied and detailed analyses revealed hybridization and gene flow between them, suggesting a single northern species with a complex phylogeographic history. No apparent hybridization was observed among lineages south of the GIF Ridge, inferring the existence of two more species. Differences in proteomic profiles between the three putative species were minimal, implying an ongoing or recent speciation process. Population differentiation was high, even among closely associated populations, and higher in mitochondrial COI than nuclear ddRAD loci. Gene flow is apparently male-biased, leading to hybrid zones and instances of complete exchange of the local nuclear genome through immigrating males. This study did not confirm the existence of subspecies defined by male characters, which probably instead refer to different male developmental stages.

A new species of Expanathura (Crustacea: Isopoda: Anthuroidea) from Iriomote Island, Japan, with a note on male polymorphism.

We describe Expanathura monile sp. nov. from Iriomote Island, Okinawa, Japan, northwestern Pacific. This species resembles E. collaris, E. macronesia, and E. haddae in having very broad uropodal rami, a broad telson, and a moderately shortened pleopod 1 endopod, but differs from these three species in the female having (1) the head with a dorsal V-shaped band of brown pigmentation, (2) the antennular flagellum with 34 articles, (3) antennal peduncular article 2 with an outer triangular projection, (4) the uropodal endopod longer than wide, (5) the uropodal exopod with an acute triangular projection, (6) an oval telson, and (7) the posterolateral setae on the telson similar in length. We present a revised key to known species in the genus Expanathura and briefly discuss the male polymorphism observed in E. monile.

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

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

Evolutionary transition to XY sex chromosomes associated with Y-linked duplication of a male hormone gene in a terrestrial isopod.

Sex chromosomes are highly variable in some taxonomic groups, but the evolutionary mechanisms underlying this diversity are not well understood. In terrestrial isopod crustaceans, evolutionary turnovers in sex chromosomes are frequent, possibly caused by Wolbachia, a vertically-transmitted endosymbiont causing male-to-female sex reversal. Here, we use surgical manipulations and genetic crosses, plus genome sequencing, to examine sex chromosomes in the terrestrial isopod Trachelipus rathkei. Although an earlier cytogenetics study suggested a ZZ/ZW sex chromosome system in this species, we surprisingly find multiple lines of evidence that in our study population, sex is determined by an XX/XY system. Consistent with a recent evolutionary origin for this XX/XY system, the putative male-specific region of the genome is small. The genome shows evidence of Y-linked duplications of the gene encoding the androgenic gland hormone, a major component of male sexual differentiation in isopods. Our analyses also uncover sequences horizontally acquired from past Wolbachia infections, consistent with the hypothesis that Wolbachia may have interfered with the evolution of sex determination in T. rathkei. Overall, these results provide evidence for the co-occurrence of multiple sex chromosome systems within T. rathkei, further highlighting the relevance of terrestrial isopods as models for the study of sex chromosome evolution.