Distribution and evolution of repeated sequences in genomes of Triatominae (Hemiptera-Reduviidae) inferred from genomic in situ hybridization

The subfamily Triatominae, vectors of Chagas disease, comprises 140 species characterized by a highly homogeneous chromosome number. We analyzed the chromosomal distribution and evolution of repeated sequences in Triatominae genomes by Genomic in situ Hybridization using Triatoma delpontei and Triatoma infestans genomic DNAs as probes. Hybridizations were performed on their own chromosomes and on nine species included in six genera from the two main tribes: Triatomini and Rhodniini. Genomic probes clearly generate two different hybridization patterns, dispersed or accumulated in specific regions or chromosomes. The three used probes generate the same hybridization pattern in each species. However, these patterns are species-specific. In closely related species, the probes strongly hybridized in the autosomal heterochromatic regions, resembling C-banding and DAPI patterns. However, in more distant species these co-localizations are not observed. The heterochromatic Y chromosome is constituted by highly repeated sequences, which is conserved among 10 species of Triatomini tribe suggesting be an ancestral character for this group. However, the Y chromosome in Rhodniini tribe is markedly different, supporting the early evolutionary dichotomy between both tribes. In some species, sex chromosomes and autosomes shared repeated sequences, suggesting meiotic chromatin exchanges among these heterologous chromosomes. Our GISH analyses enabled us to acquire not only reliable information about autosomal repeated sequences distribution but also an insight into sex chromosome evolution in Triatominae. Furthermore, the differentiation obtained by GISH might be a valuable marker to establish phylogenetic relationships and to test the controversial origin of the Triatominae subfamily.

The evolution of sex ratio distorter suppression affects a 25 cM genomic region in the butterfly Hypolimnas bolina

Symbionts that distort their host's sex ratio by favouring the production and survival of females are common in arthropods. Their presence produces intense Fisherian selection to return the sex ratio to parity, typified by the rapid spread of host ‘suppressor’ loci that restore male survival/development. In this study, we investigated the genomic impact of a selective event of this kind in the butterfly Hypolimnas bolina. Through linkage mapping, we first identified a genomic region that was necessary for males to survive Wolbachia-induced male-killing. We then investigated the genomic impact of the rapid spread of suppression, which converted the Samoan population of this butterfly from a 100∶1 female-biased sex ratio in 2001 to a 1∶1 sex ratio by 2006. Models of this process revealed the potential for a chromosome-wide effect. To measure the impact of this episode of selection directly, the pattern of genetic variation before and after the spread of suppression was compared. Changes in allele frequencies were observed over a 25 cM region surrounding the suppressor locus, with a reduction in overall diversity observed at loci that co-segregate with the suppressor. These changes exceeded those expected from drift and occurred alongside the generation of linkage disequilibrium. The presence of novel allelic variants in 2006 suggests that the suppressor was likely to have been introduced via immigration rather than through de novo mutation. In addition, further sampling in 2010 indicated that many of the introduced variants were lost or had declined in frequency since 2006. We hypothesize that this loss may have resulted from a period of purifying selection, removing deleterious material that introgressed during the initial sweep. Our observations of the impact of suppression of sex ratio distorting activity reveal a very wide genomic imprint, reflecting its status as one of the strongest selective forces in nature.

The sex ratio of the offspring produced by an individual can be an evolutionary battleground. In many arthropod species, maternally inherited microbes selectively kill male hosts, and the host may in turn evolve strategies to restore the production or survival of males. When males are rare, the intensity of selection on the host may be extreme. We recently observed one such episode, in which the population sex ratio of the butterfly Hypolimnas bolina shifted from 100 females per male to near parity, through the evolution of a suppressor gene. In our current study, we investigate the hypothesis that the strength of selection in this case was so strong that the genomic impact would go well beyond the suppressor gene itself. After mapping the location of the suppressor within the genome of H. bolina, we examined changes in genetic variation at sites on the same chromosome as the suppressor. We show that a broad region of the genome was affected by the spread of the suppressor. Our data also suggest that the selection may have been sufficiently strong to introduce deleterious material into the population, which was later purged by selection.

Transcription activator-like effector nuclease (TALEN)-mediated female-specific sterility in the silkworm, Bombyx mori

Engineering sex-specific sterility is critical for developing transgene-based sterile insect technology. Targeted genome engineering achieved by customized zinc-finger nuclease, transcription activator-like effector nuclease (TALEN) or clustered, regularly interspaced, short palindromic repeats/Cas9 systems has been exploited extensively in a variety of model organisms; however, screening mutated individuals without a detectable phenotype is still challenging. In addition, genetically recessive mutations only detectable in homozygotes make the experiments time-consuming. In the present study, we model a novel genetic system in the silkworm, Bombyx mori, that results in female-specific sterility by combining transgenesis with TALEN technologies. This system induces sex-specific sterility at a high efficiency by targeting the female-specific exon of the B. mori doublesex (Bmdsx) gene, which has sex-specific splicing isoforms regulating somatic sexual development. Transgenic animals co-expressing TALEN left and right arms targeting the female-specific Bmdsx exon resulted in somatic mutations and female mutants lost fecundity because of lack of egg storage and abnormal external genitalia. The wild-type sexual dimorphism of abdominal segment was not evident in mutant females. In contrast, there were no deleterious effects in mutant male moths. The current somatic TALEN technologies provide a promising approach for future insect functional genetics, thus providing the basis for the development of attractive genetic alternatives for insect population management.

Transcriptional changes induced by in vivo exposure to pentachlorophenol (PCP) in Chironomus riparius (Diptera) aquatic larvae

Pentachlorophenol (PCP) has been extensively used worldwide as a pesticide and biocide and is frequently detected in the aquatic environment. In the present work, the toxicity of PCP was investigated in Chironomus riparius aquatic larvae. The effects following short- and long-term exposures were evaluated at the molecular level by analyzing changes in the transcriptional profile of different endocrine genes, as well as in genes involved in the stress response and detoxification. Interestingly, although no differences were found after 12- and 24-h treatments, at 96-h exposures PCP was able to induce significant increases in transcripts from the ecdysone receptor gene (EcR), the early ecdysone-inducible E74 gene, the estrogen-related receptor gene (ERR), the Hsp70 gene and the CYP4G gene. In contrast, the Hsp27 gene appeared to be downregulated, while the ultraspiracle gene (usp) (insect ortholog of the retinoid X receptor) was not altered in any of the conditions assayed. Moreover, Glutathione-S-Transferase (GST) activity was not affected. The results obtained show the ability of PCP to modulate transcription of different biomarker genes from important cellular metabolic activities, which could be useful in genomic approaches to monitoring. In particular, the significant upregulation of hormonal genes represents the first evidence at the genomic level of the potential endocrine disruptive effects of PCP on aquatic invertebrates.

Accidental genetic engineers: horizontal sequence transfer from parasitoid wasps to their Lepidopteran hosts

We show here that 105 regions in two Lepidoptera genomes appear to derive from horizontally transferred wasp DNA. We experimentally verified the presence of two of these sequences in a diverse set of silkworm (Bombyx mori) genomes. We hypothesize that these horizontal transfers are made possible by the unusual strategy many parasitoid wasps employ of injecting hosts with endosymbiotic polydnaviruses to minimize the host's defense response. Because these virus-like particles deliver wasp DNA to the cells of the host, there has been much interest in whether genetic information can be permanently transferred from the wasp to the host. Two transferred sequences code for a BEN domain, known to be associated with polydnaviruses and transcriptional regulation. These findings represent the first documented cases of horizontal transfer of genes between two organisms by a polydnavirus. This presents an interesting evolutionary paradigm in which host species can acquire new sequences from parasitoid wasps that attack them. Hymenoptera and Lepidoptera diverged ∼300 MYA, making this type of event a source of novel sequences for recipient species. Unlike many other cases of horizontal transfer between two eukaryote species, these sequence transfers can be explained without the need to invoke the sequences 'hitchhiking' on a third organism (e.g. retrovirus) capable of independent reproduction. The cellular machinery necessary for the transfer is contained entirely in the wasp genome. The work presented here is the first such discovery of what is likely to be a broader phenomenon among species affected by these wasps.

Polymorphic chromosomal inversions in Anopheles moucheti, a major malaria vector in Central Africa

Anopheles moucheti Evans (Diptera: Culicidae) is a major vector of malaria in forested areas of Central Africa. However, few genetic tools are available for this species. The present study represents the first attempt to characterize chromosomes in An. moucheti females collected in Cameroon. Ovarian nurse cells contained polytene chromosomes, which were suitable for standard cytogenetic applications. The presence of three polymorphic chromosomal inversions in An. moucheti was revealed. Two of these inversions were located on the 2R chromosome arm. The homology between the 2R chromosome arms of An. moucheti and Anopheles gambiae Giles was established by fluorescent in situ hybridization of six An. gambiae genic sequences. Mapping of the probes on chromosomes of An. moucheti detected substantial gene order reshuffling between the two species. The presence of polytene chromosomes and polymorphic inversions in An. moucheti provides a new basis for further population genetic, taxonomic and ecological studies of this neglected malaria vector.

An evaluation of the ecological relationship between Drosophila species and their parasitoid wasps as an opportunity for horizontal transposon transfer

Evidences of horizontal transfer, the exchange of genetic material between reproductively isolated species, have accumulated over the last decades, including for multicellular eukaryotic organisms. However, the mechanisms and ecological relationships that promote such phenomenon is still poorly known. Host-parasite interaction is one type of relationship usually pointed in the literature that could potentially increase the probability of the horizontal transfer between species, because the species involved in such relationships are generally in close contact. Transposable elements, which are well-known genomic parasites, are DNA entities that tend to be involved in horizontal transfer due to their ability to mobilize between different genomic locations. Using Drosophila species and their parasitoid wasps as a host-parasite model, we evaluated the hypothesis that horizontal transposon transfers (HTTs) are more frequent in this set of species than in species that do not exhibit a close ecological and phylogenetic relationship. For this purpose, we sequenced two sets of species using a metagenomic and single-species genomic sampling approach through next-generation DNA sequencing. The first set was composed of five generalist Drosophila (D. maculifrons, D. bandeirantorum, D. polymorpha, D. mercatorum and D. willistoni) species and their associated parasitoid wasps, whereas the second set was composed of D. incompta, which is a flower specialist species, and its parasitoid wasp. We did not find strong evidence of HTT in the two sets of Drosophila and wasp parasites. However, at least five cases of HTT were observed between the generalist and specialist Drosophila species. Moreover, we detected an HT event involving a Wolbachia lineage between generalist and specialist species, indicating that these endosymbiotic bacteria could play a role as HTT vectors. In summary, our results do not support the hypothesis of prevalent HTT between species with a host-parasite relationship, at least for the studied wasp-Drosophila pairs. Moreover, it suggests that other mechanisms or parasites are involved in promoting HTT between Drosophila species as the Wolbachia endosymbiotic bacteria.

Heritable genome editing with CRISPR/Cas9 in the silkworm, Bombyx mori

We report the establishment of an efficient and heritable gene mutagenesis method in the silkworm Bombyx mori using modified type II clustered regularly interspaced short palindromic repeats (CRISPR) with an associated protein (Cas9) system. Using four loci Bm-ok, BmKMO, BmTH, and Bmtan as candidates, we proved that genome alterations at specific sites could be induced by direct microinjection of specific guide RNA and Cas9-mRNA into silkworm embryos. Mutation frequencies of 16.7-35.0% were observed in the injected generation, and DNA fragments deletions were also noted. Bm-ok mosaic mutants were used to test for mutant heritability due to the easily determined translucent epidermal phenotype of Bm-ok-disrupted cells. Two crossing strategies were used. In the first, injected Bm-ok moths were crossed with wild-type moths, and a 28.6% frequency of germline mutation transmission was observed. In the second strategy, two Bm-ok mosaic mutant moths were crossed with each other, and 93.6% of the offsprings appeared mutations in both alleles of Bm-ok gene (compound heterozygous). In summary, the CRISPR/Cas9 system can act as a highly specific and heritable gene-editing tool in Bombyx mori.

The expression and phylogenetics of the Inhibitor Cysteine Knot peptide OCLP1 in the honey bee Apis mellifera

Small cysteine-rich peptides have diverse functions in insects including antimicrobial defense, phenoloxidase activity regulation, and toxic inhibition of ion channels of prey or predator. We combined bioinformatics and measurements of transcript abundance to start characterizing AmOCLP1, a recently discovered Inhibitor Cysteine Knot peptide in the honey bee Apis mellifera. We found that the genomes of ants, bees, and the wasp Nasonia vitripennis encode orthologous sequences indicating that OCLP1 is a conserved peptide and not unique to the honey bee. Search of available EST libraries and quantitative real time PCR analyses indicate that the transcript of AmOCLP1 is ubiquitous with expression in life stages ranging from embryos to adults and in all tested tissues. In worker honey bees AmOCLP1 expression was not associated with age or task and did not show clear enrichment in any of the tested tissues. There was however a consistent trend toward higher transcript levels in the abdomen of foragers relative to levels in the head or thorax, and compared to levels in the abdomen of younger worker bees. By contrast, in drones AmOCLP1 transcript levels appeared higher in the head relative to the abdomen. Finer analyses of the head and abdomen indicated that the AmOCLP1 transcript is not enriched in the stinger and the associated venom sac or in cephalic exocrine glands. The evolutionary conservation in the Hymenoptera, the ubiquitous expression, and the lack of enrichment in the venom gland, stinger, exocrine glands, and the brain are not consistent with the hypotheses that OCLP1 is a secreted honeybee toxin or an endotoxin acting in the central nervous system. Rather we hypothesize that OCLP1 is a conserved antimicrobial or phenoloxidase inhibitor peptide.

Effect of Incomplete Lineage Sorting On Tree-Reconciliation-Based Inference of Gene Duplication

In the tree reconciliation approach to infer the duplication history of a gene family, the gene (family) tree is compared to the corresponding species tree. Incomplete lineage sorting (ILS) gives rise to stochastic variation in the topology of a gene tree and hence likely introduces false duplication events when a tree reconciliation method is used. We quantify the effect of ILS on gene duplication inference in a species tree in terms of the expected number of false duplication events inferred from reconciling a random gene tree, which occurs with a probability predicted in coalescent theory, and the species tree. We computationally examine the relationship between the effect of ILS on duplication inference in a species tree and its topological parameters. Our findings suggest that ILS may cause non-negligible bias on duplication inference, particularly on an asymmetric species tree. Hence, when gene duplication is inferred via tree reconciliation or any other approach that takes gene tree topology into account, the ILS-induced bias should be examined cautiously.

Genome-guided transcript assembly by integrative analysis of RNA sequence data

The identification of full length transcripts entirely from short-read RNA sequencing data (RNA-seq) remains a challenge in the annotation of genomes. Here we describe an automated pipeline for genome annotation that integrates RNA-seq and gene-boundary data sets, which we call Generalized RNA Integration Tool, or GRIT. Applying GRIT to Drosophila melanogaster short-read RNA-seq, cap analysis of gene expression (CAGE) and poly(A)-site-seq data collected for the modENCODE project, we recovered the vast majority of previously annotated transcripts and doubled the total number of transcripts cataloged. We found that 20% of protein coding genes encode multiple protein-localization signals and that, in 20-d-old adult fly heads, genes with multiple polyadenylation sites are more common than genes with alternative splicing or alternative promoters. GRIT demonstrates 30% higher precision and recall than the most widely used transcript assembly tools. GRIT will facilitate the automated generation of high-quality genome annotations without the need for extensive manual annotation.

Drosophila epigenome reorganization during oocyte differentiation and early embryogenesis

In sexually reproducing organisms, propagation of the species relies on specialized haploid cells (gametes) produced by germ cells. During their development in the adult germline, the female and male gametes undergo a complex differentiation process that requires transcriptional regulation and chromatin reorganization. After fertilization, the gametes then go through extensive epigenetic reprogramming, which resets the cells to a totipotent state essential for the development of the embryo. Several histone modifications characterize distinct developmental stages of gamete formation and early embryonic development, but it is unknown whether these modifications have any physiological role. Furthermore, accumulating evidence suggests that environmentally induced chromatin changes can be inherited, yet the mechanisms underlying zygotic inheritance of the gamete epigenome remain unclear. This review gives a brief overview of the mechanisms of transgenerational epigenetic inheritance and examines the function of epigenetics during oogenesis and early embryogenesis with a focus on histone posttranslational modifications.

A deterministic analysis of genome integrity during neoplastic growth in Drosophila

The development of cancer has been associated with the gradual acquisition of genetic alterations leading to a progressive increase in malignancy. In various cancer types this process is enabled and accelerated by genome instability. While genome sequencing-based analysis of tumor genomes becomes increasingly a standard procedure in human cancer research, the potential necessity of genome instability for tumorigenesis in Drosophila melanogaster has, to our knowledge, never been determined at DNA sequence level. Therefore, we induced formation of tumors by depletion of the Drosophila tumor suppressor Polyhomeotic and subjected them to genome sequencing. To achieve a highly resolved delineation of the genome structure we developed the Deterministic Structural Variation Detection (DSVD) algorithm, which identifies structural variations (SVs) with high accuracy and at single base resolution. The employment of long overlapping paired-end reads enables DSVD to perform a deterministic, i.e. fragment size distribution independent, identification of a large size spectrum of SVs. Application of DSVD and other algorithms to our sequencing data reveals substantial genetic variation with respect to the reference genome reflecting temporal separation of the reference and laboratory strains. The majority of SVs, constituted by small insertions/deletions, is potentially caused by erroneous replication or transposition of mobile elements. Nevertheless, the tumor did not depict a loss of genome integrity compared to the control. Altogether, our results demonstrate that genome stability is not affected inevitably during sustained tumor growth in Drosophila implying that tumorigenesis, in this model organism, can occur irrespective of genome instability and the accumulation of specific genetic alterations.

Comparison of ultra-conserved elements in drosophilids and vertebrates

Metazoan genomes contain many ultra-conserved elements (UCEs), long sequences identical between distant species. In this study we identified UCEs in drosophilid and vertebrate species with a similar level of phylogenetic divergence measured at protein-coding regions, and demonstrated that both the length and number of UCEs are larger in vertebrates. The proportion of non-exonic UCEs declines in distant drosophilids whilst an opposite trend was observed in vertebrates. We generated a set of 2,126 Sophophora UCEs by merging elements identified in several drosophila species and compared these to the eutherian UCEs identified in placental mammals. In contrast to vertebrates, the Sophophora UCEs are depleted around transcription start sites. Analysis of 52,954 P-element, piggyBac and Minos insertions in the D. melanogaster genome revealed depletion of the P-element and piggyBac insertions in and around the Sophophora UCEs. We examined eleven fly strains with transposon insertions into the intergenic UCEs and identified associated phenotypes in five strains. Four insertions behave as recessive lethals, and in one case we observed a suppression of the marker gene within the transgene, presumably by silenced chromatin around the integration site. To confirm the lethality is caused by integration of transposons we performed a phenotype rescue experiment for two stocks and demonstrated that the excision of the transposons from the intergenic UCEs restores viability. Sequencing of DNA after the transposon excision in one fly strain with the restored viability revealed a 47 bp insertion at the original transposon integration site suggesting that the nature of the mutation is important for the appearance of the phenotype. Our results suggest that the UCEs in flies and vertebrates have both common and distinct features, and demonstrate that a significant proportion of intergenic drosophila UCEs are sensitive to disruption.

Decoding the ubiquitin-mediated pathway of arthropod disease vectors

Protein regulation by ubiquitin has been extensively described in model organisms. However, characterization of the ubiquitin machinery in disease vectors remains mostly unknown. This fundamental gap in knowledge presents a concern because new therapeutics are needed to control vector-borne diseases, and targeting the ubiquitin machinery as a means for disease intervention has been already adopted in the clinic. In this study, we employed a bioinformatics approach to uncover the ubiquitin-mediated pathway in the genomes of Anopheles gambiae, Aedes aegypti, Culex quinquefasciatus, Ixodes scapularis, Pediculus humanus and Rhodnius prolixus. We observed that (1) disease vectors encode a lower percentage of ubiquitin-related genes when compared to Drosophila melanogaster, Mus musculus and Homo sapiens but not Saccharomyces cerevisiae; (2) overall, there are more proteins categorized as E3 ubiquitin ligases when compared to E2-conjugating or E1-activating enzymes; (3) the ubiquitin machinery within the three mosquito genomes is highly similar; (4) ubiquitin genes are more than doubled in the Chagas disease vector (R. prolixus) when compared to other arthropod vectors; (5) the deer tick I. scapularis and the body louse (P. humanus) genomes carry low numbers of E1-activating enzymes and HECT-type E3 ubiquitin ligases; (6) R. prolixus have low numbers of RING-type E3 ubiquitin ligases; and (7) C. quinquefasciatus present elevated numbers of predicted F-box E3 ubiquitin ligases, JAB and UCH deubiquitinases. Taken together, these findings provide novel opportunities to study the interaction between a pathogen and an arthropod vector.

A SNP based high-density linkage map of Apis cerana reveals a high recombination rate similar to Apis mellifera

Background

The Eastern honey bee, Apis cerana Fabricius, is distributed in southern and eastern Asia, from India and China to Korea and Japan and southeast to the Moluccas. This species is also widely kept for honey production besides Apis mellifera. Apis cerana is also a model organism for studying social behavior, caste determination, mating biology, sexual selection, and host-parasite interactions. Few resources are available for molecular research in this species, and a linkage map was never constructed. A linkage map is a prerequisite for quantitative trait loci mapping and for analyzing genome structure. We used the Chinese honey bee, Apis cerana cerana to construct the first linkage map in the Eastern honey bee.

Results

F2 workers (N = 103) were genotyped for 126,990 single nucleotide polymorphisms (SNPs). After filtering low quality and those not passing the Mendel test, we obtained 3,000 SNPs, 1,535 of these were informative and used to construct a linkage map. The preliminary map contains 19 linkage groups, we then mapped the 19 linkage groups to 16 chromosomes by comparing the markers to the genome of A. mellfiera. The final map contains 16 linkage groups with a total of 1,535 markers. The total genetic distance is 3,942.7 centimorgans (cM) with the largest linkage group (180 loci) measuring 574.5 cM. Average marker interval for all markers across the 16 linkage groups is 2.6 cM.

Conclusion

We constructed a high density linkage map for A. c. cerana with 1,535 markers. Because the map is based on SNP markers, it will enable easier and faster genotyping assays than randomly amplified polymorphic DNA or microsatellite based maps used in A. mellifera.

Optimized position weight matrices in prediction of novel putative binding sites for transcription factors in the Drosophila melanogaster genome

Position weight matrices (PWMs) have become a tool of choice for the identification of transcription factor binding sites in DNA sequences. DNA-binding proteins often show degeneracy in their binding requirement and thus the overall binding specificity of many proteins is unknown and remains an active area of research. Although existing PWMs are more reliable predictors than consensus string matching, they generally result in a high number of false positive hits. Our previous study introduced a promising approach to PWM refinement in which known motifs are used to computationally mine putative binding sites directly from aligned promoter regions using composition of similar sites. In the present study, we extended this technique originally tested on single examples of transcription factors (TFs) and showed its capability to optimize PWM performance to predict new binding sites in the fruit fly genome. We propose refined PWMs in mono- and dinucleotide versions similarly computed for a large variety of transcription factors of Drosophila melanogaster. Along with the addition of many auxiliary sites the optimization includes variation of the PWM motif length, the binding sites location on the promoters and the PWM score threshold. To assess the predictive performance of the refined PWMs we compared them to conventional TRANSFAC and JASPAR sources. The results have been verified using performed tests and literature review. Overall, the refined PWMs containing putative sites derived from real promoter content processed using optimized parameters had better general accuracy than conventional PWMs.

Genomic treasure troves: complete genome sequencing of herbarium and insect museum specimens

Unlocking the vast genomic diversity stored in natural history collections would create unprecedented opportunities for genome-scale evolutionary, phylogenetic, domestication and population genomic studies. Many researchers have been discouraged from using historical specimens in molecular studies because of both generally limited success of DNA extraction and the challenges associated with PCR-amplifying highly degraded DNA. In today's next-generation sequencing (NGS) world, opportunities and prospects for historical DNA have changed dramatically, as most NGS methods are actually designed for taking short fragmented DNA molecules as templates. Here we show that using a standard multiplex and paired-end Illumina sequencing approach, genome-scale sequence data can be generated reliably from dry-preserved plant, fungal and insect specimens collected up to 115 years ago, and with minimal destructive sampling. Using a reference-based assembly approach, we were able to produce the entire nuclear genome of a 43-year-old Arabidopsis thaliana (Brassicaceae) herbarium specimen with high and uniform sequence coverage. Nuclear genome sequences of three fungal specimens of 22–82 years of age (Agaricus bisporus, Laccaria bicolor, Pleurotus ostreatus) were generated with 81.4–97.9% exome coverage. Complete organellar genome sequences were assembled for all specimens. Using de novo assembly we retrieved between 16.2–71.0% of coding sequence regions, and hence remain somewhat cautious about prospects for de novo genome assembly from historical specimens. Non-target sequence contaminations were observed in 2 of our insect museum specimens. We anticipate that future museum genomics projects will perhaps not generate entire genome sequences in all cases (our specimens contained relatively small and low-complexity genomes), but at least generating vital comparative genomic data for testing (phylo)genetic, demographic and genetic hypotheses, that become increasingly more horizontal. Furthermore, NGS of historical DNA enables recovering crucial genetic information from old type specimens that to date have remained mostly unutilized and, thus, opens up a new frontier for taxonomic research as well.

Widespread selection across coding and noncoding DNA in the pea aphid genome

Genome-wide patterns of diversity and selection are critical measures for understanding how evolution has shaped the genome. Yet, these population genomic estimates are available for only a limited number of model organisms. Here we focus on the population genomics of the pea aphid (Acyrthosiphon pisum). The pea aphid is an emerging model system that exhibits a range of intriguing biological traits not present in classic model systems. We performed low-coverage genome resequencing of 21 clonal pea aphid lines collected from alfalfa host plants in North America to characterize genome-wide patterns of diversity and selection. We observed an excess of low-frequency polymorphisms throughout coding and noncoding DNA, which we suggest is the result of a founding event and subsequent population expansion in North America. Most gene regions showed lower levels of Tajima's D than synonymous sites, suggesting that the majority of the genome is not evolving neutrally but rather exhibits significant constraint. Furthermore, we used the pea aphid's unique manner of X-chromosome inheritance to assign genomic scaffolds to either autosomes or the X chromosome. Comparing autosomal vs. X-linked sequence variation, we discovered that autosomal genes show an excess of low frequency variants indicating that purifying selection acts more efficiently on the X chromosome. Overall, our results provide a critical first step in characterizing the genetic diversity and evolutionary pressures on an aphid genome.

The pea aphid (Acyrthosiphon pisum) genome encodes two divergent early developmental programs

The pea aphid (Acyrthosiphon pisum) can reproduce either sexually or asexually (parthenogenetically), giving rise, in each case, to almost identical adults. These two modes of reproduction are accompanied by differences in ovarian morphology and the developmental environment of the offspring, with sexual forms producing eggs that are laid, whereas asexual development occurs within the mother. Here we examine the effect each mode of reproduction has on the expression of key maternal and axis patterning genes; orthodenticle (otd), hunchback (hb), caudal (cad) and nanos (nos). We show that three of these genes (Ap-hb, Ap-otd and Ap-cad) are expressed differently between the sexually and asexually produced oocytes and embryos of the pea aphid. We also show, using immunohistochemistry and cytoskeletal inhibitors, that Ap-hb RNA is localized differently between sexually and asexually produced oocytes, and that this is likely due to differences in the 3' untranslated regions of the RNA. Furthermore, Ap-hb and Ap-otd have extensive expression domains in early sexually produced embryos, but are not expressed at equivalent stages in asexually produced embryos. These differences in expression likely correspond with substantial changes in the gene regulatory networks controlling early development in the pea aphid. These data imply that in the evolution of parthenogenesis a new program has evolved to control the development of asexually produced embryos, whilst retaining the existing, sexual, developmental program. The patterns of modification of these developmental processes mirror the changes that we see in developmental processes between species, in that early acting pathways in development are less constrained, and evolve faster, than later ones. We suggest that the evolution of the novel asexual development pathway in aphids is not a simple modification of an ancestral system, but the evolution of two very different developmental mechanisms occurring within a single species.

Complex genetic architecture of cardiac disease in a wild type inbred strain of Drosophila melanogaster

Natural populations of the fruit fly, Drosophila melanogaster, segregate genetic variation that leads to cardiac disease phenotypes. One nearly isogenic line from a North Carolina peach orchard, WE70, is shown to harbor two genetically distinct heart phenotypes: elevated incidence of arrhythmias, and a dramatically constricted heart diameter in both diastole and systole, with resemblance to restrictive cardiomyopathy in humans. Assuming the source to be rare variants of large effect, we performed Bulked Segregant Analysis using genomic DNA hybridization to Affymetrix chips to detect single feature polymorphisms, but found that the mutant phenotypes are more likely to have a polygenic basis. Further mapping efforts revealed a complex architecture wherein the constricted cardiomyopathy phenotype was observed in individual whole chromosome substitution lines, implying that variants on both major autosomes are sufficient to produce the phenotype. A panel of 170 Recombinant Inbred Lines (RIL) was generated, and a small subset of mutant lines selected, but these each complemented both whole chromosome substitutions, implying a non-additive (epistatic) contribution to the “disease” phenotype. Low coverage whole genome sequencing was also used to attempt to map chromosomal regions contributing to both the cardiomyopathy and arrhythmia, but a polygenic architecture had to be again inferred to be most likely. These results show that an apparently simple rare phenotype can have a complex genetic basis that would be refractory to mapping by deep sequencing in pedigrees. We present this as a cautionary tale regarding assumptions related to attempts to map new disease mutations on the assumption that probands carry a single causal mutation.

Transposition-driven genomic heterogeneity in the Drosophila brain

Recent studies in mammals have documented the neural expression and mobility of retrotransposons and have suggested that neural genomes are diverse mosaics. We found that transposition occurs among memory-relevant neurons in the Drosophila brain. Cell type-specific gene expression profiling revealed that transposon expression is more abundant in mushroom body (MB) αβ neurons than in neighboring MB neurons. The Piwi-interacting RNA (piRNA) proteins Aubergine and Argonaute 3, known to suppress transposons in the fly germline, are expressed in the brain and appear less abundant in αβ MB neurons. Loss of piRNA proteins correlates with elevated transposon expression in the brain. Paired-end deep sequencing identified more than 200 de novo transposon insertions in αβ neurons, including insertions into memory-relevant loci. Our observations indicate that genomic heterogeneity is a conserved feature of the brain.

Comparative metatranscriptomic signatures of wood and paper feeding in the gut of the termite Reticulitermes flavipes (Isoptera: Rhinotermitidae)

Termites are highly eusocial insects that thrive on recalcitrant materials like wood and soil and thus play important roles in global carbon recycling and also in damaging wooden structures. Termites, such as Reticulitermes flavipes (Rhinotermitidae), owe their success to their ability to extract nutrients from lignocellulose (a major component of wood) with the help of gut-dwelling symbionts. With the aim to gain new insights into this enzymatic process we provided R. flavipes with a complex lignocellulose (wood) or pure cellulose (paper) diet and followed the resulting differential gene expression on a custom oligonucleotide-microarray platform. We identified a set of expressed sequence tags (ESTs) with differential abundance between the two diet treatments and demonstrated the source (host/symbiont) of these genes, providing novel information on termite nutritional symbiosis. Our results reveal: (1) the majority of responsive wood- and paper-abundant ESTs are from host and symbionts, respectively; (2) distinct pathways are associated with lignocellulose and cellulose feeding in both host and symbionts; and (3) sets of diet-responsive ESTs encode putative digestive and wood-related detoxification enzymes. Thus, this study illuminates the dynamics of termite nutritional symbiosis and reveals a pool of genes as potential targets for termite control and functional studies of termite-symbiont interactions.

The complete mitochondrial genome of spittlebug Paphnutius ruficeps (Insecta: Hemiptera: Cercopidae) with a fairly short putative control region

The mitochondrial genome of the spittlebug Paphnutius ruficeps is a double-strand DNA circular molecule of 14,841 bp with a total A and T content of 73.8%. It is one of the shortest genomes among published hemipteran mitogenomes and encodes 13 protein-coding genes, 2 ribosome RNA genes and 22 transfer RNA (tRNA) genes. The gene order is consistent with the hypothesized ancestral arthropod genome arrangement. Most of the protein-coding genes use ATG as start and TAA as stop codon. The codons show an evident bias toward the nucleotides T and A at the third codon position and the most commonly used codons contain more A and T than their synonymous ones. The anticodons of the 22 tRNA genes are identical to those of the mitogenome of Philaenus spumarius, another studied spittlebug. All the tRNAs could be folded into traditional clover leaf secondary structures. The putative control region (traditionally called A + T-rich region) is the main non-coding part of the mitogenome. The AT content of this region (74.5%) is not significantly higher than that of the total mitogenome (73.8%) and slightly lower than that of the N-chain protein-coding genes (75.3%). The absence of repeat sequences as well as its short length is the most obvious characteristics of the mitochondrial genome of Paphnutius ruficeps compared with those of other published hemipteran species.

Draft genome of the mountain pine beetle, Dendroctonus ponderosae Hopkins, a major forest pest

BACKGROUND

The mountain pine beetle, Dendroctonus ponderosae Hopkins, is the most serious insect pest of western North American pine forests. A recent outbreak destroyed more than 15 million hectares of pine forests, with major environmental effects on forest health, and economic effects on the forest industry. The outbreak has in part been driven by climate change, and will contribute to increased carbon emissions through decaying forests.

RESULTS

We developed a genome sequence resource for the mountain pine beetle to better understand the unique aspects of this insect's biology. A draft de novo genome sequence was assembled from paired-end, short-read sequences from an individual field-collected male pupa, and scaffolded using mate-paired, short-read genomic sequences from pooled field-collected pupae, paired-end short-insert whole-transcriptome shotgun sequencing reads of mRNA from adult beetle tissues, and paired-end Sanger EST sequences from various life stages. We describe the cytochrome P450, glutathione S-transferase, and plant cell wall-degrading enzyme gene families important to the survival of the mountain pine beetle in its harsh and nutrient-poor host environment, and examine genome-wide single-nucleotide polymorphism variation. A horizontally transferred bacterial sucrose-6-phosphate hydrolase was evident in the genome, and its tissue-specific transcription suggests a functional role for this beetle.

CONCLUSIONS

Despite Coleoptera being the largest insect order with over 400,000 described species, including many agricultural and forest pest species, this is only the second genome sequence reported in Coleoptera, and will provide an important resource for the Curculionoidea and other insects.