Detecting selection-induced departures from Hardy-Weinberg proportions

Chromosomal inversions harbour excess mutational load in the coral, Acropora kenti, on the Great Barrier Reef

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future.

Evolution in response to climate in the native and introduced ranges of a globally distributed plant

The extent to which species can adapt to spatiotemporal climatic variation in their native and introduced ranges remains unresolved. To address this, we examined how clines in cyanogenesis (hydrogen cyanide [HCN] production-an antiherbivore defense associated with decreased tolerance to freezing) have shifted in response to climatic variation in space and time over a 60-year period in both the native and introduced ranges of Trifolium repens. HCN production is a polymorphic trait controlled by variation at two Mendelian loci (Ac and Li). Using phenotypic assays, we estimated within-population frequencies of HCN production and dominant alleles at both loci (i.e., Ac and Li) from 10,575 plants sampled from 131 populations on five continents, and then compared these frequencies to those from historical data collected in the 1950s. There were no clear relationships between changes in the frequency of HCN production, Ac, or Li and changes in temperature between contemporary and historical samples. We did detect evidence of continued evolution to temperature gradients in the introduced range, whereby the slope of contemporary clines for HCN and Ac in relation to winter temperature became steeper than historical clines and more similar to native clines. These results suggest that cyanogenesis clines show no clear changes through time in response to global warming, but introduced populations continue to adapt to their contemporary environments.

Commonly used Hardy-Weinberg equilibrium filtering schemes impact population structure inferences using RADseq data

Reduced representation sequencing (RRS) is a widely used method to assay the diversity of genetic loci across the genome of an organism. The dominant class of RRS approaches assay loci associated with restriction sites within the genome (restriction site associated DNA sequencing, or RADseq). RADseq is frequently applied to non‐model organisms since it enables population genetic studies without relying on well‐characterized reference genomes. However, RADseq requires the use of many bioinformatic filters to ensure the quality of genotyping calls. These filters can have direct impacts on population genetic inference, and therefore require careful consideration. One widely used filtering approach is the removal of loci that do not conform to expectations of Hardy–Weinberg equilibrium (HWE). Despite being widely used, we show that this filtering approach is rarely described in sufficient detail to enable replication. Furthermore, through analyses of in silico and empirical data sets we show that some of the most widely used HWE filtering approaches dramatically impact inference of population structure. In particular, the removal of loci exhibiting departures from HWE after pooling across samples significantly reduces the degree of inferred population structure within a data set (despite this approach being widely used). Based on these results, we provide recommendations for best practice regarding the implementation of HWE filtering for RADseq data sets.

Discrepancies of RET gene and risk of differentiated thyroid carcinoma

BACKGROUND

Somatic variations in rearranged during transfection (RET) proto-oncogene acts to influence Thyroid cancer (TC) in a low penetrance manner, but their effects tend to vary between different populations.

OBJECTIVE

This case-control study was aimed to evaluate effect of RET G691S, S904S and L769L single nucleotide polymorphisms (SNPs) on the risk for differentiated thyroid carcinoma (DTC).

METHODS

A total of 180 patients and 220 controls were genotyped by Polymerase chain reaction - restriction fragment length polymorphism (PCR-RFLP). Di-Deoxy Sanger sequencing was performed on 100 samples with variations and 20 wild samples for each amplified exon. In addition, In Silico tools were used to evaluate structural and functional impact of individual SNPs in disease progression.

RESULTS

In RET G691S/L769L/S904S SNPs, frequency of variant genotypes in DTC cases was 61.1%, 54.4% and 76.6% as compared to 45.9%, 43.6% and 89.09% in controls respectively (P⩽ 0.05). In Silico analysis revealed that different protein formed due to G691S substitution decreases the stability of 3D structure of protein. The RET G691S and L769L SNP followed "Dominant" but RET S904S SNP confirmed an "Additive" mode of inheritance.

CONCLUSION

RET G691S/L769L/S904S SNPs are significantly associated with DTC with G691S SNP declining the stability of final protein product.

Genetic diversity and population structure in divergent German cattle selection lines on the basis of milk protein polymorphisms

The aim of this study was to analyze the genetic structure of the casein cluster in eight selection lines of the Holstein Friesian (HF), German Simmental (GS) and German Black Pied cattle (“Deutsches Schwarzbuntes Niederungsrind”, DSN) breeds. A total of 2962 milk samples were typed at αs1-casein (αs1-CN), β-casein (β-CN), αs2-casein (αs2-CN) and κ-casein (κ-CN) loci using isoelectric focusing. The number of alleles per locus ranged from one (αs2-CN) to five (β-CN), and the average expected heterozygosity and polymorphic information content of all loci were 0.33 and 0.27, respectively. The unrooted dendrogram revealed that the selection lines of the endangered DSN breed were clearly separated from the HF and GS breeds due to their predominance of the β-CN A1 allele and the comprehensive haplotype BA1A (in the abbreviation of αs1-β-κ-CN). Temporal changes in allele distributions indicated decreasing genetic diversity at the casein loci, explaining the moderate level of genetic differentiation among selection lines (7.1 %). The variability of the casein should be exploited in future using breeding programs to select genetic lines for specific protein production in bovine milk but also to preserve biodiversity.

The search for sexually antagonistic genes: Practical insights from studies of local adaptation and statistical genomics

Sexually antagonistic (SA) genetic variation-in which alleles favored in one sex are disfavored in the other-is predicted to be common and has been documented in several animal and plant populations, yet we currently know little about its pervasiveness among species or its population genetic basis. Recent applications of genomics in studies of SA genetic variation have highlighted considerable methodological challenges to the identification and characterization of SA genes, raising questions about the feasibility of genomic approaches for inferring SA selection. The related fields of local adaptation and statistical genomics have previously dealt with similar challenges, and lessons from these disciplines can therefore help overcome current difficulties in applying genomics to study SA genetic variation. Here, we integrate theoretical and analytical concepts from local adaptation and statistical genomics research-including F ST and F IS statistics, genome-wide association studies, pedigree analyses, reciprocal transplant studies, and evolve-and-resequence experiments-to evaluate methods for identifying SA genes and genome-wide signals of SA genetic variation. We begin by developing theoretical models for between-sex F ST and F IS, including explicit null distributions for each statistic, and using them to critically evaluate putative multilocus signals of sex-specific selection in previously published datasets. We then highlight new statistics that address some of the limitations of F ST and F IS, along with applications of more direct approaches for characterizing SA genetic variation, which incorporate explicit fitness measurements. We finish by presenting practical guidelines for the validation and evolutionary analysis of candidate SA genes and discussing promising empirical systems for future work.

HLA-A, -B, -DRB1 and -DQB1 polymorphisms among Iraqi Arabs

In this report, HLA polymorphisms (A, B, DRB1 and DQB1 loci) were determined in 149 unrelated Iraqi Arab potential bone marrow and kidney donors. Molecular genotyping was carried out by polymerase chain reaction followed by specific oligonucleotide probe hybridizations. Data were analyzed by Arlequin software. HLA-A, -B and -DRB1 genotype frequencies were significantly deviated from Hardy-Weinberg equilibrium, while HLA-DQB1 frequencies showed no deviation. A*03, B*35, DRB1*11 and DQB1*02 were the most frequent allele groups, while A*02-B*07-DRB1*04-DQB1*03 was the most frequent haplotype. HLA data are available in the Allele Frequencies Net Database (AFND: 3680) under the population name "Iraq Arabs".

Predicting the strength of urban-rural clines in a Mendelian polymorphism along a latitudinal gradient

Cities are emerging as models for addressing the fundamental question of whether populations evolve in parallel to similar environments. Here, we examine the environmental factors that drive the evolution of parallel urban‐rural clines in a Mendelian trait—the cyanogenic antiherbivore defense of white clover (Trifolium repens). Previous work suggested urban‐rural gradients in frost and snow depth could drive the evolution of reduced hydrogen cyanide (HCN) frequencies in urban populations. Here, we sampled over 700 urban and rural clover populations across 16 cities along a latitudinal transect in eastern North America. In each population, we quantified changes in the frequency of genotypes that produce HCN, and in a subset of the cities we estimated the frequency of the alleles at the two genes (CYP79D15 and Li) that epistatically interact to produce HCN. We then tested the hypothesis that cold climatic conditions are necessary for the evolution of cyanogenesis clines by comparing the strength of clines among cities located along a latitudinal gradient of winter temperature and frost exposure. Overall, half of the cities exhibited urban‐rural clines in the frequency of HCN, whereby urban populations evolved lower HCN frequencies. Clines did not evolve in cities with the lowest temperatures and greatest snowfall, supporting the hypothesis that snow buffers plants against winter frost and constrains the formation of clines. By contrast, the strongest clines occurred in the warmest cities where snow and frost are rare, suggesting that alternative selective agents are maintaining clines in warmer cities. Some clines were driven by evolution at only CYP79D15, consistent with stronger and more consistent selection on this locus than on Li. Together, our results demonstrate that urban environments often select for similar phenotypes, but different selective agents and targets underlie the evolutionary response in different cities.

Sex linkage of the skeletal muscle sodium channel gene (SCN4A) explains apparent deviations from Hardy-Weinberg equilibrium of tetrodotoxin-resistance alleles in garter snakes (Thamnophis sirtalis)

The arms race between tetrodotoxin-bearing Pacific newts (Taricha) and their garter snake predators (Thamnophis) in western North America has become a classic example of coevolution, shedding light on predator-prey dynamics, the molecular basis of adaptation, and patterns of convergent evolution. Newts are defended by tetrodotoxin (TTX), a neurotoxin that binds to voltage-gated sodium channels (Na_v proteins), arresting electrical activity in nerves and muscles and paralyzing would-be predators. However, populations of the common garter snake (T. sirtalis) have overcome this defense, largely through polymorphism at the locus SCN4A, which renders the encoded protein (Na_v1.4) less vulnerable to TTX. Previous work suggests that SCN4A commonly shows extreme deviations from Hardy-Weinberg equilibrium (HWE) in these populations, which has been interpreted as the result of intense selection imposed by newts. Here we show that much of this apparent deviation can be attributed to sex linkage of SCN4A. Using genomic data and quantitative PCR, we show that SCN4A is on the Z chromosome in Thamnophis and other advanced snakes. Taking Z-linkage into account, we find that most apparent deviations from HWE can be explained by female hemizygosity rather than low heterozygosity. Sex linkage can affect mutation rates, selection, and drift, and our results suggest that Z-linkage of SCN4A may make significant contributions to the overall dynamics of the coevolutionary arms race between newts and snakes.

microRNA-27a and microRNA-146a SNP in cerebral malaria

Background

During Plasmodium falciparum infection, microRNA expression alters in brain tissue of mice with cerebral malaria compared to noninfected controls. MicroRNA regulates gene expression post‐transcriptionally to influence biological processes. Cerebral malaria pathology caused mainly by the immunological disorder. We hypothesize that single‐nucleotide polymorphism in a microRNA influences microRNA biogenesis or target gene recognition and altering susceptibility to cerebral malaria.

Methods

We performed a literature search based on immunological mechanism and applied microRNA‐related single‐nucleotide polymorphisms database to examine candidate microRNA SNPs possibly responsible for cerebral malaria. MicroRNA‐27a and microRNA‐146a are supposed to involve in cerebral malaria pathology. To assess the relationship of microRNA SNP to cerebral malaria outcome, we performed TaqMan Genotyping Assays in 110 cerebral malaria and 207 uncomplicated malaria cases for three candidate microRNA SNPs (rs895819 of microRNA‐27a, rs57095329 and rs2910164 of microRNA‐146a).

Results

Our study detected no significant difference in genotype and allele frequency of individual microRNA SNPs as well as in haplotypes of microRNA‐146a between these two groups of malaria patients in Thailand. Hardy–Weinberg disequilibrium of rs57095329 in the cerebral malaria group showed a heterozygous excess which might be due to natural selection.

Conclusion

Our data supported that the candidate microRNA SNPs have no major role to develop cerebral malaria.

Signatures of selection in mammalian clock genes with coding trinucleotide repeats: Implications for studying the genomics of high-pace adaptation

Climate change is predicted to affect the reproductive ecology of wildlife; however, we have yet to understand if and how species can adapt to the rapid pace of change. Clock genes are functional genes likely critical for adaptation to shifting seasonal conditions through shifts in timing cues. Many of these genes contain coding trinucleotide repeats, which offer the potential for higher rates of change than single nucleotide polymorphisms (SNPs) at coding sites, and, thus, may translate to faster rates of adaptation in changing environments. We characterized repeats in 22 clock genes across all annotated mammal species and evaluated the potential for selection on repeat motifs in three clock genes (NR1D1,CLOCK, and PER1) in three congeneric species pairs with different latitudinal range limits: Canada lynx and bobcat (Lynx canadensis and L. rufus), northern and southern flying squirrels (Glaucomys sabrinus and G. volans), and white‐footed and deer mouse (Peromyscus leucopus and P. maniculatus). Signatures of positive selection were found in both the interspecific comparison of Canada lynx and bobcat, and intraspecific analyses in Canada lynx. Northern and southern flying squirrels showed differing frequencies at common CLOCK alleles and a signature of balancing selection. Regional excess homozygosity was found in the deer mouse at PER1 suggesting disruptive selection, and further analyses suggested balancing selection in the white‐footed mouse. These preliminary signatures of selection and the presence of trinucleotide repeats within many clock genes warrant further consideration of the importance of candidate gene motifs for adaptation to climate change.

Testing for Hardy-Weinberg proportions: have we lost the plot?

Testing for Hardy-Weinberg proportions (HWP) is routine in almost all genetic studies of natural populations, but many researchers do not demonstrate a full understanding of the purposes of these tests or how to interpret the results. Common problems include a lack of understanding of statistical power and the difference between statistical significance and biological significance, how to interpret results of multiple tests, and how to distinguish between various factors that can cause statistically significant departures. In this perspective, which focuses on analysis of genetic data for nonmodel species, I 1) review factors that can cause departures from HWP at individual loci and linkage disequilibrium (LD) at pairs of loci; 2) discuss commonly used tests for HWP and LD, with an emphasis on multiple-testing issues; 3) show how to distinguish among possible causes of departures from HWP; and 4) outline some simple steps to follow when significant test results are found. Finally, I 5) identify some issues that merit particular attention as we move into an era in which analysis of genomics-scale datasets for nonmodel species is commonplace.

Toy models for macroevolutionary patterns and trends

Many models have been used to simplify and operationalize the subtle but complex mechanisms of biological evolution. Toy models are gross simplifications that nevertheless attempt to retain major essential features of evolution, bridging the gap between empirical reality and formal theoretical understanding. In this paper, we examine thirteen models which describe evolution that also qualify as such toy models, including the tree of life, branching processes, adaptive ratchets, fitness landscapes, and the role of nonlinear avalanches in evolutionary dynamics. Such toy models are intended to capture features such as evolutionary trends, coupled evolutionary dynamics of phenotype and genotype, adaptive change, branching, and evolutionary transience. The models discussed herein are applied to specific evolutionary contexts in various ways that simplify the complexity inherent in evolving populations. While toy models are overly simplistic, they also provide sufficient dynamics for capturing the fundamental mechanism(s) of evolution. Toy models might also be used to aid in high-throughput data analysis and the understanding of cultural evolutionary trends. This paper should serve as an introductory guide to the toy modeling of evolutionary complexity.

Searching for the bull's eye: agents and targets of selection vary among geographically disparate cyanogenesis clines in white clover (Trifolium repens L.)

The recurrent evolution of adaptive clines within a species can be used to elucidate the selective factors and genetic responses that underlie adaptation. White clover is polymorphic for cyanogenesis (HCN release with tissue damage), and climate-associated cyanogenesis clines have evolved throughout the native and introduced species range. This polymorphism arises through two independently segregating Mendelian polymorphisms for the presence/absence of two required components: cyanogenic glucosides and their hydrolyzing enzyme linamarase. Cyanogenesis is commonly thought to function in herbivore defense; however, the individual cyanogenic components may also serve other physiological functions. To test whether cyanogenesis clines have evolved in response to the same selective pressures acting on the same genetic targets, we examined cyanogenesis cline shape and its environmental correlates in three world regions: southern New Zealand, the central United States and the US Pacific Northwest. For some regional comparisons, cline shapes are remarkably similar despite large differences in the spatial scales over which clines occur (40-1600 km). However, we also find evidence for major differences in both the agents and targets of selection among the sampled clines. Variation in cyanogenesis frequency is best predicted using a combination of minimum winter temperature and aridity variables. Together, our results provide evidence that recurrent adaptive clines do not necessarily reflect shared adaptive processes.

No association of catechol-O-methyltransferase polymorphisms with schizophrenia in the Han Chinese population

AIMS

Genetics play a major role in the etiology of schizophrenia (SZ). Catechol-O-methyltransferase (COMT) is one of the promising candidate genes for SZ. A nonsynonymous single-nucleotide polymorphism (SNP), rs4680, causing a Valine (Val) to Methionine (Met) substitution, has been widely studied in relation to psychiatric phenotypes, including SZ, but with conflicting results. We conducted a two-stage study to examine the association of COMT polymorphisms with SZ in the Han Chinese population.

RESULTS

Association analysis of nine SNPs in 768 patients and 1348 controls failed to detect any positive markers or haplotypes. Then, we tested rs4680 in a validation sample of 963 patients and 992 controls, and no significant association was observed, but the cases significantly deviated from Hardy-Weinberg equilibrium (p=5.7e-4). There was no association of rs4680 with SZ in the combined sample (n=4071, p=0.110, odds ratio=1.08).

CONCLUSIONS

Our results do not support the association of COMT with SZ in the Han Chinese population.

Association between immunoglobulin GM and KM genotypes and placental malaria in HIV-1 negative and positive women in western Kenya

Immunoglobulin (Ig) GM and KM allotypes, genetic markers of γ and κ chains, are associated with humoral immune responsiveness. Previous studies have shown the relationships between GM6-carrying haplotypes and susceptibility to malaria infection in children and adults; however, the role of the genetic markers in placental malaria (PM) infection and PM with HIV co-infection during pregnancy has not been investigated. We examined the relationship between the gene polymorphisms of Ig GM6 and KM allotypes and the risk of PM infection in pregnant women with known HIV status. DNA samples from 728 pregnant women were genotyped for GM6 and KM alleles using polymerase chain reaction-restriction fragment length polymorphism method. Individual GM6 and KM genotypes and the combined GM6 and KM genotypes were assessed in relation to PM in HIV-1 negative and positive women, respectively. There was no significant effect of individual GM6 and KM genotypes on the risk of PM infection in HIV-1 negative and positive women. However, the combination of homozygosity for GM6(+) and KM3 was associated with decreased risk of PM (adjusted OR, 0.25; 95% CI, 0.08-0.8; P = 0.019) in HIV-1 negative women while in HIV-1 positive women the combination of GM6(+/-) with either KM1-3 or KM1 was associated with increased risk of PM infection (adjusted OR, 2.10; 95% CI, 1.18-3.73; P = 0.011). Hardy-Weinberg Equilibrium (HWE) tests further showed an overall significant positive F(is) (indication of deficit in heterozygotes) for GM6 while there was no deviation for KM genotype frequency from HWE in the same population. These findings suggest that the combination of homozygous GM6(+) and KM3 may protect against PM in HIV-1 negative women while the HIV-1 positive women with heterozygous GM6(+/-) combined with KM1-3 or KM1 may be more susceptible to PM infection. The deficit in heterozygotes for GM6 further suggests that GM6 could be under selection likely by malaria infection.

Large-scale MHC class II genotyping of a wild lemur population by next generation sequencing

The critical role of major histocompatibility complex (MHC) genes in disease resistance, along with their putative function in sexual selection, reproduction and chemical ecology, make them an important genetic system in evolutionary ecology. Studying selective pressures acting on MHC genes in the wild nevertheless requires population-wide genotyping, which has long been challenging because of their extensive polymorphism. Here, we report on large-scale genotyping of the MHC class II loci of the grey mouse lemur (Microcebus murinus) from a wild population in western Madagascar. The second exons from MHC-DRB and -DQB of 772 and 672 individuals were sequenced, respectively, using a 454 sequencing platform, generating more than 800,000 reads. Sequence analysis, through a stepwise variant validation procedure, allowed reliable typing of more than 600 individuals. The quality of our genotyping was evaluated through three independent methods, namely genotyping the same individuals by both cloning and 454 sequencing, running duplicates, and comparing parent-offspring dyads; each displaying very high accuracy. A total of 61 (including 20 new) and 60 (including 53 new) alleles were detected at DRB and DQB genes, respectively. Both loci were non-duplicated, in tight linkage disequilibrium and in Hardy-Weinberg equilibrium, despite the fact that sequence analysis revealed clear evidence of historical selection. Our results highlight the potential of 454 sequencing technology in attempts to investigate patterns of selection shaping MHC variation in contemporary populations. The power of this approach will nevertheless be conditional upon strict quality control of the genotyping data.

Evolutionary history and adaptation from high-coverage whole-genome sequences of diverse African hunter-gatherers

To reconstruct modern human evolutionary history and identify loci that have shaped hunter-gatherer adaptation, we sequenced the whole genomes of five individuals in each of three different hunter-gatherer populations at > 60× coverage: Pygmies from Cameroon and Khoesan-speaking Hadza and Sandawe from Tanzania. We identify 13.4 million variants, substantially increasing the set of known human variation. We found evidence of archaic introgression in all three populations, and the distribution of time to most recent common ancestors from these regions is similar to that observed for introgressed regions in Europeans. Additionally, we identify numerous loci that harbor signatures of local adaptation, including genes involved in immunity, metabolism, olfactory and taste perception, reproduction, and wound healing. Within the Pygmy population, we identify multiple highly differentiated loci that play a role in growth and anterior pituitary function and are associated with height.