Extensive hybridization and introgression between Melastoma candidum and M. sanguineum

Long non-coding RNA-mediated competing endogenous RNA regulatory network during flower development and color formation in Melastoma candidum

M. candidum, an evergreen shrubby flower known for its superior adaptation ability in South China, has gained increased attention in garden applications. However, scant attention has been paid to its flower development and color formation process at the non-coding RNA level. To fill this gap, we conducted a comprehensive analysis based on long non-coding RNA sequencing (lncRNA-seq), RNA-seq, small RNA sequencing (sRNA-seq), and widely targeted metabolome detection of three different flower developmental stages of M. candidum. After differentially expressed lncRNAs (DElncRNAs), differentially expressed mRNAs (DEmRNAs), differentially expressed microRNAs (DEmiRNAs), and differentially synthesized metabolites (DSmets) analyses between the different flower developmental stages, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted to identify some key genes and metabolites in flavonoid, flavone, anthocyanin, carotenoid, and alkaloid-related GO terms and biosynthetic pathways. Three direct-acting models, including antisense-acting, cis-acting, and trans-acting between lncRNAs and mRNAs, were detected to illustrate the direct function of lncRNAs on target genes during flower development and color formation. Based on the competitive endogenous RNA (ceRNA) regulatory theory, we constructed a lncRNA-mediated regulatory network composed of DElncRNAs, DEmiRNAs, DEmRNAs, and DSmets to elucidate the indirect role of lncRNAs in the flower development and color formation of M. candidum. By utilizing correlation analyses between DERNAs and DSmets within the ceRNA regulatory network, alongside verification trials of the ceRNA regulatory mechanism, the study successfully illustrated the significance of lncRNAs in flower development and color formation process. This research provides a foundation for improving and regulating flower color at the lncRNA level in M. candidum, and sheds light on the potential applications of non-coding RNA in studies of flower development.

Genomic evidence reveals high genetic diversity in a narrowly distributed species and natural hybridization risk with a widespread species in the genus Geodorum

BACKGROUND

Understanding genetic diversity is a core issue in conservation genetics. However, previous genetic diversity evaluations of narrowly distributed species have rarely used closely related widespread species as a reference. Furthermore, identifying natural hybridization signals between narrowly and widely distributed sympatric species is of great importance for the development of species conservation programs.

METHODS

In this study, population genotyping by sequencing (GBS) was performed for a narrowly distributed species, Geodorum eulophioides (endemic and endangered in Southwest China), and a widespread species, G. densiflorum. A total of 18,490 high-quality single nucleotide polymorphisms (SNPs) were identified at the whole-genome level.

RESULTS

The results showed that the nucleotide diversity and heterozygosity of G. eulophioides were significantly higher than those of G. densiflorum, confirming that narrowly distributed species can still preserve high genetic diversity. Consistent with taxonomic boundaries, all sampled individuals from the two species were divided into two genetic clusters and showed high genetic differentiation between species. However, in a sympatric population, a few G. eulophioides individuals were detected with genetic components from G. densiflorum, suggesting potential interspecific natural hybridization. This hypothesis was supported by Treemix analysis and hand-hybridization trials. Invasion of the habitat of G. eulophioides invasion by G. densiflorum under anthropogenic disturbance may be the main factor causing interspecific hybridization.

CONCLUSIONS

Therefore, reducing or avoiding habitat disturbance is a key measure to protect the G. eulophioides populations. This study provides valuable information for future conservation programs for narrowly distributed species.

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

Chromosomal-level genome assembly of Melastoma candidum provides insights into trichome evolution

Melastoma, consisting of ~100 species diversified in tropical Asia and Oceania in the past 1-2 million years, represents an excellent example of rapid speciation in flowering plants. Trichomes on hypanthia, twigs and leaves vary markedly among species of this genus and are the most important diagnostic traits for species identification. These traits also play critical roles in contributing to differential adaptation of these species to their own habitats. Here we sequenced the genome of M. candidum, a common, erect-growing species from southern China, with the aim to provide genomic insights into trichome evolution in this genus. We generated a high-quality, chromosome-level genome assembly of M. candidum, with the genome size of 256.2 Mb and protein-coding gene number of 40,938. The gene families specific to, and significantly expanded in Melastoma are enriched for GO terms related to trichome initiation and differentiation. We provide evidence that Melastoma and its sister genus Osbeckia have undergone two whole genome duplications (WGDs) after the triplication event (γ) shared by all core eudicots. Preferential retention of trichome development-related transcription factor genes such as C2H2, bHLH, HD-ZIP, WRKY, and MYB after both WGDs might provide raw materials for trichome evolution and thus contribute to rapid species diversification in Melastoma. Our study provides candidate transcription factor genes related to trichome evolution in Melastoma, which can be used to evolutionary and functional studies of trichome diversification among species of this genus.

Novel Approaches for Species Concepts and Delimitation in Polyploids and Hybrids

Hybridization and polyploidization are important processes for plant evolution. However, classification of hybrid or polyploid species has been notoriously difficult because of the complexity of processes and different evolutionary scenarios that do not fit with classical species concepts. Polyploid complexes are formed via combinations of allopolyploidy, autopolyploidy and homoploid hybridization with persisting sexual reproduction, resulting in many discrete lineages that have been classified as species. Polyploid complexes with facultative apomixis result in complicated net-work like clusters, or rarely in agamospecies. Various case studies illustrate the problems that apply to traditional species concepts to hybrids and polyploids. Conceptual progress can be made if lineage formation is accepted as an inevitable consequence of meiotic sex, which is established already in the first eukaryotes as a DNA restoration tool. The turnaround of the viewpoint that sex forms species as lineages helps to overcome traditional thinking of species as "units". Lineage formation and self-sustainability is the prerequisite for speciation and can also be applied to hybrids and polyploids. Species delimitation is aided by the improved recognition of lineages via various novel -omics methods, by understanding meiosis functions, and by recognizing functional phenotypes by considering morphological-physiological-ecological adaptations.

Flower-visiting insects of genus Melastoma (Myrtales: Melastomataceae) at the Fushan Botanical Garden, Taiwan

Background

We investigated the diversity and behaviour of insects that visit flowers of four native Melastoma (Family Melastomataceae) species of Taiwan and a horticultural hybrid Melastoma species at the Fushan Botanical Garden, Taiwan biweekly from May to August 2020. Visits of flower-visiting insects were classified into seven behavioural categories, based on the insects' behaviour and positions on the flower. The data are further assigned into four insect-flower interactions, namely pollination, herbivory, commensalism and neutralism. Our goal is to provide baseline data of insect-plant interactions of Melastoma, which is a common, but understudied plant genus in the country.

New information

A total of 1,289 visits to flowers were recorded by at least 63 insect morphospecies belonging to seven orders. The number of insect species recorded per Melastoma species ranged from 9 to 39. Visiting, sonication and passing were the three most frequently recorded types of behaviour, collectively accounting for 90.2% (n = 1,240) of the total observations. Pollination was the most dominant insect-flower interaction, accounting for 70.2% of the total observations, followed by neutralism (20.0%), herbivory (6.3%) and commensalism (3.5%). Sweat bees of the genera Lasioglossum and Maculonomia (Hymenoptera: Halictidae) are considered key pollinators to Melastoma species in Fushan Botanical Garden, based on their high number of visits and sonication behaviour. Our study provides the first list of insects that visit the flowers of all Taiwan's known Melastoma species and description of their interactions with the plants.

Lectotypification of the name Melastoma candidum f. albiflorum and its taxonomic status

A nomenclatural and taxonomic treatment of the name Melastoma candidum f. albiflorum (Melastomataceae) is presented. A lectotype is designated for this name, with an updated morphological description based on fresh material. The name Melastoma candidum f. albiflorum is proposed as a heterotypic synonym of Melastoma candidum.

Development and characterization of genomic microsatellite markers in the tree species, Rhodoleia championii, R. parvipetala, and R. forrestii (Hamamelidaceae)

Rhodoleia Champion ex Hooker is one of the most primitive relict genera of Hamamelidaceae, a key family exploited to understand the origin and early evolution of flowering plants. Genomic simple sequence repeats (SSRs) were developed for R. championii to perform genetic diversity, phylogeographical structure or even systematic evolution studies of the genus. Among the 278,743 contigs (105,758,242 bps) de novo assembled from the low-coverage whole genome sequencing of R. championii, a total of 9106 SSRs were detected in 8370 contigs, and SSR primer pairs were successfully designed for 6677 SSRs. Among the 110 selected primer pairs, 41 were amplified successfully in the preliminary test of SSR screening. Further amplification of these 41 primer pairs across the 122 individuals collected from six populations of the three Rhodoleia species showed that 32 and 40 SSR markers can be amplified in Vietnam and Jinping populations of R. parvipetala, 41, 33, and 41 SSR markers in Boluo, Hongkong and Xinyi populations of R. championii, 25 SSR markers in Fugong population of R. forrestii, and 20 SSR markers demonstrated to be polymorphic across the three species. Genetic analysis for these 20 polymorphic SSRs showed that Allele number (A) ranged from four to 13 and polymorphic information content (PIC) ranged from 0.479 to 0.876 across the three species. At the population level, observed heterozygosity (H_O) ranged from 0.000 to 1.000, and expected heterozygosity (H_E) ranged from 0.091 to 0.851. In the present study, we provided the first whole-genome sequencing database for the species R. championii, identified ample SSR loci with designed primers, and revealed that 20 of the 110 selected SSRs were polymorphic across three Rhodoleia species. These provide valuable resources for future studies on genetic study, species delimitation, phylogeography, and conservation of this genus.

De novo Assembly and Characterization of the Floral Transcriptomes of Two Varieties of Melastoma malabathricum

Melastoma malabathricum is an important medicinal and landscape plant that is globally distributed in temperate and subtropical regions. However, available genomic information for the entire Melastomataceae family is notably limited. In view of the application potential of floral parts in secondary metabolite extraction, we characterized for the first time the floral transcriptomes of two key M. malabathricum varieties, purple variety and white variety. Our transcriptome assembly generated 52,498 and 49,380 unigenes with an N50 of 1,906 and 1,929 bases for the purple and white varieties, respectively. Comparative analysis of two transcriptomes demonstrated that they are highly similar but also highlighted genes that are presumably lineage specific, which explains the phenotypes of each variety. Additionally, a shared transcriptional signature across the floral developmental stages was identified in both M. malabathricum varieties; this signature included pathways related to secondary metabolite synthesis, plant hormone signaling and production, energy homeostasis and nutrient assimilation pathways, and cellular proliferation. The expression levels of flavonoid accumulation and candidate flavonoid biosynthesis-related genes in M. malabathricum flower development stages validated the transcriptome findings. The transcriptome data presented in this study will serve as a valuable resource for future work on the exploitation of M. malabathricum and other related species. The gene expression dynamics during flower development will facilitate the discovery of lineage-specific genes associated with phenotypic characteristics and will elucidate the mechanism of the ontogeny of individual flower types.

Molecular identification of natural hybridization between Melastoma malabathricum and Melastoama beccarianum in Sarawak, Malaysia

Hybridization is very common in flowering plants and it plays a significant role in plant evolution and adaptation. Melastoma L. (Melastomataceae) comprises about 80-90 species in tropical Asia and Oceania, among which 41 species occur in Borneo. Natural hybridization is frequently reported in Melastoma in China, but so far there have been no confirmed cases of hybridization in Southeast Asia (including Borneo), where most species occur. Here, we identified a case of natural hybridization between Melastoma malabathricum L. and Melastoma beccarianum Cogn. in Sarawak, Malaysia, by using sequence data of three nuclear genes and one chloroplast intergenic spacer. Melastoma malabathricum is the most widespread species of this genus, occurring in almost the whole range of this genus, while M. beccarianum is a local species endemic to northern Borneo. Our results showed that natural hybridization and introgression occur between M. malabathricum and M. beccarianum, and the introgression was asymmetrical, mainly from M. malabathricum to M. beccarianum. As adaptive traits can be transferred by introgression, our study suggests that natural hybridization should be a significant mechanism for the evolution and adaptation of Melastoma in Southeast Asia. However, introgression from the common species M. malabathricum to the relatively rare species M. beccarianum may cause the decline of M. beccarianum, incurring conservation concern. With a large number of species of Melastoma and almost year-around flowering in Southeast Asia, more cases of natural hybridization are expected to be found and identified in near future.

Comparative transcriptomics sheds light on differential adaptation and species diversification between two Melastoma species and their F1 hybrid

Variation in gene expression has been shown to promote adaptive divergence, and can lead to speciation. The plant genus Melastoma, thought to have diversified through adaptive radiation, provides an excellent model for the study of gene expressional changes during adaptive differentiation and following interspecific hybridization. In this study, we performed RNA-seq on M. candidum, M. sanguineum and their F1 hybrid, to investigate the role of gene expression in species diversification within the genus. Reference transcriptomes were assembled using combined data from both parental species, resulting in 50 519 and 48 120 transcripts for the leaf and flower petal, after removing redundancy. Differential expression analysis uncovered 3793 and 2116 differentially expressed (DE) transcripts, most of which are between M. candidum and M. sanguineum. Differential expression was observed for genes related to light responses, as well as genes that regulate the development of leaf trichomes, a trait that among others is thought to protect plants against sunlight, suggesting the differential adaptation of the species to sunlight intensity. The analysis of positively selected genes between the two species also revealed possible differential adaptation to other abiotic stresses such as drought and temperature. In the hybrid, almost all possible modes of expression were observed at the DE transcripts, although at most transcripts, the expression levels were similar to that of either parent instead of being intermediate. A small number of transgressively expressed transcripts that matched genes known to promote plant growth and adaptation to stresses in new environments were also found, possibly explaining the vigour observed in the hybrid. The findings in this study provided insights into the role of gene expression in the diversification of Melastoma, which we believe is an important example for more cross-taxa comparisons in the future.

The specific status of Melastoma kudoi (Melastomataceae, Melastomeae)

BACKGROUND

Melastoma has undergone rapid species radiation during the last one million years, and circumscription of some species in the genus has remained controversial. Melastoma kudoi, an erect species narrowly endemic to central Taiwan was previously treated as a synonym of M. intermedium, a semicreeping hybrid between the erect species M. candidum and the creeping M. dodecandrum, making its identity questionable. We addressed this question based on molecular and morphological data.

RESULTS

Phylogenetic analyses based on nrITS sequence data revealed that M. kudoi is most closely related to M. dodecandrum. Further analyses of six nuclear genes (cam, chi, gapC, gbss, tpi and vr) and two chloroplast markers (trnL-trnF and psbA) showed that M. kudoi is well diverged from its close relatives. Morphologically, it is also easily distinguished from related species by its erect habit, center-positioned stigma, and spreading, basally enlarged hairs on the hypanthium.

CONCLUSIONS

Both molecular phylogenetic and morphological data suggest that M. kudoi is well separated from M. intermedium, M. dodecandrum, and O. scaberrima, and should be treated as a distinct species. Taxonomic treatment and detailed description of M. kudoi are provided.

Microwave-Assisted Extraction of Phenolic Compounds from Melastoma sanguineum Fruit: Optimization and Identification

A microwave-assisted extraction (MAE) technology optimized by response surface methodology (RSM) was established to extract phenolic compounds from the fruit of Melastoma sanguineum. The effects of solvent composition, ratio of solvent to material, temperature, time and microwave power on phenol yield were evaluated in single-factor tests. The three parameters exerting main impacts on phenol yield were further optimized by RSM. Under optimal extraction conditions (31.33% ethanol, solvent/material ratio of 32.21 mL/g, 52.24 °C, 45 min and 500 W), the total phenolic content was 39.02 ± 0.73 mg gallic acid equivalent (GAE)/g dry weight (DW). This MAE method performed better in comparison with two conventional methods, those being maceration (25.79 ± 1.03 mg GAE/g DW) and Soxhlet extraction (18.40 ± 1.34 mg GAE/g DW), using lower process temperature, shorter irradiation time, and lower organic solvent consumption. In addition, five flavonoids (epicatechin gallate, epicatechin, rutin, pigallocatechin and quercetin) and two phenolic acids (protocatechuic acid and chlorogenic acid) in the extract were identified and quantified using UPLC-MS/MS.

Molecular evidence of hybridization in sympatric populations of the Enantia jethys complex (Lepidoptera: Pieridae)

Hybridization events are frequently demonstrated in natural butterfly populations. One interesting butterfly complex species is the Enantia jethys complex that has been studied for over a century; many debates exist regarding the species composition of this complex. Currently, three species that live sympatrically in the Gulf slope of Mexico (Enantia jethys, E. mazai, and E. albania) are recognized in this complex (based on morphological and molecular studies). Where these species live in sympatry, some cases of interspecific mating have been observed, suggesting hybridization events. Considering this, we employed a multilocus approach (analyses of mitochondrial and nuclear sequences: COI, RpS5, and Wg; and nuclear dominant markers: inter-simple sequence repeat (ISSRs) to study hybridization in sympatric populations from Veracruz, Mexico. Genetic diversity parameters were determined for all molecular markers, and species identification was assessed by different methods such as analyses of molecular variance (AMOVA), clustering, principal coordinate analysis (PCoA), gene flow, and PhiPT parameters. ISSR molecular markers were used for a more profound study of hybridization process. Although species of the Enantia jethys complex have a low dispersal capacity, we observed high genetic diversity, probably reflecting a high density of individuals locally. ISSR markers provided evidence of a contemporary hybridization process, detecting a high number of hybrids (from 17% to 53%) with significant differences in genetic diversity. Furthermore, a directional pattern of hybridization was observed from E. albania to other species. Phylogenetic study through DNA sequencing confirmed the existence of three clades corresponding to the three species previously recognized by morphological and molecular studies. This study underlines the importance of assessing hybridization in evolutionary studies, by tracing the lineage separation process that leads to the origin of new species. Our research demonstrates that hybridization processes have a high occurrence in natural populations.

Evidence for continual hybridization rather than hybrid speciation between Ligularia duciformis and L. paradoxa (Asteraceae)

Background

Hybrids possess phenotypic traits that are often intermediate between their parental taxa, which commonly serves as evidence of hybridization in morphological analyses. Natural hybridization has been shown to occur frequently in Ligularia (Asteraceae). In a previous study, Ligularia ×maoniushanensis was demonstrated as a natural hybrid species between L. duciformis and L. paradoxa based on morphological and reproductive traits.

Methods

We used three chloroplast (cpDNA) fragments (psbA-trnH, trnL-rpl32 and trnQ-5′rps16), the nuclear ribosomal internal transcribed spacer (nrITS), and co-dominant SSR and dominant ISSR markers to study natural hybridization between L. duciformis and L. paradoxa growing sympatrically in two locations. Parental taxa were inferred using network analyses of cpDNA and nrITS haplotypes. Admixture among individuals was examined using the Bayesian clustering programs STRUCTURE and NewHybrids based on the SSR and ISSR data; and potential introgression in the SSR loci was assessed using the INTROGRESS package.

Results

The putative parental species were clearly distinguished from other sympatric Ligularia species by nrITS data, and L. ×maoniushanensis individuals were confirmed to be the hybrid offspring of L. duciformis and L. paradoxa. Moreover, introgression was detected among several individuals morphologically identified as L. duciformis or L. paradoxa. Analyses of the cpDNA data revealed primarily unidirectional hybridization between L. duciformis and L. paradoxa, with L. paradoxa as the maternal parent in Mt. Maoniu, whereas bidirectional but asymmetrical hybridization was inferred to occur in Heihai Lake. The STRUCTURE analyses based on the SSR data detected two distinct clusters among the three taxa. The NewHybrids analyses showed that individuals circumscribed as L. ×maoniushanensis were dominated by early- and later-generation and backcrossing hybrids. The NewHybrids results based on the ISSR data were congruent with SSR results. In addition, introgression was detected in some SSR loci, and heterogeneity among loci was found in terms of detected patterns of introgression.

Conclusions

Our data provide strong evidence for hybridization and introgression between L. duciformis and L. paradoxa. Ligularia ×maoniushanensis was demonstrated to be of hybrid origin. Since no evident reproductive isolation was found between the two parental species, detected hybrids appear to be part of hybrid swarms resulting from frequent and ongoing gene flow, which might impede the formation of a new hybrid species.

The complete chloroplast genome sequence of Melastoma candidum (Melastomataceae)

The plant genus Melastoma is comprised of members estimated to have formed through recent species radiation. Natural hybridization among member species further complicates taxonomy within the genus. Herein, we report the complete chloroplast genome of M. candidum, assembled from partial data obtained from a parallel whole-genome Illumina paired-end sequencing effort on the species. The chloroplast genome was 156,682 bp in length, with a large single-copy (LSC) region of 86,084 bp, a small single-copy (SSC) region of 17,094 bp, separated by two inverted repeat (IR) regions of 26,752 bp each. It was predicted to contain a total of 129 genes, with an overall GC content of 37.17%. Phylogenetic analysis placed M. candidum in the same clade as species within the Melastomeae tribe of Melastomataceae.

Ultrasound-Assisted Extraction and Identification of Natural Antioxidants from the Fruit of Melastoma sanguineum Sims

The fruit of Melastoma sanguineum Sims is an edible and sweet wild fruit. In our previous study, the fruit was found to have a strong antioxidant property. In this study, an ultrasound-assisted extraction (UAE) method was developed to extract natural antioxidants from the fruit of Melastoma sanguineum Sims, and a response surface methodology was used to optimize the conditions of UAE to maximize the extraction efficiency. The influence of five independent extraction parameters (ethanol concentration, solvent/material ratio, extracting time, temperature, and ultrasound power) on the extraction efficiency were investigated using a single factor experiment, and then a central composite rotatable design was used to investigate the interaction of three key parameters. The results showed that the optimal extraction conditions were 42.98% ethanol, 28.29 mL/g solvent/material ratio, 34.29 min extracting time, 60 °C temperature, and 600 W ultrasound power. Under these conditions, the Trolox equivalent antioxidant capacity (TEAC) value of the extracts was 1074.61 ± 32.56 μmol Trolox/g dry weight (DW). Compared with conventional maceration (723.27 ± 11.61 μmol Trolox/g DW) and Soxhlet extraction methods (518.37 ± 23.23 μmol Trolox/g DW), the UAE method improved the extraction efficiency, in a shorter period of time. In addition, epicatechin gallate, epicatechin, rutin, epigallocatechin, protocatechuic acid, chlorogenic acid, and quercetin, were identified and quantified in the fruit extracts of Melastoma sanguineum Sims by UPLC-MS/MS.

Similar Morphologies but Different Origins: Hybrid Status of Two More Semi-creeping Taxa of Melastoma

Inferring the origins of hybrid taxa based on morphology alone is difficult because morphologically similar hybrids can arise from hybridization between different populations of the same parental species or be produced by hybridization of different parental species. In this study, we investigated the origins of two semi-creeping taxa in Melastoma, which are morphologically similar to a natural hybrid, M. intermedium, by sequencing a chloroplast intergenic spacer, nuclear ribosomal internal transcribed spacer and two low-copy nuclear genes (tpi and cam) in these taxa and their putative parental species. Our sequence analysis provides compelling evidence for the hybrid status of the two semi-creeping taxa: one originating from hybridization between M. dodecandrum and M. malabathricum, and the other between M. dodecandrum and M. normale. The origins of these hybrids are therefore clearly different from M. intermedium, and morphological similarity for the three hybrids is most likely due to their origins from hybridization between the same creeping species M. dodecandrum and a different erect species in each of the three cases. We also observed low rate of introgression from M. normale to M. dodecandrum, and genetic exchange between them may transfer adaptive traits to M. dodecandrum. Rare occurrence of these two hybrids may be due to small range overlaps between parental species in one case, and different flowering periods between parental species in the other.

The complex evolutionary history of big-eared horseshoe bats (Rhinolophus macrotis complex): insights from genetic, morphological and acoustic data

Palaeoclimatic oscillations and different landscapes frequently result in complex population-level structure or the evolution of cryptic species. Elucidating the potential mechanisms is vital to understanding speciation events. However, such complex evolutionary patterns have rarely been reported in bats. In China, the Rhinolophus macrotis complex contains a large form and a small form, suggesting the existence of a cryptic bat species. Our field surveys found these two sibling species have a continuous and widespread distribution with partial sympatry. However, their evolutionary history has received little attention. Here, we used extensive sampling, morphological and acoustic data, as well as different genetic markers to investigate their evolutionary history. Genetic analyses revealed discordance between the mitochondrial and nuclear data. Mitochondrial data identified three reciprocally monophyletic lineages: one representing all small forms from Southwest China, and the other two containing all large forms from Central and Southeast China, respectively. The large form showed paraphyly with respect to the small form. However, clustering analyses of microsatellite and Chd1 gene sequences support two divergent clusters separating the large form and the small form. Moreover, morphological and acoustic analyses were consistent with nuclear data. This unusual pattern in the R. macrotis complex might be accounted for by palaeoclimatic oscillations, shared ancestral polymorphism and/or interspecific hybridization.