Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow lsa Template on Twitter
Research Article
Transparent Process
Open Access

Distinct mechanisms mediate X chromosome dosage compensation in Anopheles and Drosophila

View ORCID ProfileClaudia Isabelle Keller Valsecchi  Correspondence email, View ORCID ProfileEric Marois, View ORCID ProfileM Felicia Basilicata, Plamen Georgiev, View ORCID ProfileAsifa Akhtar  Correspondence email
Claudia Isabelle Keller Valsecchi
1Max Planck Institute of Immunobiology and Epigenetics, Freiburg, Germany
3Institute of Molecular Biology (IMB), Mainz, Germany
Roles: Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Project administration, Writing—original draft, review, and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Claudia Isabelle Keller Valsecchi
  • For correspondence: c.keller@imb-mainz.de
Eric Marois
2Université de Strasbourg, Centre National de la Recherche Scientifique (CNRS) UPR9022, Inserm U1257, Strasbourg, France
Roles: Conceptualization, Resources, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Eric Marois
M Felicia Basilicata
1Max Planck Institute of Immunobiology and Epigenetics, Freiburg, Germany
3Institute of Molecular Biology (IMB), Mainz, Germany
Roles: Formal analysis, Supervision, Validation, Investigation, Visualization, Methodology, Project administration, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for M Felicia Basilicata
Plamen Georgiev
1Max Planck Institute of Immunobiology and Epigenetics, Freiburg, Germany
Roles: Conceptualization, Resources, Investigation, Methodology, Project administration, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Asifa Akhtar
1Max Planck Institute of Immunobiology and Epigenetics, Freiburg, Germany
Roles: Conceptualization, Resources, Supervision, Funding acquisition, Methodology, Project administration, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Asifa Akhtar
  • For correspondence: akhtar@ie-freiburg.mpg.de
Published 15 July 2021. DOI: 10.26508/lsa.202000996
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

Sex chromosomes induce potentially deleterious gene expression imbalances that are frequently corrected by dosage compensation (DC). Three distinct molecular strategies to achieve DC have been previously described in nematodes, fruit flies, and mammals. Is this a consequence of distinct genomes, functional or ecological constraints, or random initial commitment to an evolutionary trajectory? Here, we study DC in the malaria mosquito Anopheles gambiae. The Anopheles and Drosophila X chromosomes evolved independently but share a high degree of homology. We find that Anopheles achieves DC by a mechanism distinct from the Drosophila MSL complex–histone H4 lysine 16 acetylation pathway. CRISPR knockout of Anopheles msl-2 leads to embryonic lethality in both sexes. Transcriptome analyses indicate that this phenotype is not a consequence of defective X chromosome DC. By immunofluorescence and ChIP, H4K16ac does not preferentially enrich on the male X. Instead, the mosquito MSL pathway regulates conserved developmental genes. We conclude that a novel mechanism confers X chromosome up-regulation in Anopheles. Our findings highlight the pluralism of gene-dosage buffering mechanisms even under similar genomic and functional constraints.

Introduction

Chromosome-wide copy number alterations are typically linked to developmental failure and diseases such as cancer. Heteromorphic sex chromosomes are a notable exception. To account for potentially deleterious expression imbalances in XY males in comparison with XX females, X chromosomes are often subjected to dosage compensation (DC) (Samata & Akhtar, 2018). In Dipterans such as fruit flies or mosquitos, sex chromosomes have formed multiple times (Vicoso & Bachtrog, 2015) and show high evolutionary turnover (Herpin & Schartl, 2015). RNA-seq studies revealed that X chromosome DC in multiple Dipteran genera occurs by up-regulation of the single male X (Rose et al, 2016). Owing to their short generation times and related, yet rapidly evolving genomes, Dipteran species are excellent models to understand how cells co-opt pathways to keep chromosome-wide gene dosage alterations in check. Dipterans have also been intensively studied as some of its species are prominent disease vectors, for example, the malaria-transmitting mosquito Anopheles gambiae. Understanding the processes responsible for sex differentiation and regulation of heteromorphic sex chromosomes is highly relevant as these factors are among the most promising targets for gene-drive strategies aiming at preventing transmission of infectious diseases by vector control (Kyrou et al, 2018).

DC in the fruit fly Drosophila melanogaster (Dmel) is mediated by association of the male-specific lethal (MSL) complex with the X chromosome in males. Via deposition of histone H4 lysine 16 acetylation (H4K16ac), the MSL complex enhances transcription to achieve approximately twofold up-regulated expression of X-linked genes (Lucchesi & Kuroda, 2015; Samata & Akhtar, 2018). Loss of MSL complex member genes results in male-specific lethality, whereas females remain unaffected (Belote & Lucchesi, 1980). Key to establishing MSL-DC is its subunit MSL2, which is male-specifically expressed and recruits the complex to the X chromosome (Bashaw & Baker, 1995). H4K16ac enrichment on the X as well as other characteristic features of the MSL-DC system are conserved in other drosophilids, for example, the distantly related Drosophila virilis or Drosophila busckii (Alekseyenko et al, 2013; Renschler et al, 2019). The MSL complex is also present in mammals, but its in vivo functions remain poorly characterized and mammalian DC is mediated by X chromosome inactivation in females via an entirely different pathway (Keller & Akhtar, 2015).

Like Drosophila, the malaria mosquito A. gambiae (Agam) achieves DC by hypertranscription of the male X (Deitz et al, 2018). The Dmel and Agam X chromosomes evolved independently from the same ancestral autosome (Zdobnov et al, 2002). It remains unclear whether X up-regulation in mosquitos is achieved by the same molecular machinery as in drosophilids or whether different mechanisms may arise in these contexts. This is of particular interest, as similar functional constraints and gene contents have shaped the acquisition of DC on newly evolving sex chromosomes in these lineages. Hence, comparing Agam with Dmel offers an excellent opportunity to study whether the independent evolution of sex chromosomes is accompanied by the acquisition of the same DC mechanism.

Here, we analyzed whether the MSL complex–H4K16ac pathway confers X chromosome–wide DC in Agam. Expression analyses indicate that unlike in Drosophila, the MSL complex subunits are expressed in both males and females in Anopheles. Loss of Anopheles msl-2 affects both sexes, where zygotic KO induces an early embryonic lethal phenotype and germ line KO is results in sterility. Transcriptome analyses indicate that msl-2 KO embryos do not display a global perturbation of X-linked gene expression. Furthermore, ChIP-seq and immunofluorescence show that H4K16ac is neither enriched on the male mosquito X chromosome nor displays a sexually dimorphic pattern. Instead, the MSL–H4K16 pathway appears to be involved in regulating developmental genes. Our study highlights that X chromosome DC can be achieved by entirely different molecular strategies even in very similar genomic contexts.

Results

MSL complex members are expressed in both male and female A. gambiae

Cytogenetic and genomic studies suggest that the ancestral dipteran karyotype consisted of six chromosomal pairs commonly referred to as Muller elements A-F (Vicoso & Bachtrog, 2015). Two independent evolutionary events resulted in Muller element A becoming sex-linked in Drosophila and Anopheles (Fig 1A) (Landeen & Presgraves, 2013; Vicoso & Bachtrog, 2015). The X chromosomes of these species exhibit considerable homology (Fig 1B [Zdobnov et al, 2002]), with around half of the Agam X-linked genes sharing an annotated X-linked orthologue in Dmel (Fig 1C).

Figure 1
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1 MSL complex members are expressed in both male and female Anopheles gambiae.

(A) Schematic representation of the evolution of the X chromosomes in the indicated species from the ancient Dipteran karyotype consisting of Muller elements A-F. Note that the mosquito Aedes aegyptii does not contain a differentiated sex chromosome. (B) Circos plot of pairs of synteny blocks that were generated over genomic regions of at least 100 kb among Drosophila melanogaster and Anopheles gambiae genomes. Synteny blocks on the Dmel X chromosome are shaded in orange, autosomes in grey. (C) Bar plot illustrating the chromosomal location of 1:1 orthologues of all A. gambiae genes in the D. melanogaster genome. If a given A. gambiae gene does not contain an annotated D. melanogaster orthologue, the bar is colored in black, X-linked genes in orange and genes on the autosomal arms in grey. Orthologues were obtained from VectorBase. (D) Amino acid sequence conservation of the individual domains of D. melanogaster MSL proteins obtained by BLASTP in A. gambiae. Detailed results are reported in Table S1. (E) RT-qPCR analysis of polyA+ RNA levels of the indicated genes in adult A. gambiae females and males. The mRNA level of each gene was normalized to Rps4 and expressed relative to the expression level in males. The bar plot represents the mean ± SEM with overlaid data points reflecting one biological replicate. Rps4, RpL32, and RPL19 are control genes, AGAP000634 and GPRMTH2 are representative X-linked and autosomal genes, respectively. Significance was evaluated by a Wilcoxon rank sum test, *P-value < 0.05, **P-value < 0.01, details in Table S1. (F) Line plot representing the mean log2 expression levels of msl genes during A. gambiae embryogenesis. The data from Goltsev et al (2009) was obtained from VectorBase. (G) Cropped immunoblots of male and female extracts of msl-2::HA (endogenously tagged) D. melanogaster heads. RNA Polymerase 2 (POL2) serves as loading control. (H), as in (G), Tubulin serves as loading control. The asterisk indicates the major isoform of MSL3 that is only expressed in males, the band detected at 50 kD is unspecific. The band below 50 kD accumulates only in females. (I) Left: schematic representation of the polysome profiling strategy to assess the protein expression of msls. Extracts are fractionated on a sucrose gradient to separate mRNA that associates with actively translating ribosomes from free/not-translated RNA. Right: bar plot displaying the relative enrichment of msls and controls (not translated: U6 snRNA; translated: Rps4 mRNA) in the polysomal fractions. The data for all fractions individually is reported in Fig S1F. The barplot represents the mean ± SEM with overlaid data points reflecting one polysome profiling experiment/biological replicate.

To elucidate whether the MSL complex mediates DC in Agam, we first analyzed the amino acid sequence conservation of the core MSL subunits MSL1, MSL2, MSL3, and MOF, as well as the associated RNA helicase MLE (Table S1 and Fig 1D; note that the high evolutionary turnover of noncoding RNAs hinders the computational identification of putative orthologues of Dmel roX1 and roX2 [Quinn et al, 2016]). The highest degree of conservation was found for MLE (57% identity) and MOF (61% identity), the latter of which displays a shortened N-terminus in Anopheles compared to Drosophila. MSL3 displays a similar length and more than 50% similarity for both the chromo- and MRG domains. The RING and CXC domains of MSL2 are also conserved, but the flexible linker regions adjacent to these structured domains are substantially extended in Anopheles. The DNA binding CXC domain of Agam MSL2 appears more similar to mammalian than Drosophila MSL2 (Fig S1A). There is no annotated Agam msl-1 gene and extensive BLAST searches using the full-length protein failed to identify a bona-fide orthologue of MSL1. Because MSL1 is largely unstructured (Hallacli et al, 2012), the overall sequence level conservation may be rather low. Indeed, individual domain searches revealed hits for the coiled-coil and PEHE regions, albeit at rather low fidelity (Table S1; note that even mammalian and Drosophila MSL1 show <20% amino acid similarity despite providing the same scaffolding function to the complex).

Table S1 Drosophila melanogaster MSL protein sequences and results of BLAST searches in Anopheles gambiae, RNA-seq mapping statistics, and statistical test results.

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Analysis of MSL complex components in male and female Anopheles gambiae.

(A) Sequence alignment of the MSL2 RING (top) and CXC (bottom) domains in Drosophila melanogaster (dm), Drosophila virilis (dv), A. gambiae, and Homo sapiens. Residues that are conserved in all four species are in green, for the CXC domain, residues that are conserved in any of the three other species are additionally shaded with a green background. (B) Bar plot of the mean log2 expression level and STDEV of msl genes, which were obtained from RNA-seq data of carcasses (somatic) and reproductive tissue of adult males and females published by Papa et al (2017). Padj and classification/bias values (legend) were obtained from Papa et al (2017). (C) Bar plot of the mean log2 expression level of msl genes in adult D. melanogaster males and females of the indicated age. The data were obtained from modEncode and downloaded from FlyBase. (D) RT-qPCR analysis of RNA levels of the indicated genes in adult D. melanogaster females and males. The RNA level of each gene was normalized to RpL32 and expressed relative to the expression level in males. The bar plot represents the mean ± SEM with overlaid data points reflecting one biological replicate. -RA and -RB indicate two isoforms of msl-2 and msl-1, respectively, as annotated in FlyBase. Significance was evaluated by a Wilcoxon rank sum test, *P-value < 0.05, details in Table S1. (E) Line plot of the mean expression level (rpkm) during development of D. melanogaster. The data were obtained from modEncode and downloaded from FlyBase. (F) Line plot displaying the mean ± SEM of n = 3 replicates of the relative RNA enrichment of msl genes and controls (snRNA U6, not translated; Rps4 mRNA, translated) in each of the polysome-profile fractions (left: females; right: males). Fraction 1 corresponds to the top of the gradient, Fraction 8 to the bottom. Large particles such as ribosomes and polysomes sediment in the high % sucrose-containing fractions (4–8).

Assessing the expression level of the core MSL complex members by RT-qPCR from adult Anopheles, we found a modestly biased expression of mof and mle in males (Fig 1E, statistical analyses in Table S1). Since entire adult mosquitos contain the gonads, which could be introducing a bias, we also analyzed published RNA-seq data (Papa et al, 2017) of dissected carcasses (soma) and reproductive tissues (germline) (Fig S1B; modEncode and RT-qPCR data from Dmel as a comparison in Fig S1C and D). In this RNA-seq dataset, Agam msls were all classified as not sex-biased. This discrepancy could be a consequence of differences in sensitivity (RT-qPCR versus RNA-seq), or alternatively, other tissues with male-biased expression present in the sample analyzed by us.

Next, we turned our attention to the expression pattern of msls during Anopheles embryogenesis (mixed-sex data from Goltsev et al [2009]). msl-3 and mof mRNA levels drop along development, whereas mle and msl-2 appear continuously expressed (Fig 1F). In Drosophila, msl-2 and mle expression exhibits a different pattern during embryogenesis, with a pronounced induction shortly after zygotic genome activation (ZGA) ([Lott et al, 2011; Samata et al, 2020], Fig S1E).

In Drosophila, both sexes express MOF and MLE proteins (Fig 1G and H [Kotlikova et al, 2006; Conrad et al, 2012]), and at least in some tissues, MSL1 and MSL3 are not only present in males, but also in females (Chlamydas et al, 2016; McCarthy et al, 2019 Preprint). However, posttranscriptional mechanisms prevent MSL2 protein expression in females thereby providing the DC function specifically to males (Beckmann et al, 2005). As antibodies against Agam MSL2 are unavailable, we assessed its protein production by polysome profiling. This revealed that unlike Drosophila, MSL2 protein is produced in both male and female mosquitos (Figs 1I and S1F).

msl-2 KO in A. gambiae triggers sex-independent phenotypes

Our observation that unlike Drosophila, MSL2 is expressed in both male and female mosquitos prompted us to investigate whether it is involved in sex-specific functions or in a more general role in gene expression regulation. We generated transgenic Agam for conditional KO of the Agam msl-2 gene by CRISPR/Cas9 (Fig 2A, see the Materials and Methods section). We crossed vasa-Cas9 mothers with msl-2 gRNA transgenic fathers, where Cas9 protein is maternally deposited into the fertilized egg. This results in msl-2 KO shortly after the onset of zygotic transcription of the paternally provided gRNA transgenes (Fig 2A). Compared to the controls, which completed embryonic development within around 42 h followed by hatching (Goltsev et al, 2007), the msl-2 KO progeny displayed an arrest in early embryogenesis (Figs 2B and S2A). To analyze cell death, we performed TUNEL stainings in embryos collected at 6, 16 and 24 h after egg laying. TUNEL-stained cells could be detected in msl-2 KO embryos around 16 h after egg laying and this became much more pronounced after 24 h of development (Figs 2C and S2B). Compared to Dmel, where MSL2 loss results in male-specific lethality beginning to manifest at the third instar larval stage (after around 5 d of development), this severe early embryonic phenotype affecting both sexes was unexpected. However, it was in line with expression analyses by RT-qPCR (see above), where both msl-2 and msl-3 are highly transcribed at ZGA, which occurs around 2–4 h after egg laying (Figs 2A and S2C). RNA levels of early patterning genes such as hunchback, elovl or eve as well as msl-3 or the X-linked AGAP000634 appeared rather unaffected in msl-2 KO Agam (see below for transcriptome-wide analyses). RT–PCR from these mixed-sex populations collected at different time points after fertilization revealed that there were both males (assessed by expression of Yob [Krzywinska et al, 2016]) and females (assessed by female-spliced Dsx [Kyrou et al, 2018]) among the transcriptionally active population (Fig 2D). This was validated by single-embryo genotyping by qPCR, where we detected approximately equal numbers of males and females in the F1 embryos across different developmental time points (Fig 2E). Therefore, the proportion between the two sexes is unchanged by msl-2 KO throughout embryogenesis until their arrest, indicating that both sexes are equally affected upon loss of MSL2.

Figure 2
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2 msl-2 knockout in Anopheles gambiae triggers sex-independent phenotypes distinct from Drosophila.

(A) Top: schematic representation of the A. gambiae msl-2 gene and position of the four gRNAs that were used to generate gRNA-expressing transgenic mosquitos. Bottom: RT-qPCR analysis of msl-2 RNA levels in F1 progeny of vasa-Cas9 mothers crossed with control (left, grey bars) or msl-2 gRNA (right, blue bars) transgenic fathers. Each bar represents the progeny of one cross collected at the indicated time points after egg laying. (B) as in (A), representative pictures of control or msl-2 gRNA-expressing transgenic F1 progeny embryos collected at the indicated time points after egg laying. The pictures were obtained with a 40× magnification on an inverted microscope. (C) as in (A), representative TUNEL staining and corresponding bright field (BF) images of msl-2 gRNA-expressing transgenic F1 progeny embryos collected at the indicated time points after egg laying. The control embryos did not show a staining (Fig S2B). The pictures were cropped and rotated to orient them properly according to anterior-posterior and dorsal–ventral axis (uncropped pictures in Fig S2B). Scale bar = 20 μm. (D) cropped agarose gel of semi quantitative RT–PCR products obtained from mixed-sex F1 progeny obtained from crosses set as in (A). The RT was primed with gene-specific primers for the male-specifically expressed Yob, the primers for Dsx detect the splice variant exclusively present in females (upper band, arrow) or in both sexes (lower band). ctrl1 and ctrl2 RNAs serve as loading controls (see primer list). (E) Result of individual F1 embryo genotyping (cross as in [A]) as assessed by the expression of male-specific Yob by RT-qPCR. The single embryos were collected at the indicated time points after egg laying. Rps4 was amplified as a control expressed in both males and females (not shown in figure). (F) Cropped immunoblots of extracts obtained from mixed-sex F1 embryos (cross as in [A]). RPB3 serves as a loading control. (G) representative pictures of dissected abdominal tissues of female F1 progeny obtained from crosses of vasa-Cas9 fathers with control (ctrl) or msl-2 gRNA-expressing mothers. The resulting msl-2 KO is induced in the F1 germline tissue, stars point towards the ovaries. Blood feeding triggers egg laying in ctrl, but not msl-2 KO embryos, where ovaries are completely atrophied. Scale bar = 500 μm. (H) Immunofluorescence for H4K16ac (green) in female F1 ovaries obtained from crosses as in (G). The pictures are single z-planes with scale bar = 50 μm (left and middle panel) or 10 μm (right panel) with DAPI shown in blue. Note the gross overall size and difference in tissue morphology. The arrows highlight the grape-shaped germ line tissue (see ctrl), which are highly atrophied (middle) or completely absent (right) in the KO females.

Figure S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S2. Analysis of MSL complex components in male and female Anopheles gambiae.

(A) as in Fig. 2B; representative pictures of control or msl-2 gRNA-expressing transgenic F1 progeny embryos collected at 46 h after egg laying. The pictures were obtained with a 40× magnification on a cell culture microscope. (B) as in Fig. 2C, representative TUNEL staining with pictures of msl-2 gRNA-expressing transgenic F1 progeny embryos collected at the indicated time points after egg laying. The control embryos did not show a staining. Scale bar = 20 μm. (C) RT-qPCR analysis of the RNA levels of the indicated genes in F1 progeny of vasa-Cas9 mothers crossed with control (left, grey bars) or msl-2 gRNA (right, blue bars) transgenic fathers. Each bar represents the progeny of one cross collected at the indicated time points (hours) after egg laying. (D) as in Fig. 2H; representative pictures taken with a stereomicroscope of control or msl-2 gRNA-expressing transgenic male F1 progeny testis (vasa-Cas9 fathers crossed with gRNA-expressing transgenic mothers, scale bar = 320 μm).

In Drosophila, the loss of msl-2 is linked to a bulk decrease of its downstream epigenetic mark, H4K16ac. When comparing Anopheles wild-type and msl-2 KO embryo extracts by Western blot, we found no scorable differences in H4K16ac levels (Fig 2F). In summary, these findings point towards an essential function for MSL2 in both male and female mosquitos early in embryogenesis that is distinct from Drosophila.

We next crossed msl-2 gRNA-expressing mothers with vasa-Cas9 transgenic fathers, which induces a KO in germ line tissues of the F1 progeny. When F1 males were crossed to wild-type females, normal numbers of eggs were laid but those did not hatch, indicating that the F1 males were sterile (pictures of testis in Fig S2D). F1 females displayed a drastic atrophy of the ovaries and failed to lay eggs even after two consecutive blood meals, whereas in the control, a single blood feeding reliably triggered ovary development and egg laying (Fig 2G). This phenotype is again distinct from Drosophila, where msl-2 KO females are fertile (Belote & Lucchesi, 1980). We performed immunofluorescence stainings of H4K16ac in female F1 KO ovaries (Fig 2H; note the different scale bars and overall differences in organ size in control compared with KO). A range in severity of phenotypes was observable. Most females displayed a complete absence of the germ line tissue, where only the wild-type somatic ovary epithelium was present. Some individuals exhibited a defective ovarian tissue that still contained an oocyte. Overall, the tissue architecture and morphology were drastically compromised upon msl-2 KO in the female germ line (Videos 1 and 2). We conclude that MSL2 is essential for faithful progression of embryonic and germ line development in both male and female mosquitos.

Video 1

Z-stacks of Immunofluorescence of control Anopheles gambiae F1 ovary stained for H4K16ac (green) and DAPI (blue).Download video

Video 2

Z-stacks of Immunofluorescence of msl-2 KO Anopheles gambiae F1 ovary stained for H4K16ac (green) and DAPI (blue).Download video

MSL2 confers gene-by-gene regulation on all chromosomes

The early embryonic and pronounced germ line msl-2 KO phenotypes could be a consequence of an additional function of the MSL complex in mosquitos compared with Drosophila, which could mask a phenotype related to regulating the Anopheles X chromosome. To analyze this, we characterized the transcriptome alterations upon loss of msl-2 in Anopheles embryos, in particular with regards to the chromosomal location. We analyzed two different time points: 6–7 h after egg laying with n = 3 biological replicates and 15 h after egg laying with n = 2 biological replicates (Fig S3A; we also generated a single dataset at 24 h after egg laying, which was only considered for qualitative analyses). Because it is not possible to phenotypically separate sexes at this stage, we performed RNA-seq experiments from mixed embryo populations. To validate that this approach per se permits capturing effects on DC of the X, we performed an equivalent experiment from a mixed male and female Dmel embryo population (Fig 3A). We scored differentially expressed (DE) genes using DESeq2 with a false-discovery rate (FDR) cutoff of <0.05 in control versus msl-2 KO (Fig S3B and Tables S2 and S3). As expected from the role of Dmel MSL2 in regulating the X chromosome, we made the following observations in the Drosophila dataset: First, X-linked genes are significantly overrepresented among all down-regulated genes in the msl-2 KO embryos (Fig 3B, 717 of 1,304 genes, Fisher’s exact test P = 1.65645 × 10−273). Second, genes on the X chromosome are globally down-regulated, irrespective of whether a gene is scored as DE or not (Fig 3C, statistical analyses in Table S1). Furthermore, autosomes are moderately up-regulated. Third, the fold change/extent of the down-regulation occurs in a range that is consistent with an approximately twofold, DC-like effect (Fig 3A; log2FC of −0.415 assuming a 50:50 male:female population, where males are affected by twofold and females are unaffected; grey bar). Indeed, the majority of the downregulated genes (1,103 of 1,304) are reduced to 25–75% of control levels. We conclude that X-linked perturbations due to defective DC can be reliably scored from such a mixed-population experiment.

Table S2 List of differentially expressed genes (DESeq2) in Anopheles gambiae msl-2 KO embryos (15 h after egg laying).

Table S3 List of differentially expressed genes (DESeq2) in Drosophila melanogaster msl-2 KO embryos (overnight collection).

Figure S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S3. Transcriptome alterations upon loss of msl-2 in Anopheles gambiae and Drosophila melanogaster.

(A) RNA-seq experiment in D. melanogaster. The heat map illustrates the euclidean distances between all D. melanogaster RNA-seq samples and replicates. (B) RNA-seq experiment in D. melanogaster. Heat maps of the Top75 (by Padj) down-regulated genes in homozygous msl-2Δ/msl-2Δ null mutant compared with the heterozygous control msl-2Δ/CyO-GFP embryos. (C) RNA-seq experiment in A. gambiae. The heat map illustrates the euclidean distances between all A. gambiae RNA-seq samples and replicates. Note that the different time points cluster together, indicating that the effect of transcriptome alterations along embryogenesis are larger than the effect of msl-2 KO. (D) RNA-seq experiment in A. gambiae; Heat maps of all down-regulated (top) and up-regulated (bottom) genes in 6–7 h F1 A. gambiae embryos obtained from vasa-Cas9 mothers crossed with control (referred to as WT) or msl-2 gRNA (gMSL2) transgenic fathers. The D. melanogaster (Dmel) orthologue and function in Drosophila are described. The asterisk (*) denotes genes that are directly adjacent to the msl-2 gene and are likely disrupted as a direct consequence of the gRNAs cutting the locus. NA, no orthologue annotated.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3. Transcriptome alterations upon loss of msl-2 in Drosophila melanogaster and Anopheles gambiae embryos.

(A) RNA-seq experiments were conducted on Drosophila melanogaster F1 embryos obtained from crosses of msl-2Δ/CyO-GFP fathers with msl-2Δ/CyO-GFP mothers. Non-fluorescent msl-2Δ/msl-2Δ null mutant embryos (abbreviated as msl-2) were compared with fluorescent msl-2Δ/CyO-GFP heterozygous controls. MA-Plots showing normalized counts versus the log2(fold change) obtained with DESeq2. Differentially expressed (DE) genes (false-discovery rate < 0.05) are colored in red and the number of DE genes is reported in each panel. The grey-shaded bar indicates the expected log2FC position of down-regulated genes with a dosage compensation-like effect affected in only males but not females in a 50:50 mixed-sex population. See Table S1 for mapping statistics. (B) RNA-seq as in (A), Bar plot for the number of DE genes with respect to their location on the chromosomal arms. (C) RNA-seq as in (A), Density plots of the log2(fold change) obtained with DESeq2 for genes on each of the indicated chromosomal arms. All genes were taken into account for the analyses, irrespective of whether they are scored as DE or not. (D) RNA-seq experiments were conducted on A. gambiae F1 embryos collected at 15 h after egg laying obtained from crosses of vasa-Cas9 mothers with control or msl-2 gRNA-expressing fathers. MA-Plots showing normalized counts versus the log2(fold change) obtained with DESeq2. DE genes (false-discovery rate < 0.05) are colored in red and the number of DE genes is reported in each panel. See Table S1 for mapping statistics. (E) RNA-seq as in (D), Bar plot for the number of DE genes with respect to their location on the chromosomal arms. (F), RNA-seq as in (D), Density plots of the log2(fold change) obtained with DESeq2 for genes on each of the indicated chromosomal arms. All genes were taken into account for the analyses, irrespective of whether they are scored as DE or not.

For Anopheles, we first inspected the earlier dataset around 6–7 h after egg laying, and found few gene expression changes in msl-2 KO embryos (Fig S3C and D). A total of five significantly down-regulated genes (FDR < 0.05) included msl-2 itself, three transcripts directly adjacent to msl-2, and the orthologue of Dmel mira, an actin/myosin binding protein. These limited changes suggest that the aforementioned phenotypes (Fig 2) likely arise from loss of msl-2, rather than from Cas9 off-target effects. Because ZGA occurs after around 3 h of development (Goltsev et al, 2007), this also implies that ZGA and early embryo development occurs faithfully without msl-2. We turned our attention to the 15 h after egg laying time point, where we found 566 down-regulated transcripts, as well as 360 up-regulated transcripts (Fig 3D). Distinct from Drosophila (orange bars in Fig 3B, Fisher’s exact test in Dmel P = 1.65645 × 10−273, see above), the fraction of X-linked, down-regulated genes in Anopheles tends to only reflect the overall proportion of X chromosomal genes in the genome (orange bars in Fig 3E, n = 62 of 566, Fisher’s exact test P = 0.007221). In addition, the extent of the perturbations was much more pronounced compared with Drosophila, as the majority of down-regulated targets (523 of 566) were affected by at least 75% or more (Tables S2 and S3 and Fig 3D). An overall down-regulation of the X chromosome or an up-regulation of autosomes was also not detectable in the Anopheles dataset (Fig 3F and Table S1). We conclude that loss of MSL2 in Anopheles embryos leads to global changes in gene expression affecting all chromosomes (for functional analyses on MSL2-evoked transcriptome alterations not related to X-chromosome–wide DC see below).

Tissue-specific distribution of histone H4 lysine 16 acetylation in A. gambiae

The induction of a subnuclear territory enriched for an epigenetic mark in one of the sexes (H4K16ac in male Dmel, histone H3 lysine 27 trimethylation [H3K27me3] in female mammals, and histone H4 lysine 20 monomethylation [H4K20me1] in hermaphrodite Caenorhabditis elegans) is a characteristic feature of all known DC systems (Wells et al, 2012). To analyze this, we performed immunostainings in adult male and female mosquitos and characterized the tissue-specific enrichment of H4K16ac, in particular with regards to the presence of a X chromosome territory (Fig 4A). H4K16ac exhibited homogeneous staining of the entire nucleus in male midguts (panel 4) and Malpighian tubules (panel 3). The testis (panels 1 and 2) showed a staining pattern where H4K16ac levels were varying across the different cellular stages of sperm differentiation. The early germ cells appeared enriched, wereas the more mature, transcriptionally inert stages (e.g., spermatocytes) were depleted of H4K16ac. This cell type–specific pattern of H4K16ac in the testis was antagonized by the presence of the repressive mark H3K27me3 (Fig S4A). In females, the midguts and Malpighian tubules displayed a similar staining pattern for H4K16ac as in males (Fig 4B, panels 2 and 3; Fig 4C and D for a comparison of Agam to the X chromosome territory in Dmel). The somatic cells of the ovaries showed a high staining intensity and the oocyte itself was H4K16ac-positive. Similar to males, H3K27me3 was homogeneously staining the female nucleus and showed a tissue-specific distribution, where H4K16ac high cells displayed low levels of H3K27me3 and vice versa (Fig S4B). However, despite extensive analyses, we could not identify any cell type in adult Anopheles that showed a nuclear substructure akin to a putative X chromosomal territory that is marked by H4K16ac.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4. Tissue-specific distribution of Histone H4 Lysine 16 acetylation in Anopheles gambiae.

(A) Immunofluorescence for H4K16ac (green) in adult male A. gambiae organs. The areas 1 and 2 highlight the testis, 3 the Malpighian tubules, and 4 the gut. The pictures are orthogonal projections of a z-stack with scale bar = 100 μm. The panels 1–4 on the right are single z-planes, where the highlighted areas were imaged at a higher magnification, scale bar = 20 μm. DAPI is shown in blue. (B) as in (A), but in adult female A. gambiae organs. The areas 1 highlight the ovary, 2 the Malpighian tubules and 3 the gut. The image labelled “zoom” shows a closeup of the ovarian area highlighted with the solid white rectangle (same image cropped and zoomed). (C) close-up of nuclei in Malpighian tubules (polyploid) or guts (diploid) that are highlighted with a solid white rectangle in (A, B). (D) Immunofluorescence for H4K16ac (green) in male and female Drosophila salivary glands (polyploid) and guts (diploid). The pictures are single z-planes with scale bar = 20 μm and DAPI shown in blue. The bottom panel shows a close-up of one nucleus in each males and females (same image zoomed).

Figure S4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S4. Immunofluorescence of H3K27me3 in male and female Anopheles gambiae organs.

(A) Immunofluorescence for H3K27me3 (red) in adult male A. gambiae organs. The pictures are orthogonal projections of a z-stack with scale bar = 20 μm. DAPI is shown in blue, the phalloidin-staining in grey. In mammals, but not A. gambiae, H3K27me3 is enriched on the inactive female X chromosome and forms a territory. (B) as in (A), but in adult female A. gambiae organs.

Genome-wide profiles reveal H4K16ac on active genes in mosquitos

The fact that H4K16ac is not associated with a territory in Anopheles could be a consequence of a less pronounced X chromosome enrichment and hence, immunofluorescence not being sensitive enough to score it. We therefore generated high-resolution MNase-ChIP-seq profiles of H4K16ac, as well as a total H3 control from adult male and female mosquitos. We segmented the genome into 2 kb bins and analyzed the Input normalized ChIP enrichment on each of the chromosomal arms (Fig 5A). Although this approach readily revealed H4K16ac enrichment on the Dmel male X chromosome (Valsecchi et al, 2018), both Anopheles males and females looked rather similar to Drosophila females. Next, we used an unsupervised clustering approach to identify different enrichment patterns at male in comparison with female transcription start sites (TSSs) (Fig 5B). In Drosophila, this approach revealed a Cluster 1, which displays a markedly higher and broad enrichment of H4K16ac in males, but not females. Cluster 1 consisted exclusively of X-chromosomal genes (Fig 5C). Cluster 2 shows a TSS-proximal H4K16ac signal in both males and females. In contrast to Drosophila, the clustering approach did not reveal differences at TSS between male and female mosquitos (Fig 5B) and/or a predominant enrichment of X-linked genes in any of the clusters (Fig 5C). Overall, the H4K16ac enrichment on genes was highly correlated between Agam males and females, but did not show the striking enrichment as on Dmel male X-linked genes (Fig S5A and B; see Fig S5C for Dmel H4K16ac data from different tissues). In line with H4K16ac promoting transcription, the level of H4K16ac on Anopheles genes (low to high by ChIP enrichment in quartiles Q1-Q4) was correlating with low to high read counts at the RNA expression level in RNA-seq (Fig S5D). Agam genes that showed the highest levels of H4K16ac were involved in biological processes broadly involving cellular metabolism and provided molecular functions such as “RNA binding,” “signaling receptor activity” or “molecular transducer activity” (Fig S6A). Our findings were corroborated by visual inspection of the ChIP data in the genome browser (Fig 5D and E). Compared to the H3 control, H4K16ac in Anopheles appeared clearly enriched on individual, active genes on both X and autosomes. This indicates that TSS-proximal H4K16ac promotes gene expression in a similar fashion in Agam males and females and thereby serving a more general role in gene regulation, instead of sex-specific functions.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5. Genome-wide profiles reveal H4K16ac on active genes, but no male-specific X chromosome enrichment.

(A) ChIP-seq profiles for total Histone H3 and H4K16ac in Anopheles gambiae were generated from dissected adult midguts, in Drosophila melanogaster from L3 third instar larvae. Box plots show the mean log2(fold change) ChIP versus Input enrichments per 1 kb bin on each of the chromosomal arms. The biological replicates were merged (also see the Materials and Methods section). (B) as in (A), three unsupervised K-means clusters were generated from the log2(fold change) H4K16ac ChIP versus Input data. The transcription start site of each Agam or Dmel gene served as a reference point, while plotting the mean enrichment profile ± 1 kb. (C) as in (B), the bar plots show the chromosomal location of the genes present in clusters 1–3 compared to the overall number of genes on each of the chromosomal arms. (D) as in (A), genome browser snapshots of the normalized ChIP enrichment in Agam males and females on example regions of the X chromosome (top) and autosome (bottom). (E) as in (A), genome browser snapshots of the normalized ChIP enrichment in Dmel males and females on an example region on the X chromosome.

Figure S5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S5. H4K16ac ChIP-seq analyses in Anopheles gambiae compared with Drosophila melanogaster.

(A) Scatterplot comparing the normalized H4K16ac ChIP-seq (log2FC versus Input) enrichment on each A. gambiae gene in males and females. The bottom plot shows normalized total H3 ChIP-seq (log2FC versus Input) enrichment. The ChIP-seq data were generated from adult mosquitos. The scatterplot is represented as a smoothed color density representation obtained through a (2D) kernel density estimate (smoothScatter function in R). (B) as in (A) but comparing the normalized H4K16ac enrichment (top) or H3K36me3 enrichment (bottom) in D. melanogaster males and females. The data were published in Valsecchi et al (2018). The left plot shows all genes. In the right plot, X-chromosomal genes were excluded. Note the cloud on the bottom right, which show a much higher enrichment in males compared to females and corresponds to X-linked genes. H3K36me3 serves as a control, which does not show any differential pattern between males and females. (C) IGV genome browser snapshots of Drosophila H4K16ac ChIP-seq performed with two different fragmentation methods (sonication versus micrococcal nuclease [MNase] digestion) and in various different stages and genotypes. Data were from S2 cells (Straub et al, 2013), adults (modEncode Consortium et al, 2010), embryos (Samata et al, 2020), and third instar larvae (Valsecchi et al, 2018). (D) Density plots of the RNA expression level of A. gambiae genes grouped by their H4K16ac ChIP-seq enrichment in quartiles 1 (low, blue) to 4 (high, red). RNA expression levels by normalized log10 (counts) in RNA-seq were published by Papa et al (2017). Genes with low H4K16ac levels (Q1 and Q2) are rather lowly expressed, whereas Q3 and Q4 (high H4K16ac) are highly expressed.

Figure S6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S6. Features of transcripts that are altered upon loss of msl-2 in Anopheles gambiae.

(A) Gene Ontology-term analysis of the genes with the highest H4K16ac levels (Q4 quartile) in A. gambiae. The top 15 biological processes (left) and molecular functions (right) are reported with the false-discovery rate, which is color coded with the most significant terms from red to white. (B) RNA-seq as in Fig 3D and n = 926 msl-2 DE genes in A. gambiae embryos (msl-2Δ versus control, false-discovery rate < 0.05) were analyzed for distinctive gene features using the ShinyGo (v0.61) gene-set enrichment tool (Ge et al, 2020).

A. gambiae MSL2 regulates conserved developmental genes

The aforementioned observations indicate that the MSL complex-H4K16ac pathway is unlikely to mediate chromosome-wide DC in Agam. We were therefore interested in which functions and gene-regulatory networks MSL2 might be controlling. To this end, we analyzed our RNA-seq dataset for the enrichment of distinctive features among the DE genes. This revealed that DE genes in msl-2 KO embryos tend to be significantly longer than average in the Anopheles genome (e.g., at the level of coding sequence, transcript length, gene span as well as the number of exons, Figs 6A and S6B). Furthermore, we found a significantly higher fraction of Drosophila orthologues in the msl-2 KO DE group than expected (Fisher’s test P = 6.34 × 10−18) indicating that the Anopheles MSL-regulated genes tend to be conserved among Dipterans (Fig 6B, left panel). There was no difference between the expected and DE gene groups at the sequence level (% identity of the orthologues) (Fig 6B, right panel). We next analyzed the features of the proteins encoded by msl-2 KO DE genes using the STRING database of known and predicted protein–protein interactions (Fig 6C) (Szklarczyk et al, 2018). We find that the proteins encoded by the Top400 DE genes (msl-2 KO versus control) display significantly more protein–protein interactions than expected, indicating that they are highly interactive. If the targets regulated by MSL2 tend to frequently engage in protein–protein interactions and are part of multisubunit complexes, this may explain why their expression levels need to be tightly controlled by chromatin-regulatory pathways such as the H4K16ac–MSL pathway: poorly controlled expression could lead to stoichiometry imbalances, proteotoxicity and aggregation (Brennan et al, 2019).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6. Transcriptome analyses reveal a role for Anopheles MSL2 in regulating conserved developmental-regulatory genes.

(A) RNA-seq as in Fig 3D and n = 926 msl-2 DE genes in A. gambiae embryos (msl-2 KO versus control, false-discovery rate < 0.05, 15 h after egg laying) were analyzed for overrepresentation concerning the numbers of exons versus the expected distribution using the ShinyGo (v0.61) gene-set enrichment tool (Ge et al, 2020). (B) RNA-seq as in Fig 3D, left: msl-2 DE genes or all A. gambiae genes (expected) were analyzed for the presence of an Ensembl BioMart orthologue (black: no orthologue, light blue: low confidence orthologue, blue: high confidence orthologue). The statistical significance was evaluated by comparing underrepresentation in the “no orthologue” group (DE, n = 204 genes without orthologue) versus expected using a Fisher’s exact test. Right: msl-2 DE genes or all A. gambiae genes (expected) were analyzed for % of sequence identity in the target Dmel orthologue. The statistical significance was evaluated by comparing the % sequence identity for all orthologues (n = 971) versus expected using a Wilcoxon test. (C) RNA-seq as in Fig 3D, Protein–protein interaction network and statistical significance obtained on STRING (v11.0) (Szklarczyk et al, 2018) using the Top400 (by Padj) DE genes. (D) RNA-seq as in Fig 3D, enriched functional protein domains and statistical significance obtained from Pfam using the ShinyGo (v0.61) gene-set enrichment tool (Ge et al, 2020). (E) RNA-seq as in Fig 3D, Gene Ontology (GO) analysis of biological processes enrichment of msl-2 DE genes. Enriched GO Terms were visualized with REVIGO. (F) RNA-seq as in Fig 3D, GO analysis of cellular components enrichment of DE genes upon msl-2 KO in A. gambiae. (G) RNA-seq as in Fig 3D, Genome browser snapshots of the normalized read coverage in control (top, black) and msl-2 KO embryos (bottom, blue) at the A. gambie HOX cluster on chromosome 2R. The two tracks represent two different biological replicates.

Analysis of the functional domains of MSL2 targets with the Pfam database revealed an enrichment of “homeobox domain,” which is a well-known module of developmental transcription factors, but also various other domains mediating processes occurring at the cellular surface (Fig 6D). We also performed Gene Ontology term analyses (Fig 6E and F), where we found misregulated biological processes such as “multicellular organism development,” “tissue development,” “tissue morphogenesis,” or “cell differentiation” and cellular components such as “plasma membrane,” “cell periphery,” “extracellular space,” or “cell junction.” A prominent example in the down-regulated genes group was the A. gambiae HOX cluster (Fig 6G). HOX genes encode homeodomain-containing transcription factors that provide specificity and positioning information along the anterior-posterior axis and are conserved master regulators for setting up animal body plans (Gehring et al, 2009). Other misregulated genes include SPZ6 (chromosome 2L), a Spätzle-like cytokine involved in signalling and embryo patterning, or the OSIRIS family members Osi7, Osi9, and Osi14 (chromosome 2R), a conserved dosage-sensitive gene family encoding putative transmembrane receptors. Collectively, this indicates that the MSL complex fine-tunes the expression of conserved genes that participate in highly intertwined developmental-regulatory networks orchestrating Anopheles embryogenesis. A perturbation of this function could provide a likely explanation for the early arrest of msl-2 KO embryos (Fig 2).

Discussion

In this study, we have explored the MSL complex–H4K16ac pathway and its involvement in X chromosome DC in the malaria mosquito A. gambiae. We find similarities, but also notable differences to the related Dipteran D. melanogaster, where these factors have been most extensively characterized. Zygotic msl-2 KO Anopheles embryos fail to develop beyond the segmentation stage. In RNA-seq, we detect misexpression of key developmental-regulatory genes, including the Hox cluster. Genes affected upon msl-2 KO tend to contain more exons, display a higher degree of conservation and participate in more protein–protein interactions than expected. Minor changes in the intertwined gene-regulatory networks operating in a highly precise manner during development can cumulatively lead to severe phenotypes, such as the ones observed by us and therefore explain why they require the fine-tuning by the MSL complex. Indeed, the fine-tuning of developmental-regulatory genes appears to be a common theme for MSL function in insects up to mammals (Valsecchi et al, 2018).

However, the time point of phenotype manifestation at only few hours of embryonic development in Anopheles is very distinct from the effect of MSL loss in Drosophila and mammals. In fruit flies, loss of msls results in lethality occurring at the third instar larval stage (after around 5 d of development) (Belote & Lucchesi, 1980). The early Anopheles phenotype is also very different from the consequences of MSL3 loss in humans, which is associated with global developmental delay post birth affecting both male and female individuals (Basilicata et al, 2018). The reasons for these various time points of phenotype manifestation will need to be carefully dissected in future studies. They could, for example, be related to differences in maternally deposited factors that might allow faithful embryogenesis despite a zygotic loss, the presence of species-specific isoforms and paralogues or differences in stage-specific gene-regulatory networks in a given species.

In Anopheles, MSL2 is expressed in both males and females, and zygotic or germ line KO affects both sexes. Germ line–specific loss of msl-2 results in sterility that is accompanied by severe atrophy of germline tissue. This is different from Drosophila, where MSL2 is male-specifically expressed, zygotic loss of msl-2 is known to only affect males and KO females are fully fertile (Beckmann et al, 2005). If the fine-tuning of developmental genes appears to be the “ancient” function of the MSL complex, how can it be achieved in Drosophila females, if the component that provides the chromatin targeting function to the complex (MSL2) is not present? We note that other MSL proteins, for example, MOF, can bind nucleic acids (Akhtar et al, 2000), making it possible that a subcomplex devoid of MSL2 (Samata et al, 2020) could be taking on these roles. In addition, other protein or ncRNA components could substitute for MSL2. One potential candidate is the protein CLAMP (Soruco et al, 2013), which aids recruitment of the MSL complex and is expressed in both sexes. Last but not least, some studies conversely suggest that MSL2 could indeed be expressed to some degree in females (Cheetham & Brand, 2018). Future studies will have to dissect the stage- and tissue-specific compositions and functions of the MSL complex in Drosophila.

The second intriguing question that arises is by which mechanism the MSL complex became co-opted to specifically target the male X chromosome in drosophilids. MSL2 is a DNA binding protein (Fauth et al, 2010) but the DNA binding CXC domain of MSL2 is dispensable for DC in vivo (Tikhonova et al, 2019). However, MSL2 also binds roX RNAs (Ilik et al, 2013), which are not present in species other than Drosophila (Quinn et al, 2016). Indeed, by virtue of a condensate-forming mechanism, MSL2-roX interaction governs stable association with the X (Valsecchi et al, 2021). The evolution of this species-specific RNA–protein interaction module therefore provides targeting to the X chromosome in Drosophila, but probably not other dipterans.

Another noteworthy aspect of our study is that X chromosome up-regulation in mosquitos is unlikely to be provided by the MSL-H4K16ac pathway. In our RNA-seq, loss of msls in Drosophila is associated with global X chromosome down-regulation. However, an alternative model of inverse DC proposes that loss of msl genes induces autosomal up-regulation, instead of X down-regulation (Birchler, 2016). Indeed, we do observe a moderate up-regulation in our Drosophila data, which would be consistent with this idea. Whereas our datasets do not allow us to unambiguously distinguish between these two models for Drosophila, our Anopheles dataset does not provide any evidence for skewed expression of an entire chromosome, neither the X nor autosomes. In addition, our ChIP-seq profiles show neither evidence of global X chromosome enrichment nor a sexually dimorphic pattern on a gene-by-gene level. Furthermore, H4K16ac does not form the typical X chromosome “territory” in immunofluorescence indicating that the deposition of this mark is not a universal feature for DC in Dipterans. This is reminiscent of another dipteran, Sciara, where mle and H4K16ac appear to play no role in DC (Ruiz et al, 2000). However, different from Anopheles, whose X chromosome is directly comparable to Drosophila, Sciara possesses a peculiar sex chromosome system undergoing imprinted elimination. The anole lizard, a reptile, displays enrichment of H4K16ac on the X chromosome of males compared to females (Marin et al, 2017). A study in monarch butterflies revealed a dichotomy of epigenetic marks on the neoZ chromosome, where both up-regulation (by H4K16ac) and dampening (by H4K20me1) DC mechanisms seem to operate concurrently in the same species (Gu et al, 2019). Yet in both cases, it remains unclear whether H4K16ac is deposited by the MSL complex or another enzyme. Multiple mechanisms have been proposed to be modulating gene-dosage imbalance ranging from transcription, RNA stability as well as protein turnover (Disteche, 2016). In light of the multitude of mechanisms involved in X chromosome regulation in different species, it will be interesting to further explore the epigenome of A. gambiae (Ruiz et al, 2021). Nevertheless, our data suggests that gene content appears not to be the main driver that governs the molecular pathway by which X chromosome DC is achieved. Elegant studies have demonstrated that Dipteran sex chromosomes show a high evolutionary turnover (Vicoso & Bachtrog, 2015). Our data adds on these findings and highlights that the complexity of dosage balancing mechanisms is probably much higher than previously anticipated. This fascinating diversity opens exciting opportunities for future evolutionary and mechanistic studies in Dipteran and other non-model organisms.

Materials and Methods

Mosquito rearing and transgenesis

A. gambiae mosquitos were maintained in standard insectary conditions (26–28°C, 75–80% humidity, and 12-/12-h light/dark cycle). Loss of function of msl-2 (AGAP000553) was achieved by crossing transgenic mosquitos expressing Cas9 to transgenic partners expressing four different guide RNAs targeting exons of msl-2. The 20-nucleotide msl-2–specific sequences in the gRNAs were GATGAGTCCGCTGGGTGGTG, GAGCGTGAGCGTGGAGGAGG, GCGCAAGAGCCTCAGCTGTG, GTGGTAATGTTCACGATGCG. Each gRNA was cloned under the control of the AGAP013557 U6 promoter in a cloning vector compatible with Golden Gate Cloning. The four gRNA modules were then assembled by Golden Gate Cloning in transgenesis plasmid pDSARN (Volohonsky et al, 2015). The full sequence of the resulting plasmid, pDSARN-4xMSL2gRNA can be provided upon request. Transgenic mosquitos were generated by inserting this plasmid in the Agam attP site of docking line X1 (Volohonsky et al, 2015). The resulting transgenic line was called gMSL2 and made homozygous by COPAS sorting of neonate larvae carrying two copies of the DsRedNLS fluorescence transgenesis marker. As crossing partners of gMSL2 mosquitos, two different Cas9 parental lines were tested, both expressing Cas9 inserted in the X1 locus under the control of the germ line vasa promoter (Papathanos et al, 2009): one (YC9) encodes native Cas9 from Streptococcus pyogenes and carries a YFP transgenesis fluorescent marker; the other (eSpC9) is a Cas9 mutant less prone to off-target cutting (Slaymaker et al, 2016) and carries a DsRedNLS fluorescent marker. We initially attempted to generate heritable msl-2 mutations by crossing Cas9-expressing to gRNA-expressing mosquitos to recover mutant progeny as described in Dong et al (2018). However, the lethality and sterility of F1 crosses prevented this approach. Instead, we produced msl-2 knockout embryos by crossing Cas9 females with gRNA males. In this case, maternal deposition of active Cas9 in the oocyte and zygotic expression of gRNAs from the paternal copy of the U6-gRNA transgene resulted in 99% and 95% lethality using YC9 and eSpC9, respectively. The 1–5% escaper larvae presumably resulted from residual MSL2 function in embryos with incomplete mutagenesis or mutations preserving the gene’s function. Most of these larvae died during development. Crossing Cas9 males to gRNA females produced normal F1 progeny, but CRISPR-Cas9 activity in the developing germ line resulted in atrophy of the male and female gonad and sterility of these F1 mosquitos. Note that when using vasa-Cas9 fathers for crosses, we usually observe no significant contribution of paternal Cas9 protein via the sperm of the father.

Polysome profiling

Flash-frozen male and female mosquitos were homogenized in 6 mM Na2HPO4, 4 mM NaH2PO4, 1% NP-40, 150 mM KCl, 5 mM MgCl2, 1 mM EGTA, 50 mM NaF, 0.1 mM Na3VO4, 1 mM PMSF, 1× complete protease inhibitors, and 0.1 mg/ml cycloheximide and the extract was clarified by centrifugation. 140 OD (260 nm) of the extract were layered on a 15–60% sucrose gradient followed by centrifugation at 39,000 rpm for 3 h in a SW40 Rotor. After collecting the fractions, RNA was isolated from the fractions using phenol-chloroform extraction followed by isopropanol precipitation. RNA recovery was determined with a Qubit RNA assay. cDNA was synthesized from 100 ng RNA using the GoScript Reverse Transcription System (Promega) with random primers and subsequently quantified by qRT-PCR. The data were analyzed as described in Keller et al (2012) calculating the RNA enrichment relative to the total amount of RNA in a given fraction.

Protein extracts and Western blots

Extracts from adult Drosophila were prepared by collecting 10 adult male or female flies aged 12–24 h after eclosion, which were decapitated, and heads were homogenized in 50 μl 1× Roti-Load. Anopheles extracts were prepared by bead-beating entire mosquitos in 1× Roti-Load. Samples were incubated for 5 min at 70°C and 10–20 μl used for separation by SDS–PAGE. Proteins were transferred overnight at 60 mA or 1 h 11 min at 111 V to polyvinylidene fluoride (PVDF) membranes using a Bio-Rad Wet Tank Blotting System in Tris-Glycine Transfer-Buffer with 10% methanol. The membranes were blocked for 30 min in 5% milk in PBS-0.1% Tween before incubation with primary antibodies in 0.5% milk PBS-0.1% Tween. Primary antibodies were used at the following dilutions: anti-HA (Mouse, #MMS-101P, 1:5,000; Covance), anti-MOF (Rabbit, 1:2,000 [Raja et al, 2010]), anti-MLE (Rat, 1:2,000 [Ilik et al, 2013]), anti-RNA Pol2 (Mouse clone 4H8, 101,307, 1:10,000; Active Motif), anti-MSL3 (Rat, 1:1,000 [Raja et al, 2010]), anti-MSL1 (Rabbit, 1:3,000 [Raja et al, 2010]), anti-Tubulin (Mouse, 44928, 1:5,000; Abcam), anti-RPB3 (Rabbit, 1:2,000 [Valsecchi et al, 2018]), anti-H4K16ac (Rabbit, 07-329, 1:3,000; Merck Milipore), anti–acetyl-Histone H3 (Rabbit, 06-599; Merck Milipore), anti-H3 (Mouse, 39763 1:5,000; Active Motif), and anti-H4 (61521, 1:1,000; Active Motif). Secondary antibodies (1:10,000) used were anti-mouse IgG HRP (NXA931), anti-rat IgG HRP (NA935V), and anti-rabbit IgG HRP (NA934) from Sigma-Aldrich. Blots were developed using Lumi-Light Western Blotting substrate (12015200001; Roche; Sigma-Aldrich) and imaged on a ChemiDoc XRS+ (Bio-Rad).

RNA expression analysis

Samples were homogenized in TRIzol Reagent (15596026; Thermo Fisher Scientific) followed by RNA purification using a Direct-zol RNA MiniPrep kit or Direct-zol RNA Microprep (R2050 or R2062; Zymo Research) according to the manufacturer’s instruction. cDNA was synthesized with the GoScript Reverse Transcription System (Promega) with random or oligo dT primers as indicated according to the manufacturer’s instructions. For stranded polyA+ mRNA-Seq Library Preparation, the TruSeq stranded mRNA sample preparation kit (RS-122-2101; Illumina) was used.

Quantitative real-time PCR

qPCR was performed on a Roche LightCycler II using FastStart Universal SYBR Green Master (04913914001; Roche) in a 7-μl reaction at 300 nM final primer concentration. Cycling conditions as recommended by the manufacturer were applied. We corrected for primer efficiency using serial dilutions.

Semiquantitative RT–PCR

Semiquantitative RT–PCR was conducted with a one-step RT–PCR kit (210210; QIAGEN), where the reverse transcription with a gene-specific primer and subsequent PCR is conducted in a single reaction. 2.5 ng of RNA template was used in a total reaction mix of 10 μl as of the manufacturers’ instructions. The cycling was conducted as follows: 30 min at 50°C, 15 min at 95°C (1 cycle); 40 s at 94°C, 40 s at 55°C, and 60 s at 72°C (28 cycles); 10 min at 72°C (1 cycle). Products were separated by 1.6% agarose gel electrophoresis in Tris-Borate buffer and stained with SYBR Safe–DNA Gel Stain (S33102; Thermo Fisher Scientific). Gels were scanned using a Typhoon FLA scanner (GE Healthcare Life Sciences).

View this table:
  • View inline
  • View popup
Table 1.

RT-qPCR primers

Drosophila rearing conditions

D. melanogaster were reared on a corn flour–molasses fruit fly medium (1 liter water, 12 g agar–agar threads, 18 g bakery yeast, 10 g soya flour, 80 g corn flour, 22 g molasses, 80 g malt extract, 2.4 g 4-hydroxibenzoic acid methylester [Nipagin], and 6.25 ml propionic acid) at 25°C, 70% relative humidity, and 12-h dark/12-h light cycle. All experiments were conducted at 25°C unless otherwise indicated. White-eyed (w1118) Oregon-R or endogenously tagged msl-2::HA (Valsecchi et al, 2018) was used as a wild-type strain.

Immunofluorescence and TUNEL staining

Tissues were stained according to standard procedures. In brief, tissues were fixed in 4% formaldehyde in PEM (0.1 M PIPES (pH 6.9), 1 mM EGTA, and 1 mM MgCl2) for 20 min and washed three times with PBS. Samples were blocked for 1 h rocking with freshly prepared 0.5% BSA, 0.2% Triton X-100 in 1× PBS solution followed by overnight incubation with primary antibody (anti-H4K16ac [Rabbit, 07-329, 1:200; Merck Millipore], anti-H3K27me3 [Mouse, C15200181-50, 1:300; Diagenode]; staining validated also with anti-H3K27me3 [Rabbit, 9733, 1:250; Cell Signalling]). After Alexa fluorophore–labelled secondary antibody incubation (1:300), samples were thoroughly washed with PT (0.2% Triton X-100 in 1× PBS) before dissection of organs and mounting. Slides were mounted using ProLong Gold Antifade Mountant with DAPI (P36935; Thermo Fisher Scientific) and imaged using a Zeiss LSM880 Airyscan Microscope. For the TUNEL staining, embryos were dechorionated before fixation in 4% PFA. Embryos were washed three times for 5 min in PBS + 0.5% Triton X-100, followed by another three washes in PBS + 0.2% Tween-20. This was followed by TUNEL staining with the In Situ Cell Death Detection Kit, TMR red (12156792910; Sigma-Aldrich) following the procedures described in the manual. The embryos were mounted in ProLong Diamond Antifade Mountant (P36961; Thermo Fisher Scientific) and images obtained using a fluorescence spinning disc confocal microscope, VisiScope 5 Elements (Visitron Systems GmbH), which is based on a Ti-2E (Nikon) stand and equipped with a spinning disc unit (CSU-W1, 50 μm pinhole; Yokogawa). The set-up was controlled by the VisiView 5.0 software and images were acquired with a 20× air objective (CSI Plan Apo Lambda, NA 0.75; Nikon) and a sCMOS camera (BSI; Photometrics). 3D stacks of images were recorded for each sample, the exposure time was set to 200 ms for each channel. Maximum intensity projections of z-stacks and image processing were performed in Fiji (version 2.0.0-rc-68/1.52 h). The immunofluorescence stainings have been replicated in n independent experiments: H4K16ac (n = 4 for both males and females), H3K27me3 (n = 4 males, n = 2 females for Cell Signalling Antibody; n = 2 for both males and females for mouse Diagenode antibody), TUNEL (n = 2).

ChIP

ChIP was performed as described in Valsecchi et al (2018). In brief, midguts of A. gambiae were dissected in PBS and fixed in 0.2% formaldehyde for 15 min at RT. The fixation was stopped by addition of Glycine and the fixative was removed. Midguts were homogenized in buffer G1 (5 mM Hepes [pH 7.4], 20 mM KCl, 10 mM MgCl2, 0.5% Tween-20, and 20% Glycerol containing 1 × Complete Protease Inhibitors) at 4°C. The homogenate was cleared through a sucrose cushion (G1 + 30% sucrose). Nuclei were washed once in (10 mM Tris [pH 8.0], 5 mM MgCl2, 1 mM CaCl2, 10 mM NaCl, 0.1% IGEPAL CA-630, and 0.25 M Sucrose) and then resuspended in Aline buffer (20 mM Tris [pH 8.0], 5 mM MgCl2, 1 mM CaCl2, 10 mM NaCl, 1% Triton X-100, and 0.25 M Sucrose). This was followed by MNase digestion at 25°C for 15 min typically using 0.5–1 μl Micrococcal Nuclease (M0247; New England Biolabs, amounts optimized depending on the batch). 10× High Salt Buffer (200 mM EDTA and 4 M NaCl) was added to stop the MNase reaction before treatment in a Bioruptor Pico (10 cycles, 30 s ON/OFF). After clarification by centrifugation (10 min at 12,000g), the supernatant was used for immunoprecipitation. 1 μl of anti–acetyl-Histone H4 (Lys16) Antibody, 07-329; Merck Millipore or 0.5 μl of Histone H3 antibody (mAb), MABI 0301 (39763; Active Motif) were used for immunoprecipitations followed by washes, elution by reverse cross-linking, RNase A, and Proteinase K treatment according to standard procedures. DNA was purified by phenol–chloroform extraction followed by ethanol precipitation. Libraries were made with a NEBNext UltraTM II DNA Library Prep Kit for Illumina (E7645; New England Biolabs) according to the manufacturer’s instructions.

Processing of RNA-seq and ChIP-seq datasets

For quality control, we used FastQC version: v0.11.3. For RNA-seq, paired-end reads were mapped with RNA STAR (Galaxy version 2.5.2b-0) to AgamP4_ensembl or dm6_ensembl with default parameters. Mapping statistics are reported in Table S1. The raw read counts per gene were obtained with featureCounts (1.4.6.p5). This was followed by differential expression analysis with DESeq2, R version 3.0.2. DESeq2 normalizes the read counts with a linear model, where counts for genes are modeled using a negative binomial distribution with fitted mean and a gene-specific dispersion parameter (details in Love et al [2014]). Table S1 reports statistical analyses for the Drosophila RNA-seq dataset for both the conventional DESeq2 standard normalization as well as where the data are normalized to the rRNA reads (note that rRNA species are typically present in RNA-seq data despite the polyA enrichment in the kit for mRNA-Seq). For the latter analysis, we also used DESeq2 and specified the rRNA genes read counts to be used for size factor estimation. For ChIP-seq, paired-end reads were mapped to AgamP4_ensembl with bowtie2 (Galaxy version v2.2.6.2) using default parameters. To generate coverage box plots for H4K16ac or H3 ChIP on each of the chromosomal arms, we used deepTools bamCompare v2.4.2 to calculate the readCount normalized log2 ratio versus the Input per 1 kb bin. To generate readCount normalized coverage bigwig files, we calculated the log2 ratio versus the Input using a bin size of 5 bp with deepTools bamCompare v2.4.2. Duplicate reads and reads with a mapping quality <10 were removed, the X chromosome was ignored for scaling. The bigwig files were then used to generate k-means clusters with deepTools plotHeatmap v2.4.2, calculation of the mean enrichment for each gene using multiBigwigSummary v2.4.2 and analysis with the IGV genome browser. GOterm analyses were conducted with the PANTHER Overrepresentation Test and Fisher analyses with Bonferroni correction on http://geneontology.org/ and visualized using REVIGO (http://revigo.irb.hr/r). Circos plots visualizing the homology between A. gambiae and D. melanogaster chromosomes were obtained with SyMAP Version v4.2 (default parameters). Previously published data were Drosophila H4K16ac ChIP-seq from (modEncode Consortium et al, 2010; Straub et al, 2013; Valsecchi et al, 2018; Samata et al, 2020), A. gambiae RNA expression data from Papa et al (2017) and Goltsev et al (2009). D. melanogaster modEncode RNA-seq data were obtained from FlyBase.

Statistics

All statistics were calculated with R Studio (v1.3.959). In the box plots the line that divides the box into two parts represents the median, box bottom, and top edges represent interquartile ranges (IQR, 0.25th to 0.75th quartile [Q1-Q3]), whiskers represent Q1 − 1.5*IQR (bottom), Q3 + 1.5*IQR (top). Bar plots represent the mean ± SEM with overlaid data points representing independent experiments. Results were considered significant at FDR (P-value) below 0.05. ns, not significant; NA, not analyzed; SEM, standard error of the mean. Unless otherwise indicated statistical tests were two-sided.

Protein sequences and alignments

Protein sequences were retrieved from FlyBase (http://flybase.org) and VectorBase (https://www.vectorbase.org/). Protein alignments were created on Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clustalo/). Results of BLAST searches (BLASTP, 2.2.29+) are reported in Table S1. Lists of 1:1 orthologues were obtained via the Biomart tool on VectorBase.

Bioinformatic and web resources

Bowtie2 (https://github.com/BenLangmead/bowtie2) (Langmead & Salzberg, 2012); deepTools2 (https://deeptools.readthedocs.io/en/latest/) (Ramírez et al, 2016); Galaxy (https://github.com/bgruening/galaxytools) (Grüning et al, 2017); IGV (https://software.broadinstitute.org/software/igv/) (Robinson et al, 2011); R (https://www.r-project.org); DESeq2 (http://bioconductor.org/packages/DESeq2/); GoTerm analysis (http://geneontology.org/) (The Gene Ontology Consortium, 2019); REVIGO (http://revigo.irb.hr/r) (Supek et al, 2011); STRING (https://string-db.org/) (Szklarczyk et al, 2018); ShinyGO (http://bioinformatics.sdstate.edu/go/) (Ge et al, 2020).

Data Availability

RNA-seq and ChIP-seq data have been deposited to the Gene Expression Omnibus under the accession number GSE153780. All other relevant data supporting the key findings of this study are available within the article and its supplementary files or from the corresponding authors upon request.

Acknowledgements

We thank the Max Planck Institute of Immunobiology and Epigenetics facilities and the IMB Microscopy Core Facility for support. We thank Marc Bayer and Sarah Toscano for contributing initial analyses to set up the project, as well as Devon Ryan and Gina Renschler for help with data processing and bioinformatic analyses. We thank Gina Renschler, Maria Shvedunova, and Agata Kalita for their input on the manuscript. CI Keller Valsecchi was supported by a Human Frontier Science Program long-term fellowship (000233/2014-L). This study was supported by the German Research Foundation (DFG) under Germany’s Excellence Strategy (CIBSS–EXC-2189–Project ID 390939984). Mosquito production and insectarium operation were supported by Agence Nationale de la Recherche grant #ANR-11-EQPX-0022. This work was supported by the DFG under the CRC 992 (project A02) and CRC 1381 (project B3) awarded to A Akhtar.

Author Contributions

  • CI Keller Valsecchi: conceptualization, data curation, formal analysis, supervision, funding acquisition, validation, investigation, visualization, methodology, project administration, and writing—original draft, review, and editing.

  • E Marois: conceptualization, resources, funding acquisition, validation, investigation, visualization, methodology, and writing—review and editing.

  • MF Basilicata: formal analysis, supervision, validation, investigation, visualization, methodology, project administration, and writing—review and editing.

  • P Georgiev: conceptualization, resources, investigation, methodology, project administration, and writing—review and editing.

  • A Akhtar: conceptualization, resources, supervision, funding acquisition, methodology, project administration, and writing—review and editing.

Conflict of Interest Statement

The authors declare that they have no conflict of interest.

  • Received December 16, 2020.
  • Revision received June 27, 2021.
  • Accepted June 28, 2021.
  • © 2021 Keller Valsecchi et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Akhtar A,
    2. Zink D,
    3. Becker PB
    (2000) Chromodomains are protein-RNA interaction modules. Nature 407: 405–409. doi:10.1038/35030169
    OpenUrlCrossRefPubMed
  2. ↵
    1. Alekseyenko AA,
    2. Ellison CE,
    3. Gorchakov AA,
    4. Zhou Q,
    5. Kaiser VB,
    6. Toda N,
    7. Walton Z,
    8. Peng S,
    9. Park PJ,
    10. Bachtrog D, et al.
    (2013) Conservation and de novo acquisition of dosage compensation on newly evolved sex chromosomes in Drosophila. Genes Dev 27: 853–858. doi:10.1101/gad.215426.113
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Bashaw GJ,
    2. Baker BS
    (1995) The msl-2 dosage compensation gene of Drosophila encodes a putative DNA-binding protein whose expression is sex specifically regulated by Sex-lethal. Development 121: 3245–3258. doi:10.1242/dev.121.10.3245
    OpenUrlAbstract
  4. ↵
    1. Basilicata MF,
    2. Bruel AL,
    3. Semplicio G,
    4. Valsecchi CIK,
    5. Aktaş T,
    6. Duffourd Y,
    7. Rumpf T,
    8. Morton J,
    9. Bache I,
    10. Szymanski WG, et al.
    (2018) De novo mutations in MSL3 cause an X-linked syndrome marked by impaired histone H4 lysine 16 acetylation. Nat Genet 50: 1442–1451. doi:10.1038/s41588-018-0220-y
    OpenUrlCrossRef
  5. ↵
    1. Beckmann K,
    2. Grskovic M,
    3. Gebauer F,
    4. Hentze MW
    (2005) A dual inhibitory mechanism restricts msl-2 mRNA translation for dosage compensation in Drosophila. Cell 122: 529–540. doi:10.1016/j.cell.2005.06.011
    OpenUrlCrossRefPubMed
  6. ↵
    1. Belote JM,
    2. Lucchesi JC
    (1980) Male-specific lethal mutations of Drosophila melanogaster. Genetics 96: 165–186. doi:10.1093/genetics/96.1.165
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Birchler JA
    (2016) Parallel universes for models of X chromosome dosage compensation in Drosophila: A review. Cytogenet Genome Res 148: 52–67. doi:10.1159/000445924
    OpenUrlCrossRef
  8. ↵
    1. Brennan CM,
    2. Vaites LP,
    3. Wells JN,
    4. Santaguida S,
    5. Paulo JA,
    6. Storchova Z,
    7. Harper JW,
    8. Marsh JA,
    9. Amon A
    (2019) Protein aggregation mediates stoichiometry of protein complexes in aneuploid cells. Genes Dev 33: 1031–1047. doi:10.1101/gad.327494.119
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Cheetham SW,
    2. Brand AH
    (2018) RNA-DamID reveals cell-type-specific binding of roX RNAs at chromatin-entry sites. Nat Struct Mol Biol 25: 109–114. doi:10.1038/s41594-017-0006-4
    OpenUrlCrossRef
  10. ↵
    1. Chlamydas S,
    2. Holz H,
    3. Samata M,
    4. Chelmicki T,
    5. Georgiev P,
    6. Pelechano V,
    7. Dündar F,
    8. Dasmeh P,
    9. Mittler G,
    10. Cadete FT, et al.
    (2016) Functional interplay between MSL1 and CDK7 controls RNA polymerase II Ser5 phosphorylation. Nat Struct Mol Biol 23: 580–589. doi:10.1038/nsmb.3233
    OpenUrlCrossRef
  11. ↵
    1. Conrad T,
    2. Cavalli FM,
    3. Holz H,
    4. Hallacli E,
    5. Kind J,
    6. Ilik I,
    7. Vaquerizas JM,
    8. Luscombe NM,
    9. Akhtar A
    (2012) The MOF chromobarrel domain controls genome-wide H4K16 acetylation and spreading of the MSL complex. Dev Cell 22: 610–624. doi:10.1016/j.devcel.2011.12.016
    OpenUrlCrossRefPubMed
  12. ↵
    1. Deitz KC,
    2. Takken W,
    3. Slotman MA
    (2018) The effect of hybridization on dosage compensation in member species of the Anopheles gambiae species complex. Genome Biol Evol 10: 1663–1672. doi:10.1093/gbe/evy108
    OpenUrlCrossRef
  13. ↵
    1. Disteche CM
    (2016) Dosage compensation of the sex chromosomes and autosomes. Semin Cell Dev Biol 56: 9–18. doi:10.1016/j.semcdb.2016.04.013
    OpenUrlCrossRefPubMed
  14. ↵
    1. Dong Y,
    2. Simões ML,
    3. Marois E,
    4. Dimopoulos G
    (2018) CRISPR/Cas9 -mediated gene knockout of Anopheles gambiae FREP1 suppresses malaria parasite infection. PLoS Pathog 14: e1006898. doi:10.1371/journal.ppat.1006898
    OpenUrlCrossRefPubMed
  15. ↵
    1. Fauth T,
    2. Müller-Planitz F,
    3. König C,
    4. Straub T,
    5. Becker PB
    (2010) The DNA binding CXC domain of MSL2 is required for faithful targeting the Dosage Compensation Complex to the X chromosome. Nucleic Acids Res 38: 3209–3221. doi:10.1093/nar/gkq026
    OpenUrlCrossRefPubMed
  16. ↵
    1. Ge SX,
    2. Jung D,
    3. Yao R
    (2020) ShinyGO: A graphical gene-set enrichment tool for animals and plants. Bioinformatics 36: 2628–2629. doi:10.1093/bioinformatics/btz931
    OpenUrlCrossRefPubMed
  17. ↵
    1. Gehring WJ,
    2. Kloter U,
    3. Suga H
    (2009) Evolution of the Hox gene complex from an evolutionary ground state. Curr Top Dev Biol 88: 35–61. doi:10.1016/S0070-2153(09)88002-2
    OpenUrlCrossRefPubMed
  18. ↵
    1. Goltsev Y,
    2. Fuse N,
    3. Frasch M,
    4. Zinzen RP,
    5. Lanzaro G,
    6. Levine M
    (2007) Evolution of the dorsal-ventral patterning network in the mosquito, Anopheles gambiae. Development 134: 2415–2424. doi:10.1242/dev.02863
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Goltsev Y,
    2. Rezende GL,
    3. Vranizan K,
    4. Lanzaro G,
    5. Valle D,
    6. Levine M
    (2009) Developmental and evolutionary basis for drought tolerance of the Anopheles gambiae embryo. Dev Biol 330: 462–470. doi:10.1016/j.ydbio.2009.02.038
    OpenUrlCrossRefPubMed
  20. ↵
    1. Grüning BA,
    2. Fallmann J,
    3. Yusuf D,
    4. Will S,
    5. Erxleben A,
    6. Eggenhofer F,
    7. Houwaart T,
    8. Batut B,
    9. Videm P,
    10. Bagnacani A, et al.
    (2017) The RNA workbench: Best practices for RNA and high-throughput sequencing bioinformatics in Galaxy. Nucleic Acids Res 45: W560–W566. doi:10.1093/nar/gkx409
    OpenUrlCrossRef
  21. ↵
    1. Gu L,
    2. Reilly PF,
    3. Lewis JJ,
    4. Reed RD,
    5. Andolfatto P,
    6. Walters JR
    (2019) Dichotomy of dosage compensation along the neo Z chromosome of the monarch butterfly. Curr Biol 29: 4071–4077.e3. doi:10.1016/j.cub.2019.09.056
    OpenUrlCrossRef
  22. ↵
    1. Hallacli E,
    2. Lipp M,
    3. Georgiev P,
    4. Spielman C,
    5. Cusack S,
    6. Akhtar A,
    7. Kadlec J
    (2012) Msl1-mediated dimerization of the dosage compensation complex is essential for male X-chromosome regulation in Drosophila. Mol Cell 48: 587–600. doi:10.1016/j.molcel.2012.09.014
    OpenUrlCrossRefPubMed
  23. ↵
    1. Herpin A,
    2. Schartl M
    (2015) Plasticity of gene-regulatory networks controlling sex determination: Of masters, slaves, usual suspects, newcomers, and usurpators. EMBO Rep 16: 1260–1274. doi:10.15252/embr.201540667
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Ilik IA,
    2. Quinn JJ,
    3. Georgiev P,
    4. Tavares-Cadete F,
    5. Maticzka D,
    6. Toscano S,
    7. Wan Y,
    8. Spitale RC,
    9. Luscombe N,
    10. Backofen R, et al.
    (2013) Tandem stem-loops in roX RNAs act together to mediate X chromosome dosage compensation in Drosophila. Mol Cell 51: 156–173. doi:10.1016/j.molcel.2013.07.001
    OpenUrlCrossRefPubMed
  25. ↵
    1. Keller C,
    2. Adaixo R,
    3. Stunnenberg R,
    4. Woolcock KJ,
    5. Hiller S,
    6. Bühler M
    (2012) HP1(Swi6) mediates the recognition and destruction of heterochromatic RNA transcripts. Mol Cell 47: 215–227. doi:10.1016/j.molcel.2012.05.009
    OpenUrlCrossRefPubMed
  26. ↵
    1. Keller CI,
    2. Akhtar A
    (2015) The MSL complex: Juggling RNA-protein interactions for dosage compensation and beyond. Curr Opin Genet Dev 31: 1–11. doi:10.1016/j.gde.2015.03.007
    OpenUrlCrossRefPubMed
  27. ↵
    1. Kotlikova IV,
    2. Demakova OV,
    3. Semeshin VF,
    4. Shloma VV,
    5. Boldyreva LV,
    6. Kuroda MI,
    7. Zhimulev IF
    (2006) The Drosophila dosage compensation complex binds to polytene chromosomes independently of developmental changes in transcription. Genetics 172: 963–974. doi:10.1534/genetics.105.045286
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Krzywinska E,
    2. Dennison NJ,
    3. Lycett GJ,
    4. Krzywinski J
    (2016) A maleness gene in the malaria mosquito Anopheles gambiae. Science 353: 67–69. doi:10.1126/science.aaf5605
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Kyrou K,
    2. Hammond AM,
    3. Galizi R,
    4. Kranjc N,
    5. Burt A,
    6. Beaghton AK,
    7. Nolan T,
    8. Crisanti A
    (2018) A CRISPR-Cas9 gene drive targeting doublesex causes complete population suppression in caged Anopheles gambiae mosquitoes. Nat Biotechnol 36: 1062–1066. doi:10.1038/nbt.4245
    OpenUrlCrossRefPubMed
  30. ↵
    1. Landeen EL,
    2. Presgraves DC
    (2013) Evolution: From autosomes to sex chromosomes--and back. Curr Biol 23: R848–R850. doi:10.1016/j.cub.2013.08.021
    OpenUrlCrossRef
  31. ↵
    1. Langmead B,
    2. Salzberg SL
    (2012) Fast gapped-read alignment with Bowtie 2. Nat Methods 9: 357–359. doi:10.1038/nmeth.1923
    OpenUrlCrossRefPubMed
  32. ↵
    1. Lott SE,
    2. Villalta JE,
    3. Schroth GP,
    4. Luo S,
    5. Tonkin LA,
    6. Eisen MB
    (2011) Noncanonical compensation of zygotic X transcription in early Drosophila melanogaster development revealed through single-embryo RNA-seq. PLoS Biol 9: e1000590. doi:10.1371/journal.pbio.1000590
    OpenUrlCrossRefPubMed
  33. ↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15: 550. doi:10.1186/s13059-014-0550-8
    OpenUrlCrossRefPubMed
  34. ↵
    1. Lucchesi JC,
    2. Kuroda MI
    (2015) Dosage compensation in Drosophila. Cold Spring Harb Perspect Biol 7: a019398. doi:10.1101/cshperspect.a019398
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Marin R,
    2. Cortez D,
    3. Lamanna F,
    4. Pradeepa MM,
    5. Leushkin E,
    6. Julien P,
    7. Liechti A,
    8. Halbert J,
    9. Brüning T,
    10. Mössinger K, et al.
    (2017) Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage. Genome Res 27: 1974–1987. doi:10.1101/gr.223727.117
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. McCarthy A,
    2. Sarkar K,
    3. Martin ET,
    4. Upadhyay M,
    5. James JR,
    6. Lin JM,
    7. Jang S,
    8. Williams ND,
    9. Forni PE,
    10. Buszczak M, et al.
    (2019) MSL3 coordinates a transcriptional and translational meiotic program in female Drosophila. BioRxiv doi:10.1101/2019.12.18.879874 (Preprint posted December 19, 2019).
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Roy S,
    2. Ernst J,
    3. Kharchenko PV,
    4. Kheradpour P,
    5. Negre N,
    6. Eaton ML,
    7. Landolin JM,
    8. Bristow CA,
    9. Ma L,
    10. Lin MF, et al., modEncode Consortium
    , (2010) Identification of functional elements and regulatory circuits by Drosophila modENCODE. Science 330: 1787–1797. doi:10.1126/science.1198374
    OpenUrlCrossRefPubMed
  38. ↵
    1. Papa F,
    2. Windbichler N,
    3. Waterhouse RM,
    4. Cagnetti A,
    5. D’Amato R,
    6. Persampieri T,
    7. Lawniczak MKN,
    8. Nolan T,
    9. Papathanos PA
    (2017) Rapid evolution of female-biased genes among four species of Anopheles malaria mosquitoes. Genome Res 27: 1536–1548. doi:10.1101/gr.217216.116
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Papathanos PA,
    2. Windbichler N,
    3. Menichelli M,
    4. Burt A,
    5. Crisanti A
    (2009) The vasa regulatory region mediates germline expression and maternal transmission of proteins in the malaria mosquito Anopheles gambiae: A versatile tool for genetic control strategies. BMC Mol Biol 10: 65. doi:10.1186/1471-2199-10-65
    OpenUrlCrossRefPubMed
  40. ↵
    1. Quinn JJ,
    2. Zhang QC,
    3. Georgiev P,
    4. Ilik IA,
    5. Akhtar A,
    6. Chang HY
    (2016) Rapid evolutionary turnover underlies conserved lncRNA-genome interactions. Genes Dev 30: 191–207. doi:10.1101/gad.272187.115
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Raja SJ,
    2. Charapitsa I,
    3. Conrad T,
    4. Vaquerizas JM,
    5. Gebhardt P,
    6. Holz H,
    7. Kadlec J,
    8. Fraterman S,
    9. Luscombe NM,
    10. Akhtar A
    (2010) The nonspecific lethal complex is a transcriptional regulator in Drosophila. Mol Cell 38: 827–841. doi:10.1016/j.molcel.2010.05.021
    OpenUrlCrossRefPubMed
  42. ↵
    1. Ramírez F,
    2. Ryan DP,
    3. Grüning B,
    4. Bhardwaj V,
    5. Kilpert F,
    6. Richter AS,
    7. Heyne S,
    8. Dündar F,
    9. Manke T
    (2016) deepTools2: A next generation web server for deep-sequencing data analysis. Nucleic Acids Res 44: W160–W165. doi:10.1093/nar/gkw257
    OpenUrlCrossRefPubMed
  43. ↵
    1. Renschler G,
    2. Richard G,
    3. Valsecchi CK,
    4. Toscano S,
    5. Arrigoni L,
    6. Ramírez F,
    7. Akhtar A
    (2019) Hi-C guided assemblies reveal conserved regulatory topologies on X and autosomes despite extensive genome shuffling. Genes Dev 33: 1591–1612. doi:10.1101/gad.328971.119
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Robinson JT,
    2. Thorvaldsdóttir H,
    3. Winckler W,
    4. Guttman M,
    5. Lander ES,
    6. Getz G,
    7. Mesirov JP
    (2011) Integrative genomics viewer. Nat Biotechnol 29: 24–26. doi:10.1038/nbt.1754
    OpenUrlCrossRefPubMed
  45. ↵
    1. Rose G,
    2. Krzywinska E,
    3. Kim J,
    4. Revuelta L,
    5. Ferretti L,
    6. Krzywinski J
    (2016) Dosage compensation in the african malaria mosquito Anopheles gambiae. Genome Biol Evol 8: 411–425. doi:10.1093/gbe/evw004
    OpenUrlCrossRefPubMed
  46. ↵
    1. Ruiz JL,
    2. Ranford-Cartwright LC,
    3. Gómez-Díaz E
    (2020) The regulatory genome of the malaria vector Anopheles gambiae: Integrating chromatin accessibility and gene expression. NAR Genom Bioinform 3: lqaa113. doi:10.1093/nargab/lqaa113
    OpenUrlCrossRef
  47. ↵
    1. Ruiz MF,
    2. Esteban MR,
    3. Doñoro C,
    4. Goday C,
    5. Sánchez L
    (2000) Evolution of dosage compensation in Diptera: The gene maleless implements dosage compensation in Drosophila (Brachycera suborder) but its homolog in Sciara (Nematocera suborder) appears to play no role in dosage compensation. Genetics 156: 1853–1865. doi:10.1093/genetics/156.4.1853
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Samata M,
    2. Akhtar A
    (2018) Dosage compensation of the X chromosome: A complex epigenetic assignment involving chromatin regulators and long noncoding RNAs. Annu Rev Biochem 87: 323–350. doi:10.1146/annurev-biochem-062917-011816
    OpenUrlCrossRef
  49. ↵
    1. Samata M,
    2. Alexiadis A,
    3. Richard G,
    4. Georgiev P,
    5. Nuebler J,
    6. Kulkarni T,
    7. Renschler G,
    8. Basilicata MF,
    9. Zenk FL,
    10. Shvedunova M, et al.
    (2020) Intergenerationally maintained histone H4 lysine 16 acetylation is instructive for future gene activation. Cell 82: 127–144.e23. doi:10.1016/j.cell.2020.05.026
    OpenUrlCrossRef
  50. ↵
    1. Slaymaker IM,
    2. Gao L,
    3. Zetsche B,
    4. Scott DA,
    5. Yan WX,
    6. Zhang F
    (2016) Rationally engineered Cas9 nucleases with improved specificity. Science 351: 84–88. doi:10.1126/science.aad5227
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Soruco MM,
    2. Chery J,
    3. Bishop EP,
    4. Siggers T,
    5. Tolstorukov MY,
    6. Leydon AR,
    7. Sugden AU,
    8. Goebel K,
    9. Feng J,
    10. Xia P, et al.
    (2013) The CLAMP protein links the MSL complex to the X chromosome during Drosophila dosage compensation. Genes Dev 27: 1551–1556. doi:10.1101/gad.214585.113
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Straub T,
    2. Zabel A,
    3. Gilfillan GD,
    4. Feller C,
    5. Becker PB
    (2013) Different chromatin interfaces of the Drosophila dosage compensation complex revealed by high-shear ChIP-seq. Genome Res 23: 473–485. doi:10.1101/gr.146407.112
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Supek F,
    2. Bošnjak M,
    3. Škunca N,
    4. Šmuc T
    (2011) REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One 6: e21800. doi:10.1371/journal.pone.0021800
    OpenUrlCrossRefPubMed
  54. ↵
    1. Szklarczyk D,
    2. Gable AL,
    3. Lyon D,
    4. Junge A,
    5. Wyder S,
    6. Huerta-Cepas J,
    7. Simonovic M,
    8. Doncheva NT,
    9. Morris JH,
    10. Bork P, et al.
    (2018) STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res 47: D607–D613. doi:10.1093/nar/gky1131
    OpenUrlCrossRef
  55. ↵
    1. The Gene Ontology Consortium
    (2019) The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res 47: D330–D338. doi:10.1093/nar/gky1055
    OpenUrlCrossRefPubMed
  56. ↵
    1. Tikhonova E,
    2. Fedotova A,
    3. Bonchuk A,
    4. Mogila V,
    5. Larschan EN,
    6. Georgiev P,
    7. Maksimenko O
    (2019) The simultaneous interaction of MSL2 with CLAMP and DNA provides redundancy in the initiation of dosage compensation in Drosophila males. Development 146: dev179663. doi:10.1242/dev.179663
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Valsecchi CIK,
    2. Basilicata MF,
    3. Georgiev P,
    4. Gaub A,
    5. Seyfferth J,
    6. Kulkarni T,
    7. Panhale A,
    8. Semplicio G,
    9. Manjunath V,
    10. Holz H, et al.
    (2021) RNA nucleation by MSL2 induces selective X chromosome compartmentalization. Nature 589: 137–142. doi:10.1038/s41586-020-2935-z
    OpenUrlCrossRef
  58. ↵
    1. Valsecchi CIK,
    2. Basilicata MF,
    3. Semplicio G,
    4. Georgiev P,
    5. Gutierrez NM,
    6. Akhtar A
    (2018) Facultative dosage compensation of developmental genes on autosomes in Drosophila and mouse embryonic stem cells. Nat Commun 9: 3626. doi:10.1038/s41467-018-05642-2
    OpenUrlCrossRef
  59. ↵
    1. Vicoso B,
    2. Bachtrog D
    (2015) Numerous transitions of sex chromosomes in Diptera. PLoS Biol 13: e1002078. doi:10.1371/journal.pbio.1002078
    OpenUrlCrossRefPubMed
  60. ↵
    1. Volohonsky G,
    2. Terenzi O,
    3. Soichot J,
    4. Naujoks DA,
    5. Nolan T,
    6. Windbichler N,
    7. Kapps D,
    8. Smidler AL,
    9. Vittu A,
    10. Costa G, et al.
    (2015) Tools for Anopheles gambiae transgenesis. G3 (Bethesda) 5: 1151–1163. doi:10.1534/g3.115.016808
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Wells MB,
    2. Csankovszki G,
    3. Custer LM
    (2012) Finding a balance: How diverse dosage compensation strategies modify histone h4 to regulate transcription. Genet Res Int 2012: 795069. doi:10.1155/2012/795069
    OpenUrlCrossRefPubMed
  62. ↵
    1. Zdobnov EM,
    2. von Mering C,
    3. Letunic I,
    4. Torrents D,
    5. Suyama M,
    6. Copley RR,
    7. Christophides GK,
    8. Thomasova D,
    9. Holt RA,
    10. Subramanian GM, et al.
    (2002) Comparative genome and proteome analysis of Anopheles gambiae and Drosophila melanogaster. Science 298: 149–159. doi:10.1126/science.1077061
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Distinct mechanisms mediate X chromosome dosage compensation in Anopheles and Drosophila
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
X chromosome dosage compensation in Anopheles
Claudia Isabelle Keller Valsecchi, Eric Marois, M Felicia Basilicata, Plamen Georgiev, Asifa Akhtar
Life Science Alliance Jul 2021, 4 (9) e202000996; DOI: 10.26508/lsa.202000996

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
X chromosome dosage compensation in Anopheles
Claudia Isabelle Keller Valsecchi, Eric Marois, M Felicia Basilicata, Plamen Georgiev, Asifa Akhtar
Life Science Alliance Jul 2021, 4 (9) e202000996; DOI: 10.26508/lsa.202000996
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
Issue Cover

In this Issue

Volume 4, No. 9
September 2021
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • ChIP
    • Data Availability
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Chromatin & Epigenetics
  • Evolution
  • Genetics, Gene Therapy & Genetic Disease

Related Articles

  • No related articles found.

Cited By...

  • No citing articles found.
  • Google Scholar

More in this TOC Section

  • NLRP3 controls ATM activation
  • RIF1 stabilizes replication domains
  • Chromatin context-dependent interferon response
Show more Research Article

Similar Articles

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2023 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.