A new analytical approach to landscape genetic modelling: least-cost transect analysis and linear mixed models

Landscape genetics aims to assess the effect of the landscape on intraspecific genetic structure. To quantify interdeme landscape structure, landscape genetics primarily uses landscape resistance surfaces (RSs) and least-cost paths or straight-line transects. However, both approaches have drawbacks. Parameterization of RSs is a subjective process, and least-cost paths represent a single migration route. A transect-based approach might oversimplify migration patterns by assuming rectilinear migration. To overcome these limitations, we combined these two methods in a new landscape genetic approach: least-cost transect analysis (LCTA). Habitat-matrix RSs were used to create least-cost paths, which were subsequently buffered to form transects in which the abundance of several landscape elements was quantified. To maintain objectivity, this analysis was repeated so that each landscape element was in turn regarded as migration habitat. The relationship between explanatory variables and genetic distances was then assessed following a mixed modelling approach to account for the nonindependence of values in distance matrices. Subsequently, the best fitting model was selected using the statistic. We applied LCTA and the mixed modelling approach to an empirical genetic dataset on the endangered damselfly, Coenagrion mercuriale. We compared the results to those obtained from traditional least-cost, effective and resistance distance analysis. We showed that LCTA is an objective approach that identifies both the most probable migration habitat and landscape elements that either inhibit or facilitate gene flow. Although we believe the statistical approach to be an improvement for the analysis of distance matrices in landscape genetics, more stringent testing is needed.     

QTL and quantitative genetic analysis of beak morphology reveals patterns of standing genetic variation in an Estrildid finch

The intra- and interspecific diversity of avian beak morphologies is one of the most compelling examples for the power of natural selection acting on a morphological trait. The development and diversification of the beak have also become a textbook example for evolutionary developmental biology, and variation in expression levels of several genes is known to causally affect beak shape. However, until now, no genomic polymorphisms have been identified, which are related to beak morphology in birds. QTL mapping does reveal the location of causal polymorphisms, albeit with poor spatial resolution. Here, we estimate heritability and genetic correlations for beak length, depth and width and perform a QTL linkage analysis for these traits based on 1404 informative single-nucleotide polymorphisms genotyped in a four-generation pedigree of 992 captive zebra finches (Taeniopygia guttata). Beak size, relative to body size, was sexually dimorphic (larger in males). Heritability estimates ranged from 0.47 for beak length to 0.74 for beak width. QTL mapping revealed four to five regions of significant or suggestive genome-wide linkage for each of the three beak dimensions (nine different regions in total). Eight out of 11 genes known to influence beak morphology are located in these nine peak regions. Five QTL do not cover known candidates demonstrating that yet unknown genes or regulatory elements may influence beak morphology in the zebra finch.     

The probability of genetic parallelism and convergence in natural populations

Genomic and genetic methods allow investigation of how frequently the same genes are used by different populations during adaptive evolution, yielding insights into the predictability of evolution at the genetic level. We estimated the probability of gene reuse in parallel and convergent phenotypic evolution in nature using data from published studies. The estimates are surprisingly high, with mean probabilities of 0.32 for genetic mapping studies and 0.55 for candidate gene studies. The probability declines with increasing age of the common ancestor of compared taxa, from about 0.8 for young nodes to 0.1–0.4 for the oldest nodes in our study. Probability of gene reuse is higher when populations begin from the same ancestor (genetic parallelism) than when they begin from divergent ancestors (genetic convergence). Our estimates are broadly consistent with genomic estimates of gene reuse during repeated adaptation to similar environments, but most genomic studies lack data on phenotypic traits affected. Frequent reuse of the same genes during repeated phenotypic evolution suggests that strong biases and constraints affect adaptive evolution, resulting in changes at a relatively small subset of available genes. Declines in the probability of gene reuse with increasing age suggest that these biases diverge with time.     

mmod: an R library for the calculation of population differentiation statistics

mmod is a library for the R programming language that allows the calculation of the population differentiation measures D_est, G″_ST and φ′_ST. R provides a powerful environment in which to conduct and record population genetic analyses but, at present, no R libraries provide functions for the calculation of these statistics from standard population genetic files. In addition to the calculation of differentiation measures, mmod can produce parametric bootstrap and jackknife samples of data sets for further analysis. By integrating with and complimenting the existing libraries adegenet and pegas , mmod extends the power of R as a population genetic platform.     

Comparison of diverse affinity based high-abundance protein depletion strategies for improved bio-marker discovery in oral fluids

Oral fluids (OF) have been suggested as a source of biomarkers for oral and systemic health, but as with other bio-fluids, the presence of high-abundance proteins interferes with the detection of lower-abundance biomarkers. Here, we compared the performance of four depletion treatments: triple depletion (TD) of amylases, albumins and immunoglobulin G; multiple depletion (MD) of amylases and a panel of 20 proteins, a combination of the two (EMD) and combinatorial peptide ligand library based depletion termed CPLL. TD, MD, EMD and CPLL removed 76%, 83%, 85% and 94% of total proteins, respectively, coupled with increased low abundance protein detection and narrowed dynamic range. 2-DE revealed that all depletion pretreatments successfully clarified areas hampered by high-abundance proteins; however, EMD and CPLL exposed the highest number of proteins. Quantitative MS of EMD samples relative to none treated samples indicated that most of downregulated proteins (>90%) were EMD target proteins. In conclusion, a multiple step EMD and CPLL depletion approaches bring about the highest number of protein detection ability and the best hampered-area clearance. As CPLL requires at least 10 fold more protein starting material, we suggest EMD pretreatment as a new detection tool in instances of low protein starting material.     

Age‐dependent skewing of X chromosome inactivation appears delayed in centenarians’ offspring. Is there a role for allelic imbalance in Healthy Aging and Longevity?

Recently, it has been proposed that age-related X chromosome inactivation (XCI) skewing can clinically result in late-onset X-linked disorders. This observation leads to hypothesize that age-related skewed XCI might also influence lifespan in women. To investigate this issue, we employed a new experimental model of longevity and healthy aging including 55 female centenarians, 40 of their offspring, 33 age-matched offspring of both non-long-lived parents and 41 young women. Peripheral blood DNA from 169 females was screened for heterozygosity at the HUMARA locus. We confirmed that skewing of XCI is an age-dependent phenomenon. However, skewed XCI was significantly less severe and frequent in centenarians' offspring [degree of skewing (DS) = 0.16 ± 0.02] compared to age-matched offspring of both non-long-lived parents (DS = 0.24 ± 0.02) (P < 0.05). A second goal was to assess whether changes in XCI pattern could be a consequence of loss of methylation on X chromosome. Using a methylation array evaluating 1085 CpG sites across X chromosome and eleven CpG sites located at HUMARA locus, no differences in methylation levels and profiles emerged between all groups analysed, thus suggesting that age-associated epigenetic changes could not influence HUMARA results. In conclusion, the results presented herein highlight for the first time an interesting link between skewing of XCI and healthy aging and longevity. We speculate that the allelic imbalance produced by XCI skewing may compromise the cooperative and compensatory organization occurring between the two cell populations that make up the female mosaic.     

Species delimitation and digit number in a North African skink

Delimitation of species is an important and controversial area within evolutionary biology. Many species boundaries have been defined using morphological data. New genetic approaches now offer more objective evaluation and assessment of the reliability of morphological variation as an indicator that speciation has occurred. We examined geographic variation in morphology of the continuously distributed skink Chalcides mionecton from Morocco and used Bayesian analyses of nuclear and mitochondrial DNA (mtDNA) loci to examine: (i) their concordance with morphological patterns, (ii) support for species delimitation, (iii) timing of speciation, and (iv) levels of gene flow between species. Four digit individuals were found at sites between Cap Rhir (in the south) and the northern extreme of the range, whereas five‐digit individuals were found in two disjunct areas: (i) south of Cap Rhir and (ii) the north of the range where they were often syntopic with four‐digit individuals. The pattern of variation in generalized body dimensions was largely concordant with that in digit number, suggesting two general morphotypes. Bayesian analyses of population structure showed that individuals from sites south of Cap Rhir formed one genetic cluster, but that northern four‐ and five‐digit individuals clustered together. Statistical support for delimitation of these genetic clusters into two species was provided by a recent Bayesian method. Phylogenetic–coalescent dating with external time calibrations indicates that speciation was relatively recent, with a 95% posterior interval of 0.46–2.66 mya. This postdates equivalent phylogenetic dating estimates of sequence divergence by approximately 1 Ma. Statistical analyses of a small number of independent loci provide important insights into the history of the speciation process in C. mionecton and support delimitation of populations into two species with distributions that are spatially discordant with patterns of morphological variation.     

Improbable but true: the invasive inbreeding ambrosia beetle Xylosandrus morigerus has generalist genotypes

The wide distribution and dominance of invasive inbreeding species in many forest ecosystems seems paradoxical in face of their limited genetic variation. Successful establishment of invasive species in new areas is nevertheless facilitated by clonal reproduction: parthenogenesis, regular self-fertilization, and regular inbreeding. The success of clonal lineages in variable environments has been explained by two models, the frozen niche variation (FNV) model and the general-purpose genotype (GPG) model. We tested these models on a widely distributed forest pest that has been recently established in Costa Rica-the sibling-mating ambrosia beetle Xylosandrus morigerus. Two deeply diverged mitochondrial haplotypes coexist at multiple sites in Costa Rica. We find that these two haplotypes do not differ in their associations with ecological factors. Overall the two haplotypes showed complete overlap in their resource utilization; both genotypes have broad niches, supporting the GPG model. Thus, probable or not, our findings suggest that X. morigerus is a true ecological generalist. Clonal aspects of reproduction coupled with broad niches are doubtless important factors in the successful colonization of new habitats in distant regions.     

Effective/census population size ratio estimation: a compendium and appraisal

With an ecological-evolutionary perspective increasingly applied toward the conservation and management of endangered or exploited species, the genetic estimation of effective population size (N_e) has proliferated. Based on a comprehensive analysis of empirical literature from the past two decades, we asked: (i) how often do studies link N_e to the adult census population size (N)? (ii) To what extent is N_e correctly linked to N? (iii) How readily is uncertainty accounted for in both N_e and N when quantifying N_e/N ratios? and (iv) how frequently and to what degree might errors in the estimation of N_e or N affect inferences of N_e/N ratios? We found that only 20% of available N_e estimates (508 of 2617; 233 studies) explicitly attempted to link N_e and N; of these, only 31% (160 of 508) correctly linked N_e and N. Moreover, only 7% (41 of 508) of N_e/N ratios (correctly linked or not) reported confidence intervals for both N_e and N; for those cases where confidence intervals were reported for N_e only, 31% of N_e/N ratios overlapped with 1, of which more than half also reached below N_e/N = 0.01. Uncertainty in N_e/N ratios thus sometimes spanned at least two orders of magnitude. We conclude that the estimation of N_e/N ratios in natural populations could be significantly improved, discuss several options for doing so, and briefly outline some future research directions.     

Intracoastal shipping drives patterns of regional population expansion by an invasive marine invertebrate

Understanding the factors contributing to expansion of nonnative populations is a critical step toward accurate risk assessment and effective management of biological invasions. Nevertheless, few studies have attempted explicitly to test hypotheses regarding factors driving invasive spread by seeking correlations between patterns of vector movement and patterns of genetic connectivity. Herein, we describe such an attempt for the invasive tunicate Styela clava in the northeastern Pacific. We utilized microsatellite data to estimate gene flow between samples collected throughout the known range of S. clava in the region, and assessed correlation of these estimates with patterns of intracoastal commercial vessel traffic. Our results suggest that recent shipping patterns have contributed to the contemporary distribution of genetic variation. However, the analysis also indicates that other factors—including a complex invasion history and the influence of other vectors—have partially obscured genetic patterns associated with intracoastal population expansion.