Genomic data revealed inbreeding despite a geographically connected stable effective population size since the Holocene in the protected Formosan Long-Arm Scarab beetle, Cheirotonus formosanus.

Biodiversity conservation is a top priority in the face of global environmental change, and the practical restoration of biodiversity has emerged as a key objective. Nevertheless, the question of how to effectively contribute to biodiversity restoration and identify suitable systems for such efforts continues to present major challenges. By using genome-wide SNP data, our study revealed that populations from different mountain ranges of the Formosan Long-Arm Scarab beetle, a flagship species that receives strict protection, exhibited a single genetic cluster with no subdivision. Additionally, our result implied an association between the demographic history and historical fluctuations in climate and environmental conditions. Furthermore, we showed that, despite a stable and moderately sized effective population over recent history, all the individuals we studied exhibited signs of genetic inbreeding. We argued that the current practice of protecting the species as one evolutionarily significant unit remains the best conservation plan and that recent habitat change may have led to the pattern of significant inbreeding. We closed by emphasizing the importance of conservation genetic studies in guiding policy decisions and highlighting the potential of genomic data for identifying ideal empirical systems for genetic rescue, or assisted gene flow studies.

Systematics and evolutionary dynamics of insect-fern interactions in the specialized fern-spore feeding Cuprininae (Lepidoptera, Stathmopodidae).

Fern-spore-feeding (FSF) is rare and found in only four families of Lepidoptera. Stathmopodidae is the most speciose family that contains FSF species, and its subfamily Cuprininae exclusively specializes on FSF. However, three species of Stathmopodinae also specialize on FSF. To better understand the evolutionary history of FSF and, more generally, the significance of specialization on a peculiar host, a phylogenetic and taxonomic revision for this group is necessary. We reconstructed the most comprehensive molecular phylogeny, including one mitochondrial and four nuclear genes, of Stathmopodidae to date, including 137 samples representing 62 species, with a particular focus on the FSF subfamily, Cuprininae, including 33 species (41% of named species) from 6 of the 7 Cuprininae genera. Species from two other subfamilies, Stathmopodinae and Atkinsoniinae, were also included. We found that FSF evolved only once in Stathmopodidae and that the previous hypothesis of multiple origins of FSF was misled by inadequate taxonomy. Moreover, we showed that (1) speciation/extinction rates do not differ significantly between FSF and non-FSF groups and that (2) oligophage is the ancestral character state in Cuprininae. We further revealed that a faster rate of accumulating specialists over time, and thus a higher number of specialists, was achieved by a higher transition rate from oligophagages to specialists compared to the transition rate in the opposite direction. We finish by describing three new genera, Trigonodagen. nov., Petalagen. nov., and Pediformisgen. nov., and revalidating five genera: Cuprina, Calicotis, Thylacosceles, Actinoscelis, Thylacosceloides in Cuprininae, and we provide an updated taxonomic key to genera and a revised global checklist of Cuprininae.

The genomes of Hercules beetles reveal putative adaptive loci and distinct demographic histories in pristine North American forests.

Beetles, despite their remarkable biodiversity and a long history of research, remain lacking in reference genomes annotated with structural variations in loci of adaptive significance. We sequenced and assembled high-quality chromosome-level genomes of four Hercules beetles which exhibit divergence in male horn size and shape and body colouration. The four Hercules beetle genomes were assembled to 11 pseudo-chromosomes, where the three genomes assembled using Nanopore data (Dynastes grantii, D. hyllus and D. tityus) were mapped to the genome assembled using PacBio + Hi-C data (D. maya). We demonstrated a striking similarity in genome structure among the four species. This conservative genome structure may be attributed to our use of the D. maya assembly as the reference; however, it is worth noting that such a conservative genome structure is a recurring phenomenon among scarab beetles. We further identified homologues of nine and three candidate-gene families that may be associated with the evolution of horn structure and body colouration respectively. Structural variations in Scr and Ebony2 were detected and discussed for their putative impacts on generating morphological diversity in beetles. We also reconstructed the demographic histories of the four Hercules beetles using heterozygosity information from the diploid genomes. We found that the demographic histories of the beetles closely recapitulated historical changes in suitable forest habitats driven by climate shifts.

Testing the efficacy of different molecular tools for parasite conservation genetics: a case study using horsehair worms (Phylum: Nematomorpha).

In recent years, parasite conservation has become a globally significant issue. Because of this, there is a need for standardized methods for inferring population status and possible cryptic diversity. However, given the lack of molecular data for some groups, it is challenging to establish procedures for genetic diversity estimation. Therefore, universal tools, such as double-digest restriction-site-associated DNA sequencing (ddRADseq), could be useful when conducting conservation genetic studies on rarely studied parasites. Here, we generated a ddRADseq dataset that includes all 3 described Taiwanese horsehair worms (Phylum: Nematomorpha), possibly one of the most understudied animal groups. Additionally, we produced data for a fragment of the cytochrome c oxidase subunit I (COXI) for the said species. We used the COXI dataset in combination with previously published sequences of the same locus for inferring the effective population size (Ne) trends and possible population genetic structure.We found that a larger and geographically broader sample size combined with more sequenced loci resulted in a better estimation of changes in Ne. We were able to detect demographic changes associated with Pleistocene events in all the species. Furthermore, the ddRADseq dataset for Chordodes formosanus did not reveal a genetic structure based on geography, implying a great dispersal ability, possibly due to its hosts. We showed that different molecular tools can be used to reveal genetic structure and demographic history at different historical times and geographical scales, which can help with conservation genetic studies in rarely studied parasites.

Genetic Assessment of Fertile F1 Hybrids between Two Hercules Beetles, Dynastes maya Hardy and D. grantii Horn (Scarabaeidae).

Although hybridization may complicate taxonomic practices, it can be common between animal species. Animal hybridization not only can help with generating phenotypic and species diversity in nature, but also with understanding the genetic and genomic basis of phenotypic evolution in the laboratory. We assessed the genetic composition of captive bred F1 hybrids between two Hercules beetle species using mitochondrial CO1 and nuclear loci from a double-digest restriction-site associated DNA sequencing (ddRADseq) library. We showed that the F1 hybrids were genetically clustered with samples from the maternal species, D. grantii, based on CO1 data. Nuclear genome data, on the other hand, clearly showed that the F1 individuals were genetically intermediate between D. maya, the paternal species, and D. grantii, based on a principal component analysis. Our results also revealed that sampling design may have a major impact on the inferred genetic structure and hybrid individuals using ddRADseq data sets. We discuss the importance and potential from studying the genomics of this hybrid progeny in terms of understanding the origin and maintenance of both intraspecific and interspecific phenotypic divergence and convergence.

Modeling the geographical distributions of Chordodes formosanus and its mantis hosts in Taiwan, with considerations for their niche overlaps.

Species distribution models (SDMs) have conventionally been used for evaluating the distribution of individual species, but they can also be used, through comparing different SDMs, to evaluate the geographic similarity between taxa. In this study, we used a parasite and host system to infer the geographic overlaps between species with tight biological interaction, for example, parasites and their obligate host. Specifically, we used the horsehair worm Chordodes formosanus and its three mantis hosts to study the extent of niche overlap. We retrieved presence points for the host species and the parasite, and then we built SDMs with MaxEnt implemented in ENMeval using selected bioclim variables (based on variance inflation factor values) at 30s scale. The models showed that the hosts and parasite do not occur in the high elevation areas in Taiwan, which is expected based on their biology. Interestingly, the predicted parasite distribution included areas without collection records, implying local extinction or sampling bias. We subsequently evaluated niche overlap between hosts and the parasite according to five similarity indices (Schoener's D, I statistic, relative rank, Pearson correlation coefficient, and the rank correlation coefficient rho). Our models showed a high similarity of SDM predictions between hosts and the parasite. There were differences among metrics for which host shared the highest similarity with the parasite, but the majority of the results indicated that the Japanese boxing mantis had the highest niche similarity with the horsehair worm. The choice of the niche overlap metric to use can uncover information on the parasite's ecology, which can be important for endangered species. SDMs are reliable tools for host and parasite conservation management and could help improve our understanding of parasite biology and ecology.

Long-read Sequencing Data Reveals Dynamic Evolution of Mitochondrial Genome Size and the Phylogenetic Utility of Mitochondrial DNA in Hercules Beetles (Dynastes; Scarabaeidae).

The evolutionary dynamics and phylogenetic utility of mitochondrial genomes (mitogenomes) have been of particular interest to systematists and evolutionary biologists. However, certain mitochondrial features, such as the molecular evolution of the control region in insects, remain poorly explored due to technological constraints. Using a combination of long- and short-read sequencing data, we assembled ten complete mitogenomes from ten Hercules beetles. We found large-sized mitogenomes (from 24 to 28 kb), which are among the largest in insects. The variation in genome size can be attributed to copy-number evolution of tandem repeats in the control region. Furthermore, one type of tandem repeat was found flanking the conserved sequence block in the control region. Importantly, such variation, which made up around 30% of the size of the mitogenome, may only become detectable should long-read sequencing technology be applied. We also found that, although different mitochondrial loci often inferred different phylogenetic histories, none of the mitochondrial loci statistically reject a concatenated mitochondrial phylogeny, supporting the hypothesis that all mitochondrial loci share a single genealogical history. We on the other hand reported statistical support for mito-nuclear phylogenetic discordance in 50% of mitochondrial loci. We argue that long-read DNA sequencing should become a standard application in the rapidly growing field of mitogenome sequencing. Furthermore, mitochondrial gene trees may differ even though they share a common genealogical history, and ND loci could be better candidates for phylogenetics than the commonly used COX1.

Interpreting phylogenetic conflict: Hybridization in the most speciose genus of lichen-forming fungi.

While advances in sequencing technologies have been invaluable for understanding evolutionary relationships, increasingly large genomic data sets may result in conflicting evolutionary signals that are often caused by biological processes, including hybridization. Hybridization has been detected in a variety of organisms, influencing evolutionary processes such as generating reproductive barriers and mixing standing genetic variation. Here, we investigate the potential role of hybridization in the diversification of the most speciose genus of lichen-forming fungi, Xanthoparmelia. As Xanthoparmelia is projected to have gone through recent, rapid diversification, this genus is particularly suitable for investigating and interpreting the origins of phylogenomic conflict. Focusing on a clade of Xanthoparmelia largely restricted to the Holarctic region, we used a genome skimming approach to generate 962 single-copy gene regions representing over 2 Mbp of the mycobiont genome. From this genome-scale dataset, we inferred evolutionary relationships using both concatenation and coalescent-based species tree approaches. We also used three independent tests for hybridization. Although different species tree reconstruction methods recovered largely consistent and well-supported trees, there was widespread incongruence among individual gene trees. Despite challenges in differentiating hybridization from ILS in situations of recent rapid radiations, our genome-wide analyses detected multiple potential hybridization events in the Holarctic clade, suggesting one possible source of trait variability in this hyperdiverse genus. This study highlights the value in using a pluralistic approach for characterizing genome-scale conflict, even in groups with well-resolved phylogenies, while highlighting current challenges in detecting the specific impacts of hybridization.

Extensive sampling of Saccharomyces cerevisiae in Taiwan reveals ecology and evolution of predomesticated lineages.

The ecology and genetic diversity of the model yeast Saccharomyces cerevisiae before human domestication remain poorly understood. Taiwan is regarded as part of this yeast's geographic birthplace, where the most divergent natural lineage was discovered. Here, we extensively sampled the broadleaf forests across this continental island to probe the ancestral species' diversity. We found that S. cerevisiae is distributed ubiquitously at low abundance in the forests. Whole-genome sequencing of 121 isolates revealed nine distinct lineages that diverged from Asian lineages during the Pleistocene, when a transient continental shelf land bridge connected Taiwan to other major landmasses. Three lineages are endemic to Taiwan and six are widespread in Asia, making this region a focal biodiversity hotspot. Both ancient and recent admixture events were detected between the natural lineages, and a genetic ancestry component associated with isolates from fruits was detected in most admixed isolates. Collectively, Taiwanese isolates harbor genetic diversity comparable to that of the whole Asia continent, and different lineages have coexisted at a fine spatial scale even on the same tree. Patterns of variations within each lineage revealed that S. cerevisiae is highly clonal and predominantly reproduces asexually in nature. We identified different selection patterns shaping the coding sequences of natural lineages and found fewer gene family expansion and contractions that contrast with domesticated lineages. This study establishes that S. cerevisiae has rich natural diversity sheltered from human influences, making it a powerful model system in microbial ecology.

First detection of Colletotrichumfructicola (Ascomycota) on horsehair worms (Nematomorpha).

Fungal members of Colletotrichum (Ascomycota) were found to be associated with Chordodesformosanus, one of the three currently known horsehair worm (Nematomorpha) species in Taiwan. The fungi were identified as Colletotrichumfructicola, which is mostly known as a plant pathogen, through the use of the nuclear ribosomal internal transcribed spacer and partial large subunit (nrITS + nrLSU) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) DNA sequences. To our knowledge, this report represents both the first records for Colletotrichum associated with hairworms and for fungi on Nematomorpha. These findings expand the knowledge on the ecological relationships of both clades.

Evolution of adult male horn developmental phenotypes and character displacement in Xylotrupes beetles (Scarabaeidae).

Character displacement that leads to divergent phenotypes between sympatric species has been hypothesized to facilitate coexistence and promote the accumulation of biodiversity. However, there are alternative evolutionary mechanisms that may also lead to the evolution of phenotypic divergence between sympatric species; one of the mechanisms is evolutionary contingency. We studied the evolution of the presence and absence of a major male horn phenotype, which may have ecological implications for promoting coexistence between sympatric beetles, across geographic populations from different Xylotrupes beetles. By using a previously published phylogeny with 80 Xylotrupes taxa, we estimated the transition rates between the two phenotypic states (i.e., presence vs. absence of a major male phenotype). Based on the estimated transition rates, we then simulated possible phenotypic outcomes between sympatric species. We found that sympatric species were equally likely to evolve the same versus distinct phenotypic states based on the estimated transition rates given the phylogeny. The empirically observed number of sympatric species showing different phenotypic states can be explained by evolutionary contingency alone. We discussed the importance of applying phylogenetic comparative methods when studying phenotypic evolution and more generally to investigate the effect of stochastic processes before making deterministic inferences.

The Origins of Coca: Museum Genomics Reveals Multiple Independent Domestications from Progenitor Erythroxylum gracilipes.

Coca is the natural source of cocaine as well as a sacred and medicinal plant farmed by South American Amerindians and mestizos. The coca crop comprises four closely related varieties classified into two species (Amazonian and Huánuco varieties within Erythroxylum coca Lam., and Colombian and Trujillo varieties within Erythroxylum novogranatense (D. Morris) Hieron.) but our understanding of the domestication and evolutionary history of these taxa is nominal. In this study, we use genomic data from natural history collections to estimate the geographic origins and genetic diversity of this economically and culturally important crop in the context of its wild relatives. Our phylogeographic analyses clearly demonstrate the four varieties of coca comprise two or three exclusive groups nested within the diverse lineages of the widespread, wild species Erythroxylum gracilipes; establishing a new and robust hypothesis of domestication wherein coca originated two or three times from this wild progenitor. The Colombian and Trujillo coca varieties are descended from a single, ancient domestication event in northwestern South America. Huánuco coca was domesticated more recently, possibly in southeastern Peru. Amazonian coca either shares a common domesticated ancestor with Huánuco coca, or it was the product of a third and most recent independent domestication event in the western Amazon basin. This chronology of coca domestication reveals different Holocene peoples in South America were able to independently transform the same natural resource to serve their needs; in this case, a workaday stimulant. [Erythroxylum; Erythroxylaceae; Holocene; Museomics; Neotropics; phylogeography; plant domestication; target-sequence capture.].

Macroevolutionary integration of phenotypes within and across ant worker castes.

Phenotypic traits are often integrated into evolutionary modules: sets of organismal parts that evolve together. In social insect colonies, the concepts of integration and modularity apply to sets of traits both within and among functionally and phenotypically differentiated castes. On macroevolutionary timescales, patterns of integration and modularity within and across castes can be clues to the selective and ecological factors shaping their evolution and diversification. We develop a set of hypotheses describing contrasting patterns of worker integration and apply this framework in a broad (246 species) comparative analysis of major and minor worker evolution in the hyperdiverse ant genus Pheidole. Using geometric morphometrics in a phylogenetic framework, we inferred fast and tightly integrated evolution of mesosoma shape between major and minor workers, but slower and more independent evolution of head shape between the two worker castes. Thus, Pheidole workers are evolving as a mixture of intracaste and intercaste integration and rate heterogeneity. The decoupling of homologous traits across worker castes may represent an important process facilitating the rise of social complexity.

Is population subdivision different from speciation? From phylogeography to species delimitation.

Species-level diversity and the underlying mechanisms that lead to the formation of new species, that is, speciation, have often been confounded with intraspecific diversity and population subdivision. The delineation between intraspecific and interspecific divergence processes has received much less attention than species delimitation. The ramifications of confounding speciation and population subdivision are that the term speciation has been used to describe many different biological divergence processes, rendering the results, or inferences, between studies incomparable. Phylogeographic studies have advanced our understanding of how spatial variation in the pattern of biodiversity can begin, become structured, and persist through time. Studies of species delimitation have further provided statistical and model-based approaches to determine the phylogeographic entities that merit species status. However, without a proper understanding and delineation between the processes that generate and maintain intraspecific and interspecific diversity in a study system, the delimitation of species may still not be biologically and evolutionarily relevant. I argue that variation in the continuity of the divergence process among biological systems could be a key factor leading to the enduring contention in delineating divergence patterns, or species delimitation, meriting future comparative studies to help us better understand the nature of biological species.

Socially Parasitic Ants Evolve a Mosaic of Host-Matching and Parasitic Morphological Traits.

A basic expectation of evolution by natural selection is that species morphologies will adapt to their ecological niche. In social organisms, this may include selective pressure from the social environment. Many non-ant parasites of ant colonies are known to mimic the morphology of their host species, often in striking fashion [1, 2], indicating there is selection on parasite morphology to match the host (Batesian and/or Wasmannian mimicry [3]). However, ants that parasitize other ant societies are usually closely related to their hosts (Emery's rule) [4-8] and expected to be similar due to common ancestry, making any kind of mimicry difficult to detect [9]. Here, we investigate the diversification of the hyperdiverse ant genus Pheidole in Madagascar, including the evolution of 13 putative social parasite species within a broader radiation of over 100 ant species on the island. We find that the parasitic species are monophyletic and that their associated hosts are spread across the Malagasy Pheidole radiation. This provides an opportunity to test for selection on morphological similarity and divergence between parasites and hosts. Using X-ray microtomography and both linear measurements and three-dimensional (3D) geometric morphometrics, we show that ant social parasite worker morphologies feature a mix of "host-matching" and "parasitic" traits, where the former converge on the host phenotype and the latter diverge from typical Pheidole phenotypes to match a common parasitic syndrome. This finding highlights the role of social context in shaping the evolution of phenotypes and raises questions about the role of morphological sensing in nestmate recognition.

Genome-scale data reveal the role of hybridization in lichen-forming fungi.

Advancements in molecular genetics have revealed that hybridization may be common among plants, animals, and fungi, playing a role in evolutionary dynamics and speciation. While hybridization has been well-documented in pathogenic fungi, the effects of these processes on speciation in fungal lineages with different life histories and ecological niches are largely unexplored. Here we investigated the potential influence of hybridization on the emergence of morphologically and reproductively distinct asexual lichens. We focused on vagrant forms (growing obligately unattached to substrates) within a clade of rock-dwelling, sexually reproducing species in the Rhizoplaca melanophthalma (Lecanoraceae, Ascomycota) species complex. We used phylogenomic data from both mitochondrial and nuclear genomes to infer evolutionary relationships and potential patterns of introgression. We observed multiple instances of discordance between the mitochondrial and nuclear trees, including the clade comprising the asexual vagrant species R. arbuscula, R. haydenii, R. idahoensis, and a closely related rock-dwelling lineage. Despite well-supported phylogenies, we recovered strong evidence of a reticulated evolutionary history using a network approach that incorporates both incomplete lineage sorting and hybridization. These data suggest that the rock-dwelling western North American subalpine endemic R. shushanii is potentially the result of a hybrid speciation event, and introgression may have also played a role in other taxa, including vagrant species R. arbuscula, R. haydenii and R. idahoensis. We discuss the potential roles of hybridization in terms of generating asexuality and novel morphological traits in lichens. Furthermore, our results highlight the need for additional study of reticulate phylogenies when investigating species boundaries and evolutionary history, even in cases with well-supported topologies inferred from genome-scale data.

Accelerated diversifications in three diverse families of morphologically complex lichen-forming fungi link to major historical events.

Historical mass extinction events had major impacts on biodiversity patterns. The most recent and intensively studied event is the Cretaceous - Paleogene (K-Pg) boundary (ca. 66 million years ago [MYA]). However, the factors that may have impacted diversification dynamics vary across lineages. We investigated the macroevolutionary dynamics with a specific focus on the impact of major historical events such as the K-Pg mass extinction event on two major subclasses - Lecanoromycetidae and Ostropomycetidae - of lichen-forming fungi and tested whether variation in the rate of diversification can be associated with the evolution of a specific trait state - macrolichen. Our results reveal accelerated diversification events in three families of morphologically complex lichen-forming fungi - Cladoniaceae, Parmeliaceae, and Peltigeraceae - which are from the subclass Lecanoromycetidae and mostly composed of macrolichens, those that form three dimensional structures. Our RTT plot result for the subclass Lecanoromycetidae also reveals accelerated diversification. Changes in diversification rates occurred around the transition between Mesozoic and Cenozoic eras and was likely related to the K-Pg mass extinction event. The phylogenetic positions for rate increases estimated based on marginal shift probability are, however, scattered from 100 to 40 MYA preventing us from making explicit inference. Although we reveal that the phenotypic state of macrolichens is associated with a higher diversification rate than microlichens, we also show that the evolution of macrolichens predated the K-Pg event. Furthermore, the association between macrolichens and increased diversification is not universal and can be explained, in part, by phylogenetic relatedness. By investigating the macroevolutionary dynamics of lichen-forming fungi our study provides a new empirical system suitable to test the effect of major historical event on shaping biodiversity patterns and to investigate why changes in biodiversity patterns are not in concordance across clades. Our results imply that multiple historical events during the transition from Mesozoic to Cenozoic eras, including the K-Pg mass extinction event, impacted the evolutionary dynamics in lichen-forming fungi. However, future studies focusing on individual lichen-forming fungal families are required to ascertain whether diversification rates are associated with growth form and certain geological events.

Multiple historical processes obscure phylogenetic relationships in a taxonomically difficult group (Lobariaceae, Ascomycota).

In the age of next-generation sequencing, the number of loci available for phylogenetic analyses has increased by orders of magnitude. But despite this dramatic increase in the amount of data, some phylogenomic studies have revealed rampant gene-tree discordance that can be caused by many historical processes, such as rapid diversification, gene duplication, or reticulate evolution. We used a target enrichment approach to sample 400 single-copy nuclear genes and estimate the phylogenetic relationships of 13 genera in the lichen-forming family Lobariaceae to address the effect of data type (nucleotides and amino acids) and phylogenetic reconstruction method (concatenation and species tree approaches). Furthermore, we examined datasets for evidence of historical processes, such as rapid diversification and reticulate evolution. We found incongruence associated with sequence data types (nucleotide vs. amino acid sequences) and with different methods of phylogenetic reconstruction (species tree vs. concatenation). The resulting phylogenetic trees provided evidence for rapid and reticulate evolution based on extremely short branches in the backbone of the phylogenies. The observed rapid and reticulate diversifications may explain conflicts among gene trees and the challenges to resolving evolutionary relationships. Based on divergence times, the diversification at the backbone occurred near the Cretaceous-Paleogene (K-Pg) boundary (65 Mya) which is consistent with other rapid diversifications in the tree of life. Although some phylogenetic relationships within the Lobariaceae family remain with low support, even with our powerful phylogenomic dataset of up to 376 genes, our use of target-capturing data allowed for the novel exploration of the mechanisms underlying phylogenetic and systematic incongruence.

Holocene Population Decline and Conservation Implication for the Western Hercules Beetle, Dynastes grantii (Coleoptera, Scarabaeidae).

The Western Hercules beetle (Dynastes grantii) is endemic to the highland forest habitats of southwestern United States and northern Mexico. The habitats harbor many endemic species, but are being threatened by rapid climate change and urban development. In this study, the genetic structure of D. grantii populations from southwestern United States was investigated. Specifically, genomic data from double-digest restriction-site-associated DNA sequencing libraries were utilized to test whether geographically distant populations from the Mogollon Rim (Arizona [N = 12 individuals] and New Mexico [N = 10 individuals]) are genetically structured. The study also estimated the effective population size of the Mogollon Rim populations based on genetic diversity. The results indicated that the 2 geographic populations from the Mogollon Rim were not genetically structured. A population size reduction was detected since the end of the last glacial period, which coincided with a reduction of forest habitat in the study area. The results implied that the connectivity and the size of highland forest habitats in the Mogollon Rim could have been the major factors shaping the population genetic structure and demographic history of D. grantii. The Western Hercules beetle could be a useful flagship species for local natural history education and to promote the conservation of highland forest habitats.

Characterisation of the complete mitochondrial genome of Lucanus chengyuani (Coleoptera: Lucanidae).

We sequenced and assembled the complete mitochondrial genome of Lucanus chengyuani, from the Alishan, Chiayi County, Taiwan. The length of the complete mitogenome of L. chengyuani is 16,926 bp and the mitogenome contains 13 protein-coding, 22 tRNA, and 2 rDNA genes. Nucleotide compositions of the whole mitogenome of L. chengyuani are 38.37% for A, 27.96% for T, 23.03% for C, and 10.637% for G. The AT and GC skewness of mitogenome sequence are 0.157 and -0.368, showing the A-skew and C-skew. The reconstructed phylogenetic relationships of 9 Lucanidae species based on 13 mitochondrial protein-coding genes are highly supported. The clade including Neolucanus maximus and Odontolabis cuvera is sister to the rest of the stag beetle clades, which contains L. chengyuani and L. mazama. Mitogenomic data from this study will provide useful information for further studies for the population genetics, speciation, biogeography, and conservation of L. chengyuani in the future.