Diversification in tropics and subtropics following the mid-Miocene climate change: A case study of the spider genus Nesticella

A targeted gene phylogenetic framework to investigate diversification in the highly diverse yet geographically restricted red devil spiders (Araneae, Dysderidae)

The family Dysderidae is a highly diverse group of nocturnal ground-dwelling and active-hunter spiders. Dysderids are mostly restricted to the Western Palearctic, and particularly rich and abundant around the Mediterranean region. Interestingly, the distribution of species richness among its 24 genera and three subfamilies is highly biased-80% of its 644 documented species belong to just two genera, Dysdera (326) and Harpactea (211). Dysderidae provides an excellent study case for evolutionary and ecological research. It includes cases of trophic specialization, which are uncommon among spiders, and exhibit other remarkable biological (e.g. holocentric chromosomes), behavioural (e.g. cryptic female choice), evolutionary (e.g. adaptive radiation) and ecological features (e.g. recurrent colonization of the subterranean environment). The lack of a quantitative hypothesis on its phylogenetic structure has hampered its potential as a testing ground for evolutionary, biogeographical and ecological hypotheses. Here, we present the results of a target, multi-locus phylogenetic analysis, using mitochondrial (cox1, 16s and 12s) and nuclear genes (h3, 28s and 18s), of the most exhaustive taxonomic sample within Dysderidae (104 spp.) to date and across related families (Synspermiata) (83 spp.). We estimate divergence times using a combination of fossil and biogeographic node calibrations and use this timeline to identify shifts in diversification rates. Our results support the monophyly of the Dysderidae subfamilies Rhodinae and Dysderinae but reject Harpacteinae as currently defined. Moreover, the clades recovered within Harpacteinae do not support its current taxonomy. The origin of the family most likely post-dated the break-up of Pangea, and cave colonization may be older than previously considered. After correcting for the taxonomic artefacts, we identified a significant shift in diversification rates at the base of the genus Dysdera. Although the unique coexistence of specialist and generalist diets within the lineage could be suggested as the potential driver for the rate acceleration, further quantitative analyses would be necessary to test this hypothesis.

Environmental specificity of karst cave habitats evidenced by diverse symbiotic bacteria in Opiliones

BACKGROUND

Karst caves serve as natural laboratories, providing organisms with extreme and constant conditions that promote isolation, resulting in a genetic relationship and living environment that is significantly different from those outside the cave. However, research on cave creatures, especially Opiliones, remains scarce, with most studies focused on water, soil, and cave sediments.

RESULTS

The structure of symbiotic bacteria in different caves were compared, revealing significant differences. Based on the alpha and beta diversity, symbiotic bacteria abundance and diversity in the cave were similar, but the structure of symbiotic bacteria differed inside and outside the cave. Microorganisms in the cave play an important role in material cycling and energy flow, particularly in the nitrogen cycle. Although microbial diversity varies inside and outside the cave, Opiliones in Beijing caves and Hainan Island exhibited a strong similarity, indicating that the two environments share commonalities.

CONCLUSIONS

The karst cave environment possesses high microbial diversity and there are noticeable differences among different caves. Different habitats lead to significant differences in the symbiotic bacteria in Opiliones inside and outside the cave, and cave microorganisms have made efforts to adapt to extreme environments. The similarity in symbiotic bacteria community structure suggests a potential similarity in host environments, providing an explanation for the appearance of Sinonychia martensi in caves in the north.

Sustainable Development versus Extractivist Deforestation in Tropical, Subtropical, and Boreal Forest Ecosystems: Repercussions and Controversies about the Mother Tree and the Mycorrhizal Network Hypothesis

This research reviews the phenomenon of extractive deforestation as a possible trigger for cascade reactions that could affect part of the forest ecosystem and its biodiversity (surface, aerial, and underground) in tropical, subtropical, and boreal forests. The controversy and disparities in criteria generated in the international scientific community around the hypothesis of a possible link between "mother trees" and mycorrhizal networks in coopetition for nutrients, nitrogen, and carbon are analyzed. The objective is to promote awareness to generate more scientific knowledge about the eventual impacts of forest extraction. Public policies are emphasized as crucial mediators for balanced sustainable development. Currently, the effects of extractive deforestation on forest ecosystems are poorly understood, which requires caution and forest protection. Continued research to increase our knowledge in molecular biology is advocated to understand the adaptation of biological organisms to the new conditions of the ecosystem both in the face of extractive deforestation and reforestation. The environmental impacts of extractive deforestation, such as the loss of biodiversity, soil degradation, altered water cycles, and the contribution of climate change, remain largely unknown. Long-term and high-quality research is essential to ensure forest sustainability and the preservation of biodiversity for future generations.

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.

Integrative taxonomic revision of the genera Nesticella and Howaia in Japan with the description of five new species (Araneae, Nesticidae, Nesticellini)

The Japanese species of the genera Nesticella Lehtinen & Saaristo, 1980 and Howaia Lehtinen & Saaristo, 1980 (Araneae, Nesticidae, Nesticellini) are revised using an integrative taxonomic approach. Their morphology, phylogenetic position within the genera, assignment to species groups, and distribution in mainland Japan and the Ryukyu islands are discussed herein. A phylogenetic and species delimitation analysis was conducted to confirm the boundaries between the putative species that were initially discriminated based on external and genital morphology. As a result of the present study, six species of Nesticella and three species of Howaia are proposed for the fauna of Japan based on the combined results of morphology and molecular analyses. Male and female of the previously known species H.mogera (Yaginuma, 1972), N.brevipes (Yaginuma, 1970), and N.okinawaensis (Yaginuma, 1979) are redescribed and illustrated using type specimens or specimens collected at the type locality. Nesticellaterrestris (Yaginuma, 1970) is resurrected as a valid species and distinguished from the closely related N.brevipes based on morphological and molecular evidence. The male of N.terrestris is described for the first time. We herein describe three new species of Nesticella and two new species of Howaia from different islands in the Ryukyu Archipelago, southwest Japan: N.insulanasp. nov. (♂, ♀) from Yonaguni-jima Island, N.occultasp. nov. (♀) from Ishigaki-jima Island, N.silvicolasp. nov. (♂, ♀) from Yakushima Island, H.albasp. nov. (♂, ♀) from Miyako-jima Island and H.subterraneasp. nov. (♀) from Okinoerabu-jima Island. Nesticellaoccultasp. nov., H.albasp. nov. and H.subterraneasp. nov. dwell exclusively in caves and show extensive morphological adaptation to subterranean life.

Geological events and climate change drive diversification and speciation of mute cicadas in eastern continental Asia

The poor mobility of nymphs living underground, usually for many years and the weak flying ability of adults make cicadas unique for evolutionary biology and bio-geographical study. Cicadas of the genus Karenia are unusual in Cicadidae in lacking the timbals that produce sound. Population differentiation, genetic structure, dispersal and evolutionary history of the eastern Asian mute cicada Karenia caelatata were investigated based on morphological, acoustic and molecular data. The results reveal a high level of genetic differentiation in this species. Six independent clades with nearly unique sets of haplotypes corresponding to geographically isolated populations are recognized. Genetic and geographic distances are significantly correlated among lineages. The phenotypic differentiation is generally consistent with the high levels of genetic divergence across populations. Results of ecological niche modeling suggest that the potential distribution range of this mountain-habitat specialist during the Last Glacial Maximum was broader than its current range, indicating this species had benefited from the climate change during the early Pleistocene in southern China. Geological events such as orogeny in Southwest China and Pleistocene climate oscillations have driven the differentiation and divergence of this species, and basins, plains and rivers function as natural "barriers" to block the gene flow. Besides significant genetic divergence being found among clades, the populations occurring in the Wuyi Mountains and the Hengduan Mountains are significantly different in the calling song structure from other populations. This may have resulted from significant population differentiation and subsequent adaptation of related populations. We conclude that ecological differences in habitats, coupled with geographical isolation, have driven population divergence and allopatric speciation. This study provides a plausible example of incipient speciation in Cicadidae and improves understanding of population differentiation, acoustic signal diversification and phylogeographic relationships of this unusual cicada species. It informs future studies on population differentiation, speciation and phylogeography of other mountain-habitat insects in the East Asian continent.

New Insights into the Variation and Admixture of the Cave-Dwelling Spider Trogloneta yunnanensis in South China Karst

Subterranean karst caves can contain unexpected biodiversity, but few studies related to spider population genetics have been conducted in the karst area of Southern China. In this study, we investigated the population genetic structure of Trogloneta yunnanensis (Song & Zhu, 1994) based on 73 spider samples from six underground populations in South China Karst. Population genetic structure analysis showed a clear divergence (FST > 0.9 and Nm < 0.05) among populations according to mitochondrial genes. The phylogenetic gene tree constructed by BI and ML methods recovered six geographic clades. Divergence time estimation indicated that the divergence of these six populations can be traced back to the late Pleistocene. We supposed that the geographic isolation led to the extreme population structure. According to this study and previous studies about troglobites living in this region, the subterranean habitats of the Yunnan-Guizhou Plateau may contain many organisms with similar genetic structures. The subterranean biodiversity in the karst area of Southern China needs to be re-evaluated and protected.

A molecular phylogeny of the European nesticid spiders (Nesticidae, Araneae): Implications for their systematics and biogeography

Nesticidae is a small family of spiders with a worldwide distribution that includes 15 genera and 272 described species. Seven genera and 56 species are known from Europe, distributed from the Iberian Peninsula to the Caucasus and the Ural Mountains. Most of these European species are cave dwellers and many of them are troglobites. In this study we present the first molecular phylogeny of the family Nesticidae in Europe with a wide geographical sampling across the continent. In our analysis the European nesticid fauna is well represented, including six genera and 40 of the 56 currently accepted species including the type species of all sampled genera. We have included in the analysis representatives of the North American and Asian fauna to test the monophyly of the European species and the phylogenetic relationships of European lineages. Phylogenetic relationships were reconstructed using maximum likelihood and Bayesian inference. As part of our Bayesian analyses, we also dated the phylogeny using two approaches, one based only on fossil calibrations and one that included an additional biogeographical constraint. Our results show paraphyly of the European nesticids with respect to the Asian and North American taxa. We recover four main lineages within Europe. These four European lineages and all European genera have 100% bootstrap support and high posterior probability support in the BEAST2 analysis. The Typhlonesticus lineage is the earliest branching clade present in Europe and includes seven species, the five currently accepted species plus T. parvus from Bosnia and Herzegovina and T. silvestrii from western North America. The Eastern lineage includes the genus Aituaria and is the sister group of the Asian genera Nesticella and Wraios. The Domitius lineage is likely the sister group of the Central European lineage and spreads over the Iberian and Italian peninsulas. Finally, the Central European lineage includes three genera: Kryptonesticus, distributed from the karstic massifs of the Balkan Peninsula to Turkey, Nesticus with a single synanthropic species N. cellulanus and Carpathonesticus, exclusive to the Carpathian Mountains. With the exception of the genus Typhlonesticus, all European genera show an allopatric distribution (except for the two European synanthropic species). The results obtained in this study together with the revision of the original descriptions, redescriptions, and illustrations, lead us to propose 11 nomenclatural changes (new combinations) concerning the genera Typhlonesticus, Nesticus and Carpathonesticus.

A new species of Paracortina from a Vietnamese cave, with remarkable secondary sexual characters in males (Callipodida, Paracortinidae)

A new millipede species, Paracortina kyrangsp. nov., is described from a cave in Cao Bang Province, northern Vietnam. The new species is diagnosed by having an extraordinarily long projection on the head of males, reduced eyes, a gonocoxite with two processes, a long and slender gonotelopodite with two long, clavate prefemoroidal processes densely covered with long macrosetae apically, and with a distal, reverse, short spine on mesal side, and a rather sinuous distal part of the telopodite. This is the third species of the genus that is known from Vietnam. A brief comparison of some secondary sexual characters is made.

Biotic colonization of subtropical East Asian caves through time

Caves are home to unique and fragile biotas with high levels of endemism. However, little is known about how the biotic colonization of caves has developed over time, especially in caves from middle and low latitudes. Subtropical East Asia holds the world's largest karst landform with numerous ancient caves, which harbor a high diversity of cave-dwelling organisms and are regarded as a biodiversity hotspot. Here, we assess the temporal dynamics of biotic colonization of subtropical East Asian caves through a multi-taxon analysis with representatives of green plants, animals, and fungi. We then investigate the consequences of paleonviromental changes on the colonization dynamics of these caves in combination with reconstructions of vegetation, temperature, and precipitation. We discover that 88% of cave colonization events occurred after the Oligocene-Miocene boundary, and organisms from the surrounding forest were a major source for subtropical East Asian cave biodiversity. Biotic colonization of subtropical East Asian caves during the Neogene was subject to periods of acceleration and decrease, in conjunction with large-scale, seasonal climatic changes and evolution of local forests. This study highlights the long-term evolutionary interaction between surface and cave biotas; our climate-vegetation-relict model proposed for the subtropical East Asian cave biota may help explain the evolutionary origins of other mid-latitude subterranean biotas.

Immigration dynamics of tropical and subtropical Southeast Asian limestone karst floras

Ex situ origins and dispersal of taxa have played important roles in the assembly of island-like biodiversity hotspots. Insular limestone karsts in Southeast Asia are hotspots of biodiversity and endemism, but the immigration processes of their unique floras are still poorly known. Here, we used Gesneriaceae as a proxy to investigate the immigration dynamics of tropical and subtropical Southeast Asian karst floras. We present the most comprehensive phylogenetic analysis of the Old World gesneriads to date based on twelve loci. By estimating divergence times and reconstructing ancestral states (habitat, soil type and range), we found that immigration into subtropical Southeast Asian karst floras first occurred in the Early Miocene, with two peaks in the Early-Middle Miocene and the Pliocene-Early Pleistocene, whereas immigration into tropical Southeast Asian karsts initiated in the Late Eocene, with two peaks in the Late Oligocene and the Late Miocene. We also discover that Southeast Asian karst biodiversity comprises immigrant pre-adapted lineages and descendants from local acid soil ancestors, although niche shift from acid soil to karst in tropical Southeast Asian islands was lacking. This study advances our understanding of the historical assembly of Southeast Asian karst floras.

Ancient volcanos as species pumps: A case study of freshwater amphipods in Northeast Asia

Volcano-tectonic processes have been viewed as primary drivers in the formation of present-day diversity. Volcanos associated with mountain uplifts drive allopatric speciation through vicariance and may impact the surrounding areas like species pump or species attractor. However, the application of these hypotheses to aquatic fauna has rarely been tested explicitly. We tested these hypotheses in the Changbai Mountains (Mts), which are one of the most typical, active volcanic ranges in Northeast (NE) Asia with a long and turbulent geological history. The Gammarus nekkensis species complex of amphipod crustaceans, widely distributed throughout NE Asia with poor dispersal abilities and a long evolutionary history, is a suitable model for testing hypotheses of species pump or species attractor. Phylogenetic and ancestral range reconstructions demonstrated that the studied amphipod originated from the Changbai Mts ~27 Ma and diverged into eastern (Clade I) and western (Clade II) clades, which corresponds well with the initial volcanic eruption of the Changbai Mts in the Late Oligocene. The subsequent diversifications of subclades CI-3, CII-1a and CII-2a were probably driven by second and third eruptions of the Changbai Mts during the Miocene. In particular, the Changbai lineages had spread to the Russian Far East multiple times since the Early Miocene, and widely colonized the region during the Pleistocene. Our discoveries suggest that the ancient volcanos of the Changbai Mts act as species pumps in NE Asia, resulted in burst of diversification around the Changbai Mts and subsequent dispersals into adjacent regions.

The mid-Miocene Zhangpu biota reveals an outstandingly rich rainforest biome in East Asia

During the Mid-Miocene Climatic Optimum [MMCO, ~14 to 17 million years (Ma) ago], global temperatures were similar to predicted temperatures for the coming century. Limited megathermal paleoclimatic and fossil data are known from this period, despite its potential as an analog for future climate conditions. Here, we report a rich middle Miocene rainforest biome, the Zhangpu biota (~14.7 Ma ago), based on material preserved in amber and associated sedimentary rocks from southeastern China. The record shows that the mid-Miocene rainforest reached at least 24.2°N and was more widespread than previously estimated. Our results not only highlight the role of tropical rainforests acting as evolutionary museums for biodiversity at the generic level but also suggest that the MMCO probably strongly shaped the East Asian biota via the northern expansion of the megathermal rainforest biome. The Zhangpu biota provides an ideal snapshot for biodiversity redistribution during global warming.

Global Diversification of Anelosimus Spiders Driven by Long-Distance Overwater Dispersal and Neogene Climate Oscillations

Vicariance and dispersal events, combined with intricate global climatic history, have left an imprint on the spatiotemporal distribution and diversity of many organisms. Anelosimus cobweb spiders (Theridiidae), are organisms ranging in behavior from solitary to highly social, with a cosmopolitan distribution in temperate to tropical areas. Their evolutionary history and the discontinuous distribution of species richness suggest that 1) long-distance overwater dispersal and 2) climate change during the Neogene (23-2.6 Ma), may be major factors in explaining their distribution and diversification. Here, we test these hypotheses, and explicitly test if global Miocene/Pliocene climatic cooling in the last 8 Ma affected Anelosimus radiation in parallel in South America and Madagascar. To do so, we investigate the phylogeny and spatiotemporal biogeography of Anelosimus through a culmination of a 20-year comprehensive global sampling at the species level (69 species, including 84% of the known 75 species worldwide, represented by 268 individuals) using nucleotide data from seven loci (5.5 kb). Our results strongly support the monophyly of Anelosimus with an Oligocene ($\sim $30 Ma) South American origin. Major clades on other continents originate via multiple, long-distance dispersal events, of solitary or subsocial-but not social-lineages, from the Americas. These intercontinental dispersals were to Africa, Madagascar (twice), and SE Asia/Australasia. The early diversification of Anelosimus spiders coincides with a sudden thermal increase in the late Oligocene ($\sim $27-25 Ma), though no causal connection can be made. Our results, however, strongly support the hypothesis that global Neogene climatic cooling in the last 8 Ma drove Anelosimus radiation in parallel in South America and Madagascar, offering a rare empirical evidence for diversification of a socially diverse group driven by an interplay between long-distance dispersal and global Neogene climatic changes. [Cobweb spiders; diversification; global biogeography; long-distance dispersal; molecular phylogenetics; neogene climate changes; sociality; vicariance.].

Climatic and topographic changes since the Miocene influenced the diversification and biogeography of the tent tortoise (Psammobates tentorius) species complex in Southern Africa

BACKGROUND

Climatic and topographic changes function as key drivers in shaping genetic structure and cladogenic radiation in many organisms. Southern Africa has an exceptionally diverse tortoise fauna, harbouring one-third of the world's tortoise genera. The distribution of Psammobates tentorius (Kuhl, 1820) covers two of the 25 biodiversity hotspots in the world, the Succulent Karoo and Cape Floristic Region. The highly diverged P. tentorius represents an excellent model species for exploring biogeographic and radiation patterns of reptiles in Southern Africa.

RESULTS

We investigated genetic structure and radiation patterns against temporal and spatial dimensions since the Miocene in the Psammobates tentorius species complex, using multiple types of DNA markers and niche modelling analyses. Cladogenesis in P. tentorius started in the late Miocene (11.63-5.33 Ma) when populations dispersed from north to south to form two geographically isolated groups. The northern group diverged into a clade north of the Orange River (OR), followed by the splitting of the group south of the OR into a western and an interior clade. The latter divergence corresponded to the intensification of the cold Benguela current, which caused western aridification and rainfall seasonality. In the south, tectonic uplift and subsequent exhumation, together with climatic fluctuations seemed responsible for radiations among the four southern clades since the late Miocene. We found that each clade occurred in a habitat shaped by different climatic parameters, and that the niches differed substantially among the clades of the northern group but were similar among clades of the southern group.

CONCLUSION

Climatic shifts, and biome and geographic changes were possibly the three major driving forces shaping cladogenesis and genetic structure in Southern African tortoise species. Our results revealed that the cladogenesis of the P. tentorius species complex was probably shaped by environmental cooling, biome shifts and topographic uplift in Southern Africa since the late Miocene. The Last Glacial Maximum (LGM) may have impacted the distribution of P. tentorius substantially. We found the taxonomic diversify of the P. tentorius species complex to be highest in the Greater Cape Floristic Region. All seven clades discovered warrant conservation attention, particularly Ptt-B-Ptr, Ptt-A and Pv-A.

Biogeography of land snail genus Acusta (Gastropoda: Camaenidae): Diversification on East Asian islands

East Asia has highly diverse and endemic biota due to its complex geological and climatic history and its diversified topography. The continental and insular distributions of land snail genus Acusta in East Asia provide a good opportunity to compare the evolutionary processes in this group under different biogeographical conditions. In this study, we inferred the evolutionary history of the land snail genus Acusta by a molecular phylogeny and investigated how the palaeogeographic events shaped species diversity and the distribution of the Acusta genus within the island arc. A concatenated dataset generated from sequences of one nuclear (ITS2) and two mitochondrial (16S, COI) gene fragments, include most of nominal taxa of the genus, four related species and one outgroup. We constructed the phylogeny and the evolutionary history of the genus through maximum-likelihood and Bayesian inference methods, using a Bayesian molecular clock and ancestral range estimation. Our results suggested that currently recognized species in Acusta are polyphyletic. The traditionally accepted concept of the affinity of Acusta and Bradybaena is not supported. The hypothesis of colonization via land bridges during the Pleistocene glaciations for the biota of East Asian islands is not supported. Instead, the origin and diversification of the genus Acusta was dated to the late Miocene-Pliocene from an area around North and Northeast China to South China and East Asian islands Three major evolutionary lineages were identified. Two of the major lineages demonstrate distinct evolutionary histories, as sympatric speciation is the major speciation process for the continental clade, while the insular clade originated from founder events. Taiwan functioned as an important source of diversification for species on the East Asian islands possibly through passive dispersal of different mechanisms. The sea level fluctuations caused by the Pleistocene glacial cycles play a role in the subsequent dispersion and diversification of species of the continental clade, such as the more recent range expansion of A. redfieldi from South China to Taiwan and Japan.

Phylogeography of the cicada Platypleura hilpa in subtropical and tropical East Asia based on mitochondrial and nuclear genes and microsatellite markers

Climate change and geographical events play key roles in driving population genetic structure of organisms, but different scenarios were suggested for species occurring in tropical and subtropical areas. We investigated the population genetic structure, diversity and demographic history of the cicada Platypleura hilpa which occurs in coastal areas of southern China and northeastern Indo-Burma, and analysed the potential impact of climate change and geological events on its evolutionary history. Our data imply that P. hilpa comprises five main lineages with nearly unique sets of haplotypes and distinct geographic distributions. A major split of the lineages occurred in the Pleistocene. Geographic distance and geomorphic barriers serve as the primary factors shaping the genetic population structure, and several climatic factors are associated with the divergence. The potential range during the Last Glacial Maximum has apparently increased in south China and the exposed South China Sea Shelf. The Pleistocene sea-level fluctuations had profound effects on the regional genetic structure. The Beibu Gulf represents a more important geographic barrier than the Qiongzhou Strait in blocking gene flow among populations. These results contribute to a better understanding of the pressure climatic change and geographical events impose on insects in coastal areas of East Asia.

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.