Transcriptome sequencing and phylogenetic analysis of four species of luminescent beetles

Eyeless cave-dwelling Leptonetela spiders still rely on light

Subterranean animals living in perpetual darkness may maintain photoresponse. However, the evolutionary processes behind the conflict between eye loss and maintenance of the photoresponse remain largely unknown. We used Leptonetela spiders to investigate the driving forces behind the maintenance of the photoresponse in cave-dwelling spiders. Our behavioral experiments showed that all eyeless/reduced-eyed cave-dwelling species retained photophobic response and that they had substantially decreased survival at cave entrances due to weak drought resistance. The transcriptomic analysis demonstrated that nearly all phototransduction pathway genes were present and that all tested phototransduction pathway genes were subjected to strong functional constraints in cave-dwelling species. Our results suggest that cave-dwelling eyeless spiders still use light and that light detection likely plays a role in avoiding the cave entrance habitat. This study confirms that some eyeless subterranean animals have retained their photosensitivity due to natural selection and provides a case of mismatch between phenotype and genotype or physiological function in a long-term evolutionary process.

Rhagophthalmidae Olivier, 1907 (Coleoptera, Elateroidea): described genera and species, current problems, and prospects for the bioluminescent and paedomorphic beetle lineage

Rhagophthalmidae are a small beetle family known from the eastern Palaearctic and Oriental realms. Rhagophthalmidae are closely related to railroad worms (Phengodidae) and fireflies (Lampyridae) with which they share highly modified paedomorphic females and the ability to emit light. Currently, Rhagophthalmidae include 66 species classified in the following 12 genera: Bicladodrilus Pic, 1921 (two spp.), Bicladum Pic, 1921 (two spp.), Dioptoma Pascoe, 1860 (two spp.), Diplocladon Gorham, 1883 (two spp.), Dodecatoma Westwood, 1849 (eight spp.), Falsophrixothrix Pic, 1937 (six spp.), Haplocladon Gorham, 1883 (two spp.), Menghuoius Kawashima, 2000 (three spp.), Mimoochotyra Pic, 1937 (one sp.), Monodrilus Pic, 1921 (two spp. in two subgenera), Pseudothilmanus Pic, 1918 (two spp.), and Rhagophthalmus Motschulsky, 1854 (34 spp.). The replacement name Haplocladongorhami Kundrata, nom. nov. is proposed for Diplocladonhasseltii Gorham, 1883b (described in subgenus Haplocladon) which is preoccupied by Diplocladonhasseltii Gorham, 1883a. The genus Reductodrilus Pic, 1943 is tentatively placed in Lampyridae: Ototretinae. Lectotypes are designated for Pseudothilmanusalatus Pic, 1918 and P.marginalis Pic, 1918. Interestingly, in the eastern part of their distribution, Rhagophthalmidae have remained within the boundaries of the Sunda Shelf and the Philippines demarcated by the Wallace Line, which separates the Oriental and Australasian realms. This study is intended to be a first step towards a comprehensive revision of the group on both genus and species levels. Additionally, critical problems and prospects for rhagophthalmid research are briefly discussed.

Morphological and molecular characterization of variation in common bean (Phaseolus vulgaris L.) germplasm from Azad Jammu and Kashmir, Pakistan

Phaseolus vulgaris, an essential food and source of protein, is cultivated across the world. This study was carried out to investigate the diversity and population structure of 34 P. vulgaris landrace accessions collected from the Azad Jammu and Kashmir (AJ&K) regions of Pakistan. The samples were analyzed both morphologically and using genetic variation identified through RNA sequencing. Our results indicated that most genetic variation occurs among local accessions, with little genetic variation occurring between geographical regions. In addition, the accessions fell into two major genetic groups. Morphological analysis revealed that these two genetic groups differ in a number of quantitative traits, including seed length, seed width, and seed weight. One accession, DUD-11, appears to be a mixture of the two major groups genetically as well as morphologically. Among the other accessions, DUD-8, RWK-2, and NGD-1 depicted particularly high seed weight along with higher seed length, seed width, and seed yield per plant. We suggest focusing on these accessions in future breeding programs. More generally, our results provide baseline data that will be useful for crop improvement and effective cultivation practices in Pakistan.

The complete mitochondrial genome of the firefly Curtos fulvocapitalis (Coleoptera: Lampyridae)

We report the complete mitochondrial genome of the firefly Curtos fulvocapitalis Jeng et Yang 1998. The circular genome is 16,398 bp and has a base composition of A (42.21%), C (11.22%), G (7.73%), and T (38.85%). Our sequence is similar to those of other metazoans, which contain 13 protein-coding genes. All 13 protein-coding genes were initiated by the ATN (ATT, ATA, and ATG) codon. Nine protein-coding genes stopped with TAA or TAG codons, and the other four genes had an incomplete termination codon, a single T. We sequenced the mitochondrial genome of fireflies to analyze phylogenetic relationships and determine the evolution of their flashing signals.

Characterization of the First Complete Mitochondrial Genome of Cyphonocerinae (Coleoptera: Lampyridae) with Implications for Phylogeny and Evolution of Fireflies

Complete mitochondrial genomes are valuable resources for phylogenetics in insects. The Cyphonoceridae represents an important lineage of fireflies. However, no complete mitogenome is available until now. Here, the first complete mitochondrial genome from this subfamily was reported, with Cyphonocerus sanguineus klapperichi as a representative. The mitogenome of C. sanguineus klapperichi was conserved in the structure and comparable to that of others in size and A+T content. Nucleotide composition was A+T-biased, and all genes exhibited a positive AT-skew and negative GC-skew. Two types of tandem repeat sequence units were present in the control region (136 bp × 2; 171 bp × 2 + 9 bp). For reconstruction of Lampyridae's phylogeny, three different datasets were analyzed by both maximum likelihood (ML) and Bayesian inference (BI) methods. As a result, the same topology was produced by both ML analysis of 13 protein-coding genes and 2rRNA and BI analysis of 37 genes. The results indicated that Lampyridae, Lampyrinae, Luciolinae (excluding Emeia) were monophyletic, but Ototretinae was paraphyletic, of which Stenocladius was recovered as the sister taxon to all others, while Drilaster was more closely related to Cyphonocerinae; Phturinae + Emeia were included in a monophyletic clade, which comprised sister groups with Lampyridae. Vesta was deeply rooted in the Luciolinae.

Phylogenetic analysis of eight species of Anomopoda based on transcriptomic and mitochondrial DNA sequences

Anomopoda is the widespread planktonic microcrustacean, which plays a crucial role in aquatic ecosystem. There are few studies about the evolutionary relationships among various Anomopoda basing on molecular data. In the present study, phylogenetic analysis of eight Anomopoda was carried out. Firstly, the culture system was developed to breed cladocerans. By using this system, eight species (Daphnia magna, D. pulex, D. sinensis, Ceriodaphnia reticulata, Moina micrura, Scapholeberis kingi, Simocephalus vetulus and Eurycercus lamellatus) were purified and cultured stably in the laboratory. Then, transcriptomic sequences and partial mitochondrial DNA sequences were both used to reconstruct the phylogenetic tree among 8 species. Transcriptomic sequences were sequenced on Illumina Hiseq 2500 platform. After assembly and annotation, transcriptomic sequences were spliced together and aligned for phylogenetic analysis. Basing on the orthologous genes derived from transcriptomic sequences, the phylogenetic analysis showed that 4 genera of Daphniidae were clustered into one group, and among the 4 genera, Ceriodaphnia was closer to Daphnia than Simocephalus, while Scapholeberis was farthest from other species. In addition, Eurycercidae was closer to Daphniidae than Moinidae. The phylogenetic trees based on both 12S rRNA and 16S rRNA sequences were similar with that based on transcriptomic sequences. Meanwhile, the phylogenetic tree based on 16S rRNA sequences was more suitable than that based on 12S rRNA sequences. These results suggested that the phylogenetic analysis basing on the transcriptomic sequences was available in cladocerans, which will help us to effectively understand the phylogenetic relationships among various cladocerans.

Inferring the Phylogenetic Positions of Two Fig Wasp Subfamilies of Epichrysomallinae and Sycophaginae Using Transcriptomes and Mitochondrial Data

Fig wasps are a group of insects (Hymenoptera: Chalcidoidea) that live in the compact syconia of fig trees (Moraceae: Ficus). Accurate classification and phylogenetic results are very important for studies of fig wasps, but the taxonomic statuses of some fig wasps, especially the non-pollinating subfamilies are difficult to determine, such as Epichrysomallinae and Sycophaginae. To resolve the taxonomic statuses of Epichrysomallinae and Sycophaginae, we obtained transcriptomes and mitochondrial genome (mitogenome) data for four species of fig wasps. These newly added data were combined with the data of 13 wasps (data on 11 fig wasp species were from our laboratory and two wasp species were download from NCBI). Based on the transcriptome and genome data, we obtained 145 single-copy orthologous (SCO) genes in 17 wasp species, and based on mitogenome data, we obtained 13 mitochondrial protein-coding genes (PCGs) for each of the 17 wasp species. Ultimately, we used 145 SCO genes, 13 mitochondrial PCGs and combined SCO genes and mitochondrial genes data to reconstruct the phylogenies of fig wasps using both maximum likelihood (ML) and Bayesian inference (BI) analyses. Our results suggest that both Epichrysomallinae and Sycophaginae are more closely related to Agaonidae with a high statistical support.

The genetic variations in the mitochondrial genomes of three Luciolinae fireflies

This paper studied the complete mitochondrial genomes of three fireflies, Pygoluciola qingyu, Emeia pseudosauteri and Abscondita terminalis. We discussed the variations in the mitochondrial genomes of samples of each firefly from different populations. The mitochondrial genomes of Abs. terminalis and P. qingyu are very stable among their different populations. The mitochondrial genome of E. pseudosauteri shows some variations among the different populations, especially in the COI sequence.

Complete mitochondrial genome of the Korean endemic firefly, Luciola unmunsana (Coleoptera: Lampyridae)

In this study, we announce the complete mitochondrial genome (mitogenome) of the Korean endemic firefly, Luciola unmunsana Doi, 1931. The full-length circular genome was 15,858 bp, with 77.94% A/T content. It contained the typical set of 37 metazoan genes: 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and 2 ribosomal RNA (rRNA) genes, as well as an A + T-rich region. The gene arrangement of the species is identical to that of the ancestral arrangement found in the majority of insects. The maximum-likelihood tree, built using all PCGs and two rRNAs via randomized accelerated maximum likelihood (RAxML) showed that L. unmunsana was grouped as a sister to L. curtithorax with the highest nodal support. However, another Luciola species clustered with the Aquatica species such that the genus Luciola was a non-monophyletic group. Therefore, more sampling is required to clarify the phylogeny of Luciola.

Identification of potential candidate genes involved in the sex determination cascade in an aquatic firefly, Sclerotia aquatilis (Coleoptera, Lampyridae)

Sexual differentiation, dimorphism, and courtship behavior are the downstream developmental programs of the sex determination cascade. The sex determination cascade in arthropods often involves key genes, transformer (tra), doublesex (dsx), transformer-2 (tra2), and fruitless (fru). These genes are conserved among insect taxa; however, they have never been reported in fireflies. In this study, the candidate genes for these key genes were identified for the first time in an aquatic firefly, Sclerotia aquatilis using transcriptome analysis. A comparative protein-protein interaction (PPI) network of sex determination cascade was reconstructed for S. aquatilis based on a network of a model insect, Drosophila melanogaster. Subsequently, a sex determination cascade in S. aquatilis was proposed based on the amino acid sequence structures and expression profiles of these candidates. This study describes the first efforts toward understanding the molecular control of sex determination cascade in fireflies.

Trehalase Gene as a Molecular Signature of Dietary Diversification in Mammals

Diet is a key factor in determining and structuring animal diversity and adaptive radiations. The mammalian fossil record preserves phenotypic evidence of many dietary shifts, whereas genetic changes followed by dietary diversification in mammals remain largely unknown. To test whether living mammals preserve molecular evidence of dietary shifts, we examined the trehalase gene (Treh), which encodes an enzyme capable of digesting trehalose from insect blood, in bats and other mammals with diverse diets. Bats represent the largest dietary radiation among all mammalian orders, with independent origins of frugivory, nectarivory, carnivory, omnivory, and even sanguivory in an otherwise insectivorous clade. We found that Treh has been inactivated in unrelated bat lineages that independently radiated into noninsectivorous niches. Consistently, purifying selection has been markedly relaxed in noninsectivorous bats compared with their insectivorous relatives. Enzymatic assays of intestinal trehalase in bats suggest that trehalase activity tends to be lost or markedly reduced in noninsectivorous bats compared with their insectivorous relatives. Furthermore, our survey of Treh in 119 mammal species, which represent a deeper evolutionary timeframe, additionally identified a number of other independent losses of Treh in noninsectivorous species, recapitulating the evolutionary pattern that we found in bats. These results document a molecular record of dietary diversification in mammals, and suggest that such molecular signatures of dietary shifts would help us understand both historical and modern changes of animal diets.

Convergent reduction of V1R genes in subterranean rodents

BACKGROUND

Vomeronasal type 1 receptor genes (V1Rs) are expected to detect intraspecific pheromones. It is believed that rodents rely heavily on pheromonal communication mediated by V1Rs, but pheromonal signals are thought to be confined in subterranean rodents that live in underground burrows. Thus, subterranean rodents may show a contrasting mode of V1R evolution compared with their superterranean relatives.

RESULTS

We examined the V1R evolution in subterranean rodents by analyzing currently available genomes of 24 rodents, including 19 superterranean and 5 subterranean species from three independent lineages. We identified a lower number of putatively functional V1R genes in each subterranean rodent (a range of 22-40) compared with superterranean species (a range of 63-221). After correcting phylogenetic inertia, the positive correlation remains significant between the small V1R repertoire size and the subterranean lifestyle. To test whether V1Rs have been relaxed from functional constraints in subterranean rodents, we sequenced 22 intact V1Rs in 29 individuals of one subterranean rodent (Spalax galili) from two soil populations, which have been proposed to undergo incipient speciation. We found 12 of the 22 V1Rs to show significant genetic differentiations between the two natural populations, indicative of diversifying selection.

CONCLUSION

Our study demonstrates convergent reduction of V1Rs in subterranean rodents from three independent lineages. Meanwhile, it is noteworthy that most V1Rs in the two Spalax populations are under diversifying selection rather than relaxed selection, suggesting that functional constraints on these genes may have retained in some subterranean species.

Phylogenetic analysis provides insights into the evolution of Asian fireflies and adult bioluminescence

Fireflies are one of the best-known examples of luminescent organisms. The limited geographic distribution and rarity of some firefly genera have hindered molecular phylogenetic analysis, resulting in uncertainty in regard to firefly phylogeny. Here, using genome skimming next-generation sequencing, we sequenced 23 Asian firefly species from 15 genera (Lampyridae: 14; Rhagophthalmidae: one) and assembled their mitochondrial genomes (mitogenomes) and nuclear ribosomal DNA (rDNA) repeat unit. The mitogenomes (including 15 mitochondrial genes: COX1-3, ATP6&8, ND1-6&4L, CYTB, 12S, and 16S) were recovered for almost all 23 species; furthermore, three regions of the nuclear rDNA repeat unit (18S, 28S, and 5.8S) were recovered for 22 out of the 23 species. The mitogenomes of 11 genera and 22 species as well as the complete rDNA from 22 species are reported here for the first time. Combined with previously published sequences of mitochondrial and rDNA coding regions, 166 species (170 populations with four overlapping in Lampyridae) were included in the current analyses. We selected different species groups and coding regions to infer phylogenies, and then employed tree certainty (TC) and internode certainty (IC) to quantify any phylogenetic incongruence. Phylogenetic analysis of 18 coding regions (15 mitochondrial genes and three regions of the nuclear rDNA repeat unit) from different species groups showed that the 144-species selection group (excluding 22 species outside Lampyridae) had relatively high TC (101.39). Further phylogenetic analysis of the 144 species using different coding regions indicated that the phylogeny of the 13 coding regions (10 mitochondrial genes: COX1-2, ATP6&8, ND1, ND4-5, CYTB, 12S and 16S; three rDNA regions: 18S, 5.8S, and 28S) demonstrated higher TC (103.02) than the phylogenies based on the 18 coding regions (TC = 101.39), conserved-regions (c-regions, i.e., 12S, 16S, COX1, 18S, and 28S) (TC = 95.11), or conserved-sites (c-sites, TC = 92.31) for the mitochondrial genes. In contrast, the c-sites strengthened the deeper nodes of the 144-species phylogeny compared to the c-regions. All of the 144-species phylogenies using different coding regions (except the c-regions) consistently recovered the monophyly of each of the three luminous families and their combination (Lampyridae, Rhagophthalmidae, and Phengodidae) with high IC support. Our phylogenetic analyses clarified the position of firefly genera Lamprigera, Vesta, Stenocladius, Pyrocoelia, Diaphanes, Abscondita, Pygoluciola, Emeia, Pristolycus, and Menghuoius. We also inferred the evolutionary pattern of adult bioluminescence in Lampyridae based on the phylogenies of 166 and 144 species. Our data suggest that the common ancestor of Lampyridae possessed adult bioluminescence, with a higher loss rate than gain rate of bioluminescence during its lineage evolution. Our results provide insight into Asian firefly phylogeny, and also enrich mitogenome and rDNA data resources for further study.

Reconstruction of insect hormone pathways in an aquatic firefly, Sclerotia aquatilis (Coleoptera: Lampyridae), using RNA-seq

Insect hormones: ecdysteroids and juvenile hormones have crucial functions during the regulation of different developmental pathways in insects. Insect metamorphosis is one of the primary pathways regulated by these hormones. The insect hormone biosynthetic pathway is conserved among arthropods, including insects, with some variations in the form of hormones used among each group of insects. In this study, the candidate genes involved in the insect hormone pathways and their functional roles were assessed in an aquatic firefly, Sclerotia aquatilis using a high-throughput RNA sequencing technique. Illumina next-generation sequencing (NGS) was used to generate transcriptome data for the different developmental stages (i.e., larva, pupa, and adult) of S. aquatilis. A total of 82,022 unigenes were generated across all different developmental stages. Functional annotation was performed for each gene, based on multiple biological databases, generating 46,230 unigenes. These unigenes were subsequently mapped using KEGG pathways. Accordingly, 221 protein-encoding genes involved in the insect hormone pathways were identified, including, JHAMT, CYP15A1, JHE, and Halloween family genes. Twenty potential gene candidates associated with the biosynthetic and degradation pathways for insect hormones were subjected to real-time PCR, reverse transcriptase PCR (RT-PCR) and sequencing analyses. The real-time PCR results showed similar expression patterns as those observed for transcriptome expression profiles for most of the examined genes. RT-PCR and Sanger sequencing confirmed the expressed coding sequences of these gene candidates. This study is the first to examine firefly insect hormone pathways, facilitating a better understanding of firefly growth and development.

Phylogenomic analyses and divergence time estimation of Elateroidea (Coleoptera) based on RNA-Seq data

Bioluminescence, the emission of visible light in a living organism, is an intriguing phenomenon observed in different species and environments. In terrestrial organisms, the bioluminescence is observed mainly in beetles of the Elateroidea superfamily (Coleoptera). Several phylogenetic studies have been used different strategies to propose a scenario for the origin and evolution of the bioluminescence within this group, however some of them showed incongruences, mainly about the relationship of the bioluminescent families. In order to increase the number of molecular markers available for Elateroidea species and to propose a more accurate phylogeny, with high supported topology, we employed the Next-Generation Sequencing (NGS) methodology to perform the RNA-Seq analysis of luminescent (Elateridae, Phengodidae, Rhagophthalmidae, and Lampyridae) and non-luminescent (Cantharidae) species of Neotropical beetles. We used the RNA-Seq data to construct a calibrated phylogeny of Elateroidea superfamily using a large number of nuclear molecular markers. The results indicate Lampyridae and Phengodidae/Rhagophthalmidae as sister-groups, suggesting that the bioluminescence evolved later in Elateridae than other families (Lampyridae, Phengodidae, and Rhagophthalmidae), and indicating the Upper Cretaceous as the period for the main diversification of Elateroidea bioluminescent species.

The transcriptome analysis of Protaetia brevitarsis Lewis larvae

Larvae of the pest Protaetia brevitarsis are used to treat infections in traditional Chinese medicine. However, genomic information about this non-model species is currently lacking. To better understand the fundamental biology of this non-model species, its transcriptome was obtained using next generation sequencing and then analyzed. A total of 7.62 Gb of clean reads were obtained, which were assembled into 169,087 transcripts corresponding to 142,000 annotated unigenes. These unigenes were functionally classified according to Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. A total of 41,921 unigenes were assigned to 56 GO terms, 21,454 unigenes were divided among 26 KOG categories, and 16,368 unigenes were assigned to 32 KEGG pathways. In addition, 19,144 simple sequence repeats (SSRs) were identified. Furthermore, several kinds of natural antimicrobial peptides and proteins, 4 histones with potential antimicrobial activity, and 41 potential antimicrobial peptide sequences were identified. These data are the first reported whole transcriptome sequence of P. brevitarsis larvae, which represents a valuable genomic resource for studying this species, thus promoting the utilization of its medical potential.

Genome sequences identify three families of Coleoptera as morphologically derived click beetles (Elateridae)

Plastoceridae Crowson, 1972, Drilidae Blanchard, 1845 and Omalisidae Lacordaire, 1857 (Elateroidea) are families of the Coleoptera with obscure phylogenetic relationships and modified morphology showing neotenic traits such as soft bodies, reduced wing cases and larviform females. We shotgun sequenced genomes of Plastocerus, Drilus and Omalisus and incorporated them into data matrices of 66 and 4202 single-copy nuclear genes representing Elateroidea. Phylogenetic analyses indicate their terminal positions within the broadly defined well-sclerotized and fully metamorphosed Elateridae and thus Omalisidae should now be considered as Omalisinae stat. nov. in Elateridae Leach, 1815. The results support multiple independent origins of incomplete metamorphosis in Elateridae and indicate the parallel evolution of morphological and ecological traits. Unlike other neotenic elateroids derived from the supposedly pre-adapted aposematically coloured and unpalatable soft-bodied elateroids, such as fireflies (Lampyridae) and net-winged beetles (Lycidae), omalisids and drilids evolved from well-sclerotized click beetles. These findings suggest sudden morphological shifts through incomplete metamorphosis, with important implications for macroevolution, including reduced speciation rate and high extinction risk in unstable habitats. Precise phylogenetic placement is necessary for studies of the molecular mechanisms of ontogenetic shifts leading to profoundly changed morphology.

Long-read sequence assembly of the firefly Pyrocoelia pectoralis genome

Background

Fireflies are a family of insects within the beetle order Coleoptera, or winged beetles, and they are one of the most well-known and loved insect species because of their bioluminescence. However, the firefly is in danger of extinction because of the massive destruction of its living environment. In order to improve the understanding of fireflies and protect them effectively, we sequenced the whole genome of the terrestrial firefly Pyrocoelia pectoralis.

Findings

Here, we developed a highly reliable genome resource for the terrestrial firefly Pyrocoelia pectoralis (E. Oliv., 1883; Coleoptera: Lampyridae) using single molecule real time (SMRT) sequencing on the PacBio Sequel platform. In total, 57.8 Gb of long reads were generated and assembled into a 760.4-Mb genome, which is close to the estimated genome size and covered 98.7% complete and 0.7% partial insect Benchmarking Universal Single-Copy Orthologs. The k-mer analysis showed that this genome is highly heterozygous. However, our long-read assembly demonstrates continuousness with a contig N50 length of 3.04 Mb and the longest contig length of 13.69 Mb. Furthermore, 135 589 SSRs and 341 Mb of repeat sequences were detected. A total of 23 092 genes were predicted; 88.44% of genes were annotated with one or more related functions.

Conclusions

We assembled a high-quality firefly genome, which will not only provide insights into the conservation and biodiversity of fireflies, but also provide a wealth of information to study the mechanisms of their sexual communication, bio-luminescence, and evolution.

In silico identification and characterization of a diverse subset of conserved microRNAs in bioenergy crop Arundo donax L

MicroRNAs (miRNAs) are small non-coding RNA molecules involved in the post-transcriptional regulation of gene expression in plants. Arundo donax L. is a perennial C₃ grass considered one of the most promising bioenergy crops. Despite its relevance, many fundamental aspects of its biology still remain to be elucidated. In the present study we carried out the first in silico mining and tissue-specific characterization of microRNAs and their putative targets in A. donax. We identified a total of 141 miRNAs belonging to 14 families along with the corresponding primary miRNAs, precursor miRNAs and a total of 462 high-confidence predicted targets and novel target sites were validated by 5′-race. Gene Ontology functional annotation showed that miRNA targets are constituted mainly by transcription factors, but three of the newly validated targets are enzymes involved in novel functions like RNA editing, acyl lipid metabolism and post-Golgi trafficking. Folding variability of pre-miRNA loops and phylogenetic analyses indicate variable selective pressure acting on the different miRNA families. The set of miRNAs identified in this study will pave the road to further miRNA research in Arundo donax and contribute towards a better understanding of miRNA-mediated gene regulatory processes in other bioenergy crops.

Genome sequencing of Rhinorhipus Lawrence exposes an early branch of the Coleoptera

Background

Rhinorhipidae Lawrence, 1988 is an enigmatic beetle family represented by a single species, Rhinorhipus tamborinensis Lawrence, 1988, from Australia, with poorly established affinities near the superfamily Elateroidea (click beetles, soldier beetles and fireflies) or the more inclusive series (infraorder) Elateriformia. Its evolutionary position may inform the basal relationships of the suborder Polyphaga, the largest clade of Coleoptera.

Results

We analyzed four densely sampled DNA datasets of major coleopteran lineages for mitogenomes, rRNA genes and single copy nuclear genes. Additionally, genome sequencing was used for incorporation of R. tamborinensis into a set of 4220 orthologs for 24 terminals representing 12 polyphagan superfamilies. Topologies differed to various degrees, but all consistently refute the proposed placement of Rhinorhipidae in Elateroidea and instead indicate either sister relationships with other Elateriformia, frequently together with Nosodendridae, another divergent small family hitherto placed in Derodontoidea, or in an isolated position among the deepest lineages of Polyphaga. The phylogenomic analyses recovered Rhinorhipus in a sister position to all other Elateriformia composed of five superfamilies. Therefore, we erect the new superfamily Rhinorhipoidea Lawrence, 1988, stat. Nov., with the type-family Rhinorhipidae. The origins of the Rhinorhipidae were dated to the Upper Triassic/Lower Jurassic at the very early phase of polyphagan diversification.

Conclusions

Thus, Rhinorhipidae adds another example to several recently recognized ancient relict lineages which are interspersed within contemporaneous hugely species-rich lineages of Coleoptera.

Electronic supplementary material

The online version of this article (10.1186/s12983-018-0262-0) contains supplementary material, which is available to authorized users.

Transcriptome sequencing and differential gene expression analysis of the schistosome-transmitting snail Oncomelania hupensis inhabiting hilly and marshland regions

The freshwater snail Oncomelania hupensis is the unique intermediate host of the blood fluke Schistosoma japonicum, which is the major cause of schistosomiasis. The snail inhabits two contrasting environments: the hilly and marshland regions. The hilly snails are smaller in size and have the typical smooth shell, whereas the marshland snails are larger and possess the ribbed shell. To reveal the differences in gene expression between the hilly and marshland snails, a total of six snails, three per environment, were individually examined by RNA sequencing technology. All paired-end reads were assembled into contigs from which 34,760 unigenes were predicted. Based on single nucleotide polymorphisms, principal component analysis and neighbor-joining clustering revealed two distinct clusters of hilly and marshland snails. Analysis of expression changes between environments showed that upregulated genes relating to immunity and development were enriched in hilly snails, while those associated with reproduction were over-represented in marshland snails. Eight differentially expressed genes between the two types of snails were validated by qRT-PCR. Our study identified candidate genes that could be targets for future functional studies, and provided a link between expression profiling and ecological adaptation of the snail that may have implications for schistosomiasis control.

The complete mitochondrial genome sequence and phylogenetic analysis of Luciola lateralis, one of the most famous firefly in Japan (Coleoptera: Lampyridae)

We will report the complete mitochondrial genome sequence of Japanese firefly ‘Heike Botaru’, Luciola lateralis (Coleoptera: Lampyridae). Total length of this mitogenome was 16,719 bp and the composition of each base was A (42.50%), C (9.01%), G (14.16%), T (34.33%), respectively. The obtained sequence fulfils general mitogenome composition of metazoan (13 protein coding sequences (CDSs), 22 tRNA genes, two rRNA subunits, and an AT-rich region). From the phylogenetic tree analysis using 25 kinds of insect mitogenome including firefly family was found that L. lateralis is the closest to the genus Aquatica.