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.

The mosquito family Culicidae is divided into 2 subfamilies named the Culicinae and Anophelinae. Nix, the dominant male-determining factor, has only been found in the culicines Aedes aegypti and Aedes albopictus, 2 important arboviral vectors that belong to the subgenus Stegomyia. Here we performed sex-specific whole-genome sequencing and RNAseq of divergent mosquito species and explored additional male-inclusive datasets to investigate the distribution of Nix. Except for the Culex genus, Nix homologs were found in all species surveyed from the Culicinae subfamily, including 12 additional species from 3 highly divergent tribes comprising 4 genera, suggesting Nix originated at least 133 to 165 million years ago (MYA). Heterologous expression of 1 of 3 divergent Nix open reading frames (ORFs) in Ae. aegypti resulted in partial masculinization of genetic females as evidenced by morphology and doublesex splicing. Phylogenetic analysis suggests Nix is related to femaleless (fle), a recently described intermediate sex-determining factor found exclusively in anopheline mosquitoes. Nix from all species has a conserved structure, including 3 RNA-recognition motifs (RRMs), as does fle. However, Nix has evolved at a much faster rate than fle. The RRM3 of both Nix and fle are distantly related to the single RRM of a widely distributed and conserved splicing factor transformer-2 (tra2). The RRM3-based phylogenetic analysis suggests this domain in Nix and fle may have evolved from tra2 or a tra2-related gene in a common ancestor of mosquitoes. Our results provide insights into the evolution of sex determination in mosquitoes and will inform broad applications of mosquito-control strategies based on manipulating sex ratios toward nonbiting males.

Aedes aegypti is a major vector of arboviruses that cause dengue, chikungunya, yellow fever, and Zika. Although recent success in reverse genetics has facilitated rapid progress in basic and applied research, integration of forward genetics with modern technologies remains challenging in this important species, as up to 47% of its chromosome is refractory to genetic mapping due to extremely low rate of recombination. Here, we report the development of a marker-assisted mapping strategy to readily screen for and genotype only the rare but informative recombinants, drastically increasing both the resolution and signal-to-noise ratio. Using marker-assisted mapping, we mapped a transgene that was inserted in a >100-Mb recombination desert and a sex-linked spontaneous red-eye (re) mutation just outside the region. We subsequently determined, by CRISPR/Cas9-mediated knockout, that cardinal is the causal gene of re, which is the first forward genetic identification of a causal gene in Ae. aegypti. The identification of the causal gene of the sex-linked re mutation provides the molecular foundation for using gene editing to develop versatile and stable genetic sexing methods. To facilitate genome-wide forward genetics in Ae. aegypti, we generated and compiled a number of lines with markers throughout the genome. Thus, by overcoming the challenges presented by the vast recombination deserts and the scarcity of markers, we have shown that effective forward genetic analysis is increasingly feasible in this important arboviral vector species.

Skeletal muscle from meat-producing livestock such as cattle is a major source of food for humans. To improve skeletal muscle growth efficiency or quality in cattle, it is necessary to understand the genetic and physiological mechanisms that govern skeletal muscle composition, development, and growth. Satellite cells are the myogenic progenitor cells in postnatal skeletal muscle. In this study we analyzed the composition of bovine satellite cells with single-cell RNA sequencing (scRNA-seq). We isolated satellite cells from a 2-week-old male calf, cultured them in growth medium for a week, and performed scRNA-seq using the 10x Genomics platform. Deep sequencing of two scRNA-seq libraries constructed from cultured bovine satellite cells yielded 860 million reads. Cell calling analyses revealed that these reads were sequenced from 19,096 individual cells. Clustering analyses indicated that these reads represented 15 cell clusters that differed in gene expression profile. Based on the enriched expression of markers of satellite cells (PAX7 and PAX3), markers of myoblasts (MYOD1, MYF5), and markers of differentiated myoblasts or myocytes (MYOG), three clusters were determined to be satellite cells, two clusters myoblasts, and two clusters myocytes. Gene ontology and trajectory inference analyses indicated that cells in these myogenic clusters differed in proliferation rate and differentiation stage. Two of the remaining clusters were enriched with PDGFRA, a marker of fibro-adipogenic (FAP) cells, the progenitor cells for intramuscular fat, and are therefore considered to be FAP cells. Gene ontology analyses indicated active lipogenesis in one of these two clusters. The identity of the remaining six clusters could not be defined. Overall, the results of this study support the hypothesis that bovine satellite cells are composed of subpopulations that differ in transcriptional and myogenic state. The results of this study also support the hypothesis that intramuscular fat in cattle originates from fibro-adipogenic cells.

Chromosome level assemblies are accumulating in various taxonomic groups including mosquitoes. However, even in the few reference-quality mosquito assemblies, a significant portion of the heterochromatic regions including telomeres remain unresolved. Here we produce ade novoassembly of the New World malaria mosquito,Anopheles albimanusby integrating Oxford Nanopore sequencing, Illumina, Hi-C and optical mapping. This 172.6 Mbps female assembly, which we call AalbS3, is obtained by scaffolding polished large contigs (contig N50 = 13.7 Mbps) into three chromosomes. All chromosome arms end with telomeric repeats, which is the first in mosquito assemblies and represents a significant step toward the completion of a genome assembly. These telomeres consist of tandem repeats of a novel 30-32 bp Telomeric Repeat Unit (TRU) and are confirmed by analyzing the termini of long reads and through both chromosomalin situhybridization and aBal31 sensitivity assay. The AalbS3 assembly included previously uncharacterized centromeric and rDNA clusters and more than doubled the content of transposable elements and other repetitive sequences. This telomere-to-telomere assembly, although still containing gaps, represents a significant step toward resolving biologically important but previously hidden genomic components. The comparison of different scaffolding methods will also inform future efforts to obtain reference-quality genomes for other mosquito species.

A dominant male-determining locus (M-locus) establishes the male sex (M/m) in the yellow fever mosquito, Aedes aegypti. Nix, a gene in the M-locus, was shown to be a male-determining factor (M factor) as somatic knockout of Nix led to feminized males (M/m) while transient expression of Nix resulted in partially masculin-ized females (m/m), with male reproductive organs but retained female antennae. It was not clear whether any of the other 29 genes in the 1.3-Mb M-locus are also needed for complete sex -conversion. Here, we report the generation of multiple transgenic lines that express Nix under the control of its own promoter. Ge-netic and molecular analyses of these lines provided insights un-attainable from previous transient experiments. We show that the Nix transgene alone, in the absence of the M-locus, was sufficient to convert females into males with all male-specific sexually dimorphic features and male-like gene expression. The converted m/m males are flightless, unable to perform the nuptial flight required for mat-ing. However, they were able to father sex-converted progeny when presented with cold-anesthetized wild-type females. We show that myo-sex, a myosin heavy-chain gene also in the M-locus, was re-quired for male flight as knockout of myo-sex rendered wild-type males flightless. We also show that Nix-mediated female-to-male conversion was 100% penetrant and stable over many genera-tions. Therefore, Nix has great potential for developing mosquito control strategies to reduce vector populations by female-to-male sex conversion, or to aid in a sterile insect technique that requires releasing only non-biting males.

A dominant and hemizygous male-determining locus (M locus) establishes the male sex (M/m) in the yellow fever mosquito, Aedes aegypti . Nix is a male-determining factor (M factor) in the M locus and its transient expression in females (m/m) results in partial masculinization. Here, we show that the Nix transgene alone was sufficient to convert females into fertile males, which continued to produce sex-converted progeny in subsequent generations. However, assisted mating with wild-type females was necessary, as the converted m/m males could not fly. Knockout of myo-sex , a myosin heavy chain gene and the only other protein-coding gene reported in the M locus, rendered wild-type males flightless. Thus, Nix alone converts female Ae. aegypti to fertile males and myo-sex is required for male flight. Only female Ae. aegypti mosquitoes bite and transmit disease-causing viruses. Nix -mediated female-to-male conversion is 100% penetrant and stable over many generations, indicating great potential for mosquito control.

We previously showed that Guy1, a primary signal expressed from the Y chromosome, is a strong candidate for a male-determining factor that confers female-specific lethality in Anopheles stephensi (Criscione et al., 2016). Here, we present evidence that Guy1 increases X gene expression in Guy1-transgenic females from two independent lines, providing a mechanism underlying the Guy1-conferred female lethality. The median level gene expression (MGE) of X-linked genes is significantly higher than autosomal genes in Guy1-transgenic females while there is no significant difference in MGE between X and autosomal genes in wild-type females. Furthermore, Guy1 significantly upregulates at least 40% of the 996 genes across the X chromosome in transgenic females. Guy1-conferred female-specific lethality is remarkably stable and completely penetrant. These findings indicate that Guy1 regulates dosage compensation in An. stephensi and components of dosage compensation may be explored to develop novel strategies to control mosquito-borne diseases.

A dominant and hemizygous male-determining locus (M locus) establishes the male sex (M/m) in the yellow fever mosquito,Aedes aegypti.Nixis a male-determining factor (M factor) in the M locus and its transient expression in females (m/m) results in partial masculinization. Here, we show that theNixtransgene alone was sufficient to convert females into fertile males, which continued to produce sex-converted progeny in subsequent generations. However, assisted mating with wild-type females was necessary, as the converted m/m males could not fly. Knockout ofmyo-sex, a myosin heavy chain gene and the only other protein-coding gene reported in the M locus, rendered wild-type males flightless. Thus,Nixalone converts femaleAe. aegyptito fertile males andmyo-sexis required for male flight. Only femaleAe. aegyptimosquitoes bite and transmit disease-causing viruses.Nix-mediated female-to-male conversion is 100% penetrant and stable over many generations, indicating great potential for mosquito control.

Ticks transmit more pathogens to humans and animals than any other arthropod. We describe the 2.1 Gbp nuclear genome of the tick, Ixodes scapularis (Say), which vectors pathogens that cause Lyme disease, human granulocytic anaplasmosis, babesiosis and other diseases. The large genome reflects accumulation of repetitive DNA, new lineages of retro-transposons, and gene architecture patterns resembling ancient metazoans rather than pancrustaceans. Annotation of scaffolds representing ∼57% of the genome, reveals 20,486 protein-coding genes and expansions of gene families associated with tick–host interactions. We report insights from genome analyses into parasitic processes unique to ticks, including host ‘questing’, prolonged feeding, cuticle synthesis, blood meal concentration, novel methods of haemoglobin digestion, haem detoxification, vitellogenesis and prolonged off-host survival. We identify proteins associated with the agent of human granulocytic anaplasmosis, an emerging disease, and the encephalitis-causing Langat virus, and a population structure correlated to life-history traits and transmission of the Lyme disease agent.

Despite their importance in sexual differentiation and reproduction, Y chromosome genes are rarely described because they reside in repeat-rich regions that are difficult to study. Here, we show that Guy1, a unique Y chromosome gene of a major urban malaria mosquito Anopheles stephensi, confers 100% female lethality when placed on the autosomes. We show that the small GUY1 protein (56 amino acids in length) causes female lethality and that males carrying the transgene are reproductively more competitive than their non-transgenic siblings under laboratory conditions. The GUY1 protein is a primary signal from the Y chromosome that affects embryonic development in a sex-specific manner. Our results have demonstrated, for the first time in mosquitoes, the feasibility of stable transgenic manipulation of sex ratios using an endogenous gene from the male-determining chromosome. These results provide insights into the elusive M factor and suggest exciting opportunities to reduce mosquito populations and disease transmission.

Sex determination in the mosquito Aedes aegypti is governed by a dominant male-determining factor (M factor) located within a Y chromosome-like region called the M locus. Here, we show that an M-locus gene, Nix, functions as an M factor in A. aegypti. Nix exhibits persistent M linkage and early embryonic expression, two characteristics required of an M factor. Nix knockout with clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 resulted in largely feminized genetic males and the production of female isoforms of two key regulators of sexual differentiation: doublesex and fruitless. Ectopic expression of Nix resulted in genetic females with nearly complete male genitalia. Thus, Nix is both required and sufficient to initiate male development. This study provides a foundation for mosquito control strategies that convert female mosquitoes into harmless males.

Complete dosage compensation refers to hyperexpression of the entire X or Z chromosome in organisms with heterogametic sex chromosomes (XY male or ZW female) in order to compensate for having only one copy of the X or Z chromosome. Recent analyses suggest that complete dosage compensation, as in Drosophila melanogaster, may not be the norm. There has been no systematic study focusing on dosage compensation in mosquitoes. However, analysis of dosage compensation in Anopheles mosquitoes provides opportunities for evolutionary insights, as the X chromosome of Anopheles and that of its Dipteran relative, D. melanogaster formed independently from the same ancestral chromosome. Furthermore, Culicinae mosquitoes, including the Aedes genus, have homomorphic sex-determining chromosomes, negating the need for dosage compensation. Thus, Culicinae genes provide a rare phylogenetic context to investigate dosage compensation in Anopheles mosquitoes. Here, we performed RNA-seq analysis of male and female samples of the Asian malaria mosquito Anopheles stephensi and the yellow fever mosquito Aedes aegypti. Autosomal and X-linked genes in An. stephensi showed very similar levels of expression in both males and females, indicating complete dosage compensation. The uniformity of average expression levels of autosomal and X-linked genes remained when An. stephensi gene expression was normalized by that of their Ae. aegypti orthologs, strengthening the finding of complete dosage compensation in Anopheles. In addition, we comparatively analyzed the differentially expressed genes between adult males and adult females in both species, investigated sex-biased gene chromosomal distribution patterns in An. stephensi and provided three examples where gene duplications may have enabled the acquisition of sex-specific expression during mosquito evolution.

Background Anopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range. Results Here, we report the genome sequence and annotation of the Indian strain of the type form of An. stephensi . The 221 Mb genome assembly represents more than 92% of the entire genome and was produced using a combination of 454, Illumina, and PacBio sequencing. Physical mapping assigned 62% of the genome onto chromosomes, enabling chromosome-based analysis. Comparisons between An. stephensi and An. gambiae reveal that the rate of gene order reshuffling on the X chromosome was three times higher than that on the autosomes. An. stephensi has more heterochromatin in pericentric regions but less repetitive DNA in chromosome arms than An. gambiae . We also identify a number of Y-chromosome contigs and BACs. Interspersed repeats constitute 7.1% of the assembled genome while LTR retrotransposons alone comprise more than 49% of the Y contigs. RNA-seq analyses provide new insights into mosquito innate immunity, development, and sexual dimorphism. Conclusions The genome analysis described in this manuscript provides a resource and platform for fundamental and translational research into a major urban malaria vector. Chromosome-based investigations provide unique perspectives on Anopheles chromosome evolution. RNA-seq analysis and studies of immunity genes offer new insights into mosquito biology and mosquito-parasite interactions.

Anopheles stephensi is a principal vector of urban malaria on the Indian subcontinent and an emerging model for molecular and genetic studies of mosquito biology. To enhance our understanding of female mosquito reproduction, and to develop new tools for basic research and for genetic strategies to control mosquito-borne infectious diseases, we identified 79 genes that displayed previtellogenic germline-specific expression based on RNA-Seq data generated from 11 life stage–specific and sex-specific samples. Analysis of this gene set provided insights into the biology and evolution of female reproduction. Promoters from two of these candidates, vitellogenin receptor and nanos, were used in independent transgenic cassettes for the expression of artificial microRNAs against suspected mosquito maternal-effect genes, discontinuous actin hexagon and myd88. We show these promoters have early germline-specific expression and demonstrate 73% and 42% knockdown of myd88 and discontinuous actin hexagon mRNA in ovaries 48 hr after blood meal, respectively. Additionally, we demonstrate maternal-specific delivery of mRNA and protein to progeny embryos. We discuss the application of this system of maternal delivery of mRNA/miRNA/protein in research on mosquito reproduction and embryonic development, and for the development of a gene drive system based on maternal-effect dominant embryonic arrest.

Background: Y chromosomes are responsible for the initiation of male development, male fertility, and other male-related functions in diverse species. However, Y genes are rarely characterized outside a few model species due to the arduous nature of studying the repeat-rich Y.Results: The chromosome quotient (CQ) is a novel approach to systematically discover Y chromosome genes. In the CQ method, genomic DNA from males and females is sequenced independently and aligned to candidate reference sequences. The female to male ratio of the number of alignments to a reference sequence, a parameter called the chromosome quotient ( CQ), is used to determine whether the sequence is Y-linked. Using the CQ method, we successfully identified known Y sequences from Homo sapiens and Drosophila melanogaster. The CQ method facilitated the discovery of Y chromosome sequences from the malaria mosquitoes Anopheles stephensi and An. gambiae. Comparisons to transcriptome sequence data with blastn led to the discovery of six Anopheles Y genes, three from each species. All six genes are expressed in the early embryo. Two of the three An. stephensi Y genes were recently acquired from the autosomes or the X. Although An. stephensi and An. gambiae belong to the same subgenus, we found no evidence of Y genes shared between the species.Conclusions: The CQ method can reliably identify Y chromosome sequences using the ratio of alignments from male and female sequence data. The CQ method is widely applicable to species with fragmented genome assemblies produced from next-generation sequencing data. Analysis of the six Y genes characterized in this study indicates rapid Y chromosome evolution between An. stephensi and An. gambiae. The Anopheles Y genes discovered by the CQ method provide unique markers for population and phylogenetic analysis, and opportunities for novel mosquito control measures through the manipulation of sexual dimorphism and fertility.

Horizontal transfer of genetic material between complex organisms often involves transposable elements (TEs). For example, a DNA transposon mariner has been shown to undergo horizontal transfer between different orders of insects and between different phyla of animals. Here we report the discovery and characterization of an ITmD37D transposon, MJ1, in Anopheles sinensis. We show that some MJ1 elements in Aedes aegypti and An. sinensis contain intact open reading frames and share nearly 99% nucleotide identity over the entire transposon, which is unexpectedly high given that these two genera had diverged 145-200 million years ago. Chromosomal hybridization and TE-display showed that MJ1 copy number is low in An. sinensis. Among 24 mosquito species surveyed, MJ1 is only found in Ae. aegypti and the hyrcanus group of anopheline mosquitoes to which An. sinensis belongs. Phylogenetic analysis is consistent with horizontal transfer and provides the basis for inference of its timing and direction. Although report of horizontal transfer of DNA transposons between higher eukaryotes is accumulating, our analysis is one of a small number of cases in which horizontal transfer of nearly identical TEs among highly divergent species has been thoroughly investigated and strongly supported. Horizontal transfer involving mosquitoes is of particular interest because there are ongoing investigations of the possibility of spreading pathogen-resistant genes into mosquito populations to control malaria and other infectious diseases. The initial indication of horizontal transfer of MJ1 came from comparisons between a 0.4x coverage An. sinensis 454 sequence database and available TEs in mosquito genomes. Therefore we have shown that it is feasible to use low coverage sequencing to systematically uncover horizontal transfer events. Expanding such efforts across a wide range of species will generate novel insights into the relative frequency of horizontal transfer of different TEs and provide the evolutionary context of these lateral transfer events.

Culex quinquefasciatus (the southern house mosquito) is an important mosquito vector of viruses such as West Nile virus and St. Louis encephalitis virus, as well as of nematodes that cause lymphatic filariasis. C. quinquefasciatus is one species within the Culex pipiens species complex and can be found throughout tropical and temperate climates of the world. The ability of C. quinquefasciatus to take blood meals from birds, livestock, and humans contributes to its ability to vector pathogens between species. Here, we describe the genomic sequence of C. quinquefasciatus : Its repertoire of 18,883 protein-coding genes is 22% larger than that of Aedes aegypti and 52% larger than that of Anopheles gambiae with multiple gene-family expansions, including olfactory and gustatory receptors, salivary gland genes, and genes associated with xenobiotic detoxification.

Background SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE ( SINE200 ) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms. Methods A SINE -PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis , Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed. Results Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted evolution and/or mixtures of M and S IGS-sequences among the arrays of single chromatids. Conclusion The approach utilized allowed to develop new easy-to-use co-dominant markers for the analysis of genetic differentiation between M and S-forms and opens new perspectives in the study of the speciation process ongoing within A. gambiae .