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.

The hindgut microbiota of coconut rhinoceros beetles (Oryctes rhinoceros) in relation to their geographical populations.

The coconut rhinoceros beetle (CRB, Oryctes rhinoceros) is a palm tree pest capable of rapidly expanding its population in new territories. Previous studies identified a digestive symbiosis between CRB and its gut microbes. However, no research compared the genetic variation of CRBs with their hindgut microbiota on a global scale. This study aims to investigate the genetic divergence of CRB and the compositional variation of CRB's microbiota across different geographical locations, and explore the association between them and their predicted functional profiles and environmental data. The research reveals a distinct and consistent microbial community within local populations, but it varies across different geographical populations. The microbial functional profiles linked to the production of digestive enzymes, including cellulases and ligninases, are nonetheless globally conserved. This suggests that CRBs employ specific mechanisms to select and maintain microbes with functional benefits, contributing to host adaptability, stress tolerance, and fitness. The CRB microbial communities did not appear to recapitulate the genetic variation of their hosts. Rather than depend on obligate symbionts, CRBs seem to establish similar digestive associations with whatever environmentally acquired microbes are available wherever they are, aiding them in successfully establishing after invading a new location.IMPORTANCECoconut rhinoceros beetles (CRBs) are notorious pests on Arecaceae plants, posing destructive threats to countries highly reliant on coconut, oil palm, and date palm as economic crops. In the last century, CRBs have rapidly expanded their presence to territories that were once free of these beetles. The United States, for instance, has officially designated CRBs as invasive and alien pests. Given their remarkable ability to swiftly adapt to new environments, their gut microbes may play a crucial role in this process. While the microbiota of CRBs vary depending on geographical location, these beetles consistently exhibit a functionally identical digestive association with locally acquired microbes. This underscores the significance of CRB-microbe association in shaping the adaptive strategies of this agricultural pest.

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.

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.

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.

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.

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.].

Different species or genetically divergent populations? Integrative species delimitation of the Primulina hochiensis complex from isolated karst habitats.

To consistently and objectively delineate species-level divergence from population subdivision has been a challenge in systematics. This is particularly evident in naturally fragmented and allopatric systems in which small population size often leads to extreme population structuring. Here we evaluated the robustness of the species delimitation methods implemented in BEAST, BPP, and iBPP in the Primulina hochiensis complex comprising four described and one candidate species (five taxa in total) distributed in karst landscapes of southern China. We analyzed levels of molecular and morphological divergence among species using multilocus sequence data (nine chloroplast loci and 10 nuclear loci), and morphological data (16 quantitative and 12 qualitative traits), for 124 individuals from 25 populations of the complex. Independent analyses of cpDNA and nDNA sequence data revealed high levels of genetic differentiation among the five taxa. Both BPP and iBPP delimited five candidate species, which correspond to the five genetic clusters recovered with population structure analysis. In contrast, morphological differences among populations were more limited, so that results from principal component analysis (PCA) recovered only three distinct clusters. We ruled out the possibility of morphologically cryptic species because reciprocally monophyletic groups were not supported among the morphologically un-differentiated taxa. Our results represent a case where extreme population genetic structuring leads to oversplit of species diversity by molecular data using the multispecies coalescent (MSC) methods. The observed congruence across multiple analyses corroborates the recognition of a new species P. lianpingensis and indicates its sister species relationship with P. yingdeensis. This study highlights the dangers of violating model assumption and the importance of incorporating multiple evidence into species delimitation of a particular system.

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.

Oligocene origin and drivers of diversification in the genus Sticta (Lobariaceae, Ascomycota).

A major challenge to evolutionary biologists is to understand how biodiversity is distributed through space and time and across the tree of life. Diversification of organisms is influenced by many factors that act at different times and geographic locations but it is still not clear which have a significant impact and how drivers interact. To study diversification, we chose the lichen genus Sticta, by sampling through most of the global range and producing a time tree. We estimate that Sticta originated about 30 million years ago, but biogoegraphic analysis was unclear in estimating the origin of the genus. Furthermore, we investigated the effect of dispersal ability finding that Sticta has a high dispersal rate, as collections from Hawaii showed that divergent lineages colonized the islands at least four times. Symbiont interactions were investigated using BiSSE to understand if green-algal or cyanobacterial symbiont interactions influenced diversification, only to find that the positive results were driven almost completely by Type I error. On the other hand, another BiSSE analysis found that an association with Andean tectonic activity increases the speciation rate of species.

The Species versus Subspecies Conundrum: Quantitative Delimitation from Integrating Multiple Data Types within a Single Bayesian Approach in Hercules Beetles.

With the recent attention and focus on quantitative methods for species delimitation, an overlooked but equally important issue regards what has actually been delimited. This study investigates the apparent arbitrariness of some taxonomic distinctions, and in particular how species and subspecies are assigned. Specifically, we use a recently developed Bayesian model-based approach to show that in the Hercules beetles (genus Dynastes) there is no statistical difference in the probability that putative taxa represent different species, irrespective of whether they were given species or subspecies designations. By considering multiple data types, as opposed to relying exclusively on genetic data alone, we also show that both previously recognized species and subspecies represent a variety of points along the speciation spectrum (i.e., previously recognized species are not systematically further along the continuum than subspecies). For example, based on evolutionary models of divergence, some taxa are statistically distinguishable on more than one axis of differentiation (e.g., along both phenotypic and genetic dimensions), whereas other taxa can only be delimited statistically from a single data type. Because both phenotypic and genetic data are analyzed in a common Bayesian framework, our study provides a framework for investigating whether disagreements in species boundaries among data types reflect (i) actual discordance with the actual history of lineage splitting, or instead (ii) differences among data types in the amount of time required for differentiation to become apparent among the delimited taxa. We discuss what the answers to these questions imply about what characters are used to delimit species, as well as the diverse processes involved in the origin and maintenance of species boundaries. With this in mind, we then reflect more generally on how quantitative methods for species delimitation are used to assign taxonomic status.

Parapatric genetic introgression and phenotypic assimilation: testing conditions for introgression between Hercules beetles (Dynastes, Dynastinae).

The prevalence and consequences of genetic introgression between species have been intensively debated. I used Hercules beetles as examples to test for conditions that may be associated with the occurrence of introgression. RADseq data were used to reconstruct the species tree and history of introgression between Hercules beetles. Image data from museum specimens were used to investigate the phenotypic similarity of two adaptive traits between species from two distinct climatic realms (Nearctic vs. Neotropical). Genetic introgression was identified between Hercules beetles living in geographic proximity (parapatric). Phylogenetic relatedness and phenotypic similarity did not predict nor preclude genetic introgression between species. Phenotypic assimilation in body coloration was evident between distantly related Hercules beetles codistributed in Central America, where directional introgression was also statistically supported from the putative donor to receiver lineages. The number of introgressed loci was significantly higher between species with than without phenotypic similarity. I discuss the implications of recent studies on adaptive genetic introgression by providing supporting evidence from the Hercules beetle system.

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.

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.

Lineage-specific late pleistocene expansion of an endemic subtropical gossamer-wing damselfly, Euphaea formosa, in Taiwan.

BACKGROUND: Pleistocene glacial oscillations have significantly affected the historical population dynamics of temperate taxa. However, the general effects of recent climatic changes on the evolutionary history and genetic structure of extant subtropical species remain poorly understood. In the present study, phylogeographic and historical demographic analyses based on mitochondrial and nuclear DNA sequences were used. The aim was to investigate whether Pleistocene climatic cycles, paleo-drainages or mountain vicariance of Taiwan shaped the evolutionary diversification of a subtropical gossamer-wing damselfly, Euphaea formosa. RESULTS: E. formosa populations originated in the middle Pleistocene period (0.3 Mya) and consisted of two evolutionarily independent lineages. It is likely that they derived from the Pleistocene paleo-drainages of northern and southern Minjiang, or alternatively by divergence within Taiwan. The ancestral North-central lineage colonized northwestern Taiwan first and maintained a slowly growing population throughout much of the early to middle Pleistocene period. The ancestral widespread lineage reached central-southern Taiwan and experienced a spatial and demographic expansion into eastern Taiwan. This expansion began approximately 30,000 years ago in the Holocene interglacial period. The ancestral southern expansion into eastern Taiwan indicates that the central mountain range (CMR) formed a barrier to east-west expansion. However, E. formosa populations in the three major biogeographic regions (East, South, and North-Central) exhibit no significant genetic partitions, suggesting that river drainages and mountains did not form strong geographical barriers against gene flow among extant populations. CONCLUSIONS: The present study implies that the antiquity of E. formosa's colonization is associated with its high dispersal ability and larval tolerance to the late Pleistocene dry grasslands. The effect of late Pleistocene climatic changes on the subtropical damselfly's historical demography is lineage-specific, depending predominantly on its colonization history and geography. It is proposed that the Riss and Würm glaciations in the late Pleistocene period had a greater impact on the evolutionary diversification of subtropical insular species than the last glacial maximum (LGM).

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.

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.

Diversification in subtropical mountains: phylogeography, Pleistocene demographic expansion, and evolution of polyphenic mandibles in Taiwanese stag beetle, Lucanus formosanus.

Pleistocene glacial oscillations have had profound impacts on the historical population dynamics of extant species. However, the genetic consequences of past climatic changes depend largely on the latitude and topography of the regions in question. This study investigates the effect of Pleistocene glacial periods and the Central Mountain Range on the phylogeography, historical demography, and phenotypic differentiation of a montane forest-dwelling stag beetle, Lucanus formosanus (Coleoptera: Lucanidae), which exhibits extensive mandible variations across mountain ranges in subtropical Taiwan. Analyses of mitochondrial (cox1) and nuclear (wg) loci reveal that L. formosanus originated nearly 1.6 million years ago (Mya) in the early Pleistocene period and consisted of geographically overlapping Alishan and Widespread clades. A drastic population expansion starting approximately 0.2 Mya in the Widespread clade likely resulted from altitudinal range shift of the temperate forests, which was closely tied to the arrival of the Riss glacial period in the late Middle Pleistocene. A ring-like pattern of historical gene flow among neighboring populations in the vicinity of the Central Mountain Range indicates that the mountains constitute a strong vicariant barrier to the east-west gene flow of L. formosanus populations. A geographic cline of decreasing mandible size from central to north and south, and onto southeast of Taiwan is inconsistent with the low overall phylogeographic structures. The degree of mandible variation does not correlate with the expected pattern of neutral evolution, indicating that the evolutionary diversification of this morphological weapon is most likely subject to sexual or natural selection. We hypothesize that the adaptive evolution of mandibles in L. formosanus is shaped largely by the habitat heterogeneity.