Speciation and gene flow across an elevational gradient in New Guinea kingfishers

Amazonian rivers are leaky barriers to gene flow in forest understory birds

Ever since Alfred Russel Wallace's nineteenth-century observation that related terrestrial species are often separated on opposing riverbanks, major Amazonian rivers have been recognized as key drivers of speciation. However, rivers are dynamic entities whose widths and courses may vary through time. It thus remains unknown how effective rivers are at reducing gene flow and promoting speciation over long timescales. We fit demographic models to genomic sequences to reconstruct the history of gene flow in three pairs of avian taxa fully separated by different Amazonian rivers, and whose geographic ranges do not make contact in headwater regions where rivers may cease to be barriers. Models with gene flow were best fit but still supported an initial period without any gene flow, which ranged from 187 000 to over 959 000 years, suggesting that rivers are capable of initiating speciation through long stretches of allopatric divergence. Allopatry was followed by either bursts or prolonged episodes of gene flow that retarded genomic differentiation but did not fully homogenize populations. Our results support Amazonian rivers as key barriers that promoted speciation and the build-up of species richness, but they also suggest that river barriers are often leaky, with genomic divergence accumulating slowly owing to episodes of substantial gene flow.

Biological clocks as age estimation markers in animals: a systematic review and meta-analysis

Various biological attributes associated with individual fitness in animals change predictably over the lifespan of an organism. Therefore, the study of animal ecology and the work of conservationists frequently relies upon the ability to assign animals to functionally relevant age classes to model population fitness. Several approaches have been applied to determining individual age and, while these methods have proved useful, they are not without limitations and often lack standardisation or are only applicable to specific species. For these reasons, scientists have explored the potential use of biological clocks towards creating a universal age-determination method. Two biological clocks, tooth layer annulation and otolith layering have found universal appeal. Both methods are highly invasive and most appropriate for post-mortem age-at-death estimation. More recently, attributes of cellular ageing previously explored in humans have been adapted to studying ageing in animals for the use of less-invasive molecular methods for determining age. Here, we review two such methods, assessment of methylation and telomere length, describing (i) what they are, (ii) how they change with age, and providing (iii) a summary and meta-analysis of studies that have explored their utility in animal age determination. We found that both attributes have been studied across multiple vertebrate classes, however, telomere studies were used before methylation studies and telomere length has been modelled in nearly twice as many studies. Telomere length studies included in the review often related changes to stress responses and illustrated that telomere length is sensitive to environmental and social stressors and, in the absence of repair mechanisms such as telomerase or alternative lengthening modes, lacks the ability to recover. Methylation studies, however, while also detecting sensitivity to stressors and toxins, illustrated the ability to recover from such stresses after a period of accelerated ageing, likely due to constitutive expression or reactivation of repair enzymes such as DNA methyl transferases. We also found that both studied attributes have parentally heritable features, but the mode of inheritance differs among taxa and may relate to heterogamy. Our meta-analysis included more than 40 species in common for methylation and telomere length, although both analyses included at least 60 age-estimation models. We found that methylation outperforms telomere length in terms of predictive power evidenced from effect sizes (more than double that observed for telomeres) and smaller prediction intervals. Both methods produced age correlation models using similar sample sizes and were able to classify individuals into young, middle, or old age classes with high accuracy. Our review and meta-analysis illustrate that both methods are well suited to studying age in animals and do not suffer significantly from variation due to differences in the lifespan of the species, genome size, karyotype, or tissue type but rather that quantitative method, patterns of inheritance, and environmental factors should be the main considerations. Thus, provided that complex factors affecting the measured trait can be accounted for, both methylation and telomere length are promising targets to develop as biomarkers for age determination in animals.

Time trees and clock genes: a systematic review and comparative analysis of contemporary avian migration genetics

Timing is a crucial aspect for survival and reproduction in seasonal environments leading to carefully scheduled annual programs of migration in many species. But what are the exact mechanisms through which birds (class: Aves) can keep track of time, anticipate seasonal changes, and adapt their behaviour? One proposed mechanism regulating annual behaviour is the circadian clock, controlled by a highly conserved set of genes, collectively called 'clock genes' which are well established in controlling the daily rhythmicity of physiology and behaviour. Due to diverse migration patterns observed within and among species, in a seemingly endogenously programmed manner, the field of migration genetics has sought and tested several candidate genes within the clock circuitry that may underlie the observed differences in breeding and migration behaviour. Among others, length polymorphisms within genes such as Clock and Adcyap1 have been hypothesised to play a putative role, although association and fitness studies in various species have yielded mixed results. To contextualise the existing body of data, here we conducted a systematic review of all published studies relating polymorphisms in clock genes to seasonality in a phylogenetically and taxonomically informed manner. This was complemented by a standardised comparative re-analysis of candidate gene polymorphisms of 76 bird species, of which 58 are migrants and 18 are residents, along with population genetics analyses for 40 species with available allele data. We tested genetic diversity estimates, used Mantel tests for spatial genetic analyses, and evaluated relationships between candidate gene allele length and population averages for geographic range (breeding- and non-breeding latitude), migration distance, timing of migration, taxonomic relationships, and divergence times. Our combined analysis provided evidence (i) of a putative association between Clock gene variation and autumn migration as well as a putative association between Adcyap1 gene variation and spring migration in migratory species; (ii) that these candidate genes are not diagnostic markers to distinguish migratory from sedentary birds; and (iii) of correlated variability in both genes with divergence time, potentially reflecting ancestrally inherited genotypes rather than contemporary changes driven by selection. These findings highlight a tentative association between these candidate genes and migration attributes as well as genetic constraints on evolutionary adaptation.

PAReTT: A Python Package for the Automated Retrieval and Management of Divergence Time Data from the TimeTree Resource for Downstream Analyses

Evolutionary processes happen gradually over time and are, thus, considered time dependent. In addition, several evolutionary processes are either adaptations to local habitats or changing habitats, otherwise restricted thereby. Since evolutionary processes driving speciation take place within the landscape of environmental and temporal bounds, several published studies have aimed at providing accurate, fossil-calibrated, estimates of the divergence times of both extant and extinct species. Correct calibration is critical towards attributing evolutionary adaptations and speciation both to the time and paleogeography that contributed to it. Data from more than 4000 studies and nearly 1,50,000 species are available from a central TimeTree resource and provide opportunities of retrieving divergence times, evolutionary timelines, and time trees in various formats for most vertebrates. These data greatly enhance the ability of researchers to investigate evolution. However, there is limited functionality when studying lists of species that require batch retrieval. To overcome this, a PYTHON package termed Python-Automated Retrieval of TimeTree data (PAReTT) was created to facilitate a biologist-friendly interaction with the TimeTree resource. Here, we illustrate the use of the package through three examples that includes the use of timeline data, time tree data, and divergence time data. Furthermore, PAReTT was previously used in a meta-analysis of candidate genes to illustrate the relationship between divergence times and candidate genes of migration. The PAReTT package is available for download from GitHub or as a pre-compiled Windows executable, with extensive documentation on the package available on GitHub wiki pages regarding dependencies, installation, and implementation of the various functions.

Role of paleoclimatic and paleohydrological processes in lineage divergence in freshwater organisms: A snippet from lentic genus Pila

The Indian subcontinent is extremely diverse in terms of its flora and fauna. However, only a handful of studies have aimed to understand the diversity of freshwater invertebrates using multiple lines of evidence in recent times. Here we aimed to estimate the cryptic diversity of two widespread freshwater snail species within the genus Pila (Röding, 1798) and uncover the processes behind lineage diversification in these species. We sequenced mitochondrial and nuclear markers from a comprehensive sampling of specimens from different river basins in India. We implemented an integrative taxonomy approach to delimit the lineages in these groups, employing phylogenetic, geometric morphometric and niche modelling-based methods. Then, we investigated the drivers of lineage divergence in these species using population genetic tools in conjunction with divergence time estimation. We found that both species consist of several genetically and ecologically distinct lineages. The genetic data showed that several of these lineages are restricted to a single or a few river basins. The divergence time estimation analyses indicated that the time frame of divergence within the species coincided with paleohydrological and paleoclimatic events in the Miocene. The diversification was primarily driven by allopatric isolation into different river basins. To conclude, the study sheds light on the complex interaction between the habitat preference of the species and the environment in shaping the diversification patterns in this group.

Repeated divergence of amphibians and reptiles across an elevational gradient in northern Madagascar

How environmental factors shape patterns of biotic diversity in tropical ecosystems is an active field of research, but studies examining the possibility of ecological speciation in terrestrial tropical ecosystems are scarce. We use the isolated rainforest herpetofauna on the Montagne d'Ambre (Amber Mountain) massif in northern Madagascar as a model to explore elevational divergence at the level of populations and communities. Based on intensive sampling and DNA barcoding of amphibians and reptiles along a transect ranging from ca. 470-1470 m above sea level (a.s.l.), we assessed a main peak in species richness at an elevation of ca. 1000 m a.s.l. with 41 species. The proportion of local endemics was highest (about 1/3) at elevations >1100 m a.s.l. Two species of chameleons (Brookesia tuberculata, Calumma linotum) and two species of frogs (Mantidactylus bellyi, M. ambony) studied in depth by newly developed microsatellite markers showed genetic divergence up the slope of the mountain, some quite strong, others very weak, but in each case with genetic breaks between 1100 and 1270 m a.s.l. Genetic clusters were found in transect sections significantly differing in bioclimate and herpetological community composition. A decrease in body size was detected in several species with increasing elevation. The studied rainforest amphibians and reptiles show concordant population genetic differentiation across elevation along with morphological and niche differentiation. Whether this parapatric or microallopatric differentiation will suffice for the completion of speciation is, however, unclear, and available phylogeographic evidence rather suggests that a complex interplay between ecological and allopatric divergence processes is involved in generating the extraordinary species diversity of Madagascar's biota. Our study reveals concordant patterns of diversification among main elevational bands, but suggests that these adaptational processes are only part of the complex of processes leading to species formation, among which geographical isolation is probably also important.

Multiple concordant cytonuclear divergences and potential hybrid speciation within a species complex in Asia

Various environmental factors impact the distribution, population structure, demography and evolutionary trajectory of a bird species, leading to genetic and morphological divergences between populations across its distribution. The Paradoxornis webbianus species complex is found throughout much of East Asia, where its geographically distinct populations exhibit dramatic morphological variation. This has resulted in a hotly debated taxonomy. This study intended to identify genetic divergence patterns and their underlying contributing factors for this species complex. We collected 243 birds, whose data was combined with those available in GenBank to perform phylogeographic analyses using one mitochondrial and six nuclear loci. Six mitochondrial clades were observed in the species complex, while individual-based Bayesian clustering using nuclear markers showed multiple congruent breaks. Overall, the six molecular lineages could be recognized as independent species under the lineage species concept in view of genetic divergence, clade-specific morphological changes and distribution: P. webbianus, P. w. bulomachus, P. alphonsianus, P. a. ganluoensis, P. brunneus brunneus and P. b. ricketti. The estimated divergence times range from 0.46 to 3.36 million years ago, suggesting it was likely impacted by paleoclimatic changes. Interestingly, P. alphonsianus carries two divergent mitochondrial lineages shared with P. webbianus and P. a. ganluoensis, respectively, and analyses based on nuclear loci found a similar pattern. We discussed the various hypotheses for this pattern and argued that P. alphonsianus was likely the result of hybridization between P. webbianus and P. a. ganluoensis. Further data on genome, transcriptome and breeding ecology are needed to address the hypothesis of hybrid speciation and its underlying mechanisms.

Fine-scale population genetic structure and barriers to gene flow in a widespread seabird (Ardenna pacifica)

Pelagic seabirds are highly mobile, reducing opportunities for population isolation that might promote differentiation and speciation. At the same time, many seabirds are philopatric, and their tendency to return to their natal islands to breed might reduce gene flow sufficiently to permit local adaptation and differentiation. To test the net impact of these competing processes, estimates of differentiation and gene flow based on comprehensive geographical sampling are required. We leveraged diverse source material to achieve comprehensive geographical sampling in a widespread seabird, the Wedge-tailed Shearwater (Ardenna pacifica). Using data from sequence capture and high-throughput sequencing of 2402 loci containing 20 780 single nucleotide polymorphisms, we tested for population differentiation and gene flow among breeding areas. We found little evidence of deep divergences within A. pacifica but were able to resolve fine-scale differentiation across island groups. This differentiation was sufficient to assign individuals sampled away from breeding areas to their likely source populations. Estimated effective migration surfaces revealed reduced migration between the Indian Ocean and Pacific Ocean, presumably owing to land barriers, and across the equatorial Pacific Ocean, perhaps associated with differences in breeding schedule. Our results reveal that, despite their mobility, pelagic seabirds can exhibit fine-scale population differentiation and reduced gene flow among ocean basins.

Further evidence from common garden rearing experiments of heritable traits separating lean and siscowet lake charr (Salvelinus namaycush) ecotypes

Genetic evidence of selection for complex and polygenically regulated phenotypes can easily become masked by neutral population genetic structure and phenotypic plasticity. Without direct evidence of genotype‐phenotype associations it can be difficult to conclude to what degree a phenotype is heritable or a product of environment. Common garden laboratory studies control for environmental stochasticity and help to determine the mechanism that regulate traits. Here we assess lipid content, growth, weight, and length variation in full and hybrid F₁ crosses of deep and shallow water sympatric lake charr ecotypes reared for nine years in a common garden experiment. Redundancy analysis (RDA) and quantitative‐trait‐loci (QTL) genomic scans are used to identify associations between genotypes at 19,714 single nucleotide polymorphisms (SNPs) aligned to the lake charr genome and individual phenotypes to determine the role that genetic inheritance plays in ecotype phenotypic diversity. Lipid content, growth, length, and weight differed significantly among lake charr crosses throughout the experiment suggesting that pedigree plays a large roll in lake charr development. Polygenic scores of 15 SNPs putatively associated with lipid content and/or condition factor indicated that ecotype distinguishing traits are polygenically regulated and additive. A QTL identified on chromosome 38 contained >200 genes, some of which were associated with lipid metabolism and growth, demonstrating the complex nature of ecotype diversity. The results of our common garden study further indicate that lake charr ecotypes observed in nature are predetermined at birth and that ecotypes differ fundamentally in lipid metabolism and growth.

Population diversification in the frog Mantidactylus bellyi on an isolated massif in northern Madagascar based on genetic, morphological, bioacoustic and ecological evidence

In the processes that give rise to new species, changes first occur at the population level. But with the continuous nature of the divergence process, change in biological properties delimiting the shift from “individuals of divergent populations” towards “individuals of distinct species”, as well as abiotic factors driving the change, remain largely ambivalent. Here we study diversification processes at the population level in a semi-aquatic frog, Mantidactylus (Brygoomantis) bellyi, across the diverse vegetation types of Montagne d’Ambre National Park (MANP), Madagascar. Genetic diversity was assessed with seven newly developed microsatellite markers as well as mitochondrial DNA sequences and concordance with patterns of ecological, morphological, and bioacoustic divergence evaluated. We found M. bellyi lacking mitochondrial differentiation within MANP, while microsatellite datasets partitioned them into three highly differentiated, geographically separated subpopulations (with indications for up to five subpopulations). The molecular grouping–primarily clustering individuals by geographic proximity–was coincident with differences in mean depth and width of waters, suggesting a possible role of fluvial characteristics in genetic exchange in this stream-breeding species. Genetic clustering not consistent with differences in call properties, except for dominant call frequencies under the two-subpopulations model. Morphological divergence was mostly consistent with the genetic clustering; subpopulations strongly differed by their snout-vent length, with individuals from high-elevation subpopulations smaller than those from populations below 1000 m above sea level. These results exemplify how mountains and environmental conditions might primarily shape genetic and morphological divergence in frog populations, without strongly affecting their calls.

The formation of avian montane diversity across barriers and along elevational gradients

Tropical mountains harbor exceptional concentrations of Earth's biodiversity. In topographically complex landscapes, montane species typically inhabit multiple mountainous regions, but are absent in intervening lowland environments. Here we report a comparative analysis of genome-wide DNA polymorphism data for population pairs from eighteen Indo-Pacific bird species from the Moluccan islands of Buru and Seram and from across the island of New Guinea. We test how barrier strength and relative elevational distribution predict population differentiation, rates of historical gene flow, and changes in effective population sizes through time. We find population differentiation to be consistently and positively correlated with barrier strength and a species' altitudinal floor. Additionally, we find that Pleistocene climate oscillations have had a dramatic influence on the demographics of all species but were most pronounced in regions of smaller geographic area. Surprisingly, even the most divergent taxon pairs at the highest elevations experience gene flow across barriers, implying that dispersal between montane regions is important for the formation of montane assemblages.

Landscape and climatic variations shaped secondary contacts amid barn owls of the Western Palearctic

The combined actions of climatic variations and landscape barriers shape the history of natural populations. When organisms follow their shifting niches, obstacles in the landscape can lead to the splitting of populations, on which evolution will then act independently. When two such populations are reunited, secondary contact occurs in a broad range of admixture patterns, from narrow hybrid zones to the complete dissolution of lineages. A previous study suggested that barn owls colonized the Western Palearctic after the last glaciation in a ring-like fashion around the Mediterranean Sea, and conjectured an admixture zone in the Balkans. Here, we take advantage of whole-genome sequences of 94 individuals across the Western Palearctic to reveal the complex history of the species in the region using observational and modeling approaches. Even though our results confirm that two distinct lineages colonized the region, one in Europe and one in the Levant, they suggest that it predates the last glaciation and identify a secondary contact zone between the two in Anatolia. We also show that barn owls recolonized Europe after the glaciation from two distinct glacial refugia: a previously identified western one in Iberia and a new eastern one in Italy. Both glacial lineages now communicate via eastern Europe, in a wide and permeable contact zone. This complex history of populations enlightens the taxonomy of Tyto alba in the region, highlights the key role played by mountain ranges and large water bodies as barriers and illustrates the power of population genomics in uncovering intricate demographic patterns.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.

Evolutionary conservatism will limit responses to climate change in the tropics

Rapid species turnover in tropical mountains has fascinated biologists for centuries. A popular explanation for this heightened beta diversity is that climatic stability at low latitudes promotes the evolution of narrow thermal tolerance ranges, leading to local adaptation, evolutionary divergence and parapatric speciation along elevational gradients. However, an emerging consensus from research spanning phylogenetics, biogeography and behavioural ecology is that this process rarely, if ever, occurs. Instead, closely related species typically occupy a similar elevational niche, while species with divergent elevational niches tend to be more distantly related. These results suggest populations have responded to past environmental change not by adapting and diverging in place, but instead by shifting their distributions to tightly track climate over time. We argue that tropical species are likely to respond similarly to ongoing and future climate warming, an inference supported by evidence from recent range shifts. In the absence of widespread in situ adaptation to new climate regimes by tropical taxa, conservation planning should prioritize protecting large swaths of habitat to facilitate movement.

Distinguishing genomic homogenization from parapatric speciation in an elevationally replacing pair of Ramphocelus tanagers

Geographically connected species pairs with weakly differentiated genomes could either represent cases of genomic homogenization in progress or of incipient parapatric speciation. Discriminating between these processes is difficult because intermediate stages of either may produce weakly differentiated genomes that diverge at few locations. We used coalescent modelling applied to a genome-wide sample of SNPs to discriminate between speciation with gene flow and genomic homogenization in two phenotypically distinct but genomically weakly diverged species of elevationally replacing Ramphocelus tanagers, forming a hybrid zone in the Andean foothills. We found overwhelming support for a model of genomic homogenization following secondary contact. Simulating under this model suggested that our species pair was differentiated (F_ST  = 0.30) at secondary contact but that most of the genome has rapidly homogenized during 254 Ky of high gene flow towards the present (F_ST  = 0.02). Despite extensive genome-wide homogenization, plumage remains distinctive with a narrower than expected geographic cline width, indicating divergent selection on colour. We found two SNPs significantly associated with plumage colour, which retain moderately high F_ST . We conclude that the majority of the genome has fused, but that divergent selection on select loci probably maintains the geographically structured colour differences between these incipient species.

The latitudinal taxonomy gradient

Emerging large-scale datasets coupled with statistical advances have provided new insights into the processes that generate the latitudinal diversity gradient (LDG). But many of these studies run into an old, if often underappreciated, problem: The interpretation of the data critically depends on the consistent application of criteria to define what constitutes a species. This is particularly pernicious for the LDG because good species have been easier to recognize in temperate than in tropical regions. We provide evidence that this latitudinal taxonomy gradient exists, discuss how this potentially impacts inferences about latitudinal variation in ecoevolutionary processes such as population differentiation and speciation, and provide a roadmap for how to mitigate taxonomic biases in the study of biodiversity patterns.

Mitogenomes Reveal Multiple Colonization of Mountains by Rattus in Sundaland

Tropical mountains are cradles of biodiversity and endemism. Sundaland, tropical Southeast Asia, hosts 3 species of Rattus endemic to elevations above 2000 m with an apparent convergence in external morphology: Rattus korinchi and R. hoogerwerfi from Sumatra, and R. baluensis from Borneo. A fourth one, R. tiomanicus, is restricted to lowland elevations across the whole region. The origins of these endemics are little known due to the absence of a robust phylogenetic framework. We use complete mitochondrial genomes from the 3 high altitude Rattus, and several related species to determine their relationships, date divergences, reconstruct their history of colonization, and test for selection on the mitochondrial DNA. We show that mountain colonization happened independently in Borneo (<390 Kya) and Sumatra (~1.38 Mya), likely from lowland lineages. The origin of the Bornean endemic R. baluensis is very recent and its genetic diversity is nested within the diversity of R. tiomanicus. We found weak evidence of positive selection in the high-elevation lineages and attributed the greater nonsynonymous mutations on these branches (specially R. baluensis) to lesser purifying selection having acted on the terminal branches in the phylogeny.

Museomics identifies genetic erosion in two butterfly species across the 20th century in Finland

Erosion of biodiversity generated by anthropogenic activities has been studied for decades and in many areas at the species level, using taxa monitoring. In contrast, genetic erosion within species has rarely been tracked, and is often studied by inferring past population dynamics from contemporaneous estimators. An alternative to such inferences is the direct examination of past genes, by analysing museum collection specimens. While providing direct access to genetic variation over time, historical DNA is usually not optimally preserved, and it is necessary to apply genotyping methods based on hybridization-capture to unravel past genetic variation. In this study, we apply such a method (i.e., HyRAD), to large time series of two butterfly species in Finland, and present a new bioinformatic pipeline, namely PopHyRAD, that standardizes and optimizes the analysis of HyRAD data at the within-species level. In the localities for which the data retrieved have sufficient power to accurately examine genetic dynamics through time, we show that genetic erosion has increased across the last 100 years, as revealed by signatures of allele extinctions and heterozygosity decreases, despite local variations. In one of the two butterflies (Erebia embla), isolation by distance also increased through time, revealing the effect of greater habitat fragmentation over time.