Parallel duplication and partial subfunctionalization of β-catenin/armadillo during insect evolution

Global Gene Expression Analysis Reveals Complex Cuticle Organization of the Tribolium Compound Eye

The red flour beetle Tribolium castaneum is a resource-rich model for genomic and developmental studies. To extend previous studies on Tribolium eye development, we produced transcriptomes for normal-eyed and eye-depleted heads of pupae and adults to identify differentially transcript-enriched (DE) genes in the visual system. Unexpectedly, cuticle-related genes were the largest functional class in the pupal compound eye DE gene population, indicating differential enrichment in three distinct cuticle components: clear lens facet cuticle, highly melanized cuticle of the ocular diaphragm, which surrounds the Tribolium compound eye for internal fortification, and newly identified facet margins of the tanned cuticle, possibly enhancing external fortification. Phylogenetic, linkage, and high-throughput gene knockdown data suggest that most cuticle proteins (CPs) expressed in the Tribolium compound eye stem from the deployment of ancient CP genes. Consistent with this, TcasCPR15, which we identified as the major lens CP gene in Tribolium, is a beetle-specific but pleiotropic paralog of the ancient CPR RR-2 CP gene family. The less abundant yet most likely even more lens-specific TcasCP63 is a member of a sprawling family of noncanonical CP genes, documenting a role of local gene family expansions in the emergence of the Tribolium compound eye CP repertoire. Comparisons with Drosophila and the mosquito Anopheles gambiae reveal a steady turnover of lens-enriched CP genes during insect evolution.

Utilizing a chromosomal-length genome assembly to annotate the Wnt signaling pathway in the Asian citrus psyllid, Diaphorina citri

The Asian citrus psyllid, Diaphorina citri, is an insect vector that transmits Candidatus Liberibacter asiaticus, the causal agent of the Huanglongbing (HLB), or citrus greening disease. This disease has devastated Florida's citrus industry, and threatens California's industry as well as other citrus producing regions around the world. To find novel solutions to the disease, a better understanding of the vector is needed. The D. citri genome has been used to identify and characterize genes involved in Wnt signaling pathways. Wnt signaling is utilized for many important biological processes in metazoans, such as patterning and tissue generation. Curation based on RNA sequencing data and sequence homology confirms 24 Wnt signaling genes within the D. citri genome, including homologs for beta-catenin, Frizzled receptors, and seven Wnt-ligands. Through phylogenetic analysis, we classify D. citri Wnt ligands as Wg/Wnt1, Wnt5, Wnt6, Wnt7, Wnt10, Wnt11, and WntA. The D. citri version 3.0 genome with chromosomal length scaffolds reveals a conserved Wnt1-Wnt6-Wnt10 gene cluster with a gene configuration like that in Drosophila melanogaster. These findings provide greater insight into the evolutionary history of D. citri and Wnt signaling in this important hemipteran vector. Manual annotation was essential for identifying high quality gene models. These gene models can be used to develop molecular systems, such as CRISPR and RNAi, which target and control psyllid populations to manage the spread of HLB. Manual annotation of Wnt signaling pathways was done as part of a collaborative community annotation project.

Recurrent sequence evolution after independent gene duplication

Background

Convergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa.

Results

We develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3.

Conclusions

The presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

New genomic data and analyses challenge the traditional vision of animal epithelium evolution

Background

The emergence of epithelia was the foundation of metazoan expansion. Epithelial tissues are a hallmark of metazoans deeply rooted in the evolution of their complex developmental morphogenesis processes. However, studies on the epithelial features of non-bilaterians are still sparse and it remains unclear whether the last common metazoan ancestor possessed a fully functional epithelial toolkit or if it was acquired later during metazoan evolution.

Results

To investigate the early evolution of animal epithelia, we sequenced the genome and transcriptomes of two new sponge species to characterize epithelial markers such as the E-cadherin complex and the polarity complexes for all classes (Calcarea, Demospongiae, Hexactinellida, Homoscleromorpha) of sponges (phylum Porifera) and compare them with their homologues in Placozoa and in Ctenophora. We found that Placozoa and most sponges possess orthologues of all essential genes encoding proteins characteristic of bilaterian epithelial cells, as well as their conserved interaction domains. In stark contrast, we found that ctenophores lack several major polarity complex components such as the Crumbs complex and Scribble. Furthermore, the E-cadherin ctenophore orthologue exhibits a divergent cytoplasmic domain making it unlikely to interact with its canonical cytoplasmic partners.

Conclusions

These unexpected findings challenge the current evolutionary paradigm on the emergence of epithelia. Altogether, our results raise doubt on the homology of protein complexes and structures involved in cell polarity and adhesive-type junctions between Ctenophora and Bilateria epithelia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4715-9) contains supplementary material, which is available to authorized users.

A C-terminally truncated form of β-catenin acts as a novel regulator of Wnt/β-catenin signaling in planarians

β-Catenin, the core element of the Wnt/β-catenin pathway, is a multifunctional and evolutionarily conserved protein which performs essential roles in a variety of developmental and homeostatic processes. Despite its crucial roles, the mechanisms that control its context-specific functions in time and space remain largely unknown. The Wnt/β-catenin pathway has been extensively studied in planarians, flatworms with the ability to regenerate and remodel the whole body, providing a ‘whole animal’ developmental framework to approach this question. Here we identify a C-terminally truncated β-catenin (β-catenin4), generated by gene duplication, that is required for planarian photoreceptor cell specification. Our results indicate that the role of β-catenin4 is to modulate the activity of β-catenin1, the planarian β-catenin involved in Wnt signal transduction in the nucleus, mediated by the transcription factor TCF-2. This inhibitory form of β-catenin, expressed in specific cell types, would provide a novel mechanism to modulate nuclear β-catenin signaling levels. Genomic searches and in vitro analysis suggest that the existence of a C-terminally truncated form of β-catenin could be an evolutionarily conserved mechanism to achieve a fine-tuned regulation of Wnt/β-catenin signaling in specific cellular contexts.

The Wnt signaling pathway is essential for proper intercellular communication in every developmental process since it controls basic cellular events as cell fate or proliferation. The key element of the Wnt signaling is β-catenin, which controls the transcription of multiple genes in the Wnt receiving cell. A main level of regulation of the Wnt/β-catenin signaling occurs in the cytoplasm, where β-catenin protein levels depend on the activity of the β-catenin destruction complex. However, once it reaches the nucleus, β-catenin transcriptional activity requires a fine-tuned regulation to enable the multiple context-specific responses that it performs. These nuclear mechanisms that regulate the Wnt/β-catenin signaling remain poorly understood. Here we report the existence of C-terminal truncated forms of β-catenin in planarians (β-cat3 and 4), which, in vitro, do not show transactivation activity and compete with the canonical planarian β-catenin (β-cat1), thus acting as competitor inhibitors. Functional analyses in planarians indicate that β-cat4 acts as a negative regulator of β-cat1 during planarian eye photoreceptor specification. We provide evidence to suggest that this novel mechanism for the regulation of nuclear β-catenin activity could be conserved across animal evolution.

Evolution of selenophosphate synthetases: emergence and relocation of function through independent duplications and recurrent subfunctionalization

Selenoproteins are proteins that incorporate selenocysteine (Sec), a nonstandard amino acid encoded by UGA, normally a stop codon. Sec synthesis requires the enzyme Selenophosphate synthetase (SPS or SelD), conserved in all prokaryotic and eukaryotic genomes encoding selenoproteins. Here, we study the evolutionary history of SPS genes, providing a map of selenoprotein function spanning the whole tree of life. SPS is itself a selenoprotein in many species, although functionally equivalent homologs that replace the Sec site with cysteine (Cys) are common. Many metazoans, however, possess SPS genes with substitutions other than Sec or Cys (collectively referred to as SPS1). Using complementation assays in fly mutants, we show that these genes share a common function, which appears to be distinct from the synthesis of selenophosphate carried out by the Sec- and Cys- SPS genes (termed SPS2), and unrelated to Sec synthesis. We show here that SPS1 genes originated through a number of independent gene duplications from an ancestral metazoan selenoprotein SPS2 gene that most likely already carried the SPS1 function. Thus, in SPS genes, parallel duplications and subsequent convergent subfunctionalization have resulted in the segregation to different loci of functions initially carried by a single gene. This evolutionary history constitutes a remarkable example of emergence and evolution of gene function, which we have been able to trace thanks to the singular features of SPS genes, wherein the amino acid at a single site determines unequivocally protein function and is intertwined to the evolutionary fate of the entire selenoproteome.

Faster evolving Drosophila paralogs lose expression rate and ubiquity and accumulate more non-synonymous SNPs

Background

Duplicated genes can indefinately persist in genomes if either both copies retain the original function due to dosage benefit (gene conservation), or one of the copies assumes a novel function (neofunctionalization), or both copies become required to perform the function previously accomplished by a single copy (subfunctionalization), or through a combination of these mechanisms. Different models of duplication retention imply different predictions about substitution rates in the coding portion of paralogs and about asymmetry of these rates.

Results

We analyse sequence evolution asymmetry in paralogs present in 12 Drosophila genomes using the nearest non-duplicated orthologous outgroup as a reference. Those paralogs present in D. melanogaster are analysed in conjunction with the asymmetry of expression rate and ubiquity and of segregating non-synonymous polymorphisms in the same paralogs. Paralogs accumulate substitutions, on average, faster than their nearest singleton orthologs. The distribution of paralogs’ substitution rate asymmetry is overdispersed relative to that of orthologous clades, containing disproportionally more unusually symmetric and unusually asymmetric clades. We show that paralogs are more asymmetric in: a) clades orthologous to highly constrained singleton genes; b) genes with high expression level; c) genes with ubiquitous expression and d) non-tandem duplications. We further demonstrate that, in each asymmetrically evolving pair of paralogs, the faster evolving member of the pair tends to have lower average expression rate, lower expression uniformity and higher frequency of non-synonymous SNPs than its slower evolving counterpart.

Conclusions

Our findings are consistent with the hypothesis that many duplications in Drosophila are retained despite stabilising selection being more relaxed in one of the paralogs than in the other, suggesting a widespread unfinished pseudogenization. This phenomenon is likely to make detection of neo- and subfunctionalization signatures difficult, as these models of duplication retention also predict asymmetries in substitution rates and expression profiles.

Reviewers

This article has been reviewed by Dr. Jia Zeng (nominated by Dr. I. King Jordan), Dr. Fyodor Kondrashov and Dr. Yuri Wolf.

High plasticity in epithelial morphogenesis during insect dorsal closure

Insect embryos complete the outer form of the body via dorsal closure (DC) of the epidermal flanks, replacing the transient extraembryonic (EE) tissue. Cell shape changes and morphogenetic behavior are well characterized for DC in Drosophila, but these data represent a single species with a secondarily reduced EE component (the amnioserosa) that is not representative across the insects. Here, we examine DC in the red flour beetle, Tribolium castaneum, providing the first detailed, functional analysis of DC in an insect with complete EE tissues (distinct amnion and serosa). Surprisingly, we find that differences between Drosophila and Tribolium DC are not restricted to the EE tissue, but also encompass the dorsal epidermis, which differs in cellular architecture and method of final closure (zippering). We then experimentally manipulated EE tissue complement via RNAi for Tc-zen1, allowing us to eliminate the serosa and still examine viable DC in a system with a single EE tissue (the amnion). We find that the EE domain is particularly plastic in morphogenetic behavior and tissue structure. In contrast, embryonic features and overall kinetics are robust to Tc-zen1(RNAi) manipulation in Tribolium and conserved with a more distantly related insect, but remain substantially different from Drosophila. Although correct DC is essential, plasticity and regulative, compensatory capacity have permitted DC to evolve within the insects. Thus, DC does not represent a strong developmental constraint on the nature of EE development, a property that may have contributed to the reduction of the EE component in the fly lineage.

The evolutionary origin of epithelial cell-cell adhesion mechanisms

A simple epithelium forms a barrier between the outside and the inside of an organism, and is the first organized multicellular tissue found in evolution. We examine the relationship between the evolution of epithelia and specialized cell-cell adhesion proteins comprising the classical cadherin/β-catenin/α-catenin complex (CCC). A review of the divergent functional properties of the CCC in metazoans and non-metazoans, and an updated phylogenetic coverage of the CCC using recent genomic data reveal: (1) The core CCC likely originated before the last common ancestor of unikonts and their closest bikont sister taxa. (2) Formation of the CCC may have constrained sequence evolution of the classical cadherin cytoplasmic domain and β-catenin in metazoa. (3) The α-catenin-binding domain in β-catenin appears to be the favored mutation site for disrupting β-catenin function in the CCC. (4) The ancestral function of the α/β-catenin heterodimer appears to be an actin-binding module. In some metazoan groups, more complex functions of α-catenin were gained by sequence divergence in the non-actin-binding (N-, M-) domains. (5) Allosteric regulation of α-catenin may have evolved for more complex regulation of the actin cytoskeleton.

New insights into the evolution of metazoan cadherins and catenins

E-Cadherin and β-catenin are the best studied representatives of the superfamilies of transmembrane cadherins and intracellular armadillo catenins, respectively. However, in over 600 million years of multicellular animal evolution, these two superfamilies have diversified remarkably both structurally and functionally. Although their basic building blocks, respectively, the cadherin repeat domain and the armadillo repeat domain, predate metazoans, the specific and complex domain compositions of the different family members and their functional roles in cell adhesion and signaling appear to be key features for the emergence of multicellular animal life. Basal animals such as placozoans and sponges have a limited number of distinct cadherins and catenins. The origin of vertebrates, in particular, coincided with a large increase in the number of cadherins and armadillo proteins, including modern "classical" cadherins, protocadherins, and plakophilins. Also, α-catenins increased. This chapter introduces the many different family members and describes the putative evolutionary relationships between them.

Asymmetrically expressed axin required for anterior development in Tribolium

Canonical Wnt signaling has been implicated in an AP axis polarizing mechanism in most animals, despite limited evidence from arthropods. In the long-germ insect, Drosophila, Wnt signaling is not required for global AP patterning, but in short-germ insects including Tribolium castaneum, loss of Wnt signaling affects development of segments in the growth zone but not those defined in the blastoderm. To determine the effects of ectopic Wnt signaling, we analyzed the expression and function of axin, which encodes a highly conserved negative regulator of the pathway. We found Tc-axin transcripts maternally localized to the anterior pole in freshly laid eggs. Expression spread toward the posterior pole during the early cleavage stages, becoming ubiquitous by the time the germ rudiment formed. Tc-axin RNAi produced progeny phenotypes that ranged from mildly affected embryos with cuticles displaying a graded loss of anterior structures, to defective embryos that condensed at the posterior pole in the absence of serosa. Altered expression domains of several blastodermal markers indicated anterior expansion of posterior fates. Analysis of other canonical Wnt pathway components and the expansion of Tc-caudal expression, a Wnt target, suggest that the effects of Tc-axin depletion are mediated through this pathway and that Wnt signaling must be inhibited for proper anterior development in Tribolium. These studies provide unique evidence that canonical Wnt signaling must be carefully regulated along the AP axis in an arthropod, and support an ancestral role for Wnt activity in defining AP polarity and patterning in metazoan development.