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.

Piwi interacts with chromatin at nuclear pores and promiscuously binds nuclear transcripts in Drosophila ovarian somatic cells.

Piwi in a complex with Piwi-interacting RNAs (piRNAs) triggers transcriptional silencing of transposable elements (TEs) in Drosophila ovaries, thus ensuring genome stability. To do this, Piwi must scan the nascent transcripts of genes and TEs for complementarity to piRNAs. The mechanism of this scanning is currently unknown. Here we report the DamID-seq mapping of multiple Piwi-interacting chromosomal domains in somatic cells of Drosophila ovaries. These domains significantly overlap with genomic regions tethered to Nuclear Pore Complexes (NPCs). Accordingly, Piwi was coimmunoprecipitated with the component of NPCs Elys and with the Xmas-2 subunit of RNA transcription and export complex, known to interact with NPCs. However, only a small Piwi fraction has transient access to DNA at nuclear pores. Importantly, although 36% of the protein-coding genes overlap with Piwi-interacting domains and RNA-immunoprecipitation results demonstrate promiscuous Piwi binding to numerous genic and TE nuclear transcripts, according to available data Piwi does not silence these genes, likely due to the absence of perfect base-pairing between piRNAs and their transcripts.

Extended disordered regions of ribosome-associated NAC proteins paralogs belong only to the germline in Drosophila melanogaster.

The nascent polypeptide-associated complex (NAC) consisting of α- and ß-subunits is an essential ribosome-associated protein conserved in eukaryotes. NAC is a ubiquitously expressed co-translational regulator of nascent protein folding and sorting providing for homeostasis of cellular proteins. Here we report on discovering the germline-specific NACαß paralogs (gNACs), whose ß-subunits, non-distinguishable by ordinary immunodetection, are encoded by five highly homologous gene copies, while the α-subunit is encoded by a single αNAC gene. The gNAC expression is detected in the primordial embryonic and adult gonads via immunostaining. The germline-specific α and ß subunits differ from the ubiquitously expressed paralogs by the extended intrinsically disordered regions (IDRs) acquired at the N- and C-termini of the coding regions, predicted to be phosphorylated. The presence of distinct phosphorylated isoforms of gNAC-ß subunits is confirmed by comparing of their profiles by 2D-isoeletrofocusing resolution before and after phosphatase treatment of testis ribosomes. We revealed that the predicted S/T sites of phosphorylation in the individual orthologous IDRs of gNAC-ß sequences of Drosophila species are positionally conserved despite these disordered regions are drastically different. We propose the IDR-dependent molecular crowding and specific coordination of NAC and other proteostasis regulatory factors at the ribosomes of germinal cells. Our findings imply that there may be a functional crosstalk between the germinal and ubiquitous α- and ß-subunits based on assessing their depletion effects on the fly viability and gonad development.

Evidence for post-transcriptional regulation of clustered microRNAs in Drosophila.

BACKGROUND: MicroRNAs (miRNA) are short 21-23nt RNAs capable of inhibiting translation of complementary target messenger RNAs. Almost half of D. melanogaster miRNA genes are grouped in genomic clusters. RESULTS: The peculiarities of the expression of clustered miRNAs were studied using publicly available libraries of sequenced small RNAs from different Drosophila tissues. We have shown that although miRNAs from almost all clusters have similar tissue expression profiles (coordinated clusters), some clusters contain miRNAs with uncoordinated expression profiles. The predicted transcription start sites (TSSs) of such clusters are located upstream of the first miRNA, but no TSSs are found within the clusters. The expression profiles of miR and miR* sequences in uncoordinated clustered miRNAs do not correlate while their profiles from the coordinated clustered miRNAs are similar. CONCLUSIONS: The presence of exclusively upstream promoters in miRNA clusters containing uncoordinated miRNAs means that the clusters are transcribed as single transcription units. The difference of tissue expression profiles of uncoordinated miRNAs and the corresponding miRs* suggests a post-transcriptional regulation of their processing or stability.

Drosophila Model for the Analysis of Genesis of LIM-kinase 1-Dependent Williams-Beuren Syndrome Cognitive Phenotypes: INDELs, Transposable Elements of the Tc1/Mariner Superfamily and MicroRNAs.

Genomic disorders, the syndromes with multiple manifestations, may occur sporadically due to unequal recombination in chromosomal regions with specific architecture. Therefore, each patient may carry an individual structural variant of DNA sequence (SV) with small insertions and deletions (INDELs) sometimes less than 10 bp. The transposable elements of the Tc1/mariner superfamily are often associated with hotspots for homologous recombination involved in human genetic disorders, such as Williams Beuren Syndromes (WBS) with LIM-kinase 1-dependent cognitive defects. The Drosophila melanogaster mutant agnts3 has unusual architecture of the agnostic locus harboring LIMK1: it is a hotspot of chromosome breaks, ectopic contacts, underreplication, and recombination. Here, we present the analysis of LIMK1-containing locus sequencing data in agnts3 and three D. melanogaster wild-type strains-Canton-S, Berlin, and Oregon-R. We found multiple strain-specific SVs, namely, single base changes and small INDEls. The specific feature of agnts3 is 28 bp A/T-rich insertion in intron 1 of LIMK1 and the insertion of mobile S-element from Tc1/mariner superfamily residing ~460 bp downstream LIMK1 3'UTR. Neither of SVs leads to amino acid substitutions in agnts3 LIMK1. However, they apparently affect the nucleosome distribution, non-canonical DNA structure formation and transcriptional factors binding. Interestingly, the overall expression of miRNAs including the biomarkers for human neurological diseases, is drastically reduced in agnts3 relative to the wild-type strains. Thus, LIMK1 DNA structure per se, as well as the pronounced changes in total miRNAs profile, probably lead to LIMK1 dysregulation and complex behavioral dysfunctions observed in agnts3 making this mutant a simple plausible Drosophila model for WBS.

RNA helicase Spn-E is required to maintain Aub and AGO3 protein levels for piRNA silencing in the germline of Drosophila.

Germline-specific RNA helicase Spindle-E (Spn-E) is known to be essential for piRNA silencing in Drosophila that takes place mainly in the perinuclear nuage granules. Loss-of-function spn-E mutations lead to tandem Stellate genes derepression in the testes and retrotransposon mobilization in the ovaries. However, Spn-E functions in the piRNA pathway are still obscure. Analysis of total library of short RNAs from the testes of spn-E heterozygous flies revealed the presence of abundant piRNA ping-pong pairs originating from Su(Ste) transcripts. The abundance of these ping-pong pairs were sharply reduced in the library from the testes of spn-E mutants. Thus we found that ping-pong mechanism contributed to Su(Ste) piRNA generation in the testes. The lack of Spn-E caused a significant drop of protein levels of key ping-pong participants, Aubergine (Aub) and AGO3 proteins of PIWI subfamily, in the germline of both males and females, but did not disrupt of their assembly in nuage granules. We found that observed decline of the protein expression was not caused by suppression of aub and ago3 transcription as well as total transcription, indicating possible contribution of Spn-E to post-transcriptional regulation.

Expansion and evolution of the X-linked testis specific multigene families in the melanogaster species subgroup.

The testis specific X-linked genes whose evolution is traced here in the melanogaster species subgroup are thought to undergo fast rate of diversification. The CK2ßtes and NACßtes gene families encode the diverged regulatory ß-subunits of protein kinase CK2 and the homologs of ß-subunit of nascent peptide associated complex, respectively. We annotated the CK2ßtes-like genes related to CK2ßtes family in the D. simulans and D. sechellia genomes. The ancestor CK2ßtes-like genes preserved in D. simulans and D. sechellia are considered to be intermediates in the emergence of the D. melanogaster specific Stellate genes related to the CK2ßtes family. The CK2ßtes-like genes are more similar to the unique autosomal CK2ßtes gene than to Stellates, taking into account their peculiarities of polymorphism. The formation of a variant the CK2ßtes gene Stellate in D. melanogaster as a result of illegitimate recombination between a NACßtes promoter and a distinct polymorphic variant of CK2ßtes-like ancestor copy was traced. We found a close nonrandom proximity between the dispersed defective copies of DINE-1 transposons, the members of Helitron family, and the CK2ßtes and NACßtes genes, suggesting an involvement of DINE-1 elements in duplication and amplification of these genes.

A novel organelle, the piNG-body, in the nuage of Drosophila male germ cells is associated with piRNA-mediated gene silencing.

Proteins of the PIWI subfamily Aub and AGO3 associated with the germline-specific perinuclear granules (nuage) are involved in the silencing of retrotransposons and other selfish repetitive elements in the Drosophila genome. PIWI proteins and their 25- to 30-nt PIWI-interacting RNA (piRNAs) are considered as key participants of the piRNA pathway. Using immunostaining, we found a large, nuage-associated organelle in the testes, the piNG-body (piRNA nuage giant body), which was significantly more massive than an ordinary nuage granule. This body contains known ovarian nuage proteins, including Vasa, Aub, AGO3, Tud, Spn-E, Bel, Squ, and Cuff, as well as AGO1, the key component of the microRNA pathway. piNG-bodies emerge at the primary spermatocyte stage of spermatogenesis during the period of active transcription. Aub, Vasa, and Tud are located at the periphery of the piNG-body, whereas AGO3 is found in its core. Mutational analysis revealed that Vasa, Aub, and AGO3 were crucial for both the maintenance of the piNG-body structure and the silencing of selfish Stellate repeats. The piNG-body destruction caused by csul mutations that abolish specific posttranslational symmetrical arginine methylation of PIWI proteins is accompanied by strong derepression of Stellate genes known to be silenced via the piRNA pathway.