Long-Term Balancing Selection in LAD1 Maintains a Missense Trans-Species Polymorphism in Humans, Chimpanzees, and Bonobos

Parallel Selection in Domesticated Atlantic Salmon from Divergent Founders Including on Whole-Genome Duplication-derived Homeologous Regions

Domestication and artificial selection for desirable traits have driven significant phenotypic changes and left detectable genomic footprints in farmed animals. Since the 1960s, intensive breeding has led to the rapid domestication of Atlantic salmon (Salmo salar), with multiple independent events that make it a valuable model for studying early domestication stages and the parallel evolution of populations of different origins subjected to similar selection pressures. Some aquatic species, including Atlantic salmon, have undergone whole-genome duplication (WGD), raising the possibility that genetic redundancy resulting from WGD has contributed to adaptation in captive environments, as seen in plants. Here, we examined the genomic responses to domestication in Atlantic salmon, focusing on potential signatures of parallel selection, including those associated with WGD. Candidate genomic regions under selection were identified by comparing whole-genome sequences from aquaculture and wild populations across 2 independently domesticated lineages (Western Norway and North America) using a genome-wide scan that combined 3 statistical methods: allele frequencies (FST), site frequency (Tajima's D), and haplotype differentiation (XP-EHH). These analyses revealed shared selective sweeps on identical SNPs in major histocompatibility complex (MHC) genes across aquaculture populations. This suggests that a combination of long-term balancing selection and recent human-induced selection has shaped MHC gene evolution in domesticated salmon. Additionally, we observed selective sweeps on a small number of gene pairs in homeologous regions originating from WGD, offering insights into how historical genome duplication events may intersect with recent selection pressures in aquaculture species.

Deep genetic substructure within bonobos

Establishing the genetic and geographic structure of populations is fundamental, both to understand their evolutionary past and preserve their future. Nevertheless, the patterns of genetic population structure are unknown for most endangered species. This is the case for bonobos (Pan paniscus), which, together with chimpanzees (Pan troglodytes), are humans' closest living relatives. Chimpanzees live across equatorial Africa and are classified into four subspecies,1 with some genetic population substructure even within subspecies. Conversely, bonobos live exclusively in the Democratic Republic of Congo and are considered a homogeneous group with low genetic diversity,2 despite some population structure inferred from mtDNA. Nevertheless, mtDNA aside, their genetic structure remains unknown, hampering our understanding of the species and conservation efforts. Mapping bonobo genetic diversity in space is, however, challenging because, being endangered, only non-invasive sampling is possible for wild individuals. Here, we jointly analyze the exomes and mtDNA from 20 wild-born bonobos, the whole genomes of 10 captive bonobos, and the mtDNA of 136 wild individuals. We identify three genetically distinct bonobo groups of inferred Central, Western, and Far-Western geographic origin within the bonobo range. We estimate the split time between the central and western populations to be ∼145,000 years ago and genetic differentiation to be in the order of that of the closest chimpanzee subspecies. Furthermore, our estimated long-term Ne for Far-West (∼3,000) is among the lowest estimated for any great ape lineage. Our results highlight the need to attend to the bonobo substructure, both in terms of research and conservation.

Long-term balancing selection for pathogen resistance maintains trans-species polymorphisms in a planktonic crustacean

Balancing selection is an evolutionary process that maintains genetic polymorphisms at selected loci and strongly reduces the likelihood of allele fixation. When allelic polymorphisms that predate speciation events are maintained independently in the resulting lineages, a pattern of trans-species polymorphisms may occur. Trans-species polymorphisms have been identified for loci related to mating systems and the MHC, but they are generally rare. Trans-species polymorphisms in disease loci are believed to be a consequence of long-term host-parasite coevolution by balancing selection, the so-called Red Queen dynamics. Here we scan the genomes of three crustaceans with a divergence of over 15 million years and identify 11 genes containing identical-by-descent trans-species polymorphisms with the same polymorphisms in all three species. Four of these genes display molecular footprints of balancing selection and have a function related to immunity. Three of them are located in or close to loci involved in resistance to a virulent bacterial pathogen, Pasteuria, with which the Daphnia host is known to coevolve. This provides rare evidence of trans-species polymorphisms for loci known to be functionally relevant in interactions with a widespread and highly specific parasite. These findings support the theory that specific antagonistic coevolution is able to maintain genetic diversity over millions of years.

Temporal challenges in detecting balancing selection from population genomic data

The role of balancing selection in maintaining genetic variation remains an open question in population genetics. Recent years have seen numerous studies identifying candidate loci potentially experiencing balancing selection, most predominantly in human populations. There are however numerous alternative evolutionary processes that may leave similar patterns of variation, thereby potentially confounding inference, and the expected signatures of balancing selection additionally change in a temporal fashion. Here we use forward-in-time simulations to quantify expected statistical power to detect balancing selection using both site frequency spectrum- and linkage disequilibrium-based methods under a variety of evolutionarily realistic null models. We find that whilst site frequency spectrum-based methods have little power immediately after a balanced mutation begins segregating, power increases with time since the introduction of the balanced allele. Conversely, linkage disequilibrium-based methods have considerable power whilst the allele is young, and power dissipates rapidly as the time since introduction increases. Taken together, this suggests that site frequency spectrum-based methods are most effective at detecting long-term balancing selection (>25N generations since the introduction of the balanced allele) whilst linkage disequilibrium-based methods are effective over much shorter timescales (<1N generations), thereby leaving a large time frame over which current methods have little power to detect the action of balancing selection. Finally, we investigate the extent to which alternative evolutionary processes may mimic these patterns, and demonstrate the need for caution in attempting to distinguish the signatures of balancing selection from those of both neutral processes (e.g. population structure and admixture) as well as of alternative selective processes (e.g. partial selective sweeps).

The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

Bw4 ligand and direct T-cell receptor binding induced selection on HLA A and B alleles

Introduction

The HLA region is the hallmark of balancing selection, argued to be driven by the pressure to present a wide variety of viral epitopes. As such selection on the peptide-binding positions has been proposed to drive HLA population genetics. MHC molecules also directly binds to the T-Cell Receptor and killer cell immunoglobulin-like receptors (KIR).

Methods

We here combine the HLA allele frequencies in over six-million Hematopoietic Stem Cells (HSC) donors with a novel machine-learning-based method to predict allele frequency.

Results

We show for the first time that allele frequency can be predicted from their sequences. This prediction yields a natural measure for selection. The strongest selection is affecting KIR binding regions, followed by the peptide-binding cleft. The selection from the direct interaction with the KIR and TCR is centered on positively charged residues (mainly Arginine), and some positions in the peptide-binding cleft are not associated with the allele frequency, especially Tyrosine residues.

Discussion

These results suggest that the balancing selection for peptide presentation is combined with a positive selection for KIR and TCR binding.

Anchoring Filament Protein Ladinin-1 is an Immunosuppressive Microenvironment and Cold Tumor Correlated Prognosticator in Lung Adenocarcinoma

Anchoring filament protein ladinin-1 (LAD1) codes for an anchor filament protein in the basement membrane. Here, we have aimed to determine its potential role in LUAD. According to the comprehensive analyses conducted in this study, we studied the expression, prognostic significance, function, methylation, copy number variations, and the immune cell infiltration of LAD1 in LUAD. A higher level of LAD1 gene expression was observed in the LUAD tumor tissues compared to the normal lung tissues (p < 0.001). Furthermore, the multivariate analysis indicated that a higher LAD1 gene expression level was the independent prognostic factor. Additionally, the DNA methylation level of the LAD1 was inversely linked to its expression (p < 0.001). We noted that the patients affected due to LAD1 hypomethylation showed a very low overall survival rate compared to the patients with a higher LAD1 methylation score (p < 0.05). Moreover, the results of the immunity analysis indicated that the LAD1 expression might be inversely linked to the immune cell infiltration degree, expression of the infiltrated immune cells, and the PD-L1 levels. Lastly, we supplemented some verification to increase the rigor of the study. The results suggested that high expression of LAD1 may be related to cold tumors. Hence, this indirectly reflects that the immunotherapy effect of LUAD patients with high LAD1 expression might be worse. Based on the role played by the LAD1 in the tumor immune microenvironment, it can be considered a potential biomarker for predicting the immunotherapy response to LUAD.

LAD1 promotes malignant progression by diminishing ubiquitin-dependent degradation of vimentin in gastric cancer

BACKGROUND

Ladinin-1 (LAD1), an anchoring filament protein, has been associated with several cancer types, including cancers of the colon, lungs, and breast. However, it is still unclear how and why LAD1 causes gastric cancer (GC).

METHODS

Multiple in vitro and in vivo, functional gains and loss experiments were carried out in the current study to confirm the function of LAD1. Mass spectrometry was used to find the proteins that interact with LAD1. Immunoprecipitation analyses revealed the mechanism of LAD1 involved in promoting aggressiveness.

RESULTS

The results revealed that the LAD1 was overexpressed in GC tissues, and participants with increased LAD1 expression exhibited poorer disease-free survival (DFS) and overall survival (OS). Functionally, LAD1 promotes cellular invasion, migration, proliferation, and chemoresistance in vivo and in vitro in the subcutaneous patient-and cell-derived xenograft (PDX and CDX) tumor models. Mechanistically, LAD1 competitively bound to Vimentin, preventing it from interacting with the E3 ubiquitin ligase macrophage erythroblast attacher (MAEA), which led to a reduction in K48-linked ubiquitination of Vimentin and an increase in Vimentin protein levels in GC cells.

CONCLUSIONS

In conclusion, the current investigation indicated that LAD1 has been predicted as a possible prognostic biomarker and therapeutic target for GC due to its ability to suppress Vimentin-MAEA interaction.

ZEB1 Is Regulated by K811 Acetylation to Promote Stability, NuRD Complex Interactions, EMT, and NSCLC Metastasis

Epithelial-to-mesenchymal transition results in loss of specialized epithelial cell contacts and acquisition of mesenchymal invasive capacity. The transcription repressor zinc finger E-box-binding homeobox 1 (ZEB1) binds to E-boxes of gene promoter regions to suppress the expression of epithelial genes. ZEB1 has inconsistent molecular weights, which have been attributed to posttranslational modifications (PTM). We performed mass spectrometry and identified K811 acetylation as a novel PTM in ZEB1. To define the role of ZEB1 acetylation in regulating function, we generated ZEB1 acetyl-mimetic (K811Q) and acetyl-deficient (K811R) mutant-expressing non-small cell lung cancer cell lines (NSCLC). We demonstrate that the K811R ZEB1 (125 kDa) has a shorter protein half-life than wild-type (WT) ZEB1 and K811Q ZEB1 (∼225 kDa), suggesting that lack of ZEB1 acetylation in the lower molecular weight form affects protein stability. Further, the acetylated form of ZEB1 recruits the nucleosome remodeling and deacetylase (NuRD) complex to bind the promoter of its target genes mir200c-141 and SEMA3F. RNA-sequencing revealed that WT ZEB1 and K811Q ZEB1 downregulate the expression of epithelial genes to promote lung adenocarcinoma invasion and metastasis, whereas the K811R ZEB1 does not. Our findings establish that the K811 acetylation promotes ZEB1 protein stability, interaction with other protein complexes, and subsequent invasion/metastasis of lung adenocarcinoma via epithelial-to-mesenchymal transition.

IMPLICATIONS

The molecular mechanisms by which ZEB1 is regulated by K811 acetylation to promote protein stability, NuRD complex and promoter interactions, and function are relevant to the development of treatment strategies to prevent and treat metastasis in patients with NSCLC.

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

Prognostic roles of a novel basement membranes-related gene signature in lung adenocarcinoma

Background: The basement membranes (BMs) are involved in tumor progression, while few comprehensive analyses to date are performed on the role of BM-related gene signatures in lung adenocarcinoma (LUAD). Thus, we aimed to develop a novel prognostic model in LUAD based on BMs-related gene profiling. Methods: The LUAD BMs-related gene profiling and corresponding clinicopathological data were obtained from the basement membrane BASE, The Cancer Genome Atlas (TCGA) and gene expression omnibus (GEO) databases. The Cox regression and least absolute shrinkage and selection operator (LASSO) methods were used to construct a BMs-based risk signature. The concordance index (C-index), receiver operating characteristic (ROC), and calibration curves were generated to evaluate the nomogram. The GSE72094 dataset was used to validate prediction of the signature. The differences in functional enrichment, immune infiltration, and drug sensitivity analyses were compared based on risk score. Results: In TCGA training cohort, 10 BMs-related genes were found, (e.g., ACAN, ADAMTS15, ADAMTS8, BCAN, etc). The signal signature based on these 10 genes was categorized into high- and low-risk groups regarding survival differences (p < 0.001). Multivariable analysis showed that the signature of combined 10 BMs-related genes was an independent prognostic predictor. Such a prognostic value of BMs-based signature in validation cohort of the GSE72094 were further verified. The GEO verification, C-index, and ROC curve showed that the nomogram had accurate prediction performance. The functional analysis suggested that BMs were mainly enriched in extracellular matrix-receptor (ECM-receptor) interaction. Moreover, the BMs-based model was correlated with immune checkpoint. Conclusion: This study identified BMs-based risk signature genes and demonstrated their ability to predict prognosis and guide personalized treatment of patients with LUAD.

The Tumorigenic Effect of the High Expression of Ladinin-1 in Lung Adenocarcinoma and Its Potential as a Therapeutic Target

The oncogenic role of Ladinin-1 (LAD1), an anchoring filament protein, is largely unknown. In this study, we conducted a series of studies on the oncogenic role of LAD1 in lung adenocarcinoma (LUAD). Firstly, we analyzed the aberrant expression of LAD1 in LUAD and its correlation with patient survival, tumor immune infiltration, and the activation of cancer signaling pathways. Furthermore, the relationship between LAD1 expression and K-Ras and EGF signaling activation, tumor cell proliferation, migration, and colony formation was studied by gene knockout/knockout methods. We found that LAD1 was frequently overexpressed in LUAD, and high LAD1 expression predicts a poor prognosis. LAD1 exhibits promoter hypomethylation in LUAD, which may contribute to its mRNA upregulation. Single-sample gene set enrichment analysis (ssGSEA) showed that acquired immunity was negatively correlated with LAD1 expression, which was verified by the downregulated GO terms of "Immunoglobulin receptor binding" and "Immunoglobulin complex circulating" in the LAD1 high-expression group through Gene Set Variation Analysis (GSVA). Notably, the Ras-dependent signature was the most activated signaling in the LAD1 high-expression group, and the phosphorylation of downstream effectors, such as ERK and c-jun, was strongly inhibited by LAD1 deficiency. Moreover, we demonstrated that LAD1 depletion significantly inhibited the proliferation, migration, and cell-cycle progression of LUAD cells and promoted sensitivity to Gefitinib, K-Ras inhibitor, and paclitaxel treatments. We also confirmed that LAD1 deficiency remarkably retarded tumor growth in the xenograft model. Conclusively, LAD1 is a critical prognostic biomarker for LUAD and has potential as an intervention target.

Ladinin 1 Shortens Survival via Promoting Proliferation and Enhancing Invasiveness in Lung Adenocarcinoma

Lung cancer is one of the deadliest cancers worldwide, including in Taiwan. The poor prognosis of the advanced lung cancer lies in delayed diagnosis and non-druggable targets. It is worth paying more attention to these ongoing issues. Public databases and an in-house cohort were used for validation. The KM plotter was utilized to discover the clinical significance. GSEA and GSVA were adopted for a functional pathway survey. Molecular biological methods, including proliferation, migration, and the EMT methods, were used for verification. Based on public databases, the increased expression of Ladinin 1 (LAD1) was presented in tumor and metastatic sites. Furthermore, an in-house cohort revealed a higher intensity of LAD1 in tumor rather than in normal parts. The greater the expression of LAD1 was, the shorter the duration of lung adenocarcinoma (LUAD) patient survival. Moreover, the association of B3GNT3 with LAD1 affected the survival of LUAD patients. Functional analyses using GSEA and GSVA revealed the associations with survival, migration, invasion, and EMT. Biologic functions supported the roles of LAD1 in proliferation via the cell cycle and migration in EMT. This study reveals that LAD1 plays a major role in regulating proliferation and migration in lung cancer and impacts survival in LUAD. It is worth investing in further studies and in the development of drugs targeting LAD1.

Landscape of circular RNAs in different types of lung cancer and an emerging role in therapeutic resistance (Review)

Lung cancer is one of the most common malignant tumor types and the leading cause of cancer‑associated death worldwide. Different types of lung cancer exhibit differences in terms of pathophysiology and pathogenesis, and also treatment and prognosis. Accumulating evidence has indicated that circular RNAs (circRNAs) are abnormally expressed among different types of lung cancer and confer important biological functions in progression and prognosis. However, studies comparing different circRNAs in lung cancer subtypes are scarce. Furthermore, circRNAs have an important role in drug resistance and are related to clinicopathological features in lung cancer. Summaries of the association of circRNAs with drug resistance are also scarce in the literature. The present study outlined the biological functions of circRNAs and focused on discriminating differential circRNA patterns and mechanisms in three different types of lung cancer. The emerging roles of circRNAs in the resistance to chemotherapy, targeted therapy, radiotherapy and immunotherapy were also highlighted. Understanding these aspects of circRNAs sheds light on novel physiological and pathophysiological processes of lung cancer and suggests the application of circRNAs as biomarkers for diagnosis and prognosis, as well as therapeutic resistance.

DNA Methylome and Transcriptome Study of Triterpenoid CDDO in TPA-Mediated Skin Carcinogenesis Model

Overexposure to ultraviolet radiation and environmental carcinogens drive skin cancer development through redox imbalance and gene mutation. Antioxidants such as triterpenoids have exhibited anti-oxidative and anti-inflammatory potentials to alleviate skin carcinogenesis. This study investigated the methylome and transcriptome altered by tumor promoter 12-O-tetradecanoylphorbol-13-acetate (TPA) or TPA with 2-cyano 2,3-dioxoolean-1,9-dien-28-oic acid (CDDO). The results show that CDDO blocks TPA-induced transformation dose dependently. Several differential expressed genes (DEGs) involved in skin cell transformation, while counteracted by CDDO, were revealed by differential expression analysis including Lyl1, Lad1, and Dennd2d. In CpG methylomic profiles, the differentially methylated regions (DMRs) in the promoter region altered by TPA while showing the opposite methylation status in the CDDO treatment group were identified. The correlation between DNA methylation and RNA expression has been established and DMRs showing inverse correlation were further studied as potential therapeutic targets. From the CpG methylome and transcriptome results, CDDO significantly restored gene expression of NAD(P)H:quinone oxidoreductase 1 (Nqo1) inhibited by TPA by decreasing their promoter CpG methylation. Ingenuity Pathways Analysis (IPA) shows that CDDO neutralized the effect of TPA through modulating cell cycles, cell migration, and inflammatory and immune response regulatory pathways. Notably, Tumor Necrosis Factor Receptor 2 (TNFR2) signaling was significantly downregulated by CDDO potentially contributing to prevention of TPA-induced cell transformation. Overall, incorporating the transcriptome, CpG methylome, and signaling pathway network, we reveal potential therapeutic targets and pathways by which CDDO could reverse TPA-induced carcinogenesis. The results could be useful for future human study and targets development for skin cancer.

The roles of balancing selection and recombination in the evolution of rattlesnake venom

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.

The immunogenetic impact of European colonization in the Americas

The introduction of pathogens originating from Eurasia into the Americas during early European contact has been associated with high mortality rates among Indigenous peoples, likely contributing to their historical and precipitous population decline. However, the biological impacts of imported infectious diseases and resulting epidemics, especially in terms of pathogenic effects on the Indigenous immunity, remain poorly understood and highly contentious to this day. Here, we examine multidisciplinary evidence underpinning colonization-related immune genetic change, providing contextualization from anthropological studies, paleomicrobiological evidence of contrasting host-pathogen coevolutionary histories, and the timings of disease emergence. We further summarize current studies examining genetic signals reflecting post-contact Indigenous population bottlenecks, admixture with European and other populations, and the putative effects of natural selection, with a focus on ancient DNA studies and immunity-related findings. Considering current genetic evidence, together with a population genetics theoretical approach, we show that post-contact Indigenous immune adaptation, possibly influenced by selection exerted by introduced pathogens, is highly complex and likely to be affected by multifactorial causes. Disentangling putative adaptive signals from those of genetic drift thus remains a significant challenge, highlighting the need for the implementation of population genetic approaches that model the short time spans and complex demographic histories under consideration. This review adds to current understandings of post-contact immunity evolution in Indigenous peoples of America, with important implications for bettering our understanding of human adaptation in the face of emerging infectious diseases.

Diverse functions associate with non-coding polymorphisms shared between humans and chimpanzees

Background

Long-term balancing selection (LTBS) can maintain allelic variation at a locus over millions of years and through speciation events. Variants shared between species in the state of identity-by-descent, hereafter “trans-species polymorphisms”, can result from LTBS, often due to host–pathogen interactions. For instance, the major histocompatibility complex (MHC) locus contains TSPs present across primates. Several hundred candidate LTBS regions have been identified in humans and chimpanzees; however, because many are in non-protein-coding regions of the genome, the functions and potential adaptive roles for most remain unknown.

Results

We integrated diverse genomic annotations to explore the functions of 60 previously identified regions with multiple shared polymorphisms (SPs) between humans and chimpanzees, including 19 with strong evidence of LTBS. We analyzed genome-wide functional assays, expression quantitative trait loci (eQTL), genome-wide association studies (GWAS), and phenome-wide association studies (PheWAS) for all the regions. We identify functional annotations for 59 regions, including 58 with evidence of gene regulatory function from GTEx or functional genomics data and 19 with evidence of trait association from GWAS or PheWAS. As expected, the SPs associate in humans with many immune system phenotypes, including response to pathogens, but we also find associations with a range of other phenotypes, including body size, alcohol intake, cognitive performance, risk-taking behavior, and urate levels.

Conclusions

The diversity of traits associated with non-coding regions with multiple SPs support previous hypotheses that functions beyond the immune system are likely subject to LTBS. Furthermore, several of these trait associations provide support and candidate genetic loci for previous hypothesis about behavioral diversity in human and chimpanzee populations, such as the importance of variation in risk sensitivity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02020-x.

Similarity-Based Analysis of Allele Frequency Distribution among Multiple Populations Identifies Adaptive Genomic Structural Variants

Structural variants have a considerable impact on human genomic diversity. However, their evolutionary history remains mostly unexplored. Here, we developed a new method to identify potentially adaptive structural variants based on a similarity-based analysis that incorporates genotype frequency data from 26 populations simultaneously. Using this method, we analyzed 57,629 structural variants and identified 576 structural variants that show unusual population differentiation. Of these putatively adaptive structural variants, we further showed that 24 variants are multiallelic and overlap with coding sequences, and 20 variants are significantly associated with GWAS traits. Closer inspection of the haplotypic variation associated with these putatively adaptive and functional structural variants reveals deviations from neutral expectations due to: 1) population differentiation of rapidly evolving multiallelic variants, 2) incomplete sweeps, and 3) recent population-specific negative selection. Overall, our study provides new methodological insights, documents hundreds of putatively adaptive variants, and introduces evolutionary models that may better explain the complex evolution of structural variants.

BalLeRMix +: mixture model approaches for robust joint identification of both positive selection and long-term balancing selection

SUMMARY

The growing availability of genomewide polymorphism data has fueled interest in detecting diverse selective processes affecting population diversity. However, no model-based approaches exist to jointly detect and distinguish the two complementary processes of balancing and positive selection. We extend the BalLeRMix  B-statistic framework described in Cheng and DeGiorgio (2020) for detecting balancing selection and present BalLeRMix+, which implements five B statistic extensions based on mixture models to robustly identify both types of selection. BalLeRMix+ is implemented in Python and computes the composite likelihood ratios and associated model parameters for each genomic test position.

AVAILABILITY AND IMPLEMENTATION

BalLeRMix+ is freely available at https://github.com/bioXiaoheng/BallerMixPlus.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

High expression of ladinin-1 (LAD1) predicts adverse outcomes: a new candidate docetaxel resistance gene for prostatic cancer (PCa)

Docetaxel resistance is one of the major obstacles that undermine the treatment outcome of PCa. Exploring molecular mechanisms associated with docetaxel resistance could provide insights into the formulation of novel strategies enhancing the efficacy of PCa treatment. Ladinin-1 (LAD1) is an anchoring filament protein in basement membranes, which contributes to the association of the epithelial cells with the underlying mesenchyme. LAD1 has been implicated in the progression of different cancers. However, its role in PCa remains to be investigated. In the present study, we found that LAD1 was highly expressed in docetaxel-resistant PCa cells, while its expression was significantly suppressed in tumor samples after docetaxel treatment. Moreover, the expression level of LAD1 in PCa tissues was significantly higher than that of normal tissue, and high expression level of LAD1 was significantly associated with adverse outcomes of PCa patients. Finally, high expression of LAD1 in PCa tissue was also correlated with the expression level of genes involving in tumor cell proliferation and invasive behaviors. Collectively, our data suggest that LAD1 may serve as a potential prognostic factor in PCa patients.

Balancing selection maintains ancient polymorphisms at conserved enhancers for the olfactory receptor genes of a Chinese marine fish

The study of balancing selection, as a selective force maintaining adaptive genetic variation in gene pools longer than expected by drift, is currently experiencing renewed interest due to the increased availability of new data, methods of analysis, and case studies. In this investigation, evidence of balancing selection operating on conserved enhancers of the olfactory receptor (OR) genes is presented for the Chinese sleeper (Bostrychus sinensis), a coastal marine fish that is emerging as a model species for evolutionary studies in the Northwest Pacific marginal seas. Coupled with tests for Gene Ontology enrichment and transcription factor binding, population genomic data allow for the identification of an OR cluster in the sleeper with a downstream flanking region containing three enhancers that are conserved with human and other fish species. Phylogenetic and population genetic analyses indicate that the enhancers are under balancing selection as evidenced by their translineage polymorphisms, excess common alleles, and increased within-group diversities. Age comparisons between the translineage polymorphisms and most recent common ancestors of neutral genealogies substantiate that the former are old, and thus, due to ancient balancing selection. The survival and reproduction of vertebrates depend on their sense of smell, and thereby, on their ORs. In addition to locus duplication and allelic variation of structural genes, this study highlights a third mechanism by which receptor diversity can be achieved for detecting and responding to the huge variety of environmental odorants (i.e., by balancing selection acting on OR gene expression through their enhancer variability).

Cooperation, involuntary defection, and Polymorphic Evolutionary Games

The subject of involuntary defection in the context of cycles of interaction approach to direct reciprocal cooperation was introduced some time ago (J. Theor. Biol., 242: 873-879). Current work is motivated by the subsequent accumulation of empirical evidence and the advances in the methodology of evolutionary games. In recent decades it become clear that individuals in many animal species vary consistently in their behavioral responses to specific challenges-animal personality. Moreover, these differences have a hereditary component. Finally, investigations into the effects of neuropeptides on behavior suggest that the variations in animal personalities involve polymorphisms based on non-Mendelian heritability within the neuropeptide signaling systems. The last observation suggests that animal personalities can be productively analyzed via Polymorphic Evolutionary Games, which allow us to add genetic model(s) to standard (phenotypic) evolutionary games. Such an analysis of reciprocal cooperation is the subject of the current paper. The results indicate that there is a marked difference between models that assume Mendelian vs. non-Mendelian inheritance of the pertinent traits. Monomorphic and polymorphic-Mendelian models predict pure-strategy (single phenotype) ESS, whereas the non-Mendelian genetic model predicts a mixed strategy ESS exhibiting all three phenotypes.

Distinguishing between recent balancing selection and incomplete sweep using deep neural networks

Balancing selection is an important adaptive mechanism underpinning a wide range of phenotypes. Despite its relevance, the detection of recent balancing selection from genomic data is challenging as its signatures are qualitatively similar to those left by ongoing positive selection. In this study, we developed and implemented two deep neural networks and tested their performance to predict loci under recent selection, either due to balancing selection or incomplete sweep, from population genomic data. Specifically, we generated forward-in-time simulations to train and test an artificial neural network (ANN) and a convolutional neural network (CNN). ANN received as input multiple summary statistics calculated on the locus of interest, while CNN was applied directly on the matrix of haplotypes. We found that both architectures have high accuracy to identify loci under recent selection. CNN generally outperformed ANN to distinguish between signals of balancing selection and incomplete sweep and was less affected by incorrect training data. We deployed both trained networks on neutral genomic regions in European populations and demonstrated a lower false-positive rate for CNN than ANN. We finally deployed CNN within the MEFV gene region and identified several common variants predicted to be under incomplete sweep in a European population. Notably, two of these variants are functional changes and could modulate susceptibility to familial Mediterranean fever, possibly as a consequence of past adaptation to pathogens. In conclusion, deep neural networks were able to characterize signals of selection on intermediate frequency variants, an analysis currently inaccessible by commonly used strategies.

The inflated significance of neutral genetic diversity in conservation genetics

The current rate of species extinction is rapidly approaching unprecedented highs, and life on Earth presently faces a sixth mass extinction event driven by anthropogenic activity, climate change, and ecological collapse. The field of conservation genetics aims at preserving species by using their levels of genetic diversity, usually measured as neutral genome-wide diversity, as a barometer for evaluating population health and extinction risk. A fundamental assumption is that higher levels of genetic diversity lead to an increase in fitness and long-term survival of a species. Here, we argue against the perceived importance of neutral genetic diversity for the conservation of wild populations and species. We demonstrate that no simple general relationship exists between neutral genetic diversity and the risk of species extinction. Instead, a better understanding of the properties of functional genetic diversity, demographic history, and ecological relationships is necessary for developing and implementing effective conservation genetic strategies.

The Pan social brain: An evolutionary history of neurochemical receptor genes and their potential impact on sociocognitive differences

Humans have unique cognitive capacities that, compared with apes, are not only simply expressed as a higher level of general intelligence, but also as a quantitative difference in sociocognitive skills. Humans' closest living relatives, bonobos (Pan paniscus), and chimpanzees (Pan troglodytes), show key between-species differences in social cognition despite their close phylogenetic relatedness, with bonobos arguably showing greater similarities to humans. To better understand the evolution of these traits, we investigate the neurochemical mechanisms underlying sociocognitive skills by focusing on variation in genes encoding proteins with well-documented roles in mammalian social cognition: the receptors for vasopressin (AVPR1A), oxytocin (OXTR), serotonin (HTR1A), and dopamine (DRD2). Although these genes have been well studied in humans, little is known about variation in these genes that may underlie differences in social behavior and cognition in apes. We comparatively analyzed sequence data for 33 bonobos and 57 chimpanzees, together with orthologous sequence data for other apes. In all four genes, we describe genetic variants that alter the amino acid sequence of the respective receptors, raising the possibility that ligand binding or signal transduction may be impacted. Overall, bonobos show 57% more fixed substitutions than chimpanzees compared with the ancestral Pan lineage. Chimpanzees, show 31% more polymorphic coding variation, in line with their larger historical effective population size estimates and current wider distribution. An extensive literature review comparing allelic changes in Pan with known human behavioral variants revealed evidence of homologous evolution in bonobos and humans (OXTR rs4686301(T) and rs237897(A)), while humans and chimpanzees shared OXTR rs2228485(A), DRD2 rs6277(A), and DRD2 rs11214613(A) to the exclusion of bonobos. Our results offer the first in-depth comparison of neurochemical receptor gene variation in Pan and put forward new variants for future behavior-genotype association studies in apes, which can increase our understanding of the evolution of social cognition in modern humans.

LAD1 expression is associated with the metastatic potential of colorectal cancer cells

BACKGROUND

Anchoring filament protein ladinin-1 (LAD1) was related to the aggressive progression of breast, lung, laryngeal and thyroid cancers. However, the association of LAD1 with colorectal cancer remained unknown. Here, to determine the relationship of LAD1 with colorectal cancer progression, we explored the effect of LAD1 loss on the malignant features of colorectal cancer cells.

METHODS

We constructed LAD1-depleted cell lines and examined the effect of LAD1 deficiency on the phenotypic and molecular features of colorectal cancer cells in vitro. The function of LAD1 in metastasis in vivo was examined by establishing a spleen-to-liver metastasis mouse model. LAD1 protein expression in colorectal cancer patient specimens was assessed by immunohistochemistry of tumor microarrays.

RESULTS

We found that LAD1 was abundant in most colorectal cancer cells. In addition, high expression of LAD1 significantly correlated with poor patient outcome. LAD1 depletion inhibited the migration and invasion of two different colorectal cancer cell lines, SW620 and Caco-2, without affecting their proliferation. In addition, LAD1 loss led to defects in liver metastasis of SW620 cells in the mouse model. Immunohistochemistry of colorectal cancer tissues revealed LAD1 enrichment in metastatic tissues compared to that in primary tumor and normal tissues.

CONCLUSION

These results suggest that LAD1 expression is associated with the metastatic progression of colorectal cancer by promoting the migration and invasion of cancer cells.

Genomic insights on the contribution of balancing selection and local adaptation to the long-term survival of a widespread living fossil tree, Cercidiphyllum japonicum

'Living fossils' are testimonies of long-term sustained ecological success, but how demographic history and natural selection contributed to their survival, resilience, and persistence in the face of Quaternary climate fluctuations remains unclear. To better understand the interplay between demographic history and selection in shaping genomic diversity and evolution of such organisms, we assembled the whole genome of Cercidiphyllum japonicum, a widespread East Asian Tertiary relict tree, and resequenced 99 individuals of C. japonicum and its sister species, Cercidiphyllum magnificum (Central Japan). We dated this speciation event to the mid-Miocene, and the intraspecific lineage divergence of C. japonicum (China vs Japan) to the Early Pliocene. Throughout climatic upheavals of the late Tertiary/Quaternary, population bottlenecks greatly reduced the genetic diversity of C. japonicum. However, this polymorphism loss was likely counteracted by, first, long-term balancing selection at multiple chromosomal and heterozygous gene regions, potentially reflecting overdominance, and, second, selective sweeps at stress response and growth-related genes likely involved in local adaptation. Our findings contribute to a better understanding of how living fossils have survived climatic upheaval and maintained an extensive geographic range; that is, both types of selection could be major factors contributing to the species' survival, resilience, and persistence.

Flexible Mixture Model Approaches That Accommodate Footprint Size Variability for Robust Detection of Balancing Selection

Long-term balancing selection typically leaves narrow footprints of increased genetic diversity, and therefore most detection approaches only achieve optimal performances when sufficiently small genomic regions (i.e., windows) are examined. Such methods are sensitive to window sizes and suffer substantial losses in power when windows are large. Here, we employ mixture models to construct a set of five composite likelihood ratio test statistics, which we collectively term B statistics. These statistics are agnostic to window sizes and can operate on diverse forms of input data. Through simulations, we show that they exhibit comparable power to the best-performing current methods, and retain substantially high power regardless of window sizes. They also display considerable robustness to high mutation rates and uneven recombination landscapes, as well as an array of other common confounding scenarios. Moreover, we applied a specific version of the B statistics, termed B2, to a human population-genomic data set and recovered many top candidates from prior studies, including the then-uncharacterized STPG2 and CCDC169-SOHLH2, both of which are related to gamete functions. We further applied B2 on a bonobo population-genomic data set. In addition to the MHC-DQ genes, we uncovered several novel candidate genes, such as KLRD1, involved in viral defense, and SCN9A, associated with pain perception. Finally, we show that our methods can be extended to account for multiallelic balancing selection and integrated the set of statistics into open-source software named BalLeRMix for future applications by the scientific community.

Binding affinities of 438 HLA proteins to complete proteomes of seven pandemic viruses and distributions of strongest and weakest HLA peptide binders in populations worldwide

We report detailed peptide‐binding affinities between 438 HLA Class I and Class II proteins and complete proteomes of seven pandemic human viruses, including coronaviruses, influenza viruses and HIV‐1. We contrast these affinities with HLA allele frequencies across hundreds of human populations worldwide. Statistical modelling shows that peptide‐binding affinities classified into four distinct categories depend on the HLA locus but that the type of virus is only a weak predictor, except in the case of HIV‐1. Among the strong HLA binders (IC₅₀ ≤ 50), we uncovered 16 alleles (the top ones being A*02:02, B*15:03 and DRB1*01:02) binding more than 1% of peptides derived from all viruses, 9 (top ones including HLA‐A*68:01, B*15:25, C*03:02 and DRB1*07:01) binding all viruses except HIV‐1, and 15 (top ones A*02:01 and C*14:02) only binding coronaviruses. The frequencies of strongest and weakest HLA peptide binders differ significantly among populations from different geographic regions. In particular, Indigenous peoples of America show both higher frequencies of strongest and lower frequencies of weakest HLA binders. As many HLA proteins are found to be strong binders of peptides derived from distinct viral families, and are hence promiscuous (or generalist), we discuss this result in relation to possible signatures of natural selection on HLA promiscuous alleles due to past pathogenic infections. Our findings are highly relevant for both evolutionary genetics and the development of vaccine therapies. However they should not lead to forget that individual resistance and vulnerability to diseases go beyond the sole HLA allelic affinity and depend on multiple, complex and often unknown biological, environmental and other variables.

Evolutionary and population (epi)genetics of immunity to infection

Immune response is one of the functions that have been more strongly targeted by natural selection during human evolution. The evolutionary genetic dissection of the immune system has greatly helped to distinguish genes and functions that are essential, redundant or advantageous for human survival. It is also becoming increasingly clear that admixture between early Eurasians with now-extinct hominins such as Neanderthals or Denisovans, or admixture between modern human populations, can be beneficial for human adaptation to pathogen pressures. In this review, we discuss how the integration of population genetics with functional genomics in diverse human populations can inform about the changes in immune functions related to major lifestyle transitions (e.g., from hunting and gathering to farming), the action of natural selection to the evolution of the immune system, and the history of past epidemics. We also highlight the need of expanding the characterization of the immune system to a larger array of human populations-particularly neglected human groups historically exposed to different pathogen pressures-to fully capture the relative contribution of genetic, epigenetic, and environmental factors to immune response variation in humans.

A review of research in primate sanctuaries

While non-human primate studies have long been conducted in laboratories, and more recently at zoological parks, sanctuaries are increasingly considered a viable setting for research. Accredited sanctuaries in non-range countries house thousands of primates formerly used as subjects of medical research, trained performers or personal pets. In range countries, however, sanctuaries typically house orphaned primates confiscated from illegal poaching and the bushmeat and pet trafficking trades. Although the primary mission of these sanctuaries is to rescue and rehabilitate residents, many of these organizations are increasingly willing to participate in non-invasive research. Notably, from a scientific standpoint, most sanctuaries provide potential advantages over traditional settings, such as large, naturalistic physical and social environments which may result in more relevant models of primates' free-ranging wild counterparts than other captive settings. As a result, an impressive scope of research in the fields of primate behaviour, cognition, veterinary science, genetics and physiology have been studied in sanctuaries. In this review, we examine the range and form of research that has been conducted at accredited sanctuaries around the world. We also describe the potential challenges of sanctuary-based work and the considerations that external researchers may face when deciding to collaborate with primate sanctuaries on their research projects.

Genetic susceptibility to severe childhood asthma and rhinovirus-C maintained by balancing selection in humans for 150 000 years

Selective pressures imposed by pathogens have varied among human populations throughout their evolution, leading to marked inter-population differences at some genes mediating susceptibility to infectious and immune-related diseases. Here, we investigated the evolutionary history of a common polymorphism resulting in a Y529 versus C529 change in the cadherin related family member 3 (CDHR3) receptor which underlies variable susceptibility to rhinovirus-C infection and is associated with severe childhood asthma. The protective variant is the derived allele and is found at high frequency worldwide (69-95%). We detected genome-wide significant signatures of natural selection consistent with a rapid increase of the haplotypes carrying the allele, suggesting that non-neutral processes have acted on this locus across all human populations. However, the allele has not fixed in any population despite multiple lines of evidence suggesting that the mutation predates human migrations out of Africa. Using an approximate Bayesian computation method, we estimate the age of the mutation while explicitly accounting for past demography and positive or frequency-dependent balancing selection. Our analyses indicate a single emergence of the mutation in anatomically modern humans ~150 000 years ago and indicate that balancing selection has maintained the beneficial allele at high equilibrium frequencies worldwide. Apart from the well-known cases of the MHC and ABO genes, this study provides the first evidence that negative frequency-dependent selection plausibly acted on a human disease susceptibility locus, a form of balancing selection compatible with typical transmission dynamics of communicable respiratory viruses that might exploit CDHR3.

Polymorphic Evolutionary Games and Non-Mendelian Genetics

Classical evolutionary game theory (EGT) focuses on competition among phenotypes while assuming asexual transmission of these phenotypes to the next generation. However, phenotypic selection and sexual recombination are not necessarily mutually reinforcing in populations with sexual reproduction. In particular, it has been long known that some of the evolutionarily stable strategies derived by EGT methods cannot be achieved by sexually reproducing, real-world, populations. Thus, the recently formulated polymorphic evolutionary game theory (PEGT), which adds underlying genetics and sexual reproduction to evolutionary games, has the potential to revolutionize game theoretical modeling of coevolutionary processes. To illustrate the advantages of PEGT over classical EGT, I analyze two of the best known EGT models: Hawk/Retaliator/Dove and Defector/Tit-for-Tat/Altruist, by PEGT methods. I show that if one admits non-Mendelian genetics-common in heritable behavior, both of these games exhibit the properties of moderated aggression and conditional cooperation as components of population-level polymorphisms.

CREBBP and WDR 24 Identified as Candidate Genes for Quantitative Variation in Red-Brown Plumage Colouration in the Chicken

Plumage colouration in birds is important for a plethora of reasons, ranging from camouflage, sexual signalling, and species recognition. The genes underlying colour variation have been vital in understanding how genes can affect a phenotype. Multiple genes have been identified that affect plumage variation, but research has principally focused on major-effect genes (such as those causing albinism, barring, and the like), rather than the smaller effect modifier loci that more subtly influence colour. By utilising a domestic × wild advanced intercross with a combination of classical QTL mapping of red colouration as a quantitative trait and a targeted genetical genomics approach, we have identified five separate candidate genes (CREBBP, WDR24, ARL8A, PHLDA3, LAD1) that putatively influence quantitative variation in red-brown colouration in chickens. By treating colour as a quantitative rather than qualitative trait, we have identified both QTL and genes of small effect. Such small effect loci are potentially far more prevalent in wild populations, and can therefore potentially be highly relevant to colour evolution.

Genetic diversity of CHC22 clathrin impacts its function in glucose metabolism

CHC22 clathrin plays a key role in intracellular membrane traffic of the insulin-responsive glucose transporter GLUT4 in humans. We performed population genetic and phylogenetic analyses of the CHC22-encoding CLTCL1 gene, revealing independent gene loss in at least two vertebrate lineages, after arising from gene duplication. All vertebrates retained the paralogous CLTC gene encoding CHC17 clathrin, which mediates endocytosis. For vertebrates retaining CLTCL1, strong evidence for purifying selection supports CHC22 functionality. All human populations maintained two high frequency CLTCL1 allelic variants, encoding either methionine or valine at position 1316. Functional studies indicated that CHC22-V1316, which is more frequent in farming populations than in hunter-gatherers, has different cellular dynamics than M1316-CHC22 and is less effective at controlling GLUT4 membrane traffic, altering its insulin-regulated response. These analyses suggest that ancestral human dietary change influenced selection of allotypes that affect CHC22's role in metabolism and have potential to differentially influence the human insulin response.

Influence of Recombination and GC-biased Gene Conversion on the Adaptive and Nonadaptive Substitution Rate in Mammals versus Birds

Recombination is expected to affect functional sequence evolution in several ways. On the one hand, recombination is thought to improve the efficiency of multilocus selection by dissipating linkage disequilibrium. On the other hand, natural selection can be counteracted by recombination-associated transmission distorters such as GC-biased gene conversion (gBGC), which tends to promote G and C alleles irrespective of their fitness effect in high-recombining regions. It has been suggested that gBGC might impact coding sequence evolution in vertebrates, and particularly the ratio of nonsynonymous to synonymous substitution rates (dN/dS). However, distinctive gBGC patterns have been reported in mammals and birds, maybe reflecting the documented contrasts in evolutionary dynamics of recombination rate between these two taxa. Here, we explore how recombination and gBGC affect coding sequence evolution in mammals and birds by analyzing proteome-wide data in six species of Galloanserae (fowls) and six species of catarrhine primates. We estimated the dN/dS ratio and rates of adaptive and nonadaptive evolution in bins of genes of increasing recombination rate, separately analyzing AT → GC, GC → AT, and G ↔ C/A ↔ T mutations. We show that in both taxa, recombination and gBGC entail a decrease in dN/dS. Our analysis indicates that recombination enhances the efficiency of purifying selection by lowering Hill-Robertson effects, whereas gBGC leads to an overestimation of the adaptive rate of AT → GC mutations. Finally, we report a mutagenic effect of recombination, which is independent of gBGC.

Long-term balancing selection drives evolution of immunity genes in Capsella

Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.

Detection of Shared Balancing Selection in the Absence of Trans-Species Polymorphism

Trans-species polymorphism has been widely used as a key sign of long-term balancing selection across multiple species. However, such sites are often rare in the genome and could result from mutational processes or technical artifacts. Few methods are yet available to specifically detect footprints of trans-species balancing selection without using trans-species polymorphic sites. In this study, we develop summary- and model-based approaches that are each specifically tailored to uncover regions of long-term balancing selection shared by a set of species by using genomic patterns of intraspecific polymorphism and interspecific fixed differences. We demonstrate that our trans-species statistics have substantially higher power than single-species approaches to detect footprints of trans-species balancing selection, and are robust to those that do not affect all tested species. We further apply our model-based methods to human and chimpanzee whole-genome sequencing data. In addition to the previously established major histocompatibility complex and malaria resistance-associated FREM3/GYPE regions, we also find outstanding genomic regions involved in barrier integrity and innate immunity, such as the GRIK1/CLDN17 intergenic region, and the SLC35F1 and ABCA13 genes. Our findings not only echo the significance of pathogen defense but also reveal novel candidates in maintaining balanced polymorphisms across human and chimpanzee lineages. Finally, we show that these trans-species statistics can be applied to and work well for an arbitrary number of species, and integrate them into open-source software packages for ease of use by the scientific community.

Signatures of Long-Term Balancing Selection in Human Genomes

Balancing selection maintains advantageous diversity in populations through various mechanisms. Although extensively explored from a theoretical perspective, an empirical understanding of its prevalence and targets lags behind our knowledge of positive selection. Here, we describe the Non-central Deviation (NCD), a simple yet powerful statistic to detect long-term balancing selection (LTBS) that quantifies how close frequencies are to expectations under LTBS, and provides the basis for a neutrality test. NCD can be applied to a single locus or genomic data, and can be implemented considering only polymorphisms (NCD1) or also considering fixed differences with respect to an outgroup (NCD2) species. Incorporating fixed differences improves power, and NCD2 has higher power to detect LTBS in humans under different frequencies of the balanced allele(s) than other available methods. Applied to genome-wide data from African and European human populations, in both cases using chimpanzee as an outgroup, NCD2 shows that, albeit not prevalent, LTBS affects a sizable portion of the genome: ∼0.6% of analyzed genomic windows and 0.8% of analyzed positions. Significant windows (P < 0.0001) contain 1.6% of SNPs in the genome, which disproportionally fall within exons and change protein sequence, but are not enriched in putatively regulatory sites. These windows overlap ∼8% of the protein-coding genes, and these have larger number of transcripts than expected by chance even after controlling for gene length. Our catalog includes known targets of LTBS but a majority of them (90%) are novel. As expected, immune-related genes are among those with the strongest signatures, although most candidates are involved in other biological functions, suggesting that LTBS potentially influences diverse human phenotypes.

Gene buddies: linked balanced polymorphisms reinforce each other even in the absence of epistasis

The fates of genetic polymorphisms maintained by balancing selection depend on evolutionary dynamics at linked sites. While coevolution across linked, epigenetically-interacting loci has been extensively explored, such supergenes may be relatively rare. However, genes harboring adaptive variation can occur in close physical proximity while generating independent effects on fitness. Here, I present a model in which two linked loci without epistasis are both under balancing selection for unrelated reasons. Using forward-time simulations, I show that recombination rate strongly influences the retention of adaptive polymorphism, especially for intermediate selection coefficients. A locus is more likely to retain adaptive variation if it is closely linked to another locus under balancing selection, even if the two loci have no interaction. Thus, two linked polymorphisms can both be retained indefinitely even when they would both be lost to drift if unlinked. While these results may be intuitive, they have important implications for genetic architecture: clusters of mutually reinforcing genes may underlie phenotypic variation in natural populations, and such genes cannot be assumed to be functionally associated. Future studies that measure selection coefficients and recombination rates among closely linked genes will be fruitful for characterizing the extent of this phenomenon.

Detecting Long-Term Balancing Selection Using Allele Frequency Correlation

Balancing selection occurs when multiple alleles are maintained in a population, which can result in their preservation over long evolutionary time periods. A characteristic signature of this long-term balancing selection is an excess number of intermediate frequency polymorphisms near the balanced variant. However, the expected distribution of allele frequencies at these loci has not been extensively detailed, and therefore existing summary statistic methods do not explicitly take it into account. Using simulations, we show that new mutations which arise in close proximity to a site targeted by balancing selection accumulate at frequencies nearly identical to that of the balanced allele. In order to scan the genome for balancing selection, we propose a new summary statistic, β, which detects these clusters of alleles at similar frequencies. Simulation studies show that compared with existing summary statistics, our measure has improved power to detect balancing selection, and is reasonably powered in non-equilibrium demographic models and under a range of recombination and mutation rates. We compute β on 1000 Genomes Project data to identify loci potentially subjected to long-term balancing selection in humans. We report two balanced haplotypes-localized to the genes WFS1 and CADM2-that are strongly linked to association signals for complex traits. Our approach is computationally efficient and applicable to species that lack appropriate outgroup sequences, allowing for well-powered analysis of selection in the wide variety of species for which population data are rapidly being generated.

A genomic perspective on HLA evolution

Several decades of research have convincingly shown that classical human leukocyte antigen (HLA) loci bear signatures of natural selection. Despite this conclusion, many questions remain regarding the type of selective regime acting on these loci, the time frame at which selection acts, and the functional connections between genetic variability and natural selection. In this review, we argue that genomic datasets, in particular those generated by next-generation sequencing (NGS) at the population scale, are transforming our understanding of HLA evolution. We show that genomewide data can be used to perform robust and powerful tests for selection, capable of identifying both positive and balancing selection at HLA genes. Importantly, these tests have shown that natural selection can be identified at both recent and ancient timescales. We discuss how findings from genomewide association studies impact the evolutionary study of HLA genes, and how genomic data can be used to survey adaptive change involving interaction at multiple loci. We discuss the methodological developments which are necessary to correctly interpret genomic analyses involving the HLA region. These developments include adapting the NGS analysis framework so as to deal with the highly polymorphic HLA data, as well as developing tools and theory to search for signatures of selection, quantify differentiation, and measure admixture within the HLA region. Finally, we show that high throughput analysis of molecular phenotypes for HLA genes-namely transcription levels-is now a feasible approach and can add another dimension to the study of genetic variation.

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans

Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

Long-term balancing selection contributes to adaptation in Arabidopsis and its relatives

BACKGROUND

In contrast to positive selection, which reduces genetic variation by fixing beneficial alleles, balancing selection maintains genetic variation within a population or species and plays crucial roles in adaptation in diverse organisms. However, which genes, genome-wide, are under balancing selection and the extent to which these genes are involved in adaptation are largely unknown.

RESULTS

We performed a genome-wide scan for genes under balancing selection across two plant species, Arabidopsis thaliana and its relative Capsella rubella, which diverged about 8 million generations ago. Among hundreds of genes with shared coding-region polymorphisms, we find evidence for long-term balancing selection in five genes: AT1G35220, AT2G16570, AT4G29360, AT5G38460, and AT5G44000. These genes are involved in the response to biotic and abiotic stress and other fundamental biochemical processes. More intriguingly, for these genes, we detected significant ecological diversification between the two haplotype groups, suggesting that balancing selection has been very important for adaptation.

CONCLUSIONS

Our results indicate that beyond the well-known S-locus genes and resistance genes, many loci are under balancing selection. These genes are mostly correlated with resistance to stress or other fundamental functions and likely play a more important role in adaptation to diverse habitats than previously thought.

Molecular characterisation, genetic variability and detection of a functional polymorphism influencing the promoter activity of OXT gene in goat and sheep

The purpose of the study described in this Research Communication was to report the full characterisation of the goat and sheep oxytocin-neurophysin I gene (OXT), their promoters and amino acid sequences. Using the genomic DNA as template, we sequenced and compared the whole OXT gene (3 exons), plus 958/960 nucleotides at the 5' flanking region and 478/477 nucleotides at the 3' flanking region, in 46 sheep and 24 goats belonging to different breeds/genetic types reared in Italy, Greece and Germany. The comparison of the obtained sequences showed a high degree of genetic variability at these loci. In particular, we focused on the SNP g.438T > C as possible example of trans-specific polymorphism. This SNP alters a putative binding site of the transcription factor Oct-1. The set-up of a luciferase assay confirmed that the C variant of this SNP negatively affects the promoter activity of the sheep OXT gene. The results of this study suggest that the SNP g.438T > C might be useful to promote association studies with traits/physiological processes controlled by this hormone.

Selective Landscapes in newt Immune Genes Inferred from Patterns of Nucleotide Variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome.

The psoriasis-associated deletion of late cornified envelope genes LCE3B and LCE3C has been maintained under balancing selection since Human Denisovan divergence

BACKGROUND

A common, 32kb deletion of LCE3B and LCE3C genes is strongly associated with psoriasis. We recently found that this deletion is ancient, predating Human-Denisovan divergence. However, it was not clear why negative selection has not removed this deletion from the population.

RESULTS

Here, we show that the haplotype block that harbors the deletion (i) retains high allele frequency among extant and ancient human populations; (ii) harbors unusually high nucleotide variation (π, P < 4.1 × 10^-3); (iii) contains an excess of intermediate frequency variants (Tajima's D, P < 3.9 × 10^-3); and (iv) has an unusually long time to coalescence to the most recent common ancestor (TSel, 0.1 quantile).

CONCLUSIONS

Our results are most parsimonious with the scenario where the LCE3BC deletion has evolved under balancing selection in humans. More broadly, this is consistent with the hypothesis that a balance between autoimmunity and natural vaccination through increased exposure to pathogens maintains this deletion in humans.