Molecular phylogenetics of the lizard genus Microlophus (squamata:tropiduridae): aligning and retrieving indel signal from nuclear introns

Diversification of the African legless skinks in the subfamily Acontinae (Family Scincidae)

Cladogenic diversification is often explained by referring to climatic oscillations and geomorphic shifts that cause allopatric speciation. In this regard, southern Africa retains a high level of landscape heterogeneity in vegetation, geology, and rainfall patterns. The legless skink subfamily Acontinae occurs broadly across the southern African subcontinent and therefore provides an ideal model group for investigating biogeographic patterns associated with the region. A robust phylogenetic study of the Acontinae with comprehensive coverage and adequate sampling of each taxon has been lacking up until now, resulting in unresolved questions regarding the subfamily's biogeography and evolution. In this study, we used multi-locus genetic markers (three mitochondrial and two nuclear) with comprehensive taxon coverage (all currently recognized Acontinae species) and adequate sampling (multiple specimens for most taxa) of each taxon to infer a phylogeny for the subfamily. The phylogeny retrieved four well-supported clades in Acontias and supported the monophyly of Typhlosaurus. Following the General Lineage Concept (GLC), many long-standing phylogenetic enigmas within Acontias occidentalis and the A. kgalagadi, A. lineatus and A. meleagris species complexes, and within Typhlosaurus were resolved. Our species delimitation analyses suggest the existence of hidden taxa in the A. occidentalis, A. cregoi and A. meleagris species groups, but also suggest that some currently recognized species in the A. lineatus and A. meleagris species groups, and within Typhlosaurus, should be synonymised. We also possibly encountered "ghost introgression" in A. occidentalis. Our inferred species tree revealed a signal of gene flow, which implies possible cross-over in some groups. Fossil evidence calibration dating results showed that the divergence between Typhlosaurus and Acontias was likely influenced by cooling and increasing aridity along the southwest coast in the mid-Oligocene caused by the opening of the Drake Passage. Further cladogenesis observed in Typhlosaurus and Acontias was likely influenced by Miocene cooling, expansion of open habitat, uplifting of the eastern Great Escarpment (GE), and variation in rainfall patterns, together with the effect of the warm Agulhas Current since the early Miocene, the development of the cold Benguela Current since the late Miocene, and their co-effects. The biogeographic pattern of the Acontinae bears close resemblance to that of other herpetofauna (e.g., rain frogs and African vipers) in southern Africa.

Chromosome-level genome assembly of the bar-headed goose (Anser indicus)

Bar-headed geese (Anser indicus) are adaptable to plateau environments. In this study, we sequenced and assembled a high-quality chromosome-level genome of the bar-headed goose using PacBio long reads and Hi-C technique, and generated 115.73 Gb of Illumina short-reads and 95.89 Gb of PacBio long-reads. The assembled bar-headed goose genome, with a contig N50 of 5.734 Mb and a scaffold N50 of 65.77 Mb, is 1.129 Gb in length and includes 33 chromosomes and 451 fragments. BUSCO assessment yielded a completeness score of 94.4%. In total, 15,376 protein-coding genes were predicted, of which 94.95% had homologs in protein databases. We identified 78 positively selected genes (PSGs) in the bar-headed goose genome, which were mainly enriched in calcium ion and ATP-binding. This bar-headed goose genome will be an important resource for increasing our understanding regarding the genetic basis of adaptation to life at a high altitude.

The reference genome and transcriptome of the limestone langur, Trachypithecus leucocephalus, reveal expansion of genes related to alkali tolerance

Background

Trachypithecus leucocephalus, the white-headed langur, is a critically endangered primate that is endemic to the karst mountains in the southern Guangxi province of China. Studying the genomic and transcriptomic mechanisms underlying its local adaptation could help explain its persistence within a highly specialized ecological niche.

Results

In this study, we used PacBio sequencing and optical assembly and Hi-C analysis to create a high-quality de novo assembly of the T. leucocephalus genome. Annotation and functional enrichment revealed many genes involved in metabolism, transport, and homeostasis, and almost all of the positively selected genes were related to mineral ion binding. The transcriptomes of 12 tissues from three T. leucocephalus individuals showed that the great majority of genes involved in mineral absorption and calcium signaling were expressed, and their gene families were significantly expanded. For example, FTH1 primarily functions in iron storage and had 20 expanded copies.

Conclusions

These results increase our understanding of the evolution of alkali tolerance and other traits necessary for the persistence of T. leucocephalus within an ecologically unique limestone karst environment.

Geometric morphometric analysis of the head of Microlophus atacamensis (Tropiduridae) in a latitudinal gradient

BACKGROUND

Clinal variation is defined as gradual variation in a character associated with geographic distance among sites. Microlophus atacamensis is a medium large lizard species which inhabits the intertidal zone of northern Chile, distributed from Antofagasta (23° 39' S) to Arrayán, La Serena (29° 41' S) in a gradient which could show clinal variation. Geometric morphometrics analyzes differences in shape independent of size; information about shape allows a more complete biological interpretation than information on morphological variation. This studyevaluated clinal variation in the head shape of M.atacamensis fromfive localities (Antofagasta, Paposo, Pan de Azúcar, Caldera, and Tres Playitas) using dorsal and lateral views, comparing form variation with latitudinal distribution. The heads of adults collected were photographed in lateral and dorsal views.

RESULTS

The analysis did not find significant differences in form among the five localities, in contrast to the proposal of earlier studies, and no differences were recorded between the sexes. Possible reasons why these populations are not differentiated in the latitudinal gradient are discussed.

CONCLUSIONS

Our results show that there are no differences between the studied, among the explanations are that populationsof this species are subjected to similar environments that promote convergence of the structures analyzed.

Molecular phylogeny and new classification of the genera Eulophias and Zoarchias (PISCES, Zoarcoidei)

Morphological and osteological studies of the Zoarcoidei group have previously been undertaken, but the group (especially the genera Eulophias and Zoarchias) still remains enigmatic. Therefore, we conducted molecular phylogenetic studies on the two genera Eulophias and Zoarchias using two mitochondrial (16S rRNA and COI) and two nuclear genes (RAG2 and RNF213). Our phylogenetic analysis supported the monophyly of the suborder level of the Zoarcoidei, but rejected the previous morphology- and osteology-based classification hypotheses regarding the two genera. Conflict between mtDNA and nDNA phylogenies within the genus Eulophias implies that the genus shows a complicated relationship such as hybridization in the process of the evolutionary history. The genetic distances between the Eulophias (or Zoarchias) and other Zoarcoidei spp. were the greatest, showing different family-level affiliations. In addition, the mtDNA topology showed the two genera were clearly separated from each other as well as from the families Stichaeidae and Zoarcidae. Considering the new molecular phylogeny, we suggest a new classification for the two genera: (1) Eulophias belongs to a new family named as the Eulophiidae; (2) Zoarchias belongs to the family Neozoarcidae (sensu Radchenko et al., 2012b) rather than to Stichaeidae and Zoarcidae.

Parsimony and model-based analyses of indels in avian nuclear genes reveal congruent and incongruent phylogenetic signals

Insertion/deletion (indel) mutations, which are represented by gaps in multiple sequence alignments, have been used to examine phylogenetic hypotheses for some time. However, most analyses combine gap data with the nucleotide sequences in which they are embedded, probably because most phylogenetic datasets include few gap characters. Here, we report analyses of 12,030 gap characters from an alignment of avian nuclear genes using maximum parsimony (MP) and a simple maximum likelihood (ML) framework. Both trees were similar, and they exhibited almost all of the strongly supported relationships in the nucleotide tree, although neither gap tree supported many relationships that have proven difficult to recover in previous studies. Moreover, independent lines of evidence typically corroborated the nucleotide topology instead of the gap topology when they disagreed, although the number of conflicting nodes with high bootstrap support was limited. Filtering to remove short indels did not substantially reduce homoplasy or reduce conflict. Combined analyses of nucleotides and gaps resulted in the nucleotide topology, but with increased support, suggesting that gap data may prove most useful when analyzed in combination with nucleotide substitutions.

Iteratively refined guide trees help improving alignment and phylogenetic inference in the mushroom family Bolbitiaceae

Reconciling traditional classifications, morphology, and the phylogenetic relationships of brown-spored agaric mushrooms has proven difficult in many groups, due to extensive convergence in morphological features. Here, we address the monophyly of the Bolbitiaceae, a family with over 700 described species and examine the higher-level relationships within the family using a newly constructed multilocus dataset (ITS, nrLSU rDNA and EF1-alpha). We tested whether the fast-evolving Internal Transcribed Spacer (ITS) sequences can be accurately aligned across the family, by comparing the outcome of two iterative alignment refining approaches (an automated and a manual) and various indel-treatment strategies. We used PRANK to align sequences in both cases. Our results suggest that--although PRANK successfully evades overmatching of gapped sites, referred previously to as alignment overmatching--it infers an unrealistically high number of indel events with natively generated guide-trees. This 'alignment undermatching' could be avoided by using more rigorous (e.g. ML) guide trees. The trees inferred in this study support the monophyly of the core Bolbitiaceae, with the exclusion of Panaeolus, Agrocybe, and some of the genera formerly placed in the family. Bolbitius and Conocybe were found monophyletic, however, Pholiotina and Galerella require redefinition. The phylogeny revealed that stipe coverage type is a poor predictor of phylogenetic relationships, indicating the need for a revision of the intrageneric relationships within Conocybe.

Re-mind the gap! Insertion - deletion data reveal neglected phylogenetic potential of the nuclear ribosomal internal transcribed spacer (ITS) of fungi

Rapidly evolving, indel-rich phylogenetic markers play a pivotal role in our understanding of the relationships at multiple levels of the tree of life. There is extensive evidence that indels provide conserved phylogenetic signal, however, the range of phylogenetic depths for which gaps retain tree signal has not been investigated in detail. Here we address this question using the fungal internal transcribed spacer (ITS), which is central in many phylogenetic studies, molecular ecology, detection and identification of pathogenic and non-pathogenic species. ITS is repeatedly criticized for indel-induced alignment problems and the lack of phylogenetic resolution above species level, although these have not been critically investigated. In this study, we examined whether the inclusion of gap characters in the analyses shifts the phylogenetic utility of ITS alignments towards earlier divergences. By re-analyzing 115 published fungal ITS alignments, we found that indels are slightly more conserved than nucleotide substitutions, and when included in phylogenetic analyses, improved the resolution and branch support of phylogenies across an array of taxonomic ranges and extended the resolving power of ITS towards earlier nodes of phylogenetic trees. Our results reconcile previous contradicting evidence for the effects of data exclusion: in the case of more sophisticated indel placement, the exclusion of indel-rich regions from the analyses results in a loss of tree resolution, whereas in the case of simpler alignment methods, the exclusion of gapped sites improves it. Although the empirical datasets do not provide to measure alignment accuracy objectively, our results for the ITS region are consistent with previous simulations studies alignment algorithms. We suggest that sophisticated alignment algorithms and the inclusion of indels make the ITS region and potentially other rapidly evolving indel-rich loci valuable sources of phylogenetic information, which can be exploited at multiple taxonomic levels.

Incorporating indels as phylogenetic characters: impact for interfamilial relationships within Arctoidea (Mammalia: Carnivora)

Insertion and deletion events (indels) provide a suite of markers with enormous potential for molecular phylogenetics. Using many more indel characters than those in previous studies, we here for the first time address the impact of indel inclusion on the phylogenetic inferences of Arctoidea (Mammalia: Carnivora). Based on 6843 indel characters from 22 nuclear intron loci of 16 species of Arctoidea, our analyses demonstrate that when the indels were not taken into consideration, the monophyly of Ursidae and Pinnipedia tree and the monophyly of Pinnipedia and Musteloidea tree were both recovered, whereas inclusion of indels by using three different indel coding schemes give identical phylogenetic tree topologies supporting the monophyly of Ursidae and Pinnipedia. Our work brings new perspectives on the previously controversial placements among Arctoidea families, and provides another example demonstrating the importance of identifying and incorporating indels in the phylogenetic analyses of introns. In addition, comparison of indel incorporation methods revealed that the three indel coding methods are all advantageous over treating indels as missing data, given that incorporating indels produces consistent results across methods. This is the first report of the impact of different indel coding schemes on phylogenetic reconstruction at the family level in Carnivora, which indicates that indels should be taken into account in the future phylogenetic analyses.

Phylogenetic relationships among the colobine monkeys revisited: new insights from analyses of complete mt genomes and 44 nuclear non-coding markers

Background

Phylogenetic relationships among Asian and African colobine genera have been disputed and are not yet well established. In the present study, we revisit the contentious relationships within the Asian and African Colobinae by analyzing 44 nuclear non-coding genes (>23 kb) and mitochondrial (mt) genome sequences from 14 colobine and 4 non-colobine primates.

Principal Findings

The combined nuclear gene and the mt genome as well as the combined nuclear and mt gene analyses yielded different phylogenetic relationships among colobine genera with the exception of a monophyletic ‘odd-nosed’ group consisting of Rhinopithecus, Pygathrix and Nasalis, and a monophyletic African group consisting of Colobus and Piliocolobus. The combined nuclear data analyses supported a sister-grouping between Semnopithecus and Trachypithecus, and between Presbytis and the odd-nosed monkey group, as well as a sister-taxon association of Pygathrix and Rhinopithecus within the odd-nosed monkey group. In contrast, mt genome data analyses revealed that Semnopithecus diverged earliest among the Asian colobines and that the odd-nosed monkey group is sister to a Presbytis and Trachypithecus clade, as well as a close association of Pygathrix with Nasalis. The relationships among these genera inferred from the analyses of combined nuclear and mt genes, however, varied with the tree-building methods used. Another remarkable finding of the present study is that all of our analyses rejected the recently proposed African colobine paraphyly and hybridization hypothesis and supported reciprocal monophyly of the African and Asian groups.

Significance

The phylogenetic utility of large-scale new non-coding genes was assessed using the Colobinae as a model, We found that these markers were useful for distinguishing nodes resulting from rapid radiation episodes such as the Asian colobine radiation. None of these markers here have previously been used for colobine phylogenetic reconstruction, increasing the spectrum of molecular markers available to mammalian systematics.

Unravelling the peculiarities of island life: vicariance, dispersal and the diversification of the extinct and extant giant Galápagos tortoises

In isolated oceanic islands, colonization patterns are often interpreted as resulting from dispersal rather than vicariant events. Such inferences may not be appropriate when island associations change over time and new islands do not form in a simple linear trend. Further complexity in the phylogeography of ocean islands arises when dealing with endangered taxa as extinctions, uncertainty on the number of evolutionary 'units', and human activities can obscure the progression of colonization events. Here, we address these issues through a reconstruction of the evolutionary history of giant Galápagos tortoises, integrating DNA data from extinct and extant species with information on recent human activities and newly available geological data. Our results show that only three of the five extinct or nearly extinct species should be considered independent evolutionary units. Dispersal from mainland South America started at approximately 3.2 Ma after the emergence of the two oldest islands of San Cristobal and Española. Dispersal from older to younger islands began approximately 1.74 Ma and was followed by multiple colonizations from different sources within the archipelago. Vicariant events, spurred by island formation, coalescence, and separation, contributed to lineage diversifications on Pinzón and Floreana dating from 1.26 and 0.85 Ma. This work provides an example of how to reconstruct the history of endangered taxa in spite of extinctions and human-mediated dispersal events and highlights the need to take into account both vicariance and dispersal when dealing with organisms from islands whose associations are not simply explained by a linear emergence model.

On the phylogeny of Mustelidae subfamilies: analysis of seventeen nuclear non-coding loci and mitochondrial complete genomes

BACKGROUND

Mustelidae, as the largest and most-diverse family of order Carnivora, comprises eight subfamilies. Phylogenetic relationships among these Mustelidae subfamilies remain argumentative subjects in recent years. One of the main reasons is that the mustelids represent a typical example of rapid evolutionary radiation and recent speciation event. Prior investigation has been concentrated on the application of different mitochondrial (mt) sequence and nuclear protein-coding data, herein we employ 17 nuclear non-coding loci (>15 kb), in conjunction with mt complete genome data (>16 kb), to clarify these enigmatic problems.

RESULTS

The combined nuclear intron and mt genome analyses both robustly support that Taxidiinae diverged first, followed by Melinae. Lutrinae and Mustelinae are grouped together in all analyses with strong supports. The position of Helictidinae, however, is enigmatic because the mt genome analysis places it to the clade uniting Lutrinae and Mustelinae, whereas the nuclear intron analysis favors a novel view supporting a closer relationship of Helictidinae to Martinae. This finding emphasizes a need to add more data and include more taxa to resolve this problem. In addition, the molecular dating provides insights into the time scale of the origin and diversification of the Mustelidae subfamilies. Finally, the phylogenetic performances and limits of nuclear introns and mt genes are discussed in the context of Mustelidae phylogeny.

CONCLUSION

Our study not only brings new perspectives on the previously obscured phylogenetic relationships among Mustelidae subfamilies, but also provides another example demonstrating the effectiveness of nuclear non-coding loci for reconstructing evolutionary histories in a group that has undergone rapid bursts of speciation.

Progressive colonization and restricted gene flow shape island-dependent population structure in Galápagos marine iguanas (Amblyrhynchus cristatus)

BACKGROUND

Marine iguanas (Amblyrhynchus cristatus) inhabit the coastlines of large and small islands throughout the Galápagos archipelago, providing a rich system to study the spatial and temporal factors influencing the phylogeographic distribution and population structure of a species. Here, we analyze the microevolution of marine iguanas using the complete mitochondrial control region (CR) as well as 13 microsatellite loci representing more than 1200 individuals from 13 islands.

RESULTS

CR data show that marine iguanas occupy three general clades: one that is widely distributed across the northern archipelago, and likely spread from east to west by way of the South Equatorial current, a second that is found mostly on the older eastern and central islands, and a third that is limited to the younger northern and western islands. Generally, the CR haplotype distribution pattern supports the colonization of the archipelago from the older, eastern islands to the younger, western islands. However, there are also signatures of recurrent, historical gene flow between islands after population establishment. Bayesian cluster analysis of microsatellite genotypes indicates the existence of twenty distinct genetic clusters generally following a one-cluster-per-island pattern. However, two well-differentiated clusters were found on the easternmost island of San Cristóbal, while nine distinct and highly intermixed clusters were found on youngest, westernmost islands of Isabela and Fernandina. High mtDNA and microsatellite genetic diversity were observed for populations on Isabela and Fernandina that may be the result of a recent population expansion and founder events from multiple sources.

CONCLUSIONS

While a past genetic study based on pure FST analysis suggested that marine iguana populations display high levels of nuclear (but not mitochondrial) gene flow due to male-biased dispersal, the results of our sex-biased dispersal tests and the finding of strong genetic differentiation between islands do not support this view. Therefore, our study is a nice example of how recently developed analytical tools such as Bayesian clustering analysis and DNA sequence-based demographic analyses can overcome potential biases introduced by simply relying on FST estimates from markers with different inheritance patterns.

Island biogeography of Galápagos lava lizards (Tropiduridae: Microlophus): species diversity and colonization of the archipelago

The "lava lizards" (Microlophus) are distributed throughout the Galápagos Archipelago, and consist of radiations derived from two independent colonizations. The "Eastern Radiation" includes M. bivittatus and M. habeli endemic to San Cristobal and Marchena Islands. The "Western Radiation" includes five to seven historically recognized species distributed across almost the entire Archipelago. We combine dense geographic sampling and multilocus sequence data to estimate a phylogenetic hypothesis for the Western Radiation, to delimit species boundaries in this radiation, and to estimate a time frame for colonization events. Our phylogenetic hypothesis rejects two earlier topologies for the Western Radiation and paraphyly of M. albemarlensis, while providing strong support for single colonizations on each island. The colonization history implied by our phylogeny is consistent with general expectations of an east-to-west route predicted by the putative age of island groups, and prevailing ocean currents in the Archipelago. Additionally, combined evidence suggests that M. indefatigabilis from Santa Fe should be recognized as a full species. Finally, molecular divergence estimates suggest that the two colonization events likely occurred on the oldest existing islands, and the Western Radiation represents a recent radiation that, in most cases, has produced species that are considerably younger than the islands they inhabit.

Colonization and diversification of Galápagos terrestrial fauna: a phylogenetic and biogeographical synthesis

Remote oceanic islands have long been recognized as natural models for the study of evolutionary processes involved in diversification. Their remoteness provides opportunities for isolation and divergence of populations, which make islands remarkable settings for the study of diversification. Groups of islands may share a relatively similar geological history and comparable climate, but their inhabitants experience subtly different environments and have distinct evolutionary histories, offering the potential for comparative studies. A range of organisms have colonized the Galápagos Islands, and various lineages have radiated throughout the archipelago to form unique assemblages. This review pays particular attention to molecular phylogenetic studies of Galápagos terrestrial fauna. We find that most of the Galápagos terrestrial fauna have diversified in parallel to the geological formation of the islands. Lineages have occasionally diversified within islands, and the clearest cases occur in taxa with very low vagility and on large islands with diverse habitats. Ecology and habitat specialization appear to be critical in speciation both within and between islands. Although the number of phylogenetic studies is continuously increasing, studies of natural history, ecology, evolution and behaviour are essential to completely reveal how diversification proceeded on these islands.