A workflow of massive identification and application of intron markers using snakes as a model

Relative to the commonly used mitochondrial and nuclear protein‐coding genes, the noncoding intron sequences are a promising source of informative markers that have the potential to resolve difficult phylogenetic nodes such as rapid radiations and recent divergences. Yet many issues exist in the use of intron markers, which prevent their extensive application as conventional markers. We used the diverse group of snakes as an example to try paving the way for massive identification and application of intron markers. We performed a series of bioinformatics screenings which identified appropriate introns between single‐copy and conserved exons from two snake genomes, adding particular constraints on sequence length variability and sequence variability. A total of 1,273 candidate intron loci were retrieved. Primers for nested polymerase chain reaction (PCR) were designed for over a hundred candidates and tested in 16 snake representatives. 96 intron markers were developed that could be amplified across a broad range of snake taxa with high PCR successful rates. The markers were then applied to 49 snake samples. The large number of amplicons was subjected to next‐generation sequencing (NGS). An analytic strategy was developed to accurately recover the amplicon sequences, and approximately, 76% of the marker sequences were recovered. The average p‐distances of the intron markers at interfamily, intergenus, interspecies, and intraspecies levels were .168, .052, .015, and .004, respectively, suggesting that they were useful to study snake relationships of different evolutionary depths. A snake phylogeny was constructed with the intron markers, which produced concordant results with robust support at both interfamily and intragenus levels. The intron markers provide a convenient way to explore the signals in the noncoding regions to address the controversies on the snake tree. Our improved strategy of genome screening is effective and can be applied to other animal groups. NGS coupled with appropriate sequence processing can greatly facilitate the extensive application of molecular markers.     

Adaptive significance of food income in European snakes: body size is related to prey energetics

Feeding strategies and diet patterns have been extensively investigated in vertebrates and, more specifically, in snakes. Although it has been hypothesized that prey species may differ in terms of energy content, almost no theoretical or practical study has been carried out to determine actual nutritional values of the common prey types of wild snakes. Our model taxa were a selection of widely distributed and well known European snake species, which have all been studied in depth: approximately 76% of their diet is composed of mammals, reptiles, and insects. We therefore selected a single model species for each of these categories and proceeded with the analyses. Nutritional values were determined using a standard procedure: lizards and mice were richer in proteins than insects (crickets); insects and mice were richer in lipids than lizards, and mice and crickets have a higher energy content than lizards; lizards were rich in ashes. We then applied our experimental results to a selected sample of European terrestrial snakes (11 populations, ten species, seven genera, two families) characterized by different body size (50–160 cm total length) and reproductive strategies (oviparous versus viviparous), aiming to correlate these parameters with patterns of energy income. A direct relationship was found between body mass/body length ratio (BCI, body condition index) and meal energetics: the higher the BCI, the higher was the metabolic requirement, whereas BCI was independent of species or of reproductive system effect. Large‐sized snakes thus need a highly diversified and more energy‐rich diet than smaller snakes, supporting previous hypotheses. The simple applicability of this method could be of valuable support in further comparative research work, reducing experimental costs and stimulating further ecological, behavioural, and, possibly, phylogenetic comparisons. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 307–317.     

Comments on the evolution of the jaw adductor musculature of snakes

The aim of this study is to provide a general view of the adductor musculature of the alethinophidian snakes. The aponeurotic system present in anilioid snakes is here described as being also present in colubroid and booid snakes. Although modified in various groups, this aponeurotic system retains the same topographical pattern in the anilioids, booids and colubroids, and is thus hypothesized to be homologous. An analysis of the aponeurotic system and related muscular bundles within the alethinophidian snakes is given. A new terminology is proposed for the jaw adductor muscles where the muscles levator anguli oris and adductor mandibulae externus superficialis (proper) of snakes (sensu Lakjer, 1926; Haas, 1962) retain these names even if this fails to reflect the presumed homologies with the bundles of the same name in lizards (see Rieppel, 1988b); the fibres originating from the temporal tendon in the Anilioidea, and presumed to form a bundle of composite nature (Rieppel, 1980b), are named the M. adductor mandibulae externus temporalis (lost by the Macrostomata); the M. adductor mandibulae externus medialis is a composite muscle in the Anilioidea (Rieppel, 1980b) which give rise to two different muscles in the ‘booids’, the M. adductor mandibulae externus medialis, pars anterior and the M. adductor mandibulae externus profundus, the former being secondarily lost by the Caenophidia which retains only fibres homologues of the 3b and 3c heads of the profundus layer of lizards; the so-called M. adductor mandibular externus profundus of snakes (sensu Lackjer, 1926; Haas, 1962) is also a composite muscle in the Anilioidea (Rieppel, 1980b), in the alethinophidians it is essentially made of fibres homologous with the posterior pinnate part of the medialis layer of lizards, and is here named the M. adductor mandibulae externus medialis, pars posterior. As a result from this analysis it follows that: (1) the Macrostomata are characterized by the downward extension of the fibres forming the M. adductor mandibulae externus medialis, pars anterior and the loss of the M. adductor mandibulae externus temporalis: (2) the Xenopeltidae are set apart from the remaining macrostomatan snakes by the retention of the M. levator anguli oris and of a well developed lateral sheet of the quadrate aponeurosis; (3) the ‘booids’ form a monophyletic group comprising only the Boidae and Bolyeriidae (with the exclusion of the Xenopeltidae and Tropidophiidae) which is characterized by a differentiated M. adductor mandibulae externus medialis, pars anterior inserting on the lateral surface of the compound bone via its own aponeurosis; (4) the Tropidophiidae are set apart from all other snakes by the peculiar course of their lateral head vein; however, they belong to the Caenophidia as they show a facial carotid artery which passes dorsally to the mandibular and maxillary branches of the trigeminus; (5) a possible additional character in favour of an Acrochordoidea + Colubroidea monophyletic unit may be given by the pattern of innervation of the jaw adductor muscles in these two taxa; (6) a new interpretation of the compressor glandulae muscular complex of Atractaspis resulted in a morphologically similar pattern to that of the viperids; the phylogenetic implications of such similarity are discussed in detail.     

Identification of putative orthologous genes for the phylogenetic reconstruction of temperate woody bamboos (Poaceae: Bambusoideae)

The temperate woody bamboos (Arundinarieae) are highly diverse in morphology but lack a substantial amount of genetic variation. The taxonomy of this lineage is intractable, and the relationships within the tribe have not been well resolved. Recent studies indicated that this tribe could have a complex evolutionary history. Although phylogenetic studies of the tribe have been carried out, most of these phylogenetic reconstructions were based on plastid data, which provide lower phylogenetic resolution compared with nuclear data. In this study, we intended to identify a set of desirable nuclear genes for resolving the phylogeny of the temperate woody bamboos. Using two different methodologies, we identified 209 and 916 genes, respectively, as putative single copy orthologous genes. A total of 112 genes was successfully amplified and sequenced by next‐generation sequencing technologies in five species sampled from the tribe. As most of the genes exhibited intra‐individual allele heterozygotes, we investigated phylogenetic utility by reconstructing the phylogeny based on individual genes. Discordance among gene trees was observed and, to resolve the conflict, we performed a range of analyses using BUCKy and HybTree. While caution should be taken when inferring a phylogeny from multiple conflicting genes, our analysis indicated that 74 of the 112 investigated genes are potential markers for resolving the phylogeny of the temperate woody bamboos.     

The braincases of mosasaurs and Varanus, and the relationships of snakes

The braincase structure of the mosasaur Platecarpus is described in detail and compared to that of Varanus and snakes. The braincase of mosasaurs and Varanus is found to be closely similar in most respects other than the consequences of obliteration of the metakinetic axis in mosasaurs. Neither Varanus, nor mosasaurs, approach snakes in braincase structure. The hypothesis of a sister-group relationship of snakes and mosasauroids is discussed in the light of how hypotheses of homology, or synapomorphy, can be established on an empirical, i.e. testable and potentially falsifiable basis. The establishment of homology qua synapomorphy is recognized as a procedure involving at least two levels of generalization. The most basic level is the conjecture of similarity of constituent elements of two or more organisms. Such conjectures of similarity maintain their testability, and falsifiability, only if established by reference to topographical equivalence, or equivalence of connectivity.     

The evolutionary history of the ‘alba’ polymorphism in the butterfly subfamily Coliadinae (Lepidoptera: Pieridae)

Polymorphisms are common in the natural world and have played an important role in our understanding of how selection maintains multiple phenotypes within extant populations. Studying the evolutionary history of polymorphisms has revealed important features of this widespread form of phenotypic diversity, including its role in speciation, niche breadth, and range size. In the present study, we examined the evolutionary history of a ubiquitous colour polymorphism in the sulphur butterflies (subfamily: Coliadinae) termed the ‘alba’ polymorphism. We investigated the origin and stability of the ‘alba’ polymorphism using ancestral state reconstruction analysis. Our results indicate that the ancestor of the Coliadinae was polymorphic and that this polymorphism has undergone repeated transitions to monomorphism. Repeated loss of polymorphism suggests that the ‘alba’ polymorphism may be relatively unstable over evolutionary time. These results provide a framework for future studies on the origin and maintenance of the ‘alba’ polymorphism and guide the direction of future hypotheses. We discuss these results in light of current understandings of how the ‘alba’ polymorphism is maintained in extant populations.     

The evolutionary history of Afrocanarian blue tits inferred from genomewide SNPs

A common challenge in phylogenetic reconstruction is to find enough suitable genomic markers to reliably trace splitting events with short internodes. Here, we present phylogenetic analyses based on genomewide single‐nucleotide polymorphisms (SNPs) of an enigmatic avian radiation, the subspecies complex of Afrocanarian blue tits (Cyanistes teneriffae). The two sister species, the Eurasian blue tit (Cyanistes caeruleus) and the azure tit (Cyanistes cyanus), constituted the out‐group. We generated a large data set of SNPs for analysis of population structure and phylogeny. We also adapted our protocol to utilize degraded DNA from old museum skins from Libya. We found strong population structuring that largely confirmed subspecies monophyly and constructed a coalescent‐based phylogeny with full support at all major nodes. The results are consistent with a recent hypothesis that La Palma and Libya are relic populations of an ancient Afrocanarian blue tit, although a small data set for Libya could not resolve its position relative to La Palma. The birds on the eastern islands of Fuerteventura and Lanzarote are similar to those in Morocco. Together they constitute the sister group to the clade containing the other Canary Islands (except La Palma), in which El Hierro is sister to the three central islands. Hence, extant Canary Islands populations seem to originate from multiple independent colonization events. We also found population divergences in a key reproductive trait, viz. sperm length, which may constitute reproductive barriers between certain populations. We recommend a taxonomic revision of this polytypic species, where several subspecies should qualify for species rank.     

Molecular and morphological criteria for revision of the genus Microcoleus (Oscillatoriales,Cyanobacteria)

Ninety‐two strains of Microcoleus vaginatus (=nomenclatural‐type species of the genus Microcoleus Desmazières ex Gomont) and Phormidium autumnale Trevisan ex Gomont from a wide diversity of regions and biotopes were examined using a combination of morphological and molecular methods. Phylogenies based on the 16S rDNA and 16S‐23S ITS (partial) demonstrated that the 92 strains, together with a number of strains in GenBank, were members of a highly supported monophyletic clade of strains (Bayesian posterior probability = 1.0) distant from the species‐cluster containing the generitype of Phormidium. Similarity of the 16S rRNA gene exceeded 95.5% among all members of the Microcoleus clade, but was less than 95% between any Microcoleus strains and species outside of the clade (e.g., Phormidium sensu stricto). These findings, which are in agreement with earlier studies on these taxa, necessitate the revision of Microcoleus to include P. autumnale. Furthermore, the cluster of Phormidium species in the P. autumnale group (known as Group VII) must be moved into Microcoleus as well, and these nomenclatural transfers are included in this study. The main diacritical characters defining Microcoleus are related to the cytomorphology of trichomes, including: narrowed trichome ends, calyptra, cells shorter than wide up to more or less isodiametric, and facultative presence of sheaths. The majority of species are 4–10 μm in diameter. The possession of multiple trichomes in a common sheath is present facultatively in many but not all species.     

Adaptive radiations in butterflies: evolutionary history of the genus Erebia (Nymphalidae: Satyrinae)

We studied the speciose butterfly genus Erebia by reconstructing its phylogenetic relationships using parsimony and Bayesian approaches. We estimated times and rates of diversification for its lineages and employed a biogeographical analysis in order to reconstruct its evolutionary history. DNA sequence data from one mitochondrial gene and three nuclear genes were analyzed for a total of 74 species in Erebia. The estimated dates of origin and diversification for clades, in combination with a biogeographical analysis, suggest that the genus originated in Asian Russia and started its diversification process around 23 Myr. An important event was the dispersal of a lineage from Asia to Western Europe between 23 and 17 Myr, which allowed the radiation of most of species in the genus. The diversification pattern is consistent with a model of diversity limited by clade richness, which implies an early rapid diversification followed by deceleration due to a decrease in speciation. We argue that these characteristics of the evolutionary history of Erebia are consistent with a density‐dependent scenario, with species radiation limited by filling of niche space and reduced resources. We found that the Boeberia parmenio appears strongly supported in the genus Erebia and therefore we place Boeberia Prout, 1901 as a junior synonym of Erebia Dalman, 1816 ( syn. nov. ).     

Classification of crucigenoid algae: phylogenetic position of the reinstated genus Lemmermannia,Tetrastrum spp. Crucigenia tetrapedia,and C. lauterbornii (Trebouxiophyceae,Chlorophyta)1

The subfamily Crucigenioideae was traditionally classified within the well‐characterized family Scenedesmaceae (Chlorophyceae). Several morpho‐logical revisions and questionable taxonomic changes hampered the correct classification of crucigenoid species resulting in a high number of synonymous genera. We used a molecular approach to determine the phylogenetic position of several Tetrastrum and Crucigenia species. The molecular results were correlated with morphological and ontogenetic characters. Phylogenetic analyses of the SSU rDNA gene resolved the position of Tetrastrum heteracanthum and T. staurogeniaeforme as a new lineage within the Oocystis clade of the Trebouxiophyceae. Crucigenia tetrapedia, T. triangulare, T. punctatum, and T. komarekii were shown to be closely related to Botryococcus (Trebouxiophyceae) and were transferred to Lemmer‐mannia. Crucigenia lauterbornii was not closely related to the other Crucigenia strains, but was recovered within the Chlorella clade of the Trebouxiophyceae.     

Karyotype diversity and genome size variation in Neotropical Maxillariinae orchids

• Orchidaceae is a widely distributed plant family with very diverse vegetative and floral morphology, and such variability is also reflected in their karyotypes. However, since only a low proportion of Orchidaceae has been analysed for chromosome data, greater diversity may await to be unveiled. Here we analyse both genome size (GS) and karyotype in two subtribes recently included in the broadened Maxillariinea to detect how much chromosome and GS variation there is in these groups and to evaluate which genome rearrangements are involved in the species evolution.
• To do so, the GS (14 species), the karyotype – based on chromosome number, heterochromatic banding and 5S and 45S rDNA localisation (18 species) – was characterised and analysed along with published data using phylogenetic approaches.
• The GS presented a high phylogenetic correlation and it was related to morphological groups in Bifrenaria (larger plants – higher GS). The two largest GS found among genera were caused by different mechanisms: polyploidy in Bifrenaria tyrianthina and accumulation of repetitive DNA in Scuticaria hadwenii. The chromosome number variability was caused mainly through descending dysploidy, and x=20 was estimated as the base chromosome number.
• Combining GS and karyotype data with molecular phylogeny, our data provide a more complete scenario of the karyotype evolution in Maxillariinae orchids, allowing us to suggest, besides dysploidy, that inversions and transposable elements as two mechanisms involved in the karyotype evolution. Such karyotype modifications could be associated with niche changes that occurred during species evolution.

     

Evolution,biogeography and systematics of the western Palaearctic Zamenis ratsnakes

The phylogenetic relationships between western Palaearctic Zamenis and Rhinechis ratsnakes have been troubled, with recent estimates based on the supermatrix approach questioning their monophyly and providing contradictory results. In this study, we generated a comprehensive molecular data set for Zamenis and closely related ratsnakes to assess their phylogenetic and systematic relationships and infer their spatial and temporal modes of diversification. We obtained a fully resolved and well‐supported phylogeny, which is consistent across markers, taxon‐sets and phylogenetic methods. The close phylogenetic relationship between Rhinechis and Zamenis is well‐established. However, the early branching pattern within this clade, and the position of R. scalaris and Z. hohenackeri, remains poorly supported. The Persian ratsnake Z. persicus is sister to the Mediterranean species Z. situla, Z. longissimus and Z. lineatus, of which Z. situla is sister to a clade containing the latter two species. These results are consistent with a recent phylogenomic study on ratsnakes based on hundreds of loci. Whereas, topological tests based on our data and evidence from such phylogenomic study strongly rejected previous phylogenetic estimates based on the supermatrix approach and demonstrate that these “mega‐phylogenies”, with hundreds of taxa and high levels of missing data, have recovered inconsistent relationships with spurious nodal support. Biogeographical and molecular dating analyses suggest an origin of the ancestor of Rhinechis and Zamenis in the Aegean region with early cladogenesis during the Late Miocene associated with the Aegean arch formation and support a scenario of east‐to‐west diversification. Finally, while we have little morphological and phylogenetic evidence for the distinctiveness between Rhinechis and Zamenis, a classification of them in a single genus, and the designation of Zamenis scalaris (Schinz, 1822), reflects better their evolutionary relationships.     

Prospects on the evolutionary mitogenomics of plants: A case study on the olive family (Oleaceae)

The mitogenome is rarely used to reconstruct the evolutionary history of plants, contrary to nuclear and plastid markers. Here, we evaluate the usefulness of mitochondrial DNA for molecular evolutionary studies in Oleaceae, in which cases of cytoplasmic male sterility (CMS) and of potentially contrasted organelle inheritance are known. We compare the diversity and the evolution of mitochondrial and chloroplast genomes by focusing on the olive complex and related genera. Using high‐throughput techniques, we reconstructed complete mitogenomes (ca. 0.7 Mb) and plastomes (ca. 156 kb) for six olive accessions and one Chionanthus. A highly variable organization of mitogenomes was observed at the species level. In olive, two specific chimeric genes were identified in the mitogenome of lineage E3 and may be involved in CMS. Plastid‐derived regions (mtpt) were observed in all reconstructed mitogenomes. Through phylogenetic reconstruction, we demonstrate that multiple integrations of mtpt regions have occurred in Oleaceae, but mtpt regions shared by all members of the olive complex derive from a common ancestor. We then assembled 52 conserved mitochondrial gene regions and complete plastomes of ten additional accessions belonging to tribes Oleeae, Fontanesieae and Forsythieae. Phylogenetic congruence between topologies based on mitochondrial regions and plastomes suggests a strong disequilibrium linkage between both organellar genomes. Finally, while phylogenetic reconstruction based on plastomes fails to resolve the evolutionary history of maternal olive lineages in the Mediterranean area, their phylogenetic relationships were successfully resolved with complete mitogenomes. Overall, our study demonstrates the great potential of using mitochondrial DNA in plant phylogeographic and metagenomic studies.     

Patterns,biases and prospects in the distribution and diversity of Neotropical snakes

Motivation

We generated a novel database of Neotropical snakes (one of the world's richest herpetofauna) combining the most comprehensive, manually compiled distribution dataset with publicly available data. We assess, for the first time, the diversity patterns for all Neotropical snakes as well as sampling density and sampling biases.

Main types of variables contained

We compiled three databases of species occurrences: a dataset downloaded from the Global Biodiversity Information Facility (GBIF), a verified dataset built through taxonomic work and specialized literature, and a combined dataset comprising a cleaned version of the GBIF dataset merged with the verified dataset.

Spatial location and grain

Neotropics, Behrmann projection equivalent to 1° × 1°.

Time period

Specimens housed in museums during the last 150 years.

Major taxa studied

Squamata: Serpentes.

Software format

Geographical information system (GIS).

Results

The combined dataset provides the most comprehensive distribution database for Neotropical snakes to date. It contains 147,515 records for 886 species across 12 families, representing 74% of all species of snakes, spanning 27 countries in the Americas. Species richness and phylogenetic diversity show overall similar patterns. Amazonia is the least sampled Neotropical region, whereas most well‐sampled sites are located near large universities and scientific collections. We provide a list and updated maps of geographical distribution of all snake species surveyed.

Main conclusions

The biodiversity metrics of Neotropical snakes reflect patterns previously documented for other vertebrates, suggesting that similar factors may determine the diversity of both ectothermic and endothermic animals. We suggest conservation strategies for high‐diversity areas and sampling efforts be directed towards Amazonia and poorly known species.     

The phylogenetic distribution and morphological variation of the ‘pouch’ in female snakes

Siegel D.S., Miralles A, Trauth S.E. and Aldridge R.D. 2011. The phylogenetic distribution and morphological variation of the ‘pouch’ in female snakes. —Acta Zoologica (Stockholm) 93 : 400–408. The urodaeum of female snakes is a chamber that receives the oviducts, urinary ducts, and intestine. Previous studies have indicated that a derived region of unknown function persists between the urodaeum and oviducts in some snakes of the Colubroides. This structure was recently termed the ‘pouch’; however, it is also commonly referred to as the vaginal pouch. A broad sampling of female snakes revealed that the pouch evolved once on the branch leading to Colubroides. The presence of the pouch in Colubroides is not unambiguous, as this feature was lost in some members of the Atractaspididae and Lamprophiidae. Variation was also observed in pouch morphology between taxa of the Colubroides. Variation included the following: (1) whether the pouch was bifurcated and if so, were the arms of the pouch separated caudally, (2) whether separated pouch arms were split by a caudal bifurcation of the urodaeum or through medial septa of the pouch arms, and (3) the absence or presence of oviducal papillae invading the cranial/lateral extremity of the pouch. The variation in pouch morphology observed does not appear to be correlated with the evolutionary history of Colubroides and polymorphism within families is common.     

Renewed diversification is associated with new ecological opportunity in the Neotropical turtle ants

Ecological opportunity, defined as access to new resources free from competitors, is thought to be a catalyst for the process of adaptive radiation. Much of what we know about ecological opportunity, and the larger process of adaptive radiation, is derived from vertebrate diversification on islands. Here, we examine lineage diversification in the turtle ants (Cephalotes), a species‐rich group of ants that has diversified throughout the Neotropics. We show that crown group turtle ants originated during the Eocene (around 46 mya), coincident with global warming and the origin of many other clades. We also show a marked lineage‐wide slowdown in diversification rates in the Miocene. Contrasting this overall pattern, a species group associated with the young and seasonally harsh Chacoan biogeographic region underwent a recent burst of diversification. Subsequent analyses also indicated that there is significant phylogenetic clustering within the Chacoan region and that speciation rates are highest there. Together, these findings suggest that recent ecological opportunity, from successful colonization of novel habitat, may have facilitated renewed turtle ant diversification. Our findings highlight a central role of ecological opportunity within a successful continental radiation.