Annotated Draft Genomes of Two Caddisfly Species Plectrocnemia conspersa CURTIS and Hydropsyche tenuis NAVAS (Insecta: Trichoptera)

Members of the speciose insect order Trichoptera (caddisflies) provide important ecosystem services, for example, nutrient cycling through breaking down of organic matter. They are also of industrial interest due to their larval silk secretions. These form the basis for their diverse case-making behavior that allows them to exploit a wide range of ecological niches. Only five genomes of this order have been published thus far, with variable qualities regarding contiguity and completeness. A low-cost sequencing strategy, that is, using a single Oxford Nanopore flow cell per individual along with Illumina sequence reads was successfully used to generate high-quality genomes of two Trichoptera species, Plectrocnemia conspersa and Hydropsyche tenuis. Of the de novo assembly methods compared, assembly of low coverage Nanopore reads (∼18×) and subsequent polishing with long reads followed by Illumina short reads (∼80-170× coverage) yielded the highest genome quality both in terms of contiguity and BUSCO completeness. The presented genomes are the shortest to date and extend our knowledge of genome size across caddisfly families. The genomic region that encodes for light (L)-chain fibroin, a protein component of larval caddisfly silk was identified and compared with existing L-fibroin gene clusters. The new genomic resources presented in this paper are among the highest quality Trichoptera genomes and will increase the knowledge of this important insect order by serving as the basis for phylogenomic and comparative genomic studies.

The joy of balancers

Balancer chromosomes are multiply inverted and rearranged chromosomes that are widely used in Drosophila genetics. First described nearly 100 years ago, balancers are used extensively in stock maintenance and complex crosses. Recently, the complete molecular structures of several commonly used balancers were determined by whole-genome sequencing. This revealed a surprising amount of variation among balancers derived from a common progenitor, identified genes directly affected by inversion breakpoints, and cataloged mutations shared by balancers. These studies emphasized that it is important to choose the optimal balancer, because different inversions suppress meiotic recombination in different chromosomal regions. In this review, we provide a brief history of balancers in Drosophila, discuss how they are used today, and provide examples of unexpected recombination events involving balancers that can lead to stock breakdown.

Chromatin Architecture in the Fly: Living without CTCF/Cohesin Loop Extrusion?: Alternating Chromatin States Provide a Basis for Domain Architecture in Drosophila

The organization of the genome into topologically associated domains (TADs) appears to be a fundamental process occurring across a wide range of eukaryote organisms, and it likely plays an important role in providing an architectural foundation for gene regulation. Initial studies emphasized the remarkable parallels between TAD organization in organisms as diverse as Drosophila and mammals. However, whereas CCCTC-binding factor (CTCF)/cohesin loop extrusion is emerging as a key mechanism for the formation of mammalian topological domains, the genome organization in Drosophila appears to depend primarily on the partitioning of chromatin state domains. Recent work suggesting a fundamental conserved role of chromatin state in building domain architecture is discussed and insights into genome organization from recent studies in Drosophila are considered.

Survey of the Bradysia odoriphaga Transcriptome Using PacBio Single-Molecule Long-Read Sequencing

The damage caused by Bradysia odoriphaga is the main factor threatening the production of vegetables in the Liliaceae family. However, few genetic studies of B. odoriphaga have been conducted because of a lack of genomic resources. Many long-read sequencing technologies have been developed in the last decade; therefore, in this study, the transcriptome including all development stages of B. odoriphaga was sequenced for the first time by Pacific single-molecule long-read sequencing. Here, 39,129 isoforms were generated, and 35,645 were found to have annotation results when checked against sequences available in different databases. Overall, 18,473 isoforms were distributed in 25 various Clusters of Orthologous Groups, and 11,880 isoforms were categorized into 60 functional groups that belonged to the three main Gene Ontology classifications. Moreover, 30,610 isoforms were assigned into 44 functional categories belonging to six main Kyoto Encyclopedia of Genes and Genomes functional categories. Coding DNA sequence (CDS) prediction showed that 36,419 out of 39,129 isoforms were predicted to have CDS, and 4319 simple sequence repeats were detected in total. Finally, 266 insecticide resistance and metabolism-related isoforms were identified as candidate genes for further investigation of insecticide resistance and metabolism in B. odoriphaga.

Visualizing DNA folding and RNA in embryos at single-cell resolution

The establishment of cell types during development requires precise interactions between genes and distal regulatory sequences. We have a limited understanding of how these interactions look in three dimensions, vary across cell types in complex tissue, and relate to transcription. Here we describe optical reconstruction of chromatin architecture (ORCA), a method that can trace the DNA path in single cells with nanoscale accuracy and genomic resolution reaching two kilobases. We used ORCA to study a Hox gene cluster in cryosectioned Drosophila embryos and labelled around 30 RNA species in parallel. We identified cell-type-specific physical borders between active and Polycomb-repressed DNA, and unexpected Polycomb-independent borders. Deletion of Polycomb-independent borders led to ectopic enhancer-promoter contacts, aberrant gene expression, and developmental defects. Together, these results illustrate an approach for high-resolution, single-cell DNA domain analysis in vivo, identify domain structures that change with cell identity, and show that border elements contribute to the formation of physical domains in Drosophila.

Nanopore sequencing reads improve assembly and gene annotation of the Parochlus steinenii genome

Parochlus steinenii is a winged midge from King George Island. It is cold-tolerant and endures the harsh Antarctic winter. Previously, we reported the genome of this midge, but the genome assembly with short reads had limited contig contiguity, which reduced the completeness of the genome assembly and the annotated gene sets. Recently, assembly contiguity has been increased using nanopore technology. A number of methods for enhancing the low base quality of the assembly have been reported, including long-read (e.g. Nanopolish) or short-read (e.g. Pilon) based methods. Based on these advances, we used nanopore technologies to upgrade the draft genome sequence of P. steinenii. The final assembled genome was 145,366,448 bases in length. The contig number decreased from 9,132 to 162, and the N50 contig size increased from 36,946 to 1,989,550 bases. The BUSCO completeness of the assembly increased from 87.8 to 98.7%. Improved assembly statistics helped predict more genes from the draft genome of P. steinenii. The completeness of the predicted gene model increased from 79.5 to 92.1%, but the numbers and types of the predicted repeats were similar to those observed in the short read assembly, with the exception of long interspersed nuclear elements. In the present study, we markedly improved the P. steinenii genome assembly statistics using nanopore sequencing, but found that genome polishing with high-quality reads was essential for improving genome annotation. The number of genes predicted and the lengths of the genes were greater than before, and nanopore technology readily improved genome information.

How chromosomal rearrangements shape adaptation and speciation: Case studies in Drosophila pseudoobscura and its sibling species Drosophila persimilis

The gene arrangements of Drosophila have played a prominent role in the history of evolutionary biology from the original quantification of genetic diversity to current studies of the mechanisms for the origin and establishment of new inversion mutations within populations and their subsequent fixation between species supporting reproductive barriers. This review examines the genetic causes and consequences of inversions as recombination suppressors and the role that recombination suppression plays in establishing inversions in populations as they are involved in adaptation within heterogeneous environments. This often results in the formation of clines of gene arrangement frequencies among populations. Recombination suppression leads to the differentiation of the gene arrangements which may accelerate the accumulation of fixed genetic differences among populations. If these fixed mutations cause incompatibilities, then inversions pose important reproductive barriers between species. This review uses the evolution of inversions in Drosophila pseudoobscura and D. persimilis as a case study for how inversions originate, establish and contribute to the evolution of reproductive isolation.

Population-specific dynamics and selection patterns of transposable element insertions in European natural populations

Transposable elements (TEs) are ubiquitous sequences in genomes of virtually all species. While TEs have been investigated for several decades, only recently we have the opportunity to study their genome-wide population dynamics. Most of the studies so far have been restricted either to the analysis of the insertions annotated in the reference genome or to the analysis of a limited number of populations. Taking advantage of the European Drosophila population genomics consortium (DrosEU) sequencing data set, we have identified and measured the dynamics of TEs in a large sample of European Drosophila melanogaster natural populations. We showed that the mobilome landscape is population-specific and highly diverse depending on the TE family. In contrast with previous studies based on SNP variants, no geographical structure was observed for TE abundance or TE divergence in European populations. We further identified de novo individual insertions using two available programs and, as expected, most of the insertions were present at low frequencies. Nevertheless, we identified a subset of TEs present at high frequencies and located in genomic regions with a high recombination rate. These TEs are candidates for being the target of positive selection, although neutral processes should be discarded before reaching any conclusion on the type of selection acting on them. Finally, parallel patterns of association between the frequency of TE insertions and several geographical and temporal variables were found between European and North American populations, suggesting that TEs can be potentially implicated in the adaptation of populations across continents.

Stress response, behavior, and development are shaped by transposable element-induced mutations in Drosophila

Most of the current knowledge on the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions annotated in the reference genome of Drosophila melanogaster, in 91 samples from 60 worldwide natural populations. We identified a set of 300 polymorphic TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Besides, we identified a set of 21 fixed TEs also likely to be adaptive. Indeed, we, and others, found evidence of selection for 84 of these reference TE insertions. The analysis of the genes located nearby these 84 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affects expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well.

Artificial selection on storage protein 1 possibly contributes to increase of hatchability during silkworm domestication

Like other domesticates, the efficient utilization of nitrogen resources is also important for the only fully domesticated insect, the silkworm. Deciphering the way in which artificial selection acts on the silkworm genome to improve the utilization of nitrogen resources and to advance human-favored domestication traits, will provide clues from a unique insect model for understanding the general rules of Darwin's evolutionary theory on domestication. Storage proteins (SPs), which belong to a hemocyanin superfamily, basically serve as a source of amino acids and nitrogen during metamorphosis and reproduction in insects. In this study, through blast searching on the silkworm genome and further screening of the artificial selection signature on silkworm SPs, we discovered a candidate domestication gene, i.e., the methionine-rich storage protein 1 (SP1), which is clearly divergent from other storage proteins and exhibits increased expression in the ova of domestic silkworms. Knockout of SP1 via the CRISPR/Cas9 technique resulted in a dramatic decrease in egg hatchability, without obvious impact on egg production, which was similar to the effect in the wild silkworm compared with the domestic type. Larval development and metamorphosis were not affected by SP1 knockout. Comprehensive ova comparative transcriptomes indicated significant higher expression of genes encoding vitellogenin, chorions, and structural components in the extracellular matrix (ECM)-interaction pathway, enzymes in folate biosynthesis, and notably hormone synthesis in the domestic silkworm, compared to both the SP1 mutant and the wild silkworm. Moreover, compared with the wild silkworms, the domestic one also showed generally up-regulated expression of genes enriched in the structural constituent of ribosome and amide, as well as peptide biosynthesis. This study exemplified a novel case in which artificial selection could act directly on nitrogen resource proteins, further affecting egg nutrients and eggshell formation possibly through a hormone signaling mediated regulatory network and the activation of ribosomes, resulting in improved biosynthesis and increased hatchability during domestication. These findings shed new light on both the understanding of artificial selection and silkworm breeding from the perspective of nitrogen and amino acid resources.

Like other domesticates, nitrogen resources are also important for the only fully domesticated insect, the silkworm. Deciphering the way in which artificial selection acts on the silkworm genome to improve the utilization of nitrogen resources, thereby advancing human-favored domestication traits, will provide clues from a unique insect model for understanding the general rules of Darwin's theory on artificial selection. However, the mechanisms of domestication in the silkworm remain largely unknown. In this study, we focused on one important nitrogen resource, the storage protein (SP). We discovered that the methionine-rich storage protein 1 (SP1), which is divergent from other SPs, is the only target of artificial selection. Based on functional evidence, together with key findings from the comprehensive comparative transcriptome, we propose that artificial selection favored higher expression of SP1 in the domestic silkworm, which would influence the genes or pathways vital for egg development and eggshell formation. Artificial selection also consistently favored activated ribosome activities and improved amide and peptide biosynthesis in the ova, like what they may act in the silk gland to increase silk-cocoon yield. We highlighted a novel case in which artificial selection could directly act on a nitrogen resource protein associated with a human-desired domestication trait.

The first complete mitochondrial genome of marigold pest thrips, Neohydatothrips samayunkur (Sericothripinae) and comparative analysis

Complete mitogenomes from the order Thysanoptera are limited to representatives of the subfamily Thripinae. Therefore, in the present study, we sequenced the mitochondrial genome of Neohydatothrips samayunkur (15,295 bp), a member of subfamily Sericothripinae. The genome possesses the canonical 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs) as well as two putative control regions (CRs). The majority strand was 77.42% A + T content, and 22.58% G + C with weakly positive AT skew (0.04) and negative GC skew (-0.03). The majority of PCGs start with ATN codons as observed in other insect mitochondrial genomes. The GCG codon (Alanine) was not used in N. samayunkur. Most tRNAs have the typical cloverleaf secondary structure, however the DHU stem and loop were absent in trnV and trnS1, while the TΨC loop was absent in trnR and trnT. The two putative control regions (CR1 and CR2) show 99% sequence similarity indicated a possible duplication, and shared 57 bp repeats were identified. N. samayunkur showed extensive gene rearrangements, with 11 PCGs, 22 tRNAs, and two rRNAs translocated when compared to the ancestral insect. The gene trnL2 was separated from the 'trnL2-cox2' gene block, which is a conserved, ancestral gene order found in all previously sequenced thrips mitogenomes. Both maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees resulted in similar topologies. The phylogenetic position of N. samayunkur indicates that subfamily Sericothripinae is sister to subfamily Thripinae. More molecular data from different taxonomic groups is needed to understand thrips phylogeny and evolution.

Wolbachia Infection Associated with Increased Recombination in Drosophila

Wolbachia is a maternally-transmitted endosymbiotic bacteria that infects a large diversity of arthropod and nematode hosts. Some strains of Wolbachia are parasitic, manipulating host reproduction to benefit themselves, while other strains of Wolbachia exhibit obligate or facultative mutualisms with their host. The effects of Wolbachia on its host are many, though primarily relate to host immune and reproductive function. Here we test the hypothesis that Wolbachia infection alters the frequency of homologous recombination during meiosis. We use D. melanogaster as a model system, and survey recombination in eight wild-derived Wolbachia-infected (strain wMel) and Wolbachia-uninfected strains, controlling for genotype. We measure recombination in two intervals of the genome. Our results indicate that Wolbachia infection is associated with increased recombination in one genomic interval and not the other. The effect of Wolbachia infection on recombination is thus heterogenous across the genome. Our data also indicate a reproductive benefit of Wolbachia infection; infected females show higher fecundity than their uninfected genotypic controls. Given the prevalence of Wolbachia infection in natural populations, our findings suggest that Wolbachia infection is likely to contribute to recombination rate and fecundity variation among individuals in nature.

Transforming insect population control with precision guided sterile males with demonstration in flies

The sterile insect technique (SIT) is an environmentally safe and proven technology to suppress wild populations. To further advance its utility, a novel CRISPR-based technology termed precision guided SIT (pgSIT) is described. PgSIT mechanistically relies on a dominant genetic technology that enables simultaneous sexing and sterilization, facilitating the release of eggs into the environment ensuring only sterile adult males emerge. Importantly, for field applications, the release of eggs will eliminate burdens of manually sexing and sterilizing males, thereby reducing overall effort and increasing scalability. Here, to demonstrate efficacy, we systematically engineer multiple pgSIT systems in Drosophila which consistently give rise to 100% sterile males. Importantly, we demonstrate that pgSIT-generated sterile males are fit and competitive. Using mathematical models, we predict pgSIT will induce substantially greater population suppression than can be achieved by currently-available self-limiting suppression technologies. Taken together, pgSIT offers to potentially transform our ability to control insect agricultural pests and disease vectors.

Spatiotemporally controlled genetic perturbation for efficient large-scale studies of cell non-autonomous effects

Studies in genetic model organisms have revealed much about the development and pathology of complex tissues. Most have focused on cell-intrinsic gene functions and mechanisms. Much less is known about how transformed, or otherwise functionally disrupted, cells interact with healthy ones toward a favorable or pathological outcome. This is largely due to technical limitations. We developed new genetic tools in Drosophila melanogaster that permit efficient multiplexed gain- and loss-of-function genetic perturbations with separable spatial and temporal control. Importantly, our novel tool-set is independent of the commonly used GAL4/UAS system, freeing the latter for additional, non-autonomous, genetic manipulations; and is built into a single strain, allowing one-generation interrogation of non-autonomous effects. Altogether, our design opens up efficient genome-wide screens on any deleterious phenotype, once plasmid or genome engineering is used to place the desired miRNA(s) or ORF(s) into our genotype. Specifically, we developed tools to study extrinsic effects on neural tumor growth but the strategy presented has endless applications within and beyond neurobiology, and in other model organisms.

Genome-wide patterns of local adaptation in Western European Drosophila melanogaster natural populations

Signatures of spatially varying selection have been investigated both at the genomic and transcriptomic level in several organisms. In Drosophila melanogaster, the majority of these studies have analyzed North American and Australian populations, leading to the identification of several loci and traits under selection. However, several studies based mainly in North American populations showed evidence of admixture that likely contributed to the observed population differentiation patterns. Thus, disentangling demography from selection might be challenging when analyzing these populations. European populations could help identify loci under spatially varying selection provided that no recent admixture from African populations would have occurred. In this work, we individually sequence the genome of 42 European strains collected in populations from contrasting environments: Stockholm (Sweden) and Castellana Grotte (Southern Italy). We found low levels of population structure and no evidence of recent African admixture in these two populations. We thus look for patterns of spatially varying selection affecting individual genes and gene sets. Besides single nucleotide polymorphisms, we also investigated the role of transposable elements in local adaptation. We concluded that European populations are a good dataset to identify candidate loci under spatially varying selection. The analysis of the two populations sequenced in this work in the context of all the available D. melanogaster data allowed us to pinpoint genes and biological processes likely to be relevant for local adaptation. Identifying and analyzing populations with low levels of population structure and admixture should help to disentangle selective from non-selective forces underlying patterns of population differentiation in other species as well.

Firefly genomes illuminate parallel origins of bioluminescence in beetles

Fireflies and their luminous courtships have inspired centuries of scientific study. Today firefly luciferase is widely used in biotechnology, but the evolutionary origin of bioluminescence within beetles remains unclear. To shed light on this long-standing question, we sequenced the genomes of two firefly species that diverged over 100 million-years-ago: the North American Photinus pyralis and Japanese Aquatica lateralis. To compare bioluminescent origins, we also sequenced the genome of a related click beetle, the Caribbean Ignelater luminosus, with bioluminescent biochemistry near-identical to fireflies, but anatomically unique light organs, suggesting the intriguing hypothesis of parallel gains of bioluminescence. Our analyses support independent gains of bioluminescence in fireflies and click beetles, and provide new insights into the genes, chemical defenses, and symbionts that evolved alongside their luminous lifestyle.

Using genomic data to study insecticide resistance in the house fly, Musca domestica

The house fly, Musca domestica, is a major pest at livestock facilities throughout the world. Insecticides have been the most common control strategy for flies, but many populations have evolved resistance. The speed by which we are able to identify the mutations responsible for resistance has been a major challenge for the development of high throughput resistance monitoring assays as new insecticides are introduced for control. This is particularly true for mutations that cause trans regulation of a gene, which then results in resistance. In this paper we take advantage of the conserved homology of dipteran chromosomes to assign 3069 genes to chromosomes. Of these, 234 were of toxicological interest (CYPs, esterases/hydrolases, glutathione S-transferases (GSTs) and target sites). The chromosomal location of genes known from linkage analysis studies matched the location predicted by homology mapping in ten out of ten cases, indicating a high reliability of our approach. The CYPs, esterases/hydrolases and GSTs were not randomly distributed throughout the genome. They clustered on chromosomes, but the pattern was different between the CYPs, esterases/hydrolases and GSTs. Examples are provided for how the availability of the house fly genome, combined with an ability to assign genes to chromosomes, will help to accelerate research in house flies.

The origin of the odorant receptor gene family in insects

The origin of the insect odorant receptor (OR) gene family has been hypothesized to have coincided with the evolution of terrestriality in insects. Missbach et al. (2014) suggested that ORs instead evolved with an ancestral OR co-receptor (Orco) after the origin of terrestriality and the OR/Orco system is an adaptation to winged flight in insects. We investigated genomes of the Collembola, Diplura, Archaeognatha, Zygentoma, Odonata, and Ephemeroptera, and find ORs present in all insect genomes but absent from lineages predating the evolution of insects. Orco is absent only in the ancestrally wingless insect lineage Archaeognatha. Our new genome sequence of the zygentoman firebrat Thermobia domestica reveals a full OR/Orco system. We conclude that ORs evolved before winged flight, perhaps as an adaptation to terrestriality, representing a key evolutionary novelty in the ancestor of all insects, and hence a molecular synapomorphy for the Class Insecta.

Ancestral polymorphisms explain the role of chromosomal inversions in speciation

Understanding the role of chromosomal inversions in speciation is a fundamental problem in evolutionary genetics. Here, we perform a comprehensive reconstruction of the evolutionary histories of the chromosomal inversions in Drosophila persimilis and D. pseudoobscura. We provide a solution to the puzzling origins of the selfish Sex-Ratio arrangement in D. persimilis and uncover surprising patterns of phylogenetic discordance on this chromosome. These patterns show that, contrary to widely held views, all fixed chromosomal inversions between D. persimilis and D. pseudoobscura were already present in their ancestral population long before the species split. Our results suggest that patterns of higher genomic divergence and an association of reproductive isolation genes with chromosomal inversions may be a direct consequence of incomplete lineage sorting of ancestral polymorphisms. These findings force a reconsideration of the role of chromosomal inversions in speciation, not as protectors of existing hybrid incompatibilities, but as fertile grounds for their formation.

The Complete Mitochondrial Genome of Ugyops sp. (Hemiptera: Delphacidae)

The complete mitochondrial genome (mitogenome) of Ugyops sp. (Hemiptera: Delphacidae) was sequenced, making it the first determined mitogenome from the subfamily Asiracinae, the basal clade of the family Delphacidae. The mitogenome was 15,259 bp in length with A + T content of 77.65% and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region. The gene order was identical with that of the ancestral insect. The nucleotide composition analysis indicated that the whole mitogenome was strongly A-skewed (0.288) and highly C-skewed (-0.270). For PCGs on the J-strand, the AT skew was positive, and the GC skew was negative. All PCGs started with canonical ATN codons, except for cox1 and nad5, which used CTG and GTG as start codon, respectively. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnS1 (AGN), in which the dihydrouridine arm was reduced to a simple loop. The control region included a poly-T stretch downstream of the small rRNA gene (rrnS), a subregion of higher A + T content and tandemly repeated sequence near trnI. The mitogenome of Ugyops sp. could be very helpful in exploring the diversity and evolution of mitogenomes in Delphacidae.

Chromatin immunoprecipitation (ChIP) method for non-model fruit flies (Diptera: Tephritidae) and evidence of histone modifications

Interactions between DNA and proteins located in the cell nucleus play an important role in controlling physiological processes by specifying, augmenting and regulating context-specific transcription events. Chromatin immunoprecipitation (ChIP) is a widely used methodology to study DNA-protein interactions and has been successfully used in various cell types for over three decades. More recently, by combining ChIP with genomic screening technologies and Next Generation Sequencing (e.g. ChIP-seq), it has become possible to profile DNA-protein interactions (including covalent histone modifications) across entire genomes. However, the applicability of ChIP-chip and ChIP-seq has rarely been extended to non-model species because of a number of technical challenges. Here we report a method that can be used to identify genome wide covalent histone modifications in a group of non-model fruit fly species (Diptera: Tephritidae). The method was developed by testing and refining protocols that have been used in model organisms, including Drosophila melanogaster. We demonstrate that this method is suitable for a group of economically important pest fruit fly species, viz., Bactrocera dorsalis, Ceratitis capitata, Zeugodacus cucurbitae and Bactrocera tryoni. We also report an example ChIP-seq dataset for B. tryoni, providing evidence for histone modifications in the genome of a tephritid fruit fly for the first time. Since tephritids are major agricultural pests globally, this methodology will be a valuable resource to study taxa-specific evolutionary questions and to assist with pest management. It also provides a basis for researchers working with other non-model species to undertake genome wide DNA-protein interaction studies.

Measuring Sister Chromatid Cohesion Protein Genome Occupancy in Drosophila melanogaster by ChIP-seq

This chapter presents methods to conduct and analyze genome-wide chromatin immunoprecipitation of the cohesin complex and the Nipped-B cohesin loading factor in Drosophila cells using high-throughput DNA sequencing (ChIP-seq). Procedures for isolation of chromatin, immunoprecipitation, and construction of sequencing libraries for the Ion Torrent Proton high throughput sequencer are detailed, and computational methods to calculate occupancy as input-normalized fold-enrichment are described. The results obtained by ChIP-seq are compared to those obtained by ChIP-chip (genomic ChIP using tiling microarrays), and the effects of sequencing depth on the accuracy are analyzed. ChIP-seq provides similar sensitivity and reproducibility as ChIP-chip, and identifies the same broad regions of occupancy. The locations of enrichment peaks, however, can differ between ChIP-chip and ChIP-seq, and low sequencing depth can splinter broad regions of occupancy into distinct peaks.

Genome-wide identification of long non-coding RNA genes and their association with insecticide resistance and metamorphosis in diamondback moth, Plutella xylostella

Long non-coding RNA (lncRNA) is a class of noncoding RNA >200 bp in length that has essential roles in regulating a variety of biological processes. Here, we constructed a computational pipeline to identify lncRNA genes in the diamondback moth (Plutella xylostella), a major insect pest of cruciferous vegetables. In total, 3,324 lncRNAs corresponding to 2,475 loci were identified from 13 RNA-Seq datasets, including samples from parasitized, insecticide-resistant strains and different developmental stages. The identified P. xylostella lncRNAs had shorter transcripts and fewer exons than protein-coding genes. Seven out of nine randomly selected lncRNAs were validated by strand-specific RT-PCR. In total, 54-172 lncRNAs were specifically expressed in the insecticide resistant strains, among which one lncRNA was located adjacent to the sodium channel gene. In addition, 63-135 lncRNAs were specifically expressed in different developmental stages, among which three lncRNAs overlapped or were located adjacent to the metamorphosis-associated genes. These lncRNAs were either strongly or weakly co-expressed with their overlapping or neighboring mRNA genes. In summary, we identified thousands of lncRNAs and presented evidence that lncRNAs might have key roles in conferring insecticide resistance and regulating the metamorphosis development in P. xylostella.

Large-scale analysis reveals that the genome features of simple sequence repeats are generally conserved at the family level in insects

BACKGROUND

Simple sequence repeats (SSR), also called microsatellites, have been widely used as genetic markers, and have been extensively studied in some model insects. At present, the genomes of more than 100 insect species are available. However, the features of SSRs in most insect genomes remain largely unknown.

RESULTS

We identified 15.01 million SSRs across 136 insect genomes. The number of identified SSRs was positively associated with genome size in insects, but the frequency and density per megabase of genomes were not. Most insect SSRs (56.2-93.1%) were perfect (no mismatch). Imperfect (at least one mismatch) SSRs (average length 22-73 bp) were longer than perfect SSRs (16-30 bp). The most abundant insect SSRs were the di- and trinucleotide types, which accounted for 27.2% and 22.0% of all SSRs, respectively. On average, 59.1%, 36.8%, and 3.7% of insect SSRs were located in intergenic, intronic, and exonic regions, respectively. The percentages of various types of SSRs were similar among insects from the same family. However, they were dissimilar among insects from different families within orders. We carried out a phylogenetic analysis using the SSR frequencies. Species from the same family were generally clustered together in the evolutionary tree. However, insects from the same order but not in the same family did not cluster together. These results indicated that although SSRs undergo rapid expansions and contractions in different populations of the same species, the general genomic features of insect SSRs remain conserved at the family level.

CONCLUSION

Millions of insect SSRs were identified and their genome features were analyzed. Most insect SSRs were perfect and were located in intergenic regions. We presented evidence that the variance of insect SSRs accumulated after the differentiation of insect families.

Comparative Gene Mapping as a Tool to Understand the Evolution of Pest Crop Insect Chromosomes

The extent of the conservation of synteny and gene order in aphids has been previously investigated only by comparing a small subset of linkage groups between the pea aphid Acyrthosiphon pisum and a few other aphid species. Here we compared the localization of eight A. pisum scaffolds (covering more than 5 Mb and 83 genes) in respect to the Drosophila melanogaster Muller elements identifying orthologous loci spanning all the four A. pisum chromosomes. Comparison of the genetic maps revealed a conserved synteny across different loci suggesting that the study of the fruit fly Muller elements could favour the identification of chromosomal markers useful for the study of chromosomal rearrangements in aphids. A. pisum is the first aphid species to have its genome sequenced and the finding that there are several chromosomal regions in synteny between Diptera and Hemiptera indicates that the genomic tools developed in A. pisum will be broadly useful not only for the study of other aphids but also for other insect species.