Genome-Wide Detection of Selective Signatures in Chicken through High Density SNPs

Chicken is recognized as an excellent model for studies of genetic mechanism of phenotypic and genomic evolution, with large effective population size and strong human-driven selection. In the present study, we performed Extended Haplotype Homozygosity (EHH) tests to identify significant core regions employing 600K SNP Chicken chip in an F2 population of 1,534 hens, which was derived from reciprocal crosses between White Leghorn and Dongxiang chicken. Results indicated that a total of 49,151 core regions with an average length of 9.79 Kb were identified, which occupied approximately 52.15% of genome across all autosomes, and 806 significant core regions attracted us mostly. Genes in candidate regions may experience positive selection and were considered to have possible influence on beneficial economic traits. A panel of genes including AASDHPPT, GDPD5, PAR3, SOX6, GPC1 and a signal pathway of AKT1 were detected with the most extreme P-values. Further enrichment analyses indicated that these genes were associated with immune function, sensory organ development and neurogenesis, and may have experienced positive selection in chicken. Moreover, some of core regions exactly overlapped with genes excavated in our previous GWAS, suggesting that these genes have undergone positive selection may affect egg production. Findings in our study could draw a comparatively integrate genome-wide map of selection signature in the chicken genome, and would be worthy for explicating the genetic mechanisms of phenotypic diversity in poultry breeding.     

Microbial diversity--insights from population genetics

Although many environmental microbial populations are large and genetically diverse, both the level of diversity and the extent to which it is ecologically relevant remain enigmatic. Because the effective (or long-term) population size, N_e, is one of the parameters that determines population genetic diversity, tests and simulations that assume selectively neutral mutations may help to identify the processes that have shaped microbial diversity. Using ecologically important genes, tests of selective neutrality suggest that adaptive as well as non-adaptive types of selection act and that departure from neutrality may be widespread or restricted to small groups of genotypes. Population genetic simulations using population sizes between 10³ and 10⁷ suggest extremely high levels of microbial diversity in environments that sustain large populations. However, census and effective population sizes may differ considerably, and because we know nothing of the evolutionary history of environmental microbial populations, we also have no idea what N_e of environmental populations is. On the one hand, this reflects our ignorance of the microbial world. On the other hand, the tests and simulations illustrate interactions between microbial diversity and microbial population genetics that should inform our thinking in microbial ecology. Because of the different views on microbial diversity across these disciplines, such interactions are crucial if we are to understand the role of genes in microbial communities.     

Three Signatures of Adaptive Polymorphism Exemplified by Malaria-Associated Genes

Malaria has been one of the strongest selective pressures on our species. Many of the best-characterized cases of adaptive evolution in humans are in genes tied to malaria resistance. However, the complex evolutionary patterns at these genes are poorly captured by standard scans for nonneutral evolution. Here, we present three new statistical tests for selection based on population genetic patterns that are observed more than once among key malaria resistance loci. We assess these tests using forward-time evolutionary simulations and apply them to global whole-genome sequencing data from humans, and thus we show that they are effective at distinguishing selection from neutrality. Each test captures a distinct evolutionary pattern, here called Divergent Haplotypes, Repeated Shifts, and Arrested Sweeps, associated with a particular period of human prehistory. We clarify the selective signatures at known malaria-relevant genes and identify additional genes showing similar adaptive evolutionary patterns. Among our top outliers, we see a particular enrichment for genes involved in erythropoiesis and for genes previously associated with malaria resistance, consistent with a major role for malaria in shaping these patterns of genetic diversity. Polymorphisms at these genes are likely to impact resistance to malaria infection and contribute to ongoing host–parasite coevolutionary dynamics.     

Drift and selection influence geographic variation at immune loci of prairie-chickens

Previous studies of immunity in wild populations have focused primarily on genes of the major histocompatibility complex (MHC); however, studies of model species have identified additional immune-related genes that also affect fitness. In this study, we sequenced five non-MHC immune genes in six greater prairie-chicken (Tympanuchus cupido) populations that have experienced varying degrees of genetic drift as a consequence of population bottlenecks and fragmentation. We compared patterns of geographic variation at the immune genes with six neutral microsatellite markers to investigate the relative effects of selection and genetic drift. Global F(ST) outlier tests identified positive selection on just one of five immune genes (IAP-1) in one population. In contrast, at other immune genes, standardized G'(ST) values were lower than those at microsatellites for a majority of pairwise population comparisons, consistent with balancing selection or with species-wide positive or purifying selection resulting in similar haplotype frequencies across populations. The effects of genetic drift were also evident as summary statistics (e.g., Tajima's D) did not differ from neutrality for the majority of cases, and immune gene diversity (number of haplotypes per gene) was correlated positively with population size. In summary, we found that both genetic drift and selection shaped variation at the five immune genes, and the strength and type of selection varied among genes. Our results caution that neutral forces, such as drift, can make it difficult to detect current selection on genes.     

Skin pigmentation,biogeographical ancestry and admixture mapping

Ancestry informative markers (AIMs) are genetic loci showing alleles with large frequency differences between populations. AIMs can be used to estimate biogeographical ancestry at the level of the population, subgroup (e.g. cases and controls) and individual. Ancestry estimates at both the subgroup and individual level can be directly instructive regarding the genetics of the phenotypes that differ qualitatively or in frequency between populations. These estimates can provide a compelling foundation for the use of admixture mapping (AM) methods to identify the genes underlying these traits. We present details of a panel of 34 AIMs and demonstrate how such studies can proceed, by using skin pigmentation as a model phenotype. We have genotyped these markers in two population samples with primarily African ancestry, viz. African Americans from Washington D.C. and an African Caribbean sample from Britain, and in a sample of European Americans from Pennsylvania. In the two African population samples, we observed significant correlations between estimates of individual ancestry and skin pigmentation as measured by reflectometry (R(2)=0.21, P<0.0001 for the African-American sample and R(2)=0.16, P<0.0001 for the British African-Caribbean sample). These correlations confirm the validity of the ancestry estimates and also indicate the high level of population structure related to admixture, a level that characterizes these populations and that is detectable by using other tests to identify genetic structure. We have also applied two methods of admixture mapping to test for the effects of three candidate genes (TYR, OCA2, MC1R) on pigmentation. We show that TYR and OCA2 have measurable effects on skin pigmentation differences between the west African and west European parental populations. This work indicates that it is possible to estimate the individual ancestry of a person based on DNA analysis with a reasonable number of well-defined genetic markers. The implications and applications of ancestry estimates in biomedical research are discussed.     

整合冠心病风险SNPs功能及分子网络特征筛选疾病易感基因

目的:冠心病(Coronary Heart Disease,CHD)是一种由多因素(遗传因素、环境因素以及它们之间的相互作用)引起的复杂疾病。本文从遗传因素和分子互作模式识别新的冠心病易感基因。方法:结合冠心病群体遗传SNPs数据和PPI数据,通过群体遗传数据的风险评估、功能SNPs的判定和PPI网络基因的分类,以功能SNPs属性、网络拓扑属性和基因功能属性为特征,利用两步分类的方法筛选新的冠心病易感基因。结果:获得了69个新的冠心病易感基因,其中43个被文献证实与冠心病的发生发展密切相关,且识别的新的易感基因注释的KEGG通路中有很多是已知的易感基因所没有注释到的,如MAPK signaling pathway,Calcium signaling pathway,Focal adhesion和Chemokine signaling pathway等,其中Chemokine signaling pathway被证实是CHD发展的关键通路。结论:应用本文提出的整合筛选策略,能识别与冠心病相关的新的易感基因,可为冠心病的预防、诊断和治疗提供新的研究方向。     

A Drosophila overexpression screen for modifiers of Rho signalling in cytokinesis

To identify genes that modulate Rho signalling during cytokinesis we tested the effect of overexpressing a set of 2190 genes on an eye phenotype caused by defective Rho activation. The resulting 112 modifier loci fell into three main classes: cell cycle genes, signalling effectors and metabolic enzymes. We developed a further series of genetic tests to refine the interactors into those most likely to modify Rho signalling during cytokinesis. In addition to a number of genes previously implicated in the Rho pathway during cytokinesis, we identified four novel primary candidates: cdc14, Pitslre, PDK1 and thread/diap1. cdc14 orthologs have, however, been implicated in cytokinesis in other organisms, as have molecules related to Thread/Diap1. The identification of several modifiers that are genetically redundant paralogs highlights the ability of overexpression screens to identify genes that are refractory to traditional loss-of-function approaches. Overexpression screens and sensitized phenotypes, therefore, may help identify the many factors that are expected to be involved in cytokinesis but have not been discovered by previous genetic screens.     

Genetic evidence for a distinct subtype of schizophrenia characterized by pervasive cognitive deficit

A novel phenotyping strategy in schizophrenia, targeting different neurocognitive domains, neurobehavioral features, and selected personality traits, has allowed us to identify a homogeneous familial subtype of the disease, characterized by pervasive neurocognitive deficit. Our genome scan data indicate that this subtype, which accounts for up to 50% of our sample, has a distinct genetic basis and explains linkage to chromosome 6p24 reported previously. If representative of other populations, the ratio of schizophrenia subtypes observed in our families could have a profound impact on sample heterogeneity and on the power of genetic studies to detect linkage and association. Our proposed abbreviated battery of tests should facilitate phenotype characterization for future genetic analyses and allow a focus on a crisply defined schizophrenia subtype, thus promoting a more informed search for susceptibility genes.     

A Population Genetic Signal of Polygenic Adaptation

Adaptation in response to selection on polygenic phenotypes may occur via subtle allele frequencies shifts at many loci. Current population genomic techniques are not well posed to identify such signals. In the past decade, detailed knowledge about the specific loci underlying polygenic traits has begun to emerge from genome-wide association studies (GWAS). Here we combine this knowledge from GWAS with robust population genetic modeling to identify traits that may have been influenced by local adaptation. We exploit the fact that GWAS provide an estimate of the additive effect size of many loci to estimate the mean additive genetic value for a given phenotype across many populations as simple weighted sums of allele frequencies. We use a general model of neutral genetic value drift for an arbitrary number of populations with an arbitrary relatedness structure. Based on this model, we develop methods for detecting unusually strong correlations between genetic values and specific environmental variables, as well as a generalization of comparisons to test for over-dispersion of genetic values among populations. Finally we lay out a framework to identify the individual populations or groups of populations that contribute to the signal of overdispersion. These tests have considerably greater power than their single locus equivalents due to the fact that they look for positive covariance between like effect alleles, and also significantly outperform methods that do not account for population structure. We apply our tests to the Human Genome Diversity Panel (HGDP) dataset using GWAS data for height, skin pigmentation, type 2 diabetes, body mass index, and two inflammatory bowel disease datasets. This analysis uncovers a number of putative signals of local adaptation, and we discuss the biological interpretation and caveats of these results.     

A hierarchical Bayesian model for next-generation population genomics

The demography of populations and natural selection shape genetic variation across the genome and understanding the genomic consequences of these evolutionary processes is a fundamental aim of population genetics. We have developed a hierarchical Bayesian model to quantify genome-wide population structure and identify candidate genetic regions affected by selection. This model improves on existing methods by accounting for stochastic sampling of sequences inherent in next-generation sequencing (with pooled or indexed individual samples) and by incorporating genetic distances among haplotypes in measures of genetic differentiation. Using simulations we demonstrate that this model has a low false-positive rate for classifying neutral genetic regions as selected genes (i.e., Φ(ST) outliers), but can detect recent selective sweeps, particularly when genetic regions in multiple populations are affected by selection. Nonetheless, selection affecting just a single population was difficult to detect and resulted in a high false-negative rate under certain conditions. We applied the Bayesian model to two large sets of human population genetic data. We found evidence of widespread positive and balancing selection among worldwide human populations, including many genetic regions previously thought to be under selection. Additionally, we identified novel candidate genes for selection, several of which have been linked to human diseases. This model will facilitate the population genetic analysis of a wide range of organisms on the basis of next-generation sequence data.     

Isolation by environment in White‐breasted Nuthatches (Sitta carolinensis) of the Madrean Archipelago sky islands: a landscape genomics approach

Understanding landscape processes driving patterns of population genetic differentiation and diversity has been a long‐standing focus of ecology and evolutionary biology. Gene flow may be reduced by historical, ecological or geographic factors, resulting in patterns of isolation by distance (IBD) or isolation by environment (IBE). Although IBE has been found in many natural systems, most studies investigating patterns of IBD and IBE in nature have used anonymous neutral genetic markers, precluding inference of selection mechanisms or identification of genes potentially under selection. Using landscape genomics, the simultaneous study of genomic and ecological landscapes, we investigated the processes driving population genetic patterns of White‐breasted Nuthatches (Sitta carolinensis) in sky islands (montane forest habitat islands) of the Madrean Archipelago. Using more than 4000 single nucleotide polymorphisms and multiple tests to investigate the relationship between genetic differentiation and geographic or ecological distance, we identified IBE, and a lack of IBD, among sky island populations of S. carolinensis. Using three tests to identify selection, we found 79 loci putatively under selection; of these, seven matched CDS regions in the Zebra Finch. The loci under selection were highly associated with climate extremes (maximum temperature of warmest month and minimum precipitation of driest month). These results provide evidence for IBE – disentangled from IBD – in sky island vertebrates and identify potential adaptive genetic variation.     

Microarray genotyping resource to determine population stratification in genetic association studies of complex disease

We have developed a robust microarray genotyping chip that will help advance studies in genetic epidemiology. In population-based genetic association studies of complex disease, there could be hidden genetic substructure in the study populations, resulting in false-positive associations. Such population stratification may confound efforts to identify true associations between genotype/haplotype and phenotype. Methods relying on genotyping additional null single nucleotide polymorphism (SNP) markers have been proposed, such as genomic control (GC) and structured association (SA), to correct association tests for population stratification. If there is an association of a disease with null SNPs, this suggests that there is a population subset with different genetic background plus different disease susceptibility. Genotyping over 100 null SNPs in the large numbers of patient and control DNA samples that are required in genetic association studies can be prohibitively expensive. We have therefore developed and tested a resequencing chip based on arrayed primer extension (APEX) from over 2000 DNA probe features that facilitate multiple interrogations of each SNP, providing a powerful, accurate, and economical means to simultaneously determine the genotypes at 110 null SNP loci in any individual. Based on 1141 known genotypes from other research groups, our GC SNP chip has an accuracy of 98.5%, including non-calls.     

A roadmap for the genetic analysis of renal aging

Several studies show evidence for the genetic basis of renal disease, which renders some individuals more prone than others to accelerated renal aging. Studying the genetics of renal aging can help us to identify genes involved in this process and to unravel the underlying pathways. First, this opinion article will give an overview of the phenotypes that can be observed in age‐related kidney disease. Accurate phenotyping is essential in performing genetic analysis. For kidney aging, this could include both functional and structural changes. Subsequently, this article reviews the studies that report on candidate genes associated with renal aging in humans and mice. Several loci or candidate genes have been found associated with kidney disease, but identification of the specific genetic variants involved has proven to be difficult. CUBN, UMOD, and SHROOM3 were identified by human GWAS as being associated with albuminuria, kidney function, and chronic kidney disease (CKD). These are promising examples of genes that could be involved in renal aging, and were further mechanistically evaluated in animal models. Eventually, we will provide approaches for performing genetic analysis. We should leverage the power of mouse models, as testing in humans is limited. Mouse and other animal models can be used to explain the underlying biological mechanisms of genes and loci identified by human GWAS. Furthermore, mouse models can be used to identify genetic variants associated with age‐associated histological changes, of which Far2, Wisp2, and Esrrg are examples. A new outbred mouse population with high genetic diversity will facilitate the identification of genes associated with renal aging by enabling high‐resolution genetic mapping while also allowing the control of environmental factors, and by enabling access to renal tissues at specific time points for histology, proteomics, and gene expression.     

Suppression mechanisms: themes from variations.

Suppressor analysis is a commonly used strategy to identify functional relationships between genes that might not have been revealed through other genetic or biochemical means. Many mechanisms that explain the phenomenon of genetic suppression have been described, but the wide variety of possible mechanisms can present a challenge to defining the relationship between a suppressor and the original gene. This article provides a broad framework for classifying suppression mechanisms and describes a series of genetic tests that can be applied to determine the most likely mechanism of suppression.     

Recent population decline and selection shape diversity of taxol-related genes

Taxanes are defensive metabolites produced by Taxus species (yews) and used in anticancer therapies. Despite their medical interest, patterns of natural diversity in taxane-related genes are unknown. We examined variation at five main genes of Taxus baccata in the Iberian Peninsula, a region where unique yew genetic resources are endangered. We looked at several gene features and applied complementary neutrality tests, including diversity/divergence tests, tests solely based on site frequency spectrum (SFS) and Zeng's compound tests. To account for specific demography, microsatellite data were used to infer historical changes in population size based on an Approximate Bayesian Computation (ABC) approach. Polymorphism-divergence tests pointed to positive selection for genes TBT and TAT and balancing selection for DBAT. In addition, neutrality tests based on SFS found that while a recent reduction in population size may explain most statistics' values, selection may still be in action in genes TBT and DBAT, at least in some populations. Molecular signatures on taxol genes suggest the action of frequent selective waves with different direction or intensity, possibly related to varying adaptive pressures produced by the host-enemy co-evolution on defence-related genes. Such natural selection processes may have produced taxane variants still undiscovered.     

Population Genetic Assignment of Confiscated Gopher Tortoises

Abstract As gopher tortoises (Gopherus polyphemus) increasingly become threatened throughout their range in Florida, USA, the need for management and conservation will intensify. Here we evaluate the forensic applicability of genetic assignment tests based on microsatellite genotypic data to 1) accurately assign individuals in our genetic database to the sample location or population of origin and 2) determine the origin of 6 confiscated tortoises. Overall, we could correctly assign 90% of the individuals in the database to their population of origin, but we were unable to determine the source of the confiscated tortoises. However, these individuals are unlikely to have come from any of our sampled sites and all 6 may have come from the same population. This approach can be used by law enforcement personnel to identify the origin of confiscated tortoises as well as by developers and wildlife managers to determine the genetic appropriateness of potential recipient populations when it is necessary to relocate individuals.     

Migratory Flyways May Affect Population Structure in Double-Crested Cormorants

Double-crested cormorants (Phalacrocorax auritus) recovered from a demographic bottleneck so well that they are now considered a nuisance species at breeding and wintering grounds across the United States and Canada. Management of this species could be improved by refining genetic population boundaries and assigning individuals to their natal population. Further, recent radio-telemetry data suggest the existence of Interior and Atlantic migratory flyways, which could reduce gene flow and result in substantial genetic isolation. In this study, we used 1,784 individuals collected across the eastern United States, a large panel of microsatellite markers developed for this species, and individuals banded as chicks and recaptured as adults to explore the effects of migratory flyways on population structure, quantify the genetic effects of demographic bottlenecks, and determine whether individuals could be assigned to their natal population based on genotype. We found evidence for genetic population division only along migratory flyways, no evidence of genetic bottlenecks, and mixed effectiveness of assignment tests. Our population structure findings suggest that gene flow is high across large scales; for example, individuals from New York, Minnesota, and Alabama are all in panmixia. We also found that traditional subspecies ranges may not be valid because >1 subspecies was present in single genetic populations. The lack of evidence for genetic bottlenecks also likely underscores the vagility of this species, suggesting that even during demographic bottlenecks, populations were not isolated from allelic exchange. Finally, the failure of assignment tests to consistently perform is likely due in part to imperfect a priori sampling of Atlantic and Interior chicks and the high vagility of adults. We conclude that the demographic bottleneck is not likely to have reduced genetic diversity, and that assignment tests remain unreliable for this species. We recommend double-crested cormorants be managed by flyway. Further development of genomic resources in this species could improve population subdivision resolution, improve assignment tests, and reveal further information on demographic histories. © 2020 The Wildlife Society.     

A semiparametric test to detect associations between quantitative traits and candidate genes in structured populations

MOTIVATION: Although population-based association mapping may be subject to the bias caused by population stratification, alternative methods that are robust to population stratification such as family-based linkage analysis have lower mapping resolution. Recently, various statistical methods robust to population stratification were proposed for association studies, using unrelated individuals to identify associations between candidate genes and traits of interest. The association between a candidate gene and a quantitative trait is often evaluated via a regression model with inferred population structure variables as covariates, where the residual distribution is customarily assumed to be from a symmetric and unimodal parametric family, such as a Gaussian, although this may be inappropriate for the analysis of many real-life datasets. RESULTS: In this article, we proposed a new structured association (SA) test. Our method corrects for continuous population stratification by first deriving population structure and kinship matrices through a set of random genetic markers and then modeling the relationship between trait values, genotypic scores at a candidate marker and genetic background variables through a semiparametric model, where the error distribution is modeled as a mixture of Polya trees centered around a normal family of distributions. We compared our model to the existing SA tests in terms of model fit, type I error rate, power, precision and accuracy by application to a real dataset as well as simulated datasets.     

Population genetics of Plasmodium resistance genes in Anopheles gambiae: no evidence for strong selection

Anopheles mosquitoes are the primary vectors for malaria in Africa, transmitting the disease to more than 100 million people annually. Recent functional studies have revealed mosquito genes that are crucial for Plasmodium development, but there is presently little understanding of which genes mediate vector competence in the wild, or evolve in response to parasite-mediated selection. Here, we use population genetic approaches to study the strength and mode of natural selection on a suite of mosquito immune system genes, CTL4, CTLMA2, LRIM1, and APL2 (LRRD7), which have been shown to affect Plasmodium development in functional studies. We sampled these genes from two African populations of An. gambiae s.s., along with several closely related species, and conclude that there is no evidence for either strong directional or balancing selection on these genes. We highlight a number of challenges that need to be met in order to apply population genetic tests for selection in Anopheles mosquitoes; in particular the dearth of suitable outgroup species and the potential difficulties that arise when working within a closely-related species complex.