LlamaTags: A Versatile Tool to Image Transcription Factor Dynamics in Live Embryos.

Embryonic cell fates are defined by transcription factors that are rapidly deployed, yet attempts to visualize these factors in vivo often fail because of slow fluorescent protein maturation. Here, we pioneer a protein tag, LlamaTag, which circumvents this maturation limit by binding mature fluorescent proteins, making it possible to visualize transcription factor concentration dynamics in live embryos. Implementing this approach in the fruit fly Drosophila melanogaster, we discovered stochastic bursts in the concentration of transcription factors that are correlated with bursts in transcription. We further used LlamaTags to show that the concentration of protein in a given nucleus heavily depends on transcription of that gene in neighboring nuclei; we speculate that this inter-nuclear signaling is an important mechanism for coordinating gene expression to delineate straight and sharp boundaries of gene expression. Thus, LlamaTags now make it possible to visualize the flow of information along the central dogma in live embryos.

Twist-dependent ratchet functioning downstream from Dorsal revealed using a light-inducible degron.

Graded transcription factors are pivotal regulators of embryonic patterning, but whether their role changes over time is unclear. A light-regulated protein degradation system was used to assay temporal dependence of the transcription factor Dorsal in dorsal-ventral axis patterning of Drosophila embryos. Surprisingly, the high-threshold target gene snail only requires Dorsal input early but not late when Dorsal levels peak. Instead, late snail expression can be supported by action of the Twist transcription factor, specifically, through one enhancer, sna.distal This study demonstrates that continuous input is not required for some Dorsal targets and downstream responses, such as twist, function as molecular ratchets.

The Snail signaling branch downstream of the TGF-ß/Smad3 pathway mediates Rho activation and subsequent stress fiber formation.

Cancer cells acquire malignant phenotypes through an epithelial-mesenchymal transition, which is induced by environmental factors or extracellular signaling molecules, including transforming growth factor-ß (TGF-ß). Among epithelial-mesenchymal transition-associated cell responses, cell morphological changes and cell motility are closely associated with remodeling of the actin stress fibers. Here, we examined the TGF-ß signaling pathways leading to these cell responses. Through knockdown experiments in A549 lung adenocarcinoma cells, we found that Smad3-mediated induction of Snail, but not that of Slug, is indispensable for morphological changes, stress fiber formation, and enhanced motility in cells stimulated with TGF-ß. Ectopic expression of Snail in SMAD3-knockout cells rescued the defect in morphological changes and stress fiber formation by TGF-ß, indicating that the role of Smad3 in these responses is to upregulate Snail expression. Mechanistically, Snail is required for TGF-ß-induced upregulation of Wnt5b, which in turn activates RhoA and subsequent stress fiber formation in cooperation with phosphoinositide 3-kinase. However, ectopic expression of Snail in SMAD3-knockout cells failed to rescue the defect in cell motility enhancement by TGF-ß, indicating that activation of the Smad3/Snail/Wnt5b axis is indispensable but not sufficient for enhancing cell motility; a Smad3-dependent but Snail-independent pathway to activate Rac1 is additionally required. Therefore, the Smad3-dependent pathway leading to enhanced cell motility has two branches: a Snail-dependent branch to activate RhoA and a Snail-independent branch to activate Rac1. Coordinated activation of these branches, together with activation of non-Smad signaling pathways, mediates enhanced cell motility induced by TGF-ß.

Prey Shifts Drive Venom Evolution in Cone Snails.

Venom systems are complex traits that have independently emerged multiple times in diverse plant and animal phyla. Within each venomous lineage there typically exists interspecific variation in venom composition where several factors have been proposed as drivers of variation, including phylogeny and diet. Understanding these factors is of broad biological interest and has implications for the development of antivenom therapies and venom-based drug discovery. Because of their high species richness and the presence of several major evolutionary prey shifts, venomous marine cone snails (genus Conus) provide an ideal system to investigate drivers of interspecific venom variation. Here, by analyzing the venom gland expression profiles of ∼3,000 toxin genes from 42 species of cone snail, we elucidate the role of prey-specific selection pressures in shaping venom variation. By analyzing overall venom composition and individual toxin structures, we demonstrate that the shifts from vermivory to piscivory in Conus are complemented by distinct changes in venom composition independent of phylogeny. In vivo injections of venom from piscivorous cone snails in fish further showed a higher potency compared with venom of nonpiscivores demonstrating a selective advantage. Together, our findings provide compelling evidence for the role of prey shifts in directing the venom composition of cone snails and expand our understanding of the mechanisms of venom variation and diversification.

Dynamic EMT: a multi-tool for tumor progression.

The process of epithelial-mesenchymal transition (EMT) is fundamental for embryonic morphogenesis. Cells undergoing it lose epithelial characteristics and integrity, acquire mesenchymal features, and become motile. In cancer, this program is hijacked to confer essential changes in morphology and motility that fuel invasion. In addition, EMT is increasingly understood to orchestrate a large variety of complementary cancer features, such as tumor cell stemness, tumorigenicity, resistance to therapy and adaptation to changes in the microenvironment. In this review, we summarize recent findings related to these various classical and non-classical functions, and introduce EMT as a true tumorigenic multi-tool, involved in many aspects of cancer. We suggest that therapeutic targeting of the EMT process will-if acknowledging these complexities-be a possibility to concurrently interfere with tumor progression on many levels.

Noncanonical Wnt signaling promotes colon tumor growth, chemoresistance and tumor fibroblast activation.

Colon tumors of the mesenchymal subtype have the lowest overall survival. Snail1 is essential for the acquisition of this phenotype, characterized by increased tumor stemness and invasion, and high resistance to chemotherapy. Here, we find that Snail1 expression in colon tumor cells is dependent on an autocrine noncanonical Wnt pathway. Accordingly, depletion of Ror2, the co-receptor for noncanonical Wnts such as Wnt5a, potently decreases Snail1 expression. Wnt5a, Ror2, and Snail1 participate in a self-stimulatory feedback loop since Wnt5a increases its own synthesis in a Ror2- and Snail1-dependent fashion. This Wnt5a/Ror2/Snail1 axis controls tumor invasion, chemoresistance, and formation of tumor spheres. It also stimulates TGFß synthesis; consequently, tumor cells expressing Snail1 are more efficient in activating cancer-associated fibroblasts than the corresponding controls. Ror2 downmodulation or inhibition of the Wnt5a pathway decreases Snail1 expression in primary colon tumor cells and their ability to form tumors and liver metastases. Finally, the expression of SNAI1, ROR2, and WNT5A correlates in human colon and other tumors. These results identify inhibition of the noncanonical Wnt pathway as a putative colon tumor therapy.

Deubiquitinating enzyme OTUD4 regulates metastasis in triple-negative breast cancer by stabilizing Snail1.

Metastasis is the primary cause of cancer-related deaths and remains poorly understood. Deubiquitinase OTU domain containing 4 (OTUD4) has been reported to regulate antiviral immune responses and resistance to radio- or chemo-therapies in certain cancers. However, the role of OTUD4 in cancer metastasis remain unknown. Here, we demonstrate that the depletion of OTUD4 in triple-negative breast cancer (TNBC) cells markedly suppress cell clonogenic ability, migration, invasion and cancer stem cell population in vitro as well as metastasis in vivo. Mechanistically, the tumor promoting function of OTUD4 is mainly mediated by deuiquitinating and stabilizing Snail1, one key transcriptional factor in the epithelial-mesenchymal transition. The inhibitory effect of targeting OTUD4 could be largely reversed by the reconstitution of Snail1 in OTUD4-deficient cells. Overall, our study establishes the OTUD4-Snail1 axis as an important regulatory mechanism of breast cancer metastasis and provides a rationale for potential therapeutic interventions in the treatment of TNBC.

Low PDE4A expression promoted the progression of ovarian cancer by inducing Snail nuclear translocation.

Widespread metastasis is the primary reason for the high mortality associated with ovarian cancer (OC), and effective targeted therapy for tumor aggressiveness is still insufficient in clinical practice. Therefore, it is urgent to find new targets to improve prognosis of patients. PDE4A is a cyclic nucleotide phosphodiesterase that plays a crucial role in the occurrence and development in various malignancies. Our study firstly reported the function of PDE4A in OC. Expression of PDE4A was validated through bioinformatics analysis, RT-qPCR, Western blot, and immunohistochemistry. Additionally, its impact on cell growth and motility was assessed via in vitro and in vivo experiments. PDE4A was downregulated in OC tissues compared with normal tissues and low PDE4A expression was correlated with poor clinical outcomes in OC patients. The knockdown of PDE4A significantly promoted the proliferation, migration and invasion of OC cells while overexpression of PDE4A resulted in the opposite effect. Furthermore, smaller and fewer tumor metastatic foci were observed in mice bearing PDE4A-overexpressing OVCAR3 cells. Mechanistically, downregulation of PDE4A expression can induce epithelial-mesenchymal transition (EMT) and nuclear translocation of Snail, which suggests that PDE4A plays a pivotal role in suppressing OC progression. Notably, Rolipram, the PDE4 inhibitor, mirrored the effects observed with PDE4A deletion. In summary, the downregulation of PDE4A appears to facilitate OC progression by modulating the Snail/EMT pathway, underscoring the potential of PDE4A as a therapeutic target against ovarian cancer metastasis.

A coherent FOXO3-SNAI2 feed-forward loop in autophagy.

SignificanceUnderstanding autophagy regulation is instrumental in developing therapeutic interventions for autophagy-associated disease. Here, we identified SNAI2 as a regulator of autophagy from a genome-wide screen in HeLa cells. Upon energy stress, SNAI2 is transcriptionally activated by FOXO3 and interacts with FOXO3 to form a feed-forward regulatory loop to reinforce the expression of autophagy genes. Of note, SNAI2-increased FOXO3-DNA binding abrogates CRM1-dependent FOXO3 nuclear export, illuminating a pivotal role of DNA in the nuclear retention of nucleocytoplasmic shuttling proteins. Moreover, a dFoxO-Snail feed-forward loop regulates both autophagy and cell size in Drosophila, suggesting this evolutionarily conserved regulatory loop is engaged in more physiological activities.

Cytolethal Distending Toxin Modulates Cell Differentiation and Elicits Epithelial to Mesenchymal Transition.

BACKGROUND: The bacterial genotoxin, cytolethal distending toxin (CDT), causes DNA damage in host cells, a risk factor for carcinogenesis. Previous studies have shown that CDT induces phenotypes reminiscent of epithelial to mesenchymal transition (EMT), a process involved in cancer initiation and progression. METHODS: We investigated different steps of EMT in response to Helicobacter hepaticus CDT and its active CdtB subunit using in vivo and in vitro models. RESULTS: Most of the steps of the EMT process were induced by CDT/CdtB and observed throughout the study in murine and epithelial cell culture models. CdtB induced cell-cell junction disassembly, causing individualization of cells and acquisition of a spindle-like morphology. The key transcriptional regulators of EMT (SNAIL and ZEB1) and some EMT markers were upregulated at both RNA and protein levels in response to CDT/CdtB. CdtB increased the expression and proteolytic activity of matrix metalloproteinases, as well as cell migration. A range of these results were confirmed in Helicobacter hepaticus-infected and xenograft murine models. In addition, colibactin, a genotoxic metabolite produced by Escherichia coli, induced EMT-like effects in cell culture. CONCLUSIONS: Overall, these data show that infection with genotoxin-producing bacteria elicits EMT process activation, supporting their role in tumorigenesis.

Escargot a Snail superfamily member and its multiple roles in Drosophila melanogaster development.

The Snail superfamily of transcription factors plays a crucial role in metazoan development; one of the most important vertebrate members of this family is Snai1 which is orthologous to the Drosophila melanogaster esg gene. This review offers a comprehensive examination of the roles of the esg gene in Drosophila development, covering its expression pattern and downstream targets, and draws parallels between the vertebrate Snai1 family proteins on controlling the epithelial-to-mesenchymal transition and esg. This gene regulates stemness, ploidy, and pluripontency. esg is expressed in various tissues during development, including the gut, imaginal discs, and neuroblasts. The functions of the esg include the suppression of differentiation in intestinal stem cells and the preservation of diploidy in imaginal cells. In the nervous system development, esg expression also inhibits neuroblast differentiation, thus regulating the number of neurons and the moment in development of neuronal differentiation. Loss of esg function results in diverse developmental defects, including defects in intestinal stem cell maintenance and differentiation, and alters imaginal disc and nervous system development. Expression levels of esg also play a role in regulating longevity and metabolism in adult stages. This review provides an overview of the current understanding of esg's developmental role, emphasizing cellular and tissue effects that arise from its loss of function. The insights gained may contribute to a better understanding of evolutionary conserved developmental mechanisms and certain metabolic diseases.

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.

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.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni.

BACKGROUND: Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs). RESULTS: De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages. CONCLUSIONS: The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

A conserved gastropod withdrawal circuit in Biomphalaria glabrata, an intermediate host for schistosomiasis.

Neuronal signals mediated by the biogenic amine serotonin (5-HT) underlie critical survival strategies across the animal kingdom. This investigation examined serotonin-like immunoreactive neurons in the cerebral ganglion of the panpulmonate snail Biomphalaria glabrata, a major intermediate host for the trematode parasite Schistosoma mansoni. Five neurons comprising the cerebral serotonergic F (CeSF) cluster of B. glabrata shared morphological characteristics with neurons that contribute to withdrawal behaviors in numerous heterobranch species. The largest member of this group, designated CeSF-1, projected an axon to the tentacle, a major site of threat detection. Intracellular recordings demonstrated repetitive activity and electrical coupling between the bilateral CeSF-1 cells. In semi-intact preparations, the CeSF-1 cells were not responsive to cutaneous stimuli but did respond to photic stimuli. A large FMRF-NH2-like immunoreactive neuron, termed C2, was also located on the dorsal surface of each cerebral hemiganglion near the origin of the tentacular nerve. C2 and CeSF-1 received coincident bouts of inhibitory synaptic input. Moreover, in the presence of 5-HT they both fired rhythmically and in phase. As the CeSF and C2 cells of Biomphalaria share fundamental properties with neurons that participate in withdrawal responses in Nudipleura and Euopisthobranchia, our observations support the proposal that features of this circuit are conserved in the Panpulmonata.NEW & NOTEWORTHY Neuronal signals mediated by the biogenic amine serotonin underlie critical survival strategies across the animal kingdom. This investigation identified a group of serotonergic cells in the panpulmonate snail Biomphalaria glabrata that appear to be homologous to neurons that mediate withdrawal responses in other gastropod taxa. It is proposed that an ancient withdrawal circuit has been highly conserved in three major gastropod lineages.

Accelerated mitochondrial evolution and asymmetric fitness of hybrids contribute to the persistence of Helix thessalica in the Helix pomatia range.

Interbreeding and introgression between recently diverged species is common. However, the processes that prevent these species from merging where they co-occur are not well understood. We studied the mechanisms that allowed an isolated group of populations of the snail Helix thessalica to persist within the range of the related Helix pomatia despite high gene flow. Using genomic cline analysis, we found that the nuclear gene flow between the two taxa across the mosaic hybrid zone was not different from that expected under neutral admixture, but that the exchange of mtDNA was asymmetric. Tests showed that there is relaxed selection in the mitochondrial genome of H. thessalica and that the substitution rate is elevated compared to that of H. pomatia. A lack of hybrids that combine the mtDNA of H. thessalica with a mainly (>46%) H. pomatia genomic background indicates that the nuclear-encoded mitochondrial proteins of H. pomatia are not well adapted to the more rapidly evolving proteins and RNAs encoded by the mitochondrion of H. thessalica. The presumed reduction of fitness of hybrids with the fast-evolving mtDNA of H. thessalica and a high H. pomatia ancestry, similar to 'Darwin's Corollary to Haldane's rule', resulted in a relative loss of H. pomatia nuclear ancestry compared to H. thessalica ancestry in the hybrid zone. This probably prevents the H. thessalica populations from merging quickly with the surrounding H. pomatia populations and supports the hypothesis that incompatibilities between rapidly evolving mitochondrial genes and nuclear genes contribute to speciation.

Fatty alcohols, a minor component of the tree tobacco surface wax, are associated with defence against caterpillar herbivory.

Despite decades of research resulting in a comprehensive understanding of epicuticular wax metabolism, the function of these almost ubiquitous metabolites in plant-herbivore interactions remains unresolved. In this study, we examined the effects of CRISPR-induced knockout mutations in four Nicotiana glauca (tree tobacco) wax metabolism genes. These mutations cause a wide range of changes in epicuticular wax composition, leading to altered interactions with insects and snails. Three interaction classes were examined: chewing herbivory by seven caterpillars and one snail species, phloem feeding by Myzus persicae (green peach aphid) and oviposition by Bemisia tabaci (whitefly). Although total wax load and alkane abundance did not affect caterpillar growth, a correlation across species, showed that fatty alcohols, a minor component of N. glauca surface waxes, negatively affected the growth of both a generalist caterpillar (Spodoptera littoralis) and a tobacco-feeding specialist (Manduca sexta). This negative correlation was overshadowed by the stronger effect of anabasine, a nicotine isomer, and was apparent when fatty alcohols were added to an artificial lepidopteran diet. By contrast, snails fed more on waxy leaves. Aphid reproduction and feeding activity were unaffected by wax composition but were potentially affected by altered cutin composition. Wax crystal morphology could explain the preference of B. tabaci to lay eggs on waxy wild-type plants relative to both alkane and fatty alcohol-deficient mutants. Together, our results suggest that the varied responses among herbivore classes and species are likely to be a consequence of the co-evolution that shaped the specific effects of different surface wax components in plant-herbivore interactions.

Climate change could fuel urinary schistosomiasis transmission in Africa and Europe.

The freshwater snail Bulinus truncatus is an important intermediate host for trematode parasites causing urogenital schistosomiasis, a tropical disease affecting over 150 million people. Despite its medical importance, uncertainty remains about its global distribution and the potential impacts of climate change on its future spread. Here, we investigate the distribution of B. truncatus, combining the outputs of correlative and mechanistic modelling methods to fully capitalize on both experimental and occurrence data of the species and to create a more reliable distribution forecast than ever constructed. We constructed ensemble correlative species distribution models using 273 occurrence points collected from different sources and a combination of climatic and (bio)physical environmental variables. Additionally, a mechanistic thermal suitability model was constructed, parameterized by recent life-history data obtained through extensive lab-based snail-temperature experiments and supplemented with an extensive literature review. Our findings reveal that the current suitable habitat for B. truncatus encompasses the Sahel region, the Middle East, and the Mediterranean segment of Africa, stretching from Southern Europe to Mozambique. Regions identified as suitable by both methods generally coincide with areas exhibiting high urogenital schistosomiasis prevalence. Model projections into the future suggest an overall net increase in suitable area of up to 17%. New suitable habitat is in Southern Europe, the Middle East, and large parts of Central Africa, while suitable habitat will be lost in the Sahel region. The change in snail habitat suitability may substantially increase the risk of urogenital schistosomiasis transmission in parts of Africa and Southern Europe while reducing it in the Sahel region.

Continuous variation in the shell colour of the snail Cepaea nemoralis is associated with the colour locus of the supergene.

While the shell of the land snail Cepaea nemoralis is typically classed as yellow, pink, or brown, the reality is that colour variation is continuously distributed. To further understand the origin of the continuous variation, we used crosses of C. nemoralis to compare quantitative measures of the colour with the inferred genotype of the underlying supergene locus. We also used a recently developed linkage map to find quantitative trait loci that may influence colour. The results show that the colour locus of the supergene-at around 31.385 cM on linkage group 11-is involved in determining the quantitative chromatic differences that are perceptible to human vision. We also found some evidence that variation within colour classes may be due to allelic variation at or around the supergene. There are likely other unlinked loci involved in determining colour within classes, but confirmation will require greater statistical power. Although not investigated here, environmental factors, including diet, may also impact upon variation within colour types.

Function of snail shell hairs in anti-predator defense.

The function and evolutionary background of the hairs on the shells of terrestrial gastropods is largely unknown. Many hypotheses proposed by malacologists have never been proven, and the long-held hypothesis of mechanical stability in wet environments has been rejected by recent studies. It would therefore be worthwhile to reexamine other hypotheses regarding the adaptive significance of shell hairs. We investigated the defense function of shell hairs against a specialist predator, the snail-eating firefly, in the long-haired snail Moellendorffia diminuta. The firefly larvae, which hunt snails using abdominal suckers, were unable to attach to the shell because of the shell hairs but were able to attach to the shells that had lost their hairs. About half of the hairy snails successfully defended themselves by swinging their shells and dropping firefly larvae, but most of the snails without hair failed to defend. The hairs reduce the ability of the larva to attach to the shell and increase the effectiveness of the shell-swinging defense behavior in removing the larva from the shell. As shell hairs grow longer with shell development, they may confer an advantage based on the predator's growth stage. Our findings highlight the anti-predator defense role of shell hairs in land snails, introducing a hypothesis previously overlooked in the evolutionary context of hairy snails.