Transcriptomic data from panarthropods shed new light on the evolution of insulator binding proteins in insects : Insect insulator proteins

Unique territorial and compartmental organization of chromosomes in the holocentric silkmoth

The hallmarks of chromosome organization in multicellular eukaryotes are chromosome territories (CT), chromatin compartments, and insulated domains, including topologically associated domains (TADs). Yet, most of these elements of chromosome organization are derived from analyses of a limited set of model organisms, while large eukaryotic groups, including insects, remain mostly unexplored. Here we combine Hi-C, biophysical modeling, and microscopy to characterize the 3D genome architecture of the silkworm, Bombyx mori. In contrast to other eukaryotes, B. mori chromosomes form highly separated territories. Similar to other eukaryotes, B. mori chromosomes segregate into active A and inactive B compartments, yet unlike in vertebrate systems, contacts between euchromatic A regions appear to be a strong driver of compartmentalization. Remarkably, we also identify a third compartment, called secluded S, with a unique contact pattern. Each S region shows prominent short-range self-contacts and is remarkably devoid of contacts with the rest of the chromosome, including other S regions. Compartment S hosts a unique combination of genetic and epigenetic features, localizes towards the periphery of CTs, and shows developmental plasticity. Biophysical modeling reveals that the formation of such secluded domains requires highly localized loop extrusion within them, along with a low level of extrusion in A and B. Our Hi-C data supports predicted genome-wide and localized extrusion. Such a broad, non-uniform distribution of extruders has not been seen in other organisms. Overall, our analyses support loop extrusion in insects and highlight the evolutionary plasticity of 3D genome organization, driven by a new combination of known processes.

Functional genomics in Spiralia

Our understanding of the mechanisms that modulate gene expression in animals is strongly biased by studying a handful of model species that mainly belong to three groups: Insecta, Nematoda and Vertebrata. However, over half of the animal phyla belong to Spiralia, a morphologically and ecologically diverse animal clade with many species of economic and biomedical importance. Therefore, investigating genome regulation in this group is central to uncovering ancestral and derived features in genome functioning in animals, which can also be of significant societal impact. Here, we focus on five aspects of gene expression regulation to review our current knowledge of functional genomics in Spiralia. Although some fields, such as single-cell transcriptomics, are becoming more common, the study of chromatin accessibility, DNA methylation, histone post-translational modifications and genome architecture are still in their infancy. Recent efforts to generate chromosome-scale reference genome assemblies for greater species diversity and optimise state-of-the-art approaches for emerging spiralian research systems will address the existing knowledge gaps in functional genomics in this animal group.

Evolution of 3D chromatin organization at different scales

Most animal genomes fold in 3D chromatin domains called topologically associated domains (TADs) that facilitate interactions between cis-regulatory elements (CREs) and promoters. Owing to their critical role in the control of developmental gene expression, we explore how TADs have shaped animal evolution. In the light of recent studies that profile TADs in disparate animal lineages, we discuss their phylogenetic distribution and the mechanisms that underlie their formation. We present evidence indicating that TADs are plastic entities composed of genomic strata of different ages: ancient cores are combined with newer regions and brought into extant TADs through genomic rearrangements. We highlight that newly incorporated TAD strata enable the establishment of new CRE-promoter interactions and in turn new expression patterns that can drive phenotypical innovation. We further highlight how subtle changes in chromatin folding may fine-tune the expression levels of developmental genes and hold a potential for evolutionary significance.

Evolutionary history and divergence times of Odonata (dragonflies and damselflies) revealed through transcriptomics

Dragonflies and damselflies are among the earliest flying insects with extant representatives. However, unraveling details of their long evolutionary history, such as egg laying (oviposition) strategies, is impeded by unresolved phylogenetic relationships, particularly in damselflies. Here we present a transcriptome-based phylogenetic reconstruction of Odonata, analyzing 2,980 protein-coding genes in 105 species representing nearly all the order's families. All damselfly and most dragonfly families are recovered as monophyletic. Our data suggest a sister relationship between dragonfly families of Gomphidae and Petaluridae. According to our divergence time estimates, both crown-Zygoptera and -Anisoptera arose during the late Triassic. Egg-laying with a reduced ovipositor apparently evolved in dragonflies during the late Jurassic/early Cretaceous. Lastly, we also test the impact of fossil choice and placement, particularly, of the extinct fossil species, †Triassolestodes asiaticus, and †Proterogomphus renateae on divergence time estimates. We find placement of †Proterogomphus renateae to be much more impactful than †Triassolestodes asiaticus.

Four myriapod relatives - but who are sisters? No end to debates on relationships among the four major myriapod subgroups

Background

Phylogenetic relationships among the myriapod subgroups Chilopoda, Diplopoda, Symphyla and Pauropoda are still not robustly resolved. The first phylogenomic study covering all subgroups resolved phylogenetic relationships congruently to morphological evidence but is in conflict with most previously published phylogenetic trees based on diverse molecular data. Outgroup choice and long-branch attraction effects were stated as possible explanations for these incongruencies. In this study, we addressed these issues by extending the myriapod and outgroup taxon sampling using transcriptome data.

Results

We generated new transcriptome data of 42 panarthropod species, including all four myriapod subgroups and additional outgroup taxa. Our taxon sampling was complemented by published transcriptome and genome data resulting in a supermatrix covering 59 species. We compiled two data sets, the first with a full coverage of genes per species (292 single-copy protein-coding genes), the second with a less stringent coverage (988 genes). We inferred phylogenetic relationships among myriapods using different data types, tree inference, and quartet computation approaches. Our results unambiguously support monophyletic Mandibulata and Myriapoda. Our analyses clearly showed that there is strong signal for a single unrooted topology, but a sensitivity of the position of the internal root on the choice of outgroups. However, we observe strong evidence for a clade Pauropoda+Symphyla, as well as for a clade Chilopoda+Diplopoda.

Conclusions

Our best quartet topology is incongruent with current morphological phylogenies which were supported in another phylogenomic study. AU tests and quartet mapping reject the quartet topology congruent to trees inferred with morphological characters. Moreover, quartet mapping shows that confounding signal present in the data set is sufficient to explain the weak signal for the quartet topology derived from morphological characters. Although outgroup choice affects results, our study could narrow possible trees to derivatives of a single quartet topology. For highly disputed relationships, we propose to apply a series of tests (AU and quartet mapping), since results of such tests allow to narrow down possible relationships and to rule out confounding signal.

Evolutionary History of Major Chemosensory Gene Families across Panarthropoda

Chemosensory perception is a fundamental biological process of particular relevance in basic and applied arthropod research. However, apart from insects, there is little knowledge of specific molecules involved in this system, which is restricted to a few taxa with uneven phylogenetic sampling across lineages. From an evolutionary perspective, onychophorans (velvet worms) and tardigrades (water bears) are of special interest since they represent the closest living relatives of arthropods, altogether comprising the Panarthropoda. To get insights into the evolutionary origin and diversification of the chemosensory gene repertoire in panarthropods, we sequenced the antenna- and head-specific transcriptomes of the velvet worm Euperipatoides rowelli and analyzed members of all major chemosensory families in representative genomes of onychophorans, tardigrades, and arthropods. Our results suggest that the NPC2 gene family was the only family encoding soluble proteins in the panarthropod ancestor and that onychophorans might have lost many arthropod-like chemoreceptors, including the highly conserved IR25a receptor of protostomes. On the other hand, the eutardigrade genomes lack genes encoding the DEG-ENaC and CD36-sensory neuron membrane proteins, the chemosensory members of which have been retained in arthropods; these losses might be related to lineage-specific adaptive strategies of tardigrades to survive extreme environmental conditions. Although the results of this study need to be further substantiated by an increased taxon sampling, our findings shed light on the diversification of chemosensory gene families in Panarthropoda and contribute to a better understanding of the evolution of animal chemical senses.

Overlapping but Distinct Sequences Play Roles in the Insulator and Promoter Activities of the Drosophila BEAF-Dependent scs' Insulator

Chromatin domain insulators are thought to help partition the genome into genetic units called topologically associating domains (TADs). In Drosophila, TADs are often separated by inter-TAD regions containing active housekeeping genes and associated insulator binding proteins. This raises the question of whether insulator binding proteins are involved primarily in chromosomal TAD architecture or gene activation, or if these two activities are linked. The Boundary Element-Associated Factor of 32 kDa (BEAF-32, or BEAF for short) is usually found in inter-TADs. BEAF was discovered based on binding to the scs' insulator, and is important for the insulator activity of scs' and other BEAF binding sites. There are divergent promoters in scs' with a BEAF binding site by each. Here, we dissect the scs' insulator to identify DNA sequences important for insulator and promoter activity, focusing on the half of scs' with a high affinity BEAF binding site. We find that the BEAF binding site is important for both insulator and promoter activity, as is another sequence we refer to as LS4. Aside from that, different sequences play roles in insulator and promoter activity. So while there is overlap and BEAF is important for both, insulator and promoter activity can be separated.

Promoter-Proximal Chromatin Domain Insulator Protein BEAF Mediates Local and Long-Range Communication with a Transcription Factor and Directly Activates a Housekeeping Promoter in Drosophila

BEAF (Boundary Element-Associated Factor) was originally identified as a Drosophila melanogaster chromatin domain insulator-binding protein, suggesting a role in gene regulation through chromatin organization and dynamics. Genome-wide mapping found that BEAF usually binds near transcription start sites, often of housekeeping genes, suggesting a role in promoter function. This would be a nontraditional role for an insulator-binding protein. To gain insight into molecular mechanisms of BEAF function, we identified interacting proteins using yeast two-hybrid assays. Here, we focus on the transcription factor Serendipity δ (Sry-δ). Interactions were confirmed in pull-down experiments using bacterially expressed proteins, by bimolecular fluorescence complementation, and in a genetic assay in transgenic flies. Sry-δ interacted with promoter-proximal BEAF both when bound to DNA adjacent to BEAF or > 2-kb upstream to activate a reporter gene in transient transfection experiments. The interaction between BEAF and Sry-δ was detected using both a minimal developmental promoter (y) and a housekeeping promoter (RpS12), while BEAF alone strongly activated the housekeeping promoter. These two functions for BEAF implicate it in playing a direct role in gene regulation at hundreds of BEAF-associated promoters.

Multifunctional cellulase enzymes are ancestral in Polyneoptera

Many hemimetabolous insects produce their own cellulase enzymes from the glycoside hydrolase family 9, first observed in termites and cockroaches. Phasmatodea have multiple cellulases, some of which are multifunctional and can degrade xylan or xyloglucan. To discover when these abilities evolved, we identified cellulases from the Polyneoptera sampled by the 1000 Insect Transcriptome and Evolution (1KITE) project, including all cockroach and termite transcriptomes. We hoped to identify what role enzyme substrate specificities had in the evolution of dietary specification, such as leaf-feeding or wood-feeding. Putative cellulases were identified from the transcriptomes and analysed phylogenetically. All cellulases were amplified from an exemplar set of Polyneoptera species using rapid amplification of cDNA ends PCR and heterologously expressed in an insect cell line, then tested against different polysaccharides for their digestive abilities. We identified several multifunctional xyloglucanolytic enzymes across Polyneoptera, plus a large group of cellulase-like enzymes found in nearly all insect orders with no discernible digestive ability. Multifunctional xylanolytic cellulases remain unique to Phasmatodea. The presence or absence of multifunctional enzymes does not impact dietary specification, but rather having multiple, multifunctional cellulase genes is an ancestral state for Polyneoptera and possibly Insecta. The prevalence of multifunctional cellulases in other animals demands further investigation.

Genomic data reveal high conservation but divergent evolutionary pattern of Polycomb/Trithorax group genes in arthropods

Epigenetic gene control is maintained by chromatin-associated Polycomb group (PcG) and Trithorax group (TrxG) genes, which act antagonistically via the interplay between PcG and TrxG regulation to generate silenced or permissive transcriptional states. In this study, we searched for PcG/TrxG genes in 180 arthropod genomes, covering all the sequenced arthropod genomes at the time of conducting this study, to perform a global investigation of PcG/TrxG genes in a phylogenetic frame. Results of ancestral state reconstruction analysis revealed that the ancestor of arthropod species has an almost complete repertoire of PcG/TrxG genes, and most of these genes were seldom lost above order level. The domain diversity analysis indicated that the PcG/TrxG genes show variable extent of domain structure changes; some of these changes could be associated with lineage-specific events. The likelihood ratio tests for selection pressure detected a number of PcG/TrxG genes which underwent episodic positive selection on the branch leading to the insects with holometabolous development. These results suggest that, despite their high conservation across arthropod species, different members of PcG/TrxG genes showed considerable differences in domain structure and sequence divergence in arthropod evolution. Our cross species comparisons using large-scale genomic data provide insights into divergent evolutionary pattern on highly conserved genes in arthropods.

Role of Su(Hw) zinc finger 10 and interaction with CP190 and Mod(mdg4) proteins in recruiting the Su(Hw) complex to chromatin sites in Drosophila

Su(Hw) belongs to the class of proteins that organize chromosome architecture and boundaries/insulators between regulatory domains. This protein contains a cluster of 12 zinc finger domains most of which are responsible for binding to three different modules in the consensus site. Su(Hw) forms a complex with CP190 and Mod(mdg4)-67.2 proteins that binds to well-known Drosophila insulators. To understand how Su(Hw) performs its activities and binds to specific sites in chromatin, we have examined the previously described su(Hw)^f mutation that disrupts the 10th zinc finger (ZF10) responsible for Su(Hw) binding to the upstream module. The results have shown that Su(Hw)^f loses the ability to interact with CP190 in the absence of DNA. In contrast, complete deletion of ZF10 does not prevent the interaction between Su(Hw)^Δ10 and CP190. Having studied insulator complex formation in different mutant backgrounds, we conclude that both association with CP190 and Mod(mdg4)-67.2 partners and proper organization of DNA binding site are essential for the efficient recruitment of the Su(Hw) complex to chromatin insulators.