The chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus uncovers patterns of genome evolution in mosquitoes

BACKGROUND

Understanding genome organization and evolution is important for species involved in transmission of human diseases, such as mosquitoes. Anophelinae and Culicinae subfamilies of mosquitoes show striking differences in genome sizes, sex chromosome arrangements, behavior, and ability to transmit pathogens. However, the genomic basis of these differences is not fully understood.

METHODS

In this study, we used a combination of advanced genome technologies such as Oxford Nanopore Technology sequencing, Hi-C scaffolding, Bionano, and cytogenetic mapping to develop an improved chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus.

RESULTS

We then used this assembly to annotate odorant receptors, odorant binding proteins, and transposable elements. A genomic region containing male-specific sequences on chromosome 1 and a polymorphic inversion on chromosome 3 were identified in the Cx. quinquefasciatus genome. In addition, the genome of Cx. quinquefasciatus was compared with the genomes of other mosquitoes such as malaria vectors An. coluzzi and An. albimanus, and the vector of arboviruses Ae. aegypti. Our work confirms significant expansion of the two chemosensory gene families in Cx. quinquefasciatus, as well as a significant increase and relocation of the transposable elements in both Cx. quinquefasciatus and Ae. aegypti relative to the Anophelines. Phylogenetic analysis clarifies the divergence time between the mosquito species. Our study provides new insights into chromosomal evolution in mosquitoes and finds that the X chromosome of Anophelinae and the sex-determining chromosome 1 of Culicinae have a significantly higher rate of evolution than autosomes.

CONCLUSION

The improved Cx. quinquefasciatus genome assembly uncovered new details of mosquito genome evolution and has the potential to speed up the development of novel vector control strategies.

Patterns of genetic differentiation imply distinct phylogeographic history of the mosquito species Anopheles messeae and Anopheles daciae in Eurasia

Detailed knowledge of phylogeography is important for control of mosquito species involved in the transmission of human infectious diseases. Anopheles messeae is a geographically widespread and genetically diverse dominant vector of malaria in Eurasia. A closely related species, An. daciae, was originally distinguished from An. messeae based on five nucleotide substitutions in its ribosomal DNA (rDNA). However, the patterns of phylogeographic history of these species in Eurasia remain poorly understood. Here, using internal transcribed spacer 2 (ITS2) of rDNA and karyotyping for the species identification we determined the composition of five Anopheles species in 28 locations in Eurasia. Based on the frequencies of 11 polymorphic chromosomal inversions used as genetic markers, a large-scale population genetics analysis was performed of 1932 mosquitoes identified as An. messeae, An. daciae and their hybrids. The largest genetic differences between the species were detected in the X sex chromosome suggesting a potential involvement of this chromosome in speciation. The frequencies of autosomal inversions in the same locations differed by 13%-45% between the species demonstrating a restricted gene flow between the species. Overall, An. messeae was identified as a diverse species with a more complex population structure than An. daciae. The clinal gradients in frequencies of chromosomal inversions were determined in both species implicating their possible involvement in climate adaptations. The frequencies of hybrids were low ~1% in northern Europe but high up to 50% in south-eastern populations. Thus, our study revealed critical differences in patterns of phylogeographic history between An. messeae and An. daciae in Eurasia. This knowledge will help to predict the potential of the malaria transmission in the northern territories of the continent.

The complex role of transcription factor GAGA in germline death during Drosophila spermatogenesis: transcriptomic and bioinformatic analyses

The GAGA protein (also known as GAF) is a transcription factor encoded by the Trl gene in D. melanogaster. GAGA is involved in the regulation of transcription of many genes at all stages of fly development and life. Recently, we investigated the participation of GAGA in spermatogenesis and discovered that Trl mutants experience massive degradation of germline cells in the testes. Trl underexpression induces autophagic death of spermatocytes, thereby leading to reduced testis size. Here, we aimed to determine the role of the transcription factor GAGA in the regulation of ectopic germline cell death. We investigated how Trl underexpression affects gene expression in the testes. We identified 15,993 genes in three biological replicates of our RNA-seq analysis and compared transcript levels between hypomorphic Trl ^R85/Trl ³⁶² and Oregon testes. A total of 2,437 differentially expressed genes were found, including 1,686 upregulated and 751 downregulated genes. At the transcriptional level, we detected the development of cellular stress in the Trl-mutant testes: downregulation of the genes normally expressed in the testes (indicating slowed or abrogated spermatocyte differentiation) and increased expression of metabolic and proteolysis-related genes, including stress response long noncoding RNAs. Nonetheless, in the Flybase Gene Ontology lists of genes related to cell death, autophagy, or stress, there was no enrichment with GAGA-binding sites. Furthermore, we did not identify any specific GAGA-dependent cell death pathway that could regulate spermatocyte death. Thus, our data suggest that GAGA deficiency in male germline cells leads to an imbalance of metabolic processes, impaired mitochondrial function, and cell death due to cellular stress.

New Cytogenetic Photomap and Molecular Diagnostics for the Cryptic Species of the Malaria Mosquitoes Anopheles messeae and Anopheles daciae from Eurasia

The Eurasian malaria vector Anopheles messeae is a widely spread and genetically diverse species. Five widespread polymorphic chromosomal inversions were found in natural populations of this mosquito. A cryptic species, Anopheles daciae, was differentiated from An. messeae by the presence of several nucleotide substitutions in the Internal Transcribed Spacer 2 (ITS2) region of ribosomal DNA. However, because of the absence of a high-quality reference cytogenetic map, the inversion polymorphisms in An. daciae and An. messeae remain poorly understood. Moreover, a recently determined heterogeneity in ITS2 in An. daciae questioned the accuracy of the previously used Restriction Fragment Length Polymorphism (RFLP) assay for species diagnostics. In this study, a standard-universal cytogenetic map was constructed based on orcein stained images of chromosomes from salivary glands for population studies of the chromosomal inversions that can be used for both An. messeae and An. daciae. In addition, a new ITS2-RFLP approach for species diagnostics was developed. Both methods were applied to characterize inversion polymorphism in populations of An. messeae and An. daciae from a single location in Western Siberia in Russia. The analysis demonstrates that cryptic species are remarkably different in their frequencies of chromosomal inversion variants. Our study supports previous observations that An. messeae has higher inversion polymorphism in all autosomes than the cryptic species An. daciae.

rs2072580T>A Polymorphism in the Overlapping Promoter Regions of the SART3 and ISCU Genes Associated with the Risk of Breast Cancer

We analyzed association of potentially regulatory polymorphisms (rs590352, rs11542583, rs3829202, rs207258, and rs4796672) with breast cancer. A significant association was found between this disease and rs2072580T>A (p=0.001) located in the overlapping promoter regions of the SART3 and ISCU genes. In women with AA and AT genotypes, the risk of breast cancer is higher by 6.7 times (p=0.001) and 12 times (p=0.001), respectively, in comparison with TT genotype. Under a codominant model of inheritance (AT vs AA+TT), the risk of breast cancer was increased by 4.2 times (р=0.001) for the AT genotype. Under a recessive model of inheritance (TT vs AA+TT), the risk of disease was 10-fold higher (р=0.001) for the TT genotype. It has been demonstrated that the T>A substitution affects the binding properties of transcription factors CREB1 and REST.

Potential regulatory SNPs in the ATXN7L3B and KRT15 genes are associated with gender-specific colorectal cancer risk

Aim: According to the current data, a major factor for phenotypic variation of complex traits and disease susceptibility is the cis-acting effects of noncoding variants on gene expression. Our purpose was to evaluate the association between colorectal cancer (CRC) and six single nucleotide polymorphisms identified using our original bioinformatics approach as regulatory and putatively related to CRC. Materials: One hundred and sixty CRC patients and 185 healthy controls have been genotyped for rs590352, rs2072580, rs78317230, rs3829202, rs11542583 and rs4796672. Results: Genotypes and alleles distributions of rs590352 of ATXN7L3B gene were significantly different between the male CRC subjects and controls. Significant correlation of genotype with CRC is observable for women only for the rs4796672 of KRT15 gene. Analysis of haplotypes reveals that rs2072580 of the ISCU and SART3 genes can be also associated with CRC. Conclusion: We have identified three SNPs associated with CRC risk and demonstrated a gender specificity of rs590352 and rs4796672.

Personalized Approach to Assessing mRNA Expression of Histidine-Rich Glycoprotein and Immunohistochemical Markers in Diseases of the Breast

Biopsy material of patients with malignant and benign breast diseases was examined. HRG mRNA expression was detected in 70% of cases in biopsy material obtained from patients with nonspecific invasive carcinoma and in 66.7% of cases in biopsy material of patients with benign breast diseases. Immunohistochemical analysis revealed expression of collagen II, the beta-1 integrin, and E-cadherin-markers of epithelial-mesenchymal transition. The use of RT-qPCR combined with immunohistochemical study made it possible to identify atypical cells, which can be regarded as precancerous changes, in individual patients.

Messenger RNA of the Histidine-Rich Glycoprotein in Breast Tumors

The content of mRNA of the histidine-rich glycoprotein (HRG), a potential marker of malignant neoplasia, which can be used in differential diagnosis of breast tumors, was determined in 110 breast tumor biopsy samples. The presence of HRG mRNA did not depend on the cancer type, on the preoperative treatment or its absence, as well as on the tumor progression stage and the presence of metastases.

[DNA diagnosis of mixed invasions of Opisthorchis felineus and Metorchis bilis by polymerase chain reaction]

Opisthorchiasis is one of the significant naturofocal diseases in Russia. The diagnosis of opisthorchiasis may mask some diseases, such as opisthorchiasis, metorchiasis, and clonorchiasis - biohelminthoses that are induced by various representatives of the family Opisthorchiidae-Opisthorchis felineus/O.viverrini, Metorchis bilis, and Clonorchis sinensis, respectively. Coproovoscopy and serologic methods fail to accurately define the species-specific affiliation; in this connection the identification of opisthorchids, by using DNA diagnostic techniques, becomes urgent. The present paper gives the results of development of DNA diagnosticum, which differentiates parasitic diseases induced by O. felineus and M. bills. The ribosomal RNA gene cluster fragment incorporating internal transcribed spacer 2 (ITS2) was used as a diagnostic marker. A system for diagnosing opisthorchiasis was developed as a multiplex polymerase chain analysis and tested on 37 patients infected with various species of opisthorchids.