Research on the Mechanism of Metaldehyde on Pomacea canaliculata.

Golden apple snail (Pomacea canaliculata), a major alien invasive organism in China, affects food production and poses a threat to human health. Metaldehyde is a highly effective, commonly used snail killer with low toxicity. Virulence determination, tissue section, iTRAQ and RNA interference were used to systematically study the toxicity of metaldehyde on P. canaliculata. The molluscicidal activity tests showed that metaldehyde exhibits strong toxicity against P. canaliculata. Physiological and biochemical data indicate that metaldehyde can cause damage to the gills, liver, pancreas, and kidneys of snails, also reduce the oxygen consumption rate and ammonia excretion rate of golden apple snails, and cause neurological diseases. The proteome of the gill region of the golden apple snail after exposure to metaldehyde was analyzed by using iTRAQ technology. A total of 360 differential proteins were identified, and four target proteins were screened, namely, alpha-protein kinase 1 (ALPK1), cubilin (CUBN), sodium- and chloride-dependent GABA transporter 2 (GAT2), and acetylcholinesterase (AChE). RNAi was used to target the four proteins. After the ALPK1 and CUBN protein genes were interfered with by metaldehyde treatment, it was found that the mortality rate of the golden apple snail significantly increased. However, interference of GAT2 and AChE protein genes by metaldehyde led to no significant change in the mortality rates of the snails. The histopathological observation of the gill showed that the rate of cilia shedding in the gill decreased after the interference of ALPK1 and CUBN protein genes.

Odd-Paired is Involved in Morphological Divergence of Snail-Feeding Beetles.

Body shape and size diversity and their evolutionary rates correlate with species richness at the macroevolutionary scale. However, the molecular genetic mechanisms underlying the morphological diversification across related species are poorly understood. In beetles, which account for one-fourth of the known species, adaptation to different trophic niches through morphological diversification appears to have contributed to species radiation. Here, we explored the key genes for the morphological divergence of the slender to stout body shape related to divergent feeding methods on large to small snails within the genus Carabus. We show that the zinc-finger transcription factor encoded by odd-paired (opa) controls morphological variation in the snail-feeding ground beetle Carabus blaptoides. Specifically, opa was identified as the gene underlying the slender to stout morphological difference between subspecies through genetic mapping and functional analysis via gene knockdown. Further analyses revealed that changes in opa cis-regulatory sequences likely contributed to the differences in body shape and size between C. blaptoides subspecies. Among opa cis-regulatory sequences, single nucleotide polymorphisms on the transcription factor binding sites may be associated with the morphological differences between C. blaptoides subspecies. opa was highly conserved in a wide range of taxa, especially in beetles. Therefore, opa may play an important role in adaptive morphological divergence in beetles.

Littoporins: Novel actinoporin-like proteins in caenogastropod genus Littorina.

In the course of searching for genes controlling the immune system in caenogastropod mollusks, we characterized and phylogenetically placed five new actinoporin-like cytolysins expressed in periwinkles of the genus Littorina. These newly discovered proteins, named littoporins (LitP), contain a central cytolysin/lectin domain and exhibit a predicted protein fold that is almost identical to the three-dimensional structures of actinoporins. Two of these proteins, LitP-1 and LitP-2, were found to be upregulated in L. littorea kidney tissues and immune cells in response to natural and experimental infection with the trematode Himasthla elongata, suggesting their potential role as perforins in the systemic anti-trematode immune response. The primary sequence divergence of littoporins is hypothesized to be attributed to the taxonomic range of cell membranes they can recognize and permeabilize.

Genetic diversity and population structure of Physella acuta (Gastropoda: Physidae) in Thailand using mitochondrial gene markers: COI and 16S rDNA.

Physella acuta is a freshwater snail native to North America. Understanding the phylogeography and genetic structure of P. acuta will help elucidate its evolution. In this study, we used mitochondrial (COI and 16S rDNA) and nuclear (ITS1) markers to identify the species and examine its genetic diversity, population structure, and demographic history of P. acuta in Thailand. Phylogenetic and network analyses of P. acuta in Thailand pertained to clade A, which exhibits a global distribution. Analysis of the genetic structure of the population revealed that the majority of pairwise comparisons showed no genetic dissimilarity. An isolation-by-distance test indicates no significant correlation between genetic and geographical distances among P. acuta populations, suggesting that gene flow is not restricted by distance. Demographic history and haplotype network analyses suggest a population expansion of P. acuta, as evidenced by the star-like structure detected in the median-joining network. Based on these results, we concluded that P. acuta in Thailand showed gene flow and recent population expansion. Our findings provide fundamental insights into the genetic variation of P. acuta in Thailand.

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis.

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers.

BACKGROUND: The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult. OBJECTIVE: In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi. RESULTS: After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT. CONCLUSION: The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

Genome Evolution and Introgression in the New Zealand mud Snails Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa.

We have sequenced, assembled, and analyzed the nuclear and mitochondrial genomes and transcriptomes of Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa, two prosobranch snail species native to New Zealand that together span the continuum from estuary to freshwater. These two species are the closest known relatives of the freshwater species Potamopyrgus antipodarum-a model for studying the evolution of sex, host-parasite coevolution, and biological invasiveness-and thus provide key evolutionary context for understanding its unusual biology. The P. estuarinus and P. kaitunuparaoa genomes are very similar in size and overall gene content. Comparative analyses of genome content indicate that these two species harbor a near-identical set of genes involved in meiosis and sperm functions, including seven genes with meiosis-specific functions. These results are consistent with obligate sexual reproduction in these two species and provide a framework for future analyses of P. antipodarum-a species comprising both obligately sexual and obligately asexual lineages, each separately derived from a sexual ancestor. Genome-wide multigene phylogenetic analyses indicate that P. kaitunuparaoa is likely the closest relative to P. antipodarum. We nevertheless show that there has been considerable introgression between P. estuarinus and P. kaitunuparaoa. That introgression does not extend to the mitochondrial genome, which appears to serve as a barrier to hybridization between P. estuarinus and P. kaitunuparaoa. Nuclear-encoded genes whose products function in joint mitochondrial-nuclear enzyme complexes exhibit similar patterns of nonintrogression, indicating that incompatibilities between the mitochondrial and the nuclear genome may have prevented more extensive gene flow between these two species.

Freshwater trematodes differ from marine trematodes in patterns connected with division of labor.

Background: Prior research suggests that trematode rediae, a developmental stage of trematode parasites that reproduce clonally within a snail host, show evidence of division of labor (DOL). Single-species infections often have two morphologically distinct groups: small rediae, the 'soldiers', are active, aggressive, and do not appear to reproduce; large rediae, the 'reproductives', are larger, sluggish, and full of offspring. Most data supporting DOL come from trematodes infecting marine snails, while data from freshwater trematodes are more limited and generally do not supported DOL. The shorter lifespan typical of freshwater snails may partially explain this difference: defending a short-lived host at the expense of reproduction likely provides few advantages. Here, we present data from sixty-one colonies spanning twenty species of freshwater trematode exploring morphological and behavioral patterns commonly reported from marine trematodes believed to have DOL. Methods: Trematode rediae were obtained from sixty-one infected snails collected in central Vermont, USA. A portion of the COI gene was sequenced to make tentative species identifications ('COI species'). Samples of rediae were photographed, observed, and measured to look for DOL-associated patterns including a bimodal size distribution, absence of embryos in small rediae, and pronounced appendages and enlarged pharynges (mouthparts) in small rediae. Additional rediae were used to compare activity levels and likelihood to attack heterospecific trematodes in large vs. small rediae. Results: Many of the tests for DOL-associated patterns showed mixed results, even among colonies of the same COI species. However, we note a few consistent patterns. First, small rediae of most colonies appeared capable of reproduction, and we saw no indication (admittedly based on a small sample size and possibly insufficient attack trial methodology) that small rediae were more active or aggressive. This differs from patterns reported from most marine trematodes. Second, the small rediae of most colonies had larger pharynges relative to their body size than large rediae, consistent with marine trematodes. We also observed that colonies of three sampled COI species appear to produce a group of large rediae that have distinctly large pharynges. Conclusions: We conclude that these freshwater species likely do not have a group of specialized non-reproductive soldiers because small rediae of at least some colonies in almost every species do appear to produce embryos. We cannot rule out the possibility that small rediae act as a temporary soldier caste. We are intrigued by the presence of rediae with enlarged pharynges in some species and propose that they may serve an adaptive role, possibly similar to the defensive role of small 'soldier' rediae of marine trematodes. Large-pharynx rediae have been documented in other species previously, and we encourage future efforts to study these large-pharynx rediae.

Melanella martarum sp. nov. (Gastropoda: Eulimidae): the first parasitic deep-sea snail reported for the Salas & Gomez Ridge.

Eulimidae is a highly diverse family of gastropods that are often parasites of echinoderms. They are cosmopolitan and live from the intertidal to great depths. Despite its wide geographic and bathymetric distribution, no species of Eulimidae have been reported for the Salas & Gómez Ridge to date. In this study, we describe Melanella martarum sp. nov., which was collected during the EPIC oceanographic cruise onboard RV Mirai (JAMSTEC, Japan) in 2019. Seven specimens were collected with a modified Agassiz trawl on the summit of seamount "Pearl" (Zhemchuznaya) in the Salas & Gómez Ridge (25.59°S, 89.13°W) at 545 m depth. The morphology of M. martarum sp. nov. was compared with other Melanella species reported for the area, including Chile and Rapa Nui. DNA was extracted and partial sequences of the mitochondrial genes Cytochrome Oxidase 1 (COI) and 16S rDNA, and the nuclear gene Histone 3 (H3) were sequenced. Melanella martarum sp. nov. has morphological characteristics that separate it from other species of Melanella, such as the thickness and color of the shell, and the shape of the protoconch. In addition, M. martarum sp. nov. was genetically differentiated from other Melanella spp. sequences (uncorrected p distances from 18,1-8.6% in mitochondrial COI and 16S rDNA to 3% in nuclear H3 sequences). Although there is not much molecular data available for Eulimidae, the phylogenetic analysis confirms the results obtained by morphology, placing the species found on the Salas & Gómez Ridge within the genus Melanella. The current study advances the understanding of the poorly known benthic fauna found on seamounts in the easternmost part of the Sala & Gómez ridge, a location distinguished by a high level of endemism.

Morphological and molecular characteristics of Paralecithodendrium longiforme (Digenea: Lecithodendriidae) adults and cercariae from Chinese pipistrelle bats and viviparid snails in Thailand.

This study aimed to describe the morphological and molecular characteristics of Paralecithodendrium longiforme (Digenea: Lecithodendriidae) adults and cercariae isolated in Thailand. Adult flukes were isolated from the Chinese pipistrelle bat (Hypsugo sp.), and cercariae were detected in the viviparid snail (Filopaludina martensi martensi) from Chiang Mai province. The morphological characteristics were observed and described using conventional methods, and the molecular characteristics with internal transcribed spacer 2 (ITS2) and 28S rDNA gene sequences. The adult flukes were fusiform, 0.84-0.98 mm in length, and 0.37-0.49 mm in width, and were distinguishable from other species by the presence of longitudinal uterine coils. The cercariae were nonvirgulate xiphidiocercariae, with the oral sucker bigger than the acetabulum, the tail without fin fold, a body size of 117.5-138.3 × 48.3-52.2 µm, and a tail size of 100.7-103.7 × 15.0-18.9 µm. Molecular studies revealed that the adults and cercariae shared 99.3% (ITS2) and 99.6% (28S rDNA) homology with each other. They were phylogenetically close to P. longiforme with an identity of 94.5% for ITS2 and 98.7% for 28S rDNA. This study provides new information on the natural definitive host and first intermediate host of P. longiforme in Thailand. The discovery of its cercarial stage in Filopaludina snails highlights the importance of monitoring the associated second intermediate host and prevention and control of this potentially zoonotic trematode.

Both nuclear and cytoplasmic polymorphisms are involved in genetic conflicts over male fertility in the gynodioecious snail, Physa acuta.

Gynodioecy, the coexistence of hermaphrodites with females, often reflects conflicts between cytoplasmic male sterility (CMS) genes and nuclear genes restoring male fertility. CMS is frequent in plants and has been recently discovered in one animal: the freshwater snail, Physa acuta. In this system, CMS was linked to a single divergent mitochondrial genome (D), devoid of apparent nuclear restoration. Our study uncovers a second, novel CMS-associated mitogenome (K) in Physa acuta, demonstrating an extraordinary acceleration of molecular evolution throughout the entire K mitochondrial genome, akin to the previously observed pattern in D. This suggests a pervasive occurrence of accelerated evolution in both CMS-associated lineages. Through a 17-generation introgression experiment, we further show that nuclear polymorphisms in K-mitogenome individuals contribute to the restoration of male function in natural populations. Our results underscore shared characteristics in gynodioecy between plants and animals, emphasizing the presence of multiple CMS mitotypes and cytonuclear conflicts. This reaffirms the pivotal role of mitochondria in influencing male function and in generating genomic conflicts that impact reproductive processes in animals.

Two new genera and species of the valvatiform hydrobiid snails (Caenogastropoda: Truncatelloidea: Hydrobiidae) from Morocco.

Stygobiont and crenobiont minute gastropods representing the family Hydrobiidae (Caenogastropoda: Truncatelloidea), characterized by the valvatoid low-spired shell, were collected from one spring and four wells in Bouregreg region, NW Morocco. The shells were photographed and measured; shell biometry is illustrated with principal component analysis. Penes were illustrated and described. Mitochondrial cytochrome oxidase subunit 1 (COI) and 16S rRNA, as well as nuclear 18S rRNA sequences were used to infer the phylogeny. The snails represented two genera, both new to science. Their closest relatives were taxa from the Iberian Peninsula, the rough molecular estimate of the time of divergence between these Moroccan and Iberian genera coincided with the Pliocene Flooding, which restored the Strait of Gibraltar to connect the Mediterranean Basin with Atlantic Ocean.

Plasticity of circadian and circatidal rhythms in activity and transcriptomic dynamics in a freshwater snail.

Organisms have diverse biological clocks synchronised with environmental cycles depending on their habitats. Anticipation of tidal changes has driven the evolution of circatidal rhythms in some marine species. In the freshwater snail, Semisulcospira reiniana, individuals in nontidal areas exhibit circadian rhythms, whereas those in tidal areas exhibit both circadian and circatidal rhythms. We investigated whether the circatidal rhythms are genetically determined or induced by environmental cycles. The exposure to a simulated tidal cycle did not change the intensity of circatidal rhythm in individuals in the nontidal population. However, snails in the tidal population showed different activity rhythms depending on the presence or absence of the exposure. Transcriptome analysis revealed that genes with circatidal oscillation increased due to entrainment to the tidal cycle in both populations and dominant rhythmicity was consistent with the environmental cycle. These results suggest plasticity in the endogenous rhythm in the gene expression in both populations. Note that circatidal oscillating genes were more abundant in the tidal population than in the nontidal population, suggesting that a greater number of genes are associated with circatidal clocks in the tidal population compared to the nontidal population. This increase of circatidal clock-controlled genes in the tidal population could be caused by genetic changes in the biological clock or the experience of tidal cycle in the early life stage. Our findings suggest that the plasticity of biological rhythms may have contributed to the adaptation to the tidal environment in S. reiniana.

Gastropod invasions in anthropogenically impacted impoundments in South Africa: Tracing their origins and exploring field evidence of parasite spillback and amplification.

Invasive snails are associated with ecological problems in freshwater bodies worldwide. However, their impact on the transmission of digenean infections remain underreported. In the present study, 1708 specimens representing four snail species were sampled from four impoundments in the Limpopo River system in South Africa. Gyraulus chinensis (Planorbidae), Physella acuta (Physidae) and Pseudosuccinea columella (Lymnaeidae), which are invasive, were found in all the sampling sites. In contrast, the native lymnaeid Radix natalensis occurred at only one study site. Digeneans were observed only from R. natalensis (prevalence = 49%) and Ps. columella (prevalence = 23%). Morphological and genetic analyses revealed four digeneans: Fasciola nyanzae, Orientocreadium sp., Petasiger sp. and Patagifer vioscai. Pseudosuccinea columella was infected by the four digeneans while R. natalensis harboured only Orientocreadium sp. and Petasiger sp. Partial sequences of Orientocreadium sp. from the current study differed from congeners whose DNA data are available on GenBank, by p-distances of at least 1.84 and 2.2% for 28S and the internal transcribed spacer (ITS) rDNA, respectively. Phylogenetic analyses demonstrated that the present species is sister to Orientocreadium batrachoides. Genetic and phylogenetic data based on 28S and ITS rDNA suggested that Petasiger sp. from the present study and isolates of three unidentified Petasiger spp. from Kenya, Hungary and Australia, were representatives of the same species. This is the first known report of Orientocreadium, Petasiger and Patagifer from Ps. columella. The occurrence of F. nyanzae in Ps. columella indicates spillback from R. natalensis. These findings echo the concerns raised in previous studies about the potential role of Ps. columella in the amplification of digenean diseases in its introduced range. Phylogenetic analyses of partial sequences of the cytochrome c oxidase subunit 1 mitochondrial gene (cox1) showed multiple lineages of Ps. columella in North and South America. Pseudosuccinea columella specimens from the present study belong to an invasive genotype that has spread globally and has been reported from Zimbabwe, Egypt, Portugal, Australia, Argentina, Colombia and New Mexico (USA). Physella acuta from the current study had a stronger genetic relationship with isolates from Canada and Iceland, than with isolates from other parts of Africa, suggesting several invasion routes into Africa. This is the first known DNA characterisation of G. chinensis from Africa. Phylogenetic reconstruction indicated multiple exit events of G. chinensis from Asia into Europe and Africa. South African isolates clustered in a recent branch containing isolates from the Czech Republic and Hong Kong, China. Considering the presence of invasive snails in all the sampling sites in the present study, it is necessary to investigate the factors that enhance their establishment and to monitor their effects on the native snail populations.

Diverse pathways to speciation revealed by marine snails.

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

Comparative Mitogenome Analyses of Fifteen Ramshorn Snails and Insights into the Phylogeny of Planorbidae (Gastropoda: Hygrophila).

Ramshorn snails from the family Planorbidae are important freshwater snails due to their low trophic level, and some of them act as intermediate hosts for zoonotic trematodes. There are about 250 species from 40 genera of Planorbidae, but only 14 species from 5 genera (Anisus, Biomphalaria, Bulinus, Gyraulus, and Planorbella) have sequenced complete mitochondrial genomes (mitogenomes). In this study, we sequenced and assembled a high-quality mitogenome of a ramshorn snail, Polypylis sp. TS-2018, which represented the first mitogenome of the genus. The mitogenome of Polypylis sp. TS-2018 is 13,749 bp in length, which is shorter than that of most gastropods. It contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and 2 ribosomal RNA (rRNA). We compared mitogenome characteristics, selection pressure, and gene rearrangement among all of the available mitogenomes of ramshorn snails. We found that the nonsynonymous and synonymous substitution rates (Ka/Ks) of most PCGs indicated purifying and negative selection, except for atp8 of Anisus, Biomphalaria, and Gyraulus, which indicated positive selection. We observed that transpositions and reverse transpositions occurred on 10 tRNAs and rrnS, which resulted in six gene arrangement types. We reconstructed the phylogenetic trees using the sequences of PCGs and rRNAs and strongly supported the monophyly of each genus, as well as three tribes in Planorbidae. Both the gene rearrangement and phylogenetic results suggested that Polypylis had a close relationship with Anisus and Gyraulus, while Bulinus was the sister group to all of the other genera. Our results provide useful data for further investigation of species identification, population genetics, and phylogenetics among ramshorn snails.

Many roads to success: alternative routes to building an economic shell in land snails.

Land snails exhibit an extraordinary variety of shell shapes. The way shells are constructed underlies biological and mechanical constraints that vary across gastropod clades. Here, we quantify shell geometry of the two largest groups, Stylommatophora and Cyclophoroidea, to assess the potential causes for variation in shell shape and its relative frequency. Based on micro-computed tomography scans, we estimate material efficiency through 2D and 3D generalizations of the isoperimetric ratio, quantifying the ratios between area and perimeter of whorl cross-sections (2D) and shell volume and surface (3D), respectively. We find that stylommatophorans optimize material usage through whorl overlap, which may have promoted the diversification of flat-shelled species. Cyclophoroids are bound to a circular cross-section because of their operculum; flat shells are comparatively rare. Both groups show similar solutions for tall shells, where local geometry has a smaller effect because of the double overlap between previous and current whorls. Our results suggest that material efficiency is a driving factor in the selection of shell geometry. Essentially, the evolutionary success of Stylommatophora likely roots in their higher flexibility to produce an economic shell.

Evolutionary paths to new phenotypes.

Ecological model systems inform on innovative traits in plants and animals.

The genetic basis of a recent transition to live-bearing in marine snails.

Key innovations are fundamental to biological diversification, but their genetic basis is poorly understood. A recent transition from egg-laying to live-bearing in marine snails (Littorina spp.) provides the opportunity to study the genetic architecture of an innovation that has evolved repeatedly across animals. Individuals do not cluster by reproductive mode in a genome-wide phylogeny, but local genealogical analysis revealed numerous small genomic regions where all live-bearers carry the same core haplotype. Candidate regions show evidence for live-bearer-specific positive selection and are enriched for genes that are differentially expressed between egg-laying and live-bearing reproductive systems. Ages of selective sweeps suggest that live-bearer-specific alleles accumulated over more than 200,000 generations. Our results suggest that new functions evolve through the recruitment of many alleles rather than in a single evolutionary step.

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae.

BACKGROUND: The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes. RESULTS: Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs. CONCLUSION: The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.