Heterochromatin, histone modifications, and nuclear architecture in disease vectors

Epigenetic studies in insects and the valproic acid perspective

Abstract Valproic acid in association with sodium valproate (VPA) is an important anticonvulsant drug used for decades to treat neurological disorders. VPA also acts as an epigenetic modulator by inhibiting histone deacetylases, permitting histone acetylation, affecting the DNA and histone methylation status and gene expression, and inducing chromatin remodeling. Insects represent an important animal model for studies in several areas of science. Their high phenotypic plasticity makes them alternative models for epigenetic studies. This brief review emphasizes recent reports on insect epigenetics and the contribution of studies on the VPA action in insects, including effects on epigenetic markers, extending the pharmacological understanding of the potential of this drug, and demonstrating the usefulness of insects as an alternative animal model to drug studies.

PyMEGABASE: Predicting cell-type-specific structural annotations of chromosomes using the epigenome

The folding patterns of interphase genomes in higher eukaryotes, as obtained from DNA-proximity-ligation or Hi-C experiments, are used to classify loci into structural classes called compartments and subcompartments. These structurally annotated (sub)compartments are known to exhibit specific epigenomic characteristics and cell-type-specific variations. To explore the relationship between genome structure and the epigenome, we present PyMEGABASE (PYMB), a maximum-entropy-based neural network model that predicts (sub)compartment annotations of a locus based solely on the local epigenome, such as ChIP-Seq of histone post-translational modifications. PYMB builds upon our previous model while improving robustness, capability to handle diverse inputs and user-friendly implementation. We employed PYMB to predict subcompartments for over a hundred human cell types available in ENCODE, shedding light on the links between subcompartments, cell identity, and epigenomic signals. The fact that PYMB, trained on data for human cells, can accurately predict compartments in mice suggests that the model is learning underlying physicochemical principles transferable across cell types and species. Reliable at higher resolutions (up to 5 kbp), PYMB is used to investigate compartment-specific gene expression. Not only can PYMB generate (sub)compartment information without Hi-C experiments, but its predictions are also interpretable. Analyzing PYMB's trained parameters, we explore the importance of various epigenomic marks in each subcompartment prediction. Furthermore, the predictions of the model can be used as input for OpenMiChroM software, which has been calibrated to generate three-dimensional structures of the genome. Detailed documentation of PYMB is available at https://pymegabase.readthedocs.io, including an installation guide using pip or conda, and Jupyter/Colab notebook tutorials.

Regulation of chromatin organization during animal regeneration

Activation of regeneration upon tissue damages requires the activation of many developmental genes responsible for cell proliferation, migration, differentiation, and tissue patterning. Ample evidence revealed that the regulation of chromatin organization functions as a crucial mechanism for establishing and maintaining cellular identity through precise control of gene transcription. The alteration of chromatin organization can lead to changes in chromatin accessibility and/or enhancer-promoter interactions. Like embryogenesis, each stage of tissue regeneration is accompanied by dynamic changes of chromatin organization in regeneration-responsive cells. In the past decade, many studies have been conducted to investigate the contribution of chromatin organization during regeneration in various tissues, organs, and organisms. A collection of chromatin regulators were demonstrated to play critical roles in regeneration. In this review, we will summarize the progress in the understanding of chromatin organization during regeneration in different research organisms and discuss potential common mechanisms responsible for the activation of regeneration response program.

Characterisation of the epigenetic architecture of the major malaria vector Anopheles arabiensis (Diptera: Culicidae) after treatment with epigenetic modulators and heavy metals

Anopheles arabiensis (a member of the An. gambiae species complex) is a major vector of malaria in sub-Saharan Africa. Despite its disease vector status, there is currently a paucity of epigenetic information for this species. The aim this study was therefore to analyse global epigenetic markers and their response to metal exposure in insecticide susceptible and resistant laboratory strains of An. arabiensis. This was done using commercially available epigenetic marker quantification kits. In order to validate the efficacy of the kits, several kits were assessed to determine whether changes induced by known epigenetic modulators were detectable using these platforms. The efficacy of the dosages used were determined by examining the effect of the dosages used on insecticide resistant phenotypes. Upon confirmation that the dosages used were sufficient to induce a phenotypic change, the effect on epigenetic markers was assessed. Commercial kits were used to quantify 5-methylcysteine (5-mC) and 5-hydroxymethylcysteine (5-hmC) methylation in DNA, m6A methylation in mRNA as well as Histone Acetyl Transferase (HAT) activity. There was a marked difference in the phenotypic response in adult mosquitoes of the insecticide susceptible strain compared to that of its' resistant counterpart. For males and females of the resistant strain, exposure to nucleic acid modifying drugs typically increased their tolerance to insecticides. The patterns of changes in 5-mC methylation by epigenetic modulators was congruent with previous studies which quantified by mass spectrometry. The two strains differed in methylation patterns under control conditions and responded differentially to larval metal exposure. In the resistant strain, which previously was demonstrated to show increased detoxification enzyme activity and insecticide tolerance after the same treatment, the potential increase in transcriptional activity appeared to be modulated by reduced methylation and increased HAT activity. This study suggests that the commercial epigenetic quantification kits can be used to characterise phenotypic changes in An. arabiensis, and also shows that epigenetic regulation of the response to metal exposure is regulated at the DNA as opposed to the RNA level.

Effects of sodium valproate on the chromatin of Triatoma infestans (Klug, 1834) (Hemiptera, Reduviidae) under in vitro culture conditions

Sodium valproate (VPA) is a classic anticonvulsive, a histone deacetylase inhibitor, and a chromatin remodeling inducer. When injected into specimens of Triatoma infestans, a vector of Chagas disease, VPA affects the chromatin supraorganization of chromocenter heterochromatin in only a few cells of the Malpighian tubules. To test whether this result was explained by the inaccessibility of all of the organ's cells to the drug, we investigated the nuclear phenotypes and global acetylation of lysine 9 in histone H3 (H3K9ac) in Malpighian tubules cultivated in vitro for 1-24 h in the presence of 0.05 mM-1 mM VPA. The present results revealed that the chromatin decondensation event in the chromocenter body, which was detected only under low VPA concentrations up to a 4-h treatment, was not frequent during organ culture, similar to the results for injected insects. Cultivation of T. infestans Malpighian tubules in vitro for 24 h revealed inadequate for cell preservation even in the absence of the drug. Immunofluorescence signals for H3K9ac following VPA treatment showed a slightly increased intensity in the euchromatin, but were never detected in the chromocenter bodies, except with great intensity at their periphery, where the 18S rDNA is located. In conclusion, when VPA affects the chromocenter heterochromatin in this animal cell model, it occurs through a pathway that excludes a classic global H3K9ac mark. Investigation of nonhistone proteins associated with histone methylation marks is still required to further explain the differential response of T. infestans chromatin to VPA.

A Gene-Based Method for Cytogenetic Mapping of Repeat-Rich Mosquito Genomes

Long-read sequencing technologies have opened up new avenues of research on the mosquito genome biology, enabling scientists to better understand the remarkable abilities of vectors for transmitting pathogens. Although new genome mapping technologies such as Hi-C scaffolding and optical mapping may significantly improve the quality of genomes, only cytogenetic mapping, with the help of fluorescence in situ hybridization (FISH), connects genomic scaffolds to a particular chromosome and chromosome band. This mapping approach is important for creating and validating chromosome-scale genome assemblies for mosquitoes with repeat-rich genomes, which can potentially be misassembled. In this study, we describe a new gene-based physical mapping approach that was optimized using the newly assembled Aedes albopictus genome, which is enriched with transposable elements. To avoid amplification of the repetitive DNA, 15 protein-coding gene transcripts were used for the probe design. Instead of using genomic DNA, complementary DNA was utilized as a template for development of the PCR-amplified probes for FISH. All probes were successfully amplified and mapped to specific chromosome bands. The genome-unique probes allowed to perform unambiguous mapping of genomic scaffolds to chromosome regions. The method described in detail here can be used for physical genome mapping in other insects.

Chromatin Structure and Function in Mosquitoes

The principles and function of chromatin and nuclear architecture have been extensively studied in model organisms, such as Drosophila melanogaster. However, little is known about the role of these epigenetic processes in transcriptional regulation in other insects including mosquitoes, which are major disease vectors and a worldwide threat for human health. Some of these life-threatening diseases are malaria, which is caused by protozoan parasites of the genus Plasmodium and transmitted by Anopheles mosquitoes; dengue fever, which is caused by an arbovirus mainly transmitted by Aedes aegypti; and West Nile fever, which is caused by an arbovirus transmitted by Culex spp. In this contribution, we review what is known about chromatin-associated mechanisms and the 3D genome structure in various mosquito vectors, including Anopheles, Aedes, and Culex spp. We also discuss the similarities between epigenetic mechanisms in mosquitoes and the model organism Drosophila melanogaster, and advocate that the field could benefit from the cross-application of state-of-the-art functional genomic technologies that are well-developed in the fruit fly. Uncovering the mosquito regulatory genome can lead to the discovery of unique regulatory networks associated with the parasitic life-style of these insects. It is also critical to understand the molecular interactions between the vectors and the pathogens that they transmit, which could hold the key to major breakthroughs on the fight against mosquito-borne diseases. Finally, it is clear that epigenetic mechanisms controlling mosquito environmental plasticity and evolvability are also of utmost importance, particularly in the current context of globalization and climate change.

Epigenetic Molecular Mechanisms in Insects

Insects are the largest animal group on Earth both in biomass and diversity. Their outstanding success has inspired genetics and developmental research, allowing the discovery of dynamic process explaining extreme phenotypic plasticity and canalization. Epigenetic molecular mechanisms (EMMs) are vital for several housekeeping functions in multicellular organisms, regulating developmental, ontogenetic trajectories and environmental adaptations. In Insecta, EMMs are involved in the development of extreme phenotypic divergences such as polyphenisms and eusocial castes. Here, we review the history of this research field and how the main EMMs found in insects help to understand their biological processes and diversity. EMMs in insects confer them rapid response capacity allowing insect either to change with plastic divergence or to keep constant when facing different stressors or stimuli. EMMs function both at intra as well as transgenerational scales, playing important roles in insect ecology and evolution. We discuss on how EMMs pervasive influences in Insecta require not only the control of gene expression but also the dynamic interplay of EMMs with further regulatory levels, including genetic, physiological, behavioral, and environmental among others, as was earlier proposed by the Probabilistic Epigenesis model and Developmental System Theory.

A chromosome-scale assembly of the major African malaria vector Anopheles funestus

BACKGROUND

Anopheles funestus is one of the 3 most consequential and widespread vectors of human malaria in tropical Africa. However, the lack of a high-quality reference genome has hindered the association of phenotypic traits with their genetic basis in this important mosquito.

FINDINGS

Here we present a new high-quality A. funestus reference genome (AfunF3) assembled using 240× coverage of long-read single-molecule sequencing for contigging, combined with 100× coverage of short-read Hi-C data for chromosome scaffolding. The assembled contigs total 446 Mbp of sequence and contain substantial duplication due to alternative alleles present in the sequenced pool of mosquitos from the FUMOZ colony. Using alignment and depth-of-coverage information, these contigs were deduplicated to a 211 Mbp primary assembly, which is closer to the expected haploid genome size of 250 Mbp. This primary assembly consists of 1,053 contigs organized into 3 chromosome-scale scaffolds with an N50 contig size of 632 kbp and an N50 scaffold size of 93.811 Mbp, representing a 100-fold improvement in continuity versus the current reference assembly, AfunF1.

CONCLUSION

This highly contiguous and complete A. funestus reference genome assembly will serve as an improved basis for future studies of genomic variation and organization in this important disease vector.

Histone acetylation and methylation marks in chromatin of Panstrongylus megistus (Hemiptera, Reduviidae)

Panstrongylus megistus, a potential vector of Chagas disease, currently occupies a wider geographic distribution in Brazil than Triatoma infestans, another member of the hemipteran Reduviidae family and a vector of the same disease. A small heterochromatic body (chromocenter) formed by the Y chromosome is evident in the somatic cells of P. megistus, differing in size and chromosome type contribution from the well-studied chromocenters present in T. infestans. While the overall distribution of histone epigenetic marks differ when comparing the heterochromatin and euchromatin territories in T. infestans, no similar data have been established for other hemipteran reduviids, including P. megistus. In the present work, histone acetylation and methylation marks were investigated in cells of Malpighian tubules of P. megistus 5th instar nymphs using immunocytochemical assays and compared to previously published data for T. infestans. Although similarities between these species were found regarding absence of acetylated H3K9, H4K8 and H4K16, and H3K9me and H3K9me₂ in the chromocenter, presence of these marks in euchromatin, and presence of H3K9me₃ in the chromocenter, no intimate association of acetylated H4K8 and 18S rDNA was revealed in the chromocenter of P. megistus. The elevated abundance of H3K9me₂ marks at the nuclear periphery in P. megistus cells, differing from data for T. infestans, is suggested to reflect differences in the interaction of lamina-associated chromatin domains with the nuclear lamina, methyl-transferase modulation and/or association with the last DNA endoreplication step in 5th instar nymphs, which is a matter for further investigation.

Management of arthropod vector data - Social and ecological dynamics facing the One Health perspective

Emerging infectious diseases (EIDs) are spread by direct and/or indirect contacts between a pathogen or parasite and their hosts. Arthropod vectors have evolved as excellent bloodsuckers, providing an elegant transportation mode for a wide number of infectious agents. The nature of pathogen and parasite transfer and the models used to predict how a disease might spread are magnified in complexity when an arthropod vector is part of the disease cycle. One Health is a worldwide strategy for expanding interdisciplinary collaborations and communications in all aspects of health care for humans, animals and the environment. It would benefit from a structured analysis to address vectoring of arthropod-borne diseases as a dynamic transactional process. This review focused on how arthropod vector data can be used to better model and predict zoonotic disease outbreaks. With enhanced knowledge to describe arthropod vector disease transfer, researchers will have a better understanding about how to model disease outbreaks. As public health research evolves to include more social-ecological systems, the roles of society, ecology, epidemiology, pathogen/parasite evolution and animal behavior can be better captured in the research design. Overall, because of more collaborative data collection processes on arthropod vectors, disease modeling can better predict conditions where EIDs will occur.

Comparative repeatome analysis on Triatoma infestans Andean and Non-Andean lineages, main vector of Chagas disease

Triatoma infestans is the most important Chagas disease vector in South America. Two main evolutionary lineages, named Andean and non-Andean, have been recognized by geographical distribution, phenetic and genetic characteristics. One of the main differences is the genomic size, varying over 30% in their haploid DNA content. Here we realize a genome wide analysis to compare the repetitive genome fraction (repeatome) between both lineages in order to identify the main repetitive DNA changes occurred during T. infestans differentiation process. RepeatExplorer analysis using Illumina reads showed that both lineages exhibit the same amount of non-repeat sequences, and that satellite DNA is by far the major component of repetitive DNA and the main responsible for the genome size differentiation between both lineages. We characterize 42 satellite DNA families, which are virtually all present in both lineages but with different amount in each lineage. Furthermore, chromosomal location of satellite DNA by fluorescence in situ hybridization showed that genomic variations in T. infestans are mainly due to satellite DNA families located on the heterochromatic regions. The results also show that many satDNA families are located on the euchromatic regions of the chromosomes.

Tissue-specific features of the X chromosome and nucleolus spatial dynamics in a malaria mosquito, Anopheles atroparvus

Spatial organization of chromosome territories is important for maintenance of genomic stability and regulation of gene expression. Recent studies have shown tissue-specific features of chromosome attachments to the nuclear envelope in various organisms including malaria mosquitoes. However, other spatial characteristics of nucleus organization, like volume and shape of chromosome territories, have not been studied in Anopheles. We conducted a thorough analysis of tissue-specific features of the X chromosome and nucleolus volume and shape in follicular epithelium and nurse cells of the Anopheles atroparvus ovaries using a modern open-source software. DNA of the polytene X chromosome from ovarian nurse cells was obtained by microdissection and was used as a template for amplification with degenerate oligo primers. A fluorescently labeled X chromosome painting probe was hybridized with formaldehyde-fixed ovaries of mosquitoes using a 3D-FISH method. The nucleolus was stained by immunostaining with an anti-fibrillarin antibody. The analysis was conducted with TANGO-a software for a chromosome spatial organization analysis. We show that the volume and position of the X chromosome have tissue-specific characteristics. Unlike nurse cell nuclei, the growth of follicular epithelium nuclei is not accompanied with the proportional growth of the X chromosome. However, the shape of the X chromosome does not differ between the tissues. The dynamics of the X chromosome attachment regions location is tissue-specific and it is correlated with the process of nucleus growth in follicular epithelium and nurse cells.

The Physical Genome Mapping of Anopheles albimanus Corrected Scaffold Misassemblies and Identified Interarm Rearrangements in Genus Anopheles

The genome of the Neotropical malaria vector Anopheles albimanus was sequenced as part of the 16 Anopheles Genomes Project published in 2015. The draft assembly of this species consisted of 204 scaffolds with an N50 scaffold size of 18.1 Mb and a total assembly size of 170.5 Mb. It was among the smallest genomes with the longest scaffolds in the 16 Anopheles species cluster, making An. albimanus the logical choice for anchoring the genome assembly to chromosomes. In this study, we developed a high-resolution cytogenetic photomap with completely straightened polytene chromosomes from the salivary glands of the mosquito larvae. Based on this photomap, we constructed a chromosome-based genome assembly using fluorescent in situ hybridization of PCR-amplified DNA probes. Our physical mapping, assisted by an ortholog-based bioinformatics approach, identified and corrected nine misassemblies in five large genomic scaffolds. Misassemblies mostly occurred in junctions between contigs. Our comparative analysis of scaffolds with the An. gambiae genome detected multiple genetic exchanges between pericentromeric regions of chromosomal arms caused by partial-arm translocations. The final map consists of 40 ordered genomic scaffolds and corrected fragments of misassembled scaffolds. The An. albimanus physical map comprises 98.2% of the total genome assembly and represents the most complete genome map among mosquito species. This study demonstrates that physical mapping is a powerful tool for correcting errors in draft genome assemblies and for creating chromosome-anchored reference genomes.

Chromosomal Translocations in Black Flies (Diptera: Simuliidae)-Facilitators of Adaptive Radiation?

A macrogenomic investigation of a Holarctic clade of black flies-the Simulium cholodkovskii lineage-provided a platform to explore the implications of a unique, synapomorphic whole-arm interchange in the evolution of black flies. Nearly 60 structural rearrangements were discovered in the polytene complement of the lineage, including 15 common to all 138 analyzed individuals, relative to the central sequence for the entire subgenus Simulium. Three species were represented, of which two Palearctic entities (Simulium cholodkovskii and S. decimatum) were sympatric; an absence of hybrids confirmed their reproductive isolation. A third (Nearctic) entity had nonhomologous sex chromosomes, relative to the other species, and is considered a separate species, for which the name Simulium nigricoxum is revalidated. A cytophylogeny is inferred and indicates that the two Palearctic taxa are sister species and these, in turn, are the sister group of the Nearctic species. The rise of the S. cholodkovskii lineage encompassed complex chromosomal and genomic restructuring phenomena associated with speciation in black flies, viz. expression of one and the same rearrangement as polymorphic, fixed, or sex linked in different species; taxon-specific differentiation of sex chromosomes; and reciprocal translocation of chromosome arms. The translocation is hypothesized to have occurred early in male spermatogonia, with the translocated chromosomal complement being transmitted to the X- and Y-bearing sperm during spermatogenesis, resulting in alternate disjunction of viable F1 translocation heterozygotes and the eventual formation of more viable and selectable F2 translocation homozygous progeny. Of 11 or 12 independently derived whole-arm interchanges known in the family Simuliidae, at least six are associated with subsequent speciation events, suggesting a facilitating role of translocations in adaptive radiations. The findings are discussed in the context of potential structural and functional interactions for future genomic research.

Chromosome evolution in malaria mosquitoes inferred from physically mapped genome assemblies

Polymorphic inversions in mosquitoes are distributed nonrandomly among chromosomes and are associated with ecological, behavioral, and physiological adaptations related to pathogen transmission. Despite their significance, the patterns and mechanism of genome rearrangements are not well understood. Recent sequencing and physical mapping of the genomes for 16 Anopheles mosquito species provided an opportunity to study chromosome evolution at the highest resolution. New studies revealed that fixed rearrangement accumulated [Formula: see text]3 times faster on the X chromosome than on autosomes. The highest densities of transposable elements (TEs) and satellites of different sizes have also been found on the X chromosome, suggesting a mechanism for the inversion generation. The high rate of X chromosome rearrangements is in sharp contrast with the paucity of polymorphic inversions on the X in the majority of anopheline species. This paper highlights the advances in understanding chromosome evolution in malaria vectors and discusses possible future directions in studying mechanisms and biological roles of genome rearrangements.

Increased production of piRNAs from euchromatic clusters and genes in Anopheles gambiae compared with Drosophila melanogaster

BACKGROUND

Specific genomic loci, termed Piwi-interacting RNA (piRNA) clusters, manufacture piRNAs that serve as guides for the inactivation of complementary transposable elements (TEs). The piRNA pathway has been accurately detailed in Drosophila melanogaster, while it remains poorly examined in other insects. This pathway is increasingly recognized as critical for germline development and reproduction. Understanding of the piRNA functions in mosquitoes could offer an opportunity for disease vector control by the reduction of their reproductive potential.

RESULTS

To analyze the similarities and differences in this pathway between Drosophila and mosquito, we performed an in-depth analysis of the genomic loci producing piRNAs and their targets in the African malaria vector Anopheles gambiae. We identified 187 piRNA clusters in the An. gambiae genome and 155 piRNA clusters in the D. melanogaster genome. We demonstrate that many more piRNA clusters in the mosquito compared with the fruit fly are uni-directionally transcribed and are located outside pericentromeric heterochromatin. About 11 % of the An. gambiae piRNA population map to gene transcripts. This is a noticeable increase compared with the ~6 % of the piRNA population mapped to genes in D. melanogaster. A subset of the piRNA-enriched genes in An. gambiae has functions related to reproduction and development. At least 24 and 65 % of the mapped piRNAs correspond to genomic TE sequences in An. gambiae and D. melanogaster, respectively. DNA transposons and non-LTR retrotransposons are more abundant in An. gambiae, while LTR retrotransposons are more abundant in D. melanogaster. Yet, piRNAs predominantly target LTR retrotransposons in both species, which may point to a distinct feature of these elements compared to the other classes of TEs concerning their silencing by the piRNA pathway.

CONCLUSIONS

Here, we demonstrate that piRNA-producing loci have more ubiquitous distribution in the An. gambiae genome than in the genome of D. melanogaster. Also, protein-coding genes have an increased role in production of piRNAs in the germline of this mosquito. Genes involved in germline and embryonic development of An. gambiae generate a substantial portion of piRNAs, suggesting a role of the piRNA pathway in the epigenetic regulation of the reproductive processes in the African malaria vector.