Short Divergence Time Research Articles

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

The complex Y-chromosomal history of gorillas.

Studies of the evolutionary relationships among gorilla populations using autosomal and mitochondrial sequences suggest that male-mediated gene flow may have been important in the past, but data on the Y-chromosomal relationships among the gorilla subspecies are limited. Here, we genotyped blood and noninvasively collected fecal samples from 12 captives and 257 wild male gorillas of known origin representing all four subspecies (Gorilla gorilla gorilla, G. g. diehli, G. beringei beringei, and G. b. graueri) at 10 Y-linked microsatellite loci resulting in 102 unique Y-haplotypes for 224 individuals. We found that western lowland gorilla (G. g. gorilla) haplotypes were consistently more diverse than any other subspecies for all measures of diversity and comprised several genetically distinct groups. However, these did not correspond to geographical proximity and some closely related haplotypes were found several hundred kilometers apart. Similarly, our broad sampling of eastern gorillas revealed that mountain (G. b. beringei) and Grauer's (G. b. graueri) gorilla Y-chromosomal haplotypes did not form distinct clusters. These observations suggest structure in the ancestral population with subsequent mixing of differentiated haplotypes by male dispersal for western lowland gorillas, and postisolation migration or incomplete lineage sorting due to short divergence times for eastern gorillas.

Morphological, Ecological, and Molecular Divergence of Conogethes pinicolalis from C. punctiferalis (Lepidoptera: Crambidae).

Simple SummaryCurrently, two species of Conogethes have been documented in Korea: C. punctiferalis and C. pinicolalis. However, as C. pinicolalis has long been considered as a Pinaceae-feeding type of the yellow peach moth, C. punctiferalis, studies on C. pinicolalis are limited. Therefore, in this study, the divergence of C. pinicolalis from the fruit-feeding moth C. punctiferalis was analyzed in terms of morphology, ecology, and genetics. C. pinicolalis differs from C. punctiferalis in several morphological features. Through field observation, we confirmed that pine trees are the host plants for the first generation of C. pinicolalis larvae. We successfully reared C. pinicolalis larvae to adults by providing them pine needles as a diet. Sequences of mitochondrial COI of these two species substantially diverged by an average of 5.46%. Overall nuclear EF1α-based phylogeny confirmed each species as an independent clade, but a few haplotypes of EF1α indicated incomplete lineage sorting between these two species, suggesting a short divergence time.Conogethes pinicolalis has long been considered as a Pinaceae-feeding type of the yellow peach moth, C. punctiferalis, in Korea. In this study, the divergence of C. pinicolalis from the fruit-feeding moth C. punctiferalis was analyzed in terms of morphology, ecology, and genetics. C. pinicolalis differs from C. punctiferalis in several morphological features. Through field observation, we confirmed that pine trees are the host plants for the first generation of C. pinicolalis larvae, in contrast to fruit-feeding C. punctiferalis larvae. We successfully reared C. pinicolalis larvae to adults by providing them pine needles as a diet. From a genetic perspective, the sequences of mitochondrial COI of these two species substantially diverged by an average of 5.46%; moreover, phylogenetic analysis clearly assigned each species to an independent clade. On the other hand, nuclear EF1α showed a lower sequence divergence (2.10%) than COI. Overall, EF1α-based phylogenetic analysis confirmed each species as an independent clade, but a few haplotypes of EF1α indicated incomplete lineage sorting between these two species. In conclusion, our results demonstrate that C. pinicolalis is an independent species according to general taxonomic criteria; however, analysis of the EF1α sequence revealed a short divergence time.

Confronting Sources of Systematic Error to Resolve Historically Contentious Relationships: A Case Study Using Gadiform Fishes (Teleostei, Paracanthopterygii, Gadiformes).

Reliable estimation of phylogeny is central to avoid inaccuracy in downstream macroevolutionary inferences. However, limitations exist in the implementation of concatenated and summary coalescent approaches, and Bayesian and full coalescent inference methods may not yet be feasible for computation of phylogeny using complicated models and large data sets. Here, we explored methodological (e.g., optimality criteria, character sampling, model selection) and biological (e.g., heterotachy, branch length heterogeneity) sources of systematic error that can result in biased or incorrect parameter estimates when reconstructing phylogeny by using the gadiform fishes as a model clade. Gadiformes include some of the most economically important fishes in the world (e.g., Cods, Hakes, and Rattails). Despite many attempts, a robust higher-level phylogenetic framework was lacking due to limited character and taxonomic sampling, particularly from several species-poor families that have been recalcitrant to phylogenetic placement. We compiled the first phylogenomic data set, including 14,208 loci (n}{}>2.8 M bp) from 58 species representing all recognized gadiform families, to infer a time-calibrated phylogeny for the group. Data were generated with a gene-capture approach targeting coding DNA sequences from single-copy protein-coding genes. Species-tree and concatenated maximum-likelihood (ML) analyses resolved all family-level relationships within Gadiformes. While there were a few differences between topologies produced by the DNA and the amino acid data sets, most of the historically unresolved relationships among gadiform lineages were consistently well resolved with high support in our analyses regardless of the methodological and biological approaches used. However, at deeper levels, we observed inconsistency in branch support estimates between bootstrap and gene and site coefficient factors (gCF, sCF). Despite numerous short internodes, all relationships received unequivocal bootstrap support while gCF and sCF had very little support, reflecting hidden conflict across loci. Most of the gene-tree and species-tree discordance in our study is a result of short divergence times, and consequent lack of informative characters at deep levels, rather than incomplete lineage sorting. We use this phylogeny to establish a new higher-level classification of Gadiformes as a way of clarifying the evolutionary diversification of the order. We recognize 17 families in five suborders: Bregmacerotoidei, Gadoidei, Ranicipitoidei, Merluccioidei, and Macrouroidei (including two subclades). A time-calibrated analysis using 15 fossil taxa suggests that Gadiformes evolved n}{}sim 79.5 Ma in the late Cretaceous, but that most extant lineages diverged after the Cretaceous–Paleogene (K-Pg) mass extinction (66 Ma). Our results reiterate the importance of examining phylogenomic analyses for evidence of systematic error that can emerge as a result of unsuitable modeling of biological factors and/or methodological issues, even when data sets are large and yield high support for phylogenetic relationships. [Branch length heterogeneity; Codfishes; commercial fish species; Cretaceous-Paleogene (K-Pg); heterotachy; systematic error; target enrichment.]

GISH-based comparative genomic analysis in Urochloa P. Beauv.

The genus Urochloa P. Beauv. [syn. Brachiaria (Trin.) Griseb.] comprises species of great economic relevance as forages. The genomic constitution for the allotetraploid species Urochloa brizantha (cv. Marandu) and Urochloa decumbens (cv. Basilisk) and the diploid Urochloa ruziziensis was previously proposed as BBB1B1, B1B1B2B2 and B2B2, respectively. Evidence indicates U. ruziziensis as the ancestral donor of genome B2 in U. decumbens allotetraploidy, but the origin of the genomes B and B1 is still unknown. There are diploid genotypes of U. brizantha and U. decumbens that may be potential ancestors of the tetraploids. The aim of this study was to determine the genomic constitution and relationships between genotypes of U. brizantha (2x and 4x), U. decumbens (2x and 4x) and U. ruziziensis (2x) via genomic in situ hybridization (GISH). Additionally, chromosome number and genome size were verified for the diploid genotypes. The diploids U. brizantha and U. decumbens presented 2n = 2x = 18 chromosomes and DNA content of 1.79 and 1.44pg, respectively. The GISH analysis revealed high homology between the diploids U. brizantha and U. decumbens, which suggests relatively short divergence time. The GISH using genomic probes from the diploid accessions on the tetraploid accessions' chromosomes presented similar patterns, highlighting the genome B1 present in both of the tetraploids. Based on GISH results, the genomic constitution was proposed for the diploid genotypes of U. brizantha (B1B1) and U. decumbens (B1'B1') and both were pointed as donors of genome B1 (or B1'), present in the allotetraploid genotypes.

Craniofacial Shape and Nonmetric Trait Variation in Hybrids of the Japanese Macaque (Macaca fuscata) and the Taiwanese Macaque (Macaca cyclopis)

It has become apparent that natural hybridization is far more common and may play a much greater role in evolution than has historically been recognized. The skeletal morphology of hybrid primates is notoriously variable and difficult to predict. Indeed, before the advent of genetic sequencing techniques, many wild hybrid populations went undetected. Though many species of primates are now known to hybridize naturally and are likely to have done so for millions of years, anthropogenic alterations to the environment are increasingly restricting or altering primate species ranges and contact zones and driving hybridization between populations that may otherwise never have come into contact. The case of hybridizing Japanese and Taiwanese macaques (Macaca fuscata and Macaca cyclopis) is an excellent example of this, as these two island species could not have come into contact without human interference. Here we apply 3D geometric morphometrics and nonmetric trait analysis to the crania and dentition of hybrid macaques (N = 70) and their parental species, M. fuscata (N = 57) and M. cyclopis (N = 51). The exploration of 3D shape variation identifies mildly transgressive morphology in the hybrids and a general tendency toward the M. fuscata morphotype overall, but less variability in the hybrid morphotype than has been identified in previous studies of primate hybrids. We also identify a small number of nonmetric traits that differentiate the hybrids from the parental species, although the power of these traits to distinguish between groups is weak and their relationship with hybridity is unclear. We conclude that the relatively short divergence time between the parent species is likely to help explain the observed differences in hybrid morphotype, and that further exploration of the relationship between degree of evolutionary divergence and hybrid morphology may help us to better explain and predict hybrid morphology in other taxa.

Intraspecific brain size variation between coexisting sunfish ecotypes.

Variation in spatial complexity and foraging requirements between habitats can impose different cognitive demands on animals that may influence brain size. However, the relationship between ecologically related cognitive performance and brain size is not well established. We test whether variation in relative brain size and brain region size is associated with habitat use within a population of pumpkinseed sunfish composed of different ecotypes that inhabit either the structurally complex shoreline littoral habitat or simpler open-water pelagic habitat. Sunfish using the littoral habitat have on average 8.3% larger brains than those using the pelagic habitat. We found little difference in the proportional sizes of five brain regions between ecotypes. The results suggest that cognitive demands on sunfish may be reduced in the pelagic habitat given no habitat-specific differences in body condition. They also suggest that either a short divergence time or physiological processes may constrain changes to concerted, global modifications of brain size between sunfish ecotypes.

CLADES: A classification-based machine learning method for species delimitation from population genetic data.

Species are considered to be the basic unit of ecological and evolutionary studies. As multilocus genomic data are increasingly available, there have been considerable interests in the use of DNA sequence data to delimit species. In this study, we show that machine learning can be used for species delimitation. Our method treats the species delimitation problem as a classification problem for identifying the category of a new observation on the basis of training data. Extensive simulation is first conducted over a broad range of evolutionary parameters for training purposes. Each pair of known populations is combined to form training samples with a label of "same species" or "different species". We use support vector machine (SVM) to train a classifier using a set of summary statistics computed from training samples as features. The trained classifier can classify a test sample to two outcomes: "same species" or "different species". Given multilocus genomic data of multiple related organisms or populations, our method (called CLADES) performs species delimitation by first classifying pairs of populations. CLADES then delimits species by maximizing the likelihood of species assignment for multiple populations. CLADES is evaluated through extensive simulation and also tested on real genetic data. We show that CLADES is both accurate and efficient for species delimitation when compared with existing methods. CLADES can be useful especially when existing methods have difficulty in delimitation, for example with short species divergence time and gene flow.

Swimmers, Diggers, Climbers and More, a Study of Integration Across the Mustelids’ Locomotor Apparatus (Carnivora: Mustelidae)

Phenotypic integration, defined as the coordinated co-variation of parts of an organism can be an important constraint on phenotypic diversification. Functional factors, by having an heterogeneous impact across the animal body, may reinforce the integration of some parts while causing a perturbation of the integration among other parts. The integration across the locomotor apparatus should thus reflect to a certain extent the locomotor ecology of the animal. Using the mustelids as study group, we track changes in the patterns of co-variation in species belonging to four different locomotor ecologies (terrestrial, semi-arboreal, semi-fossorial, and semi-aquatic). Our results highlight the strong overall integration in mustelid long bones. The main shape changes associated with co-variations between skeletal elements are the bone robustness and proportions of the epiphyses. The pattern of co-variation is, however, only slightly impacted by allometry. Changes in co-variation between species mostly scale with phylogenetic divergence time, except for the (Mustela putorius, M. eversmanni, M. lutreola) clade which, despite a short divergence time, presents strong differences in co-variation. Co-variation patterns differ between locomotor ecologies, but few of these variations match the hypothesis of a reduction of integration due to functional specialization. This may reflect our lack of knowledge on the functional modules in species with locomotor ecologies that differ from terrestrial locomotion rather than invalidate our a priori hypotheses.

Multilocus DNA barcoding \u2013 Species Identification with Multilocus Data

Species identification using DNA sequences, known as DNA barcoding has been widely used in many applied fields. Current barcoding methods are usually based on a single mitochondrial locus, such as cytochrome c oxidase subunit I (COI). This type of barcoding method does not always work when applied to species separated by short divergence times or that contain introgressed genes from closely related species. Herein we introduce a more effective multi-locus barcoding framework that is based on gene capture and “next-generation” sequencing. We selected 500 independent nuclear markers for ray-finned fishes and designed a three-step pipeline for multilocus DNA barcoding. We applied our method on two exemplar datasets each containing a pair of sister fish species: Siniperca chuatsi vs. Sini. kneri and Sicydium altum vs. Sicy. adelum, where the COI barcoding approach failed. Both of our empirical and simulated results demonstrated that under limited gene flow and enough separation time, we could correctly identify species using multilocus barcoding method. We anticipate that, as the cost of DNA sequencing continues to fall that our multilocus barcoding approach will eclipse existing single-locus DNA barcoding methods as a means to better understand the diversity of the living world.

Integrated species delimitation and conservation implications of an endangered weevil Pachyrhynchus sonani (Coleoptera: Curculionidae) in Green and Orchid Islands of Taiwan

AbstractOceanic islands are productive habitats for generating new species and high endemism, which is primarily due to their geographical isolation, smaller population sizes and local adaptation. However, the short divergence times and subtle morphological or ecological divergence of insular organisms may obscure species identity, so the cryptic endemism on islands may be underestimated. The endangered weevil Pachyrhynchus sonani Kôno (Coleoptera: Curculionidae: Entiminae: Pachyrhynchini) is endemic to Green Island and Orchid Island of the Taiwan‐Luzon Archipelago and displays widespread variation in coloration and host range, thus raising questions regarding its species boundaries and degree of cryptic diversity. We tested the species boundaries of P. sonani using an integrated approach that combined morphological (body size and shape, genital shape, coloration and cuticular scale), genetic (four genes and restriction site‐associated DNA sequencing, RAD‐seq) and ecological (host range and distribution) diversity. The results indicated that all the morphological datasets for male P. sonani, except for the colour spectrum, reveal overlapping but statistically significant differences between islands. In contrast, the morphology of the female P. sonani showed minimum divergence between island populations. The populations of P. sonani on the two islands were significantly different in their host ranges, and the genetic clustering and phylogenies of P. sonani established two valid evolutionary species. Integrated species delimitation combining morphological, molecular and ecological characters supported two distinct species of P. sonani from Green Island and Orchid Island. The Green Island population was described as P. jitanasaius sp.n. Chen & Lin, and it is recommended that its threatened conservation status be recognized. Our findings suggest that the inter‐island speciation of endemic organisms inhabiting both islands may be more common than previously thought, and they highlight the possibility that the cryptic diversity of small oceanic islands may still be largely underestimated.

Fact or fiction: updates on how protein-coding genes might emerge de novo from previously non-coding DNA.

Over the last few years, there has been an increasing amount of evidence for the de novo emergence of protein-coding genes, i.e. out of non-coding DNA. Here, we review the current literature and summarize the state of the field. We focus specifically on open questions and challenges in the study of de novo protein-coding genes such as the identification and verification of de novo-emerged genes. The greatest obstacle to date is the lack of high-quality genomic data with very short divergence times which could help precisely pin down the location of origin of a de novo gene. We conclude that, while there is plenty of evidence from a genetics perspective, there is a lack of functional studies of bona fide de novo genes and almost no knowledge about protein structures and how they come about during the emergence of de novo protein-coding genes. We suggest that future studies should concentrate on the functional and structural characterization of de novo protein-coding genes as well as the detailed study of the emergence of functional de novo protein-coding genes.

Hierarchical Map of Orthologous Genomic Regions Reconstructed from Two Closely Related Genomes: Cucumber Case Study.

Accurate identification of orthologous genomic regions (OGRs) between two closely related genomes is crucial for the reliable detection of genomic changes, which range from small-scale changes (e.g., single nucleotide or small nucleotides) to large-scale structural changes. Although diverse OGRs inferred at different levels have been successfully applied to address various biological questions, a limited number of studies have simultaneously integrated OGRs from different levels. Here, we report on a new approach to construct a hierarchical map of OGRs. Using different types of genomic markers, this approach was applied to two very closely related cucumber genomes [ L. var. and L. var. (Royle) Alef.]. We identified two different levels of OGRs using Mugsy (denoted as dnaOGRs) and i-ADHoRe (denoted as proOGRs). Using information regarding the anchored chromosomes of the two genomes, a third level of OGRs (denoted chrOGRs) could be built at the chromosomal level. Together, these OGRs could be organized into a hierarchical map that represented the parent-child relationships (chrOGR:proOGRs:dnaOGRs) between the two genomes. For this case study, the map consisted of seven chrOGRs, 540 proOGRs, and 22,321 dnaOGRs. Based on this map, we designed different methods to detect both small-scale and large-scale genomic changes. Surprisingly, many genomic changes were detected at each OGR level despite the very short divergence time between the two subspecies. Together, our results show that a hierarchical map of OGRs and their related genomic changes are useful resources for elucidating the diversity and evolution of cucumber genomes and phenotypes.

Identification of cichlid fishes from Lake Malawi using computer vision.

BackgroundThe explosively radiating evolution of cichlid fishes of Lake Malawi has yielded an amazing number of haplochromine species estimated as many as 500 to 800 with a surprising degree of diversity not only in color and stripe pattern but also in the shape of jaw and body among them. As these morphological diversities have been a central subject of adaptive speciation and taxonomic classification, such high diversity could serve as a foundation for automation of species identification of cichlids.Methodology/Principal FindingHere we demonstrate a method for automatic classification of the Lake Malawi cichlids based on computer vision and geometric morphometrics. For this end we developed a pipeline that integrates multiple image processing tools to automatically extract informative features of color and stripe patterns from a large set of photographic images of wild cichlids. The extracted information was evaluated by statistical classifiers Support Vector Machine and Random Forests. Both classifiers performed better when body shape information was added to the feature of color and stripe. Besides the coloration and stripe pattern, body shape variables boosted the accuracy of classification by about 10%. The programs were able to classify 594 live cichlid individuals belonging to 12 different classes (species and sexes) with an average accuracy of 78%, contrasting to a mere 42% success rate by human eyes. The variables that contributed most to the accuracy were body height and the hue of the most frequent color.ConclusionsComputer vision showed a notable performance in extracting information from the color and stripe patterns of Lake Malawi cichlids although the information was not enough for errorless species identification. Our results indicate that there appears an unavoidable difficulty in automatic species identification of cichlid fishes, which may arise from short divergence times and gene flow between closely related species.

Application of molecular genetics and geometric morphometrics to taxonomy and conservation of cave beetles in central Italy

Integration of molecular genetic techniques and geometric morphometrics represent a valuable tool in the resolution of taxonomic uncertainty and the identification of significant units for conservation. We combined mitochondrial DNA cytochrome c oxidase subunit II gene sequence data and geometric morphometric analysis to examine taxonomic status and identify units for conservation in four species of the hypogean beetle Duvalius (Coleoptera, Trechinae) using mainly museum specimens collected in central Italy. Previous taxonomic studies based on morphological traits described several subspecies often inhabiting geographically distinct caves. Phylogenetic analysis identified two well supported monophyletic lineages and a number of different clades with relatively small genetic differences, suggesting a short divergence time in line with known geological history of the study area. Geometric morphometrics, on the other hand, recovered a high level of distinctiveness among specimens. Both genetic and morphometric analyses did not entirely corroborate former taxonomic nomenclature, suggesting possible rearrangements and the definition of evolutionary significant units. Beetles of the genus Duvalius are protected by regional laws and the majority of taxa considered in this study inhabit caves located outside protected areas. Our study advocates the importance of devoting protection efforts to networks of cave ecosystems rather than single locations or species.

Molecular phylogeny of extant equids and effects of ancestral polymorphism in resolving species-level phylogenies

Short divergence times and processes such as incomplete lineage sorting and species hybridization are known to hinder the inference of species-level phylogenies due to the lack of sufficient informative genetic variation or the presence of shared but incongruent polymorphism among taxa. Extant equids (horses, zebras, and asses) are an example of a recently evolved group of mammals with an unresolved phylogeny, despite a large number of molecular studies. Previous surveys have proposed trees with rather poorly supported nodes, and the bias caused by genetic introgression or ancestral polymorphism has not been assessed. Here we studied the phylogenetic relationships of all extant species of Equidae by analyzing 22 partial mitochondrial and nuclear genes using maximum likelihood and Bayesian inferences that account for heterogeneous gene histories. We also examined genetic signatures of lineage sorting and/or genetic introgression in zebras by evaluating patterns of intraspecific genetic variation. Our study improved the resolution and support of the Equus phylogeny and in particular the controversial positions of the African wild ass (E. asinus) and mountain zebra (E. zebra): the African wild ass is placed as a sister species of the Asiatic asses and the mountain zebra as the sister taxon of Grevy’s and Burchell’s zebras. A shared polymorphism (indel) detected among zebra species in the Estrogen receptor 1 gene was likely due to incomplete lineage sorting and not genetic introgression as also indicated by other mitochondrial (Cytochrome b) and nuclear (Y chromosome and microsatellites) markers. Ancestral polymorphism in equids might have contributed to the long-standing lack of clarity in the phylogeny of this highly threatened group of mammals.

Extraordinarily rapid life-history divergence between Cryptasterina sea star species

Life history plays a critical role in governing microevolutionary processes such as gene flow and adaptation, as well as macroevolutionary processes such speciation. Here, we use multilocus phylogeographic analyses to examine a speciation event involving spectacular life-history differences between sister species of sea stars. Cryptasterina hystera has evolved a suite of derived life-history traits (including internal self-fertilization and brood protection) that differ from its sister species Cryptasterina pentagona, a gonochoric broadcast spawner. We show that these species have only been reproductively isolated for approximately 6000 years (95% highest posterior density of 905-22 628), and that this life-history change may be responsible for dramatic genetic consequences, including low nucleotide diversity, zero heterozygosity and no gene flow. The rapid divergence of these species rules out some mechanisms of isolation such as adaptation to microhabitats in sympatry, or slow divergence by genetic drift during prolonged isolation. We hypothesize that the large phenotypic differences between species relative to the short divergence time suggests that the life-history differences observed may be direct responses to disruptive selection between populations. We speculate that local environmental or demographic differences at the southern range margin are possible mechanisms of selection driving one of the fastest known marine speciation events.

Introgression between cultivars and wild populations of Momordica charantia L. (Cucurbitaceae) in Taiwan.

The landrace strains of Momordica charantia are widely cultivated vegetables throughout the tropics and subtropics, but not in Taiwan, a continental island in Southeast Asia, until a few hundred years ago. In contrast, the related wild populations with smaller fruit sizes are native to Taiwan. Because of the introduction of cultivars for agricultural purposes, these two accessions currently exhibit a sympatric or parapatric distribution in Taiwan. In this study, the cultivars and wild samples from Taiwan, India, and Korea were collected for testing of their hybridization and evolutionary patterns. The cpDNA marker showed a clear distinction between accessions of cultivars and wild populations of Taiwan and a long divergence time. In contrast, an analysis of eight selectively neutral nuclear microsatellite loci did not reveal a difference between the genetic structures of these two accessions. A relatively short divergence time and frequent but asymmetric gene flows were estimated based on the isolation-with-migration model. Historical and current introgression from cultivars to wild populations of Taiwan was also inferred using MIGRATE-n and BayesAss analyses. Our results showed that these two accessions shared abundant common ancestral polymorphisms, and the timing of the divergence and colonization of the Taiwanese wild populations is consistent with the geohistory of the Taiwan Strait land bridge of the Last Glacial Maximum (LGM). Long-term and recurrent introgression between accessions indicated the asymmetric capacity to receive foreign genes from other accessions. The modern introduction of cultivars of M. charantia during the colonization of Taiwan by the Han Chinese ethnic group enhanced the rate of gene replacement in the native populations and resulted in the loss of native genes.

No evidence for Z-chromosome rearrangements between the pied flycatcher and the collared flycatcher as judged by gene-based comparative genetic maps

Revealing the genetic basis of reproductive isolation is fundamental for understanding the speciation process. Chromosome speciation models propose a role for chromosomal rearrangements in promoting the build up of reproductive isolation between diverging populations and empirical data from several animal and plant taxa support these models. The pied flycatcher and the collared flycatcher are two closely related species that probably evolved reproductive isolation during geographical separation in Pleistocene glaciation refugia. Despite the short divergence time and current hybridization, these two species demonstrate a high degree of intrinsic post-zygotic isolation and previous studies have shown that traits involved in mate choice and hybrid viability map to the Z-chromosome. Could rearrangements of the Z-chromosome between the species explain their reproductive isolation? We developed high coverage Z-chromosome linkage maps for both species, using gene-based markers and large-scale SNP genotyping. Best order maps contained 57-62 gene markers with an estimated average density of one every 1-1.5 Mb. We estimated the recombination rates in flycatcher Z-chromosomes to 1.1-1.3 cM/Mb. A comparison of the maps of the two species revealed extensive co-linearity with no strong evidence for chromosomal rearrangements. This study does therefore not provide support the idea that sex chromosome rearrangements have caused the relatively strong post-zygotic reproductive isolation between these two Ficedula species.

Mammalian Evolution May not Be Strictly Bifurcating

The massive amount of genomic sequence data that is now available for analyzing evolutionary relationships among 31 placental mammals reduces the stochastic error in phylogenetic analyses to virtually zero. One would expect that this would make it possible to finally resolve controversial branches in the placental mammalian tree. We analyzed a 2,863,797 nucleotide-long alignment (3,364 genes) from 31 placental mammals for reconstructing their evolution. Most placental mammalian relationships were resolved, and a consensus of their evolution is emerging. However, certain branches remain difficult or virtually impossible to resolve. These branches are characterized by short divergence times in the order of 1–4 million years. Computer simulations based on parameters from the real data show that as little as about 12,500 amino acid sites could be sufficient to confidently resolve short branches as old as about 90 million years ago (Ma). Thus, the amount of sequence data should no longer be a limiting factor in resolving the relationships among placental mammals. The timing of the early radiation of placental mammals coincides with a period of climate warming some 100–80 Ma and with continental fragmentation. These global processes may have triggered the rapid diversification of placental mammals. However, the rapid radiations of certain mammalian groups complicate phylogenetic analyses, possibly due to incomplete lineage sorting and introgression. These speciation-related processes led to a mosaic genome and conflicting phylogenetic signals. Split network methods are ideal for visualizing these problematic branches and can therefore depict data conflict and possibly the true evolutionary history better than strictly bifurcating trees. Given the timing of tectonics, of placental mammalian divergences, and the fossil record, a Laurasian rather than Gondwanan origin of placental mammals seems the most parsimonious explanation.