Estimating quantitative genetic parameters in wild populations: a comparison of pedigree and genomic approaches

The estimation of quantitative genetic parameters in wild populations is generally limited by the accuracy and completeness of the available pedigree information. Using relatedness at genomewide markers can potentially remove this limitation and lead to less biased and more precise estimates. We estimated heritability, maternal genetic effects and genetic correlations for body size traits in an unmanaged long‐term study population of Soay sheep on St Kilda using three increasingly complete and accurate estimates of relatedness: (i) Pedigree 1, using observation‐derived maternal links and microsatellite‐derived paternal links; (ii) Pedigree 2, using SNP‐derived assignment of both maternity and paternity; and (iii) whole‐genome relatedness at 37 037 autosomal SNPs. In initial analyses, heritability estimates were strikingly similar for all three methods, while standard errors were systematically lower in analyses based on Pedigree 2 and genomic relatedness. Genetic correlations were generally strong, differed little between the three estimates of relatedness and the standard errors declined only very slightly with improved relatedness information. When partitioning maternal effects into separate genetic and environmental components, maternal genetic effects found in juvenile traits increased substantially across the three relatedness estimates. Heritability declined compared to parallel models where only a maternal environment effect was fitted, suggesting that maternal genetic effects are confounded with direct genetic effects and that more accurate estimates of relatedness were better able to separate maternal genetic effects from direct genetic effects. We found that the heritability captured by SNP markers asymptoted at about half the SNPs available, suggesting that denser marker panels are not necessarily required for precise and unbiased heritability estimates. Finally, we present guidelines for the use of genomic relatedness in future quantitative genetics studies in natural populations.     

Complex patterns of genetic and phenotypic divergence in an island bird and the consequences for delimiting conservation units

Substantial phenotypic and genetic variation is often found below the species level and this may be useful in quantifying biodiversity and predicting future diversification. However, relatively few studies have tested whether different aspects of intraspecific variation show congruent patterns across populations. Here, we quantify several aspects of divergence between 13 insular populations of an island endemic bird, the Vanuatu white-eye ( Zosterops flavifrons ). The components of divergence studied are mitochondrial DNA (mtDNA), nuclear DNA microsatellites and morphology. These different aspects of divergence present subtly different scenarios. For instance, an mtDNA phylogenetic tree reveals a potential cryptic species on the most southerly island in Vanuatu and considerable divergence between at least two other major phylogroups. Microsatellite loci suggest that population genetic divergence between insular populations, both between and within phylogroups, is substantial, a result that is consistent with a low level of interisland gene flow. Finally, most populations were found to be strongly morphologically divergent, but no single population was morphologically diagnosable from all others. Taken together, our results show that, although many measures of divergence are concordant in this system, the number of divergent units identified varies widely depending on the characters considered and approach used. A continuum of divergence and a degree of discordance between different characters are both to be expected under simple models of evolution, but they present problems in terms of delimiting conservation units.     

History repeats itself: genomic divergence in copepods

Press stop, erase everything from now till some arbitrary time in the past and start recording life as it evolves once again. Would you see the same tape of life playing itself over and over, or would a different story unfold every time? The late Steven Jay Gould called this experiment replaying the tape of life and argued that any replay of the tape would lead evolution down a pathway radically different from the road actually taken (Gould 1989). This thought experiment has puzzled evolutionary biologists for a long time: how repeatable are evolutionary events? And if history does indeed repeat itself, what are the factors that may help us predict the path taken? A powerful means to address these questions at a small evolutionary scale is to study closely related populations that have evolved independently, under similar environmental conditions. This is precisely what Pereira et al. ( 2016 ) set out to do using marine copepods Tigriopus californicus, and present their results in this issue of Molecular Ecology. They show that evolution can be repeatable and even partly predictable, at least at the molecular level. As expected from theory, patterns of divergence were shaped by natural selection. At the same time, strong genetic drift due to small population sizes also constrained evolution down a similar evolutionary road, and probably contributed to repeatable patterns of genomic divergence.     

Erratum to: Recent range-wide demographic expansion in a Taiwan endemic montane bird,Steere's Liocichla (Liocichla steerii)

Background  

The subtropical island of Taiwan is an area of high endemism and a complex topographic environment. Phylogeographic studies indicate that vicariance caused by Taiwan's mountains has subdivided many taxa into genetic phylogroups. We used mitochondrial DNA sequences and nuclear microsatellites to test whether the evolutionary history of an endemic montane bird, Steere's Liocichla (Liocichla steerii), fit the general vicariant paradigm for a montane organism.     

Testing alternative mechanisms of evolutionary divergence in an African rain forest passerine bird

Abstract Models of speciation in African rain forests have stressed either the role of isolation or ecological gradients. Here we contrast patterns of morphological and genetic divergence in parapatric and allopatric populations of the Little Greenbul, Andropadus virens, within different and similar habitats. We sampled 263 individuals from 18 sites and four different habitat types in Upper and Lower Guinea. We show that despite relatively high rates of gene flow among populations, A. virens has undergone significant morphological divergence across the savanna-forest ecotone and mountain-forest boundaries. These data support a central component of the divergence-with-gene-flow model of speciation by suggesting that despite large amounts of gene flow, selection is sufficiently intense to cause morphological divergence. Despite evidence of isolation based on neutral genetic markers, we find little evidence of morphological divergence in fitness-related traits between hypothesized refugial areas. Although genetic evidence suggests populations in Upper and Lower Guinea have been isolated for over 2 million years, morphological divergence appears to be driven more by habitat differences than geographic isolation and suggests that selection in parapatry may be more important than geographic isolation in causing adaptive divergence in morphology.     

Niche expansion leads to small-scale adaptive divergence along an elevation gradient in a medium-sized passerine bird

Niche expansion can lead to adaptive differentiation and speciation, but there are few examples from contemporary niche expansions about how this process is initiated. We assess the consequences of a niche expansion by Mexican jays (Aphelocoma ultramarina) along an elevation gradient. We predicted that jays at high elevation would have straighter bills adapted to feeding on pine cones, whereas jays at low elevation would have hooked bills adapted to feeding on acorns. We measured morphological and genetic variation of 95 adult jays and found significant differences in hook length between elevations in accordance with predictions, a pattern corroborated by analysis at the regional scale. Genetic results from microsatellite and mtDNA variation support phenotypic differentiation in the presence of gene flow coupled with weak, but detectable genetic differentiation between high- and low-elevation populations. These results demonstrate that niche expansion can lead to adaptive divergence despite gene flow between parapatric populations along an elevation gradient, providing information on a key precursor to ecological speciation.     

Recent behavioural and population genetic divergence of an invasive ant in a novel environment

Aim Invasive species frequently exhibit high temporal and spatial variation in abundance. Although ecological aspects undoubtedly affect this variation, genetic factors may also play a part. The invasive unicolonial yellow crazy ant Anoplolepis gracilipes exhibits considerable variation in abundance throughout its extensive distribution in Australia’s Northern Territory, where it was first detected in the 1980s. First, we aimed to determine whether A. gracilipes variation in abundance was associated with behavioural and genetic differentiation of the population and to determine whether one or more introductions occurred. Second, we investigated whether the A. gracilipes population was genetically and behaviourally heterogeneous to determine whether population divergence has occurred since introduction. Location Tropical monsoonal savanna in Arnhem Land, Northern Territory, Australia. Methods Ant abundances were assessed at 13 sites throughout the study region. We used mitochondrial DNA sequences and microsatellite molecular markers to determine population genetic structure, which we correlated with abundance. Behavioural differentiation was assayed using aggression trials and analysed together with genetic data to investigate population divergence. Results Although we found considerable variation in abundance, we found no association between population structure and differences in abundance. Our analyses suggest that A. gracilipes ants in Arnhem Land resulted from a single introduction. The population is not homogeneous, however, as aggression scores varied over both genetic and geographic distance. We also found a positive relationship between genetic and geographic distance. Main conclusions The variation in abundance in the Arnhem Land population of A. gracilipes is clearly not owing to invasion by ants from different sources. The genetic and behavioural differentiation we observed is suggestive of incipient genetic and behavioural divergence, which may be expected over time when an invasive species enters in a new environment.     

Characterizing the mechanical contribution of fiber angular distribution in connective tissue: comparison of two modeling approaches

Modeling of connective tissues often includes collagen fibers explicitly as one of the components. These fibers can be oriented in many directions; therefore, several studies have considered statistical distributions to describe the fiber arrangement. One approach to formulate a constitutive framework for distributed fibers is to express the mechanical parameters, such as strain energy and stresses, in terms of angular integrals. These integrals represent the addition of the contribution of infinitesimal fractions of fibers oriented in a given direction. This approach leads to accurate results; however, it requires lengthy calculations. Recently, the use of generalized structure tensors has been proposed to represent the angular distribution in the constitutive equations of the fibers. Although this formulation is much simpler and fewer calculations are required, such structure tensors can only be used when all the fibers are in tension and the angular distribution is small. However, the amount of error introduced in these cases of non-tensile fiber loading and large angular distributions have not been quantified. Therefore, the objective of this study is to determine the range of values of angular distribution for which acceptable differences (less than 10%) between these two formulations are obtained. It was found, analytically and numerically, that both formulations are equivalent for planar distributions under equal-biaxial stretch. The comparison also showed, for other loading conditions, that the differences decrease when the fiber distribution is very small. Differences of less than 10% were usually obtained when the fiber distribution was very low (κ ≈ 0.03; κ ranges between 0 and 1/3, for aligned and isotropic distributed fibers, respectively). This range of angular distribution greatly limits the types of tissue that can be accurately analyzed using generalized structure tensors. It is expected that the results from this study guide the selection of a proper approach to analyze a particular tissue under a particular loading condition.     

Importance of genetic drift during Pleistocene divergence as revealed by analyses of genomic variation

Determining what factors affect the structuring of genetic variation is key to deciphering the relative roles of different evolutionary processes in species differentiation. Such information is especially critical to understanding how the frequent shifts and fragmentation of species distributions during the Pleistocene translates into species differences, and why the effect of such rapid climate change on patterns of species diversity varies among taxa. Studies of mitochondrial DNA (mtDNA) have detected significant population structure in many species, including those directly impacted by the glacial cycles. Yet, understanding the ultimate consequence of such structure, as it relates to how species divergence occurs, requires demonstration that such patterns are also shared with genomic patterns of differentiation. Here we present analyses of amplified fragment length polymorphisms (AFLPs) in the montane grasshopper Melanoplus oregonensis to assess the evolutionary significance of past demographic events and associated drift-induced divergence as inferred from mtDNA. As an inhabitant of the sky islands of the northern Rocky Mountains, this species was subject to repeated and frequent shifts in species distribution in response to the many glacial cycles. Nevertheless, significant genetic structuring of M. oregonensis is evident at two different geographic and temporal scales: recent divergence associated with the recolonization of the montane meadows in individual sky islands, as well as older divergence associated with displacements into regional glacial refugia. The genomic analyses indicate that drift-induced divergence, despite the lack of long-standing geographic barriers, has significantly contributed to species divergence during the Pleistocene. Moreover, the finding that divergence associated with past demographic events involves the repartitioning of ancestral variation without significant reductions of genomic diversity has intriguing implications - namely, the further amplification of drift-induced divergence by selection.     

A range-wide genetic bottleneck overwhelms contemporary landscape factors and local abundance in shaping genetic patterns of an alpine butterfly (Lepidoptera: Pieridae: Colias behrii)

Spatial and environmental heterogeneity are major factors in structuring species distributions in alpine landscapes. These landscapes have also been affected by glacial advances and retreats, causing alpine taxa to undergo range shifts and demographic changes. These nonequilibrium population dynamics have the potential to obscure the effects of environmental factors on the distribution of genetic variation. Here, we investigate how demographic change and environmental factors influence genetic variation in the alpine butterfly Colias behrii. Data from 14 microsatellite loci provide evidence of bottlenecks in all population samples. We test several alternative models of demography using approximate Bayesian computation (ABC), with the results favouring a model in which a recent bottleneck precedes rapid population growth. Applying independent calibrations to microsatellite loci and a nuclear gene, we estimate that this bottleneck affected both northern and southern populations 531–281 years ago, coinciding with a period of global cooling. Using regression approaches, we attempt to separate the effects of population structure, geographical distance and landscape on patterns of population genetic differentiation. Only 40% of the variation in F_ST is explained by these models, with geographical distance and least‐cost distance among meadow patches selected as the best predictors. Various measures of genetic diversity within populations are also decoupled from estimates of local abundance and habitat patch characteristics. Our results demonstrate that demographic change can have a disproportionate influence on genetic diversity in alpine species, contrasting with other studies that suggest landscape features control contemporary demographic processes in high‐elevation environments.     

Latitudinal divergence in a widespread amphibian: Contrasting patterns of neutral and adaptive genomic variation

Stochastic effects from demographic processes and selection are expected to shape the distribution of genetic variation in spatially heterogeneous environments. As the amount of genetic variation is central for long‐term persistence of populations, understanding how these processes affect variation over large‐scale geographical gradients is pivotal. We investigated the distribution of neutral and putatively adaptive genetic variation, and reconstructed demographic history in the moor frog (Rana arvalis) using 136 individuals from 15 populations along a 1,700‐km latitudinal gradient from northern Germany to northern Sweden. Using double digest restriction‐site associated DNA sequencing we obtained 27,590 single nucleotide polymorphisms (SNPs), and identified differentiation outliers and SNPs associated with growing season length. The populations grouped into a southern and a northern cluster, representing two phylogeographical lineages from different post‐glacial colonization routes. Hybrid index estimation and demographic model selection showed strong support for a southern and northern lineage and evidence of gene flow between regions located on each side of a contact zone. However, patterns of past gene flow over the contact zone differed between neutral and putatively adaptive SNPs. While neutral nucleotide diversity was higher along the southern than the northern part of the gradient, nucleotide diversity in differentiation outliers showed the opposite pattern, suggesting differences in the relative strength of selection and drift along the gradient. Variation associated with growing season length decreased with latitude along the southern part of the gradient, but not along the northern part where variation was lower, suggesting stronger climate‐mediated selection in the north. Outlier SNPs included loci involved in immunity and developmental processes.     

Decoupling of genetic and phenotypic divergence in a headwater landscape

In stream organisms, the landscape affecting intraspecific genetic and phenotypic divergence is comprised of two fundamental components: the stream network and terrestrial matrix. These components are known to differentially influence genetic structure in stream species, but to our knowledge, no study has compared their effects on genetic and phenotypic divergence. We examined how the stream network and terrestrial matrix affect genetic and phenotypic divergence in two stream salamanders, Gyrinophilus porphyriticus and Eurycea bislineata, in the Hubbard Brook Watershed, New Hampshire, USA. On the basis of previous findings and differences in adult terrestriality, we predicted that genetic divergence and phenotypic divergence in body morphology would be correlated in both species, but structured primarily by distance along the stream network in G. porphyriticus, and by overland distance in E. bislineata. Surprisingly, spatial patterns of genetic and phenotypic divergence were not strongly correlated. Genetic divergence, based on amplified DNA fragment length polymorphisms, increased with absolute geographic distance between sites. Phenotypic divergence was unrelated to absolute geographic distance, but related to relative stream vs. overland distances. In G. porphyriticus, phenotypic divergence was low when sites were close by stream distance alone and high when sites were close by overland distance alone. The opposite was true for E. bislineata. These results show that small differences in life history can produce large differences in patterns of intraspecific divergence, and the limitations of landscape genetic data for inferring phenotypic divergence. Our results also underscore the importance of explicitly comparing how terrestrial and aquatic conditions affect spatial patterns of divergence in species with biphasic life cycles.     

Genetic divergence and migration patterns in a North American passerine bird: implications for evolution and conservation

Like many other migratory birds, the black-throated blue warbler (Dendroica caerulescens) shows pronounced differences in migratory behaviour and other traits between populations: birds in the southern part of the breeding range have darker plumage and migrate to the eastern Caribbean during the winter, whereas those in the north have lighter plumage and migrate to the western Caribbean. We examined the phylogeography of this species, using samples collected from northern and southern populations, to determine whether differentiation between these populations dates to the Pleistocene or earlier, or whether differences in plumage and migratory behaviour have arisen more recently. We analysed variation at 369 bp of the mitochondrial control region domain I and also at seven nuclear microsatellites. Analyses revealed considerable genetic variation, but the vast majority of this variation was found within rather than between populations, and there was little differentiation between northern and southern populations. Phylogeographic analyses revealed a very shallow phylogenetic tree, a star-like haplotype network, and a unimodal mismatch distribution, all indicative of a recent range expansion from a single refugium. Coalescent modelling approaches also indicated a recent common ancestor for the entire group of birds analysed, no split between northern and southern populations, and high levels of gene flow. These results show that Pleistocene or earlier events have played little role in creating differences between northern and southern populations, suggesting that migratory and other differences between populations have arisen very recently. The implications of these results for the evolution of migration and defining taxonomic groups for conservation efforts are discussed.     

Reverse sex-biased philopatry in a cooperative bird: genetic consequences and a social cause

The genetic structure of a group or population of organisms can profoundly influence the potential for inbreeding and, through this, can affect both dispersal strategies and mating systems. We used estimates of genetic relatedness as well as likelihood-based methods to reconstruct social group composition and examine sex biases in dispersal in a Costa Rican population of white-throated magpie-jays ( Calocitta formosa , Swainson 1827), one of the few birds suggested to have female-biased natal philopatry. We found that females within groups were more closely related than males, which is consistent with observational data indicating that males disperse upon maturity, whereas females tend to remain in their natal territories and act as helpers. In addition, males were generally unrelated to one another within groups, suggesting that males do not disperse with or towards relatives. Finally, within social groups, female helpers were less related to male than female breeders, suggesting greater male turnover within groups. This last result indicates that within the natal group, female offspring have more opportunities than males to mate with nonrelatives, which might help to explain the unusual pattern of female-biased philopatry and male-biased dispersal in this system. We suggest that the novel approach adopted here is likely to be particularly useful for short-term studies or those conducted on rare or difficult-to-observe species, as it allows one to establish general patterns of philopatry and genetic structure without the need for long-term monitoring of identifiable individuals.     

The number of breeders explains genetic connectivity in an endangered bird

Connectivity is central to ecology and evolution as it focuses on the movement of individuals or genes across landscapes. Genetic connectivity approaches aim to understand gene flow but often estimate it indirectly based on metrics of genetic differentiation, which can also be affected by other evolutionary forces such as genetic drift. Gene flow and genetic drift are driven by separate ecological mechanisms with potentially differing effects on genetic differentiation and interpretations of genetic connectivity. The ecological mechanisms contributing to gene flow and genetic drift are primarily effective dispersal, or movement followed by successful reproduction, and the number of breeders in a local population, N_b, respectively. Yet, rarely are these ecological mechanisms and genetic connectivity measured simultaneously across landscapes. We examine the roles of effective dispersal and N_b on genetic connectivity across the entire range of the endangered snail kite (Rostrhamus sociabilis plumbeus), between 2006–2015. We find that both N_b and effective dispersal are important predictors of genetic connectivity across this landscape, but that N_b has a 3 × stronger effect on genetic connectivity. Furthermore, N_b is positively correlated with heterozygosity and allelic richness within patches, suggesting a potentially important role of genetic drift, in addition to gene flow, on genetic connectivity. These results emphasize that conservation efforts should focus on not only between‐patch processes of movement but also within‐patch processes regarding habitat quality and local population size for increasing genetic connectivity.     

Comparative biology comes into bloom: genomic and genetic comparison of flowering pathways in rice and Arabidopsis

Huge advances in plant biology are possible now that we have the complete genome sequences of several flowering plants. Now, genomes can be comprehensively compared and map-based cloning can be performed more easily. Association study is emerging as a powerful method for the functional identification of genes and molecular genetics has begun to reveal the basis of plant diversity. Taking the flowering pathways as an example, we discuss the potential of several approaches to comparative biology.