A bioinformatics workflow for detecting signatures of selection in genomic data

Shared rare genetic variants in multiplex autism families suggest a social memory gene under selection

Autism spectrum disorder (ASD) affects up to 1 in 59 children, and is one of the most common neurodevelopmental disorders. Recent genomic studies have highlighted the role of rare variants in ASD. This study aimed to identify genes affected by rare variants shared by siblings with ASD and validate the function of a candidate gene FRRS1L. By integrating the whole genome sequencing data of 866 multiplex families from the Hartwell Foundation's Autism Research and Technology Initiative and Autism Speaks MSSNG project, we identified rare variants shared by two or more siblings with ASD. Using shared rare variants (SRVs), we selected candidate genes for ASD. Gene prioritization by evolutionary features and expression alterations on autism identified FRRS1L in two families, including one with impaired social behaviors. One variant in this family was 6 bp away from human-specific trinucleotide fixation. Additionally, CRISPR/Cas9 experiments demonstrated downregulation by a family variant and upregulation by a fixed site. Population genetics further demonstrated that upregulation of this gene has been favored during human evolution. Various mouse behavioral tests showed that Frrs1l knockout specifically impairs social novelty recognition without altering other behavioral phenotypes. Furthermore, we constructed humanized mice by introducing human sequences into a mouse genome. These knockin mice showed improved abilities to retain social memory over time. The results of our population genetic and evolutionary analyses, behavior experiments, and genome editing propose a molecular mechanism that may confer a selective advantage through social memory enhancement and may cause autism-related social impairment when disrupted in humans. These findings highlight the role of FRRS1L, the AMPA receptor subunit, in social behavior and evolution.

Selection signatures for high altitude adaptation in livestock: A review

High altitude adapted livestock species (cattle, yak, goat, sheep, and horse) has critical role in the human socioeconomic sphere and acts as good source of animal source products including milk, meat, and leather, among other things. These species sustain production and reproduction even in harsh environments on account of adaptation resulting from continued evolution of beneficial traits. Selection pressure leads to various adaptive strategies in livestock whose footprints are evident at the different genomic sites as the "Selection Signature". Scrutiny of these signatures provides us crucial insight into the evolutionary process and domestication of livestock adapted to diverse climatic conditions. These signatures have the potential to change the sphere of animal breeding and further usher the selection programmes in right direction. Technological revolution and recent strides made in genomic studies has opened the routes for the identification of selection signatures. Numerous statistical approaches and bioinformatics tools have been developed to detect the selection signature. Consequently, studies across years have identified candidate genes under selection region found associated with numerous traits which have a say in adaptation to high-altitude environment. This makes it pertinent to have a better understanding about the selection signature, the ways to identify and how to utilize them for betterment of livestock populations as well as farmers. This review takes a closer look into the general concept, various methodologies, and bioinformatics tools commonly employed in selection signature studies and summarize the results of recent selection signature studies related to high-altitude adaptation in various livestock species. This review will serve as an informative and useful insight for researchers and students in the field of animal breeding and evolutionary biology.

Assessment of genetic diversity, inbreeding and collection completeness of Jersey bulls in the US National Animal Germplasm Program

Genomic selection and extensive use of a few elite bulls through artificial insemination are leading to reduced genetic diversity in Jersey cattle. Conservation of genetic diversity through gene banks can protect a breed's genetic diversity and genetic gain, ensuring continued genetic advancement in the future. The availability of genomic information in the US National Animal Germplasm Program (NAGP) facilitates characterization of Jersey bulls in the germplasm collection. Therefore, in this study, we compared the genetic diversity and inbreeding between Jersey bulls in the NAGP and the national cooperator database (NCD). The NCD is maintained and curated by the Council on Dairy Cattle Breeding (CDCB). We found the genetic diversity to be marginally higher in NAGP (Ho = 0.34 ± 0.17) relative to the NCD population (Ho = 0.33 ± 0.16). The average pedigree and genomic inbreeding (FPED, FGRM, FROH > 2Mb) were similar between the groups, with estimates of 7.6% with FPED, 11.07% with FGRM and 20.13% with FROH > 2Mb. An increasing trend in inbreeding was detected, and a significantly higher level of inbreeding was estimated among the older bulls in the NAGP collection, suggesting an overrepresentation of the genetics from elite bulls. Results from principal component analyses (PCA) provided evidence that the NAGP collection is representative of the genetic variation found in the NCD population and a broad majority of the loci segregating (98.2%) in the NCD population were also segregating in the NAGP. Ward's clustering was used to assess collection completeness of Jerseys in the NAGP by comparison with top 1000 sires of bulls, top 1000 sires of cow, and bulls with high Lifetime Net Merit (NM$). All the clusters were represented in the NAGP suggesting that most of the genetic diversity in the US Jersey population is represented in the NAGP and confirmed the PCA results. The decade of birth was the major driver grouping bulls into clusters, suggesting the importance of selection over time. Selection signature analysis between the historic bulls in the NAGP with the newer bulls, born in the decade after implementation of genomic selection, identified selection for milk production, fat and protein yield, fertility, health, and reproductive traits. Cluster analysis revealed that the NAGP has captured allele frequency changes over time associated with selection, validating the strategy of repeated sampling and suggests that the continuation of a repeated sampling policy is essential for the germplasm collection to maintain its future utility. While NAGP should continue to collect bulls that have large influence on the population due to selection, care should be taken to include the entire breadth of bulls, including low merit bulls.

Landscape genomics reveals regions associated with adaptive phenotypic and genetic variation in Ethiopian indigenous chickens

Climate change is a threat to sustainable livestock production and livelihoods in the tropics. It has adverse impacts on feed and water availability, disease prevalence, production, environmental temperature, and biodiversity. Unravelling the drivers of local adaptation and understanding the underlying genetic variation in random mating indigenous livestock populations informs the design of genetic improvement programmes that aim to increase productivity and resilience. In the present study, we combined environmental, genomic, and phenotypic information of Ethiopian indigenous chickens to investigate their environmental adaptability. Through a hybrid sampling strategy, we captured wide biological and ecological variabilities across the country. Our environmental dataset comprised mean values of 34 climatic, vegetation and soil variables collected over a thirty-year period for 260 geolocations. Our biological dataset included whole genome sequences and quantitative measurements (on eight traits) from 513 individuals, representing 26 chicken populations spread along 4 elevational gradients (6-7 populations per gradient). We performed signatures of selection analyses ([Formula: see text] and XP-EHH) to detect footprints of natural selection, and redundancy analyses (RDA) to determine genotype-environment and genotype-phenotype-associations. RDA identified 1909 outlier SNPs linked with six environmental predictors, which have the highest contributions as ecological drivers of adaptive phenotypic variation. The same method detected 2430 outlier SNPs that are associated with five traits. A large overlap has been observed between signatures of selection identified by[Formula: see text]and XP-EHH showing that both methods target similar selective sweep regions. Average genetic differences measured by [Formula: see text] are low between gradients, but XP-EHH signals are the strongest between agroecologies. Genes in the calcium signalling pathway, those associated with the hypoxia-inducible factor (HIF) transcription factors, and sports performance (GALNTL6) are under selection in high-altitude populations. Our study underscores the relevance of landscape genomics as a powerful interdisciplinary approach to dissect adaptive phenotypic and genetic variation in random mating indigenous livestock populations.

Genome-wide association study of drought tolerance traits in sugar beet germplasms at the seedling stage

Introduction: Sugar beets are an important crop for global sugar production. Intense drought and the increasing lack of water resources pose a great threat to sugar beet cultivation. It is a priority to investigate favourable germplasms and functional genes to improve the breeding of drought tolerant plants. Methods: Thus, in this study, 328 sugar beet germplasms were used in a genome-wide association study (GWAS) to identify single nucleotide polymorphism (SNP) markers and candidate genes associated with drought tolerance. Results: The results showed that under drought stress (9% PEG-6000), there were 11 significantly associated loci on chromosomes 2, 3, 5, 7, and 9 from the 108946 SNPs filtered using a mixed linear model (MLM). Genome-wide association analysis combined with qRT-PCR identified 13 genes that were significantly differentially expressed in drought-tolerant extreme materials. Discussion: These candidate genes mainly exhibited functions such as regulating sugar metabolism, maintaining internal environmental stability and participating in photosystem repair. This study provides valuable information for exploring the molecular mechanisms of drought tolerance and improvement in sugar beet.

Pan-genome of Citrullus genus highlights the extent of presence/absence variation during domestication and selection

The rich genetic diversity in Citrullus lanatus and the other six species in the Citrullus genus provides important sources in watermelon breeding. Here, we present the Citrullus genus pan-genome based on the 400 Citrullus genus resequencing data, showing that 477 Mb contigs and 6249 protein-coding genes were absent in the Citrullus lanatus reference genome. In the Citrullus genus pan-genome, there are a total of 8795 (30.5%) genes that exhibit presence/absence variations (PAVs). Presence/absence variation (PAV) analysis showed that a lot of gene PAV were selected during the domestication and improvement, such as 53 favorable genes and 40 unfavorable genes were identified during the C. mucosospermus to C. lanatus landrace domestication. We also identified 661 resistance gene analogs (RGAs) in the Citrullus genus pan-genome, which contains 90 RGAs (89 variable and 1 core gene) located on the pangenome additional contigs. By gene PAV-based GWAS, 8 gene presence/absence variations were found associated with flesh color. Finally, based on the results of gene PAV selection analysis between watermelon populations with different fruit colors, we identified four non-reference candidate genes associated with carotenoid accumulation, which had a significantly higher frequency in the white flesh. These results will provide an important source for watermelon breeding.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

Unraveling Signatures of Local Adaptation among Indigenous Groups from Mexico

Few studies have addressed how selective pressures have shaped the genetic structure of the current Native American populations, and they have mostly limited their inferences to admixed Latin American populations. Here, we searched for local adaptation signals, based on integrated haplotype scores and population branch statistics, in 325 Mexican Indigenous individuals with at least 99% Native American ancestry from five previously defined geographical regions. Although each region exhibited its own local adaptation profile, only PPARG and AJAP1, both negative regulators of the Wnt/β catenin signaling pathway, showed significant adaptation signals in all the tested regions. Several signals were found, mainly in the genes related to the metabolic processes and immune response. A pathway enrichment analysis revealed the overrepresentation of selected genes related to several biological phenotypes/conditions, such as the immune response and metabolic pathways, in agreement with previous studies, suggesting that immunological and metabolic pressures are major drivers of human adaptation. Genes related to the gut microbiome measurements were overrepresented in all the regions, highlighting the importance of studying how humans have coevolved with the microbial communities that colonize them. Our results provide a further explanation of the human evolutionary history in response to environmental pressures in this region.

Genomic insight into the population history and biological adaptations of high-altitude Tibetan highlanders in Nagqu

Tibetan, one of the largest indigenous populations living in the high-altitude region of the Tibetan Plateau (TP), has developed a suite of physiological adaptation strategies to cope with the extreme highland environment in TP. Here, we reported genome-wide SNP data from 48 Kham-speaking Nagqu Tibetans and analyzed it with published data from 1,067 individuals in 167 modern and ancient populations to characterize the detailed Tibetan subgroup history and population substructure. Overall, the patterns of allele sharing and haplotype sharing suggested (1) the relatively genetic homogeny between the studied Nagqu Tibetans and ancient Nepalese as well as present-day core Tibetans from Lhasa, Nagqu, and Shigatse; and (2) the close relationship between our studied Kham-speaking Nagqu Tibetans and Kham-speaking Chamdo Tibetans. The fitted qpAdm models showed that the studied Nagqu Tibetans could be fitted as having the main ancestry from late Neolithic upper Yellow River millet farmers and deeply diverged lineages from Southern East Asians (represented by Upper Paleolithic Guangxi_Longlin and Laos_Hoabinhian), and a non-neglectable western Steppe herder-related ancestry (∼3%). We further scanned the candidate genomic regions of natural selection for our newly generated Nagqu Tibetans and the published core Tibetans via FST, iHS, and XP-EHH tests. The genes overlapping with these regions were associated with essential human biological functions such as immune response, enzyme activity, signal transduction, skin development, and energy metabolism. Together, our results shed light on the admixture and evolutionary history of Nagqu Tibetan populations.

An Evolutionary Insight Into the Heterogeneous Severity Pattern of the SARS-CoV-2 Infection

The ongoing pandemic of COVID-19 has elaborated an idiosyncratic pattern of SARS-CoV-2-induced symptoms in the human host. Some populations have succumbed to the SARS-CoV-2 infection in large numbers during this pandemic, whereas others have shown a resilient side by manifesting only milder or no symptoms at all. This observation has relayed the onus of the heterogeneous pattern of SARS-CoV-2-induced critical illness among different populations to the host genetic factors. Here, the evolutionary route was explored and three genetic loci, i.e., rs10735079, rs2109069, and rs2236757, associated with COVID-19 were analyzed. Among the three, the risk allele A at genetic locus rs2236757 residing in the IFNAR2 gene was observed to have undergone recent positive selection in the African population.

Comparative genomics and signatures of selection in North Atlantic eels

Comparative genomic approaches can identify putative private and shared signatures of selection. We performed a comparative genomic study of North Atlantic eels, European eel (Anguilla Anguilla) and American eel (A. rostrata). The two sister species are nearly undistinguishable at the phenotypic level and despite a wide non-overlapping continental distribution, they spawn in partial sympatry in the Sargasso Sea. Taking advantage of the newly assembled and annotated genome, we used genome wide RAD sequencing data of 359 individuals retrieved from Sequence Nucleotide Archive and state-of-the-art statistic tests to identify putative genomic signatures of selection in North Atlantic eels. First, using the F_ST and XP-EHH methods, we detected apparent islands of divergence on a total of 7 chromosomes, particularly on chromosomes 6 and 10. Gene ontology analyses suggested candidate genes mainly related to energy production, development and regulation, which could reflect strong selection on traits related to eel migration and larval duration time. Gene effect prediction using SNPeff showed a high number of SNPs in noncoding regions, pointing to a possible regulatory role. Second, using the iHS method we detected shared regions under selection on a total of 11 chromosomes. Several hypotheses might account for the detection of shared islands of selection in North Atlantic eels, including parallel evolution due to adaptation to similar environments and introgression. Future comparative genomic studies will be needed to further clarify the causes and consequences of introgression, including the directionality of these introgression events.

Detecting selection using extended haplotype homozygosity (EHH)-based statistics in unphased or unpolarized data

Analysis of population genetic data often includes a search for genomic regions with signs of recent positive selection. One of such approaches involves the concept of extended haplotype homozygosity (EHH) and its associated statistics. These statistics typically require phased haplotypes, and some of them necessitate polarized variants. Here, we unify and extend previously proposed modifications to loosen these requirements. We compare the modified versions with the original ones by measuring the false discovery rate in simulated whole-genome scans and by quantifying the overlap of inferred candidate regions in empirical data. We find that phasing information is indispensable for accurate estimation of within-population statistics (for all but very large samples) and of cross-population statistics for small samples. Ancestry information, in contrast, is of lesser importance for both types of statistic. Our publicly available R package rehh incorporates the modified statistics presented here.

Evolutionary selection of alleles in the melanophilin gene that impacts on prostate organ function and cancer risk

Background and objectives

Several hundred inherited genetic variants or SNPs that alter the risk of cancer have been identified through genome-wide association studies. In populations of European ancestry, these variants are mostly present at relatively high frequencies. To gain insight into evolutionary origins, we screened a series of genes and SNPs linked to breast or prostate cancer for signatures of historical positive selection.

Methodology

We took advantage of the availability of the 1000 genome data and we performed genomic scans for positive selection in five different Caucasian populations as well as one African reference population. We then used prostate organoid cultures to provide a possible functional explanation for the interplay between the action of evolutionary forces and the disease risk association.

Results

Variants in only one gene showed genomic signatures of positive, evolutionary selection within Caucasian populations melanophilin (MLPH). Functional depletion of MLPH in prostate organoids, by CRISPR/Cas9 mutation, impacted lineage commitment of progenitor cells promoting luminal versus basal cell differentiation and on resistance to androgen deprivation.

Conclusions and implications

The MLPH variants influencing prostate cancer risk may have been historically selected for their adaptive benefit on skin pigmentation but MLPH is highly expressed in the prostate and the derivative, positively selected, alleles decrease the risk of prostate cancer. Our study suggests a potential functional mechanism via which MLPH and its genetic variants could influence risk of prostate cancer, as a serendipitous consequence of prior evolutionary benefits to another tissue.

Lay Summary

We screened a limited series of genomic variants associated with breast and prostate cancer risk for signatures of historical positive selection. Variants within the melanophilin (MLPH) gene fell into this category. Depletion of MLPH in prostate organoid cultures, suggested a potential functional mechanism for impacting on cancer risk, as a serendipitous consequence of prior evolutionary benefits to another tissue.

Genetic diversity of 'Very Important Pharmacogenes' in two South-Asian populations

OBJECTIVES

Reliable identification of population-specific variants is important for building the single nucleotide polymorphism (SNP) profile. In this study, genomic variation using allele frequency differences of pharmacologically important genes for Gujarati Indians in Houston (GIH) and Indian Telugu in the U.K. (ITU) from the 1000 Genomes Project vis-à-vis global population data was studied to understand its role in drug response.

METHODS

Joint genotyping approach was used to derive variants of GIH and ITU independently. SNPs of both these populations with significant allele frequency variation (minor allele frequency ≥ 0.05) with super-populations from the 1000 Genomes Project and gnomAD based on Chi-square distribution with p-value of ≤ 0.05 and Bonferroni's multiple adjustment tests were identified. Population stratification and fixation index analysis was carried out to understand genetic differentiation. Functional annotation of variants was carried out using SnpEff, VEP and CADD score.

RESULTS

Population stratification of VIP genes revealed four clusters viz., single cluster of GIH and ITU, one cluster each of East Asian, European, African populations and Admixed American was found to be admixed. A total of 13 SNPs belonging to ten pharmacogenes were identified to have significant allele frequency variation in both GIH and ITU populations as compared to one or more super-populations. These SNPs belong to VKORC1 (rs17708472, rs2359612, rs8050894) involved in Vitamin K cycle, cytochrome P450 isoforms CYP2C9 (rs1057910), CYP2B6 (rs3211371), CYP2A2 (rs4646425) and CYP2A4 (rs4646440); ATP-binding cassette (ABC) transporter ABCB1 (rs12720067), DPYD1 (rs12119882, rs56160474) involved in pyrimidine metabolism, methyltransferase COMT (rs9332377) and transcriptional factor NR1I2 (rs6785049). SNPs rs1544410 (VDR), rs2725264 (ABCG2), rs5215 and rs5219 (KCNJ11) share high fixation index (≥ 0.5) with either EAS/AFR populations. Missense variants rs1057910 (CYP2C9), rs1801028 (DRD2) and rs1138272 (GSTP1), rs116855232 (NUDT15); intronic variants rs1131341 (NQO1) and rs115349832 (DPYD) are identified to be 'deleterious'.

CONCLUSIONS

Analysis of SNPs pertaining to pharmacogenes in GIH and ITU populations using population structure, fixation index and allele frequency variation provides a premise for understanding the role of genetic diversity in drug response in Asian Indians.

Genome-wide scans for signatures of selection in Mangalarga Marchador horses using high-throughput SNP genotyping

BACKGROUND

The detection of signatures of selection in genomic regions provides insights into the evolutionary process, enabling discoveries regarding complex phenotypic traits. In this research, we focused on identifying genomic regions affected by different selection pressures, mainly highlighting the recent positive selection, as well as understanding the candidate genes and functional pathways associated with the signatures of selection in the Mangalarga Marchador genome. Besides, we seek to direct the discussion about genes and traits of importance in this breed, especially traits related to the type and quality of gait, temperament, conformation, and locomotor system.

RESULTS

Three different methods were used to search for signals of selection: Tajima's D (TD), the integrated haplotype score (iHS), and runs of homozygosity (ROH). The samples were composed of males (n = 62) and females (n = 130) that were initially chosen considering well-defined phenotypes for gait: picada (n = 86) and batida (n = 106). All horses were genotyped using a 670 k Axiom® Equine Genotyping Array​ (Axiom MNEC670). In total, 27, 104 (chosen), and 38 candidate genes were observed within the signatures of selection identified in TD, iHS, and ROH analyses, respectively. The genes are acting in essential biological processes. The enrichment analysis highlighted the following functions: anterior/posterior pattern for the set of genes (GLI3, HOXC9, HOXC6, HOXC5, HOXC4, HOXC13, HOXC11, and HOXC10); limb morphogenesis, skeletal system, proximal/distal pattern formation, JUN kinase activity (CCL19 and MAP3K6); and muscle stretch response (MAPK14). Other candidate genes were associated with energy metabolism, bronchodilator response, NADH regeneration, reproduction, keratinization, and the immunological system.

CONCLUSIONS

Our findings revealed evidence of signatures of selection in the MM breed that encompass genes acting on athletic performance, limb development, and energy to muscle activity, with the particular involvement of the HOX family genes. The genome of MM is marked by recent positive selection. However, Tajima's D and iHS results point also to the presence of balancing selection in specific regions of the genome.

Unveiling Comparative Genomic Trajectories of Selection and Key Candidate Genes in Egg-Type Russian White and Meat-Type White Cornish Chickens

Comparison of genomic footprints in chicken breeds with different selection history is a powerful tool in elucidating genomic regions that have been targeted by recent and more ancient selection. In the present work, we aimed at examining and comparing the trajectories of artificial selection in the genomes of the native egg-type Russian White (RW) and meat-type White Cornish (WC) breeds. Combining three different statistics (top 0.1% SNP by F_ST value at pairwise breed comparison, hapFLK analysis, and identification of ROH island shared by more than 50% of individuals), we detected 45 genomic regions under putative selection including 11 selective sweep regions, which were detected by at least two different methods. Four of such regions were breed-specific for each of RW breed (on GGA1, GGA5, GGA8, and GGA9) and WC breed (on GGA1, GGA5, GGA8, and GGA28), while three remaining regions on GGA2 (two sweeps) and GGA3 were common for both breeds. Most of identified genomic regions overlapped with known QTLs and/or candidate genes including those for body temperatures, egg productivity, and feed intake in RW chickens and those for growth, meat and carcass traits, and feed efficiency in WC chickens. These findings were concordant with the breed origin and history of their artificial selection. We determined a set of 188 prioritized candidate genes retrieved from the 11 overlapped regions of putative selection and reviewed their functions relative to phenotypic traits of interest in the two breeds. One of the RW-specific sweep regions harbored the known domestication gene, TSHR. Gene ontology and functional annotation analysis provided additional insight into a functional coherence of genes in the sweep regions. We also showed a greater candidate gene richness on microchromosomes relative to macrochromosomes in these genomic areas. Our results on the selection history of RW and WC chickens and their key candidate genes under selection serve as a profound information for further conservation of their genomic diversity and efficient breeding.

Analysis of selection signatures on the Z chromosome of bidirectional selection broiler lines for the assessment of abdominal fat content

Background

The discovery of selection signatures has enabled the identification of genomics regions under selective pressure, enhancing knowledge of evolutionary genotype-phenotypes. Sex chromosomes play an important role in species formation and evolution. Therefore, the exploration of selection signatures on sex chromosomes has important biological significance.

Results

In this study, we used the Cross Population Extend Haplotype Homozygosity Test (XPEHH), F-statistics (F_ST) and EigenGWAS to assess selection signatures on the Z chromosome in 474 broiler chickens via Illumina chicken 60 K SNP chips. SNP genotype data were downloaded from publicly available resources. We identified 17 selection regions, amongst which 1, 11 and 12 were identified by XPEHH, F_ST, and EigenGWAS, respectively. Each end of the Z chromosome appeared to undergo the highest levels of selection pressure. A total of 215 candidate genes were located in 17 selection regions, some of which mediated lipogenesis, fatty acid production, fat metabolism, and fat decomposition, including FGF10, ELOVL7, and IL6ST. Using abdominal adipose tissue expression data of the chickens, 187 candidate genes were expressed with 15 differentially expressed genes (DEGs) in fat vs. lean lines identified. Amongst the DEGs, VCAN was related to fat metabolism. GO pathway enrichment analysis and QTL annotations were performed to fully characterize the selection mechanism(s) of chicken abdominal fat content.

Conclusions

We have found some selection regions and candidate genes involving in fat metabolism on the Z chromosome. These findings enhance our understanding of sex chromosome selection signatures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00971-6.

An overview of current population genomics methods for the analysis of whole-genome resequencing data in eukaryotes

Characterizing the population history of a species and identifying loci underlying local adaptation is crucial in functional ecology, evolutionary biology, conservation and agronomy. The constant improvement of high-throughput sequencing techniques has facilitated the production of whole genome data in a wide range of species. Population genomics now provides tools to better integrate selection into a historical framework, and take into account selection when reconstructing demographic history. However, this improvement has come with a profusion of analytical tools that can confuse and discourage users. Such confusion limits the amount of information effectively retrieved from complex genomic data sets, and impairs the diffusion of the most recent analytical tools into fields such as conservation biology. It may also lead to redundancy among methods. To address these isssues, we propose an overview of more than 100 state-of-the-art methods that can deal with whole genome data. We summarize the strategies they use to infer demographic history and selection, and discuss some of their limitations. A website listing these methods is available at www.methodspopgen.com.

Genetic Landscape of Rare Autoinflammatory Disease Variants in Qatar and Middle Eastern Populations Through the Integration of Whole-Genome and Exome Datasets

Rare monogenic autoinflammatory diseases are a group of recurrent inflammatory genetic disorders caused due to genetic variants in over 37 genes. While a number of these disorders have been identified and reported in Middle Eastern populations, the carrier frequency of these genetic variants in the Middle Eastern population is not known. The availability of whole-genome and exome datasets of over 1,000 individuals from Qatar persuaded us to explore the genetic epidemiology of rare autoinflammatory genetic variants. We have systematically analyzed genetic variants in genome-scale datasets from Qatar with a compendium of variants associated with autoinflammatory diseases. The variants were systematically reclassified according to the American College of Medical Genetics and Genomics guidelines for interpretation of variant pathogenicity. Our analysis identified seven pathogenic and likely pathogenic variants with significant differences in their allele frequencies compared to the global population. The cumulative carrier frequency of these variants was found to be 2.58%. Furthermore, our analysis revealed that five genes, implicated in rare autoinflammatory diseases, were under natural selection. To the best of our knowledge, this is the first and most comprehensive study on the population-scale analysis and genetic epidemiology of genetic variants that cause rare autoinflammatory disease in Middle Eastern populations.

Genetic diversity and population structure of advanced clones selected over forty years by a potato breeding program in the USA

Knowledge regarding genetic diversity and population structure of breeding materials is essential for crop improvement. The Texas A&M University Potato Breeding Program has a collection of advanced clones selected and maintained in-vitro over a 40-year period. Little is known about its genetic makeup and usefulness for the current breeding program. In this study, 214 potato clones were genotyped with the Infinium Illumina 22 K V3 Potato Array. After filtering, a total of 10,106 single nucleotide polymorphic (SNP) markers were used for analysis. Heterozygosity varied by SNP, with an overall average of 0.59. Three groups of tetraploid clones primarily based on potato market classes, were detected using STRUCTURE software and confirmed by discriminant analysis of principal components. The highest coefficient of differentiation observed between the groups was 0.14. Signatures of selection were uncovered in genes controlling potato flesh and skin color, length of plant cycle and tuberization, and carbohydrate metabolism. A core set of 43 clones was obtained using Core Hunter 3 to develop a sub-collection that retains similar genetic diversity as the whole population, minimize redundancies, and facilitates long-term conservation of genetic resources. The comprehensive molecular characterization of our breeding clone bank collection contributes to understanding the genetic diversity of existing potato resources. This analysis could be applied to other breeding programs and assist in the selection of parents, fingerprinting, protection, and management of the breeding collections.

An Experimental Test of Adaptive Introgression in Locally Adapted Populations of Splash Pool Copepods

As species struggle to keep pace with the rapidly warming climate, adaptive introgression of beneficial alleles from closely related species or populations provides a possible avenue for rapid adaptation. We investigate the potential for adaptive introgression in the copepod, Tigriopus californicus, by hybridizing two populations with divergent heat tolerance limits. We subjected hybrids to strong heat selection for 15 generations followed by whole-genome resequencing. Utilizing a hybridize evolve and resequence (HER) technique, we can identify loci responding to heat selection via a change in allele frequency. We successfully increased the heat tolerance (measured as LT50) in selected lines, which was coupled with higher frequencies of alleles from the southern (heat tolerant) population. These repeatable changes in allele frequencies occurred on all 12 chromosomes across all independent selected lines, providing evidence that heat tolerance is polygenic. These loci contained genes with lower protein-coding sequence divergence than the genome-wide average, indicating that these loci are highly conserved between the two populations. In addition, these loci were enriched in genes that changed expression patterns between selected and control lines in response to a nonlethal heat shock. Therefore, we hypothesize that the mechanism of heat tolerance divergence is explained by differential gene expression of highly conserved genes. The HER approach offers a unique solution to identifying genetic variants contributing to polygenic traits, especially variants that might be missed through other population genomic approaches.

Selective sweep analysis reveals extensive parallel selection traits between large white and Duroc pigs

In the process of pig genetic improvement, different commercial breeds have been bred for the same purpose, improving meat production. Most of the economic traits, such as growth and fertility, have been selected similarly despite the discrepant selection pressure, which is known as parallel selection. Here, 28 whole-genome sequencing data of Danish large white pigs with an approximately 25-fold depth each were generated, resulting in about 12 million high-quality SNPs for each individual. Combined with the sequencing data of 27 Duroc and 23 European wild boars, we investigated the parallel selection of Danish large white and Duroc pigs using two complementary methods, Fst and iHS. In total, 67 candidate regions were identified as the signatures of parallel selection. The genes in candidate regions of parallel selection were mainly associated with sensory perception, growth rate, and body size. Further functional annotation suggested that the striking consistency of the terms may be caused by the polygenetic basis of quantitative traits, and revealing the complex genetic basis of parallel selection. Besides, some unique terms were enriched in population-specific selection regions, such as the limb development-related terms enriched in Duroc-specific selection regions, suggesting unique selections of breed-specific selected traits. These results will help us better understand the parallel selection process of different breeds. Moreover, we identified several potential causal SNPs that may contribute to the pig genetic breeding process.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

High Levels of Genetic Diversity within Nilo-Saharan Populations: Implications for Human Adaptation

Africa contains more human genetic variation than any other continent, but the majority of the population-scale analyses of the African peoples have focused on just two of the four major linguistic groups, the Niger-Congo and Afro-Asiatic, leaving the Nilo-Saharan and Khoisan populations under-represented. In order to assess genetic variation and signatures of selection within a Nilo-Saharan population and between the Nilo-Saharan and Niger-Congo and Afro-Asiatic, we sequenced 50 genomes from the Nilo-Saharan Lugbara population of North-West Uganda and 250 genomes from 6 previously unsequenced Niger-Congo populations. We compared these data to data from a further 16 Eurasian and African populations including the Gumuz, another putative Nilo-Saharan population from Ethiopia. Of the 21 million variants identified in the Nilo-Saharan population, 3.57 million (17%) were not represented in dbSNP and included predicted non-synonymous mutations with possible phenotypic effects. We found greater genetic differentiation between the Nilo-Saharan Lugbara and Gumuz populations than between any two Afro-Asiatic or Niger-Congo populations. F3 tests showed that Gumuz contributed a genetic component to most Niger-Congo B populations whereas Lugabara did not. We scanned the genomes of the Lugbara for evidence of selective sweeps. We found selective sweeps at four loci (SLC24A5, SNX13, TYRP1, and UVRAG) associated with skin pigmentation, three of which already have been reported to be under selection. These selective sweeps point toward adaptations to the intense UV radiation of the Sahel.

Sequencing era methods for identifying signatures of selection in the genome

Insights into genetic loci which are under selection and their functional roles contribute to increased understanding of the patterns of phenotypic variation we observe today. The availability of whole-genome sequence data, for humans and other species, provides opportunities to investigate adaptation and evolution at unprecedented resolution. Many analytical methods have been developed to interrogate these large data sets and characterize signatures of selection in the genome. We review here recently developed methods and consider the impact of increased computing power and data availability on the detection of selection signatures. Consideration of demography, recombination and other confounding factors is important, and use of a range of methods in combination is a powerful route to resolving different forms of selection in genome sequence data. Overall, a substantial improvement in methods for application to whole-genome sequencing is evident, although further work is required to develop robust and computationally efficient approaches which may increase reproducibility across studies.

Genome-Wide Selection Scan in an Arabian Peninsula Population Identifies a TNKS Haplotype Linked to Metabolic Traits and Hypertension

Despite the extreme and varying environmental conditions prevalent in the Arabian Peninsula, it has experienced several waves of human migrations following the out-of-Africa diaspora. Eventually, the inhabitants of the peninsula region adapted to the hot and dry environment. The adaptation and natural selection that shaped the extant human populations of the Arabian Peninsula region have been scarcely studied. In an attempt to explore natural selection in the region, we analyzed 662,750 variants in 583 Kuwaiti individuals. We searched for regions in the genome that display signatures of positive selection in the Kuwaiti population using an integrative approach in a conservative manner. We highlight a haplotype overlapping TNKS that showed strong signals of positive selection based on the results of the multiple selection tests conducted (integrated Haplotype Score, Cross Population Extended Haplotype Homozygosity, Population Branch Statistics, and log-likelihood ratio scores). Notably, the TNKS haplotype under selection potentially conferred a fitness advantage to the Kuwaiti ancestors for surviving in the harsh environment while posing a major health risk to present-day Kuwaitis.

Genome-Wide Analyses Reveal Footprints of Divergent Selection and Drought Adaptive Traits in Synthetic-Derived Wheats

Crop-wild introgressions have long been exploited without knowing the favorable recombination points. Synthetic hexaploid wheats are one of the most exploited genetic resources for bread wheat improvement. However, despite some QTL with major effects, much less is known about genome-wide patterns of introgressions and their effects on phenotypes. We used two genome-wide association approaches: SNP-GWAS and haplotype-GWAS to identify SNPs and haplotypes associated with productivity under water-limited conditions in a synthetic-derived wheat (SYN-DER) population. Haplotype-GWAS further enriched and identified 20 more genomic regions associated with drought adaptability that did not overlap with SNP-GWAS. Since GWAS is biased to the phenotypes in the study and may fail to detect important genetic diversity during breeding, we used five complementary analytical approaches (t-test, Tajima’s D, nucleotide diversity (π), Fst, and EigenGWAS) to identify divergent selections in SYN-DER compared to modern bread wheat. These approaches consistently pinpointed 89 ‘selective sweeps’, out of which 30 selection loci were identified on D-genome. These key selections co-localized with important functional genes of adaptive traits such as TaElf3-D1 (1D) for earliness per se (Eps), TaCKX-D1 (3D), TaGS1a (6D) and TaGS-D1 (7D) for grain size, weight and morphology, TaCwi-D1 (5D) influencing drought tolerance, and Vrn-D3 (7D) for vernalization. Furthermore, 55 SNPs and 23 haplotypes of agronomic and physiological importance such as grain yield, relative water content and thousand grain weight in SYN-DER, were among the top 5% of divergent selections contributed by synthetic hexaploid wheats. These divergent selections associated with improved agronomic performance carry new alleles that have been introduced to wheat. Our results demonstrated that GWAS and selection sweep analyses are powerful approaches for investigating favorable introgressions under strong selection pressure and the use of crop-wild hybridization to assist the improvement of wheat yield and productivity under moisture limiting environments.

Comparing signals of natural selection between three Indigenous North American populations

While many studies have highlighted human adaptations to diverse environments worldwide, genomic studies of natural selection in Indigenous populations in the Americas have been absent from this literature until very recently. Since humans first entered the Americas some 20,000 years ago, they have settled in many new environments across the continent. This diversity of environments has placed variable selective pressures on the populations living in each region, but the effects of these pressures have not been extensively studied to date. To help fill this gap, we collected genome-wide data from three Indigenous North American populations from different geographic regions of the continent (Alaska, southeastern United States, and central Mexico). We identified signals of natural selection in each population and compared signals across populations to explore the differences in selective pressures among the three regions sampled. We find evidence of adaptation to cold and high-latitude environments in Alaska, while in the southeastern United States and central Mexico, pathogenic environments seem to have created important selective pressures. This study lays the foundation for additional functional and phenotypic work on possible adaptations to varied environments during the history of population diversification in the Americas.

Demography and rapid local adaptation shape Creole cattle genome diversity in the tropics

The introduction of Iberian cattle in the Americas after Columbus' arrival imposed high selection pressures on a limited number of animals over a brief period of time. Knowledge of the genomic regions selected during this process may help in enhancing climatic resilience and sustainable animal production. We first determined taurine and indicine contributions to the genomic structure of modern Creole cattle. Second, we inferred their demographic history using approximate Bayesian computation (ABC), linkage disequilibrium (LD) and N e Slope (NeS) analysis. Third, we performed whole genome scans for selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) and population differentiation (F ST) to disentangle the genetic mechanisms involved in adaptation and phenotypic change by a rapid and major environmental transition. To tackle these questions, we combined SNP array data (~54,000 SNPs) in Creole breeds with their modern putative Iberian ancestors. Reconstruction of the population history of Creoles from the end of the 15th century indicated a major demographic expansion until the introduction of zebu and commercial breeds into the Americas ~180 years ago, coinciding with a drastic N e contraction. NeS analysis provided insights into short-term complexity in population change and depicted a decrease/expansion episode at the end of the ABC-inferred expansion, as well as several additional fluctuations in N e with the attainment of the current small N e only towards the end of the 20th century. Selection signatures for tropical adaptation pinpointed the thermoregulatory slick hair coat region, identifying a new candidate gene (GDNF), as well as novel candidate regions involved in immune function, behavioural processes, iron metabolism and adaptation to new feeding conditions. The outcomes from this study will help in future-proofing farm animal genetic resources (FAnGR) by providing molecular tools that allow selection for improved cattle performance, resilience and welfare under climate change.

Homo sapiens-Specific Binding Site Variants within Brain Exclusive Enhancers Are Subject to Accelerated Divergence across Human Population

Empirical assessments of human accelerated noncoding DNA frgaments have delineated presence of many cis-regulatory elements. Enhancers make up an important category of such accelerated cis-regulatory elements that efficiently control the spatiotemporal expression of many developmental genes. Establishing plausible reasons for accelerated enhancer sequence divergence in Homo sapiens has been termed significant in various previously published studies. This acceleration by including closely related primates and archaic human data has the potential to open up evolutionary avenues for deducing present-day brain structure. This study relied on empirically confirmed brain exclusive enhancers to avoid any misjudgments about their regulatory status and categorized among them a subset of enhancers with an exceptionally accelerated rate of lineage specific divergence in humans. In this assorted set, 13 distinct transcription factor binding sites were located that possessed unique existence in humans. Three of 13 such sites belonging to transcription factors SOX2, RUNX1/3, and FOS/JUND possessed single nucleotide variants that made them unique to H. sapiens upon comparisons with Neandertal and Denisovan orthologous sequences. These variants modifying the binding sites in modern human lineage were further substantiated as single nucleotide polymorphisms via exploiting 1000 Genomes Project Phase3 data. Long range haplotype based tests laid out evidence of positive selection to be governing in African population on two of the modern human motif modifying alleles with strongest results for SOX2 binding site. In sum, our study acknowledges acceleration in noncoding regulatory landscape of the genome and highlights functional parts within it to have undergone accelerated divergence in present-day human population.

Genome-Wide Association Study and Selection Signatures Detect Genomic Regions Associated with Seed Yield and Oil Quality in Flax

A genome-wide association study (GWAS) was performed on a set of 260 lines which belong to three different bi-parental flax mapping populations. These lines were sequenced to an averaged genome coverage of 19× using the Illumina Hi-Seq platform. Phenotypic data for 11 seed yield and oil quality traits were collected in eight year/location environments. A total of 17,288 single nucleotide polymorphisms were identified, which explained more than 80% of the phenotypic variation for days to maturity (DTM), iodine value (IOD), palmitic (PAL), stearic, linoleic (LIO) and linolenic (LIN) acid contents. Twenty-three unique genomic regions associated with 33 quantitative trait loci (QTL) for the studied traits were detected, thereby validating four genomic regions previously identified. The 33 QTL explained 48–73% of the phenotypic variation for oil content, IOD, PAL, LIO and LIN but only 8–14% for plant height, DTM and seed yield. A genome-wide selective sweep scan for selection signatures detected 114 genomic regions that accounted for 7.82% of the flax pseudomolecule and overlapped with the 11 GWAS-detected genomic regions associated with 18 QTL for 11 traits. The results demonstrate the utility of GWAS combined with selection signatures for dissection of the genetic structure of traits and for pinpointing genomic regions for breeding improvement.

Assessing genetic architecture and signatures of selection of dual purpose Gir cattle populations using genomic information

Gir is one of the main cattle breeds raised in tropical South American countries. Strong artificial selection through its domestication resulted in increased genetic differentiation among the countries in recent years. Over the years, genomic studies in Gir have become more common. However, studies of population structure and signatures of selection in divergent Gir populations are scarce and need more attention to better understand genetic differentiation, gene flow, and genetic distance. Genotypes of 173 animals selected for growth traits and 273 animals selected for milk production were used in this study. Clear genetic differentiation between beef and dairy populations was observed. Different criteria led to genetic divergence and genetic differences in allele frequencies between the two populations. Gene segregation in each population was forced by artificial selection, promoting isolation, and increasing genetic variation between them. Results showed evidence of selective forces in different regions of the genome. A total of 282 genes were detected under selection in the test population based on the fixation index (Fst), integrated haplotype score (iHS), and cross-population extend haplotype homozygosity (XP-EHH) approaches. The QTL mapping identified 35 genes associated with reproduction, milk composition, growth, meat and carcass, health, or body conformation traits. The investigation of genes and pathways showed that quantitative traits associated to fertility, milk production, beef quality, and growth were involved in the process of differentiation of these populations. These results would support further investigations of population structure and differentiation in the Gir breed.

Pharmacogenomic survey of Qatari populations using whole-genome and exome sequences

The Arabs represent one of the most genetically heterogeneous populations characterized by a high prevalence of Mendelian disorders due to consanguinity. Population-scale genomic datasets provide a unique opportunity to understand the epidemiology of variants associated with differential therapeutic response. We analyzed publicly available genomic data for 1005 Qatari individuals encompassing five subpopulations. The frequencies of known and novel pharmacogenetic variants were compared with global populations. Impact of genetic substructure on the pharmacogenetic landscape of the population was studied. We report an average of three clinically actionable pharmacogenetic variants with FDA-recommended genetic testing per Qatari individual regardless of their genetic ancestry. We observed extensive differences in the frequencies of clinically actionable pharmacogenetic variants among the Qatari subpopulations. Our analysis revealed 3579 deleterious pharmacogenetic variants potentially altering the function of 1163 genes associated with 1565 drugs. This study has thus compiled the first comprehensive landscape of pharmacogenetic variants for any Arab population.

Cardio-metabolic disease genetic risk factors among Māori and Pacific Island people in Aotearoa New Zealand: current state of knowledge and future directions

CONTEXT

Cardio-metabolic conditions in Aotearoa New Zealand (NZ) Māori and non-indigenous Polynesian (Pacific) populations have been increasing in prevalence and severity, especially over the last two decades.

OBJECTIVES

To assess knowledge on genetic and non-genetic risk factors for cardio-metabolic disease in the Māori and Pacific populations residing in Aotearoa NZ by a semi-systematic review of the PubMed database. To outline possible future directions in genetic epidemiological research with Māori and Pacific communities.

RESULTS

There have been few studies to confirm that risk factors in other populations also associate with cardio-metabolic conditions in Māori and Pacific populations. Such data are important when interventions are considered. Genetic studies have been sporadic, with no genome-wide association studies done.

CONCLUSIONS

Biomedical research with Māori and Pacific communities is important to reduce the prevalence and impact of the cardio-metabolic diseases, as precision medicine is implemented in other Aotearoa NZ populations using overseas findings. Genuine engagement with Māori and Pacific communities is needed to ensure positive outcomes for genetic studies, from data collection through to analysis and dissemination. Important is building trust, understanding by researchers of fundamental cultural concepts and implementing protocols that minimise risks and maximise benefits. Approaches that utilise information such as genealogical information and whole genome sequencing technologies will provide new insights into cardio-metabolic conditions-and new interventions for affected individuals and families.

Selection on the regulation of sympathetic nervous activity in humans and chimpanzees

Adrenergic α2C receptor (ADRA2C) is an inhibitory modulator of the sympathetic nervous system. Knockout mice for this gene show physiological and behavioural alterations that are associated with the fight-or-flight response. There is evidence of positive selection on the regulation of this gene during chicken domestication. Here, we find that the neuronal expression of ADRA2C is lower in human and chimpanzee than in other primates. On the basis of three-dimensional chromatin structure, we identified a cis-regulatory region whose DNA sequences have been significantly accelerated in human and chimpanzee. Active histone modification marks this region in rhesus macaque but not in human and chimpanzee; instead, repressive marks are enriched in various human brain samples. This region contains two neuron-restrictive silencer factor (NRSF) binding motifs, each of which harbours a polymorphism. Our genotyping and analysis of population genome data indicate that at both polymorphic sites, the derived allele has reached fixation in humans and chimpanzees but not in bonobos, whereas only the ancestral allele is present among macaques. Our CRISPR/Cas9 genome editing and reporter assays show that both derived nucleotides repress ADRA2C, most likely by increasing NRSF binding. In addition, we detected signatures of recent positive selection for lower neuronal ADRA2C expression in humans. Our findings indicate that there has been selective pressure for enhanced sympathetic nervous activity in the evolution of humans and chimpanzees.

Genome-wide characterization of genetic variants and putative regions under selection in meat and egg-type chicken lines

Background

Meat and egg-type chickens have been selected for several generations for different traits. Artificial and natural selection for different phenotypes can change frequency of genetic variants, leaving particular genomic footprints throghtout the genome. Thus, the aims of this study were to sequence 28 chickens from two Brazilian lines (meat and white egg-type) and use this information to characterize genome-wide genetic variations, identify putative regions under selection using Fst method, and find putative pathways under selection.

Results

A total of 13.93 million SNPs and 1.36 million INDELs were identified, with more variants detected from the broiler (meat-type) line. Although most were located in non-coding regions, we identified 7255 intolerant non-synonymous SNPs, 512 stopgain/loss SNPs, 1381 frameshift and 1094 non-frameshift INDELs that may alter protein functions. Genes harboring intolerant non-synonymous SNPs affected metabolic pathways related mainly to reproduction and endocrine systems in the white-egg layer line, and lipid metabolism and metabolic diseases in the broiler line. Fst analysis in sliding windows, using SNPs and INDELs separately, identified over 300 putative regions of selection overlapping with more than 250 genes. For the first time in chicken, INDEL variants were considered for selection signature analysis, showing high level of correlation in results between SNP and INDEL data. The putative regions of selection signatures revealed interesting candidate genes and pathways related to important phenotypic traits in chicken, such as lipid metabolism, growth, reproduction, and cardiac development.

Conclusions

In this study, Fst method was applied to identify high confidence putative regions under selection, providing novel insights into selection footprints that can help elucidate the functional mechanisms underlying different phenotypic traits relevant to meat and egg-type chicken lines. In addition, we generated a large catalog of line-specific and common genetic variants from a Brazilian broiler and a white egg layer line that can be used for genomic studies involving association analysis with phenotypes of economic interest to the poultry industry.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4444-0) contains supplementary material, which is available to authorized users.

Genetic epidemiology of familial Mediterranean fever through integrative analysis of whole genome and exome sequences from Middle East and North Africa

Familial Mediterranean fever (FMF), an autosomal recessive and rare autoinflammatory disease is caused by genetic mutations in the MEFV gene and is highly prevalent in the Mediterranean basin. Although the carrier frequency of specific disease variants in the MEFV gene has been reported from isolated studies, a comprehensive view of variants in the Mediterranean region has not been possible due to paucity of data. The recent availability of whole-genome and whole-exome datasets prompted us to study the genetic epidemiology of MEFV variants in the region. We assembled data from 5 datasets encompassing whole-genome and whole-exome datasets for 2115 individuals from multiple subpopulations in the region and also created a compendium for MEFV genetic variants, which were further systematically annotated as per the American College of Medical Genetics and Genomics (ACMG) guidelines. Our analysis points to significant differences in allele frequencies in the subpopulations, and the carrier frequency for MEFV genetic variants in the population to be about 8%. The MEFV gene appears to be under natural selection from our analysis. To the best of our knowledge, this is the most comprehensive study and analysis of population epidemiology of MEFV gene variants in the Middle East and North African populations.

Evidence of positive selection towards Zebuine haplotypes in the BoLA region of Brangus cattle

The Brangus breed was developed to combine the superior characteristics of both of its founder breeds, Angus and Brahman. It combines the high adaptability to tropical and subtropical environments, disease resistance, and overall hardiness of Zebu cattle with the reproductive potential and carcass quality of Angus. It is known that the major histocompatibility complex (MHC, also known as bovine leucocyte antigen: BoLA), located on chromosome 23, encodes several genes involved in the adaptive immune response and may be responsible for adaptation to harsh environments. The objective of this work was to evaluate whether the local breed ancestry percentages in the BoLA locus of a Brangus population diverged from the estimated genome-wide proportions and to identify signatures of positive selection in this genomic region. For this, 167 animals (100 Brangus, 45 Angus and 22 Brahman) were genotyped using a high-density single nucleotide polymorphism array. The local ancestry analysis showed that more than half of the haplotypes (55.0%) shared a Brahman origin. This value was significantly different from the global genome-wide proportion estimated by cluster analysis (34.7% Brahman), and the proportion expected by pedigree (37.5% Brahman). The analysis of selection signatures by genetic differentiation (F st ) and extended haplotype homozygosity-based methods (iHS and Rsb) revealed 10 and seven candidate regions, respectively. The analysis of the genes located within these candidate regions showed mainly genes involved in immune response-related pathway, while other genes and pathways were also observed (cell surface signalling pathways, membrane proteins and ion-binding proteins). Our results suggest that the BoLA region of Brangus cattle may have been enriched with Brahman haplotypes as a consequence of selection processes to promote adaptation to subtropical environments.

Genome-wide identification and evolution of the PIN-FORMED (PIN) gene family in Glycine max

Soybean (Glycine max) is one of the most important crop plants. Wild and cultivated soybean varieties have significant differences worth further investigation, such as plant morphology, seed size, and seed coat development; these characters may be related to auxin biology. The PIN gene family encodes essential transport proteins in cell-to-cell auxin transport, but little research on soybean PIN genes (GmPIN genes) has been done, especially with respect to the evolution and differences between wild and cultivated soybean. In this study, we retrieved 23 GmPIN genes from the latest updated G. max genome database; six GmPIN protein sequences were changed compared with the previous database. Based on the Plant Genome Duplication Database, 18 GmPIN genes have been involved in segment duplication. Three pairs of GmPIN genes arose after the second soybean genome duplication, and six occurred after the first genome duplication. The duplicated GmPIN genes retained similar expression patterns. All the duplicated GmPIN genes experienced purifying selection (Ka/Ks < 1) to prevent accumulation of non-synonymous mutations and thus remained more similar. In addition, we also focused on the artificial selection of the soybean PIN genes. Five artificially selected GmPIN genes were identified by comparing the genome sequence of 17 wild and 14 cultivated soybean varieties. Our research provides useful and comprehensive basic information for understanding GmPIN genes.

Identifying signatures of positive selection in pigmentation genes in two South Asian populations

OBJECTIVES

Skin pigmentation is a polygenic trait showing wide phenotypic variations among global populations. While numerous pigmentation genes have been identified to be under positive selection among European and East populations, genes contributing to phenotypic variation in skin pigmentation within and among South Asian populations are still poorly understood. The present study uses data from the Phase 3 of the 1000 genomes project focusing on two South Asian populations-GIH (Gujarati Indian from Houston, Texas) and ITU (Indian Telugu from UK), so as to decode the genetic architecture involved in adaptation to ultraviolet radiation in South Asian populations.

METHODS

Statistical tests included were (1) tests to identify deviations of the Site Frequency Spectrum (SFS) from neutral expectations (Tajima's D, Fay and Wu's H and Fu and Li's D* and F*), (2) tests focused on the identification of high-frequency haplotypes with extended linkage disequilibrium (iHS and Rsb), and (3) tests based on genetic differentiation between populations (LSBL).

RESULTS

Twenty-two pigmentation genes fall in the top 1% for at least one statistic in the GIH population, 5 of which (LYST, OCA2, SLC24A5, SLC45A2, and TYR) have been previously associated with normal variation in skin, hair, or eye color. In comparison, 17 genes fall in the top 1% for at least one statistic in the ITU population. Twelve loci which are identified as outliers in the ITU scan were also identified in the GIH population.

CONCLUSIONS

These results suggest that selection may have affected these loci broadly across the region.

SNP Discovery Using a Pangenome: Has the Single Reference Approach Become Obsolete?

Increasing evidence suggests that a single individual is insufficient to capture the genetic diversity within a species due to gene presence absence variation. In order to understand the extent to which genomic variation occurs in a species, the construction of its pangenome is necessary. The pangenome represents the complete set of genes of a species; it is composed of core genes, which are present in all individuals, and variable genes, which are present only in some individuals. Aside from variations at the gene level, single nucleotide polymorphisms (SNPs) are also an important form of genetic variation. The advent of next-generation sequencing (NGS) coupled with the heritability of SNPs make them ideal markers for genetic analysis of human, animal, and microbial data. SNPs have also been extensively used in crop genetics for association mapping, quantitative trait loci (QTL) analysis, analysis of genetic diversity, and phylogenetic analysis. This review focuses on the use of pangenomes for SNP discovery. It highlights the advantages of using a pangenome rather than a single reference for this purpose. This review also demonstrates how extra information not captured in a single reference alone can be used to provide additional support for linking genotypic data to phenotypic data.

Lack of direct evidence for natural selection at the candidate thrifty gene locus, PPARGC1A

Background

The gene PPARGC1A, in particular the Gly482Ser variant (rs8192678), had been proposed to be subject to natural selection, particularly in recent progenitors of extant Polynesian populations. Reasons include high levels of population differentiation and increased frequencies of the derived type 2 diabetes (T2D) risk 482Ser allele, and association with body mass index (BMI) in a small Tongan population. However, no direct statistical tests for selection have been applied.

Methods

Using a range of Polynesian populations (Tongan, Māori, Samoan) we re-examined evidence for association between Gly482Ser with T2D and BMI as well as gout. Using also Asian, European, and African 1000 Genome Project samples a range of statistical tests for selection (F_ST, integrated haplotype score (iHS), cross population extended haplotype homozygosity (XP-EHH), Tajima’s D and Fay and Wu’s H) were conducted on the PPARGC1A locus.

Results

No statistically significant evidence for association between Gly482Ser and any of BMI, T2D or gout was found. Population differentiation (F_ST) was smallest between Asian and Pacific populations (New Zealand Māori ≤ 0.35, Samoan ≤ 0.20). When compared to European (New Zealand Māori ≤ 0.40, Samoan ≤ 0.25) or African populations (New Zealand Māori ≤ 0.80, Samoan ≤ 0.66) this differentiation was larger. We did not find any strong evidence for departure from neutral evolution at this locus when applying any of the other statistical tests for selection. However, using the same analytical methods, we found evidence for selection in specific populations at previously identified loci, indicating that lack of selection was the most likely explanation for the lack of evidence of selection in PPARGC1A.

Conclusion

We conclude that there is no compelling evidence for selection at this locus, and that this gene should not be considered a candidate thrifty gene locus in Pacific populations. High levels of population differentiation at this locus and the reported absence of the derived 482Ser allele in some Melanesian populations, can alternatively be explained by multiple out-of-Africa migrations by ancestral progenitors, and subsequent genetic drift during colonisation of Polynesia. Intermediate 482Ser allele frequencies in extant Western Polynesian populations could therefore be due to recent admixture with Melanesian progenitors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-016-0341-z) contains supplementary material, which is available to authorized users.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

A thrifty variant in CREBRF strongly influences body mass index in Samoans

Samoans are a unique founder population with a high prevalence of obesity, making them well suited for identifying new genetic contributors to obesity. We conducted a genome-wide association study (GWAS) in 3,072 Samoans, discovered a variant, rs12513649, strongly associated with body mass index (BMI) (P = 5.3 × 10(-14)), and replicated the association in 2,102 additional Samoans (P = 1.2 × 10(-9)). Targeted sequencing identified a strongly associated missense variant, rs373863828 (p.Arg457Gln), in CREBRF (meta P = 1.4 × 10(-20)). Although this variant is extremely rare in other populations, it is common in Samoans (frequency of 0.259), with an effect size much larger than that of any other known common BMI risk variant (1.36-1.45 kg/m(2) per copy of the risk-associated allele). In comparison to wild-type CREBRF, the Arg457Gln variant when overexpressed selectively decreased energy use and increased fat storage in an adipocyte cell model. These data, in combination with evidence of positive selection of the allele encoding p.Arg457Gln, support a 'thrifty' variant hypothesis as a factor in human obesity.

Genetic epidemiology of pharmacogenetic variants in South East Asian Malays using whole-genome sequences

Expanding the scope of pharmacogenomic research by including multiple global populations is integral to building robust evidence for its clinical translation. Deep whole-genome sequencing of diverse ethnic populations provides a unique opportunity to study rare and common pharmacogenomic markers that often vary in frequency across populations. In this study, we aim to build a diverse map of pharmacogenetic variants in South East Asian (SEA) Malay population using deep whole-genome sequences of 100 healthy SEA Malay individuals. We investigated the allelic diversity of potentially deleterious pharmacogenomic variants in SEA Malay population. Our analysis revealed 227 common and 466 rare potentially functional single nucleotide variants (SNVs) in 437 pharmacogenomic genes involved in drug metabolism, transport and target genes, including 74 novel variants. This study has created one of the most comprehensive maps of pharmacogenetic markers in any population from whole genomes and will hugely benefit pharmacogenomic investigations and drug dosage recommendations in SEA Malays.

High rates of phasing errors in highly polymorphic species with low levels of linkage disequilibrium

Short read sequencing of diploid individuals does not permit the direct inference of the sequence on each of the two homologous chromosomes. Although various phasing software packages exist, they were primarily tailored for and tested on human data, which differ from other species in factors that influence phasing, such as SNP density, amounts of linkage disequilibrium (LD) and sample sizes. Despite becoming increasingly popular for other species, the reliability of phasing in non-human data has not been evaluated to a sufficient extent. We scrutinized the phasing accuracy for Drosophila melanogaster, a species with high polymorphism levels and reduced LD relative to humans. We phased two D. melanogaster populations and compared the results to the known haplotypes. The performance increased with size of the reference panel and was highest when the reference panel and phased individuals were from the same population. Full genomic SNP data and inclusion of sequence read information also improved phasing. Despite humans and Drosophila having similar switch error rates between polymorphic sites, the distances between switch errors were much shorter in Drosophila with only fragments <300-1500 bp being correctly phased with ≥95% confidence. This suggests that the higher SNP density cannot compensate for the higher recombination rate in D. melanogaster. Furthermore, we show that populations that have gone through demographic events such as bottlenecks can be phased with higher accuracy. Our results highlight that statistically phased data are particularly error prone in species with large population sizes or populations lacking suitable reference panels.

Human-specific derived alleles of CD33 and other genes protect against postreproductive cognitive decline

The individuals of most vertebrate species die when they can no longer reproduce. Humans are a rare exception, having evolved a prolonged postreproductive lifespan. Elders contribute to cooperative offspring care, assist in foraging, and communicate important ecological and cultural knowledge, increasing the survival of younger individuals. Age-related deterioration of cognitive capacity in humans compromises these benefits and also burdens the group with socially costly members. We investigated the contribution of the immunoregulatory receptor CD33 to a uniquely human postreproductive disease, Alzheimer's dementia. Surprisingly, even though selection at advanced age is expected to be weak, a CD33 allele protective against Alzheimer's disease is derived and unique to humans and favors a functional molecular state of CD33 resembling that of the chimpanzee. Thus, derived alleles may be compensatory and restore interactions altered as a consequence of human-specific brain evolution. We found several other examples of derived alleles at other human loci that protect against age-related cognitive deterioration arising from neurodegenerative disease or cerebrovascular insufficiency. Selection by inclusive fitness may be strong enough to favor alleles protecting specifically against cognitive decline in postreproductive humans. Such selection would operate by maximizing the contributions of postreproductive individuals to the fitness of younger kin.