Next Article in Journal
A Structure-Based Drug Discovery Paradigm
Previous Article in Journal
Reversal of MK-801-Induced Disruptions in Social Interactions and Working Memory with Simultaneous Administration of LY487379 and VU152100 in Mice
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 11
10.3390/ijms20112782
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Shortcomings of Phylogenetic Studies on Recent Radiated Insular Groups: A Meta-Analysis Using Cabo Verde Biodiversity

by
Maria M. Romeiras
1,2,*,
Ana Rita Pena
2,
Tiago Menezes
3,
Raquel Vasconcelos
4,
Filipa Monteiro
1,2,
Octávio S. Paulo
2 and
Mónica Moura
3
1
LEAF, Linking Landscape, Environment, Agriculture and Food, Instituto Superior de Agronomia, Universidade de Lisboa, 1349-017 Lisbon, Portugal
2
Centre for Ecology, Evolution and Environmental Changes (cE3c), Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisbon, Portugal
3
CIBIO, Research Centre in Biodiversity and Genetic Resources, Azores Group, InBIO Associate Laboratory, Universidade dos Açores, 9501-855 Ponta Delgada, Azores, Portugal
4
CIBIO, Research Centre in Biodiversity and Genetic Resources, InBIO Associate Laboratory, Universidade do Porto, 4485-661 Vairão, Portugal
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(11), 2782; https://doi.org/10.3390/ijms20112782
Submission received: 1 March 2019 / Revised: 16 May 2019 / Accepted: 4 June 2019 / Published: 6 June 2019
(This article belongs to the Section Molecular Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Over the previous decades, numerous studies focused on how oceanic islands have contributed to determine the phylogenetic relationships and times of origin and diversification of different endemic lineages. The Macaronesian Islands (i.e., Azores, Madeira, Selvagens, Canaries, and Cabo Verde), harbour biotas with exceptionally high levels of endemism. Within the region, the vascular plants and reptiles constitute two of the most important radiations. In this study we compare relevant published phylogenetic data and diversification rates retrieved within Cabo Verde endemic lineages and discuss the importance of choosing appropriate phylogeny-based methods to investigate diversification dynamics on islands. From this selective literature-based review, we summarize the software packages used in Macaronesian studies and discuss their adequacy considering the published data to obtain well-supported phylogenies in the target groups. We further debate the importance of Next Generation Sequencing (NGS), to investigate the evolutionary processes of diversification in the Macaronesian Islands. Analysis of genomic data provides phylogenetic resolution for rapidly evolving species radiations, suggesting a great potential to improve the phylogenetic signal and divergence time estimates in insular lineages. The most important Macaronesian reptile radiations provide good case-studies to compare classical phylogenetic methods with new tools, such as phylogenomics, revealing a high value for research on this hotspot area.

1. Background

During the last two decades, one of the fields that took the most benefit of the explosive growth of bioinformatic tools was phylogenetics, which uses molecular data to establish evolutionary relationships among species [1,2]. Presently, large amounts of genetic data are easily retrievable from public repositories, such as GenBank [3], which stores and organizes molecular sequences, while also making them publicly accessible and ready to use. Other biology focused databases include EMBL (European Molecular Biology Laboratory) [4], UniProt (Universal Protein resource) [5], DDBJ (DNA Data Bank of Japan) [6] and many others, all of them available online, allowing an ever-growing global network of information that continues to expand, giving researchers an almost never-ending flow of new data. However, an issue that arose from the high output of data was its processing. The past two decades have been astonishing, with more than a million-fold improvement in the rate of sequence generation [7], particularly with the introduction of Next Generation Sequencing (NGS), which is further discussed below. With such an extraordinary quantity of sequences being created, the number of bioinformatics tools required to handle such high levels of information, also grew and their analytic capacity dramatically improved. Currently, almost all systematic studies have some type of bioinformatics analysis and it is common to make use of multiple options.
An outstanding area for studies of evolution and speciation is the Macaronesian Region (i.e., Azores, Madeira, Selvages, Canaries, and Cabo Verde archipelagos; Figure 1A) that encompasses a great number of endemic lineages. These islands belong to the Mediterranean Basin biodiversity hotspot [8], which is characterized by high levels of endemism, and represents an excellent area to study evolution. The aim of this study is to review the present knowledge on the tree of life and evolution of Cabo Verde biodiversity using, as case-studies, the reptiles [9] and the vascular plants [10], which show high levels of endemism, with several range-restricted species occurring in different habitats, and thus constitute a good model system to understand patterns of diversification within this archipelago. We will debate the relevance of cutting-edge bioinformatics tools and methodologies to provide a better understanding of the evolutionary history of the recent radiated groups of the Macaronesian Region.

2. Phylogenetic Inference and Divergence Time Estimation

The number of bioinformatics tools available today represents a huge framework for each question that arises, presenting several alternative solutions. For phylogenetic analysis, this can prove to be an extraordinary time, as new knowledge, powerful software and hardware frameworks, are continuously improved [11,12]. Three major types of phylogenetic inference methods have been used to retrieve phylogenetic relationships within and between several Macaronesian lineages, namely, Maximum Parsimony (MP), Maximum Likelihood (ML) and Bayesian Inference (BI). Maximum Parsimony (MP) [13,14] has proven to be useful when dealing with recently diverged species with short branches and was commonly used as a standalone solution in earlier Macaronesian studies to resolve the phylogenetic relationships of some plant lineages (e.g., Convolvulus (Convolvulaceae) [15]; Gonosperminae (Asteraceae) [16]; Bencomia (Rosaceae) [17]; Pericallis (Asteraceae) [18]). More recent studies have followed the general trend, and usually couple MP approaches with other methods (i.e., ML and/or BI), providing better and more supported hypotheses, namely in several plant lineages: Lactuca (Asteraceae) [19]; Cistus (Cistaceae) [20]; Azorina, Campanula, Trachelium, Lobelia, Musschia, and Wahlenbergia (Campanulaceae) [21]; and endemic reptile genera: Gallotia (Lacertidae) [22]; Hemidactylus (Gekkonidae) [23].
When working with ML methods that were made popular by Fisher in the early 20th century [24], the objective is to infer an evolutionary tree, taking into consideration an evolutionary model which provides rates of nucleotide substitution. Likewise with MP methods, assessing each tree topology is impossible, thus heuristic search methods are used to obtain the most likely tree. More recently, the implementation of BI methods [25,26] became common, and these methods differ from ML by allowing a priori information to be assigned to the alignment, by estimating true probabilities and by resulting in a set of best, most probable, trees. A priori information defines an initial tree with a specified substitution model (e.g., Markov Chain Monte Carlo (MCMC)), distribution of rates across sites, branch lengths and tree topology [27,28].
Programs such as PAUP* (Phylogenetic Analysis Using Parsimony *and other methods) [29] and MEGA (Molecular Evolutionary Genetics Analysis) [30] have been widely used to conduct analyses with most of the above-mentioned phylogenetic inference methods in studies with reptile and with Macaronesian flora lineages [31,32,33,34]. Phylogenetic studies with Macaronesian reptiles [35,36,37,38,39,40,41,42,43,44,45,46,47,48,49] are summarised in Table 1, with detailed information on molecular markers and software used, analyses conducted and corresponding bibliographic reference. Although PAUP and MEGA allow estimating evolutionary models, an essential step in ML and BI analysis is the use of specific software such as jModelTest2 [50] and PartitionFinder2 [51]. Several of the selected research articles make use of these software, instead of only working with more generally purposed programs. jModelTest (and its predecessor ModelTest) seem to be more widely used, but more recent works seem to prefer using PartitionFinder instead (see more details for each lineage in Table 1).
To specifically perform ML and BI analysis with Macaronesian reptiles (Table 1) and plant lineages [19,20,21,31,32,33,34], two main software packages have also been often used: RAxML [52] and MrBayes [53]. Randomized Axelerated Maximum Likelihood (RAxML), and similar software packages such as PHYML [54] and GARLI [55], make use of ML estimations to perform phylogenetic inference; while MrBayes estimates phylogenies using a variant of the MCMC method called Metropolis-coupled MCMC. Newly generated and well-sampled phylogenies that made use of the ML algorithm implemented in RAxML and the Bayesian inference method in MrBayes, have contributed to establish the evolutionary and biogeography relationships within complex Macaronesian plant lineages, namely on Limonium, Plumbaginaceae [56] and on Campanulaceae endemic species [21].
Bayesian approaches also allow obtaining divergence time estimations [57,58], an analysis which provides an important insight into the evolutionary history of extant lineages and has been used to some extent with Macaronesian lineages (more details on published studies are provided below). An extremely important step in divergence time estimation is calibrating the clock (see column–Dating: Table 1). To calibrate internal nodes there are two main forms available, specifically, one that makes use of fossil and biogeographical (which requires some a priori knowledge of the relation between the species), and another that uses substitution rates previously calculated (secondary calibration). However, only reliable sources of secondary calibrations should be used, or a considerable error may be introduced in the analysis, as we further discuss in Section 4. As such, comparing several scenarios obtained with different calibration techniques is generally the best methodology.
Due to the complexity of the divergence time estimation methods, the choice of the software to perform these analyses is paramount. The most popular software used in this type of analysis is BEAST [59] and the more recent BEAST 2 [60], which make use of Bayesian multispecies coalescent methods to generate chronograms. However, working with high volumes of data, namely multi-individual and multilocus data, is a challenge, as it exponentially slows the analysis. To solve this setback a specific package was developed, StarBEAST or *BEAST [61], has shown to be very informative when working with rapid radiations [62]. This is of particular importance when working with insular species, as these systems are often fairly recent and under a considerate evolutionary pressure.
Other used approaches in Macaronesian phylogenetic inference studies, have been the estimation of phylogenetic networks, to test the existence of cryptic events (e.g., hybridization, recombination, horizontal gene transfer, or duplication or loss of genes) and often relate to intraspecific genealogies in the given molecular data set. The phylogenetic networks were used to produce complete community-level phylogenies and can be implemented in several software packages like TCS [63] and SplitsTree4 [64]. Particularly, TCS was used to analyse the genetic diversity for endemic plant lineages of Cabo Verde (i.e., Globularia amygdalifolia (Plantaginaceae), Cynanchum daltonii (Apocynaceae) and Umbilicus schmidtii (Crassulaceae) [65]) and of Azores (i.e., Ammi, Angelica lignescens, Azorina vidalii, Euphorbia stygiana, Pericallis malvifolia) [66] and of both archipelagos, i.e., Lactuca [19]. These methods also provided further insights into the population structure and genetic diversity of widespread plant lineages within Cabo Verde, Azores, Madeira and Canary Islands, such as the Campanulaceae [21], which highlights the usefulness of this tool for this type of assessment. On the other hand, SplitsTree4 contains methods specialized in the creation of recombination and hybridization networks, and were used in some studies with plants [67], where the obtained network revealed conflicting phylogenetic signals, such as rapid radiation or ancient hybridization events, within five Betoideae genera distributed in the Macaronesian Islands and in the western Mediterranean region.

3. Macaronesian Islands as Model Systems in Evolution

Islands can be broadly classified as one of two types: continental or oceanic. Continental islands are fragments of a larger landmass, whereas oceanic islands result from volcanoes rising above the water on or near a mid-ocean ridge. Particularly, oceanic islands, as is the case of the Macaronesian Islands, were the focus of studies that contributed to the understanding of speciation and the origins of biodiversity [68]. The terrestrial endemic lineages of the Macaronesian Islands are often characterized by occupying different habitats; showing striking morphological differences among species; frequent rarity and being present in few and small populations [65]. Over the past decades, research has focused on multiple facets of islands and how they have had a relevant role in the establishment of ecological and evolutionary theories [68,69].
Different patterns of colonization within the Macaronesian Region have been an ongoing subject of debate. Some of the islands within this region are near the mainland, with the closest being Fuerteventura in the Canary Islands, which distances only ca. 95 km west of mainland Africa (Figure 1A). This proximity may have facilitated dispersion and colonization with the influence of oceanic currents (North Atlantic and Canary currents), while north-easterly trade winds may have fostered the dissemination of ancestral species from the continent and nearby Canary Islands. For this reason, it has been proposed that frequent colonization events from the continent and close archipelagos may have contributed to the high number of endemic species in some of the less isolated Macaronesian Islands [70].

4. Phylogenetic Resolution and Divergence Time Estimation among Macaronesian Insular Groups–Case-Study Using A Meta-Analysis of Reptile Data

The wide variety of habitats within the Macaronesian Islands is home to a rich biodiversity, and over 1000 endemic plants have been identified so far. The evolutionary history of this huge diversity is not properly understood, and many lineages remain understudied, particularly in Azores and Cabo Verde. In fact, for Cabo Verde, few lineages were subject of molecular dating analyses to date, specifically: Campanula [21,71]; Diplotaxis [31]; Echium [32]; and Lotus [72], which represents about one-fourth of the total vascular plant diversity when considering lineages with more than one endemic species in the archipelago. On the other hand, only six reptile genera (i.e., Chalcides; Chioninia; Gallotia; Hemidactylus; Tarentola and Teira) occur within the Macaronesian Region (see Figure 1B), and from these, only three occur in Cabo Verde [9]. Due to the representativeness of reptile genera occurring in Cabo Verde, these taxa were selected as indicators of the present state of lineage resolution and divergence time estimation within Macaronesia. Twenty-three publications dealing with this group were selected, and a meta-analysis was conducted based on the methodology followed and results obtained. The present compilation included data from: (i) gene regions used; (ii) inference methods used to reconstruct phylogenies; (iii) completeness of the phylogenies generated. (iv) methods used to estimate the divergence times within each taxonomic group; (v) divergence times estimated.
The meta-analysis conducted revealed that although a huge number of methods and standard programs were used (Table 1; Figure 2, Figure 3 and Figure 4), there are still important knowledge gaps for some lineages. This analysis supports the idea that most of the endemic lineages are monophyletic but the huge morphologically diversity within some of these lineages was not always fully represented by the analysed taxa. On average, only 67%, 47% and 61% of the existent taxa were analysed in the selected studies for Chioninia, Hemidactylus and Tarentola, respectively (Figure 2). Furthermore, even within the most well-sampled taxonomic groups, outgroup selection was not consensual which might have introduced a bias in the obtained phylogenies.
While most of the divergence time estimations published to date were obtained with the application of robust bioinformatics methods, such as BEAST or *BEAST, it is evident that the lack of fossils for most lineages was an important impediment for the accurate calibration of molecular clocks (Figure 3 and Figure 4; only three studies used fossils). Due this shortcoming, several of the studies made use of mutation rates and/or geological information to date their phylogenies, which can lead to considerable error (Figure 4).
Geological dating in particular, can induce a significant amount of error when calibrating, as the accepted age of a given island or archipelago is often submitted to changes as new research methods and hypothesis are built. Some Cabo Verde Islands for example, have recently been found to be younger than previously thought, yet older island ages have been applied to dating methods thus promoting stochastic errors in divergence time estimations [73,74,75]. This means that studies based on the previously accepted geological ages were using an unfitting baseline for their calibrations, which could have led to incorrect estimations and conclusions. Another common mistake that can be made with this methodology is the incorrect assumption that all islands, and in particular those of volcanic origin and from the same archipelago, share a similar history throughout the eras and are able to harbour life all along their cycles, from start to present [76]. Each island has its own unique characteristics, and undergoes both long and short periods of change, as well as stationary phases, which dramatically changes their conditions and evolutionary pressures [77]. Cabo Verde archipelago is a good example of such a system, with different islands leading to diverse colonization events and endemic lineages. The glacial and interglacial periods are of special importance to the archipelago, since during the glacial periods, these regions were buffered from most of its effects (with temperatures variating very little) due to their tropical geographical and worked as a refuge for several species that eventually became extinct in neighbouring regions [78]. The rise and fall of the sea levels also plays an important part when studying evolutionary processes, and in particular with plants and reptiles, due to their influence on the available terrestrial area of the islands. During the interglacial periods the terrestrial area can increased dramatically (for example, in Cabo Verde, it was essentially doubled during the Last Glacial Maximum), which can allow for land bridges to form between neighbouring islands (e.g., Santa Luzia, Branco and Raso) and allow gene flow processes to occur, leading to a different evolutionary history of these thought-to-be isolated species. On the other hand, during glacial periods, this connection is severed, which not only leads to isolation and possible speciation events but can also induce several extinctions, as some organisms may not adapt to the new conditions. As such, when taking into account geological dates, particularly for oceanic volcanic islands, one must consider how these systems are often found to have very complex and non-linear geological histories. For non-insular systems, the use of this type of calibration is not as frequent, as dated continental geological events are usually punctual and restricted to particular regions [79]. The Qinghai-Tibetan Plateau (QTP) is a good example of a continental feature that has undergone several biological studies that used its formation as a calibration point to date their lineages [80,81]. The QTP and similar geological features are usually well-explored and dated but, like in insular systems, some care is necessary to assure that the calibrations are fitted to be used in a given phylogenetic reconstruction, and if possible, other complementary calibrations should be used.
All these factors hamper the development of well-supported hypotheses for the evolutionary history of the reptiles, a large portion of the endemic diversity found in the Macaronesia Region (see Figure 1B). Nonetheless, the vast amount of knowledge on Macaronesian reptile lineages has been continuously increasing along the years, and with the introduction of more advanced and progressive analysis methods, several knowledge landmarks have been achieved, such as the clarification of the Chioninia genus from the Mabuya superfamily [82] or the several phylogenetic studies on the Cabo Verde endemic lineages, that brought a clarification on their evolutionary story [39,49,83]. Recently, other methods opened new perspectives, namely the identification of plants, arthropods and bird species predated by different reptile taxa, using NGS methods (metabarcoding of faecal pellets), and a comparison with reptile lineage diversification occurring in different Cabo Verde Islands [84].
Over the past years, it was possible to witness a progressive development of the methodologies used in systematic studies of Macaronesian groups. Earlier studies used only a molecular method, such as MP, and in the case of plants, only one or two markers were used [15,16,17,18]. During the past decade, the combination of different methods, particularly the ML with BI analyses, and the addition of more DNA markers, allowed the development of more accurate phylogenies. For the Macaronesian reptiles, this new approach has proven to be very useful to reconstruct their evolutionary histories, as well as to clarify their taxonomic position. Possibly the most complex group of the three Cabo Verdean genera, the Chioninia, has undergone several taxonomic rearrangements and multiple analyses with sometimes disparate results. Problems with this group were mainly found in their relationship with its sister groups, which were previously assigned to a single large genus, Mabuya. Some studies place the Chioninia reptiles as close relatives to Middle-Eastern Trachylepis and some Mabuya species [85,86], while others place them alongside other African and Madagascan lineages [87]. A recent work by Karin et al. [37] was finally able to reconstruct a well-supported phylogeny for this reptile group that placed the Chioninia as a sister clade to Afro-Malagasy Trachylepis. To achieve this result, a wide sample of reptilian was used, and perhaps most importantly, several markers of both nuclear and mitochondrial loci were obtained, unlike previous works, which only used mitochondrial DNA.
Comparatively, the genera Hemidactylus and Tarentola are far less intricate, while still presenting their own unique challenges. The Hemidactylus genus is known for being one of the largest and more widespread reptile groups, with more than 100 described species, which makes it a common species to study due to its representativeness and multiple colonization events. Therefore, this is the well-studied group of the three Cabo Verde genera, both in the number of studies and its evolutionary and diversification history. Another positive feature of this reptile group is the availability of several fossil records, such as Sphaerodactylus, which were used to calibrate several Hemidactylus phylogenies [41,42,44,45]. As we previously discussed, fossil calibrations are generally very informative and precise when used properly, i.e., when they are properly justified [88,89]. The combination of these two factors (number of studies and fossil calibrations), has led to a mostly well-known phylogeny, which once again highlights the importance of the development of new methods of analysis and their use in numerous complementary works, which then lead to better and more interesting solutions.
Out of the three, the Tarentola is the most represented genus in the Macaronesia region, with some species in the Canary Islands, Cabo Verde, Madeira and Selvagens, while its origin region is thought to be the American continent. This group can be used as an example on how regionally focused studies can work in par with comprehensive ones that include all the species of a genus. Research in the Cabo Verde Tarentola species has been quite thorough, with several studies continuously improving the knowledge of their taxonomy and diversification patters [47,49]. In parallel, other studies shed some light in the evolutionary history of the genus, by hinting at a single colonization event from America to Africa [44], and by clarifying the relations between the several Mediterranean species [48]. When we combine these works, we can find a more robust solution, one that explains how and when this genus originated, how it came to Cabo Verde and how it has evolved since then, which is a collective knowledge that would not at all be possible with just one of these finds.
Even with all challenges considered, reptiles are still one of the most well-studied groups within the Macaronesian region, and their evolutionary patterns appear to be clearer and more comprehensible (Figure 2). Plants on the other hand, have a huge diversity of lineages that only recently have been studied using large-scale sampling within each archipelago and wide sets of molecular markers to trace their phylogenies, as is the case of the Plumbaginaceae [56] and Campanulaceae endemic species [21]. Furthermore, polytomies still persist frequently in plant phylogenies. For instance, the molecular phylogeny recently published by Franzke et al. [31] for the biggest plant radiation of Cabo Verde Islands (i.e., Diplotaxis, Brassicaceae) revealed that the tree is largely unresolved, and the clades within the tree show very weak support (55–69%). Therefore, particularly for plants, there is a need to perform more studies to understand the evolutionary patterns and dating analyses [90,91].

5. The Potential of Phylogenomics to Understand Evolutionary Relationships in Insular Lineages

As stated above, phylogenetic methods used to date to establish the relationships among insular lineages were mostly based on single locus or a small number of loci (i.e., nuclear and plastid genes for plants and mitochondrial genes for animals) and often the resulting trees included unresolved clades and were thus unable to fully reflect the speciation events within the Macaronesian Region.
Considering (i) the inherent drawbacks of the targeted molecular marker-based techniques, namely, the high effort required for finding a considerable number of informative loci; (ii) their little resolving power at low taxonomic levels, which ignores many population-level effects, and (iii) increasingly affordable sequencing costs resulting from the continuous advancement of next-generation sequencing (NGS) technologies, phylogenomics emerges as an promising alternative tool for studying evolutionary relationships based on comparative analysis of genome-scale data [92].
Phylogenomics can be seen as the phylogenetics branch which specializes in inferring species relationships based on whole genomes (WGS) or genome-partition approaches and is starting to be applied in evolutionary and ecological studies on Macaronesian islands [93,94]. The access to genomic data has proven to potentially alleviate previous problems of phylogenetics that resulted from stochastic error (limitation of sampling few genes) by expanding the number of characters [95]. However, genome wide analysis is not always straightforward, particularly when working with non-model organisms or large and complex genomes. On whole genomes approach, plastid sequencing is one of the most technically accessible regions of the genome, and its sequence conservation makes it a valuable region for comparative genome evolution, phylogenetic analysis and population studies [91]. Recent studies on endemic plants using insular lineages [96], and particularly from Azores archipelago (i.e., Laurus azorica (Seub.) Franco, [97]) show some promise to help resolving species relationships at phylogenetic level. From the various genome-partitioning strategies to sequence-selected subsets of the genome (reviewed in [96]), restriction site associated DNA sequencing (RAD-Seq) has emerged as powerful alternative to WGS, where is sequenced a representative portion of the genome, which allows for a greater depth of coverage and less financial costs, compared to complete genome sequencing [98,99,100]. As such, several studies found different uses for this methodology, that range from alignment and gene identification processes, to advanced genomic structuring evaluation [101]. Even though this method is broadly used, there are still issues to be taken in consideration, namely the high bioinformatics hub to detect loci involved in adaptation and to infer possible linkage disequilibrium effects in the data [102].
Another consideration when performing RAD-Seq analyses is the type of assembly to be used–de novo or referenced-based. As complete genomes are still scarcely available for most of the Macaronesian lineages, the selection of the method is usually limited, and the most viable solution is a de novo assembly [103]. The use of de novo assembling may provide interesting results, even when faced with a substantial volume of missing data, due to the higher number of loci on the analyses [104].
A scenario of rapid plant and reptile radiations is found among most of the Macaronesian lineages, with molecular differences not being correctly depicted by traditional molecular markers and phylogenetic analyses, mostly due to low accumulation of differences in the genomic regions screened between different species in different islands or even within species at intra-islands levels/different habitats [104,105]. Macaronesian Islands can also be used as model to study evolutionary processes of genotype-to-phenotype adaptations, and explore the adaptation of different taxa to certain habitat conditions, within each archipelago. Thus, when searching for adaptive traits, proper genomic tools are needed to detect population genetic divergence associated with processes of local adaptation [106]. NGS technologies can produce informative molecular markers that could address population differentiation and phylogeography issues, as well as find signatures of selection at a genomic scale [107]. High-throughput discovery of single nucleotide polymorphisms (SNPs) using RAD-seq [104] can also be effectively conducted on species without a reference genome, enabling studies in non-model organisms as was recently published for reptiles from Canary Islands [108]. By assessing diversity in traits under strong selective pressure, population genomics can be used to identify signatures of selection and thus genes of adaptive/fitness significance can be potentially recognized [109]. Also, SNPs that are linked to traits under selection, are highly valuable for identifying genetic loci/regions that contribute to phenotypic variation. As such, the genetic basis of adaptation under an island scenario could be a significant advance in insular evolutionary discipline, for example, could be explored to tackle genetic divergence between different species to assess colonization events, genetic variability among endemic taxa and identification of signatures of selection at a genomic scale linked to adaptive evolution on taxa present in the same environment (i.e., dry conditions) [109,110]. In an insular context, the occurrence of morphologically similar but ecologically diverse species on each island, or even within different habitats of the same island is fairly common, thus making insular ecosystems the perfect models to investigate niche shifts and adaptation [108].
In regard to insular evolutionary studies, RAD-Seq has proven to be very effective to contribute to clarify the relationships within closely related taxa, as well as providing a good alternative to study groups where Sanger sequencing failed to provide accurate results [93,111]. In cases where closely resembling species were involved, it is also generally proficient in providing good phylogenies, even when faced with poor genetic signals and conflict between gene trees [112]. Studies using RAD-Seq technology plants from Macaronesian Islands were able to resolve relationships in lineages that have undergone recent, rapid diversification resulting in extensive ecological and morphological diversity, as in Tolpis [95] and Micromeria [113].
Finally, other phylogenomics methods can be used to address the current shortcomings of relatively deep phylogenetic divergences (several million years) to shallow-scale divergences are based target enrichment to sequence protein-coding exons [114,115] or anchored hybrid enrichment [116]. High-throughput sequencing (HTS) techniques provide the opportunity to use thousands of loci to address phylogenetic problems. As well as these new technical capabilities created by NGS, there are also new types of problems that have to be addressed in phylogenomics, such as incomplete lineage sorting, gene duplication and loss, and horizontal transfer.

6. Final Considerations

The interdisciplinary field of bioinformatics was born from the need to develop and employ methodologies to study large amounts of biological data, namely proteins and DNA. Many scientific areas benefited greatly from the introduction of specialized and powerful tools that improved the implementation of databases, data analysis and biological interpretations. For evolutionary study fields such as phylogenetics, the knowledge jump has been substantial since the introduction of bioinformatic methods in their analysis. Making use of highly advanced software packages it is now possible to perform high-throughput DNA sequencing, and thoroughly analyse this data with several methods of phylogenetic inference, to obtain informative and well supported phylogenies. Moreover, the introduction of NGS was one of the biggest achievements in the biology field in the last couple of decades. Several biological fields were revolutionized and some only reached their true potential when high-throughput sequencing methods were finally widely available at a relatively low cost. For phylogenetic sciences, this was an revolutionary occurrence, as it was finally possible to sequence, store and analyse complete genomes, and in doing so, reconstruct more comprehensive phylogenies.
The power of phylogenetic sciences has been crucial as tools in the development and consolidation of knowledge of insular native/endemic species in Macaronesia, clarifying their classifications, understanding genetic patterns and evolutionary processes that occurred on the islands and/or archipelagos. The combined data from phylogeny and divergence time estimations are able to clarify diversification hypothesis and assign probable colonization events, and can even help discover new species, providing useful information for the field of conservation biology in insular environments. Phylogenomics are at a turning point and it is highly likely that in the next few years, new understanding on several issues regarding species diversification and adaptation will be instigated by the application of new methodologies to natural laboratories such as the Macaronesian Islands.

Author Contributions

Conceptualization, M.M.R.; methodology, M.M.R and M.M.; formal analysis, A.R.P. and M.M.; writing—original draft preparation, M.M.R., A.R.P., T.M., F.M. and M.M.; writing—review and editing, M.M.R., A.R.P., T.M., F.M., R.V., O.S.P. and M.M. All authors have approved the submitted version of this manuscript.

Funding

This research was funded by Fundação para a Ciência e Tecnologia (FCT) and Aga Khan Development Network (AKDN) under the project CVAgrobiodiversity/333111699 and in the scope of MACDIV–FCT-PTDC/BIABIC/0054/2014 (FCT). Fellowships SFRH/BPD/79913/2011 to RV and SFRH/BPD/104629/2014 to FM were funded by Portuguese National Funds through FCT, Portugal. Also to research units: UID/AGR/04129/2019 (LEAF) and UID/BIA/00329/2019 (cE3c).

Acknowledgments

The authors would like to acknowledge the supported provided by Fundação para a Ciência e Tecnologia (FCT) and Aga Khan Development Network (AKDN).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. López-Fernández, H.; Duque, P.; Henriques, S.; Vázquez, N.; dez-Riverola, F.; Vieira, C.P.; Reboiro-Jato, M.; Viera, J. Bioinformatics Protocols for Quickly Obtaining Large-Scale Data Sets for Phylogenetic Inferences. Interdiscip. Sci. Comput. Life Sci. 2018, 1–9. [Google Scholar] [CrossRef] [PubMed]
  2. Roumpeka, D.D.; Wallace, R.J.; Escalettes, F.; Fotheringham, I.; Watson, M. A review of bioinformatics tools for bio-prospecting from metagenomic sequence data. Front. Genet. 2017, 8, 23. [Google Scholar] [CrossRef] [PubMed]
  3. Clark, K.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Sayers, E.W. GenBank. Nucleic Acids Res. 2016, 44, D67–D72. [Google Scholar] [CrossRef] [PubMed]
  4. Kulikova, T.; Akhtar, R.; Aldebert, F.; Althorpe, N.; Andersson, M.; Baldwin, A.; Bates, K.; Bhattacharyya, S.; Bower, L.; Browne, P.; et al. EMBL Nucleotide Sequence Database in 2006. Nucleic Acids Res. 2007, 35. [Google Scholar] [CrossRef] [PubMed]
  5. The UniProt Consortium. Update on activities at the Universal Protein Resource (UniProt) in 2013. Nucleic Acids Res. 2013, 41, D43–D47. [Google Scholar]
  6. Kodama, Y.; Mashima, J.; Kaminuma, E.; Gojobori, T.; Ogasawara, O.; Takagi, T.; Okubo, K.; Nakamura, Y. The DNA Data Bank of Japan Launches a New Resource, the DDBJ Omics Archive of Functional Genomics Experiments. Nucleic Acids Res. 2012, 40. [Google Scholar] [CrossRef] [PubMed]
  7. Stratton, M.R.; Campbell, P.J.; Andrew, P.F. The Cancer Genome. Nature 2009, 458, 719–724. [Google Scholar] [CrossRef]
  8. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; Da Fonseca, G.A.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853. [Google Scholar] [CrossRef]
  9. Vasconcelos, R.; Brito, J.C.; Carranza, S.; Harris, D.J. Review of the distribution and conservation status of the terrestrial reptiles of the Cape Verde Islands. Oryx 2013, 47, 77–87. [Google Scholar] [CrossRef] [Green Version]
  10. Romeiras, M.M.; Catarino, S.; Gomes, I.; Fernandes, C.; Costa, J.C.; Caujapé-Castells, J.; Duarte, M.C. IUCN Red List assessment of the Cape Verde endemic flora: Towards a global strategy for plant conservation in Macaronesia. Bot. J. Linn. Soc. 2016, 180, 413–425. [Google Scholar] [CrossRef]
  11. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; Mcgettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X Version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed]
  13. Farris, J.S. Methods for Computing Wagner Trees. Syst. Biol. 1970, 19, 83–92. [Google Scholar] [CrossRef]
  14. Fitch, W.M. Toward defining the course of evolution: Minimum change for a specific tree topology. Syst. Biol. 1971, 20, 406–416. [Google Scholar] [CrossRef]
  15. Carine, M.A.; Francisco-Ortega, J.; Santos-Guerra, A.; Russell, S.J. Relationships of island and continental floras: Molecular evidence for multiple colonisations into Macaronesia and subsequent back-colonisation of the continent in Convolvulus L. Am. J. Bot. 2004, 91, 1070–1085. [Google Scholar] [CrossRef] [PubMed]
  16. Francisco-Ortega, J.; Barber, J.; Santos-Guerra, A.; Febles-Hernández, R.; Jansen, R.K. Origin and evolution of the endemic genera of Gonosperminae (Asteraceae: Anthemideae) from the Canary Islands: Evidence from nucleotide sequences of the internal transcribed spacers of the nuclear ribosomal DNA. Am. J. Bot. 2001, 88, 161–169. [Google Scholar] [CrossRef] [PubMed]
  17. Helfgott, D.M.; Franciso-Ortega, F.; Santos-Guerra, A.; Jansen, R.K.; Simpson, B.B. Biogeography and breeding system evolution of the woody Bencomia alliance (Rosaceae) in Macaronesia based on ITS sequence data. Syst. Bot. 2000, 25, 82–97. [Google Scholar] [CrossRef]
  18. Panero, J.L.; Francisco-Ortega, J.; Jansen, R.K.; Santos-Guerra, A. Molecular evidence for multiple origins of woodiness and a New World biogeographic connection of the Macaronesian Island endemic Pericallis (Asteraceae: Senecioneae). Proc. Natl. Acad. Sci. USA 1999, 96, 13886–13891. [Google Scholar] [CrossRef]
  19. Dias, E.F.; Kilian, N.; Silva, L.; Schaefer, H.; Carine, M.; Rudall, P.J.; Santos-Guerra, A.; Moura, M. Phylogeography of the Macaronesian Lettuce Species Lactuca watsoniana and L. palmensis (Asteraceae). Biochem. Genet. 2018, 56, 315. [Google Scholar] [CrossRef]
  20. Fernández-Mazuecos, M.; Vargas, P. Genetically depauperate in the continent but rich in oceanic islands: Cistus monspeliensis (Cistaceae) in the Canary Islands. PLoS ONE 2011, 6, e17172. [Google Scholar] [CrossRef] [PubMed]
  21. Menezes, T.; Romeiras, M.M.; Sequeira, M.M.; Moura, M. Phylogenetic relationships and phylogeography of relevant lineages within the complex Campanulaceae family in Macaronesia. Ecol. Evol. 2017, 8, 88–108. [Google Scholar] [CrossRef] [Green Version]
  22. Cox, S.C.; Carranza, S.; Brown, R.P. Divergence Times and Colonization of the Canary Islands by Gallotia Lizards. Mol. Phylogenet. Evol. 2010, 56, 747–757. [Google Scholar] [CrossRef] [PubMed]
  23. Carranza, S.; Arnold, E.N. review of the geckos of the genus Hemidactylus (Squamata: Gekkonidae) from Oman based on morphology, mitochondrial and nuclear data, with descriptions of eight new species. Zootaxa 2012, 3378, 1–95. [Google Scholar] [CrossRef]
  24. Aldrich, J.R.A. Fisher and the Making of Maximum Likelihood 1912–1922. Stat. Sci. 1997, 12, 162–176. [Google Scholar] [CrossRef]
  25. Huelsenbeck, J.P.; Ronquist, F.; Nielsen, R.; Bollback, J.P. Bayesian inference of phylogeny and its impact on evolutionary biology. Science 2001, 294, 2310–2314. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, Z.; Rannala, B. Molecular Phylogenetics: Principles and Practice. Nat. Rev. Genet. 2012, 13, 303–314. [Google Scholar] [CrossRef]
  27. Yang, Z.; Rannala, B. Bayesian phylogenetic inference using DNA sequences: A Markov chain Monte Carlo Method. Mol. Biol. Evol. 1997, 14, 717–724. [Google Scholar] [CrossRef] [PubMed]
  28. Li, S.; Pearl, D.; Doss, H. Phylogenetic tree reconstruction using Markov chain Monte Carlo. J. Am. Stat. Assoc. 2000, 95, 493–508. [Google Scholar] [CrossRef]
  29. Swofford, D.L. PAUP*. Phylogenetic Analysis Using Parsimony; Version 4; Sinauer Associates: Sunderland, MA, USA, 2002. [Google Scholar] [CrossRef]
  30. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Franzke, A.; Samani, B.S.; Neuffer, B.; Mummenhoff, K.; Hurka, H. Molecular evidence in Diplotaxis (Brassicaceae) suggests a Quaternary origin of the Cape Verdean flora. Plant Syst. Evol. 2017, 303, 467–479. [Google Scholar] [CrossRef]
  32. Romeiras, M.M.; Paulo, O.S.; Duarte, M.C.; Pina-Martins, F.; Cotrim, M.H.; Carine, M.A.; Pais, M.S. Origin and diversification of the genus Echium (Boraginaceae) in the Cape Verde archipelago. Taxon 2011, 60, 1375–1385. [Google Scholar] [CrossRef]
  33. Thiv, M.; Thulin, M.; Hjertson, M.; Kropf, M.; Linder, H.P. Evidence for a vicariant origin of Macaronesian-Eritreo/Arabian disjunctions in Campylanthus Roth (Plantaginaceae). Mol. Phylogenet. Evol. 2010, 54, 607–616. [Google Scholar] [CrossRef] [PubMed]
  34. Vitales, D.; Garnatje, T.; Pellicer, J.; Vallès, J.; Santos-Guerra, A.; Sanmartín, I. The explosive radiation of Cheirolophus (Asteraceae, Cardueae) in Macaronesia. BMC Evol. Biol. 2014, 14, 118. [Google Scholar] [CrossRef] [PubMed]
  35. Carranza, S.; Arnold, E.N.; Geniez, Ph.; Roca, J.; Mateo, J.A. Radiation, Multiple Dispersal and Parallelism in the Skinks, Chalcides and Sphenops (Squamata: Scincidae), with Comments on Scincus and Scincopus and the Age of the Sahara Desert. Mol. Phylogenet. Evol. 2008, 46, 1071–1094. [Google Scholar] [CrossRef] [PubMed]
  36. Crottini, A.; Madsen, O.; Poux, C.; Strauss, A.; Vieites, D.R.; Vences, M. Vertebrate Time-Tree Elucidates the Biogeographic Pattern of a Major Biotic Change around the K-T Boundary in Madagascar. Proc. Natl. Acad. Sci. USA 2012, 109, 5358–5363. [Google Scholar] [CrossRef] [PubMed]
  37. Karin, B.R.; Metallinou, M.; Weinell, J.L.; Jackman, T.R.; Bauer, A.M. Resolving the Higher-Order Phylogenetic Relationships of the Circumtropical Mabuya Group (Squamata: Scincidae): An out-of-Asia Diversification. Mol. Phylogenet. Evol. 2016, 102, 220–232. [Google Scholar] [CrossRef]
  38. Metallinou, M.; Weinell, J.L.; Karin, B.R.; Conradie, W.; Wagner, P.; Schmitz, A.; Jackman, T.R.; Bauer, A.M. A Single Origin of Extreme Matrotrophy in African Mabuyine Skinks. Biol. Lett. 2016, 12, 20160430. [Google Scholar] [CrossRef]
  39. Miralles, A.; Vasconcelos, R.; Perera, A.; Harris, D.J.; Carranza, S. An Integrative Taxonomic Revision of the Cape Verdean Skinks (Squamata, Scincidae). Zoologica Scripta 2011, 40, 16–44. [Google Scholar] [CrossRef]
  40. Carranza, S.; Harris, D.J.; Arnold, E.N.; Batista, V.; de La Vega, J.P.G. Phylogeography of the Lacertid Lizard, Psammodromus algirus, in Iberia and across the Strait of Gibraltar. J. Biogeogr. 2006, 33, 1279–1288. [Google Scholar] [CrossRef]
  41. Agarwal, I.; Bauer, A.M.; Jackman, T.R.; Karanth, K.P. Insights into Himalayan Biogeography from Geckos: A Molecular Phylogeny of Cyrtodactylus (Squamata: Gekkonidae). Mol. Phylogenet. Evol. 2014, 80, 145–155. [Google Scholar] [CrossRef]
  42. Bansal, R.; Karanth, K.P. Phylogenetic Analysis and Molecular Dating Suggest That Hemidactylus anamallensis Is Not a Member of the Hemidactylus Radiation and Has an Ancient Late Cretaceous Origin. PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [PubMed]
  43. Carranza, S.; Arnold, E.N. Systematics, Biogeography, and Evolution of Hemidactylus Geckos (Reptilia: Gekkonidae) Elucidated Using Mitochondrial DNA Sequences. Mol. Phylogenet. Evol. 2006, 38, 531–545. [Google Scholar] [CrossRef] [PubMed]
  44. Gamble, T.; Bauer, A.M.; Colli, G.R.; Greenbaum, E.; Jackman, T.R.; Vitt, L.J.; Simons, A.M. Coming to America: Multiple Origins of New World Geckos: Origins of New World Geckos. J. Evol. Biol. 2011, 24, 231–244. [Google Scholar] [CrossRef] [PubMed]
  45. Garcia-Porta, J.; Morales, H.E.; Gómez-Díaz, E.; Sindaco, R.; Carranza, S. Patterns of Diversification in Islands: A Comparative Study across Three Gecko Genera in the Socotra Archipelago. Mol. Phylogenet. Evol. 2016, 98, 288–299. [Google Scholar] [CrossRef] [PubMed]
  46. Šmíd, J.; Carranza, S.; Kratochvíl, L.; Gvoždík, V.; Nasher, A.K.; Moravec, J. Out of Arabia: A Complex Biogeographic History of Multiple Vicariance and Dispersal Events in the Gecko Genus Hemidactylus (Reptilia: Gekkonidae). PLoS ONE 2013, 8. [Google Scholar] [CrossRef]
  47. Carranza, S.; Arnold, E.N.; Mateo, J.A.; Geniez, P. Relationships and Evolution of the North African Geckos, Geckonia and Tarentola (Reptilia: Gekkonidae), Based on Mitochondrial and Nuclear DNA Sequences. Mol. Phylogenet. Evol. 2002, 23, 244–256. [Google Scholar] [CrossRef]
  48. Rato, C.; Carranza, S.; Harris, D.J. Evolutionary History of the Genus Tarentola (Gekkota: Phyllodactylidae) from the Mediterranean Basin, Estimated Using Multilocus Sequence Data. BMC Evol. Biol. 2012, 12, 14. [Google Scholar] [CrossRef] [PubMed]
  49. Vasconcelos, R.; Carranza, S.; Harris, D.J. Insight into an Island Radiation: The Tarentola Geckos of the Cape Verde Archipelago. J. Biogeogr. 2010, 37, 1047–1060. [Google Scholar] [CrossRef]
  50. Posada, D. JModelTest: Phylogenetic Model Averaging. Mol. Biol. Evol. 2008, 25, 1253–1256. [Google Scholar] [CrossRef]
  51. Lanfear, R.; Frandsen, P.B.; Wright, A.M.; Senfeld, T.; Calcott, B. PartitionFinder 2: New Methods for Selecting Partitioned Models of Evolution for Molecular and Morphological Phylogenetic Analyses. Mol. Biol. Evol. 2016, msw260. [Google Scholar] [CrossRef]
  52. Stamatakis, A. RAxML Version 8: A Tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [PubMed]
  53. Huelsenbeck, J.P.; Ronquist, F.P. MRBAYES: Bayesian Inference of Phylogenetic Trees. Bioinformatics 2001, 17, 754–755. [Google Scholar] [CrossRef] [PubMed]
  54. Guindon, S.; Delsuc, F.; Dufayard, J.F.; Gascuel, O. Estimating Maximum Likelihood Phylogenies with PhyML. Methods Mol. Biol. 2009, 537, 113–137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Zwickl, D.J. Genetic Algorithm Approaches for the Phylogenetic Analysis of Large Biological Sequence Datasets under the Maximum Likelihood Criterion. Ph.D. Thesis, The University of Texas at Austin, Austin, TX, USA, 2006. [Google Scholar]
  56. Koutroumpa, K.; Theodoridis, S.; Warren, B.H.; Jiménez, A.; Celep, F.; Doğan, M.; Romeiras, M.M.; Santos-Guerra, A.; Fernández-Palacios, J.M.; Caujapé-Castells, J.; et al. An expanded molecular phylogeny of Plumbaginaceae, with emphasis on Limonium (sea lavenders): Taxonomic implications and biogeographic considerations. Ecol. Evol. 2018, 8, 12397–12424. [Google Scholar] [CrossRef] [PubMed]
  57. Bromham, L.; Duchêne, S.; Hua, X.; Ritchie, A.M.; Duchêne, D.A.; Ho, S.Y. Bayesian molecular dating: Opening up the black box. Biol. Rev. 2018, 93, 1165–1191. [Google Scholar] [CrossRef] [PubMed]
  58. Drummond, A.J.; Bouckaert, R.R. Bayesian Evolutionary Analysis with BEAST; Cambridge University Press: Cambridge, UK, 2015. [Google Scholar]
  59. Drummond, A.J.; Rambaut, A. BEAST: Bayesian Evolutionary Analysis by Sampling Trees. BMC Evol. Biol. 2007, 7, 214. [Google Scholar] [CrossRef] [PubMed]
  60. Bouckaert, R.; Heled, J.; Kühnert, D.; Vaughan, T.G.; Wu, C.-H.; Xie, D.; Suchard, M.; Rambaut, A.; Drummond, A.J. Beast2: A Software Platform for Bayesian Evolutionary Analysis. PLoS Comput. Biol. 2013, 10, 1003537. [Google Scholar] [CrossRef]
  61. Ogilvie, H.A.; Bouckaert, R.R.; Drummond, A.J. StarBEAST2 Brings Faster Species Tree Inference and Accurate Estimates of Substitution Rates. Mol. Biol. Evol. 2017, 34, 2101–2114. [Google Scholar] [CrossRef]
  62. Belfiore, N.M.; Liu, L.; Moritz, C. Multilocus Phylogenetics of a Rapid Radiation in the Genus Thomomys (Rodentia: Geomyidae). Syst. Biol. 2008, 57, 294–310. [Google Scholar] [CrossRef]
  63. Clement, M.; Posada, D.; Crandall, K.A. TCS: A Computer Program to Estimate Gene Genealogies. Mol. Ecol. 2010, 9, 1657–1659. [Google Scholar] [CrossRef]
  64. Huson, D.H.; Bryant, D. Application of Phylogenetic Networks in Evolutionary Studies. Mol. Biol. Evol. 2006, 23, 254–267. [Google Scholar] [CrossRef] [PubMed]
  65. Romeiras, M.M.; Monteiro, F.; Duarte, M.C.; Schaefer, H.; Carine, M. Patterns of genetic diversity in three plant lineages endemic to the Cape Verde Islands. AoB Plants 2015, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Schaefer, H.; Moura, M.; Belo Maciel, M.G.; Silva, L.; Rumsey, F.J.; Carine, M.A. The Linnean shortfall in oceanic island biogeography: A case study in the Azores. J. Biogeogr. 2011, 38, 1345–1355. [Google Scholar] [CrossRef]
  67. Romeiras, M.M.; Vieira, A.; Silva, D.N.; Moura, M.; Santos-Guerra, A.; Batista, D.; Duarte, M.C.; Paulo, O. Evolutionary and Biogeographic Insights on the Macaronesian Beta-Patellifolia Species (Amaranthaceae) from a Time-Scaled Molecular Phylogeny. PLoS ONE 2016, 11, e0152456. [Google Scholar] [CrossRef] [PubMed]
  68. Whittaker, R.J.; Fernández-Palacios, J.M.; Matthews, T.J.; Borregaard, M.K.; Triantis, K.A. Island biogeography: Taking the long view of nature’s laboratories. Science 2017, 357, eaam8326. [Google Scholar] [CrossRef]
  69. Warren, B.H.; Simberloff, D.; Ricklefs, R.E.; Aguilee, R.; Condamine, F.L.; Gravel, D.; Morlon, H.; Mouquet, N.; Rosindell, J.; Casquet, J.; et al. Islands as model systems in ecology and evolution: Prospects fifty years after MacArthur-Wilson. Ecol. Lett. 2015, 18, 200–217. [Google Scholar] [CrossRef] [PubMed]
  70. Price, J.P.; Otto, R.; Menezes de Sequeira, M.; Kueffer, C.; Schaefer, H.; Caujapé-Castells, J.; Fernández-Palacios, J.M. Colonization and diversification shape species–area relationships in three Macaronesian archipelagos. J. Biogeogr. 2018, 45, 2027–2039. [Google Scholar] [CrossRef]
  71. Alarcón, M.; Roquet, C.; García-Fernández, A.; Vargas, P.; Aldasoro, J.J. Phylogenetic and phylogeographic evidence for a Pleistocene disjunction between Campanula jacobaea (Cape Verde Islands) and C. balfourii (Socotra). Mol. Phylogen. Evol. 2013, 69, 828–836. [Google Scholar] [CrossRef]
  72. Ojeda, D.I.; Santos-Guerra, A.; Oliva-Tejera, F.; Jaen-Molina, R.; Caujapé-Castells, J.; Marrero-Rodríguez, Á.; Cronk, Q. DNA barcodes successfully identified Macaronesian Lotus (Leguminosae) species within early diverged lineages of Cape Verde and mainland Africa. AoB Plants 2014, 6, plu050. [Google Scholar] [CrossRef]
  73. Ávila, S.P.; Melo, C.; Berning, B.; Cordeiro, R.; Landau, B.; Marques, C. Persististrombus Coronatus (Mollusca: Strombidae) in the Lower Pliocene of Santa Maria Island (Azores, NE Atlantic): Paleoecology, Paleoclimatology and Paleobiogeographic Implications. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2016, 441, 912–923. [Google Scholar] [CrossRef]
  74. Ramalho, R.S. Building the Cape Verde Islands, 1st ed.; Springer: Berlin, Germany, 2011; p. 207. [Google Scholar] [CrossRef]
  75. Torres, P.; Silva, L.; Serralheiro, A.; Tassinari, C.; Munhá, J. Enquadramento geocronológico pelo método K/Ar das principais sequências vulcano-estratigráficas da Ilha do Sal—Cabo Verde. Garcia de Orta Serviços Geológicos 2002, 18, 9–13. [Google Scholar]
  76. Fernández-Palacios, J.M.; Rijsdijk, K.F.; Norder, S.J.; Otto, R.; de Nascimento, L.; Fernández-Lugo, S.; Tjørve, E.; Whittaker, R.J. Towards a glacial-sensitive model of island biogeography. Global Ecol. Biogeogr. 2016, 25, 817–830. [Google Scholar] [CrossRef]
  77. Weigelt, P.; Steinbauer, M.J.; Cabral, J.S.; Kreft, H. Late Quaternary climate change shapes island biodiversity. Nature 2016, 532, 99. [Google Scholar] [CrossRef] [PubMed]
  78. Ávila, S.P.; Cordeiro, R.; Madeira, P.; Silva, L.; Medeiros, A.; Rebelo, A.C.; Melo, C.; Neto, A.I.; Haroun, R.; Monteiro, A.; et al. Global Change Impacts on Large-Scale Biogeographic Patterns of Marine Organisms on Atlantic Oceanic Islands. Mar. Pollut. Bull. 2018, 126, 101–112. [Google Scholar] [CrossRef]
  79. Shaffer, H.B.; Pauly, G.B.; Oliver, J.C.; Trenham, P.C. The molecular phylogenetics of endangerment: Cryptic variation and historical phylogeography of the California tiger salamander, Ambystoma californiense. Mol. Ecol. 2004, 13, 3033–3049. [Google Scholar] [CrossRef]
  80. Jin, Y.; Yang, Z.; Brown, R.P.; Liao, P.; Liu, N. Intraspecific lineages of the lizard Phrynocephalus putjatia from the Qinghai-Tibetan Plateau: Impact of physical events on divergence and discordance between morphology and molecular markers. Mol. Phylog. Evol. 2014, 71, 288–297. [Google Scholar] [CrossRef] [PubMed]
  81. Jin, Y.; Liu, N.; Brown, R.P. The geography and timing of genetic divergence in the lizard Phrynocephalus theobaldi on the Qinghai-Tibetan plateau. Sci. Rep. 2017, 7, 2281. [Google Scholar] [CrossRef]
  82. Arnold, E.N.; Vasconcelos, R.; Harris, D.J.; Mateo, J.A.; Carranza, S. Systematics, biogeography and evolution of the endemic Hemidactylus geckos (Reptilia, Squamata, Gekkonidae) of the Cape Verde Islands: Based on morphology and mitochondrial and nuclear DNA sequences. Zool. Scr. 2008, 37, 619–636. [Google Scholar] [CrossRef]
  83. Vasconcelos, R.; Perera, A.; Geniez, P.; Harris, D.J.; Carranza, S. An integrative taxonomic revision of the Tarentola geckos (Squamata, Phyllodactylidae) of the Cape Verde Islands. Zool. J. Linn. Soc. 2002, 164, 328–360. [Google Scholar] [CrossRef]
  84. Pinho, C.; Santos, B.; Mata, V.; Seguro, M.; Romeiras, M.M.; Lopes, R.; Vasconcelos, R. What Is the Giant Wall Gecko Having for Dinner? Conservation Genetics for Guiding Reserve Management in Cabo Verde. Genes 2018, 9, 599. [Google Scholar] [CrossRef]
  85. Carranza, S.; Arnold, E.N. Lizards (Reptilia, Scincidae) crossed the Atlantic twice. Syst. Biol. 2003, 1, 275–282. [Google Scholar] [CrossRef]
  86. Datta-Roy, A.; Singh, M.; Srinivasulu, C.; Karanth, K.P. Phylogeny of the Asian Eutropis (Squamata: Scincidae) reveals an “into India” endemic Indian radiation. Mol. Phylog. Evol. 2012, 63, 817–824. [Google Scholar] [CrossRef] [PubMed]
  87. Mausfeld, P.; Schmitz, A.; Böhme, W.; Misof, B.; Vrcibradic, D.; Rocha, C.F.D. Phylogenetic affinities of Mabuya atlantica Schmidt, 1945, endemic to the Atlantic Ocean archipelago of Fernando de Noronha (Brazil): Necessity of partitioning the genus Mabuya Fitzinger, 1826 (Scincidae: Lygosominae). Zool. Anz. 2002, 241, 281–293. [Google Scholar] [CrossRef]
  88. Ho, S.Y.W.; Phillips, M.J. Accounting for Calibration Uncertainty in Phylogenetic Estimation of Evolutionary Divergence Times. Syst. Biol. 2009, 58, 367–380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Parham, J.F.; Donoghue, P.C.; Bell, C.J.; Calway, T.D.; Head, J.J.; Holroyd, P.A.; Inoue, J.G.; Irmis, R.B.; Jouce, W.G.; Ksepka, D.T.; et al. Best practices for justifying fossil calibrations. Syst. Biol. 2011, 61, 346–359. [Google Scholar] [CrossRef] [PubMed]
  90. Smith, S.A.; Walker, J.F.; Brown, J.; Walker-Hale, N. Nested phylogenetic conflicts, combinability, and deep phylogenomics in plants. bioRxiv 2018, 371930. [Google Scholar] [CrossRef]
  91. Wickett, N.J.; Mirarab, S.; Nguyen, N.; Warnow, T.; Carpenter, E.; Matasci, N.; Ayyampalayam, N. Phylotranscriptomic Analysis of the Origin and Early Diversification of Land Plants. Proc. Natl. Acad. Sci. USA 2014, 111, E4859–E4868. [Google Scholar] [CrossRef] [PubMed]
  92. Nater, A.; Burri, R.; Kawakami, T.; Smeds, L.; Ellegren, H. Resolving evolutionary relationships in closely related species with whole-genome sequencing data. Syst. Biol. 2015, 64, 1000–1017. [Google Scholar] [CrossRef] [PubMed]
  93. Mort, M.E.; Crawford, D.J.; Kelly, J.K.; Santos-Guerra, A.; Menezes de Sequeira, M.; Moura, M.; Caujapé-Castells, J. Multiplexed-shotgun-genotyping data resolve phylogeny within a very recently derived insular lineage. Am. J. Bot. 2015, 102, 634–641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Zheng, Y.; Wiens, J.J. Combining Phylogenomic and Supermatrix Approaches, and a Time-Calibrated Phylogeny for Squamate Reptiles (Lizards and Snakes) Based on 52 Genes and 4162 Species. Mol. Phylogenet. Evol. 2016, 94, 537–547. [Google Scholar] [CrossRef] [PubMed]
  95. Zong, D.; Gan, P.; Zhou, A.; Zhang, Y.; Zou, X.; Duan, A.; Song, Y.; He, C. Plastome Sequences Help to Resolve Deep-Level Relationships of Populus in the Family Salicaceae. Front. Plant Sci. 2019, 10, 5. [Google Scholar] [CrossRef] [PubMed]
  96. Twyford, A.D.; Ness, R.W. Strategies for complete plastid genome sequencing. Mol. Ecol. Resour. 2017, 17, 858–868. [Google Scholar] [CrossRef] [PubMed]
  97. Liang, X.-J.; Hinsinger, D.D.; Elias, R.B.; Strijk, J.S. The plastome sequence of Laurus azorica (Seub.) Franco, an endemic tree species of the Azores islands. Mitochondrial DNA Part B 2019, 4, 363–365. [Google Scholar] [CrossRef]
  98. Andrews, K.R.; Good, J.M.; Miller, M.R.; Luikart, G.; Hohenlohe, P.A. Harnessing the Power of RADseq for Ecological and Evolutionary Genomics. Nat. Rev. Genet. 2016, 2, 81–92. [Google Scholar] [CrossRef] [PubMed]
  99. Catchen, J.M.; Hohenlohe, P.A.; Bernatchez, L.; Funk, W.C.; Andrews, K.R.; Allendorf, F.W. Unbroken: RADseq Remains a Powerful Tool for Understanding the Genetics of Adaptation in Natural Populations. Mol. Ecol. Resour. 2017, 17, 362–365. [Google Scholar] [CrossRef] [PubMed]
  100. McKinney, G.J.; Larson, W.A.; Seeb, L.W.; Seeb, J.E. RAD seq provides unprecedented insights into molecular ecology and evolutionary genetics: Comment on Breaking RAD by Lowry et al. Mol. Ecol. Resour. 2017, 17, 356–361. [Google Scholar] [CrossRef] [PubMed]
  101. Ree, R.H.; Hipp, A.L. Inferring Phylogenetic History from Restriction Site Associated DNA (RADseq). In Next-Generation Sequencing in Plant Systematics, 1st ed.; Hörandl, E., Appelhans, M.S., Eds.; Koeltz Scientific Books: Oberreifenberg, Germany, 2015; p. 298. [Google Scholar] [CrossRef]
  102. Lowry, D.B.; Hoban, S.; Kelley, J.L.; Lotterhos, K.E.; Reed, L.K.; Antolin, M.F.; Storfer, A. Breaking RAD: An Evaluation of the Utility of Restriction Site-Associated DNA Sequencing for Genome Scans of Adaptation. Mol. Ecol. Resour. 2017, 17, 142–152. [Google Scholar] [CrossRef] [PubMed]
  103. Rodrigues, A.S.B.; Silva, S.E.; Pina-Martins, F.; Loureiro, J.; Castro, M.; Gharbi, K.; Johnson, K.P.; Dietrich, C.H.; Borges, P.A.V.; Quartau, J.A.; et al. Assessing Genotype-Phenotype Associations in Three Dorsal Colour Morphs in the Meadow Spittlebug Philaenus Spumarius (L.) (Hemiptera: Aphrophoridae) Using Genomic and Transcriptomic Resources. BMC Genet. 2016, 1–16. [Google Scholar] [CrossRef] [PubMed]
  104. Grewe, F.; Huang, J.-P.; Leavitt, S.D.; Lumbsch, H.T. Reference-Based RADseq Resolves Robust Relationships among Closely Related Species of Lichen-Forming Fungi Using Metagenomic DNA. Sci. Rep. 2017, 7, 9884. [Google Scholar] [CrossRef]
  105. Hohenlohe, P.A.; Bassham, S.; Etter, P.D.; Stiffler, N.; Johnson, E.A.; Cresko, W.A. Population Genomics of Parallel Adaptation in Threespine Stickleback Using Sequenced RAD Tags. PLoS Genet. 2010, 6, e1000862. [Google Scholar] [CrossRef]
  106. McCluskey, B.M.; Postlethwait, J.H. Phylogeny of Zebrafish, a “Model Species,” within Danio, a ‘Model Genus. Mol. Biol. Evol. 2015, 32, 635–652. [Google Scholar] [CrossRef] [PubMed]
  107. Tripp, E.A.; Tsai, Y.H.E.; Zhuang, Y.; Dexter, K.G. RADseq Dataset with 90% Missing Data Fully Resolves Recent Radiation of Petalidium (Acanthaceae) in the Ultra-Arid Deserts of Namibia. Ecol. Evol. 2017, 7, 7920–7936. [Google Scholar] [CrossRef] [PubMed]
  108. Brown, R.P.; Paterson, S.; Risse, J. Genomic Signatures of Historical Allopatry and Ecological Divergence in an Island Lizard. Genome Biol. Evol. 2016, 8, 3618–3626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  109. Ellegren, H. Genome sequencing and population genomics in non-model organisms. Trends Ecol. Evol. 2014, 29, 51–63. [Google Scholar] [CrossRef] [PubMed]
  110. Barley, A.J.; Monnahan, P.J.; Thomson, R.C.; Grismer, L.L.; Brown, R.M. Sun skink landscape genomics: Assessing the roles of micro-evolutionary processes in shaping genetic and phenotypic diversity across a heterogeneous and fragmented landscape. Mol. Ecol. 2015, 24, 1696–1712. [Google Scholar] [CrossRef] [PubMed]
  111. Wessinger, C.A.; Freeman, C.C.; Mort, M.E.; Rausher, M.D.; Hileman, L.C. Multiplexed Shotgun Genotyping Resolves Species Relationships within the North American Genus Penstemon. Am. J. Bot. 2016, 103, 912–922. [Google Scholar] [CrossRef]
  112. Eaton, D.A.R.; Ree, R.H. Inferring Phylogeny and Introgression Using RADseq Data: An Example from Flowering Plants (Pedicularis: Orobanchaceae). Syst. Biol. 2013, 62, 689–706. [Google Scholar] [CrossRef] [PubMed]
  113. Curto, M.; Schachtler, C.; Puppo, P.; Meimberg, H. Using a new RAD-sequencing approach to study the evolution of Micromeria in the Canary Islands. Mol. Phylogenet. Evol. 2018, 119, 160–169. [Google Scholar] [CrossRef]
  114. Peñalba, J.V.; Smith, L.L.; Tonione, M.A.; Sass, C.; Hykin, S.M.; Skipwith, P.L.; Mcguire, J.A.; Bowie, R.C.K.; Moritz, C. Sequence capture using PCR-generated probes: A cost-effective method of targeted high-throughput sequencing for nonmodel organisms. Mol. Ecol. Resour. 2014, 14, 1000–1010. [Google Scholar] [CrossRef]
  115. Bragg, J.G.; Potter, S.; Bi, K.; Moritz, C. Exon capture phylogenomics: Efficacy across scales of divergence. Mol. Ecol. Resour. 2016, 16, 1059–1068. [Google Scholar] [CrossRef]
  116. Brandley, M.C.; Bragg, J.G.; Singhal, S.; Chapple, D.G.; Jennings, C.K.; Lemmon, A.R.; Lemmon, E.M.; Thompson, M.B.; Moritz, C. Evaluating the performance of anchored hybrid enrichment at the tips of the tree of life: A phylogenetic analysis of Australian Eugongylus group scincid lizards. BMC Evol. Biol. 2015, 15, 62. [Google Scholar] [CrossRef] [PubMed]
Ijms 20 02782 g001 550
Figure 1. The five archipelagos of the Macaronesian region—Azores, Madeira, Selvagens, Cabo Verde and Canary Islands—with the corresponding distance to the closest mainland point (A), and the diversity of endemic Squamata reptile taxa within this region (B): Chalcides (coeruleopunctatus, sexlineatus, simonyi and viridanus), Chioninia (delalandii, fogoensis, nicolauensis, spinalis boavistensis, spinalis maioensis, spinalis salensis, spinalis santiagoensis, spinalis spinalis, stangeri, vaillanti vaillanti and vaillanti xanthotis), Gallotia (atlantica mahoratae, atlantica atlantica, atlantica laurae, auaritae, bravoana, caesaris caesaris, caesaris gomerae, galloti eisentrauti, galloti galloti, galloti insulanagae, galloti palmae, intermedia, simonyi machadoi and stehlini), Hemidactylus (bouvieri bouvieri, bouvieri ssp., bouvieri razoensis, boavistensis and lopezjuradoi), Tarentola (angustimentalis, bischoffi, boavistensis, bocagei, boettgeri boettgeri, boettgeri hierrensis, caboverdiana caboverdiana, caboverdiana raziana, caboverdiana substituta, darwini, delalandii, gigas gigas, gigas brancoensis, fogoensis, gomerensis, maioensis, nicolauensis, protogigas protogigas, protogigas hartogi and rudis), and Teira dugesii. The three extinct taxa (Chioninia coctei, Gallotia goliath and Gallotia simonyi simonyi) were not included in this representative analysis. Photo of Tarentola maioensis by R. Vasconcelos.
Figure 1. The five archipelagos of the Macaronesian region—Azores, Madeira, Selvagens, Cabo Verde and Canary Islands—with the corresponding distance to the closest mainland point (A), and the diversity of endemic Squamata reptile taxa within this region (B): Chalcides (coeruleopunctatus, sexlineatus, simonyi and viridanus), Chioninia (delalandii, fogoensis, nicolauensis, spinalis boavistensis, spinalis maioensis, spinalis salensis, spinalis santiagoensis, spinalis spinalis, stangeri, vaillanti vaillanti and vaillanti xanthotis), Gallotia (atlantica mahoratae, atlantica atlantica, atlantica laurae, auaritae, bravoana, caesaris caesaris, caesaris gomerae, galloti eisentrauti, galloti galloti, galloti insulanagae, galloti palmae, intermedia, simonyi machadoi and stehlini), Hemidactylus (bouvieri bouvieri, bouvieri ssp., bouvieri razoensis, boavistensis and lopezjuradoi), Tarentola (angustimentalis, bischoffi, boavistensis, bocagei, boettgeri boettgeri, boettgeri hierrensis, caboverdiana caboverdiana, caboverdiana raziana, caboverdiana substituta, darwini, delalandii, gigas gigas, gigas brancoensis, fogoensis, gomerensis, maioensis, nicolauensis, protogigas protogigas, protogigas hartogi and rudis), and Teira dugesii. The three extinct taxa (Chioninia coctei, Gallotia goliath and Gallotia simonyi simonyi) were not included in this representative analysis. Photo of Tarentola maioensis by R. Vasconcelos.
Ijms 20 02782 g001
Ijms 20 02782 g002 550
Figure 2. Number of studies (columns) and percentage of studied taxa (red line) within Macaronesian reptile lineages.
Figure 2. Number of studies (columns) and percentage of studied taxa (red line) within Macaronesian reptile lineages.
Ijms 20 02782 g002
Ijms 20 02782 g003 550
Figure 3. Calibration methods used with Macaronesian reptile lineages. Total number of calibration points (columns) and number of unresolved clades (red line). The average percentage of unresolved clades is shown within the columns as the number of unresolved clades/total number of clades.
Figure 3. Calibration methods used with Macaronesian reptile lineages. Total number of calibration points (columns) and number of unresolved clades (red line). The average percentage of unresolved clades is shown within the columns as the number of unresolved clades/total number of clades.
Ijms 20 02782 g003
Ijms 20 02782 g004 550
Figure 4. Divergence times obtained within Macaronesia reptile lineages in a meta-analysis of selected studies. Horizontal bars represent average age estimation (Mya), grey bars indicate minimum and maximum ages estimated (Mya). Asterisks indicate lack of data. Horizontal bar colour indicates the calibration method used: green–mutation rate; black–fossil; brown–mutation rate and geological. The age for the Fogo and Sal islands is shown in a dashed orange line (reference for age; [73]).
Figure 4. Divergence times obtained within Macaronesia reptile lineages in a meta-analysis of selected studies. Horizontal bars represent average age estimation (Mya), grey bars indicate minimum and maximum ages estimated (Mya). Asterisks indicate lack of data. Horizontal bar colour indicates the calibration method used: green–mutation rate; black–fossil; brown–mutation rate and geological. The age for the Fogo and Sal islands is shown in a dashed orange line (reference for age; [73]).
Ijms 20 02782 g004
Table
Table 1. Phylogenetic studies with Macaronesian reptiles, with molecular markers and software used, analyses conducted and corresponding bibliographic reference.
Table 1. Phylogenetic studies with Macaronesian reptiles, with molecular markers and software used, analyses conducted and corresponding bibliographic reference.
Family/GenusMarkersAnalysesSoftwareDatingYearRef.
Scincidae
Chalcidescyt b, 12S, 16SML+BIModelTest, MrBayes, PhyMLMutation rate (12S, cyt b) + Geological2008[35]
ChalcidesRAG1, BDNFBIMrBayes, BEASTMutation rate + Geological2012[36]
Chioninia16S, ND2, BDNF, BRCA1, BRCA2, CMOS, EXPH5, KIF24, MC1R, MXRA5, RAG1ML+BIPartitionFinder, RAxML, MrBayesMutation rate (16S, ND2)2016[37]
Chioninia16S, ND2, BDNF, BRCA1, BRCA2, CMOS, EXPH5, KIF24, MC1R, MXRA5, RAG1BIBEASTMutation rate (previous works)2016[38]
Chioniniacyt b, COI, 12SML+BIjModelTest, PhyML, MrBayesMutation rate (12S, cyt b)2011[39]
Lacertidae
Gallotiacyt b, 12S, 16SML, MPModelTest, PhyML, PAUPMutation rate (12S, cyt b) + Geological2006[40]
Gallotiacyt b, 12S, 16S, COIBI, MPMrBayes, TNT, BEASTGeological2010[22]
Gekkonidae
HemidactylusND2, RAG1, PDCML+BIPartitionFinder, RAxML, MrBayes, BEASTFossils2014[41]
HemidactylusCmos, 12S, RAG1, PDCML+BIModelTest, Paup, BEASTFossils + Mutation rate2013[42]
Hemidactyluscyt b, 12SMLModelTest, Paup, PhyML, MrBayesMutation rate2006[43]
Hemidactylus12S, cyt b, cmos, ND4, MC1R, RAG2ML+BIRAxML, MrBayesMutation rate (12S, cyt b)2012[23]
HemidactylusRAG1, RAG2, Cmos, ACM4, PDCMLRAxMLFossils2011[44]
Hemidactylus12S, ACM4, cmos, RAG1, RAG2, PDCBIPartitionFinder, BEASTFossils + Geological2016[45]
Hemidactylus12S, cyt b, MC1R, cmos, RAG1, RAG2ML+BIjModelTest, BEAST, RAxML, MrBayesMutation rate (12S, cyt b)2013[46]
Phyllodactylidae
Tarentolacyt b, 12S, cmosML, MPPaupGeological2002[47]
Tarentola12S, 16S, PDC, ACM4, MC1R, RAG2ML+BIjModelTest, RAxML, MrBayesMutation rate2012[48]
Tarentolacyt b, 12SML+BIjModelTest, MrBayes, PhyMLMutation rate (previous works)2010[49]

Share and Cite

MDPI and ACS Style

Romeiras, M.M.; Pena, A.R.; Menezes, T.; Vasconcelos, R.; Monteiro, F.; Paulo, O.S.; Moura, M. Shortcomings of Phylogenetic Studies on Recent Radiated Insular Groups: A Meta-Analysis Using Cabo Verde Biodiversity. Int. J. Mol. Sci. 2019, 20, 2782. https://doi.org/10.3390/ijms20112782

AMA Style

Romeiras MM, Pena AR, Menezes T, Vasconcelos R, Monteiro F, Paulo OS, Moura M. Shortcomings of Phylogenetic Studies on Recent Radiated Insular Groups: A Meta-Analysis Using Cabo Verde Biodiversity. International Journal of Molecular Sciences. 2019; 20(11):2782. https://doi.org/10.3390/ijms20112782

Chicago/Turabian Style

Romeiras, Maria M., Ana Rita Pena, Tiago Menezes, Raquel Vasconcelos, Filipa Monteiro, Octávio S. Paulo, and Mónica Moura. 2019. "Shortcomings of Phylogenetic Studies on Recent Radiated Insular Groups: A Meta-Analysis Using Cabo Verde Biodiversity" International Journal of Molecular Sciences 20, no. 11: 2782. https://doi.org/10.3390/ijms20112782

APA Style

Romeiras, M. M., Pena, A. R., Menezes, T., Vasconcelos, R., Monteiro, F., Paulo, O. S., & Moura, M. (2019). Shortcomings of Phylogenetic Studies on Recent Radiated Insular Groups: A Meta-Analysis Using Cabo Verde Biodiversity. International Journal of Molecular Sciences, 20(11), 2782. https://doi.org/10.3390/ijms20112782

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop