Functional CDKN2A assay identifies frequent deleterious alleles misclassified as variants of uncertain significance

Pathogenic germline CDKN2A variants are associated with an increased risk of pancreatic ductal adenocarcinoma (PDAC). CDKN2A variants of uncertain significance (VUSs) are reported in up to 4.3% of patients with PDAC and result in significant uncertainty for patients and their family members as an unknown fraction are functionally deleterious, and therefore, likely pathogenic. Functional characterization of CDKN2A VUSs is needed to reclassify variants and inform clinical management. 29 germline CDKN2A VUSs previously reported in patients with PDAC or in ClinVar were evaluated using a validated in vitro cell proliferation assay. 12 of the 29 CDKN2A VUSs were functionally deleterious (11 VUSs) or potentially functionally deleterious (1 VUS) and were reclassified as likely pathogenic variants. Thus, over 40% of CDKN2A VUSs identified in patients with PDAC are functionally deleterious and likely pathogenic. When incorporating VUSs found to be functionally deleterious, and reclassified as likely pathogenic, the prevalence of pathogenic/likely pathogenic CDKN2A in patients with PDAC reported in the published literature is increased to up to 4.1% of patients, depending on family history. Therefore, CDKN2A VUSs may play a significant, unappreciated role in risk of pancreatic cancer. These findings have significant implications for the counselling and care of patients and their relatives.

Background selection under evolving recombination rates

Background selection (BGS), the effect that purifying selection exerts on sites linked to deleterious alleles, is expected to be ubiquitous across eukaryotic genomes. The effects of BGS reflect the interplay of the rates and fitness effects of deleterious mutations with recombination. A fundamental assumption of BGS models is that recombination rates are invariant over time. However, in some lineages recombination rates evolve rapidly, violating this central assumption. Here, we investigate how recombination rate evolution affects genetic variation under BGS. We show that recombination rate evolution modifies the effects of BGS in a manner similar to a localised change in the effective population size, potentially leading to an underestimation of the genome-wide effects of selection. Furthermore, we find evidence that recombination rate evolution in the ancestors of modern house mice may have impacted inferences of the genome-wide effects of selection in that species.

Synergistic epistasis of the deleterious effects of transposable elements

The replicative nature and generally deleterious effects of transposable elements (TEs) raise an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles. We investigated this population genetic signal in the likely ancestral Drosophila melanogaster population and found evidence supporting the presence of synergistic epistasis among TE insertions, especially TEs expected to exert large fitness impacts. Even though synergistic epistasis of TEs has been predicted to arise through ectopic recombination and TE-mediated epigenetic silencing mechanisms, we only found mixed support for the associated predictions. We observed signals of synergistic epistasis for a large number of TE families, which is consistent with the expectation that such epistatic interaction mainly happens among copies of the same family. Curiously, significant repulsion linkage was also found among TE insertions from different families, suggesting the possibility that synergism of TEs’ deleterious fitness effects could arise above the family level and through mechanisms similar to those of simple mutations. Our findings set the stage for investigating the prevalence and importance of epistatic interactions in the evolutionary dynamics of TEs.

A continent-wide high genetic load in African buffalo revealed by clines in the frequency of deleterious alleles, genetic hitchhiking and linkage disequilibrium

A high genetic load can negatively affect population viability and increase susceptibility to diseases and other environmental stressors. Prior microsatellite studies of two African buffalo (Syncerus caffer) populations in South Africa indicated substantial genome-wide genetic load due to high-frequency occurrence of deleterious alleles. The occurrence of these alleles, which negatively affect male body condition and bovine tuberculosis resistance, throughout most of the buffalo’s range were evaluated in this study. Using available microsatellite data (2–17 microsatellite loci) for 1676 animals from 34 localities (from 25°S to 5°N), we uncovered continent-wide frequency clines of microsatellite alleles associated with the aforementioned male traits. Frequencies decreased over a south-to-north latitude range (average per-locus Pearson r = -0.22). The frequency clines coincided with a multilocus-heterozygosity cline (adjusted R2 = 0.84), showing up to a 16% decrease in southern Africa compared to East Africa. Furthermore, continent-wide linkage disequilibrium (LD) at five linked locus pairs was detected, characterized by a high fraction of positive interlocus associations (0.66, 95% CI: 0.53, 0.77) between male-deleterious-trait-associated alleles. Our findings suggest continent-wide and genome-wide selection of male-deleterious alleles driven by an earlier observed sex-chromosomal meiotic drive system, resulting in frequency clines, reduced heterozygosity due to hitchhiking effects and extensive LD due to male-deleterious alleles co-occurring in haplotypes. The selection pressures involved must be high to prevent destruction of allele-frequency clines and haplotypes by LD decay. Since most buffalo populations are stable, these results indicate that natural mammal populations, depending on their genetic background, can withstand a high genetic load.

Review: Balancing Selection for Deleterious Alleles in Livestock

Harmful alleles can be under balancing selection due to an interplay of artificial selection for the variant in heterozygotes and purifying selection against the variant in homozygotes. These pleiotropic variants can remain at moderate to high frequency expressing an advantage for favorable traits in heterozygotes, while harmful in homozygotes. The impact on the population and selection strength depends on the consequence of the variant both in heterozygotes and homozygotes. The deleterious phenotype expressed in homozygotes can range from early lethality to a slightly lower fitness in the population. In this review, we explore a range of causative variants under balancing selection including loss-of-function variation (i.e., frameshift, stop-gained variants) and regulatory variation (affecting gene expression). We report that harmful alleles often affect orthologous genes in different species, often influencing analogous traits. The recent discoveries are mainly driven by the increasing genomic and phenotypic resources in livestock populations. However, the low frequency and sometimes subtle effects in homozygotes prevent accurate mapping of such pleiotropic variants, which requires novel strategies to discover. After discovery, the selection strategy for deleterious variants under balancing selection is under debate, as variants can contribute to the heterosis effect in crossbred animals in various livestock species, compensating for the loss in purebred animals. Nevertheless, gene-assisted selection is a useful tool to decrease the frequency of the harmful allele in the population, if desired. Together, this review marks various deleterious variants under balancing selection and describing the functional consequences at the molecular, phenotypic, and population level, providing a resource for further study.

Simulation of pedigree vs. fully-informative marker based relationships matrices in a loblolly pine breeding population

Computer simulations of breeding strategies are an essential resource for tree breeders because they allow exploratory analyses into potential long-term impacts on genetic gain and inbreeding consequences without bearing the cost, time, or resource requirements of field experiments. Previous work has modeled the potential long-term implications on inbreeding and genetic gain using random mating and phenotypic selection. Reduction in sequencing costs has enabled the use of DNA marker-based relationship matrices in addition to or in place of pedigree-based allele sharing estimates; this has been shown to provide a significant increase in the accuracy of progeny breeding value prediction. A potential pitfall of genomic selection using genetic relationship matrices is increased coancestry among selections, leading to the accumulation of deleterious alleles and inbreeding depression. We used simulation to compare the relative genetic gain and risk of inbreeding depression within a breeding program similar to loblolly pine, utilizing pedigree-based or marker-based relationships over ten generations. We saw a faster rate of purging deleterious alleles when using a genomic relationship matrix based on markers that track identity-by-descent of segments of the genome. Additionally, we observed an increase in the rate of genetic gain when using a genomic relationship matrix instead of a pedigree-based relationship matrix. While the genetic variance of populations decreased more rapidly when using genomic-based relationship matrices as opposed to pedigree-based, there appeared to be no long-term consequences on the accumulation of deleterious alleles within the simulated breeding strategy.

Genomic evidence for inbreeding depression and purging of deleterious genetic variation in Indian tigers

Increasing habitat fragmentation leads to wild populations becoming small, isolated, and threatened by inbreeding depression. However, small populations may be able to purge recessive deleterious alleles as they become expressed in homozygotes, thus reducing inbreeding depression and increasing population viability. We used whole-genome sequences from 57 tigers to estimate individual inbreeding and mutation load in a small–isolated and two large–connected populations in India. As expected, the small–isolated population had substantially higher average genomic inbreeding (FROH = 0.57) than the large–connected (FROH = 0.35 and FROH = 0.46) populations. The small–isolated population had the lowest loss-of-function mutation load, likely due to purging of highly deleterious recessive mutations. The large populations had lower missense mutation loads than the small–isolated population, but were not identical, possibly due to different demographic histories. While the number of the loss-of-function alleles in the small–isolated population was lower, these alleles were at higher frequencies and homozygosity than in the large populations. Together, our data and analyses provide evidence of 1) high mutation load, 2) purging, and 3) the highest predicted inbreeding depression, despite purging, in the small–isolated population. Frequency distributions of damaging and neutral alleles uncover genomic evidence that purifying selection has removed part of the mutation load across Indian tiger populations. These results provide genomic evidence for purifying selection in both small and large populations, but also suggest that the remaining deleterious alleles may have inbreeding-associated fitness costs. We suggest that genetic rescue from sources selected based on genome-wide differentiation could offset any possible impacts of inbreeding depression.

Transgenic tools for targeted chromosome rearrangements allow construction of balancer chromosomes in non-melanogaster Drosophila species

Perhaps the most valuable single set of resources for genetic studies of Drosophila melanogaster is the collection of multiply-inverted chromosomes commonly known as balancer chromosomes. Balancers prevent the recovery of recombination exchange products within genomic regions included in inversions and allow perpetual maintenance of deleterious alleles in living stocks and the execution of complex genetic crosses. Balancer chromosomes have been generated traditionally by exposing animals to ionizing radiation and screening for altered chromosome structure or for unusual marker segregation patterns. These approaches are tedious and unpredictable, and have failed to produce the desired products in some species. Here I describe transgenic tools that allow targeted chromosome rearrangements in Drosophila species. The key new resources are engineered reporter genes containing introns with yeast recombination sites and enhancers that drive fluorescent reporter genes in multiple body regions. These tools were used to generate a doubly-inverted chromosome 3R in D. simulans that serves as an effective balancer chromosome.

Comparable Number of Genes Having Experienced Positive Selection among Great Ape Species

Alleles that cause advantageous phenotypes with positive selection contribute to adaptive evolution. Investigations of positive selection in protein-coding genes rely on the accuracy of orthology, models, the quality of assemblies, and alignment. Here, based on the latest genome assemblies and gene annotations, we present a comparative analysis on positive selection in four great ape species and identify 211 high-confidence positively selected genes (PSGs). Even the differences in population size among these closely related great apes have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments, we found that they experienced comparable numbers of positive selection. We also uncovered that more than half of multigene families exhibited signals of positive selection, suggesting that imbalanced positive selection resulted in the functional divergence of duplicates. Moreover, at the expression level, although positive selection led to a more non-uniform pattern across tissues, the correlation between positive selection and expression patterns is diverse. Overall, this updated list of PSGs is of great significance for the further study of the phenotypic evolution in great apes.