Next Article in Journal
Impaired Replication Timing Promotes Tissue-Specific Expression of Common Fragile Sites
Next Article in Special Issue
Complex Structure of Lasiopodomys mandarinus vinogradovi Sex Chromosomes, Sex Determination, and Intraspecific Autosomal Polymorphism
Previous Article in Journal
A Novel Germline MLH1 In-Frame Deletion in a Slovenian Lynch Syndrome Family Associated with Uncommon Isolated PMS2 Loss in Tumor Tissue
Previous Article in Special Issue
Phylogenetic Analysis and Karyotype Evolution in Two Species of Core Gruiformes: Aramides cajaneus and Psophia viridis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 3
10.3390/genes11030327
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex

by
Atashi Sharma
1,
Nicholas A. Kinney
2,
Vladimir A. Timoshevskiy
1,
Maria V. Sharakhova
1,3,4 and
Igor V. Sharakhov
1,2,3,5,*
1
Department of Entomology, Virginia Polytechnic and State University, Blacksburg, VA 24061, USA
2
Genomics Bioinformatics and Computational Biology, Virginia Polytechnic and State University, Blacksburg, VA 24061, USA
3
Laboratory of Evolutionary Genomics of Insects, the Federal Research Center Institute of Cytology and Genetics, Siberian Branch of the Russian Academy of Sciences, 630090 Novosibirsk, Russia
4
Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 634050 Tomsk, Russia
5
Department of Cytology and Genetics, Tomsk State University, 634050 Tomsk, Russia
*
Author to whom correspondence should be addressed.
Genes 2020, 11(3), 327; https://doi.org/10.3390/genes11030327
Submission received: 19 February 2020 / Revised: 13 March 2020 / Accepted: 16 March 2020 / Published: 19 March 2020
(This article belongs to the Special Issue Mechanisms Driving Karyotype Evolution and Genomic Architecture)
Browse Figures
Versions Notes

Abstract

:
Heterochromatin is identified as a potential factor driving diversification of species. To understand the magnitude of heterochromatin variation within the Anopheles gambiae complex of malaria mosquitoes, we analyzed metaphase chromosomes in An. arabiensis, An. coluzzii, An. gambiae, An. merus, and An. quadriannulatus. Using fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA), a highly repetitive fraction of DNA, and heterochromatic Bacterial Artificial Chromosome (BAC) clones, we established the correspondence of pericentric heterochromatin between the metaphase and polytene X chromosomes of An. gambiae. We then developed chromosome idiograms and demonstrated that the X chromosomes exhibit qualitative differences in their pattern of heterochromatic bands and position of satellite DNA (satDNA) repeats among the sibling species with postzygotic isolation, An. arabiensis, An. merus, An. quadriannulatus, and An. coluzzii or An. gambiae. The identified differences in the size and structure of the X chromosome heterochromatin point to a possible role of repetitive DNA in speciation of mosquitoes. We found that An. coluzzii and An. gambiae, incipient species with prezygotic isolation, share variations in the relative positions of the satDNA repeats and the proximal heterochromatin band on the X chromosomes. This previously unknown genetic polymorphism in malaria mosquitoes may be caused by a differential amplification of DNA repeats or an inversion in the sex chromosome heterochromatin.

1. Introduction

The major malaria vector in Africa Anopheles gambiae has been a subject of extensive research over the past few decades. Initially described as a single species, An. gambiae was later subdivided into a complex of morphologically indistinguishable species by crossing experiments [1,2], fixed differences in polytene chromosome arrangement and distinct adaptations [3,4], differences in Intergenic Sequence (IGS) and Internal Transcribed Spacer (ITS) regions in the ribosomal DNA (rDNA) [5,6], and by whole-genome divergence [7,8,9]. Members of the An. gambiae complex differ in genetic diversity and they include major vectors, minor vector, and nonvectors [3,4,10]. The current list of the species includes An. arabiensis, An. amharicus, An. bwambae, An. coluzzii, An. fontenillei, An. gambiae, An. merus, An. melas, and An. quadriannulatus [4,7,9].
The majority of interspecies crosses in the An. gambiae complex produce fertile females and sterile males [2,11,12,13,14] in agreement with Haldane’s rule [15]. However, the most closely related species An. coluzzii and An. gambiae do not have postzygotic reproductive barriers as in laboratory conditions they readily mate and lay viable eggs that produce fertile hybrid males and females [16,17]. If the An. gambiae complex originated as recently as 526 thousand years ago, the An. coluzzii and An. gambiae lineages split from the common ancestor only ~61 thousand years ago [18]. The recently diverged species An. gambiae sensu stricto and An. coluzzii have been initially considered as the S and M molecular forms of An. gambiae [5,6] based on specific single nucleotide polymorphism (SNP) differences of ITS2 sequences [19,20]. Later studies demonstrated that the two forms also differ in genome sequence [8], gene expression [21], larval ecology [22], larval behavior in the presence of predators [23,24], adult swarming behavior [25], and adult mate recognition [17,26]. Thus, these data show evidence for ecological, behavioral, and genome-wide differentiation between An. coluzzii and An. gambiae.
Cytogenetic analysis of the polytene chromosomes banding patterns, including fixed chromosomal inversions, is an established tool for distinguishing species of malaria mosquitoes [3,4,27,28,29]. In addition, cytogenetic analysis of mitotic chromosomes demonstrated interspecific differences in sex chromosome heterochromatin between An. gambiae and An. arabiensis (named An. gambiae species A and B at that time) [30]. Staining with the Hoechst fluorescent dye has shown that the presence and brightness of X-chromosome heterochromatic bands differ between the two species [31]. More recently, heterochromatin variation of mitotic Y chromosomes among species of the An. gambiae complex has been clearly demonstrated [32]. Intraspecific polymorphism in the X and Y chromosome heterochromatin has also been observed in both natural and laboratory populations of An. gambiae [30,31], but is unclear if the previously observed variations can be related to possible differences between the more recently described incipient species An. coluzzii and An. gambiae s.s. [7]. This information can be useful because vector control to be successful must consider full spectrum of genetic and phenotypic variation within vector species. Further, it has been suggested that polymorphism in sex chromosome heterochromatin may affect fertility and sexual behavior of mosquitoes [33], but the lack of understanding of heterochromatin structure and function prevent researchers from mechanistic understanding of the phenotypic effects.
Sequencing of the An. gambiae genome provided important information regarding its organization [34]. It has been demonstrated that repetitive DNA component represent a substantial portion in of the An. gambiae genome (33%) that is higher than in Drosophila melanogaster genome where repeats make up approximately 24% of the genome [35,36]. The majority of repetitive DNA in the An. gambiae genome is tightly packed in heterochromatin around the centromeres [37]. Difficulty with sequencing of heterochromatin led to underrepresentation of heterochromatic sequences in the An. gambiae genome assembly [34]. Moreover, the repetitive nature of these sequences poses an impediment in mapping them correctly to chromosomes. Subsequent attempts to map heterochromatic genomic scaffolds led to an addition of ~16 Mb of heterochromatin to the genome of the An. gambiae PEST strain [38]. Further progress was made by predicting functions of 232 heterochromatin genes [39] and by mapping genes on the heterochromatin-euchromatin boundary of the polytene chromosome map [40]. Bioinformatics analysis of so called “unknown chromosome” or ~42 Mb of unmapped sequences in An. gambiae PEST genome demonstrated that it has characteristics of heterochromatin [39]. Although both An. gambiae and An. coluzzii genome assemblies are now available [39,41], corresponding information of their entire heterochromatin on a chromosome map is still missing. Likewise, sequencing of five other species from the An. gambiae complex also excluded the majority of heterochromatic sequences from the genome assembly [9,42].
Comparison of An. gambiae (the former S form) and An. coluzzii (the former M form) genomes identified pericentromeric autosomal and X-chromosome regions of high differentiation, termed “speciation islands,” or “islands of genomic divergence” [43,44]. These regions largely correspond to pericentromeric heterochromatin. Overlaps between the heterochromatin and the islands of genomic divergence are 91% in the X chromosome, 97% in the 2L arm, and 94% in the 3L arm [39]. Other studies emphasized that the highest genomic divergence between these nascent species occurred within the 4 Mb of mapped heterochromatin on the X chromosome [8,45,46]. A recent population genomic analysis of 765 field-collected mosquitoes across Africa has demonstrated that sequences in the pericentromeric X heterochromatin show the greatest separation between An. coluzzii and An. gambiae populations on a neighbor-joining tree [47]. The study also analyzed the An. coluzzii and An. gambiae genomes for CRISPR/Cas9 target sites; it found that 5474 genes had at least one viable target after excluding target sites with nucleotide variation. Interestingly, targetable genes are spread non-uniformly across the genome, falling predominantly in pericentromeric heterochromatic regions, where levels of variation are lower [47]. However, the full spectrum of heterochromatin variation and molecular composition of the X-chromosome heterochromatin in malaria mosquitoes remains unknown. Because repetitive DNA is underreplicated [48], polytene chromosomes have limited applications in studies of heterochromatin. In comparison, mitotic chromosomes are a promising system for evolutionary studies of heterochromatin in malaria mosquitoes.
In this study, we examined similarities and differences in heterochromatin patterns within X mitotic chromosomes among the major malaria vectors An. gambiae, An. coluzzii, An. arabiensis, minor vector An. merus, and zoophilic non-vector An. quadriannulatus. The correspondence of pericentric heterochromatin between the metaphase and polytene X chromosome was established for An. gambiae. We report quantitative differences in molecular organization of heterochromatin among members of the An. gambiae complex and highlight shared polymorphism of the heterochromatin structure among strains of the incipient species An. gambiae and An. coluzzii. The identified variation in heterochromatin is discussed with respect to the phylogenetic relationships among the species and a possible role of heterochromatin in speciation.

2. Materials and Methods

2.1. Mosquito Strains and Colony Maintenance

Laboratory colonies examined for this study were provided by the Biodefense and Emerging Infections Research Resources Repository and included colonies of An. gambiae PIMPERENA (MRA-861), An. gambiae KISUMU1 (MRA-762), An. gambiae ZANU (MRA-594), An. coluzzii SUA2La (MRA-765), An. coluzzii MOPTI (MRA-763), An. coluzzii MALI-NIH (MRA-860), An. arabiensis DONGOLA 2Ra2Rb3R (MRA-1235), An. quadriannulatus SANGWE (MRA-1155), and An. merus MAF (MRA-1156). The species and strains used in this study are presented in Table S1. All mosquitoes were reared at 27 °C, with 12:12 light:dark cycle and 70% relative humidity. Authentication of the species was done using PCR diagnostics [49,50]. Larvae were fed fish food and adult mosquitoes were fed 1% sugar water. To induce oviposition, females were fed defibrinated sheep blood (Colorado Serum Co., Denver, Colorado, USA) using artificial blood feeders. To perform interspecies crosses between An. gambiae and An. coluzzii, male and female pupae were separated to ensure virginity of adult mosquitoes. We differentiated males and females at the pupal stage using sex-specific differences in the shape of their terminalia [51]. After the emergence of adults, crossing experiments were performed by combining 30 females and 15 males in a single cage. Five days after random mating, the females were fed sheep blood. Two days later, an egg dish, covered with moist filter paper to keep the eggs from drying out, was put into the cage to obtain F1 hybrid larvae for cytogenetic analyses. Any possible colony contamination during the experiment was ruled out by verifying all the strains for their respective species using intentional mismatch primers (IMP) as described elsewhere [19]. All the strains showed the expected band size, with An. gambiae strains around 330 bp and An. coluzzii strains around 460 bp.

2.2. Chromosome Preparation

Preparations from early 4th instar larvae of lab colonies were made from leg and wing imaginal discs, as previously described [52]. Both male and female larvae were immobilized on ice for 10 min, and then dissected in a drop of cold, freshly prepared, hypotonic solution (0.075 M potassium chloride). After decapitating the head, the thorax was opened using dissecting scissors (Fine Science Tools, Foster City, CA, USA), followed by removal of the gut and fat from the body. A fresh drop of hypotonic solution was added to the preparation for 10 min, followed by fixation in a drop of modified Carnoy’s solution (ethanol:glacial acetic acid, 3:1) for 1 min. Next, a drop of freshly prepared 50% propionic acid was added, and the imaginal discs were covered with a 22 × 22 mm coverslip. After 5 min, the preparations were squashed using the flat rubber end of a pencil and dipped in liquid nitrogen until the bubbling stopped. Coverslips were removed using a sharp blade and slides were transferred to cold 50% ethanol and stored at −20 °C. After 2 h, slides were dehydrated in a series of 70%, 80%, and 100% ethanol. Preparations with the highest number of metaphase plates were chosen for in situ hybridization.

2.3. C0t DNA Preparation and DNA Probe Labeling

To identify the position of the centromere in chromosome X, we prepared the repetitive DNA fraction using a previously described method [52]. Genomic DNA was isolated from 500 g of non-bloodfed adult An. coluzzii MOPTI mosquitoes using a Qiagen Blood and Cell culture DNA Maxi kit (Qiagen, Hilden, Germany). Isolated C0t1 repetitive fraction was precipitated with isopropanol and labelled by nick-translation in 50 μl, containing 1 μg DNA, 0.05 mM each of unlabeled dATP, dCTP, and dGTP, and 0.015 mM of dTTP, 1 μl of fluorescein-dUTP, 0.05 mg/ml of BSA, 5 μl of 10× nick-translation buffer, 20 U of DNA-polymerase I, and 0.0012 U of DNase, at 15 °C for 2.5 h. Primers for 18S rDNA and satDNA repeats from An. gambiae AgY53A, AgY477-AgY53B, AgY477, and AgY53C were obtained, as previously described [14,53] (Table S2). An Immomix PCR kit (Bioline, Inc., Taunton, MA, USA) was used to label satellites by incorporating Cy3 and Cy5 fluorescently labeled nucleotides (Enzo Life Sciences, Inc., Farmingdale, NY, USA) directly into the PCR reaction. Each 25 µl PCR mix consisted of 35–40 ng genomic DNA, 0.3 U Taq polymerase, 1× PCR buffer, 200 µM each of dATP, dCTP, and dGTP, and 65 µM dTTP, and 0.5 µl Cy3-dUTP or 0.5 µl Cy5-dUTP (Enzo Life Sciences, Inc., Farmingdale, NY, USA). Thermocycling was performed using ImmoMix TM (Bioline Inc., Taunton, MA, USA) beginning with a 95 °C incubation for 10 min followed by 35 cycles of 95 °C for 30 sec, 52 °C for 30 sec, 72 °C for 45 sec, 72 °C for 5 min, and a final hold at 4 °C. Bacterial Artificial Chromosome (BAC) clones 05F01 (GenBank accession: AL142298, AL142299), 179F22 (GenBank accession: BH373300, BH373306), and 01K23 (GenBank accession: AL607293, AL607294) [34] were labelled by nick-translation in 50 μl, containing 1 μg DNA, 0.05 mM each of unlabeled dATP, dCTP, and dGTP and 0.015 mM of dTTP, 1 μl of Cy3-dUTP (or another fluorochrome), 0.05 mg/ml of BSA, 5 μl of 10× nick-translation buffer, 20 U of DNA-polymerase I, and 0.0012 U of DNase at 15 °C for 2.5 h.

2.4. Fluorescence in situ Hybridization (FISH)

Suitable slides with >10 metaphase plates were selected for FISH, which was performed as previously described [52,54]. Briefly, slides with good preparations were treated with 0.1 mg/ml RNase at 37 °C for 30 min. After washing twice with 2× saline-sodium citrate (SSC) for 5 min, slides were digested with 0.01% pepsin and 0.037% HCl solution for 5 min at 37 °C. After washing slides twice in 1× phosphate-buffered saline (PBS) for 5 min at room temperature, preparations were fixed in 3.7% formaldehyde for 10 min at RT. Slides were then washed in 1× PBS and dehydrated in a series of 70%, 80%, and 100% ethanol for 5 min at RT. Then, 10 μl of probes were mixed, added to the preparations, and incubated overnight at 37 °C. After washing slides in 1× SSC at 60 °C for 5 min, 4× SSC/NP40 solution at 37 °C for 10 min, and 1× PBS for 5 min at room temperature (RT), preparations were counterstained with a DAPI-antifade solution (Life Technologies, Carlsbad, CA, USA) and kept in the dark for at least 2 h before visualization with a fluorescence microscope.

2.5. Image Acquisition and Chromosome Measurements

Chromosome slide preparations were viewed with an Olympus BX61 fluorescence microscope (Olympus, Tokyo, Japan) using BioView software (BioView Inc., Billerica, MA, USA) at 1000× magnification. For idiogram development, metaphase plates were measured and analyzed from 5 larvae per strain. A total of 40 metaphase plates were used to measure chromosome lengths in An. coluzzii. For sibling species An. arabiensis, An. quadriannulatus, and An. merus, 11–13 metaphase plates were used to measure chromosome lengths (Table S3). Chromosome images were inverted using Adobe Photoshop CS6 (Adobe Inc., San Jose, CA, USA) and measured using the ruler tool at 1000× magnification. A statistical Tukey’s test or a nonparametric Kruskal–Wallis rank sum test followed by a Dunn’s test were performed to compare chromosome length between the species using JMP13 software (SAS Institute Inc., Cary, NC, USA) or R software (RStudio, Boston, MA, USA). Fluorescence intensities of proximal and distal X chromosome heterochromatin bands were measured using Adobe Photoshop CS6 (Adobe Inc., San Jose, CA, USA) and compared between three strains of An. gambiae and three strains of An. coluzzii. About forty measurements for each band were taken from forty sister chromatids of twenty X chromosomes for each strain. Pairwise comparisons between strains were performed using Student’s t-test.

2.6. Quantative Analysis of Fluorescent Signal Positions

Chromosomal positions of satDNA AgY477–AgY53B and AgY477, and the proximal DAPI band were identified and compared in three strains of An. gambiae and three strains of An. coluzzii. A custom MATLAB script [55] was written and used to measure the position of fluorescence peaks in individual channels [56]. The program was automated to measure the position of each satDNA and DAPI peak along the length of the X chromosome. The output provided the maximum likelihood of a combination of probe positions. The results were compared within and between species. The MATLAB program consisted of two user-guided steps followed by a third step of automated analysis. In the first step, sex chromosomes were identified by the user on mitotic slide preparations from An. gambiae and An. coluzzii strains. In the second step, the boundary of each identified sex chromosome was traced by a trained user. In an automated step, the pattern of heterochromatin and satDNA fluorescence was averaged longitudinally along the user-defined boundary. Program output displayed fluorescence intensity graphically from chromosome centromere to telomere. Peak fluorescence intensity was used to automatically infer the order of the heterochromatin and satDNA probes. User guidance was deliberately introduced to eliminate the pitfalls of full automation. In particular, user guidance was critical for accurate identification of sex chromosomes in the imaged mitotic slide preparations. Longitudinal averaging of fluorescence intensity during the program’s automated analysis was robust to the boundary of each chromosome traced by the user.

3. Results

3.1. Correspondence Between Heterochromatin of the Polytene and Mitotic X Chromosome in Anopheles gambiae

Using FISH with a C0t1 fraction of repetitive DNA, 18S rDNA, and heterochromatic BAC clones, we attempted to establish the correspondence of pericentric heterochromatin between the metaphase and polytene X chromosome in An. gambiae. The most repetitive C0t1 DNA fraction should presumably represent the centromeric satellites of the mitotic chromosomes corresponding to the most proximal regions of polytene arms. C0t analysis is the process of renaturation of single stranded DNA to its complementary sequence [57]. The analysis is based on the observation that the more repetitive DNA sequences require a shorter time to re-anneal following denaturation. DNA fractions with C0t1 values equal to 1.0 × 10−4–1.0 × 10−1 are considered as highly repetitive DNA fractions [57]. We expected C0t1 to map to narrow regions corresponding to centromeres in mitotic chromosomes. While that was the case for autosomes, in X chromosomes the C0t1 DNA probe localized to a wide region, including the entire proximal heterochromatin band adjacent to the rDNA locus (Figure 1A). Thus, hybridization of the C0t1 DNA fraction could be indicative of the centromere positions in autosomes, but we could not precisely determine the location of the centromere on the X chromosome in An. gambiae. We mapped three heterochromatic BAC clones, 05F01, 179F22, and 01K23, to the distal heterochromatin band of mitotic chromosomes in An. coluzzii (Figure 1B). According to the An. gambiae PEST genome assembly [38], these BAC clones are located in the most proximal heterochromatic region 6 in the X polytene chromosome map (Figure 1C). BAC clone 05F01 is mapped to scaffold AAAB01008973 (X: 20,698,730–20,818,594), BAC clone 01K23 is mapped to scaffold AAAB01008967 (X: 23,512,846–23569864), and BAC clone 179F22 is mapped to scaffold AAAB01008976 (X: 23,930,050–24,043,835) [38]. These results demonstrate that the assembled euchromatic part of the X chromosome occupies no more than half of the total X chromosome length (Figure 1C). Although polytene chromosomes are significantly larger than mitotic chromosomes, heterochromatin is under-replicated [48] and represented by diffuse structures without clear banding patterns [39], while heterochromatin in mitotic chromosomes is much more prevalent and detailed (Figure 1D). Because the An. gambiae PEST X chromosome genome assembly ends at the coordinate 24,393,108, our data indicate that the assembled heterochromatin largely corresponds to the distal heterochromatic band of the mitotic X chromosome. Only few copies of 18S rDNA sequences are found in the assembled X chromosome. Thus, most of the rDNA repeats and more proximal regions of the X chromosome heterochromatin are absent from the mapped genome assembly.

3.2. Variation in Heterochromatin Morphology Among Sibling Species of the An. gambiae Complex

We analyzed the heterochromatin morphology in species of the An. gambiae complex originating from different parts of Africa (Figure 2A) and diverging from each other from ~61,000 to ~526,000 years ago [18] (Figure 2B).
We studied mid-metaphase chromosomes obtained from imaginal discs of the 4th instar larvae of An. gambiae, An. coluzzii, An. arabiensis, An. merus, and An. quadriannulatus. Chromosomes at this stage provide reproducible heterochromatin patterns verified across multiple individuals. DAPI, a counter-stain that preferentially stains AT-rich DNA, was utilized to augment the natural banding patterns resulting from heterochromatin variation between and within species. To simplify the comparison, all fluorescence images were inverted, and converted to gray-scale images (Figure 3). The number and size of the X chromosome heterochromatin bands varied across the sibling species. Differences in the X chromosome heterochromatin pattern were clearly visible between all species with postzygotic reproductive isolation (Figure 3). Two relatively small heterochromatin bands can be seen in An. gambiae and An. coluzzii (Figure 3A,B). The KISUMU strain of An. gambiae often shows polymorphism in the size of the heterochromatin between different X chromosomes, even within one individual (Figure 1D and Figure 3A) Extended heterochromatin is observed in X chromosomes of An. arabiensis, An. quadriannulatus, and An. merus (Figure 3C). The pattern of fluorescence intensity of the heterochromatic bands differs among species as well. An. gambiae and An. coluzzii possess a dense heterochromatin block on the end proximal to the putative centromere. This is followed by a dull fluorescence area to which the 18S rDNA locus is mapped. A light distal band is present immediately next to the rDNA locus, marking the assembled heterochromatin in the current An. gambiae PEST genome assembly [38]. A large dark heterochromatin band is visible in the An. arabiensis X chromosome. Our results for An. gambiae and An. arabiensis corroborated with those for the natural populations of these species in a previous study [38]. In contrast, An. quadriannulatus, a zoophilic non-vector species of the An. gambiae complex, has one light and two dark heterochromatin bands. The longest X chromosome was found in An. merus, the marshy saltwater resident of the An. gambiae complex. The An. merus X chromosome contains three large heterochromatin blocks and a large rDNA locus making it comparable in length to the autosomes.

3.3. Molecular Variation of Heterochromatin Among Species of the An. gambiae Complex

We analyzed similarities and differences in the molecular content of the heterochromatin among sibling species. We performed FISH with 18S rDNA and previously identified satDNA sequences AgY53A, Ag53C, AgY477, and the genomic fragment AgY477–AgY53B that contains partial sequences of satellites AgY477 and AgY53B as well as a junction region [53] (Table S2). In An. gambiae PIMPERENA, satellites AgY53A and AgY477–AgY53B hybridized to single locations in the tip and base of the X chromosome proximal heterochromatin band, respectively (Figure 4A,B). This is in sharp contrast to An. merus, where AgY477–AgY53B hybridized with four distinct regions of the X chromosome heterochromatin (Figure 4C). Interestingly, AgY477–AgY53B hybridized differently in An. coluzzii SUA compared with An. gambiae PIMPERENA; AgY477–AgY53B is found in the base of the X chromosome proximal heterochromatin band of An. coluzzii SUA (Figure 4D). Satellite AgY477 has also been mapped to the base of the X chromosome proximal heterochromatin band of An. gambiae ZANU [32]. As described in the following results section, localization of these satellites represents shared variation among strains of An. gambiae and An. coluzzii. Unlike AgY53A and AgY477–AgY53B, satellite Ag53C is not sex chromosome-specific [53]; it hybridized with pericentromeric regions of autosomes in An. coluzzii MOPTI (Figure 4E). The 18S rDNA probe was mapped to the space between two dense heterochromatic bands on the X chromosome in An. quadriannulatus (Figure 4F). Previously, it was demonstrated that the AgY477–AgY53B sequence is immediately adjacent to the rDNA locus on the X chromosome of An. quadriannulatus [32]. Both AgY477–AgY53B and AgY477 co-localize at the very tip of the X chromosome outside the major dense heterochromatin band in An. arabiensis (Figure 4G). Hybridization of each of these satellites, together with 18S rDNA, demonstrated that, unlike in An. gambiae or An. coluzzii, AgY477 (Figure 4H) and AgY477–AgY53B (Figure 4I) are separated from the rDNA locus by the large band of dense heterochromatin in An. arabiensis.

3.4. Shared Heterochromatin Variation Between An. coluzzii and An. gambiae

To visualize the rDNA locus, we performed FISH of the 18S rDNA probe with mitotic chromosomes. Within all strains of An. gambiae and An. coluzzii, the rDNA locus was mapped between the proximal and distal heterochromatin bands (Figure S1). We noticed variation in the size of the rDNA locus both among the strains of An. gambiae and An. coluzzii and between homologous chromosomes within individual mosquitoes. The intra-strain polymorphism in the size of the rDNA locus was especially obvious in KISUMU, ZANU, and MOPTI. This polymorphism may suggest variation in the rDNA copy number.
To better resolve the relative satDNA positions among multiple mosquito strains, we generated F1 hybrids between laboratory strains of An. gambiae and An. coluzzii. F1 hybrid larvae resulting from the An. gambiae × An. coluzzii crosses were analyzed for heterochromatin morphology and satDNA location. Qualitative and quantitative differences in the pattern of heterochromatin blocks, rDNA loci, and satDNA location were clearly observed between homoeologous X chromosomes within the F1 hybrids. Different positions and sizes of the AgY477–AgY53B and AgY477 FISH signals with respect to the proximal heterochromatin band can be seen between homoeologous X chromosomes in F1 ♀An. gambiae PIMPERENA × ♂An. coluzzii MOPTI (Figure 5A). Different sizes of the proximal heterochromatin band and the rDNA locus can be seen between homoeologous X chromosomes in F1 ♀An. gambiae ZANU × ♂An. coluzzii MALI (Figure 5B). Also, different sizes of the hybridization signals from AgY477 and 18S rDNA can be seen between homologous X chromosomes in F1 ♀An. coluzzii MOPTI × ♂An. gambiae KISUMU (Figure 5C).
Using a custom MATLAB script [56], the position of the satDNA repeats with respect to the proximal heterochromatin band was mapped on the X chromosomes in An. gambiae and An. coluzzii strains. Because sequences of AgY477–AgY53B and AgY477 overlap, it is technically unfeasible to distinguish between the signals when mapping their distance from the proximal heterochromatin band. Therefore, we combined the fluorescence signals of AgY477–AgY53B and AgY477 and mapped them with respect to the fluorescence DAPI peak of the heterochromatin. The advantages of our automated analysis are twofold. First, fluorescent signals in each image possess a degree of overlap, which masks the order of fluorescent peaks from centromere to telomere; however, the order of peaks is easily detected using our automated analysis. Second, observation alone is insufficient to determine the distance between fluorescent peaks, which are generally separated by fewer than 10 pixels. Determining the distance between nearby peaks is made tractable in our pipeline; however, even this computational approach remains limited by microscopic resolution. Images of the mitotic sex chromosomes were captured using 1000× magnification; neighboring pixels at this resolution represent a distance of 0.05 μm. Thus, peaks spaced closer than 0.05 μm were prone to falsely inverted order. Nonetheless, peak orders were generally resolvable by aggregating multiple images averaged for multiple individuals of each strain. We measured the distance from satellites to the fluorescence DAPI peak and the linear order of the peaks. For example, the order of the satDNA and DAPI peaks is switched between An. gambiae PIMPERENA and An. coluzzii MOPTI (Figure 6A). Distribution patterns of signal intensities for satellites and DAPI along the X chromosomes also differ between An. gambiae PIMPERENA and An. coluzzii MOPTI (Figure 6B). Our mapping results revealed that the strains differed in the pattern of satDNA and DAPI peak relative positions and distance on the X chromosome (Figure 6C). Clustering analysis showed that An. gambiae PIMPERENA and An. coluzzii MALI depicted the same pattern of satDNA, followed by the DAPI band, while An. coluzzii MOPTI, An. coluzzii SUA, and An. gambiae ZANU clustered together showing the DAPI peak followed by the satDNA locus (Figure 6D). Interestingly, An. gambiae KISUMU had an almost 1:1 ratio of the direct and the reverse satDNA-DAPI band order indicating that this strain is maintaining high polymorphism in their X chromosome heterochromatin. Thus, our data shows the dynamism of the heterochromatin at the molecular level in the important malaria vectors. We found no evidence for clustering the strains by species, indicating that this heterochromatin variation is shared between An. coluzzii and An. gambiae either as ancestral polymorphism or by ongoing hybridization between the incipient species.
We also compared the average fluorescence intensities of the proximal and distal X chromosome heterochromatin between three strains of An. gambiae and three strains of An. coluzzii. A statistical analysis, using Student’s t-test, revealed no specific pattern when strains were compared with respect to the species. Most strains were significantly different from each other in their fluorescence intensities of the proximal band (P < 0.05), except for SUA-PIMPERENA (P = 0.06), SUA-ZANU (P = 0.16), ZANU-PIMPERENA (P = 0.61), and MALI-KISUMU (P = 0.74). MOPTI was significantly different from every other strain (Figure S2A). Comparison of fluorescence intensities of the distal heterochromatin band revealed a similar pattern for pairwise comparison between strains of both species. Most strains were significantly different from each other (P < 0.05), except for PIMPERENA-KISUMU (P = 0.07), ZANU-MALI (P = 0.11), and KISUMU-SUA (P = 0.89) (Figure S2B). Overall, the fluorescence intensities of both X chromosome heterochromatin bands vary among the strains, regardless of their species identity.

3.5. X Chromosome Idiograms for Species of the An. gambiae Complex

To construct idiograms, we measured the lengths of the chromosomes (Table S3) and conducted a statistical pairwise Tukey’s test between the sibling species. Our analysis revealed that the total X chromosome lengths were significantly different for all species pairs, except between An. quadriannulatus and An. arabiensis (P = 0.64). Due to the possibility that the X chromosome length data violated the assumption of equal variance, we repeated statistical analyses using a nonparametric Kruskal–Wallis rank sum test followed by a Dunn’s test. These analyses resulted in qualitatively identical results (Figure 7, Table S4). The polytene X chromosomes in all these species are similar in length; hence, the difference in mitotic X chromosomes can be attributed to the variation in heterochromatin.
We also conducted a statistical pairwise Tukey’s test for autosomes between the sibling species. We found that, unlike the X chromosome, the autosomal lengths do not significantly differ between the species (Figure 7, Table S5). This finding demonstrates the contrasting patterns of the sex chromosome and autosome evolution and highlights the rapids evolutionary changes in the X chromosome heterochromatin.
Based on mitotic chromosome measurements, banding patterns, and repeat mapping, we constructed X chromosome idiograms for An. gambiae, An. coluzzii, An. arabiensis, An. quadriannulatus, and An. merus. FISH with the genomic fragment AgY477–AgY53B revealed its multiple locations in the X chromosome heterochromatin bands in An. merus but a single location in the pericentric heterochromatin in other sibling species. These features were placed on respective idiograms for easy comparison of molecular features within heterochromatin among species (Figure 8). The identified differences in the X chromosome size and heterochromatin structure are consistent with An. gambiae, An. coluzzii, An. arabiensis, An. quadriannulatus grouping together while An. merus being a separate lineage in the species phylogeny [18] (Figure 2B). The An. merus X chromosome has the largest amount of heterochromatin among the X chromosomes of sibling species. Moreover, An. merus is the only species in this study with the rDNA locus located in the L arm of the X chromosome suggesting a fixed pericentric inversion that differentiate An. merus from the other species. The developed idiograms will facilitate future studies of genetic variation in repeat-rich heterochromatic regions of malaria mosquitoes from the An. gambiae complex.

4. Discussion

Mosquitoes, with their small number of diploid chromosomes, 2n = 6, present a convenient model to investigate the molecular organization of their genomes and chromosome evolution [58]. Since polytene chromosomes are well-developed in Anopheles mosquitoes, they are often preferred over mitotic chromosomes for gene mapping [40,59]. Analysis of polytene chromosomes has been useful for identifying genomic variations in euchromatin. Comparative genome mapping within the Anopheles genus has demonstrated that the euchromatic portion of the X chromosome is the fastest evolving by genomic rearrangements [42,60]. However, under-replication of heterochromatin in polytene chromosomes makes them less suitable for mapping repetitive elements of the genome. In this study, we utilized metaphase mitotic chromosomes for repeat mapping and for studying heterochromatin variation among species in the An. gambiae complex. Since the genome of An. gambiae was first published [34], polytene chromosomes have been used to improve the heterochromatin assembly and to characterize its genetic content [38,39,40]. Here, we established the correspondence of pericentric heterochromatin between the metaphase and polytene X chromosome of An. gambiae, and demonstrated that at least half of the An. gambiae X chromosome length is heterochromatic. Future long-read sequencing and assembly of the heterochromatin in malaria mosquitoes may yield important insights into structural and functional organization of the sex chromosomes, as demonstrated for Drosophila [61].
Our comparative analysis of the X chromosome pericentric heterochromatin identified qualitative differences among sibling species, An. arabiensis, An. merus, An. quadriannulatus, and An. coluzzii or An. gambiae, with postzygotic isolation. The identified fixed variations included the number and pattern of heterochromatic bands and chromosomal location of the satDNA repeats. A common emerging theme from studies performed in Drosophila suggests that heterochromatin plays an important role in hybrid incompatibility [62,63,64,65,66]. Rapidly evolving heterochromatic repeats and heterochromatin-binding proteins are considered major players in the evolution of intrinsic reproductive isolation. Reproductively isolated Drosophila species have different satDNA content and localization within heterochromatin [64,67]. Evolutionary forces lead to a great variation in the copy number and types of heterochromatin repeats between closely related species. As repetitive DNA elements replicate, mechanisms such as unequal crossover, rolling circle replication, and segmental duplication become sources of genome evolution [68,69]. Thus, our study supports observations from other organisms that sex chromosomes have a propensity to accumulate species-specific differences in heterochromatin [70].
SatDNA repeats are seen as prime candidates to trigger genome instability in interspecies hybrids [71]. They have been implicated in disruption of chromosome pairing, abnormal heterochromatin packaging, and selfish post-meiotic drive substitutions in Drosophila hybrids [67]. Disruption of chromosome pairing and synapsis is observed frequently in interspecies hybrids of various organisms. Interestingly, sex chromosomes more often than autosomes drive speciation and hybrid incompatibilities [70]. For example, pairing of X and Y chromosomes is more adversely affected than that of autosomes during the prophase I in male hybrids between Campbell’s dwarf hamster and the Djungarian hamster [72]. A recent work has clarified that the autosomes of male hybrids between these hamster species undergo pairing and recombination normally as their parental forms do, but the heterochromatic arms of the X and Y chromosomes show a high frequency of asynapsis and recombination failure [73]. It was proposed that asynapsis of heterospecific chromosomes in prophase I may provide a recurrently evolving trigger for the meiotic arrest of interspecific F1 hybrids of mice [74,75]. Our recent study showed that meiotic failures in sterile F1 ♀A. merus × ♂A. coluzzii hybrids are accompanied by the disruption of sex chromosome pairing and insufficient chromosome condensation. These cytogenetic abnormalities lead to a premature equational division and formation of diploid immotile sperm [14]. Demonstrated here, qualitative inter-species differences in heterochromatin suggest that they may be responsible for incompatibilities of sibling species from the An. gambiae complex.
Unlike species with postzygotic isolation, An. coluzzii and An. gambiae do not show fixed differences in the pattern of the X chromosome pericentric heterochromatin. An early study observed an intra-specific polymorphism of the X chromosome heterochromatin among lab strains of An. gambiae [31]. The described polymorphism was based on the Hoechst staining of the heterochromatin. Here, we studied the pattern of heterochromatic bands or location of satDNA in multiple strains of An. coluzzii and An. gambiae, incipient species with prezygotic isolation. We found that, although the strains may differ from each other, the species do share variation in the relative positions of the fluorescence intensity peaks for the satDNA repeats and the proximal DAPI-positive heterochromatin band of the X chromosomes. The mechanism behind this inversion-like genetic polymorphism in malaria mosquitoes could be a rapid differential amplification of satellite repeats or actual structural rearrangement in the sex chromosome heterochromatin. It is possible that further evolution of heterochromatin in An. coluzzii and An. gambiae will contribute to higher genetic differentiation, and will eventually lead to postzygotic reproductive isolation between some populations of these important disease vectors. Extensive studies of Drosophila species show that heterochromatin, in general, and satellite DNA, in particular, can be a source of tremendous genetic variation [76,77,78,79,80]. Heterochromatin can also modulate differential gene expression and cause variable phenotypes, including differences in immune response [81,82]. It is possible that changes in repetitive sequences can have fitness consequences upon which selection will act. A population genomics study in D. melanogaster found that satDNA repeats often show population differentiation, but the population structure inferred from overall satellite quantities does not recapitulate the expected population relationships based on the demographic history of this species [79]. Moreover, some satellites have likely been involved in antagonistic interactions, as inferred from negative correlations among them.
The rDNA locus is immediately adjacent to the heterochromatin bands of the X chromosome. Its variation in size, among and within species, suggests polymorphism of the rDNA copy number. Natural polymorphism in the number of rDNA copies has been recorded in vertebrate and invertebrate animals, both among and within species [83]. Associations between the heterochromatin or rDNA variation and organismal phenotypes have also been demonstrated. A study in Drosophila showed that heterochromatin formation prolongs lifespan and controls ribosomal RNA synthesis [84]. Another study demonstrated that rDNA copy number decreases during aging, and that this age-dependent decrease in rDNA copy number is transgenerationally heritable [85]. In humans, the number of rDNA repeats can vary from 250 to 670 copies per diploid genome [86], and the changes in rDNA copy number can affect genome-wide gene expression [87]. Natural variation in rDNA and heterochromatin may contribute to genome evolution, formation of reproductive barriers, and eventually to diversification of species [88,89,90]. Genomics studies of the heterochromatin in natural populations of An. coluzzii and An. gambiae can determine the actual levels of polymorphism and divergence in satDNA and rDNA copy number variation. These studies may also identify heterochromatin variations that correlate with specific environmental adaptations or behaviors of the malaria vectors.

5. Conclusions

We discovered a new type of shared cytogenetic polymorphism in the incipient species, An. gambiae and An. coluzzii—an inversion of the satDNA location with respect to the proximal heterochromatin band. This finding suggests mechanisms of rapid differentiation of the sex chromosome heterochromatin during evolution: genomic rearrangement or differential amplification of heterochromatic repeats. Future genome sequencing using long-read technologies will be crucial for determining the precise nature of the observed polymorphism. Our study also demonstrated X chromosome heterochromatin divergence among mosquito species with post-zygotic isolation, which manifests in species-specific localization of repetitive DNA as well as size and number of heterochromatin bands. The identified differences in heterochromatin support the basal phylogenetic position of An. merus and point to the role of heterochromatin in speciation. The differences in molecular organization of the sex chromosome heterochromatin may impair meiotic synapsis during meiosis in sterile inter-species male hybrids. We suggest that inter-species divergence of heterochromatin structure represents a cytogenetic threshold that triggers the evolutionary shift from prezygotic to postzygotic isolation. The new idiograms developed here for the sibling species will aid in future studies of heterochromatin organization and evolution in these malaria mosquitoes.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/3/327/s1, Figure S1: FISH of 18S rDNA with the X chromosomes in strains of An. gambiae and An. coluzzii. (A) An. gambiae strains: PIMPERENA, ZANU, KISUMU. (B) An. coluzzii strains: MALI, SUA, MOPTI. Scale bar = 1 µM. Red = 18S rDNA, Blue = DAPI., Figure S2: Fluorescence intensities of X chromosome heterochromatin in three strains of An. gambiae and three strains of An. coluzzii. (A) Fluorescence intensities of the proximal heterochromatin band. (B) Fluorescence intensities of the distal heterochromatin band., Table S1: Mosquito strains used in the study of the X chromosome heterochromatin., Table S2: Probes and primer sequences used for FISH., Table S3: Length of mitotic chromosomes in species of the An. gambiae complex. Table S4: Statistical analyses of the X chromosome lengths using a nonparametric Kruskal-Wallis rank sum test followed by a Dunn’s test., Table S5: Statistical analyses of the chromosome lengths using Tukey’s honestly significant difference test.

Author Contributions

Conceptualization, M.V.S., and I.V.S.; methodology, A.S., N.A.K., V.A.T., and I.V.S.; software, N.A.K.; validation, A.S., M.V.S., and I.V.S.; formal analysis, A.S., N.A.K. and I.V.S.; investigation, A.S., N.A.K., and V.A.T.; resources, I.V.S.; data curation, A.S., N.A.K, M.V.S., and I.V.S.; writing—original draft preparation, A.S., M.V.S., and I.V.S.; writing—review and editing, I.V.S.; visualization, A.S., N.A.K., V.A.T.; supervision, M.V.S., and I.V.S.; project administration, I.V.S.; funding acquisition, M.V.S. and I.V.S. All authors have read and agreed to the published version of the manuscript.

Funding

The chromosome mapping and analyses were supported by the US NIH NIAID grants R21AI099528 and R21AI135298 to I.V.S., the graduate program in Genomics Bioinformatics and Computational Biology (GBCB) of Virginia Tech to N.A.K., the Fralin Life Sciences Institute, and the USDA National Institute of Food and Agriculture Hatch project 223822 to I.V.S. The development of chromosome idiograms was supported by the Russian Science Foundation grant No. 19-14-00130 to M.V.S.

Acknowledgments

The following reagents were obtained through the Biodefense and Emerging Infections Research Resources Repository, NIAID, NIH: An. arabiensis, Strain DONGOLA 2Ra, 2Rb/b and 3R, MRA-1235, contributed by Ellen M. Dotson; An. coluzzii, Strain MALI-NIH, Eggs, MRA-860, contributed by Nora J. Besansky; An. coluzzii, Strain MOPTI, Eggs, MRA-763, contributed by Gregory C. Lanzaro; An. coluzzii, Strain SUA2La, Eggs, MRA-765, contributed by Alessandra della Torre; An. gambiae, Strain KISUMU1, Eggs, MRA-762, contributed by Vincent Corbel; An. gambiae, Strain PIMPERENA (S form), Eggs, MRA-861, contributed by Nora J. Besansky; An. gambiae, Strain ZANU, MRA-594, contributed by Hilary Ranson and Frank H. Collins; An. merus, Strain MAF, MRA-1156, contributed by Maureen Coetzee; and An. quadriannulatus, Strain SANGWE, Eggs, MRA-1155, contributed by Willem Takken. We thank Semen Bondarenko for his help with statistical analysis using R, Jiangtao Liang for useful discussions, as well as Jean Clarke, Kristin Rose, and Janet Webster for proofreading the text.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Davidson, G. The Five Mating Types in the Anopheles Gambiae Complex. Riv. Malariol. 1964, 43, 167–183. [Google Scholar] [PubMed]
  2. Davidson, G.; Paterson, H.E.; Coluzzi, M.; Mason, G.F.; Micks, D.W. The Anopheles gambiae Complex. In Genetics of Insect Vectors of Diesease; Wright, J.W., Pal, R., Eds.; Elsevier Publishing Company: Amsterdam, The Netherlands; London, UK; New York, NY, USA, 1967; pp. 211–249. [Google Scholar]
  3. Coluzzi, M.; Sabatini, A.; Petrarca, V.; Di Deco, M.A. Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans. R. Soc. Trop. Med. Hyg. 1979, 73, 483–497. [Google Scholar] [CrossRef]
  4. Coluzzi, M.; Sabatini, A.; della Torre, A.; Di Deco, M.A.; Petrarca, V. A polytene chromosome analysis of the Anopheles gambiae species complex. Science 2002, 298, 1415–1418. [Google Scholar] [CrossRef] [PubMed]
  5. Della Torre, A.; Costantini, C.; Besansky, N.J.; Caccone, A.; Petrarca, V.; Powell, J.R.; Coluzzi, M. Speciation within Anopheles gambiae—The glass is half full. Science 2002, 298, 115–117. [Google Scholar] [CrossRef]
  6. Della Torre, A.; Fanello, C.; Akogbeto, M.; Dossou-yovo, J.; Favia, G.; Petrarca, V.; Coluzzi, M. Molecular evidence of incipient speciation within Anopheles gambiae s.s. in West Africa. Insect Mol. Biol. 2001, 10, 9–18. [Google Scholar] [CrossRef]
  7. Coetzee, M.; Hunt, R.H.; Wilkerson, R.; Della Torre, A.; Coulibaly, M.B.; Besansky, N.J. Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa 2013, 3619, 246–274. [Google Scholar] [CrossRef] [Green Version]
  8. Lawniczak, M.K.; Emrich, S.J.; Holloway, A.K.; Regier, A.P.; Olson, M.; White, B.; Redmond, S.; Fulton, L.; Appelbaum, E.; Godfrey, J.; et al. Widespread divergence between incipient Anopheles gambiae species revealed by whole genome sequences. Science 2010, 330, 512–514. [Google Scholar] [CrossRef] [Green Version]
  9. Barron, M.G.; Paupy, C.; Rahola, N.; Akone-Ella, O.; Ngangue, M.F.; Wilson-Bahun, T.A.; Pombi, M.; Kengne, P.; Costantini, C.; Simard, F.; et al. A new species in the major malaria vector complex sheds light on reticulated species evolution. Sci. Rep. 2019, 9, 14753. [Google Scholar] [CrossRef] [Green Version]
  10. White, B.J.; Collins, F.H.; Besansky, N.J. Evolution of Anopheles gambiae in Relation to Humans and Malaria. Annu. Rev. Ecol. Evol. Syst. 2011, 42, 111–132. [Google Scholar] [CrossRef]
  11. Slotman, M.; Della Torre, A.; Powell, J.R. The genetics of inviability and male sterility in hybrids between Anopheles gambiae and An. arabiensis. Genetics 2004, 167, 275–287. [Google Scholar] [CrossRef] [Green Version]
  12. Bernardini, F.; Galizi, R.; Wunderlich, M.; Taxiarchi, C.; Kranjc, N.; Kyrou, K.; Hammond, A.; Nolan, T.; Lawniczak, M.N.K.; Papathanos, P.A.; et al. Cross-Species Y Chromosome Function Between Malaria Vectors of the Anopheles gambiae Species Complex. Genetics 2017, 207, 729–740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Presgraves, D.C.; Orr, H.A. Haldane’s rule in taxa lacking a hemizygous X. Science 1998, 282, 952–954. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Liang, J.; Sharakhov, I.V. Premeiotic and meiotic failures lead to hybrid male sterility in the Anopheles gambiae complex. Proc. Biol. Sci. 2019, 286, 20191080. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Haldane, J.B.S. Sex ratio and unisexual sterility in hybrid animals. J. Genet. 1922, 12, 101–109. [Google Scholar] [CrossRef] [Green Version]
  16. Diabate, A.; Dabire, R.K.; Millogo, N.; Lehmann, T. Evaluating the effect of postmating isolation between molecular forms of Anopheles gambiae (Diptera: Culicidae). J. Med. Entomol. 2007, 44, 60–64. [Google Scholar] [CrossRef]
  17. Aboagye-Antwi, F.; Alhafez, N.; Weedall, G.D.; Brothwood, J.; Kandola, S.; Paton, D.; Fofana, A.; Olohan, L.; Betancourth, M.P.; Ekechukwu, N.E.; et al. Experimental swap of Anopheles gambiae’s assortative mating preferences demonstrates key role of X-chromosome divergence island in incipient sympatric speciation. PLoS Genet. 2015, 11, e1005141. [Google Scholar] [CrossRef]
  18. Thawornwattana, Y.; Dalquen, D.; Yang, Z. Coalescent Analysis of Phylogenomic Data Confidently Resolves the Species Relationships in the Anopheles gambiae Species Complex. Mol. Biol. Evol. 2018, 35, 2512–2527. [Google Scholar] [CrossRef]
  19. Wilkins, E.E.; Howell, P.I.; Benedict, M.Q. IMP PCR primers detect single nucleotide polymorphisms for Anopheles gambiae species identification, Mopti and Savanna rDNA types, and resistance to dieldrin in Anopheles arabiensis. Malar. J. 2006, 5, 125. [Google Scholar] [CrossRef] [Green Version]
  20. Santolamazza, F.; Mancini, E.; Simard, F.; Qi, Y.; Tu, Z.; della Torre, A. Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms. Malar. J. 2008, 7, 163. [Google Scholar] [CrossRef] [Green Version]
  21. Caputo, B.; Nwakanma, D.; Caputo, F.P.; Jawara, M.; Oriero, E.C.; Hamid-Adiamoh, M.; Dia, I.; Konate, L.; Petrarca, V.; Pinto, J.; et al. Prominent intraspecific genetic divergence within Anopheles gambiae sibling species triggered by habitat discontinuities across a riverine landscape. Mol. Ecol. 2014, 23, 4574–4589. [Google Scholar] [CrossRef]
  22. Gimonneau, G.; Pombi, M.; Choisy, M.; Morand, S.; Dabire, R.K.; Simard, F. Larval habitat segregation between the molecular forms of the mosquito Anopheles gambiae in a rice field area of Burkina Faso, West Africa. Med. Vet. Entomol. 2012, 26, 9–17. [Google Scholar] [CrossRef] [PubMed]
  23. Dabire, K.R.; Diabate, A.; Agostinho, F.; Alves, F.; Manga, L.; Faye, O.; Baldet, T. Distribution of the members of Anopheles gambiae and pyrethroid knock-down resistance gene (kdr) in Guinea-Bissau, West Africa. Bull. Soc. Pathol. Exot. 2008, 101, 119–123. [Google Scholar] [PubMed]
  24. Gimonneau, G.; Pombi, M.; Dabire, R.K.; Diabate, A.; Morand, S.; Simard, F. Behavioural responses of Anopheles gambiae sensu stricto M and S molecular form larvae to an aquatic predator in Burkina Faso. Parasit. Vectors 2012, 5, 65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Tripet, F.; Toure, Y.T.; Taylor, C.E.; Norris, D.E.; Dolo, G.; Lanzaro, G.C. DNA analysis of transferred sperm reveals significant levels of gene flow between molecular forms of Anopheles gambiae. Mol. Ecol. 2001, 10, 1725–1732. [Google Scholar] [CrossRef]
  26. Cassone, B.J.; Kamdem, C.; Cheng, C.; Tan, J.C.; Hahn, M.W.; Costantini, C.; Besansky, N.J. Gene expression divergence between malaria vector sibling species Anopheles gambiae and An. coluzzii from rural and urban Yaounde Cameroon. Mol. Ecol. 2014, 23, 2242–2259. [Google Scholar] [CrossRef] [Green Version]
  27. Sharakhova, M.V.; Peery, A.; Antonio-Nkondjio, C.; Xia, A.; Ndo, C.; Awono-Ambene, P.; Simard, F.; Sharakhov, I.V. Cytogenetic analysis of Anopheles ovengensis revealed high structural divergence of chromosomes in the Anopheles nili group. Infect. Genet. Evol. 2013, 16, 341–348. [Google Scholar] [CrossRef] [Green Version]
  28. Stegnii, V.N.; Novikov, Y.M.; Kabanova, V.M. Cytogenetic analysis and distribution of the malaria mosquito Anopheles beklemishevi. Zool. ZH 1978, 57, 873–876. [Google Scholar]
  29. Stegnii, V.N.; Kabanova, V.M. Chromosome analysis of the malarial mosquitoes Anopheles atroparvus and A. maculipennis (Diptera, Culicidae). Zool. ZH 1978, 8, 207–225. [Google Scholar]
  30. Gatti, M.; Santini, G.; Pimpinelli, S.; Coluzzi, M. Fluorescence banding techniques in the identification of sibling species of the Anopheles gambiae complex. Heredity (Edinb) 1977, 38, 105–108. [Google Scholar] [CrossRef] [Green Version]
  31. Bonaccorsi, S.; Santini, G.; Gatti, M.; Pimpinelli, S.; Colluzzi, M. Intraspecific polymorphism of sex chromosome heterochromatin in two species of the Anopheles gambiae complex. Chromosoma 1980, 76, 57–64. [Google Scholar] [CrossRef]
  32. Hall, A.B.; Papathanos, P.A.; Sharma, A.; Cheng, C.; Akbari, O.S.; Assour, L.; Bergman, N.H.; Cagnetti, A.; Crisanti, A.; Dottorini, T.; et al. Radical remodeling of the Y chromosome in a recent radiation of malaria mosquitoes. Proc. Natl. Acad. Sci. USA 2016, 113, E2114–E2123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Fraccaro, M.; Tiepolo, L.; Laudani, U.; Marchi, A.; Jayakar, S.D. Y chromosome controls mating behaviour on Anopheles mosquitoes. Nature 1977, 265, 326–328. [Google Scholar] [CrossRef] [PubMed]
  34. Holt, R.A.; Subramanian, G.M.; Halpern, A.; Sutton, G.G.; Charlab, R.; Nusskern, D.R.; Wincker, P.; Clark, A.G.; Ribeiro, J.M.; Wides, R.; et al. The genome sequence of the malaria mosquito Anopheles gambiae. Science 2002, 298, 129–149. [Google Scholar] [CrossRef] [PubMed]
  35. Smith, C.D.; Shu, S.; Mungall, C.J.; Karpen, G.H. The Release 5.1 annotation of Drosophila melanogaster heterochromatin. Science 2007, 316, 1586–1591. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Manning, J.E.; Schmid, C.W.; Davidson, N. Interspersion of repetitive and nonrepetitive DNA sequences in the Drosophila melanogaster genome. Cell 1975, 4, 141–155. [Google Scholar] [CrossRef]
  37. Besansky, N.J.; Collins, F.H. The mosquito genome: Organization, evolution and manipulation. Parasitol. Today 1992, 8, 186–192. [Google Scholar] [CrossRef]
  38. Sharakhova, M.V.; Hammond, M.P.; Lobo, N.F.; Krzywinski, J.; Unger, M.F.; Hillenmeyer, M.E.; Bruggner, R.V.; Birney, E.; Collins, F.H. Update of the Anopheles gambiae PEST genome assembly. Genome Biol. 2007, 8, R5. [Google Scholar] [CrossRef] [Green Version]
  39. Sharakhova, M.V.; George, P.; Brusentsova, I.V.; Leman, S.C.; Bailey, J.A.; Smith, C.D.; Sharakhov, I.V. Genome mapping and characterization of the Anopheles gambiae heterochromatin. BMC Genom. 2010, 11, 459. [Google Scholar] [CrossRef] [Green Version]
  40. George, P.; Sharakhova, M.V.; Sharakhov, I.V. High-resolution cytogenetic map for the African malaria vector Anopheles gambiae. Insect Mol. Biol. 2010, 19, 675–682. [Google Scholar] [CrossRef] [Green Version]
  41. Kingan, S.B.; Heaton, H.; Cudini, J.; Lambert, C.C.; Baybayan, P.; Galvin, B.D.; Durbin, R.; Korlach, J.; Lawniczak, M.K.N. A High-Quality De novo Genome Assembly from a Single Mosquito Using PacBio Sequencing. Genes 2019, 10. [Google Scholar] [CrossRef] [Green Version]
  42. Neafsey, D.E.; Waterhouse, R.M.; Abai, M.R.; Aganezov, S.S.; Alekseyev, M.A.; Allen, J.E.; Amon, J.; Arca, B.; Arensburger, P.; Artemov, G.; et al. Mosquito genomics. Highly evolvable malaria vectors: The genomes of 16 Anopheles mosquitoes. Science 2015, 347, 1258522. [Google Scholar] [CrossRef] [Green Version]
  43. Turner, T.L.; Hahn, M.W.; Nuzhdin, S.V. Genomic islands of speciation in Anopheles gambiae. PLoS Biol. 2005, 3, e285. [Google Scholar] [CrossRef]
  44. White, B.J.; Cheng, C.; Simard, F.; Costantini, C.; Besansky, N.J. Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae. Mol. Ecol. 2010, 19, 925–939. [Google Scholar] [CrossRef] [Green Version]
  45. Nwakanma, D.C.; Neafsey, D.E.; Jawara, M.; Adiamoh, M.; Lund, E.; Rodrigues, A.; Loua, K.M.; Konate, L.; Sy, N.; Dia, I.; et al. Breakdown in the process of incipient speciation in Anopheles gambiae. Genetics 2013, 193, 1221–1231. [Google Scholar] [CrossRef] [Green Version]
  46. Neafsey, D.E.; Lawniczak, M.K.; Park, D.J.; Redmond, S.N.; Coulibaly, M.B.; Traore, S.F.; Sagnon, N.; Costantini, C.; Johnson, C.; Wiegand, R.C.; et al. SNP genotyping defines complex gene-flow boundaries among African malaria vector mosquitoes. Science 2010, 330, 514–517. [Google Scholar] [CrossRef] [Green Version]
  47. Miles, A.; Harding, N.J.; Botta, G.; Clarkson, C.; Antao, T.; Kozak, K.; Schrider, D.; Kern, A.; Redmond, S.; Sharakhov, I.V.; et al. Genetic diversity of the African malaria vector Anopheles gambiae. Nature 2017, 552, 96–100. [Google Scholar]
  48. Zhimulev, I.F. Polytene chromosomes, heterochromatin, and position effect variegation. Adv. Genet. 1998, 37, 1–555. [Google Scholar]
  49. Scott, J.A.; Brogdon, W.G.; Collins, F.H. Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am. J. Trop. Med. Hyg. 1993, 49, 520–529. [Google Scholar] [CrossRef] [Green Version]
  50. Fanello, C.; Santolamazza, F.; della Torre, A. Simultaneous identification of species and molecular forms of the Anopheles gambiae complex by PCR-RFLP. Med. Vet. Entomol. 2002, 16, 461–464. [Google Scholar] [CrossRef]
  51. Benedict, M.; Dotson, E.M. Methods in Anopheles Research. 2015, p. 408. Available online: https://www.beiresources.org/Portals/2/VectorResources/2016%20Methods%20in%20Anopheles%20Research%20full%20manual.pdf (accessed on 12 December 2019).
  52. Timoshevskiy, V.A.; Sharma, A.; Sharakhov, I.V.; Sharakhova, M.V. Fluorescent in situ Hybridization on Mitotic Chromosomes of Mosquitoes. J. Vis. Exp. 2012, e4215. [Google Scholar] [CrossRef] [Green Version]
  53. Krzywinski, J.; Sangare, D.; Besansky, N.J. Satellite DNA from the Y chromosome of the malaria vector Anopheles gambiae. Genetics 2005, 169, 185–196. [Google Scholar] [CrossRef] [Green Version]
  54. Sharakhova, M.V.; George, P.; Timoshevskiy, V.; Sharma, A.; Peery, A.; Sharakhov, I.V. Mosquitoes (Diptera). In Protocols for Cytogenetic Mapping of Arthropod Genomes; Sharakhov, I.V., Ed.; CRC Press; Taylor & Francis Group: Boca Raton, FL, USA, 2015; pp. 93–170. [Google Scholar]
  55. MATLAB. MATLAB:2010 Version 7.10.0 (R2010a); The MathWorks Inc.: Natick, MA, USA, 2010. [Google Scholar]
  56. Kinney, N.; Sharakhov, I.V. MATLAB Script to Measure the Position of Fluorescence Peaks on Chromosomes. Figshare. Software. Available online: https://doi.org/10.6084/m9.figshare.11782140.v1 (accessed on 31 January 2020).
  57. Trifonov, V.A.; Vorobieva, N.V.; Serdyukova, N.A.; Rens, W. FISH with and Without COT1 DNA. In Fluorescence In Situ Hybridization (FISH)—Application Guide; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar] [CrossRef]
  58. Kumar, A.; Rai, K.S. Chromosomal localization and copy number of 18S + 28S ribosomal RNA genes in evolutionarily diverse mosquitoes (Diptera, Culicidae). Hereditas 1990, 113, 277–289. [Google Scholar] [CrossRef]
  59. Grushko, O.G.; Rusakova, A.M.; Sharakhova, M.V.; Sharakhov, I.V.; Stegnii, V.N. Localization of repetitive DNA sequences in the pericentromeric heterochromatin of malarial mosquitoes of the “Anopheles maculipennis” complex. Tsitologiia 2006, 48, 240–245. [Google Scholar]
  60. Jiang, X.; Peery, A.; Hall, A.; Sharma, A.; Chen, X.G.; Waterhouse, R.M.; Komissarov, A.; Riehl, M.M.; Shouche, Y.; Sharakhova, M.V.; et al. Genome analysis of a major urban malaria vector mosquito, Anopheles stephensi. Genome Biol. 2014, 15, 459. [Google Scholar] [CrossRef]
  61. Chang, C.H.; Larracuente, A.M. Heterochromatin-Enriched Assemblies Reveal the Sequence and Organization of the Drosophila melanogaster Y Chromosome. Genetics 2019, 211, 333–348. [Google Scholar] [CrossRef] [Green Version]
  62. Brideau, N.J.; Flores, H.A.; Wang, J.; Maheshwari, S.; Wang, X.; Barbash, D.A. Two Dobzhansky-Muller genes interact to cause hybrid lethality in Drosophila. Science 2006, 314, 1292–1295. [Google Scholar] [CrossRef] [Green Version]
  63. Cattani, M.V.; Presgraves, D.C. Genetics and lineage-specific evolution of a lethal hybrid incompatibility between Drosophila mauritiana and its sibling species. Genetics 2009, 181, 1545–1555. [Google Scholar] [CrossRef] [Green Version]
  64. Satyaki, P.R.; Cuykendall, T.N.; Wei, K.H.; Brideau, N.J.; Kwak, H.; Aruna, S.; Ferree, P.M.; Ji, S.; Barbash, D.A. The Hmr and Lhr hybrid incompatibility genes suppress a broad range of heterochromatic repeats. PLoS Genet. 2014, 10, e1004240. [Google Scholar] [CrossRef] [Green Version]
  65. Cattani, M.V.; Presgraves, D.C. Incompatibility between X chromosome factor and pericentric heterochromatic region causes lethality in hybrids between Drosophila melanogaster and its sibling species. Genetics 2012, 191, 549–559. [Google Scholar] [CrossRef] [Green Version]
  66. Blum, J.A.; Bonaccorsi, S.; Marzullo, M.; Palumbo, V.; Yamashita, Y.M.; Barbash, D.A.; Gatti, M. The Hybrid Incompatibility Genes Lhr and Hmr Are Required for Sister Chromatid Detachment During Anaphase but Not for Centromere Function. Genetics 2017, 207, 1457–1472. [Google Scholar] [CrossRef] [Green Version]
  67. Ferree, P.M.; Prasad, S. How can satellite DNA divergence cause reproductive isolation? Let us count the chromosomal ways. Genet. Res. Int. 2012, 2012, 430136. [Google Scholar] [CrossRef] [Green Version]
  68. Kejnovsky, E.; Hobza, R.; Cermak, T.; Kubat, Z.; Vyskot, B. The role of repetitive DNA in structure and evolution of sex chromosomes in plants. Heredity 2009, 102, 533–541. [Google Scholar] [CrossRef] [Green Version]
  69. Matsunaga, S. Junk DNA promotes sex chromosome evolution. Heredity 2009, 102, 525–526. [Google Scholar] [CrossRef]
  70. Deakin, J.E.; Potter, S.; O’Neill, R.; Ruiz-Herrera, A.; Cioffi, M.B.; Eldridge, M.D.B.; Fukui, K.; Marshall Graves, J.A.; Griffin, D.; Grutzner, F.; et al. Chromosomics: Bridging the Gap between Genomes and Chromosomes. Genes 2019, 10. [Google Scholar] [CrossRef] [Green Version]
  71. Dion-Cote, A.M.; Barbash, D.A. Beyond speciation genes: An overview of genome stability in evolution and speciation. Curr. Opin. Genet. Dev. 2017, 47, 17–23. [Google Scholar] [CrossRef]
  72. Ishishita, S.; Tsuboi, K.; Ohishi, N.; Tsuchiya, K.; Matsuda, Y. Abnormal pairing of X and Y sex chromosomes during meiosis I in interspecific hybrids of Phodopus campbelli and P. sungorus. Sci. Rep. 2015, 5, 9435. [Google Scholar] [CrossRef] [Green Version]
  73. Bikchurina, T.I.; Tishakova, K.V.; Kizilova, E.A.; Romanenko, S.A.; Serdyukova, N.A.; Torgasheva, A.A.; Borodin, P.M. Chromosome Synapsis and Recombination in Male-Sterile and Female-Fertile Interspecies Hybrids of the Dwarf Hamsters (Phodopus, Cricetidae). Genes 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  74. Bhattacharyya, T.; Gregorova, S.; Mihola, O.; Anger, M.; Sebestova, J.; Denny, P.; Simecek, P.; Forejt, J. Mechanistic basis of infertility of mouse intersubspecific hybrids. Proc. Natl. Acad. Sci. USA 2013, 110, E468–E477. [Google Scholar] [CrossRef] [Green Version]
  75. Bhattacharyya, T.; Reifova, R.; Gregorova, S.; Simecek, P.; Gergelits, V.; Mistrik, M.; Martincova, I.; Pialek, J.; Forejt, J. X chromosome control of meiotic chromosome synapsis in mouse inter-subspecific hybrids. PLoS Genet. 2014, 10, e1004088. [Google Scholar] [CrossRef] [Green Version]
  76. Jagannathan, M.; Warsinger-Pepe, N.; Watase, G.J.; Yamashita, Y.M. Comparative Analysis of Satellite DNA in the Drosophila melanogaster Species Complex. G3 (Bethesda) 2017, 7, 693–704. [Google Scholar] [CrossRef] [Green Version]
  77. Craddock, E.M.; Gall, J.G.; Jonas, M. Hawaiian Drosophila genomes: Size variation and evolutionary expansions. Genetica 2016, 144, 107–124. [Google Scholar] [CrossRef] [PubMed]
  78. De Lima, L.G.; Svartman, M.; Kuhn, G.C.S. Dissecting the Satellite DNA Landscape in Three Cactophilic Drosophila Sequenced Genomes. G3 (Bethesda) 2017, 7, 2831–2843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Wei, K.H.; Grenier, J.K.; Barbash, D.A.; Clark, A.G. Correlated variation and population differentiation in satellite DNA abundance among lines of Drosophila melanogaster. Proc. Natl. Acad. Sci. USA 2014, 111, 18793–18798. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Talbert, P.B.; Kasinathan, S.; Henikoff, S. Simple and Complex Centromeric Satellites in Drosophila Sibling Species. Genetics 2018, 208, 977–990. [Google Scholar] [CrossRef] [Green Version]
  81. Lemos, B.; Branco, A.T.; Hartl, D.L. Epigenetic effects of polymorphic Y chromosomes modulate chromatin components, immune response, and sexual conflict. Proc. Natl. Acad. Sci. USA 2010, 107, 15826–15831. [Google Scholar] [CrossRef] [Green Version]
  82. Francisco, F.O.; Lemos, B. How do y-chromosomes modulate genome-wide epigenetic states: Genome folding, chromatin sinks, and gene expression. J. Genom. 2014, 2, 94–103. [Google Scholar] [CrossRef] [Green Version]
  83. Sochorova, J.; Garcia, S.; Galvez, F.; Symonova, R.; Kovarik, A. Evolutionary trends in animal ribosomal DNA loci: Introduction to a new online database. Chromosoma 2018, 127, 141–150. [Google Scholar] [CrossRef] [Green Version]
  84. Larson, K.; Yan, S.J.; Tsurumi, A.; Liu, J.; Zhou, J.; Gaur, K.; Guo, D.; Eickbush, T.H.; Li, W.X. Heterochromatin formation promotes longevity and represses ribosomal RNA synthesis. PLoS Genet. 2012, 8, e1002473. [Google Scholar] [CrossRef] [Green Version]
  85. Lu, K.L.; Nelson, J.O.; Watase, G.J.; Warsinger-Pepe, N.; Yamashita, Y.M. Transgenerational dynamics of rDNA copy number in Drosophila male germline stem cells. Elife 2018, 7. [Google Scholar] [CrossRef]
  86. Chestkov, I.V.; Jestkova, E.M.; Ershova, E.S.; Golimbet, V.E.; Lezheiko, T.V.; Kolesina, N.Y.; Porokhovnik, L.N.; Lyapunova, N.A.; Izhevskaya, V.L.; Kutsev, S.I.; et al. Abundance of ribosomal RNA gene copies in the genomes of schizophrenia patients. Schizophr. Res. 2018, 197, 305–314. [Google Scholar] [CrossRef]
  87. Gibbons, J.G.; Branco, A.T.; Yu, S.K.; Lemos, B. Ribosomal DNA copy number is coupled with gene expression variation and mitochondrial abundance in humans. Nat. Commun. 2014, 5, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  88. Dion-Cote, A.M.; Symonova, R.; Rab, P.; Bernatchez, L. Reproductive isolation in a nascent species pair is associated with aneuploidy in hybrid offspring. Proc. R. Soc. B Biol. Sci. 2015, 282, 20142826. [Google Scholar] [CrossRef] [PubMed]
  89. Dion-Cote, A.M.; Symonova, R.; Lamaze, F.C.; Pelikanova, S.; Rab, P.; Bernatchez, L. Standing chromosomal variation in Lake Whitefish species pairs: The role of historical contingency and relevance for speciation. Mol. Ecol. 2017, 26, 178–192. [Google Scholar] [CrossRef] [PubMed]
  90. Symonova, R. Integrative rDNAomics-Importance of the Oldest Repetitive Fraction of the Eukaryote Genome. Genes 2019, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 11 00327 g001 550
Figure 1. Comparison of the X heterochromatin structure between mitotic and polytene chromosomes. (A) Multicolor fluorescence in situ hybridization (FISH) of a C0t1 fraction of repetitive DNA and 18S rDNA with mitotic chromosomes of An. gambiae PIMPERENA; scale bar is 2 µm. (B) FISH mapping of BAC clones 05F01 (yellow), 179F22 (red), and 01K23 (green) to the distal heterochromatin block on the X chromosome of the An. coluzzii MOPTI strain. Chromosomes are counter-stained with DAPI. (C) Relative positions of heterochromatic blocks and DNA probes of the polytene and mitotic X chromosome idiograms. rDNA = 18S rDNA probe. Repetitive DNA = C0t1 fraction of repetitive DNA. C = putative centromere. (D) Relative sizes of mitotic and polytene chromosomes in the KISUMU strain of An. gambiae. Chromosomes are counter-stained with YOYO-1. Polytene chromosome physical map is from [38].
Figure 1. Comparison of the X heterochromatin structure between mitotic and polytene chromosomes. (A) Multicolor fluorescence in situ hybridization (FISH) of a C0t1 fraction of repetitive DNA and 18S rDNA with mitotic chromosomes of An. gambiae PIMPERENA; scale bar is 2 µm. (B) FISH mapping of BAC clones 05F01 (yellow), 179F22 (red), and 01K23 (green) to the distal heterochromatin block on the X chromosome of the An. coluzzii MOPTI strain. Chromosomes are counter-stained with DAPI. (C) Relative positions of heterochromatic blocks and DNA probes of the polytene and mitotic X chromosome idiograms. rDNA = 18S rDNA probe. Repetitive DNA = C0t1 fraction of repetitive DNA. C = putative centromere. (D) Relative sizes of mitotic and polytene chromosomes in the KISUMU strain of An. gambiae. Chromosomes are counter-stained with YOYO-1. Polytene chromosome physical map is from [38].
Genes 11 00327 g001
Genes 11 00327 g002 550
Figure 2. Geographical distribution of species and phylogeny of the An. gambiae complex. (A) A map of Africa with approximal distribution of species and places of origin of the laboratory strains. (B) Species phylogeny based on the X chromosome genomic sequences (redrawn from [18]). Times of species divergence in million years (Myrs) are shown at the tree nodes.
Figure 2. Geographical distribution of species and phylogeny of the An. gambiae complex. (A) A map of Africa with approximal distribution of species and places of origin of the laboratory strains. (B) Species phylogeny based on the X chromosome genomic sequences (redrawn from [18]). Times of species divergence in million years (Myrs) are shown at the tree nodes.
Genes 11 00327 g002
Genes 11 00327 g003 550
Figure 3. Metaphase karyotypes of females in species from the An. gambiae complex. (A) An. gambiae strains. (B) An. coluzzii strains. (C) An. arabiensis, An. quadriannulatus, and An. merus strains. Black and dark gray blocks correspond to compact and diffuse heterochromatin, respectively. X chromosomes are labeled. Species names are indicated in italics and strain names are capitalized.
Figure 3. Metaphase karyotypes of females in species from the An. gambiae complex. (A) An. gambiae strains. (B) An. coluzzii strains. (C) An. arabiensis, An. quadriannulatus, and An. merus strains. Black and dark gray blocks correspond to compact and diffuse heterochromatin, respectively. X chromosomes are labeled. Species names are indicated in italics and strain names are capitalized.
Genes 11 00327 g003
Genes 11 00327 g004 550
Figure 4. Mapping of repetitive DNA elements to mitotic chromosomes of species from the An. gambiae complex. (A) FISH of AgY477–AgY53B (green) and 18S rDNA (red) in An. gambiae PIMPERENA. (B) FISH of satellite AgY53A (green) and 18S rDNA (red) in An. gambiae PIMPERENA. (C) FISH of AgY477–AgY53B (green) and 18D rDNA (red) in An. merus. (D) FISH of AgY477–AgY53B (green) and 18S rDNA (red) in An. coluzzii SUA. (E) FISH of satellite Ag53C (yellow) and 18D rDNA (red) in An. coluzzii MOPTI. (F) FISH of 18S rDNA (red) in An. quadriannulatus. (G) FISH of satellite AgY477 (red) and AgY477–AgY53B (green) in An. arabiensis. (H) FISH of satellite AgY477 (red) and 18D rDNA (green) in An. arabiensis. (I) FISH of 18S rDNA (green) and AgY477–AgY53B (red) in An. arabiensis. Scale bar = 2 µm.
Figure 4. Mapping of repetitive DNA elements to mitotic chromosomes of species from the An. gambiae complex. (A) FISH of AgY477–AgY53B (green) and 18S rDNA (red) in An. gambiae PIMPERENA. (B) FISH of satellite AgY53A (green) and 18S rDNA (red) in An. gambiae PIMPERENA. (C) FISH of AgY477–AgY53B (green) and 18D rDNA (red) in An. merus. (D) FISH of AgY477–AgY53B (green) and 18S rDNA (red) in An. coluzzii SUA. (E) FISH of satellite Ag53C (yellow) and 18D rDNA (red) in An. coluzzii MOPTI. (F) FISH of 18S rDNA (red) in An. quadriannulatus. (G) FISH of satellite AgY477 (red) and AgY477–AgY53B (green) in An. arabiensis. (H) FISH of satellite AgY477 (red) and 18D rDNA (green) in An. arabiensis. (I) FISH of 18S rDNA (green) and AgY477–AgY53B (red) in An. arabiensis. Scale bar = 2 µm.
Genes 11 00327 g004
Genes 11 00327 g005 550
Figure 5. Variation in the pattern of heterochromatin blocks, rDNA loci, and satDNA location on the X chromosomes in F1 female hybrids between An. gambiae and An. coluzzii. (A) FISH of AgY477–AgY53B (green) and AgY477 (red) in F1 ♀An. gambiae PIMPERENA × ♂An. coluzzii MOPTI. (B) FISH of AgY477 (red) and 18S rDNA (green) in F1 ♀An. gambiae ZANU × ♂An. coluzzii MALI. (C) FISH of AgY477 (red) and 18S rDNA (green) in F1 ♀An. coluzzii MOPTI × ♂An. gambiae KISUMU. Scale bar = 2 µm.
Figure 5. Variation in the pattern of heterochromatin blocks, rDNA loci, and satDNA location on the X chromosomes in F1 female hybrids between An. gambiae and An. coluzzii. (A) FISH of AgY477–AgY53B (green) and AgY477 (red) in F1 ♀An. gambiae PIMPERENA × ♂An. coluzzii MOPTI. (B) FISH of AgY477 (red) and 18S rDNA (green) in F1 ♀An. gambiae ZANU × ♂An. coluzzii MALI. (C) FISH of AgY477 (red) and 18S rDNA (green) in F1 ♀An. coluzzii MOPTI × ♂An. gambiae KISUMU. Scale bar = 2 µm.
Genes 11 00327 g005
Genes 11 00327 g006 550
Figure 6. Polymorphism of the X heterochromatin structure among the An. gambiae and An. coluzzii strains. (A) Distance and order frequency of FISH signal and DAPI band intensity peaks in An. gambiae PIMPERENA and An. coluzzii MOPTI. (B) Distribution patterns of signal intensities for satellites and DAPI along the X chromosomes in An. gambiae PIMPERENA and An. coluzzii MOPTI. (C) Relative positions of satellite peak distance and width with respect to the proximal heterochromatin band for six strains of An. gambiae and An. coluzzii. (D) Clustering of the An. gambiae and An. coluzzii strains based on the relative positions of satellite peak distance and width with respect to the proximal heterochromatin band.
Figure 6. Polymorphism of the X heterochromatin structure among the An. gambiae and An. coluzzii strains. (A) Distance and order frequency of FISH signal and DAPI band intensity peaks in An. gambiae PIMPERENA and An. coluzzii MOPTI. (B) Distribution patterns of signal intensities for satellites and DAPI along the X chromosomes in An. gambiae PIMPERENA and An. coluzzii MOPTI. (C) Relative positions of satellite peak distance and width with respect to the proximal heterochromatin band for six strains of An. gambiae and An. coluzzii. (D) Clustering of the An. gambiae and An. coluzzii strains based on the relative positions of satellite peak distance and width with respect to the proximal heterochromatin band.
Genes 11 00327 g006
Genes 11 00327 g007 550
Figure 7. Lengths of the metaphase chromosomes in species from the An. gambiae complex.
Figure 7. Lengths of the metaphase chromosomes in species from the An. gambiae complex.
Genes 11 00327 g007
Genes 11 00327 g008 550
Figure 8. Idiograms of metaphase karyotypes of species from the An. gambiae complex. (A) An idiogram of metaphase karyotypes of An. gambiae and An. coluzzii. (B) Comparison of the X chromosome idiograms among species from the An. gambiae complex. Left and right arms are labeled with L and R, respectively. Black bands correspond to condensed heterochromatin. Dark gray bands correspond to diffuse heterochromatin. White areas represent the rDNA locus. Light gray areas show euchromatin. The area of constriction represents the putative centromere. Mapping of the BAC clones and satDNA repeats AgY53A and AgY53C were performed only with An. coluzzii or An. gambiae. Two different positions of AgY477–AgY53B on the X chromosome of An. gambiae and An. coluzzii represent polymorphism among strains.
Figure 8. Idiograms of metaphase karyotypes of species from the An. gambiae complex. (A) An idiogram of metaphase karyotypes of An. gambiae and An. coluzzii. (B) Comparison of the X chromosome idiograms among species from the An. gambiae complex. Left and right arms are labeled with L and R, respectively. Black bands correspond to condensed heterochromatin. Dark gray bands correspond to diffuse heterochromatin. White areas represent the rDNA locus. Light gray areas show euchromatin. The area of constriction represents the putative centromere. Mapping of the BAC clones and satDNA repeats AgY53A and AgY53C were performed only with An. coluzzii or An. gambiae. Two different positions of AgY477–AgY53B on the X chromosome of An. gambiae and An. coluzzii represent polymorphism among strains.
Genes 11 00327 g008

Share and Cite

MDPI and ACS Style

Sharma, A.; Kinney, N.A.; Timoshevskiy, V.A.; Sharakhova, M.V.; Sharakhov, I.V. Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex. Genes 2020, 11, 327. https://doi.org/10.3390/genes11030327

AMA Style

Sharma A, Kinney NA, Timoshevskiy VA, Sharakhova MV, Sharakhov IV. Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex. Genes. 2020; 11(3):327. https://doi.org/10.3390/genes11030327

Chicago/Turabian Style

Sharma, Atashi, Nicholas A. Kinney, Vladimir A. Timoshevskiy, Maria V. Sharakhova, and Igor V. Sharakhov. 2020. "Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex" Genes 11, no. 3: 327. https://doi.org/10.3390/genes11030327

APA Style

Sharma, A., Kinney, N. A., Timoshevskiy, V. A., Sharakhova, M. V., & Sharakhov, I. V. (2020). Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex. Genes, 11(3), 327. https://doi.org/10.3390/genes11030327

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop