Pleistocene Landscape Dynamics Drives Lineage Divergence of a Temperate Freshwater Fish Gobio rivuloides in Coastal Drainages of Northern China

Understanding historical processes underlying lineage distribution patterns is a primary goal of phylogeography. We selected Gobio rivuloides (Cypriniformes: Gobionidae) as a model to improve our knowledge about how intraspecific genetic divergence of freshwater fishes arises in coastal drainages of northern China via statistical analysis using cytochrome b gene. The time-calibrated phylogeny of G. rivuloides showed the divergence of two major lineages (I and II) at ~0.98 Ma (million years ago). Lineage I can be divided into two sub-lineages (I-A and I-B) with a divergence time of ~0.83 Ma. Sub-lineage I-A inhabits the Amur River, and sub-lineage I-B lives in the Luan River and Liao River. Lineage II is distributed in the Yellow River and Hai River, with close genetic relationships between the two drainages, and can be split into two sub-lineages (II-C and II-D) with a divergence time of ~0.60 Ma. Our findings indicate that the splitting of lineages and sub-lineages could be attributed to geographic isolation caused by the formation of the Bohai Sea, river capture, and the episodic hydrologic closing of a paleolake during the late Lower-Middle Pleistocene. It is also the first report we know of displaying a clear phylogeographic break for freshwater fishes across coastal drainages in northern China.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Circulating Illness and Changes in Thermometer Use Behavior: Series of Cross-sectional Analyses

BACKGROUND

Temperature-taking behaviors vary with levels of circulating infectious illness; however, little is known about how these behaviors differ by demographic characteristics. Populations with higher perceived risks of illness are more likely to adopt protective health behaviors.

OBJECTIVE

We investigated differences in temperature-taking frequency and the proportion of readings that were feverish among demographic groups (age, gender, urban/rural status) over influenza offseason; influenza season; and waves 1, 2, and 3 of the COVID-19 pandemic.

METHODS

Using data from smart thermometers collected from May 1, 2019, to February 28, 2021, across the United States, we calculated the frequency of temperature-taking and the proportion of temperature readings that were feverish. Mixed-effects negative binomial and mixed-effects logistic regression analyses were performed to identify demographic characteristics associated with temperature-taking frequency and the proportion of feverish readings, respectively. Separate models were fit over five study periods: influenza offseason (n=122,480), influenza season (n=174,191), wave 1 of COVID-19 (n=350,385), wave 2 (n=366,489), and wave 3 (n=391,578).

RESULTS

Both temperature-taking frequency and the proportion of feverish readings differed by study period (ANOVA P<.001) and were the highest during influenza season. During all periods, children aged 2-5 years and 6-11 years had significantly higher frequencies of temperature-taking than users aged 19-30 years, and children had the highest proportion of feverish readings of all age groups, after adjusting for covariates. During wave 1 of COVID-19, users over the age of 60 years had 1.79 times (95% CI 1.76-1.83) the rate of temperature-taking as users aged 19-30 years and 74% lower odds (95% CI 72%-75%) of a reading being feverish. Across all periods, men had significantly lower temperature-taking frequency and significantly higher odds of having a feverish reading compared to women. Users living in urban areas had significantly higher frequencies of temperature-taking than rural users during all periods, except wave 2 of COVID-19, and urban users had higher odds of a reading being feverish in all study periods except wave 1 of COVID-19.

CONCLUSIONS

Temperature-taking behavior and the proportion of readings that were feverish are associated with both population disease levels and individual demographic characteristics. Differences in the health behavior of temperature-taking may reflect changes in both perceived and actual illness risk. Specifically, older adults may have experienced an increase in perceived risk during the first three waves of COVID-19, leading to increased rates of temperature monitoring, even when their odds of fever were lower than those of younger adults. Men's perceived risk of circulating infectious illnesses such as influenza and COVID-19 may be lower than that of women, since men took their temperature less frequently and each temperature had a higher odds of being feverish across all study periods. Infectious disease surveillance should recognize and incorporate how behavior impacts illness monitoring and testing.

Evidence of Závora Bay as a critical site for reef manta rays, Mobula alfredi, in southern Mozambique

The largest known reef manta ray (Mobula alfredi) population in Africa has been monitored for more than 20 years at several locations on the coast of the Inhambane Province in southern Mozambique. Nonetheless, before this study, little had been reported on the population dynamics of M. alfredi from Závora, a remote bay in the region. Photographic mark-recapture was used to investigate the size and structure of M. alfredi that aggregate at "Red Sands," a reef cleaning station in Závora Bay. An 11 year photographic data set was used to identify 583 M. alfredi individuals between 2010 and 2021. More than half of M. alfredi individuals were resighted at least once, with most encounters (up to 18 for one individual) occurring during the peak sighting period in July-November each year. An even sex ratio was observed, 44% females and 50% males, with no significant difference in resightings between the sexes. Pollock's robust design population models were used to estimate annual abundance, emigration, annual apparent survival and capture probability at Red Sands from July to November over a 6 year period (2016-2021). Abundance estimates varied year to year, ranging from 35 (95% c.i. [30, 45]) up to 233 (95% c.i. [224, 249]) M. alfredi individuals. Given the seasonal affinity of M. alfredi observed at Red Sands, this study highlights the importance of understanding fine-scale site use within the larger home range of this population to develop local management strategies.

Insights into Mus musculus population structure across Eurasia revealed by whole-genome analysis

For more than 100 years, house mice (Mus musculus) have been used as a key animal model in biomedical research. House mice are genetically diverse, yet their genetic background at the global level has not been fully understood. Previous studies have suggested that they originated in South Asia and diverged into three major subspecies, almost simultaneously, approximately 110,000–500,000 years ago; however, they have spread across the world with the migration of modern humans in prehistoric and historic times (∼10,000 years ago to the present day) and have undergone secondary contact, which has complicated the genetic landscape of wild house mice. In this study, we sequenced the whole-genome sequences of 98 wild house mice collected from Eurasia, particularly East Asia, Southeast Asia, and South Asia. Although wild house mice were found to consist of three major genetic groups corresponding to the three major subspecies, individuals representing admixtures between subspecies were more prevalent in East Asia than has been previously recognized. Furthermore, several samples exhibited an incongruent pattern of genealogies between mitochondrial and autosomal genomes. Using samples that likely retained the original genetic components of subspecies with the least admixture, we estimated the pattern and timing of divergence among the subspecies. The estimated divergence time of the three subspecies was 187,000–226,000 years ago. These results will help us to understand the genetic diversity of wild mice on a global scale, and the findings will be particularly useful in future biomedical and evolutionary studies involving laboratory mice established from such wild mice.

Coping with Environmental Extremes: Population Ecology and Behavioural Adaptation of Erebia pronoe, an Alpine Butterfly Species

Simple Summary

High alpine meadows are home to numerous endemic butterfly species. A combination of climate change and changes in agricultural practices has led to a severe decline in many species. A seemingly unaffected representative of this habitat is Erebia pronoe. We studied the behaviour, resource use and population structure of this species to explain its resilience and estimate its future survival potential. This species shows pronounced protandry in combination with serial eclosion. Males were significantly more active and mobile and were also caught significantly more often than females, resulting in a pronounced shift in sex ratio in the predicted population structure. The adults use a wide range of nectar plants and establish homeranges in areas of high habitat quality. Thus, Erebia pronoe adults use a wide array of resources combined with a slight specialisation to avoid niche overlap with closely related species. The resulting ecological flexibility seems to be an adaptation to unpredictable environmental conditions, which should be the result of a long-lasting adaptation process. Moreover, the combination of opportunism and modest specialisation should also be a good basis for coping with future changes caused by climate and land-use change.

Abstract

A mark-recapture study of the nominotypical Erebia pronoe in the Alps was conducted to survey its ecological demands and characteristics. Population structure analysis revealed a combination of protandry (one-week earlier eclosion of males) and serial eclosion. Significant differences between both sexes were found in population density (males: 580/ha ± 37 SE; females: 241/ha ± 66 SE), sex-ratio (2.4) and behaviour (57.7 vs. 11.9% flying). Both sexes used a wide range of nectar plants (Asteraceae, 77.3%; Dipsacaceae, 12.3%; Gentianaceae, 9.7%). The use of nectar plants shows a non-specific spectrum, which, however, completely avoids overlap with the locally co-occurring species Erebia nivalis. Movement patterns show the establishment of homeranges, which significantly limits the migration potential. Due to its broad ecological niche, E. pronoe will probably be able to react plastically to the consequences of climate change. The formation of high population densities, the unconcerned endangerment status, the unspecific resource spectrum and the sedentary character of the species make E. pronoe a potential indicator of the quality and general resource occurrence of alpine rupicolous grasslands.

Chromosome-level genome assembly and population genetic analysis of a critically endangered rhododendron provide insights into its conservation

Rhododendrons are woody plants, famous throughout the world as having high horticultural value. However, many wild species are currently threatened with extinction. Here, we report for the first time a high-quality, chromosome-level genome of Rhododendron griersonianum, which has contributed to approximately 10% of all horticultural rhododendron varieties but which in its wild form has been evaluated as critically endangered. The final genome assembly, which has a contig N50 size of approximately 34 M and a total length of 677 M, is the highest-quality genome sequenced within the genus to date, in part due to its low heterozygosity (0.18%). Identified repeats constitute approximately 57% of the genome, and 38 280 protein-coding genes were predicted with high support. We further resequenced 31 individuals of R. griersonianum as well as 30 individuals of its widespread relative R. delavayi, and performed additional conservation genomic analysis. The results showed that R. griersonianum had lower genetic diversity (θ = 2.58e-3; π = 1.94e-3) when compared not only to R. delavayi (θ = 11.61e-3, π = 12.97e-3), but also to most other woody plants. Furthermore, three severe genetic bottlenecks were detected using both the Stairway plot and fastsimcoal2 analysis, which are thought to have occurred in the late Middle Pleistocene and the Last Glacial Maximum (LGM) period. After these bottlenecks, R. griersonianum recovered and maintained a constant effective population size (>25 000) until now. Intriguingly, R. griersonianum has accumulated significantly more deleterious mutations in the homozygous state than R. delavayi, and several deleterious mutations (e.g., in genes involved in the response to heat stress) are likely to have harmed the adaptation of this plant to its surroundings. This high-quality, chromosome-level genome and the population genomic analysis of the critically endangered R. griersonianum will provide an invaluable resource as well as insights for future study in this species to facilitate conservation and in the genus Rhododendron in general.

Phylogeography and demographic expansion in the widely distributed horned passalus beetle, Odontotaenius disjunctus (coleoptera: Passalidae)

Dynamic climatic oscillations during the Pleistocene had profound effects on the distribution of species across North America. Although the role of historical climate change on speciation remains controversial, the impact on genetic variation within species has been well documented. Analyses of mtDNA sequences from the cytochrome oxidase I gene (911 bp) for 115 individuals of Odontotaenius disjunctus was combined with ecological niche modelling (ENM) to infer the demographic and population differentiation scenarios under present and past conditions. We inferred three lineages that diverged during the Pleistocene and replace each other geographically across the eastern United States. One of these lineages traverses previously identified genetic barriers for terrestrial animals including the Mississippi and Apalachicola Rivers and the Appalachian Mountains. We observed overlapping ranges between two haplotype groups as well as a region of secondary contact associated with ecological transition zone in northern Florida. The two continental lineages depict a genetic signature of a recent population increase associated with expanding niche envelope, whereas the clade restricted to peninsular Florida shows stable populations in a shrinking niche envelope. Given the lack of ecological separation, overlapping distribution of haplogroups and the presence of secondary contact zones, the taxonomic status of these lineages must await robust testing using multilocus DNA data to assess species boundaries.

Conservation Genomics of a Threatened Rhododendron: Contrasting Patterns of Population Structure Revealed From Neutral and Selected SNPs

Though it is well-acknowledged that next generation sequencing (NGS) technologies can provide further insights into plant conservation management than traditional molecular markers, studies employing NGS to address conservation genomics and subsequent conservation strategies for threatened plants are still rare. Rhododendron is the largest genus of woody plants in China, and many species are threatened, however, to date there has been no conservation genetic research using NGS in this genus. In the present study, we investigated the conservation genetics of R. cyanocarpum, a threatened species endemic to the Cangshan Mountains in Yunnan, China, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. Due to the availability of sufficient SNPs, we were able to distinguish between neutral and putatively selected SNPs and were able to further investigate the genetic diversity, population structure, and differentiation in R. cyanocarpum, as well as make an estimation of its demographic history. A total of 6,584 SNPs were obtained, of which 5,729 were neutral (detected using Tajima’s D). In terms of the 5,729 neutral SNPs, R. cyanocarpum had a higher genetic diversity (π = 0.0702 ± 0.0017, H_e = 0.0675 ± 0.0016) than other plant species assessed using Rad-seq methods, while population differentiation (F_st from 0.0314 to 0.0452) was weak. Interestingly, contrasting patterns of population structure were revealed from all neutral and selected SNPs, with distinct genetic clusters forming for all SNPs and neutral SNPs, but no distinct subgroups for selected ones. Moreover, we were able to detect changes in effective population size (N_e) of R. cyanocarpum from 150,000 years ago, including a bottleneck event ca. 60,000 years ago, followed by recovery of N_e over a short period, and a subsequent gradual decline in N_e to date. Implications for conserving R. cyanocarpum based on these main results are then discussed.

An estimator of the Opportunity for Selection that is independent of mean fitness

Variation among individuals in number of offspring (fitness, k) sets an upper limit to the evolutionary response to selection. This constraint is quantified by Crow's Opportunity for Selection (I), which is the variance in relative fitness (I = σ² _k /(u_k )² ). Crow's I has been widely used but remains controversial because it depends on mean offspring number in a sample ( k ¯ ). Here, I used a generalized Wright-Fisher model that allows for unequal probabilities of producing offspring to evaluate behavior of Crow's I and related indices under a wide range of sampling scenarios. Analytical and numerical results are congruent and show that rescaling the sample variance (s² _k ) to its expected value at a fixed k ¯ 2 removes dependence of I on mean offspring number, but the result still depends on choice of k ¯ 2 . A new index is introduced, Δ_I = Π- E(Î_drift ) = Π- 1/ k ¯ , which makes Î independent of sample k ¯ without the need for variance rescaling. Δ_I has a straightforward interpretation as the component of variance in relative fitness that exceeds that expected under a null model of random reproductive success. Δ_I can be used to directly compare estimates of the Opportunity for Selection for samples from different studies, different sexes, and different life stages.

Spatial genetic structure and historical demography of East Asian wild boar

Pleistocene climatic fluctuations may have had a profound impact on the evolutionary history of many species. The geographical pattern of European wild boar (Sus scrofa) is clearly studied, and it was greatly influenced by ancient climatic events, especially the Last Glacial Maximum. Previous research on genetic variation has mainly focused on the origin and distribution histories of domestic pigs. However, some questions have not been answered, including those concerning the genetic diversity, geographical pattern and possible historic influence of climate on East Asian wild boar (EAWB). Employing the control region of mtDNA (511 bp), we investigated the contributions of historic climate, which possibly shaped the genetic pattern of wild boar. Given that the level of genetic diversity of wild boars is higher in East Asia than in Europe, 172 haplotypes were detected from 680 individuals. Phylogenetic analysis demonstrated the complex phylogeographic structure of EAWB. Mismatch analysis, neutrality tests and the Bayesian Skyline Plot results all retrieved signals of a rapid population expansion, which might have played an important role in driving the formation of complex spatial genetic structure. Genetic data and species distribution modelling showed that the Last Glacial Maximum had weak effect on the distribution of the EAWB. We suggest that, in shaping spatial genetic structure in East Asian, long-term gene flow and population history played more important roles than Pleistocene climate fluctuations.

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

Individual fitness and the effects of a changing climate on the cessation and length of the breeding period using a 34-year study of a temperate songbird

Studies of the phenological responses of animals to climate change typically emphasize the initiation of breeding although climatic effects on the cessation and length of the breeding period may be as or more influential of fitness. We quantified links between climate, the cessation and length of the breeding period, and individual survival and reproduction using a 34-year study of a resident song sparrow (Melospiza melodia) population subject to dramatic variation in climate. We show that the cessation and length of the breeding period varied strongly across years, and predicted female annual fecundity but not survival. Breeding period length was more influential of fecundity than initiation or cessation of breeding alone. Warmer annual temperature and drier winters and summers predicted an earlier cessation of breeding. Population density, the date breeding was initiated, a female's history of breeding success, and the number of breeding attempts initiated previously also predicted the cessation of breeding annually, indicating that climatic, population, and individual factors may interact to affect breeding phenology. Linking climate projections to our model results suggests that females will both initiate and cease breeding earlier in the future; this will have opposite effects on individual reproductive rate because breeding earlier is expected to increase fecundity, whereas ceasing breeding earlier should reduce it. Identifying factors affecting the cessation and length of the breeding period in multiparous species may be essential to predicting individual fitness and population demography. Given a rich history of studies on the initiation of breeding in free-living species, re-visiting those data to estimate climatic effects on the cessation and length of breeding should improve our ability to predict the impacts of climate change on multiparous species.

Specialized or opportunistic-how does the high mountain endemic butterfly Erebia nivalis survive in its extreme habitats?

High mountain ecosystems are a challenge for the survival of animal and plant species, which have to evolve specific adaptations to cope with the prevailing extreme conditions. The strategies to survive may reach from opportunistic to highly adapted traits. One species successfully surviving under these conditions is the here studied butterfly Erebia nivalis. In a mark-release-recapture study performed in the Hohe Tauern National Park (Austria) from 22 July to 26 August 2013, we marked 1386 individuals and recaptured 342 of these. For each capture event, we recorded the exact point of capture and various other traits (wing conditions, behavior, nectar sources). The population showed a partial protandrous demography with the minority of males emerging prior to the females, but the majority being synchronized with them. Males and females differed significantly in their behavior with males being more flight active and females nectaring and resting more. Both sexes showed preferences for the same plant species as nectar sources, but this specialization apparently is the result of a rapid individual adaptation to the locally available flowers. Estimates of the realized dispersal distances predicted a comparatively high amount of long-distance flights, especially for females. Therefore, the adaptation of Erebia nivalis to the unpredictable high mountain conditions might be a mixture of opportunism and specialized traits.

Scale-dependent portfolio effects explain growth inflation and volatility reduction in landscape demography

Population demography is central to fundamental ecology and for predicting range shifts, decline of threatened species, and spread of invasive organisms. There is a mismatch between most demographic work, carried out on few populations and at local scales, and the need to predict dynamics at landscape and regional scales. Inspired by concepts from landscape ecology and Markowitz's portfolio theory, we develop a landscape portfolio platform to quantify and predict the behavior of multiple populations, scaling up the expectation and variance of the dynamics of an ensemble of populations. We illustrate this framework using a 35-y time series on gypsy moth populations. We demonstrate the demography accumulation curve in which the collective growth of the ensemble depends on the number of local populations included, highlighting a minimum but adequate number of populations for both regional-scale persistence and cross-scale inference. The attainable set of landscape portfolios further suggests tools for regional population management for both threatened and invasive species.

Population demographic history of a temperate shrub, Rhododendron weyrichii (Ericaceae), on continental islands of Japan and South Korea

Continental islands provide opportunities for testing the effects of isolation and migration on genetic variation in plant populations. In characteristic of continental islands is that the geographic connections between these islands, which are currently distinguished by seaways, have experienced fluctuations caused by sea‐level changes due to climate oscillations during the Quaternary. Plant populations on the islands have migrated between these islands via the exposed seafloors or been isolated. Here, we examined the demographic history of a temperate shrub, Rhododendron weyrichii, which is distributed in the southwestern parts of the Japanese archipelago and on an island of South Korea, using statistical phylogeographic approaches based on the DNA sequences of two chloroplast and eight nuclear loci in samples analyzed from 18 populations on eight continental islands, and palaeodistribution modeling. Time estimates for four island populations indicate that the durations of vicariance history are different between these populations, and these events have continued since the last glacial or may have predated the last glacial. The constancy or expansion of population sizes on the Japanese islands, and in contrast a bottleneck in population size on the Korean island Jeju, suggests that these islands may have provided different conditions for sustaining populations. The result of palaeodistribution modeling indicates that the longitudinal range of the species as a whole has not changed greatly since the last glacial maximum. These results indicate that exposed seafloors during the glacial period formed both effective and ineffective migration corridors. These findings may shed light on the effects of seafloor exposure on the migration of plants distributed across continental islands.

Coral reef fish populations can persist without immigration

Determining the conditions under which populations may persist requires accurate estimates of demographic parameters, including immigration, local reproductive success, and mortality rates. In marine populations, empirical estimates of these parameters are rare, due at least in part to the pelagic dispersal stage common to most marine organisms. Here, we evaluate population persistence and turnover for a population of orange clownfish, Amphiprion percula, at Kimbe Island in Papua New Guinea. All fish in the population were sampled and genotyped on five occasions at 2-year intervals spanning eight years. The genetic data enabled estimates of reproductive success retained in the same population (reproductive success to self-recruitment), reproductive success exported to other subpopulations (reproductive success to local connectivity), and immigration and mortality rates of sub-adults and adults. Approximately 50% of the recruits were assigned to parents from the Kimbe Island population and this was stable through the sampling period. Stability in the proportion of local and immigrant settlers is likely due to: low annual mortality rates and stable egg production rates, and the short larval stages and sensory capacities of reef fish larvae. Biannual mortality rates ranged from 0.09 to 0.55 and varied significantly spatially. We used these data to parametrize a model that estimated the probability of the Kimbe Island population persisting in the absence of immigration. The Kimbe Island population was found to persist without significant immigration. Model results suggest the island population persists because the largest of the subpopulations are maintained due to having low mortality and high self-recruitment rates. Our results enable managers to appropriately target and scale actions to maximize persistence likelihood as disturbance frequencies increase.

Influence of gene flow on divergence dating - implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species-tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation-with-migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene- and species-based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.

Climate, not Aboriginal landscape burning, controlled the historical demography and distribution of fire-sensitive conifer populations across Australia

Climate and fire are the key environmental factors that shape the distribution and demography of plant populations in Australia. Because of limited palaeoecological records in this arid continent, however, it is unclear as to which factor impacted vegetation more strongly, and what were the roles of fire regime changes owing to human activity and megafaunal extinction (since ca 50 kya). To address these questions, we analysed historical genetic, demographic and distributional changes in a widespread conifer species complex that paradoxically grows in fire-prone regions, yet is very sensitive to fire. Genetic demographic analysis showed that the arid populations experienced strong bottlenecks, consistent with range contractions during the Last Glacial Maximum (ca 20 kya) predicted by species distribution models. In southern temperate regions, the population sizes were estimated to have been mostly stable, followed by some expansion coinciding with climate amelioration at the end of the last glacial period. By contrast, in the flammable tropical savannahs, where fire risk is the highest, demographic analysis failed to detect significant population bottlenecks. Collectively, these results suggest that the impact of climate change overwhelmed any modifications to fire regimes by Aboriginal landscape burning and megafaunal extinction, a finding that probably also applies to other fire-prone vegetation across Australia.

Genetic structure and demographic history of the endangered tree species Dysoxylum malabaricum (Meliaceae) in Western Ghats, India: implications for conservation in a biodiversity hotspot

The impact of fragmentation by human activities on genetic diversity of forest trees is an important concern in forest conservation, especially in tropical forests. Dysoxylum malabaricum (white cedar) is an economically important tree species, endemic to the Western Ghats, India, one of the world's eight most important biodiversity hotspots. As D. malabaricum is under pressure of disturbance and fragmentation together with overharvesting, conservation efforts are required in this species. In this study, range-wide genetic structure of twelve D. malabaricum populations was evaluated to assess the impact of human activities on genetic diversity and infer the species’ evolutionary history, using both nuclear and chloroplast (cp) DNA simple sequence repeats (SSR). As genetic diversity and population structure did not differ among seedling, juvenile and adult age classes, reproductive success among the old-growth trees and long distance seed dispersal by hornbills were suggested to contribute to maintain genetic diversity. The fixation index (F_IS) was significantly correlated with latitude, with a higher level of inbreeding in the northern populations, possibly reflecting a more severe ecosystem disturbance in those populations. Both nuclear and cpSSRs revealed northern and southern genetic groups with some discordance of their distributions; however, they did not correlate with any of the two geographic gaps known as genetic barriers to animals. Approximate Bayesian computation-based inference from nuclear SSRs suggested that population divergence occurred before the last glacial maximum. Finally we discussed the implications of these results, in particular the presence of a clear pattern of historical genetic subdivision, on conservation policies.

General triallelic frequency spectrum under demographic models with variable population size

It is becoming routine to obtain data sets on DNA sequence variation across several thousands of chromosomes, providing unprecedented opportunity to infer the underlying biological and demographic forces. Such data make it vital to study summary statistics that offer enough compression to be tractable, while preserving a great deal of information. One well-studied summary is the site frequency spectrum-the empirical distribution, across segregating sites, of the sample frequency of the derived allele. However, most previous theoretical work has assumed that each site has experienced at most one mutation event in its genealogical history, which becomes less tenable for very large sample sizes. In this work we obtain, in closed form, the predicted frequency spectrum of a site that has experienced at most two mutation events, under very general assumptions about the distribution of branch lengths in the underlying coalescent tree. Among other applications, we obtain the frequency spectrum of a triallelic site in a model of historically varying population size. We demonstrate the utility of our formulas in two settings: First, we show that triallelic sites are more sensitive to the parameters of a population that has experienced historical growth, suggesting that they will have use if they can be incorporated into demographic inference. Second, we investigate a recently proposed alternative mechanism of mutation in which the two derived alleles of a triallelic site are created simultaneously within a single individual, and we develop a test to determine whether it is responsible for the excess of triallelic sites in the human genome.

Geological and ecological factors drive cryptic speciation of yews in a biodiversity hotspot

The interplay of orographic uplift and climatic changes in the Himalaya-Hengduan Mountains region (HHM) have had a key role in speciation and population demography. To gain further insight into these processes, we investigated their effects on Taxus wallichiana by combining molecular phylogeography and species distribution modeling. Molecular data were obtained from 43 populations of T. wallichiana. Nineteen climatic variables were analyzed alongside genetic discontinuities. Species distribution modeling was carried out to predict potential past distribution ranges. Two distinct lineages were identified, which diverged c. 4.2 (2.0-6.5) million years ago (Ma), a timescale that corresponds well with the recent uplift of the Qinghai-Tibet Plateau and subsequent climatic changes of the region. Correlations with climatic variables also suggest that ecological factors may have further reinforced the separation of the two lineages. Both lineages experienced population expansion during the last glaciation. The high genetic divergence, long-term isolation and ecological differentiation suggest a scenario of cryptic speciation in T. wallichiana associated with geological and climatic changes in the HHM. Our findings also challenge the notion of general population 'contraction' during the last glaciation in the HHM.

Conservation genomics of threatened animal species

The genomics era has opened up exciting possibilities in the field of conservation biology by enabling genomic analyses of threatened species that previously were limited to model organisms. Next-generation sequencing (NGS) and the collection of genome-wide data allow for more robust studies of the demographic history of populations and adaptive variation associated with fitness and local adaptation. Genomic analyses can also advance management efforts for threatened wild and captive populations by identifying loci contributing to inbreeding depression and disease susceptibility, and predicting fitness consequences of introgression. However, the development of genomic tools in wild species still carries multiple challenges, particularly those associated with computational and sampling constraints. This review provides an overview of the most significant applications of NGS and the implications and limitations of genomic studies in conservation.