Resequencing Helminth Genomes for Population and Genetic Studies

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

Genomic signatures of selection associated with benzimidazole drug treatments in Haemonchus contortus field populations

Genome-wide methods offer a powerful approach to detect signatures of drug selection. However, limited availability of suitable reference genomes and the difficulty of obtaining field populations with well-defined, distinct drug treatment histories mean there is little information on the signatures of selection in parasitic nematodes and on how best to detect them. This study addresses these knowledge gaps by using field populations of Haemonchus contortus with well-defined benzimidazole treatment histories, leveraging a recently completed chromosomal-scale reference genome assembly. We generated a panel of 49,393 genomic markers to genotype 20 individual adult worms from each of four H. contortus populations: two from closed sheep flocks with an approximate 20 year history of frequent benzimidazole treatment, and two populations with a history of little or no treatment. Sampling occurred in the same geographical region to limit genetic differentiation and maximise the detection sensitivity. A clear signature of selection was detected on chromosome I, centred on the isotype-1 β-tubulin gene. Two additional, but weaker, signatures of selection were detected; one near the middle of chromosome I spanning 3.75 Mbp and 259 annotated genes, and one on chromosome II spanning a region of 3.3 Mbp and 206 annotated genes, including the isotype-2 β-tubulin locus. We also assessed how sensitivity was impacted by sequencing depth, worm number, and pooled versus individual worm sequence data. This study provides the first known direct genome-wide evidence for any parasitic nematode, that the isotype-1 β-tubulin gene is quantitatively the single most important benzimidazole resistance locus. It also identified two additional genomic regions that likely contain benzimidazole resistance loci of secondary importance. This study provides an experimental framework to maximise the power of genome-wide approaches to detect signatures of selection driven by anthelmintic drug treatments in field populations of parasitic nematodes.

A chromosome-level genome of the human blood fluke Schistosoma japonicum identifies the genomic basis of host-switching

The evolution and adaptation of S. japonicum, a zoonotic parasite that causes human schistosomiasis, remain unclear because of the lack of whole-genome data. We construct a chromosome-level S. japonicum genome and analyze it together with 72 samples representing six populations of the entire endemic region. We observe a Taiwan zoophilic lineage splitting from zoonotic populations ∼45,000 years ago, consistent with the divergent history of their intermediate hosts. Interestingly, we detect a severe population bottleneck in S. japonicum, largely coinciding with human history in Asia during the last glacial maximum. We identify several genomic regions underlying natural selection, including GATAD2A and Lmln, both showing remarkable differentiation among different areas. RNAi knockdown suggests association of GATAD2A with parasite development and infection in definitive hosts, while Lmln relates to the specificity of the intermediate hosts. Our study provides insights into the evolution of S. japonicum and serves as a resource for further studies.

Titanium dioxide-mediated fatty acids promoted carotenoid synthesis in Phaffia rhodozyma PR106 analyzed whole genome

Carotenoid, as good colorant and antioxidant, is widely used in the fields of food, medicine and feed. The whole genome of P. rhodozyma PR106 strain with 228.77 mg/L carotenoid (mainly included astaxanthin, β-carotene and lycopene) yield was sequenced, and the genome size was 16.18 Mb, the GC content was 47%. The genetic evolution analysis indicated that PR106 greatly changed in evolution process, and closely related to P. rhodozyma CBS7918. Under 500 mg/L titanium dioxide (TiO₂) stress, carotenoid yield of PR106 was 2.15 times that of the control for 48 h, and was 305.12 mg/L in PR106 to 72 h, interestingly, the yield of oleate, linoleate and α-linolenate also increased significantly among 51 fatty acids by targeted metabolomics analysis. TiO₂ promoted carotenoid synthesis of PR106 by forming astaxanthin esters to reduce the feedback inhibition of carotenoid synthesis. These results provided a theoretical basis for carotenoid production and development using P. rhodozyma.

'Escalibur' - a practical pipeline for the de novo-analysis of nucleotide variation in non-model eukaryotes

The revolution in genomics has enabled large‐scale population genetic investigations of a wide range of organisms, but there has been a relatively limited focus on improving analytical pipelines. To efficiently analyse large data sets, highly integrated and automated software pipelines, which are easy to use, efficient, reliable, reproducible and run in multiple computational environments, are required. A number of software workflows have been developed to handle and process such data sets for population genetic analyses, but effective, specialized pipelines for genetic and statistical analyses of nonmodel organisms are lacking. For most species, resources for variomes (sets of genetic variations found in populations of species) are not available, and/or genome assemblies are often incomplete and fragmented, complicating the selection of the most suitable reference genome when multiple assemblies are available. Additionally, the biological samples used often contain extraneous DNA from sources other than the species under investigation (e.g., microbial contamination), which needs to be removed prior to genetic analyses. For these reasons, we established a new pipeline, called Escalibur, which includes: functionalities, such as data trimming and mapping; selection of a suitable reference genome; removal of contaminating read data; recalibration of base calls; and variant‐calling. Escalibur uses a proven gatk variant caller and workflow description language (WDL), and is, therefore, a highly efficient and scalable pipeline for the genome‐wide identification of nucleotide variation in eukaryotes. This pipeline is available at https://gitlab.unimelb.edu.au/bioscience/escalibur (version 0.3‐beta) and is essentially applicable to any prokaryote or eukaryote.

Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

Deep amplicon sequencing highlights low intra-host genetic variability of Echinococcus multilocularis and high prevalence of the European-type haplotypes in coyotes and red foxes in Alberta, Canada

Echinococcus multilocularis (Em) is a zoonotic parasite considered a global emergent pathogen. Recent findings indicate that the parasite is expanding its range in North America and that European-type haplotypes are circulating in western Canada. However, genetic analyses are usually conducted only on a few parasites out of thousands of individuals within each definitive host, likely underestimating the prevalence of less common haplotypes. Moreover, mixed infections with several mtDNA haplotypes in the same host have been reported, but their relative abundance within the host was never estimated. We aimed to 1) estimate the frequency of co-infections of different Em haplotypes in coyotes (Canis latrans) and red foxes (Vulpes vulpes) from western Canada and their relative abundance within the definitive hosts, 2) detect less prevalent haplotypes by sampling a larger proportion of the parasite subpopulation per host, and 3) investigate differences in the distribution of Em haplotypes in these main definitive hosts; foxes and coyotes. We extracted DNA from ~10% of the worm subpopulation per host (20 foxes and 47 coyotes) and used deep amplicon sequencing (NGS technology) on four loci, targeting the most polymorphic regions from the mitochondrial genes cox1 (814 bp), nad1 (344 bp), and cob (387 bp). We detected the presence of mixed infections with multiple Em haplotypes and with different Echinococcus species including Em and E. granulosus s.l. genotypes G8/G10, low intraspecific diversity of Em, and a higher abundance of the European-type haplotypes in both hosts. Our results suggest a population expansion of the European over the North American strain in Alberta and a limited distribution of some European-type haplotypes. Our findings indicate that deep amplicon sequencing represents a valuable tool to characterize Em in multiple hosts, to assess the current distribution and possible origins of the European strain in North America. The potential use of next-generation sequencing technologies is particularly important to understand the patterns of geographic expansion of this parasite.

Patterns of relatedness and genetic diversity inferred from whole genome sequencing of archival blood fluke miracidia (Schistosoma japonicum)

Genomic approaches hold great promise for resolving unanswered questions about transmission patterns and responses to control efforts for schistosomiasis and other neglected tropical diseases. However, the cost of generating genomic data and the challenges associated with obtaining sufficient DNA from individual schistosome larvae (miracidia) from mammalian hosts have limited the application of genomic data for studying schistosomes and other complex macroparasites. Here, we demonstrate the feasibility of utilizing whole genome amplification and sequencing (WGS) to analyze individual archival miracidia. As an example, we sequenced whole genomes of 22 miracidia from 11 human hosts representing two villages in rural Sichuan, China, and used these data to evaluate patterns of relatedness and genetic diversity. We also down-sampled our dataset to test how lower coverage sequencing could increase the cost effectiveness of WGS while maintaining power to accurately infer relatedness. Collectively, our results illustrate that population-level WGS datasets are attainable for individual miracidia and represent a powerful tool for ultimately providing insight into overall genetic diversity, parasite relatedness, and transmission patterns for better design and evaluation of disease control efforts.

Characterization of the Complete Mitochondrial Genomes of Two Sibling Species of Parasitic Roundworms, Haemonchus contortus and Teladorsagia circumcincta

Haemonchus contortus and Teladorsagia circumcincta are among the two most pathogenic internal parasitic nematodes infecting small ruminants, such as sheep and goats, and are a global animal health issue. Accurate identification and delineation of Haemonchidae species is essential for development of diagnostic and control strategies with high resolution for Trichostrongyloidea infection in ruminants. Here, we describe in detail and compare the complete mitochondrial (mt) genomes of the New Zealand H. contortus and T. circumcincta field strains to improve our understanding of species- and strain-level evolution in these closely related roundworms. In the present study, we performed extensive comparative bioinformatics analyses on the recently sequenced complete mt genomes of the New Zealand H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP field strains. Amino acid sequences inferred from individual genes of each of the two mt genomes were compared, concatenated and subjected to phylogenetic analysis using Bayesian inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP). The AT-rich mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP are 14,001 bp (A+T content of 77.4%) and 14,081 bp (A+T content of 77.3%) in size, respectively. All 36 of the typical nematode mt genes are transcribed in the forward direction in both species and comprise of 12 protein-encoding genes (PCGs), 2 ribosomal RNA (rrn) genes, and 22 transfer RNA (trn) genes. The secondary structures for the 22 trn genes and two rrn genes differ between H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP, however the gene arrangements of both are consistent with other Trichostrongylidea sequenced to date. Comparative analyses of the complete mitochondrial nucleotide sequences, PCGs, A+T rich and non-coding repeat regions of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP further reinforces the high levels of diversity and gene flow observed among Trichostrongylidea, and supports their potential as ideal markers for strain-level identification from different hosts and geographical regions with high resolution for future studies. The complete mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP presented here provide useful novel markers for further studies of the meta-population connectivity and the genetic mechanisms driving evolution in nematode species.

The first mitochondrial genomes of endosymbiotic rhabdocoels illustrate evolutionary relaxation of atp8 and genome plasticity in flatworms

The first three mitochondrial (mt) genomes of endosymbiotic turbellarian flatworms are characterised for the rhabdocoels Graffilla buccinicola, Syndesmis echinorum and S. kurakaikina. Interspecific comparison of the three newly obtained sequences and the only previously characterised rhabdocoel, the free-living species Bothromesostoma personatum, reveals high mt genomic variability, including numerous rearrangements. The first intrageneric comparison within rhabdocoels shows that gene order is not fully conserved even between congeneric species. Atp8, until recently assumed absent in flatworms, was putatively annotated in two sequences. Selection pressure was tested in a phylogenetic framework and is shown to be significantly relaxed in this and another protein-coding gene: cox1. If present, atp8 appears highly derived in platyhelminths and its functionality needs to be addressed in future research. Our findings for the first time allude to a large degree of undiscovered (mt) genomic plasticity in rhabdocoels. It merits further attention whether this variation is correlated with a symbiotic lifestyle. Our results illustrate that this phenomenon is widespread in flatworms as a whole and not exclusive to the better-studied neodermatans.

Gene duplications, divergence and recombination shape adaptive evolution of the fish ectoparasite Gyrodactylus bullatarudis

Determining the molecular basis of parasite adaptation to its host is an important component in understanding host-parasite coevolution and the epidemiology of parasitic infections. Here, we investigate short- and long-term adaptive evolution in the eukaryotic parasite Gyrodactylus bullatarudis infecting Caribbean guppies (Poecilia reticulata), by comparing the reference genome of Tobagonian G. bullatarudis with other Platyhelminthes, and by analysing resequenced samples from local Trinidadian populations. At the macroevolutionary timescale, we observed duplication of G-protein and serine proteases genes, which are probably important in host-parasite arms races. Serine protease also showed strong evidence of ongoing, diversifying selection at the microevolutionary timescale. Furthermore, our analyses revealed that a hybridization event, involving two divergent genomes, followed by recombination has dramatically affected the genetic composition of Trinidadian populations. The recombinant genotypes invaded Trinidad and replaced local parasites in all populations. We localized more than 300 genes in regions fixed in local populations for variants of different origin, possibly due to diversifying selection pressure from local host populations. In addition, around 70 genes were localized in regions identified as heterozygous in some, but not all, individuals. This pattern is consistent with a very recent spread of recombinant parasites. Overall, our results are consistent with the idea that recombination between divergent genomes can result in particularly successful parasites.

The Genome Sequence of the Anthelmintic-Susceptible New Zealand Haemonchus contortus

Internal parasitic nematodes are a global animal health issue causing drastic losses in livestock. Here, we report a H. contortus representative draft genome to serve as a genetic resource to the scientific community and support future experimental research of molecular mechanisms in related parasites. A de novo hybrid assembly was generated from PCR-free whole genome sequence data, resulting in a chromosome-level assembly that is 465 Mb in size encoding 22,341 genes. The genome sequence presented here is consistent with the genome architecture of the existing Haemonchus species and is a valuable resource for future studies regarding population genetic structures of parasitic nematodes. Additionally, comparative pan-genomics with other species of economically important parasitic nematodes have revealed highly open genomes and strong collinearities within the phylum Nematoda.

Spirorchiidiasis in marine turtles: the current state of knowledge

Blood flukes of the family Spirorchiidae are important disease agents in marine turtles. The family is near cosmopolitan in distribution. Twenty-nine marine species across 10 genera are currently recognized, but taxonomic problems remain and it is likely that more species will be discovered. Spirorchiids infect the circulatory system, where they and their eggs cause a range of inflammatory lesions. Infection is sometimes implicated in the death of the turtle. In some regions, prevalence in stranded turtles is close to 100%. Knowledge of life cycles, important for control and epidemiological studies, has proven elusive until recently, when the first intermediate host identifications were made. Recent molecular studies of eggs and adult worms indicate that a considerable level of intrageneric and intraspecific diversity exists. The characterization of this diversity is likely to be of importance in exploring parasite taxonomy and ecology, unravelling life cycles, identifying the differential pathogenicity of genotypes and species, and developing antemortem diagnostic tools, all of which are major priorities for future spirorchiid research. Diagnosis to date has been reliant on copromicroscopy or necropsy, which both have significant limitations. The current lack of reliable antemortem diagnostic options is a roadblock to determining the true prevalence and epidemiology of spirorchiidiasis and the development of effective treatment regimes.

The population genetics of parasitic nematodes of wild animals

Parasitic nematodes are highly diverse and common, infecting virtually all animal species, and the importance of their roles in natural ecosystems is increasingly becoming apparent. How genes flow within and among populations of these parasites - their population genetics - has profound implications for the epidemiology of host infection and disease, and for the response of parasite populations to selection pressures. The population genetics of nematode parasites of wild animals may have consequences for host conservation, or influence the risk of zoonotic disease. Host movement has long been recognised as an important determinant of parasitic nematode population genetic structure, and recent research has also highlighted the importance of nematode life histories, environmental conditions, and other aspects of host ecology. Commonly, factors influencing parasitic nematode population genetics have been studied in isolation, such that an integrated view of the drivers of population genetic structure of parasitic nematodes is still lacking. Here, we seek to provide a comprehensive, broad, and integrative picture of these factors in parasitic nematodes of wild animals that will be a useful resource for investigators studying non-model parasitic nematodes in natural ecosystems. Increasingly, new methods of analysing the population genetics of nematodes are becoming available, and we consider the opportunities that these afford in resolving hitherto inaccessible questions of the population genetics of these important animals.

The first next-generation sequencing approach to the mitochondrial phylogeny of African monogenean parasites (Platyhelminthes: Gyrodactylidae and Dactylogyridae)

BACKGROUND

Monogenean flatworms are the main ectoparasites of fishes. Representatives of the species-rich families Gyrodactylidae and Dactylogyridae, especially those infecting cichlid fishes and clariid catfishes, are important parasites in African aquaculture, even more so due to the massive anthropogenic translocation of their hosts worldwide. Several questions on their evolution, such as the phylogenetic position of Macrogyrodactylus and the highly speciose Gyrodactylus, remain unresolved with available molecular markers. Also, diagnostics and population-level research would benefit from the development of higher-resolution genetic markers. We aim to offer genetic resources for work on African monogeneans by providing mitogenomic data of four species (two belonging to Gyrodactylidae, two to Dactylogyridae), and analysing their gene sequences and gene order from a phylogenetic perspective.

RESULTS

Using Illumina technology, the first four mitochondrial genomes of African monogeneans were assembled and annotated for the cichlid parasites Gyrodactylus nyanzae, Cichlidogyrus halli, Cichlidogyrus mbirizei (near-complete mitogenome) and the catfish parasite Macrogyrodactylus karibae (near-complete mitogenome). Complete nuclear ribosomal operons were also retrieved, as molecular vouchers. The start codon TTG is new for Gyrodactylus and for Dactylogyridae, as is the incomplete stop codon TA for Dactylogyridae. Especially the nad2 gene is promising for primer development. Gene order was identical for protein-coding genes and differed between the African representatives of these families only in a tRNA gene transposition. A mitochondrial phylogeny based on an alignment of nearly 12,500 bp including 12 protein-coding and two ribosomal RNA genes confirms that the Neotropical oviparous Aglaiogyrodactylus forficulatus takes a sister group position with respect to the other gyrodactylids, instead of the supposedly 'primitive' African Macrogyrodactylus. Inclusion of the African Gyrodactylus nyanzae confirms the paraphyly of Gyrodactylus. The position of the African dactylogyrid Cichlidogyrus is unresolved, although gene order suggests it is closely related to marine ancyrocephalines.

CONCLUSIONS

The amount of mitogenomic data available for gyrodactylids and dactylogyrids is increased by roughly one-third. Our study underscores the potential of mitochondrial genes and gene order in flatworm phylogenetics, and of next-generation sequencing for marker development for these non-model helminths for which few primers are available.

Low allelic diversity in vaccine candidates genes from different locations sustain hope for Fasciola hepatica immunization

Fasciola hepatica is a trematode parasite that causes fasciolosis in animals and humans. Fasciolosis is usually treated with triclabendazole, although drug-resistant parasites have been described in several geographical locations. An alternative to drug treatment would be the use of a vaccine, although vaccination studies that have been performed mainly in ruminants over the last 30 years, show high variability in the achieved protection and are not yet ready for commercialisation. Since F. hepatica exhibits a high degree of genomic polymorphism, variation in vaccine efficacy could be attributed, at least partially, to phenotypic differences in vaccine candidate sequences amongst parasites used in the challenge infections. To begin to address this issue, a collection of F. hepatica isolates from geographically dispersed regions, as well as parasites obtained from vaccination trials performed against a field isolate from Uruguay and the experimentally maintained South Gloucester isolate (Ridgeway Research, UK), were compiled to establish a F. hepatica Biobank. These collected isolates were used for the genetic analysis of several vaccine candidates that are important in host-parasite interactions and are the focus of the H2020 PARAGONE vaccine project (https://www.paragoneh2020.eu/), namely FhCL1, FhCL2, FhPrx, FhLAP and FhHDM. Our results show that F. hepatica exhibits a high level of conservation in the sequences encoding each of these proteins. The consequential low variability in these vaccine candidates amongst parasites from different geographical regions reinforces the idea that they would be suitable immunogens against liver fluke isolates worldwide.

Discovery of genomic intervals that underlie nematode responses to benzimidazoles

Parasitic nematodes impose a debilitating health and economic burden across much of the world. Nematode resistance to anthelmintic drugs threatens parasite control efforts in both human and veterinary medicine. Despite this threat, the genetic landscape of potential resistance mechanisms to these critical drugs remains largely unexplored. Here, we exploit natural variation in the model nematodes Caenorhabditis elegans and Caenorhabditis briggsae to discover quantitative trait loci (QTL) that control sensitivity to benzimidazoles widely used in human and animal medicine. High-throughput phenotyping of albendazole, fenbendazole, mebendazole, and thiabendazole responses in panels of recombinant lines led to the discovery of over 15 QTL in C. elegans and four QTL in C. briggsae associated with divergent responses to these anthelmintics. Many of these QTL are conserved across benzimidazole derivatives, but others show drug and dose specificity. We used near-isogenic lines to recapitulate and narrow the C. elegans albendazole QTL of largest effect and identified candidate variants correlated with the resistance phenotype. These QTL do not overlap with known benzimidazole target resistance genes from parasitic nematodes and present specific new leads for the discovery of novel mechanisms of nematode benzimidazole resistance. Analyses of orthologous genes reveal conservation of candidate benzimidazole resistance genes in medically important parasitic nematodes. These data provide a basis for extending these approaches to other anthelmintic drug classes and a pathway towards validating new markers for anthelmintic resistance that can be deployed to improve parasite disease control.

The treatment of roundworm (nematode) infections in both humans and animals relies on a small number of anti-parasitic drugs. Resistance to these drugs has appeared in veterinary parasite populations and is a growing concern in human medicine. A better understanding of the genetic basis for parasite drug resistance can be used to help maintain the effectiveness of anti-parasitic drugs and to slow or to prevent the spread of drug resistance in parasite populations. This goal is hampered by the experimental intractability of nematode parasites. Here, we use non-parasitic model nematodes to systematically explore responses to the critical benzimidazole class of anti-parasitic compounds. Using a quantitative genetics approach, we discovered unique genomic intervals that control drug effects, and we identified differences in the effects of these intervals across compounds and doses. We were able to narrow a major-effect genomic region associated with albendazole resistance and to establish that candidate genes discovered in our genetic mappings are largely conserved in important human and animal parasites. This work provides new leads for understanding parasite drug resistance and contributes a powerful template that can be extended to other anti-parasitic drug classes.

De novo assembly of the complex genome of Nippostrongylus brasiliensis using MinION long reads

BACKGROUND

Eukaryotic genome assembly remains a challenge in part due to the prevalence of complex DNA repeats. This is a particularly acute problem for holocentric nematodes because of the large number of satellite DNA sequences found throughout their genomes. These have been recalcitrant to most genome sequencing methods. At the same time, many nematodes are parasites and some represent a serious threat to human health. There is a pressing need for better molecular characterization of animal and plant parasitic nematodes. The advent of long-read DNA sequencing methods offers the promise of resolving complex genomes.

RESULTS

Using Nippostrongylus brasiliensis as a test case, applying improved base-calling algorithms and assembly methods, we demonstrate the feasibility of de novo genome assembly matching current community standards using only MinION long reads. In doing so, we uncovered an unexpected diversity of very long and complex DNA sequences repeated throughout the N. brasiliensis genome, including massive tandem repeats of tRNA genes.

CONCLUSION

Base-calling and assembly methods have improved sufficiently that de novo genome assembly of large complex genomes is possible using only long reads. The method has the added advantage of preserving haplotypic variants and so has the potential to be used in population analyses.

Human Interventions: Driving Forces of Mosquito Evolution

One of the most common strategies for controlling mosquito-borne diseases relies on the use of chemical pesticides to repel or kill the mosquito vector. Pesticide exposure interferes with several key biological functions in the mosquito and triggers a variety of adaptive responses whose underlying mechanisms are only partially elucidated. The availability of reference genome sequences opens up the possibility of tracking signatures of evolutionary changes, including the most recent, across the genomes of many vector species. In this review, we highlight the recent genomic changes, which contribute to the fascinating adaptation of malaria vectors of the sub-Saharan African region to intensive insecticide-based interventions. We emphasize the operational significance of detailed genomic knowledge for current monitoring and decision-making.