Genome-Wide Identification, Characterization and Expression Analysis of the Solute Carrier 6 Gene Family in Silkworm (Bombyx mori)

The solute carrier 6 (SLC6) gene family, initially known as the neurotransmitter transporters, plays vital roles in the regulation of neurotransmitter signaling, nutrient absorption and motor behavior. In this study, a total of 16 candidate genes were identified as SLC6 family gene homologs in the silkworm (Bombyx mori) genome. Spatio-temporal expression patterns of silkworm SLC6 gene transcripts indicated that these genes were highly and specifically expressed in midgut, brain and gonads; moreover, these genes were expressed primarily at the feeding stage or adult stage. Levels of expression for most midgut-specific and midgut-enriched gene transcripts were down-regulated after starvation but up-regulated after re-feeding. In addition, we observed that expression levels of these genes except for BmSLC6-15 and BmGT1 were markedly up-regulated by a juvenile hormone analog. Moreover, brain-enriched genes showed differential expression patterns during wandering and mating processes, suggesting that these genes may be involved in modulating wandering and mating behaviors. Our results improve our understanding of the expression patterns and potential physiological functions of the SLC6 gene family, and provide valuable information for the comprehensive functional analysis of the SLC6 gene family.

Identification and characteristic analysis of the catalase gene from Locusta migratoria

Catalase (CAT) is a ubiquitous antioxidant enzyme in almost all living organisms exposed to atmosphere, which involved in decomposing harmful hydrogen peroxide, into oxygen and water. In this study, a full-length cDNA (1524bp) encoding the catalase gene (LmCAT) from Locusta migratoria was cloned (accession number KT716445). The open reading frame of the LmCAT gene encoded 507 amino acids and shared 57.8%-97.8% amino acid identities with other insect CATs. The coding region was interrupted by 9 introns, while its promoter region contained 15 putative binding sites for 5 kinds of transcriptional regulation factors. For the stage-specific expression profile, LmCAT was highly expressed in the fourth-instar nymphs. For the tissue-specific expression profile, the LmCAT transcripts were highest in the fat bodies, and relatively abundant in the gastric caecum, Malpighian tubules, ovary and integument. Moreover, the result showed that quercetin could significantly induce the expression level of LmCAT. The expression of LmCAT could be silenced by RNAi, but the moralities were not significantly different between control and RNAi groups. Our results would provide valuable information for further study on the ROS regulation mechanism in insect.

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.

RNA-sequencing and whole-genome bisulfite sequencing in the hybrid offspring of two Nasonia parasitoid wasp species revealed strong cis-regulation of methylation and allele-specific expression. No gene was found to display genomic imprinting.

The relationship between methylation of genomic DNA and expression of the genes that it encodes—and how this relationship changes during evolution—has been widely studied in mammals, but remains less well understood for insects. Here we analyze the expressed mRNA transcripts and genomic DNA methylation of the hybrid offspring of a pair of Nasonia parasitoid wasp species, producing a wealth of information about the regulation of gene expression. We find that variation in DNA sequence impacts expression on the same strand (called “cis-regulation”), and that cytosine methylation state is also associated in cis with the regulatory consequences of this base alteration. We show that these wasp species lack differential expression dependent on parent-of-origin (called “genomic imprinting”), and that in the hybrids the alleles retain the methylation status of the parental species in a strong cis-regulated fashion. Transcript abundances were also largely driven in a cis-regulated manner, consistent with a correlation between methylation status and expression levels. Despite the many differences between Nasonia and mammals in the impact of genomic DNA methylation, in both groups the use of methylated cytosine has been co-opted in ways that help tune gene expression.

Spatial distribution of predicted transcription factor binding sites in Drosophila ChIP peaks

In the development of the Drosophila embryo, gene expression is directed by the sequence-specific interactions of a large network of protein transcription factors (TFs) and DNA cis-regulatory binding sites. Once the identity of the typically 8-10bp binding sites for any given TF has been determined by one of several experimental procedures, the sequences can be represented in a position weight matrix (PWM) and used to predict the location of additional TF binding sites elsewhere in the genome. Often, alignments of large (>200bp) genomic fragments that have been experimentally determined to bind the TF of interest in Chromatin Immunoprecipitation (ChIP) studies are trimmed under the assumption that the majority of the binding sites are located near the center of all the aligned fragments. In this study, ChIP/chip datasets are analyzed using the corresponding PWMs for the well-studied TFs; CAUDAL, HUNCHBACK, KNIRPS and KRUPPEL, to determine the distribution of predicted binding sites. All four TFs are critical regulators of gene expression along the anterio-posterior axis in early Drosophila development. For all four TFs, the ChIP peaks contain multiple binding sites that are broadly distributed across the genomic region represented by the peak, regardless of the prediction stringency criteria used. This result suggests that ChIP peak trimming may exclude functional binding sites from subsequent analyses.

Identification of Bari Transposons in 23 Sequenced Drosophila Genomes Reveals Novel Structural Variants, MITEs and Horizontal Transfer

Bari elements are members of the Tc1-mariner superfamily of DNA transposons, originally discovered in Drosophila melanogaster, and subsequently identified in silico in 11 sequenced Drosophila genomes and as experimentally isolated in four non-sequenced Drosophila species. Bari-like elements have been also studied for their mobility both in vivo and in vitro. We analyzed 23 Drosophila genomes and carried out a detailed characterization of the Bari elements identified, including those from the heterochromatic Bari1 cluster in D. melanogaster. We have annotated 401 copies of Bari elements classified either as putatively autonomous or inactive according to the structure of the terminal sequences and the presence of a complete transposase-coding region. Analyses of the integration sites revealed that Bari transposase prefers AT-rich sequences in which the TA target is cleaved and duplicated. Furthermore evaluation of transposon’s co-occurrence near the integration sites of Bari elements showed a non-random distribution of other transposable elements. We also unveil the existence of a putatively autonomous Bari1 variant characterized by two identical long Terminal Inverted Repeats, in D. rhopaloa. In addition, we detected MITEs related to Bari transposons in 9 species. Phylogenetic analyses based on transposase gene and the terminal sequences confirmed that Bari-like elements are distributed into three subfamilies. A few inconsistencies in Bari phylogenetic tree with respect to the Drosophila species tree could be explained by the occurrence of horizontal transfer events as also suggested by the results of dS analyses. This study further clarifies the Bari transposon’s evolutionary dynamics and increases our understanding on the Tc1-mariner elements’ biology.

Power, resolution and bias: recent advances in insect phylogeny driven by the genomic revolution

Our understanding on the phylogenetic relationships of insects has been revolutionised in the last decade by the proliferation of next generation sequencing technologies (NGS). NGS has allowed insect systematists to assemble very large molecular datasets that include both model and non-model organisms. Such datasets often include a large proportion of the total number of protein coding sequences available for phylogenetic comparison. We review some early entomological phylogenomic studies that employ a range of different data sampling protocols and analyses strategies, illustrating a fundamental renaissance in our understanding of insect evolution all driven by the genomic revolution. The analysis of phylogenomic datasets is challenging because of their size and complexity, and it is obvious that the increasing size alone does not ensure that phylogenetic signal overcomes systematic biases in the data. Biases can be due to various factors such as the method of data generation and assembly, or intrinsic biological feature of the data per se, such as similarities due to saturation or compositional heterogeneity. Such biases often cause violations in the underlying assumptions of phylogenetic models. We review some of the bioinformatics tools available and being developed to detect and minimise systematic biases in phylogenomic datasets. Phylogenomic-scale data coupled with sophisticated analyses will revolutionise our understanding of insect functional genomics. This will illuminate the relationship between the vast range of insect phenotypic diversity and underlying genetic diversity. In combination with rapidly developing methods to estimate divergence times, these analyses will also provide a compelling view of the rates and patterns of lineagenesis (birth of lineages) over the half billion years of insect evolution.

Agricultural applications of insect ecological genomics

Agricultural entomology is poised to benefit from the application of ecological genomics, particularly the fields of biofuels generation and pest control. Metagenomic methods can characterize microbial communities of termites, wood-boring beetles and livestock pests, and transcriptomic approaches reveal molecular bases behind wood-digesting capabilities of these insects, leading to potential mechanisms for biofuel generation. Genome sequences are being exploited to develop new pest control methods, identify candidate antigens to vaccinate livestock, and discover RNAi target sequences and potential non-target effects in other insects. Gene content analyses of pest genome sequences and their endosymbionts suggest metabolic interdependencies between organisms, exposing potential gene targets for insect control. Finally, genome-wide association studies and genotyping by high-throughput sequencing promise to improve management of pesticide resistance.

Genomic sources of phenotypic novelty in the evolution of eusociality in insects

Genomic resources are now available for closely related species that vary in social behavior, providing insight on the genomics of social evolution. Changes in the architecture of gene regulatory networks likely influence the evolutionary trajectory of social traits. Evolutionarily novel genes are likely important in the evolution of social diversity among insects, but it is unclear whether new genes played a driving role in the advent or elaboration of eusociality or if they were instead a result of other genomic features of eusociality. The worker phenotype appears to be the center of genetic novelty, but the mechanisms for this remain unresolved. Future studies are needed to understand how genetic novelty arises, becomes incorporated into existing gene regulatory networks, and the effects this has on the evolution of social traits in closely related social and solitary species.

Applications of genome editing in insects

Insect genome editing was first reported 1991 in Drosophila melanogaster but the technology used was not portable to other species. Not until the recent development of facile, engineered DNA endonuclease systems has gene editing become widely available to insect scientists. Most applications in insects to date have been technical in nature but this is rapidly changing. Functional genomics and genetics-based insect control efforts will be major beneficiaries of the application of contemporary gene editing technologies. Engineered endonucleases like Cas9 make it possible to create powerful and effective gene drive systems that could be used to reduce or even eradicate specific insect populations. 'Best practices' for using Cas9-based editing are beginning to emerge making it easier and more effective to design and use but gene editing technologies still require traditional means of delivery in order to introduce them into somatic and germ cells of insects-microinjection of developing embryos. This constrains the use of these technologies by insect scientists. Insects created using editing technologies challenge existing governmental regulatory structures designed to manage genetically modified organisms.

The Complete Mitochondrial Genome of Brachmia macroscopa (Lepidoptera: Gelechiidae) and Its Related Phylogenetic Analysis

The sweet potato leaf folder, Brachmia macroscopa, is an important pest in China. The complete mitogenome, which consists of 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and an A + T-rich region, was sequenced and found to be 15,394 bp in length (GeneBank no. KT354968). The gene order and orientation of the B. macroscopa mitogenome were similar to those of other sequenced lepidopteran species. All of the PCGs started with ATN as the canonical start codon except for cox1, which started with CGA. In regard to stop codons, most PCGs stopped at TAA except for cox2, which stopped at TA, and nad4, which stopped at a single T. Thirteen PCGs of the available species (33 taxa) were used to demonstrate phylogenetic relationships. The ditrysian cluster was supported as a monophyletic clade at high levels by using maximum likelihood and Bayesian methods. The apoditrysian group, covering the Gelechioidea, formed a monophyletic clade with a bootstrap value of 88% and a posterior probability of 1.00. The superfamily Gelechioidea was supported as a monophyletic lineage by a posterior probability of 1.00.

Systematic Identification and Characterization of Long Non-Coding RNAs in the Silkworm, Bombyx mori

Long noncoding RNAs (lncRNAs) are emerging as important regulators in various biological processes. However, to date, no systematic characterization of lncRNAs has been reported in the silkworm Bombyx mori. In the present study, we generated eighteen RNA-seq datasets with relatively high depth. Using an in-house designed lncRNA identification pipeline, 11,810 lncRNAs were identified for 5,556 loci. Among these lncRNAs, 474 transcripts were intronic lncRNAs (ilncRNAs), 6,250 transcripts were intergenic lncRNAs (lincRNAs), and 5,086 were natural antisense lncRNAs (lncNATs). Compared with protein-coding mRNAs, silkworm lncRNAs are shorter in terms of full length but longer in terms of exon and intron length. In addition, lncRNAs exhibit a lower level of sequence conservation, more repeat sequences overlapped and higher tissue specificity than protein-coding mRNAs in the silkworm. We found that 69 lncRNA transcripts from 33 gene loci may function as miRNA precursors, and 104 lncRNA transcripts from 72 gene loci may act as competing endogenous RNAs (ceRNAs). In total, 49.47% of all gene loci (2,749/5,556) for which lncRNAs were identified showed sex-biased expression. Co-expression network analysis resulted in 19 modules, 12 of which revealed relatively high tissue specificity. The highlighted darkgoldenrod module was specifically associated with middle and posterior silk glands, and the hub lncRNAs within this module were co-expressed with proteins involved in translation, translocation, and secretory processes, suggesting that these hub lncRNAs may function as regulators of the biosynthesis, translocation, and secretion of silk proteins. This study presents the first comprehensive genome-wide analysis of silkworm lncRNAs and provides an invaluable resource for genetic, evolutionary, and genomic studies of B. mori.

Methods for TALEN-Mediated Genomic Manipulations in Drosophila

TALEN (transcription activator-like effector nuclease) is a powerful tool for gene disruption and other genomic modifications. In the past 3 years or so, it has attracted eyes from every corner of the biological world, due to its characteristics of simplicity, high efficiency, low toxicity, and applicability across almost all species. In our lab, we first reported the TALEN-mediated gene disruption in Drosophila, and recently employed this technique in precisely modifying the Drosophila genome, such as in vivo tagging and gene correction. Here, we describe in detail the protocols and experiences in TALEN-mediated genomic modifications to share with the Drosophilists all over the world.

The Complete Mitochondrial Genome of two Tetragnatha Spiders (Araneae: Tetragnathidae): Severe Truncation of tRNAs and Novel Gene Rearrangements in Araneae

Mitogenomes can provide information for phylogenetic analysis and evolutionary biology. The Araneae is one of the largest orders of Arachnida with great economic importance. In order to develop mitogenome data for this significant group, we determined the complete mitogenomes of two long jawed spiders Tetragnatha maxillosa and T. nitens and performed the comparative analysis with previously published spider mitogenomes. The circular mitogenomes are 14578 bp long with A+T content of 74.5% in T. maxillosa and 14639 bp long with A+T content of 74.3% in T. nitens, respectively. Both the mitogenomes contain a standard set of 37 genes and an A+T-rich region with the same gene orientation as the other spider mitogenomes, with the exception of the different gene order by the rearrangement of two tRNAs (trnW and trnG). Most of the tRNAs lose TΨC arm stems and have unpaired amino acid acceptor arms. As interesting features, both trnS(AGN) and trnS(UCN) lack the dihydrouracil (DHU) arm and long tandem repeat units are presented in the A+T-rich region of both the spider mitogenomes. The phylogenetic relationships of 23 spider mitogenomes based on the concatenated nucleotides sequences of 13 protein-coding genes indicated that the mitogenome sequences could be useful in resolving higher-level relationship of Araneae. The molecular information acquired from the results of this study should be very useful for future researches on mitogenomic evolution and genetic diversities in spiders.

Identification, Validation, and Application of Molecular Diagnostics for Insecticide Resistance in Malaria Vectors

Insecticide resistance is a major obstacle to control of Anopheles malaria mosquitoes in sub-Saharan Africa and requires an improved understanding of the underlying mechanisms. Efforts to discover resistance genes and DNA markers have been dominated by candidate gene and quantitative trait locus studies of laboratory strains, but with greater availability of genome sequences a shift toward field-based agnostic discovery is anticipated. Mechanisms evolve continually to produce elevated resistance yielding multiplicative diagnostic markers, co-screening of which can give high predictive value. With a shift toward prospective analyses, identification and screening of resistance marker panels will boost monitoring and programmatic decision making.

Phylogenomic analyses resolve an ancient trichotomy at the base of Ischyropsalidoidea (Arachnida, Opiliones) despite high levels of gene tree conflict and unequal minority resolution frequencies

Phylogenetic resolution of ancient rapid radiations has remained problematic despite major advances in statistical approaches and DNA sequencing technologies. Here we report on a combined phylogenetic approach utilizing transcriptome data in conjunction with Sanger sequence data to investigate a tandem of ancient divergences in the harvestmen superfamily Ischyropsalidoidea (Arachnida, Opiliones, Dyspnoi). We rely on Sanger sequences to resolve nodes within and between closely related genera, and use RNA-seq data from a subset of taxa to resolve a short and ancient internal branch. We use several analytical approaches to explore this succession of ancient diversification events, including concatenated and coalescent-based analyses and maximum likelihood gene trees for each locus. We evaluate the robustness of phylogenetic inferences using a randomized locus sub-sampling approach, and find congruence across these methods despite considerable incongruence across gene trees. Incongruent gene trees are not recovered in frequencies expected from a simple multispecies coalescent model, and we reject incomplete lineage sorting as the sole contributor to gene tree conflict. Using these approaches we attain robust support for higher-level phylogenetic relationships within Ischyropsalidoidea.

The Complete Mitochondrial Genome of Delia antiqua and Its Implications in Dipteran Phylogenetics

Delia antiqua is a major underground agricultural pest widely distributed in Asia, Europe and North America. In this study, we sequenced and annotated the complete mitochondrial genome of this species, which is the first report of complete mitochondrial genome in the family Anthomyiidae. This genome is a double-stranded circular molecule with a length of 16,141 bp and an A+T content of 78.5%. It contains 37 genes (13 protein-coding genes, 22 tRNAs and 2 rRNAs) and a non-coding A+T rich region or control region. The mitochondrial genome of Delia antiqua presents a clear bias in nucleotide composition with a positive AT-skew and a negative GC-skew. All of the 13 protein-coding genes use ATN as an initiation codon except for the COI gene that starts with ATCA. Most protein-coding genes have complete termination codons but COII and ND5 that have the incomplete termination codon T. This bias is reflected in both codon usage and amino acid composition. The protein-coding genes in the D. antiqua mitochondrial genome prefer to use the codon UUA (Leu). All of the tRNAs have the typical clover-leaf structure, except for tRNA^Ser(AGN) that does not contain the dihydrouridine (DHU) arm like in many other insects. There are 7 mismatches with U-U in the tRNAs. The location and structure of the two rRNAs are conservative and stable when compared with other insects. The control region between 12S rRNA and tRNA^Ile has the highest A+T content of 93.7% in the D. antiqua mitochondrial genome. The control region includes three kinds of special regions, two highly conserved poly-T stretches, a (TA)_n stretch and several G(A)_nT structures considered important elements related to replication and transcription. The nucleotide sequences of 13 protein-coding genes are used to construct the phylogenetics of 26 representative Dipteran species. Both maximum likelihood and Bayesian inference analyses suggest a closer relationship of D. antiqua in Anthomyiidae with Calliphoridae, Calliphoridae is a paraphyly, and both Oestroidea and Muscoidea are polyphyletic.

Transcriptomic Analysis of the Anterior Silk Gland in the Domestic Silkworm (Bombyx mori) - Insight into the Mechanism of Silk Formation and Spinning

Silk proteins are synthesized in the middle and posterior silk glands of silkworms, then transit into the anterior of the silk gland, where the silk fibers are produced, stored and processed. The mechanism of formation and spinning of the silk fibers has not been fully elucidated, and transcriptome analyses specific to the anterior silk gland have not been reported. In the present study, we explored gene expression profiles in five regions of silk gland samples using the RNA-Seq method. As a result, there were 959,979,570 raw reads obtained, of which 583,068,172 reads were mapped to the silkworm genome. A total of 7419 genes were found to be expressed in terms of reads per kilobase of exon model per million mapped reads ≥ 5 in at least one sample. The gene numbers and expression levels of the expressed genes differed between these regions. The differentially expressed genes were analyzed, and 282 genes were detected as up-regulated in the anterior silk gland, compared with the other parts. Functions of these genes were addressed using the gene ontology and Kyoto Encyclopedia of Genes and Genomes databases, and seven key pathways were enriched. It suggested that the ion transportation, energy metabolism, protease inhibitors and cuticle proteins played essential roles in the process of silk formation and spinning in the anterior silk gland. In addition, 210 genes were found differently expressed between males and females, which should help to elucidate the mechanism of the quality difference in silk fibers from male and female silkworms.

The Asian Rice Gall Midge (Orseolia oryzae) Mitogenome Has Evolved Novel Gene Boundaries and Tandem Repeats That Distinguish Its Biotypes

The complete mitochondrial genome of the Asian rice gall midge, Orseolia oryzae (Diptera; Cecidomyiidae) was sequenced, annotated and analysed in the present study. The circular genome is 15,286 bp with 13 protein-coding genes, 22 tRNAs and 2 ribosomal RNA genes, and a 578 bp non-coding control region. All protein coding genes used conventional start codons and terminated with a complete stop codon. The genome presented many unusual features: (1) rearrangement in the order of tRNAs as well as protein coding genes; (2) truncation and unusual secondary structures of tRNAs; (3) presence of two different repeat elements in separate non-coding regions; (4) presence of one pseudo-tRNA gene; (5) inversion of the rRNA genes; (6) higher percentage of non-coding regions when compared with other insect mitogenomes. Rearrangements of the tRNAs and protein coding genes are explained on the basis of tandem duplication and random loss model and why intramitochondrial recombination is a better model for explaining rearrangements in the O. oryzae mitochondrial genome is discussed. Furthermore, we evaluated the number of iterations of the tandem repeat elements found in the mitogenome. This led to the identification of genetic markers capable of differentiating rice gall midge biotypes and the two Orseolia species investigated.

Comparative mitogenomic analysis of the superfamily Pentatomoidea (Insecta: Hemiptera: Heteroptera) and phylogenetic implications

BACKGROUND

Insect mitochondrial genomes (mitogenomes) are the most extensively used genetic marker for evolutionary and population genetics studies of insects. The Pentatomoidea superfamily is economically important and the largest superfamily within Pentatomomorpha with over 7,000 species. To better understand the diversity and evolution of pentatomoid species, we sequenced and annotated the mitogenomes of Eurydema gebleri and Rubiconia intermedia, and present the first comparative analysis of the 11 pentatomoid mitogenomes that have been sequenced to date.

RESULTS

We obtained the complete mitogenome of Eurydema gebleri (16,005 bp) and a nearly complete mitogenome of Rubiconia intermedia (14,967 bp). Our results show that gene content, gene arrangement, base composition, codon usage, and mitochondrial transcription termination factor sequences are highly conserved in pentatomoid species, especially for species in the same family. Evolutionary rate analyses of protein-coding genes reveal that the highest and lowest rates are found in atp8 and cox1 and distinctive evolutionary patterns are significantly correlated with the G + C content of genes. We inferred the secondary structures for two rRNA genes for eleven pentatomoid species, and identify some conserved motifs of RNA structures in Pentatomidea. All tRNA genes in pentatomoid mitogenomes have a canonical cloverleaf secondary structure, except for two tRNAs (trnS1 and trnV) which appear to lack the dihydrouridine arm. Regions that are A + T-rich have several distinct characteristics (e.g. size variation and abundant tandem repeats), and have potential as species or population level molecular markers. Phylogenetic analyses based on mitogenomic data strongly support the monophyly of Pentatomoidea, and the estimated phylogenetic relationships are: (Urostylididae + (Plataspidae + (Pentatomidae + (Cydnidae + (Dinidoridae + Tessaratomidae))))).

CONCLUSIONS

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the superfamily Pentatomoidea. Mitogenomes can be effectively used to resolve phylogenetic relationships of pentatomomorphan insects at various taxonomic levels. Sequencing more mitogenomes at various taxonomic levels, particularly from closely related species, will improve the annotation accuracy of mitochondrial genes, as well as greatly enhance our understanding of mitogenomic evolution and phylogenetic relationships in pentatomoids.

Social evolution. Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

Transgenerational inheritance of diet-induced genome rearrangements in Drosophila

Ribosomal RNA gene (rDNA) copy number variation modulates heterochromatin formation and influences the expression of a large fraction of the Drosophila genome. This discovery, along with the link between rDNA, aging, and disease, high-lights the importance of understanding how natural rDNA copy number variation arises. Pursuing the relationship between rDNA expression and stability, we have discovered that increased dietary yeast concentration, emulating periods of dietary excess during life, results in somatic rDNA instability and copy number reduction. Modulation of Insulin/TOR signaling produces similar results, indicating a role for known nutrient sensing signaling pathways in this process. Furthermore, adults fed elevated dietary yeast concentrations produce offspring with fewer rDNA copies demonstrating that these effects also occur in the germline, and are transgenerationally heritable. This finding explains one source of natural rDNA copy number variation revealing a clear long-term consequence of diet.

Improved statistical methods enable greater sensitivity in rhythm detection for genome-wide data

Robust methods for identifying patterns of expression in genome-wide data are important for generating hypotheses regarding gene function. To this end, several analytic methods have been developed for detecting periodic patterns. We improve one such method, JTK_CYCLE, by explicitly calculating the null distribution such that it accounts for multiple hypothesis testing and by including non-sinusoidal reference waveforms. We term this method empirical JTK_CYCLE with asymmetry search, and we compare its performance to JTK_CYCLE with Bonferroni and Benjamini-Hochberg multiple hypothesis testing correction, as well as to five other methods: cyclohedron test, address reduction, stable persistence, ANOVA, and F24. We find that ANOVA, F24, and JTK_CYCLE consistently outperform the other three methods when data are limited and noisy; empirical JTK_CYCLE with asymmetry search gives the greatest sensitivity while controlling for the false discovery rate. Our analysis also provides insight into experimental design and we find that, for a fixed number of samples, better sensitivity and specificity are achieved with higher numbers of replicates than with higher sampling density. Application of the methods to detecting circadian rhythms in a metadataset of microarrays that quantify time-dependent gene expression in whole heads of Drosophila melanogaster reveals annotations that are enriched among genes with highly asymmetric waveforms. These include a wide range of oxidation reduction and metabolic genes, as well as genes with transcripts that have multiple splice forms.

Culicoides Latreille (Diptera: Ceratopogonidae) taxonomy: current challenges and future directions

Culicoides Latreille biting midges (Diptera: Ceratopogonidae) cause a significant biting nuisance to humans, livestock and equines, and are the biological vectors of a range of internationally important pathogens of both veterinary and medical importance. Despite their economic significance, the delimitation and identification of species and evolutionary relationships between species within this genus remains at best problematic. To date no phylogenetic study has attempted to validate the subgeneric classification of the genus and the monophyly of many of the subgenera remains doubtful. Many informal species groupings are also known to exist but few are adequately described, further complicating accurate identification. Recent contributions to Culicoides taxonomy at the species level have revealed a high correlation between morphological and molecular analyses although molecular analyses are revealing the existence of cryptic species. This review considers the methods for studying the systematics of Culicoides using both morphological and genetic techniques, with a view to understanding the factors limiting our current understanding of Culicoides biology and hence arbovirus epidemiology. In addition, we examine the global status of Culicoides identification, highlighting areas that are poorly addressed, including the potential implementation of emerging technologies.

Environmentally selected aphid variants in clonality context display differential patterns of methylation in the genome

Heritability of acquired phenotypic traits is an adaptive evolutionary process that appears more complex than the basic allele selection guided by environmental pressure. In insects, the trans-generational transmission of epigenetic marks in clonal and/or sexual species is poorly documented. Aphids were used as a model to explore this feature because their asexual phase generates a stochastic and/or environment-oriented repertoire of variants. The a priori unchanged genome in clonal individuals prompts us to hypothesize whether covalent methyl DNA marks might be associated to the phenotypic variability and fitness selection. The full differential transcriptome between two environmentally selected clonal variants that originated from the same founder mother was mapped on the entire genomic scaffolds, in parallel with the methyl cytosine distribution. Data suggest that the assortments of heavily methylated DNA sites are distinct in these two clonal phenotypes. This might constitute an epigenetic mechanism that confers the robust adaptation of insect species to various environments involving clonal reproduction.