Comparative studies of quantitative trait and neutral marker divergence: a meta-analysis

Comparative studies of quantitative genetic and neutral marker differentiation have provided means for assessing the relative roles of natural selection and random genetic drift in explaining among-population divergence. This information can be useful for our fundamental understanding of population differentiation, as well as for identifying management units in conservation biology. Here, we provide comprehensive review and meta-analysis of the empirical studies that have compared quantitative genetic (Q(ST)) and neutral marker (F(ST)) differentiation among natural populations. Our analyses confirm the conclusion from previous reviews - based on ca. 100% more data - that the Q(ST) values are on average higher than F(ST) values [mean difference 0.12 (SD 0.27)] suggesting a predominant role for natural selection as a cause of differentiation in quantitative traits. However, although the influence of trait (life history, morphological and behavioural) and marker type (e.g. microsatellites and allozymes) on the variance of the difference between Q(ST) and F(ST) is small, there is much heterogeneity in the data attributable to variation between specific studies and traits. The latter is understandable as there is no reason to expect that natural selection would be acting in similar fashion on all populations and traits (except for fitness itself). We also found evidence to suggest that Q(ST) and F(ST) values across studies are positively correlated, but the significance of this finding remains unclear. We discuss these results in the context of utility of the Q(ST)-F(ST) comparisons as a tool for inferring natural selection, as well as associated methodological and interpretational problems involved with individual and meta-analytic studies.     

multi‐dice: r package for comparative population genomic inference under hierarchical co‐demographic models of independent single‐population size changes

Population genetic data from multiple taxa can address comparative phylogeographic questions about community‐scale response to environmental shifts, and a useful strategy to this end is to employ hierarchical co‐demographic models that directly test multi‐taxa hypotheses within a single, unified analysis. This approach has been applied to classical phylogeographic data sets such as mitochondrial barcodes as well as reduced‐genome polymorphism data sets that can yield 10,000s of SNPs, produced by emergent technologies such as RAD‐seq and GBS. A strategy for the latter had been accomplished by adapting the site frequency spectrum to a novel summarization of population genomic data across multiple taxa called the aggregate site frequency spectrum (aSFS), which potentially can be deployed under various inferential frameworks including approximate Bayesian computation, random forest and composite likelihood optimization. Here, we introduce the r package multi‐dice , a wrapper program that exploits existing simulation software for flexible execution of hierarchical model‐based inference using the aSFS, which is derived from reduced genome data, as well as mitochondrial data. We validate several novel software features such as applying alternative inferential frameworks, enforcing a minimal threshold of time surrounding co‐demographic pulses and specifying flexible hyperprior distributions. In sum, multi‐dice provides comparative analysis within the familiar R environment while allowing a high degree of user customization, and will thus serve as a tool for comparative phylogeography and population genomics.     

A framework for power and sensitivity analyses for quantitative genetic studies of natural populations, and case studies in Soay sheep (Ovis aries)

Studies of the quantitative genetics of natural populations have contributed greatly to evolutionary biology in recent years. However, while pedigree data required are often uncertain (i.e. incomplete and partly erroneous) and limited, means to evaluate the effects of such uncertainties have not been developed. We have therefore developed a general framework for power and sensitivity analyses of such studies. We propose that researchers first generate a set of pedigree data that they wish to use in a quantitative genetic study, as well as data regarding errors that occur in that pedigree. This pedigree is then permuted using the data regarding errors to generate hypothetical 'true' and 'assumed' pedigrees that differ so as to mimic pedigree errors that might occur in the study system under consideration. Phenotypic data are then simulated across the true pedigree (according to user-defined genetic and environmental covariance structures), before being analysed with standard quantitative genetic techniques in conjunction with the 'assumed' pedigree data. To illustrate this approach, we conducted power and sensitivity analyses in a well-known study of Soay sheep (Ovis aries). We found that, although the estimation of simple genetic (co)variance structures is fairly robust to pedigree errors, some potentially serious biases were detected under more complex scenarios involving maternal effects. Power analyses also showed that this study system provides high power to detect heritabilities as low as about 0.09. Given this range of results, we suggest that such power and sensitivity analyses could greatly complement empirical studies, and we provide the computer program PEDANTICS to aid in their application.     

Rapid maximum likelihood ancestral state reconstruction of continuous characters: A rerooting‐free algorithm

Ancestral state reconstruction is a method used to study the evolutionary trajectories of quantitative characters on phylogenies. Although efficient methods for univariate ancestral state reconstruction under a Brownian motion model have been described for at least 25 years, to date no generalization has been described to allow more complex evolutionary models, such as multivariate trait evolution, non‐Brownian models, missing data, and within‐species variation. Furthermore, even for simple univariate Brownian motion models, most phylogenetic comparative R packages compute ancestral states via inefficient tree rerooting and full tree traversals at each tree node, making ancestral state reconstruction extremely time‐consuming for large phylogenies. Here, a computationally efficient method for fast maximum likelihood ancestral state reconstruction of continuous characters is described. The algorithm has linear complexity relative to the number of species and outperforms the fastest existing R implementations by several orders of magnitude. The described algorithm is capable of performing ancestral state reconstruction on a 1,000,000‐species phylogeny in fewer than 2 s using a standard laptop, whereas the next fastest R implementation would take several days to complete. The method is generalizable to more complex evolutionary models, such as phylogenetic regression, within‐species variation, non‐Brownian evolutionary models, and multivariate trait evolution. Because this method enables fast repeated computations on phylogenies of virtually any size, implementation of the described algorithm can drastically alleviate the computational burden of many otherwise prohibitively time‐consuming tasks requiring reconstruction of ancestral states, such as phylogenetic imputation of missing data, bootstrapping procedures, Expectation‐Maximization algorithms, and Bayesian estimation. The described ancestral state reconstruction algorithm is implemented in the Rphylopars functions anc.recon and phylopars.     

RAD‐sequencing for estimating genomic relatedness matrix‐based heritability in the wild: A case study in roe deer

Estimating the evolutionary potential of quantitative traits and reliably predicting responses to selection in wild populations are important challenges in evolutionary biology. The genomic revolution has opened up opportunities for measuring relatedness among individuals with precision, enabling pedigree‐free estimation of trait heritabilities in wild populations. However, until now, most quantitative genetic studies based on a genomic relatedness matrix (GRM) have focused on long‐term monitored populations for which traditional pedigrees were also available, and have often had access to knowledge of genome sequence and variability. Here, we investigated the potential of RAD‐sequencing for estimating heritability in a free‐ranging roe deer (Capreolous capreolus) population for which no prior genomic resources were available. We propose a step‐by‐step analytical framework to optimize the quality and quantity of the genomic data and explore the impact of the single nucleotide polymorphism (SNP) calling and filtering processes on the GRM structure and GRM‐based heritability estimates. As expected, our results show that sequence coverage strongly affects the number of recovered loci, the genotyping error rate and the amount of missing data. Ultimately, this had little effect on heritability estimates and their standard errors, provided that the GRM was built from a minimum number of loci (above 7,000). Genomic relatedness matrix‐based heritability estimates thus appear robust to a moderate level of genotyping errors in the SNP data set. We also showed that quality filters, such as the removal of low‐frequency variants, affect the relatedness structure of the GRM, generating lower h² estimates. Our work illustrates the huge potential of RAD‐sequencing for estimating GRM‐based heritability in virtually any natural population.     

A quantitative review of MHC‐based mating preference: the role of diversity and dissimilarity

Sexual selection hypotheses stipulate that the major histocompatibility complex genes (MHC) constitute a key molecular underpinning for mate choice in vertebrates. The last four decades saw growing empirical literature on the role of MHC diversity and dissimilarity in mate choice for a wide range of vertebrate animals, but with mixed support for its significance in natural populations. Using formal phylogenetic meta‐analysis and meta‐regression techniques, we quantitatively review the existing literature on MHC‐dependent mating preferences in nonhuman vertebrates with a focus on the role of MHC diversity and dissimilarity. Overall, we found small, statistically nonsignificant, average effect sizes for both diversity‐ and dissimilarity‐based mate choice (r = 0.113 and 0.064, respectively). Importantly, however, meta‐regression models revealed statistically significant support regarding female choice for diversity, and choice for dissimilarity (regardless of choosy sex) only when dissimilarity is characterized across multiple loci. Little difference was found among vertebrate taxa; however, the lack of statistical power meant statistically significant effects were limited to some taxa. We found little sign of publication bias; thus, our results are likely to be robust. In light of our quantitative assessment, methodological improvements and fruitful future avenues of research are highlighted.     

BETYdb: a yield,trait, and ecosystem service database applied to second‐generation bioenergy feedstock production

Increasing demand for sustainable energy has led to research and development on the cultivation of diverse plant species for biomass production. To support the research and development required to domesticate and cultivate crops for bioenergy, we developed the Biofuel Ecophysiological Traits and Yields database (BETYdb). BETYdb is a centralized open‐access repository that facilitates organization, discovery, and exchange of information about plant traits, crop yields, and ecosystem functions. BETYdb provides user interfaces to simplify storage and discovery as well as programming interfaces that support automated and reproducible scientific workflows. Presently, BETYdb contains over forty thousand observations of plant traits, biomass yields, and ecosystem dynamics collected from the published articles and ongoing field studies. Over half of these records represent fewer than ten genera that have been intensively evaluated for biomass production, while the other half represent over two thousand plant species reflecting research on new crops, unmanaged ecosystems, and land use transitions associated with bioenergy. BETYdb has been accessed over twenty‐five thousand times and is used in the fields of bioenergy and ecosystem ecology to quantify yield potential and ecosystem functioning of crops and unmanaged systems under present and future climates. Here, we summarize the database contents and illustrate its applications. We show its utility in a new analysis that confirms that Miscanthus is twice as productive as switchgrass over a much wider range of environmental and management conditions than covered in previous analyses. We compare traits related to carbon uptake and water use of these species with each other and with two coppice shrubs, poplar and willow. These examples, along with a growing body of published research that used BETYdb, illustrate the scope of research supported through this open‐access database.     

Three novel quantitative trait loci for skin thickness in swine identified by linkage and genome‐wide association studies

Skin is the largest organ in the pig body and plays a key role in protecting the body against pathogens and excessive water loss. Deciphering the genetic basis of swine skin thickness would enrich our knowledge about the skin. To identify the loci for porcine skin thickness, we first performed a genome scan with 194 microsatellite markers in a White Duroc × Erhualian F₂ intercross. We identified three genome‐wide significant QTL on pig chromosomes (SSC) 4, 7 and 15 using linkage analysis. The most significant QTL was found on SSC7 with a small confidence interval of ~5 cM, explaining 23.9 percent of phenotypic variance. Further, we conducted a genome‐wide association study (GWAS) using Illumina PorcineSNP60 Beadchips for the F₂ pedigree and a population of Chinese Sutai pigs. We confirmed significant QTL in the F₂ pedigree and replicated QTL on SSC15 in Chinese Sutai pigs. A meta‐analysis of GWASs on both populations detected a genomic region associated with skin thickness on SSC4. GWAS results were generally consistent with QTL mapping. Identical‐by‐descent analysis defined QTL on SSC7 in a 683‐kb region harboring an interesting candidate gene: HMGA1. On SSC15, the linkage disequilibrium analysis showed a haplotype block of 2.20 Mb that likely harbors the gene responsible for skin thickness. Our findings provide novel insights into the genetic basis of swine skin thickness, which would benefit further understanding of porcine skin function.     

Psychotherapies for depression in low‐ and middle‐income countries: a meta‐analysis

Most psychotherapies for depression have been developed in high‐income Western countries of North America, Europe and Australia. A growing number of randomized trials have examined the effects of these treatments in non‐Western countries. We conducted a meta‐analysis of these studies to examine whether these psychotherapies are effective and to compare their effects between studies from Western and non‐Western countries. We conducted systematic searches in bibliographical databases and included 253 randomized controlled trials, of which 32 were conducted in non‐Western countries. The effects of psychotherapies in non‐Western countries were large (g=1.10; 95% CI: 0.91‐1.30), with high heterogeneity (I²=90; 95% CI: 87‐92). After adjustment for publication bias, the effect size dropped to g=0.73 (95% CI: 0.51‐0.96). Subgroup analyses did not indicate that adaptation to the local situation was associated with the effect size. Comparisons with the studies in Western countries showed that the effects of the therapies were significantly larger in non‐Western countries, also after adjusting for characteristics of the participants, the treatments and the studies. These larger effect sizes in non‐Western countries may reflect true differences indicating that therapies are indeed more effective; or may be explained by the care‐as‐usual control conditions in non‐Western countries, often indicating that no care was available; or may be the result of the relative low quality of many trials in the field. This study suggests that psychotherapies that were developed in Western countries may or may not be more effective in non‐Western countries, but they are probably no less effective and can therefore also be used in these latter countries.     

Evolutionary conservation,diversity and specificity of LTR‐retrotransposons in flowering plants: insights from genome‐wide analysis and multi‐specific comparison

The availability of complete or nearly complete genome sequences from several plant species permits detailed discovery and cross‐species comparison of transposable elements (TEs) at the whole genome level. We initially investigated 510 long terminal repeat‐retrotransposon (LTR‐RT) families comprising 32 370 elements in soybean (Glycine max (L.) Merr.). Approximately 87% of these elements were located in recombination‐suppressed pericentromeric regions, where the ratio (1.26) of solo LTRs to intact elements (S/I) is significantly lower than that of chromosome arms (1.62). Further analysis revealed a significant positive correlation between S/I and LTR sizes, indicating that larger LTRs facilitate solo LTR formation. Phylogenetic analysis revealed seven Copia and five Gypsy evolutionary lineages that were present before the divergence of eudicot and monocot species, but the scales and timeframes within which they proliferated vary dramatically across families, lineages and species, and notably, a Copia lineage has been lost in soybean. Analysis of the physical association of LTR‐RTs with centromere satellite repeats identified two putative centromere retrotransposon (CR) families of soybean, which were grouped into the CR (e.g. CRR and CRM) lineage found in grasses, indicating that the ‘functional specification’ of CR pre‐dates the bifurcation of eudicots and monocots. However, a number of families of the CR lineage are not concentrated in centromeres, suggesting that their CR roles may now be defunct. Our data also suggest that the envelope‐like genes in the putative Copia retrovirus‐like family are probably derived from the Gypsy retrovirus‐like lineage, and thus we propose the hypothesis of a single ancient origin of envelope‐like genes in flowering plants.     

CTLA4 and generalized vitiligo: two genetic association studies and a meta‐analysis of published data

Several lines of evidence have implicated the gene encoding cytotoxic T lymphocyte antigen 4 (CTLA4) in susceptibility to various autoimmune diseases. However, published studies of genetic association between CTLA4 polymorphisms and vitiligo have yielded conflicting results. Here, we describe two new genetic association studies of CTLA4 single‐nucleotide polymorphisms (SNPs) and generalized vitiligo in two independent Romanian Caucasian (CEU) case‐control cohorts. The first study, of SNPs rs1863800, rs231806, rs231775, rs3087243, rs11571302, rs11571297, and rs10932037, showed no allelic, genotypic, or haplotypic association with generalized vitiligo. The second study, of SNP rs231775, likewise showed no significant association. To enhance statistical power over that of any individual study, we carried out a meta‐analysis that incorporated these two new studies and all other published genetic association studies of CTLA4 SNPs and vitiligo in CEU populations. While there was no association with vitiligo overall, the meta‐analysis showed significant association of SNP rs231775 in that subgroup of vitiligo patients who also had other concomitant autoimmune diseases. Similarly, there was near‐significant association in this same patient subgroup with several other CTLA4 SNPs that are in linkage disequilibrium with rs231775. Our results indicate that the association of CTLA4 with vitiligo is weak, and indeed may be secondary, driven by primary genetic association of CTLA4 with other autoimmune diseases that are epidemiologically associated with vitiligo.     

The dilemma of Fisherian sexual selection: Mate choice for indirect benefits despite rarity and overall weakness of trait‐preference genetic correlation

Fisher's mechanism of sexual selection is a fundamental element of evolutionary theory. In it nonrandom mate choice causes a genetic covariance between a male trait and female preference for that trait and thereby generates a positive feedback process sustaining accelerated coevolution of the trait and preference. Numerous theoretical models of Fisher's mechanism have confirmed its mathematical underpinnings, yet biologists have often failed to find evidence for trait‐preference genetic correlation in populations in which the mechanism was expected to function. We undertook a survey of the literature to conduct a formal meta‐analysis probing the incidence and strength of trait‐preference correlation among animal species. Our meta‐analysis found significant positive genetic correlations in fewer than 20% of the species studied and an overall weighted correlation that is slightly positive. Importantly, a significant positive correlation was not found in any thorough study that included multiple subgroups. We discuss several ways in which the dynamic, multivariate nature of mate choice may reduce the trait‐preference genetic correlation predicted by Fisher's mechanism. We then entertain the possibilities that Fisherian‐like processes sometimes function without genetic correlation, and that mate choice may persist in a population as long as genetic correlation, and therefore Fisher's mechanism, occurs intermittently.     

Age and size at maturity: A quantitative review of diet‐induced reaction norms in insects

Optimality models predict that diet‐induced bivariate reaction norms for age and size at maturity can have diverse shapes, with the slope varying from negative to positive. To evaluate these predictions, we perform a quantitative review of relevant data, using a literature‐derived database of body sizes and development times for over 200 insect species. We show that bivariate reaction norms with a negative slope prevail in nearly all taxonomic and ecological categories of insects as well as in some other ectotherm taxa with comparable life histories (arachnids and amphibians). In insects, positive slopes are largely limited to species, which feed on discrete resource items, parasitoids in particular. By contrast, with virtually no meaningful exceptions, herbivorous and predatory insects display reaction norms with a negative slope. This is consistent with the idea that predictable resource depletion, a scenario selecting for positively sloped reaction norms, is not frequent for these insects. Another source of such selection—a positive correlation between resource levels and juvenile mortality rates—should similarly be rare among insects. Positive slopes can also be predicted by models which integrate life‐history evolution and population dynamics. As bottom‐up regulation is not common in most insect groups, such models may not be most appropriate for insects.     

Tardive dyskinesia risk with first‐ and second‐generation antipsychotics in comparative randomized controlled trials: a meta‐analysis

Tardive dyskinesia (TD) risk with D2/serotonin receptor antagonists or D2 receptor partial agonists (second‐generation antipsychotics, SGAs) is considered significantly lower than with D2 antagonists (first‐generation antipsychotics, FGAs). As some reports questioned this notion, we meta‐analyzed randomized controlled studies (RCTs) to estimate the risk ratio (RR) and annualized rate ratio (RaR) of TD comparing SGAs vs. FGAs and SGAs vs. SGAs. Additionally, we calculated raw and annualized pooled TD rates for each antipsychotic. Data from 57 head‐to‐head RCTs, including 32 FGA and 86 SGA arms, were meta‐analyzed, yielding 32 FGA‐SGA pairs and 35 SGA‐SGA pairs. The annualized TD incidence across FGA arms was 6.5% (95% CI: 5.3‐7.8%) vs. 2.6% (95% CI: 2.0‐3.1%) across SGA arms. TD risk and annualized rates were lower with SGAs vs. FGAs (RR=0.47, 95% CI: 0.39‐0.57, p<0.0001, k=28; RaR=0.35, 95% CI: 0.28‐0.45, p<0.0001, number‐needed‐to‐treat, NNT=20). Meta‐regression showed no FGA dose effect on FGA‐SGA comparisons (Z=?1.03, p=0.30). FGA‐SGA TD RaRs differed by SGA comparator (Q=21.8, df=7, p=0.003), with a significant advantage of olanzapine and aripiprazole over other non‐clozapine SGAs in exploratory pairwise comparisons. SGA‐SGA comparisons confirmed the olanzapine advantage vs. non‐clozapine SGAs (RaR=0.66, 95% CI: 0.49‐0.88, p=0.006, k=17, NNT=100). This meta‐analysis confirms a clinically meaningfully lower TD risk with SGAs vs. FGAs, which is not driven by high dose FGA comparators, and documents significant differences with respect to this risk between individual SGAs.