BACKGROUND: Wilson disease (WD) is a rare autosomal recessive disorder of copper metabolism caused due to mutations in the copper transporter ATP7B. There is often a striking variability of clinical manifestations among patients with ATP7B mutations, including in siblings. This phenomenon may be caused by individual differences in copper accumulation in hepatocytes and intolerance to copper toxicity as governed by genetic variations in copper metabolism genes acting as modifier loci to the disease.
OBJECTIVE: To elucidate the genetic basis of striking clinical heterogeneity among two siblings of two families with WD.
METHODS: The disease diagnosis and subsequent clinical examinations were performed by expert clinicians. The younger siblings in both families presented with early neurological manifestations at a younger age than their older siblings. Interestingly, only the younger siblings were reported to have had hepatic manifestations. Exome sequencing of all the four individuals was performed to understand their heterogeneous phenotypic outcomes.
RESULTS: Genetic screening revealed no difference in the ATP7B variant spectrum between the siblings of each family. However, the siblings of both the families were found to harbor mutually exclusive pathogenic variants in suspected modifier genes implicated in copper metabolism and/or other neurological and hepatic disorders having overlapping symptoms with WD, viz., CFTR, PPARG, ABCB11, ATP7A, CYP2D6, mTOR, TOR1A, and CP, which can potentially explain their differential clinical phenotypes.
CONCLUSIONS: Clinical heterogeneity between siblings with WD with the same ATP7B mutation profile may be attributed to the presence of different pathogenic variants in potential modifier genes.

Leber congenital amaurosis (LCA) and early-onset retinal degeneration (EORD) are inherited retinal diseases (IRD) characterized by early-onset vision impairment. Herein, we studied 15 Saudi families by whole exome sequencing (WES) and run-of-homozygosity (ROH) detection via AutoMap in 12/15 consanguineous families. This revealed (likely) pathogenic variants in 11/15 families (73%). A potential founder variant was found in RPGRIP1. Homozygous pathogenic variants were identified in known IRD genes (ATF6, CRB1, CABP4, RDH12, RIMS2, RPGRIP1, SPATA7). We established genotype-driven clinical reclassifications for ATF6, CABP4, and RIMS2. Specifically, we observed isolated IRD in the individual with the novel RIMS2 variant, and we found a retina-enriched RIMS2 isoform conserved but not annotated in mouse. The latter illustrates potential different phenotypic consequences of pathogenic variants depending on the particular tissue/cell-type specific isoforms they affect. Lastly, a compound heterozygous genotype in GUCY2D in one non-consanguineous family was demonstrated, and homozygous variants in novel candidate genes ATG2B and RUFY3 were found in the two remaining consanguineous families. Reporting these genes will allow to validate them in other IRD cohorts. Finally, the missing heritability of the two unsolved IRD cases may be attributed to variants in non-coding regions or structural variants that remained undetected, warranting future WGS studies.

To assess the likely pathogenic/pathogenic (LP/P) variants rates in Mendelian dementia genes and the moderate-to-strong risk factors rates in patients with Alzheimer disease (AD).
We included 700 patients in a prospective study and performed exome sequencing. A panel of 28 Mendelian and 6 risk-factor genes was interpreted and returned to patients. We built a framework for risk variant interpretation and risk gradation and assessed the detection rates among early-onset AD (EOAD, age of onset (AOO) ≤65 years, n = 608) depending on AOO and pedigree structure and late-onset AD (66 < AOO < 75, n = 92).
Twenty-one patients carried a LP/P variant in a Mendelian gene (all with EOAD, 3.4%), 20 of 21 affected APP, PSEN1, or PSEN2. LP/P variant detection rates in EOAD ranged from 1.7% to 11.6% based on AOO and pedigree structure. Risk factors were found in 69.5% of the remaining 679 patients, including 83 (12.2%) being heterozygotes for rare risk variants, in decreasing order of frequency, in TREM2, ABCA7, ATP8B4, SORL1, and ABCA1, including 5 heterozygotes for multiple rare risk variants, suggesting non-monogenic inheritance, even in some autosomal-dominant-like pedigrees.
We suggest that genetic screening should be proposed to all EOAD patients and should no longer be prioritized based on pedigree structure.

Sarcopenic obesity (SO) is an age-related disease characterized by the coexistence of excessive adiposity and low muscle mass or function. Although obesity and sarcopenia are heritable conditions, the genetic determinants of SO have not been fully understood. We conducted a large-scale exome-wide association analysis of SO in a sequenced sample of 2 887 cases and 113 284 controls and an imputed sample of 4 003 cases and 161 990 controls in the UK Biobank cohort. Single-variant association analysis identified one locus 1q41 (lead SNP rs1417066, LYPLAL1-AS1, odds ratio [OR] = 1.15, 95% confidence interval [CI] = [1.11-1.19], p = 1.75 × 10-14) that was significantly associated with SO at the exome-wide significance level (p < 1 × 10-8). Colocalization analysis in the Genotype-Tissue Expression expression quantitative trait locus database showed that LYPLAL1-AS1 was colocalized with SO in multiple musculoskeletal-related tissues. Gene-based burden test of rare loss-of-function variants identified 5 genes at the gene-wise significance level (p < 4.3 × 10-6): PDE3B (OR = 2.48, p = 1.10 × 10-6), MYOZ3 (OR = 25.49, p = 1.41 × 10-7), SLC15A3 (OR = 4.75, p = 6.82 × 10-7), RNF130 (OR = 25.83, p = 4.07 × 10-6), and TNK2 (OR = 4.25, p = 8.75 × 10-8). Overall, our study uncovered the genetic effects of both common and rare variants on SO susceptibility, expanded existing knowledge of the genetic architecture of SO, and improved understanding of the genetic mechanisms underlying SO.

UNASSIGNED: Perioperative anaphylaxis (POA) can lead to severe consequences. Identifying clinical risk factors and genetic loci associated with POA through pre-prescription screening may help reduce its incidence.
UNASSIGNED: Using univariate regression and covariate-adjusted multivariate regression, we retrospectively analyzed the association between clinical characteristics and POA in 72 POA patients and 72 non-POA individuals. The discovery study of whole-exome association relied on whole-exome sequencing of 73 POA cases and 1339 healthy individuals. A replication study involving an independent set of 16 POA cases and 1339 healthy individuals confirmed this association. The accurate typing of human leucocyte antigen through exome sequencing (ATHLATES) algorithm and the whole-exome sequencing data were used for genotyping the human leucocyte antigen G (HLA-G) of 73 POA patients. The HLA-G of 16 POA cases and 122 non-POA patients were genotyped through Sanger sequencing. We used Fisher\'s exact probability method to compare the allele and carrier frequencies between POA patients and healthy individuals or non-POA patients. A Pc (P/Bonferroni correction coefficient) < 0.05 represents statistical significance.
UNASSIGNED: Regression analysis identified female sex, an unconfirmed food allergy label, and a history of prior surgery as clinical variables associated with POA. The whole-exome association discovery study identified a strong signal in the major histocompatibility complex region on chromosome 6, with the rs1130356 being the most significant locus (P = 1.5E-10, OR = 3.4, 95% CI = 2.4-4.9). The replication study verified the association between the rs1130356-T allele and POA cases (P = 1.0E-6, OR = 6.3, 95% CI = 3.1-12.7). Compared with non-POA patients, HLA-G∗01:01 (Pc = 2.4E-4, OR = 2.4, 95% CI = 1.6-3.6) was significantly enriched, while HLA-G∗01:04 (Pc = 1.2E-6, OR = 0.3, 95% CI = 0.2-0.5) was lessened in POA patients.
UNASSIGNED: Our study suggested an association between POA and the risk factors of female sex, an unconfirmed food allergy label, and prior surgery. HLA-G, located in the human leucocyte antigen (HLA) region, may act as a surrogate genetic marker for POA. This suggests a causal relationship between this specific genomic region and POA. Our findings shed light on the contribution of human exome genetic variants to the susceptibility to POA.

Type 2 diabetes (T2D) is a multifactorial disease involving genetic and environmental components. Several genome-wide association studies (GWAS) have been conducted to decipher potential genetic aberrations promoting the onset of this metabolic disorder. These GWAS have identified over 400 associated variants, mostly in the intronic or intergenic regions. Recently, a growing number of exome genotyping or exome sequencing experiments have identified coding variants associated with T2D. Such studies were mainly conducted in European populations, and the few candidate-gene replication studies in North African populations revealed inconsistent results. In the present study, we aimed to discover the coding genetic etiology of T2D in the Tunisian population.
We carried out a pilot Exome Wide Association Study (EWAS) on 50 Tunisian individuals. Single variant analysis was performed as implemented in PLINK on potentially deleterious coding variants. Subsequently, we applied gene-based and gene-set analyses using MAGMA software to identify genes and pathways associated with T2D. Potential signals were further replicated in an existing large in-silico dataset, involving up to 177116 European individuals.
Our analysis revealed, for the first time, promising associations between T2D and variations in MYORG gene, implicated in the skeletal muscle fiber development. Gene-set analysis identified two candidate pathways having nominal associations with T2D in our study samples, namely the positive regulation of neuron apoptotic process and the regulation of mucus secretion. These two pathways are implicated in the neurogenerative alterations and in the inflammatory mechanisms of metabolic diseases. In addition, replication analysis revealed nominal associations of the regulation of beta-cell development and the regulation of peptidase activity pathways with T2D, both in the Tunisian subjects and in the European in-silico dataset.
The present study is the first EWAS to investigate the impact of single genetic variants and their aggregate effects on T2D risk in Africa. The promising disease markers, revealed by our pilot EWAS, will promote the understanding of the T2D pathophysiology in North Africa as well as the discovery of potential treatments.

Rare disease genomic testing is a complex process involving various resources. Accurate resource estimation is required for informed prioritization and reimbursement decisions. This study aims to analyze the costs and cost drivers of clinical genomic testing.
Based on genomic sequencing workflows we microcosted limited virtual panel analysis on exome sequencing backbone, proband and trio exome, and genome testing for proband and trio analysis in 2023 Australian Dollars ($). Deterministic and probabilistic sensitivity analyses were undertaken.
Panel testing costs AUD $2373 ($733-$6166), and exome sequencing costs $2823 ($802-$7206) and $5670 ($2006-$11,539) for proband and trio analysis, respectively. Genome sequencing costs $4840 ($2153-$9890) and $11,589 ($5842-$16,562) for proband and trio analysis. The most expensive cost component of genomic testing was sequencing (36.9%-69.4% of total cost), with labor accounting for 27.1%-63.2% of total cost.
We provide a comprehensive analysis of rare disease genomic testing costs, for a range of clinical testing types and contexts. This information will accurately inform economic evaluations of rare disease genomic testing and decision making on policy settings that assist with implementation, such as genomic testing reimbursement.

BACKGROUND: Rare variants are believed to play a substantial role in the genetic architecture of mental disorders, particularly in coding regions. However, limited evidence supports the impact of rare variants on anxiety.
METHODS: Using whole-exome sequencing data from 200,643 participants in the UK Biobank, we investigated the contribution of rare variants to anxiety. Firstly, we computed genetic risk score (GRS) of anxiety utilizing genotype data and summary data from a genome-wide association study (GWAS) on anxiety disorder. Subsequently, we identified individuals within the lowest 50% GRS, a subgroup more likely to carry pathogenic rare variants. Within this subgroup, we classified individuals with the highest 10% 7-item Generalized Anxiety Disorder scale (GAD-7) score as cases (N = 1869), and those with the lowest 10% GAD-7 score were designated as controls (N = 1869). Finally, we conducted gene-based burden tests and single-variant association analyses to assess the relationship between rare variants and anxiety.
RESULTS: Totally, 47,800 variants with MAF ≤0.01 were annotated as non-benign coding variants, consisting of 42,698 nonsynonymous SNVs, 489 nonframeshift substitution, 236 frameshift substitution, 617 stop-gain and 40 stop-loss variants. After variation aggregation, 5066 genes were included in gene-based association analysis. Totally, 11 candidate genes were detected in burden test, such as RNF123 (PBonferroni adjusted = 3.40 × 10-6), MOAP1(PBonferroni adjusted = 4.35 × 10-4), CCDC110 (PBonferroni adjusted = 5.83 × 10-4). Single-variant test detected 9 rare variants, such as rs35726701(RNF123)(PBonferroni adjusted = 3.16 × 10-10) and rs16942615(CAMTA2) (PBonferroni adjusted = 4.04 × 10-4). Notably, RNF123, CCDC110, DNAH2, and CSKMT gene were identified in both tests.
CONCLUSIONS: Our study identified novel candidate genes for anxiety in protein-coding regions, revealing the contribution of rare variants to anxiety.

OBJECTIVE: The purpose of this study was to understand the impact of population diversity and geographic variation on tumor mutation burden (TMB) scores across cancers and its implication on stratification of patients for immune checkpoint inhibitor (ICI) therapy.
METHODS: This retrospective study used whole-exome sequencing (WES) to profile 1,233 Indian patients with cancer across 30 different cancer types and to estimate their TMB scores. A WES-based pipeline was adopted, along with an indigenously developed strategy for arriving at true somatic mutations. A robust unsupervised machine learning approach was used to understand the distribution of TMB scores across different populations and within the population.
RESULTS: The results of the study showed a biphasic distribution of TMB scores in most cancers, with different threshold scores across cancer types. Patients with cancer in India had higher TMB scores compared with the Caucasian patients. We also observed that the TMB score value at 90th percentile (predicting high efficacy to ICI) was high in four different cancer types (sarcoma, ovary, head and neck, and breast) in the Indian cohort as compared with The Cancer Genome Atlas or public cohort. However, in lung and colorectal cancers, the TMB score distribution was similar between the two population cohorts.
CONCLUSIONS: The findings of this study indicate that it is crucial to benchmark both cancer-specific and population-specific TMB distributions to establish a TMB threshold for each cancer in various populations. Additional prospective studies on much larger population across different cancers are warranted to validate this observation to become the standard of care.

An increased risk of neurocognitive deficits, anxiety, and depression has been reported in childhood cancer survivors.
We analyzed associations of neurocognitive deficits, as well as anxiety and depression, with common and rare genetic variants derived from whole-exome sequencing data of acute lymphoblastic leukemia (ALL) survivors from the PETALE cohort. In addition, significant associations were assessed using stratified and multivariable analyses. Next, top-ranking common associations were analyzed in an independent SJLIFE replication cohort of ALL survivors.
Significant associations were identified in the entire discovery cohort (N = 229) between the AK8 gene and changes in neurocognitive function, whereas PTPRZ1, MUC16, TNRC6C-AS1 were associated with anxiety. Following stratification according to sex, the ZNF382 gene was linked to a neurocognitive deficit in males, whereas APOL2 and C6orf165 were associated with anxiety and EXO5 with depression. Following stratification according to prognostic risk groups, the modulatory effect of rare variants on depression was additionally found in the CYP2W1 and PCMTD1 genes. In the replication SJLIFE cohort (N = 688), the male-specific association in the ZNF382 gene was not significant; however, a P value<0.05 was observed when the entire SJLIFE cohort was analyzed. ZNF382 was significant in males in the combined cohorts as shown by meta-analyses as well as the depression-associated gene EXO5.
Further research is needed to confirm whether the current findings, along with other known risk factors, may be valuable in identifying patients at increased risk of these long-term complications.
Our results suggest that specific genes may be related to increased neuropsychological consequences.