Using variant databases for variant prioritization and to detect erroneous genotype-phenotype associations

Canopy2: tumor phylogeny inference by bulk DNA and single-cell RNA sequencing

Tumors are comprised of a mixture of distinct cell populations that differ in terms of genetic makeup and function. Such heterogeneity plays a role in the development of drug resistance and the ineffectiveness of targeted cancer therapies. Insight into this complexity can be obtained through the construction of a phylogenetic tree, which illustrates the evolutionary lineage of tumor cells as they acquire mutations over time. We propose Canopy2, a Bayesian framework that uses single nucleotide variants derived from bulk DNA and single-cell RNA sequencing to infer tumor phylogeny and conduct mutational profiling of tumor subpopulations. Canopy2 uses Markov chain Monte Carlo methods to sample from a joint probability distribution involving a mixture of binomial and beta-binomial distributions, specifically chosen to account for the sparsity and stochasticity of the single-cell data. Canopy2 demystifies the sources of zeros in the single-cell data and separates zeros categorized as non-cancerous (cells without mutations), stochastic (mutations not expressed due to bursting), and technical (expressed mutations not picked up by sequencing). Simulations demonstrate that Canopy2 consistently outperforms competing methods and reconstructs the clonal tree with high fidelity, even in situations involving low sequencing depth, poor single-cell yield, and highly-advanced and polyclonal tumors. We further assess the performance of Canopy2 through application to breast cancer and glioblastoma data, benchmarking against existing methods. Canopy2 is an open-source R package available at https://github.com/annweideman/canopy2.

Classification of feline hypertrophic cardiomyopathy-associated gene variants according to the American College of Medical Genetics and Genomics guidelines

Introduction

The correct labeling of a genetic variant as pathogenic is important as breeding decisions based on incorrect DNA tests can lead to the unwarranted exclusion of animals, potentially compromising the long-term health of a population. In human medicine, the American college of Medical Genetics (ACMG) guidelines provide a framework for variant classification. This study aims to apply these guidelines to six genetic variants associated with hypertrophic cardiomyopathy (HCM) in certain cat breeds and to propose a modified criterion for variant classification.

Methods

Genetic samples were sourced from five cat breeds: Maine Coon, Sphynx, Ragdoll, Devon Rex, and British Short- and Longhair. Allele frequencies were determined, and in the subset with phenotypes available, odds ratios to determine the association with HCM were calculated. In silico evaluation followed with joint evidence and data from other publications assisting in the classification of each variant.

Results

Two variants, MYBPC3:c.91G > C [A31P] and MYBPC3:c.2453C > T [R818W], were designated as pathogenic. One variant, MYH7:c.5647G > A [E1883K], was found likely pathogenic, while the remaining three were labeled as variants of unknown significance.

Discussion

Routine genetic testing is advised solely for the MYBPC3:c.91G > C [A31P] in the Maine Coon and MYBPC3:c.2453C > T [R818W] in the Ragdoll breed. The human ACMG guidelines serve as a suitable foundational tool to ascertain which variants to include; however, refining them for application in veterinary medicine might be beneficial.

Identification of a Novel Idiopathic Epilepsy Risk Locus and a Variant in the CCDC85A Gene in the Dutch Partridge Dog

(1) Idiopathic epilepsy (IE) is thought to have a genetic cause in several dog breeds. However, only two causal variants have been identified to date, and few risk loci are known. No genetic studies have been conducted on IE in the Dutch partridge dog (DPD), and little has been reported on the epileptic phenotype in this breed. (2) Owner-filled questionnaires and diagnostic investigations were used to characterize IE in the DPD. A genome-wide association study (GWAS) involving 16 cases and 43 controls was performed, followed by sequencing of the coding sequence and splice site regions of a candidate gene within the associated region. Subsequent whole-exome sequencing (WES) of one family (including one IE-affected dog, both parents, and an IE-free sibling) was performed. (3) IE in the DPD has a broad range in terms of age at onset, frequency, and duration of epileptic seizures. Most dogs showed focal epileptic seizures evolving into generalized seizures. A new risk locus on chromosome 12 (BICF2G630119560; p_raw = 4.4 × 10^-7; p_adj = 0.043) was identified through GWAS. Sequencing of the GRIK2 candidate gene revealed no variants of interest. No WES variants were located within the associated GWAS region. However, a variant in CCDC85A (chromosome 10; XM_038680630.1: c.689C > T) was discovered, and dogs homozygous for the variant (T/T) had an increased risk of developing IE (OR: 6.0; 95% CI: 1.6-22.6). This variant was identified as likely pathogenic according to ACMG guidelines. (4) Further research is necessary before the risk locus or CCDC85A variant can be used for breeding decisions.

The TNNT2:c.95-108G>A variant is common in Maine Coons and shows no association with hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is a common and potentially fatal heart disease in many cat breeds. An intronic variant in TNNT2, c.95-108G>A, was recently reported as the cause of HCM in the Maine Coon. The aim of this study was to determine this variant's allele frequency in different populations and its possible association with HCM. Based on 160 Maine Coon samples collected in Belgium, Italy, Sweden and the USA, the variant's allele frequency was estimated to be 0.32. Analysis of the 99 Lives feline whole genome sequencing database showed that the TNNT2 variant also occurs in other breeds, as well as mixed-breed cats. Comparison of 31 affected and 58 healthy cats did not reveal significantly increased odds for HCM in homozygotes. Based on the combined evidence and in agreement with the standards and guidelines for the interpretation of sequence variants, this variant is currently classified as a variant of unknown significance and should not be used for breeding decisions regarding HCM.

Driving mosaicism: somatic variants in reference population databases and effect on variant interpretation in rare genetic disease

BACKGROUND

Genetic variation databases provide invaluable information on the presence and frequency of genetic variants in the 'untargeted' human population, aggregated with the primary goal to facilitate the interpretation of clinically important variants. The presence of somatic variants in such databases can affect variant assessment in undiagnosed rare disease (RD) patients. Previously, the impact of somatic mosaicism was only considered in relation to two Mendelian disease-associated genes. Here, we expand the analyses to identify additional mosaicism-prone genes in blood-derived reference population databases.

RESULTS

To identify additional mosaicism-prone genes relevant to RDs, we focused on known/previously established ClinVar pathogenic and likely pathogenic single-nucleotide variants, residing in genes associated with early onset, severe autosomal dominant diseases. We asked whether any of these variants are present in a higher-than-expected frequency in the reference population databases and whether there is evidence of somatic origin (i.e., allelic imbalance) rather than germline heterozygosity (~ half of the reads supporting alternative allele). The mosaicism-prone genes identified were further categorized according to the processes they are involved in. Beyond the previously reported ASXL1 and DNMT3A, we identified 7 additional autosomal dominant RD-associated genes with known pathogenic single-nucleotide variants present in the reference population databases and good evidence of allelic imbalance: BRAF, CBL, FGFR3, IDH2, KRAS, PTPN11 and SETBP1. From this group of 9 genes, the majority (n = 7) was important for hematopoiesis. In addition, 4 of these genes were involved in cell proliferation. Further assessment of the known 156 hematopoietic genes led to identification of 48 genes (21 not yet associated with RDs) with at least some evidence of mosaicism detectable in reference population databases.

CONCLUSIONS

These results stress the importance of considering genes involved in hematopoiesis and cell proliferation when interpreting the presence and frequency of genetic variants in blood-derived reference population databases, both public and private. This is especially important when considering new variants of uncertain significance in known hematopoietic/cell proliferation RD genes and future novel gene-disease associations involving this class of genes.

Truncating SLC12A6 variants cause different clinical phenotypes in humans and dogs

Clinical, pathological, and genetic findings of a primary hereditary ataxia found in a Malinois dog family are described and compared with its human counterpart. Based on the family history and the phenotype/genotype relationships already described in humans and dogs, a causal variant was expected to be found in KCNJ10. Rather surprisingly, whole-exome sequencing identified the SLC12A6 NC_006612.3(XM_014109414.2): c.178_181delinsCATCTCACTCAT (p.(Met60Hisfs*14)) truncating variant. This loss-of-function variant perfectly segregated within the affected Malinois family in an autosomal recessive way and was not found in 562 additional reference dogs from 18 different breeds, including Malinois. In humans, SLC12A6 variants cause "agenesis of the corpus callosum with peripheral neuropathy" (ACCPN, alias Andermann syndrome), owing to a dysfunction of this K⁺-Cl^- cotransporter. However, depending on the variant (including truncating variants), different clinical features are observed within ACCPN. The variant in dogs encodes the shortest isoform described so far and its resultant phenotype is quite different from humans, as no signs of peripheral neuropathy, agenesis of the corpus callosum nor obvious mental retardation have been observed in dogs. On the other hand, progressive spinocerebellar ataxia, which is the most important feature of the canine phenotype, hindlimb paresis, and myokymia-like muscle contractions have not been described in humans with ACCPN so far. As this is the first report of a naturally occurring disease-causing SLC12A6 variant in a non-human species, the canine model will be highly valuable to better understand the complex molecular pathophysiology of SLC12A6-related neurological disorders and to evaluate novel treatment strategies.

Comprehensive characterization of horse genome variation by whole-genome sequencing of 88 horses

Whole-genome sequencing studies are vital to gain a thorough understanding of genomic variation. Here, we summarize the results of a whole-genome sequencing study comprising 88 horses and ponies from diverse breeds at 19.1× average coverage. The paired-end reads were mapped to the current EquCab3.0 horse reference genome assembly, and we identified approximately 23.5 million single nucleotide variants and 2.3 million short indel variants. Our dataset included at least 7 million variants that were not previously reported. On average, each individual horse genome carried ∼5.7 million single nucleotides and 0.8 million small indel variants with respect to the reference genome assembly. The variants were functionally annotated. We provide two examples for potentially deleterious recessive alleles that were identified in a heterozygous state in individual genome sequences. Appropriate management of such deleterious recessive alleles in horse breeding programs should help to improve fertility and reduce the prevalence of heritable diseases. This comprehensive dataset has been made publicly available, will represent a valuable resource for future horse genetic studies and supports the goal of accelerating the rates of genetic gain in domestic horse.