Post-feeding transcriptomics reveals essential genes expressed in the midgut of the desert locust.

The digestive tract constitutes an important interface between an animal's internal and external environment. In insects, available gut transcriptome studies are mostly exploratory or look at changes upon infection or upon exposure to xenobiotics, mainly performed in species belonging to holometabolan orders, such as Diptera, Lepidoptera or Coleoptera. By contrast, studies focusing on gene expression changes after food uptake and during digestion are underrepresented. We have therefore compared the gene expression profiles in the midgut of the desert locust, Schistocerca gregaria, between three different time points after feeding, i.e., 24 h (no active digestion), 10 min (the initial stage of feeding), and 2 h (active food digestion). The observed gene expression profiles were consistent with the polyphagous herbivorous lifestyle of this hemimetabolan (orthopteran) species. Our study reveals the upregulation of 576 genes 2 h post-feeding. These are mostly predicted to be associated with digestive physiology, such as genes encoding putative digestive enzymes or nutrient transporters, as well as genes putatively involved in immunity or in xenobiotic metabolism. The 10 min time point represented an intermediate condition, suggesting that the S. gregaria midgut can react rapidly at the transcriptional level to the presence of food. Additionally, our study demonstrated the critical importance of two transcripts that exhibited a significant upregulation 2 h post-feeding: the vacuolar-type H(+)-ATPase and the sterol transporter Niemann-Pick 1b protein, which upon RNAi-induced knockdown resulted in a marked increase in mortality. Their vital role and accessibility via the midgut lumen may make the encoded proteins promising insecticidal target candidates, considering that the desert locust is infamous for its huge migrating swarms that can devastate the agricultural production in large areas of Northern Africa, the Middle East, and South Asia. In conclusion, the transcriptome datasets presented here will provide a useful and promising resource for studying the midgut physiology of S. gregaria, a socio-economically important pest species.

Sexual repurposing of juvenile aposematism in locusts.

Adaptive plasticity requires an integrated suite of functional responses to environmental variation, which can include social communication across life stages. Desert locusts (Schistocerca gregaria) exhibit an extreme example of phenotypic plasticity called phase polyphenism, in which a suite of behavioral and morphological traits differ according to local population density. Male and female juveniles developing at low population densities exhibit green- or sand-colored background-matching camouflage, while at high densities they show contrasting yellow and black aposematic patterning that deters predators. The predominant background colors of these phenotypes (green/sand/yellow) all depend on expression of the carotenoid-binding "Yellow Protein" (YP). Gregarious (high-density) adults of both sexes are initially pinkish, before a YP-mediated yellowing reoccurs upon sexual maturation. Yellow color is especially prominent in gregarious males, but the reason for this difference has been unknown since phase polyphenism was first described in 1921. Here, we use RNA interference to show that gregarious male yellowing acts as an intrasexual warning signal, which forms a multimodal signal with the antiaphrodisiac pheromone phenylacetonitrile (PAN) to prevent mistaken sexual harassment from other males during scramble mating in a swarm. Socially mediated reexpression of YP thus adaptively repurposes a juvenile signal that deters predators into an adult signal that deters undesirable mates. These findings reveal a previously underappreciated sexual dimension to locust phase polyphenism, and promote locusts as a model for investigating the relative contributions of natural versus sexual selection in the evolution of phenotypic plasticity.

Prothoracicostatic Activity of the Ecdysis-Regulating Neuropeptide Crustacean Cardioactive Peptide (CCAP) in the Desert Locust.

Accurate control of innate behaviors associated with developmental transitions requires functional integration of hormonal and neural signals. Insect molting is regulated by a set of neuropeptides, which trigger periodic pulses in ecdysteroid hormone titers and coordinate shedding of the old cuticle during ecdysis. In the current study, we demonstrate that crustacean cardioactive peptide (CCAP), a structurally conserved neuropeptide described to induce the ecdysis motor program, also exhibits a previously unknown prothoracicostatic activity to regulate ecdysteroid production in the desert locust, Schistocerca gregaria. We identified the locust genes encoding the CCAP precursor and three G protein-coupled receptors that are activated by CCAP with EC50 values in the (sub)nanomolar range. Spatiotemporal expression profiles of the receptors revealed expression in the prothoracic glands, the endocrine organs where ecdysteroidogenesis occurs. RNAi-mediated knockdown of CCAP precursor or receptors resulted in significantly elevated transcript levels of several Halloween genes, which encode ecdysteroid biosynthesis enzymes, and in elevated ecdysteroid levels one day prior to ecdysis. Moreover, prothoracic gland explants exhibited decreased secretion of ecdysteroids in the presence of CCAP. Our results unequivocally identify CCAP as the first prothoracicostatic peptide discovered in a hemimetabolan species and reveal the existence of an intricate interplay between CCAP signaling and ecdysteroidogenesis.

Role of peptide hormones in insect gut physiology.

Nutrient uptake and digestion are essential for optimal growth and development. In insects, these processes are regulated by the gut-brain axis, which is a neurohumoral communication system for maintaining gut homeostasis. The insect gut is a complex organ consisting of three distinct structures, denominated foregut, midgut and hindgut, each with their specific specializations. These specializations are tightly regulated by the interplay of several neuropeptides: a versatile group of signalling molecules involved in a multitude of processes including gut physiology. Neuropeptides take part in the regulation of gut processes ranging from digestive enzyme release to muscle activity and satiety. Some neuropeptide mimetics are a promising strategy for ecological pest management. This review focuses on a selection of neuropeptides that are well-known for their role in gut physiology, and neuropeptides for which the mode of action is yet to be unravelled.

Can BRET-based biosensors be used to characterize G-protein mediated signaling pathways of an insect GPCR, the Schistocerca gregaria CRF-related diuretic hormone receptor?

G protein-coupled receptors (GPCRs) are membrane-bound receptors that are considered prime candidates for the development of novel insect pest management strategies. However, the molecular signaling properties of insect GPCRs remain poorly understood. In fact, most studies on insect GPCR signaling are limited to analysis of fluctuations in the secondary messenger molecules calcium (Ca2+) and/or cyclic adenosine monophosphate (cAMP). In the current study, we characterized a corticotropin-releasing factor-related diuretic hormone (CRF-DH) receptor of the desert locust, Schistocerca gregaria. This Schgr-CRF-DHR is mainly expressed in the nervous system and in brain-associated endocrine organs. The neuropeptide Schgr-CRF-DH induced Ca2+-dependent aequorin-based bioluminescent responses in CHO cells co-expressing this receptor with the promiscuous Gα16 protein. Furthermore, when co-expressed with the cAMP-dependent bioluminescence resonance energy transfer (BRET)-based CAMYEL biosensor in HEK293T cells, this receptor elicited dose-dependent agonist-induced responses with an EC50 in the nanomolar range (4.02 nM). In addition, we tested if vertebrate BRET-based G protein biosensors, can also be used to detect direct Gα protein subunit activation by an insect GPCR. Therefore, we analyzed ten different human BRET-based G protein biosensors, representing members of all four Gα protein subfamilies; Gαs, Gαi/o, Gαq/11 and Gα12/13. Our data demonstrate that stimulation of Schgr-CRF-DHR by Schgr-CRF-DH can dose-dependently activate Gαi/o and Gαs biosensors, while no significant effects were observed with the Gαq/11 and Gα12/13 biosensors. Our study paves the way for future biosensor-based studies to analyze the signaling properties of insect GPCRs in both fundamental science and applied research contexts.

First draft genome assembly of the desert locust, Schistocerca gregaria.

Background: At the time of publication, the most devastating desert locust crisis in decades is affecting East Africa, the Arabian Peninsula and South-West Asia. The situation is extremely alarming in East Africa, where Kenya, Ethiopia and Somalia face an unprecedented threat to food security and livelihoods. Most of the time, however, locusts do not occur in swarms, but live as relatively harmless solitary insects. The phenotypically distinct solitarious and gregarious locust phases differ markedly in many aspects of behaviour, physiology and morphology, making them an excellent model to study how environmental factors shape behaviour and development. A better understanding of the extreme phenotypic plasticity in desert locusts will offer new, more environmentally sustainable ways of fighting devastating swarms. Methods: High molecular weight DNA derived from two adult males was used for Mate Pair and Paired End Illumina sequencing and PacBio sequencing. A reliable reference genome of Schistocerca gregaria was assembled using the ABySS pipeline, scaffolding was improved using LINKS. Results: In total, 1,316 Gb Illumina reads and 112 Gb PacBio reads were produced and assembled. The resulting draft genome consists of 8,817,834,205 bp organised in 955,015 scaffolds with an N50 of 157,705 bp, making the desert locust genome the largest insect genome sequenced and assembled to date. In total, 18,815 protein-encoding genes are predicted in the desert locust genome, of which 13,646 (72.53%) obtained at least one functional assignment based on similarity to known proteins. Conclusions: The desert locust genome data will contribute greatly to studies of phenotypic plasticity, physiology, neurobiology, molecular ecology, evolutionary genetics and comparative genomics, and will promote the desert locust's use as a model system. The data will also facilitate the development of novel, more sustainable strategies for preventing or combating swarms of these infamous insects.

Peptides in insect oogenesis.

The physiological control of reproduction in insects depends on a combination of environmental and internal cues. In the adult stage, insects become sexually mature and generate gametes. In females, the latter process is designated as oogenesis. Peptides are a versatile class of extracellular signalling molecules that regulate many processes, including oogenesis. At present, the best documented physiological control mechanism of insect oogenesis is the insulin-related peptide signalling pathway. It regulates different stages of the process and provides a functional link between nutritional status and reproduction. Several other peptides have been shown to exert gonadoregulatory activities, but in most cases their exact mode of action still has to be unravelled and their effects on oogenesis could be direct or indirect. Some regulatory peptides, such as the Drosophila sex peptide, are being transferred from the male to the female during the mating process.

Pleiotropic and novel phenotypes in the Drosophila gut caused by mutation of drop-dead.

Normal gut function is vital for animal survival, and deviations from such function can contribute to malnutrition, inflammation, increased susceptibility to pathogens, diabetes, neurodegenerative diseases, and cancer. In the fruit fly Drosophila melanogaster, mutation of the gene drop-dead (drd) results in defective gut function, as measured by enlargement of the crop and reduced food movement through the gut, and drd mutation also causes the unrelated phenotypes of neurodegeneration, early adult lethality and female sterility. In the current work, adult drd mutant flies are also shown to lack the peritrophic matrix (PM), an extracellular barrier that lines the lumen of the midgut and is found in many insects including flies, mosquitos and termites. The use of a drd-gal4 construct to drive a GFP reporter in late pupae and adults revealed drd expression in the anterior cardia, which is the site of PM synthesis in Drosophila. Moreover, the ability of drd knockdown or rescue with several gal4 drivers to recapitulate or rescue the gut phenotypes (lack of a PM, reduced defecation, and reduced adult survival 10-40 days post-eclosion) was correlated to the level of expression of each driver in the anterior cardia. Surprisingly, however, knocking down drd expression only in adult flies, which has previously been shown not to affect survival, eliminated the PM without reducing defecation rate. These results demonstrate that drd mutant flies have a novel phenotype, the absence of a PM, which is functionally separable from the previously described gut dysfunction observed in these flies. As the first mutant Drosophila strain reported to lack a PM, drd mutants will be a useful tool for studying the synthesis of this structure.

Molecular cloning and characterization of the SIFamide precursor and receptor in a hymenopteran insect, Bombus terrestris.

SIFamides (SIFa) are a family of neuropeptides that are highly conserved among arthropods. In insects, this peptide is mainly expressed in four medial interneurons in the pars intercerebralis and affects sexual behavior, sleep regulation and pupal mortality. Furthermore, an influence on the hatching rate has been observed. The first SIFa receptor (SIFR) was pharmacologically characterized in Drosophila melanogaster and is homologous to the vertebrate gonadotropin-inhibitory hormone (GnIH) receptor (NPFFR). In this study, we pharmacologically characterized the SIFR of the buff-tailed bumblebee Bombus terrestris. We demonstrated an intracellular increase in calcium ions and cyclic AMP (cAMP) upon ligand binding with an EC50 value in the picomolar and nanomolar range, respectively. In addition, we studied the agonistic properties of a range of related and modified peptides. By means of quantitative real time PCR (qPCR), we examined the relative transcript levels of Bomte-SIFa and Bomte-SIFR in a variety of tissues.

Molecular cloning and characterization of the allatotropin precursor and receptor in the desert locust, Schistocerca gregaria.

Allatotropins (ATs) are pleiotropic neuropeptides initially isolated from the tobacco hornworm, Manduca sexta. In 2008, the first receptor for AT-like peptides (ATR) was characterized in Bombyx mori. Since then, ATRs have also been characterized in M. sexta, Tribolium castaneum, Aedes aegypti and Bombus terrestris. These receptors show sequence similarity to vertebrate orexin (ORX) receptors. When generating an EST-database of the desert locust (Schistocerca gregaria) central nervous system, we found cDNA sequences encoding the Schgr-AT precursor and a fragment of its putative receptor. This receptor cDNA has now been completed and functionally expressed in mammalian cell lines. Activation of this receptor, designated as Schgr-ATR, by Schgr-AT caused an increase in intracellular calcium ions, as well as cyclic AMP (cAMP), with an EC50 value in the nanomolar range. In addition, the transcript distribution of both the Schgr-AT precursor and Schgr-ATR was investigated by means of quantitative real-time PCR. Moreover, we found more evidence for the myotropic and allatostimulatory actions of Schgr-AT in the desert locust. These data are discussed and situated in a broader context by comparison with literature data on AT and ATR in insects.

Pharmacological and signalling properties of a D2-like dopamine receptor (Dop3) in Tribolium castaneum.

Dopamine is an important neurotransmitter in the central nervous system of vertebrates and invertebrates. Despite their evolutionary distance, striking parallels exist between deuterostomian and protostomian dopaminergic systems. In both, signalling is achieved via a complement of functionally distinct dopamine receptors. In this study, we investigated the sequence, pharmacology and tissue distribution of a D2-like dopamine receptor from the red flour beetle Tribolium castaneum (TricaDop3) and compared it with related G protein-coupled receptors in other invertebrate species. The TricaDop3 receptor-encoding cDNA shows considerable sequence similarity with members of the Dop3 receptor class. Real time qRT-PCR showed high expression in both the central brain and the optic lobes, consistent with the role of dopamine as neurotransmitter. Activation of TricaDop3 expressed in mammalian cells increased intracellular Ca(2+) signalling and decreased NKH-477 (a forskolin analogue)-stimulated cyclic AMP levels in a dose-dependent manner. We studied the pharmacological profile of the TricaDop3 receptor and demonstrated that the synthetic vertebrate dopamine receptor agonists, 2 - amino- 6,7 - dihydroxy - 1,2,3,4 - tetrahydronaphthalene hydrobromide (6,7-ADTN) and bromocriptine acted as agonists. Methysergide was the most potent of the antagonists tested and showed competitive inhibition in the presence of dopamine. This study offers important information on the Dop3 receptor from Tribolium castaneum that will facilitate functional analyses of dopamine receptors in insects and other invertebrates.

Mode of action of allatostatins in the regulation of juvenile hormone biosynthesis in the cockroach, Diploptera punctata.

The FGLamide allatostatins (FGL/ASTs) are a family of neuropeptides with pleiotropic functions, including the inhibition of juvenile hormone (JH) biosynthesis, vitellogenesis and muscle contraction. In the cockroach, Diploptera punctata, thirteen FGLa/ASTs and one allatostatin receptor (AstR) have been identified. However, the mode of action of ASTs in regulation of JH biosynthesis remains unclear. Here, we determined the tissue distribution of Dippu-AstR. And we expressed Dippu-AstR in vertebrate cell lines, and activated the receptor with the Dippu-ASTs. Our results show that all thirteen ASTs activated Dippu-AstR in a dose dependent manner, albeit with different potencies. Functional analysis of AstR in multiple cell lines demonstrated that activation of the AstR receptor resulted in elevated levels of Ca(2+) and cAMP, which suggests that Dippu-AstR can act through the Gαq and Gαs protein pathways. The study on the target of AST action reveals that FGL/AST affects JH biosynthesis prior to the entry of acetyl-CoA into the JH biosynthetic pathway.

Lipophorins can adhere to dsRNA, bacteria and fungi present in the hemolymph of the desert locust: a role as general scavenger for pathogens in the open body cavity.

Desert locusts are characterized by a highly sensitive and effective RNA interference (RNAi) response. Moreover, delivery of dsRNA into the open body cavity will elicit potent silencing effects throughout the body. On the other hand, many other insect species, such as Bombyx mori and Drosophila melanogaster, lack the ability to efficiently spread the RNAi-signal. In this study, we demonstrated that, in the serum of the desert locust, lipophorins adhere to dsRNA-fragments. Lipophorins can be subdivided into high density and low density lipophorins (HDLp and LDLp), according to their buoyant density, and we showed that both types of lipophorins can interact with dsRNA fragments. Furthermore, in the presence of (gram-positive) bacteria or fungi, LDLp induce the formation of pathogen aggregates, while no clear aggregation effects were detected in the presence of HDLp.

Characterisation of a functional allatotropin receptor in the bumblebee, Bombus terrestris (Hymenoptera, Apidae).

Allatotropins (ATs) are multifunctional neuropeptides initially isolated from the tobacco hornworm, Manduca sexta, where they were found to stimulate juvenile hormone synthesis and release from the corpora allata. ATs have been found in a wide range of insects, but appear to be absent in Drosophila. The first AT receptor (ATR) was characterised in 2008 in the lepidopteran Bombyx mori. Since then ATRs have been characterised in Coleoptera and Diptera and in 2012, an AT precursor gene was identified in hymenopteran species. ATRs show large sequence and structural similarity to vertebrate orexin receptors (OXR). Also, AT in insects and orexin in vertebrates show some overlap in functions, including modulation of feeding behaviour and reproduction. The goal of this study was to identify a functional ATR in a hymenopteran species. We used ATRs (insect sequences) and OXRs (vertebrate sequences) to search the genome of the bumblebee, Bombus terrestris. Two receptors (XP_003402490 and XP_003394933) with resemblance to ATRs and OXRs were found. Phylogenetic analysis provided the first indication that XP_003402490 was more closely related to ATRs than XP_003394933. We investigated the transcript level distribution of both receptors and the AT precursor gene by means of quantitative real-time reverse transcriptase PCR. XP_003402490 displayed a tissue distribution comparable with ATRs in other species, with high transcript levels in the male accessory glands. After pharmacological characterisation, it appeared that XP_003402490 is indeed a functional ATR. Activation of the receptor causes an increase in intracellular calcium and cyclic AMP levels with an EC50 value in the low nanomolar to picomolar range. XP_003394933 remains an orphan receptor.

Pharmacological characterization of a 5-HT1-type serotonin receptor in the red flour beetle, Tribolium castaneum.

Serotonin (5-hydroxytryptamine, 5-HT) is known for its key role in modulating diverse physiological processes and behaviors by binding various 5-HT receptors. However, a lack of pharmacological knowledge impedes studies on invertebrate 5-HT receptors. Moreover, pharmacological information is urgently needed in order to establish a reliable classification system for invertebrate 5-HT receptors. In this study we report on the molecular cloning and pharmacological characterization of a 5-HT1 receptor from the red flour beetle, Tribolium castaneum (Trica5-HT1). The Trica5-HT1 receptor encoding cDNA shows considerable sequence similarity with members of the 5-HT1 receptor class. Real time PCR showed high expression in the brain (without optic lobes) and the optic lobes, consistent with the role of 5-HT as neurotransmitter. Activation of Trica5-HT1 in mammalian cells decreased NKH-477-stimulated cyclic AMP levels in a dose-dependent manner, but did not influence intracellular Ca(2+) signaling. We studied the pharmacological profile of the 5-HT1 receptor and demonstrated that α-methylserotonin, 5-methoxytryptamine and 5-carboxamidotryptamine acted as agonists. Prazosin, methiothepin and methysergide were the most potent antagonists and showed competitive inhibition in presence of 5-HT. This study offers important information on a 5-HT1 receptor from T. castaneum facilitating functional research of 5-HT receptors in insects and other invertebrates. The pharmacological profiles may contribute to establish a reliable classification scheme for invertebrate 5-HT receptors.

Sex peptides and MIPs can activate the same G protein-coupled receptor.

In many animal species, copulation elicits a number of physiological and behavioral changes in the female partner. In Drosophila melanogaster, the main molecular effector of these physiological responses has been identified as sex peptide (SP). The sex peptide receptor (SPR) has been characterized and recently, its activation by Drosophila myoinhibiting peptides (MIPs)-in addition to SP-has been demonstrated. The myoinhibiting peptides are members of a conserved peptide family, also known as B-type allatostatins, which generally feature the C-terminal motif -WX6Wamide.

Phylogenetic distribution of protease inhibitors of the Kazal-family within the Arthropoda.

In mammalian pancreatic cells, the pancreatic secretory trypsin inhibitor (PSTI) belonging to the Kazal-family prevents the premature activation of digestive enzymes and thus plays an important role in a protective mechanism against tissue destruction by autophagy. Although a similar protective mechanism exists in Arthropoda, the distribution of these inhibitors in this phylum remains obscure. A comprehensive in silico search of nucleotide databases, revealed the presence of members of the Kazal-family in the four major subphyla of the Arthropoda. Especially in the Hexapoda and the Crustacea these inhibitors are widespread, while in the Chelicerata and Myriapoda only a few Kazal-like protease inhibitors were found. A sequence alignment of inhibitors retrieved in the digestive system of insects revealed a conservation of the PSTI characteristics and strong resemblance to vertebrate PSTI. A phylogenetic analysis of these inhibitors showed that they generally cluster according to their order. The results of this data mining study provide new evidence for the existence of an ancient protective mechanism in metazoan digestive systems. Kazal-like inhibitors, which play an important protective role in the pancreas of vertebrates, also seem to be present in Arthropoda.

Critical role for protein kinase A in the acquisition of gregarious behavior in the desert locust.

The mechanisms that integrate genetic and environmental information to coordinate the expression of complex phenotypes are little understood. We investigated the role of two protein kinases (PKs) in the population density-dependent transition to gregarious behavior that underlies swarm formation in desert locusts: the foraging gene product, a cGMP-dependent PK (PKG) implicated in switching between alternative group-related behaviors in several animal species; and cAMP-dependent PK (PKA), a signal transduction protein with a preeminent role in different forms of learning. Solitarious locusts acquire key behavioral characters of the swarming gregarious phase within just 1 to 4 h of forced crowding. Injecting the PKA inhibitor KT5720 before crowding prevented this transition, whereas injecting KT5823, an inhibitor of PKG, did not. Neither drug altered the behavior of long-term gregarious locusts. RNAi against foraging effectively reduced its expression in the central nervous system, but this did not prevent gregarization upon crowding. By contrast, solitarious locusts with an RNAi-induced reduction in PKA catalytic subunit C1 expression behaved less gregariously after crowding, and RNAi against the inhibitory R1 subunit promoted more extensive gregarization following a brief crowding period. A central role of PKA is congruent with the recent discovery that serotonin mediates gregarization in locusts and with findings in vertebrates that similarly implicate PKA in the capacity to cope with adverse life events. Our results show that PKA has been coopted into effecting the wide-ranging transformation from solitarious to gregarious behavior, with PKA-mediated behavioral plasticity resulting in an environmentally driven reorganization of a complex phenotype.

Functional analysis of a pancreatic secretory trypsin inhibitor-like protein in insects: silencing effects resemble the human pancreatic autodigestion phenotype.

INTRODUCTION: In mammalian pancreatic cells, the pancreatic secretory trypsin inhibitor (PSTI) prevents the premature activation of digestive enzymes and thus plays an important role in a protective mechanism against tissue destruction by autophagy, a process which may ultimately cause diseases such as pancreatitis and pancreatic cancer. Insects, however, lack a pancreas and so far no PSTI-like peptides are functionally characterized. RESULTS: In several insect species protease inhibitors that structurally resemble the mammalian PSTI were predicted in silico. A putative PSTI-like protein (LmPSTI) was cloned and sequenced in the African migratory locust, Locusta migratoria. For the first time the expression of an insect derived PSTI-like inhibitor was shown to be restricted to the digestive enzyme-producing organs in insects (midgut and caeca). LmPSTI was produced via a bacterial expression system and was found to be a potent inhibitor of bovine trypsin as well as endogenous locust gut enzymes. In the caeca, RNAi-mediated knockdown of LmPSTI resulted in a significantly upregulated expression (2-fold) of locust ATG8 transcripts (an ubiquitin-like protein crucial for autophagosome formation). These findings were confirmed by an ultrastructural study on caeca, revealing the presence of autophagy-related structures in RNAi-treated animals. CONCLUSION: The results of this study lead us to believe that LmPSTI plays an important role in controlling the proteolytic activity in the digestive system of L. migratoria. These findings provide new evidence for the existence of an ancient protective mechanism in metazoan digestive systems and open new perspectives for the study of autophagy-related diseases in the digestive tract.