Genetic architecture and adaptations of Nunavik Inuit.

The Canadian Inuit have a distinct population background that may entail particular implications for the health of its individuals. However, the number of genetic studies examining this Inuit population is limited, and much remains to be discovered in regard to its genetic characteristics. In this study, we generated whole-exome sequences and genomewide genotypes for 170 Nunavik Inuit, a small and isolated founder population of Canadian Arctic indigenous people. Our study revealed the genetic background of Nunavik Inuit to be distinct from any known present-day population. The majority of Nunavik Inuit show little evidence of gene flow from European or present-day Native American peoples, and Inuit living around Hudson Bay are genetically distinct from those around Ungava Bay. We also inferred that Nunavik Inuit have a small effective population size of 3,000 and likely split from Greenlandic Inuit ∼10.5 kya. Nunavik Inuit went through a bottleneck at approximately the same time and might have admixed with a population related to the Paleo-Eskimos. Our study highlights population-specific genomic signatures in coding regions that show adaptations unique to Nunavik Inuit, particularly in pathways involving fatty acid metabolism and cellular adhesion (CPNE7, ICAM5, STAT2, and RAF1). Subsequent analyses in selection footprints and the risk of intracranial aneurysms (IAs) in Nunavik Inuit revealed an exonic variant under weak negative selection to be significantly associated with IA (rs77470587; P = 4.6 × 10-8).

RNF213 Is Associated with Intracranial Aneurysms in the French-Canadian Population.

Intracranial aneurysms (IAs) are the result of focal weakness in the artery wall and have a complex genetic makeup. To date, genome-wide association and sequencing studies have had limited success in identifying IA risk factors. Distinct populations, such as the French-Canadian (FC) population, have increased IA prevalence. In our study, we used exome sequencing to prioritize risk variants in a discovery cohort of six FC families affected by IA, and the analysis revealed an increased variation burden for ring finger protein 213 (RNF213). We resequenced RNF213 in a larger FC validation cohort, and association tests on further identified variants supported our findings (SKAT-O, p = 0.006). RNF213 belongs to the AAA+ protein family, and two variants (p.Arg2438Cys and p.Ala2826Thr) unique to affected FC individuals were found to have increased ATPase activity, which could lead to increased risk of IA by elevating angiogenic activities. Common SNPs in RNF213 were also extracted from the NeuroX SNP-chip genotype data, comprising 257 FC IA-affected and 1,988 control individuals. We discovered that the non-ancestral allele of rs6565666 was significantly associated with the affected individuals (p = 0.03), and it appeared as though the frequency of the risk allele had changed through genetic drift. Although RNF213 is a risk factor for moyamoya disease in East Asians, we demonstrated that it might also be a risk factor for IA in the FC population. It therefore appears that the function of RNF213 can be differently altered to predispose distinct populations to dissimilar neurovascular conditions, highlighting the importance of a population's background in genetic studies of heterogeneous disease.

Multiple Pairwise Analysis of Non-homologous Centromere Coupling Reveals Preferential Chromosome Size-Dependent Interactions and a Role for Bouquet Formation in Establishing the Interaction Pattern.

During meiosis, chromosomes undergo a homology search in order to locate their homolog to form stable pairs and exchange genetic material. Early in prophase, chromosomes associate in mostly non-homologous pairs, tethered only at their centromeres. This phenomenon, conserved through higher eukaryotes, is termed centromere coupling in budding yeast. Both initiation of recombination and the presence of homologs are dispensable for centromere coupling (occurring in spo11 mutants and haploids induced to undergo meiosis) but the presence of the synaptonemal complex (SC) protein Zip1 is required. The nature and mechanism of coupling have yet to be elucidated. Here we present the first pairwise analysis of centromere coupling in an effort to uncover underlying rules that may exist within these non-homologous interactions. We designed a novel chromosome conformation capture (3C)-based assay to detect all possible interactions between non-homologous yeast centromeres during early meiosis. Using this variant of 3C-qPCR, we found a size-dependent interaction pattern, in which chromosomes assort preferentially with chromosomes of similar sizes, in haploid and diploid spo11 cells, but not in a coupling-defective mutant (spo11 zip1 haploid and diploid yeast). This pattern is also observed in wild-type diploids early in meiosis but disappears as meiosis progresses and homologous chromosomes pair. We found no evidence to support the notion that ancestral centromere homology plays a role in pattern establishment in S. cerevisiae post-genome duplication. Moreover, we found a role for the meiotic bouquet in establishing the size dependence of centromere coupling, as abolishing bouquet (using the bouquet-defective spo11 ndj1 mutant) reduces it. Coupling in spo11 ndj1 rather follows telomere clustering preferences. We propose that a chromosome size preference for centromere coupling helps establish efficient homolog recognition.

Exome sequencing of sporadic childhood-onset schizophrenia suggests the contribution of X-linked genes in males.

Childhood-onset schizophrenia (COS) is a rare and severe form of schizophrenia, defined as having an onset before the age of 13. The male COS cases have a slightly younger age of onset than female cases. They also present with a higher rate of comorbid developmental disorders. These sex differences are not explained by the frequency of chromosomal abnormalities, and the contribution of other forms of genetic variations remains unestablished. Using a whole-exome sequencing approach, we examined 12 COS trios where the unaffected parents had an affected male child. The sequencing data enabled us to test if the hemizygous variants, transmitted from the unaffected carrying mother, could mediate the phenotype (X-linked recessive inheritance model). Our results revealed that affected children have a significantly greater number of X-linked rare variants than their unaffected fathers. The variants identified in the male probands were mostly found in genes previously linked to other neuropsychiatric diseases like autism, intellectual disability, and epilepsy, including LUZP4, PCDH19, RPS6KA3, and OPHN1. The level of expression of the genes was assessed at different developmental periods in normal brain using the BrainSpan database. This approach revealed that some of them were expressed earlier in males than in females, consistent with the younger age of onset in male COS. In conclusion, this article suggests that X-linked genes might play a role in the pathophysiology of COS. Candidate genes detailed here could explain the higher level of comorbidities and the earlier age of onset observed in a subset of the male COS cases.

Efficacy, safety, and comparison of sonic hedgehog inhibitors in basal cell carcinomas: A systematic review and meta-analysis.

BACKGROUND: Sonic hedgehog inhibitors (SHHis) provide an additional treatment option for basal cell carcinomas (BCCs), especially for metastatic or locally advanced BCC. However, studies have been heterogeneous and lacked direct comparisons between molecules. OBJECTIVE: To determine the efficacy and safety of the class of molecules SHHi for treating BCC and to compare them individually. METHODS: We performed a PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses)-compliant systematic review of studies followed by a meta-analysis. RESULTS: Eighteen articles were included in our meta-analysis; 16 articles were combined for efficacy and 16 for safety. In locally advanced BCC, overall response rates (ORRs) were similar for vismodegib and sonidegib (69% vs 57%, respectively) but not complete response rates (31% vs 3%, respectively). In metastatic disease, the ORR of vismodegib was 2.7-fold higher than the ORR of sonidegib (39% vs 15%, respectively). For side effects affecting a majority of patients, prevalences for muscle spasms (67.1%), dysgeusia (54.1%), and alopecia (57.7%) were in similar proportions for sonidegib and vismodegib. Patients receiving sonidegib experienced more upper gastrointestinal distress than patients receiving vismodegib. CONCLUSION: SHHis induce a partial response to locally advanced BCC disease. Side effects are common, similar across molecules, associated with high discontinuation rates, and warrant discussion beforehand.

Whole-exome sequencing in undiagnosed genetic diseases: interpreting 119 trios.

PURPOSE: Despite the recognized clinical value of exome-based diagnostics, methods for comprehensive genomic interpretation remain immature. Diagnoses are based on known or presumed pathogenic variants in genes already associated with a similar phenotype. Here, we extend this paradigm by evaluating novel bioinformatics approaches to aid identification of new gene-disease associations. METHODS: We analyzed 119 trios to identify both diagnostic genotypes in known genes and candidate genotypes in novel genes. We considered qualifying genotypes based on their population frequency and in silico predicted effects we also characterized the patterns of genotypes enriched among this collection of patients. RESULTS: We obtained a genetic diagnosis for 29 (24%) of our patients. We showed that patients carried an excess of damaging de novo mutations in intolerant genes, particularly those shown to be essential in mice (P = 3.4 × 10(-8)). This enrichment is only partially explained by mutations found in known disease-causing genes. CONCLUSION: This work indicates that the application of appropriate bioinformatics analyses to clinical sequence data can also help implicate novel disease genes and suggest expanded phenotypes for known disease genes. These analyses further suggest that some cases resolved by whole-exome sequencing will have direct therapeutic implications.

The utility of the traditional medical genetics diagnostic evaluation in the context of next-generation sequencing for undiagnosed genetic disorders.

PURPOSE: The purpose of this study was to assess the diagnostic yield of the traditional, comprehensive clinical evaluation and targeted genetic testing, within a general genetics clinic. These data are critically needed to develop clinically and economically grounded diagnostic algorithms that consider presenting phenotype, traditional genetics testing, and the emerging role of next-generation sequencing (whole-exome/genome sequencing). METHODS: We retrospectively analyzed a cohort of 500 unselected consecutive patients who received traditional genetic diagnostic evaluations at a tertiary medical center. We calculated the diagnosis rate, number of visits to diagnosis, genetic tests, and the cost of testing. RESULTS: Thirty-nine patients were determined to not have a genetic disorder; 212 of the remaining 461 (46%) received a genetic diagnosis, and 72% of these were diagnosed on the first visit. The cost per subsequent successful genetic diagnosis was estimated at $25,000. CONCLUSION: Almost half of the patients were diagnosed using the traditional approach, most at the initial visit. For those remaining undiagnosed, next-generation sequencing may be clinically and economically beneficial. Estimating a 50% success rate for next-generation sequencing in undiagnosed genetic disorders, its application after the first clinical visit could result in a higher rate of genetic diagnosis at a considerable cost savings per successful diagnosis.

Deep resequencing of 17 glutamate system genes identifies rare variants in DISC1 and GRIN2B affecting risk of opioid dependence.

The N-methyl-D-aspartate (NMDA) glutamate receptors play important roles in the pathophysiology of substance dependence (SD), but no strong genetic evidence has associated common variants in NMDAR-related genes to SD. We hypothesized that rare variants (RVs) with minor allele frequency <1% in the NMDAR-related genes might exert large effects on SD risk. We sequenced 34 544 bp of coding and flanking intronic regions of 17 genes involved in the NMDA system in 760 subjects, all with co-occurring alcohol dependence, cocaine dependence and opioid dependence (OD), and 760 healthy control subjects. One hundred percent of the target regions were sequenced at >1000× coverage. We identified 454 variants, including 380 RVs. Based on case-control allele count differences, we genotyped 11 exonic RVs in 6751 additional subjects, and the 1520 subjects from the sequencing stage for validation. All alleles of the 11 RVs called in the sequencing stage were confirmed. We found a statistically significant association of the 11 RVs with OD in African Americans (P = 0.00080). Results from gene-based association tests showed that the association signal derived mostly from DISC1 (P = 0.0010) and GRIN2B (P = 0.00085). DISC1 is a well-validated schizophrenia risk gene. This is the first demonstration that RVs affect the risk of OD and the first demonstration of biological convergence of schizophrenia and OD risk-via DISC1.

Efficacy and Safety of Sonic Hedgehog Inhibitors in Basal Cell Carcinomas: An Updated Systematic Review and Meta-analysis (2009-2022).

BACKGROUND: Basal cell carcinoma (BCC) of the skin is the most common form of skin cancer in the United States. In life-threatening, advanced BCC, sonic hedgehog inhibitors (SSHis) remain a pre-eminent treatment option for locally advanced BCC and metastatic BCC. OBJECTIVE: In this updated systematic review and meta-analysis, we aimed to better characterize the efficacy and safety of SSHis by including final updates from pivotal clinical trials and additional new recent studies. METHODS: An electronic database search was performed for articles including clinical trials, prospective case series, and retrospective medical record reviews on human subjects. Overall response rates (ORRs) and complete response rates (CRRs) were the primary outcomes. For safety assessment, the prevalence of the following adverse effects was analyzed: muscle spasms, dysgeusia, alopecia, weight loss, fatigue, nausea, myalgias, vomiting, skin squamous cell carcinoma, increased creatine kinase, diarrhea, decreased appetite, and amenorrhea. Analyses were performed using R statistical software. Data were pooled using linear models with fixed effects meta-analysis for primary analyses, along with 95% confidence intervals (CIs) and p-values. Intermolecular differences were calculated using Fisher's exact test. RESULTS: A total of 22 studies (N = 2384 patients) were included in the meta-analysis: 19 studies assessing both efficacy and safety, 2 studies assessing safety only, and 1 study assessing efficacy only. Overall, the pooled ORR for all patients was 64.9% (95% CI 48.2-81.6%), implicating there is at least a partial response (z = 7.60, p < 0.0001) in most patients receiving SSHis. The ORR for vismodegib was 68.5% and 50.1% for sonidegib. The most common adverse effects for vismodegib and sonidegib were muscle spasms (70.5% and 61.0%, respectively), dysgeusia (58.4% and 48.6%, respectively), and alopecia (59.9% and 51.1%, respectively). Patients were likely to experience weight loss (35.1%, p < 0.0001) from vismodegib. Alternatively, patients taking sonidegib experienced more nausea, diarrhea, increased creatine kinase levels, and decreased appetite compared with those receiving vismodegib. CONCLUSION: SSHis are an effective treatment for advanced BCC disease. Given the high discontinuation rates, management of patient expectations is warranted for compliance and achieving long-term efficacy. It is essential to stay updated with the latest discoveries on the efficacy and safety of SSHis.

Gene-Environment Analyses in a UK Biobank Skin Cancer Cohort Identifies Important SNPs in DNA Repair Genes That May Help Prognosticate Disease Risk.

BACKGROUND: Despite well-established relationships between sun exposure and skin cancer pathogenesis/progression, specific gene-environment interactions in at-risk individuals remain poorly-understood. METHODS: We leveraged a UK Biobank cohort of basal cell carcinoma (BCC, n = 17,221), cutaneous squamous cell carcinoma (cSCC, n = 2,331), melanoma in situ (M-is, n = 1,158), invasive melanoma (M-inv, n = 3,798), and healthy controls (n = 448,164) to quantify the synergistic involvement of genetic and environmental factors influencing disease risk. We surveyed 8,798 SNPs from 190 DNA repair genes, and 11 demographic/behavioral risk factors. RESULTS: Clinical analysis identified darker skin (RR = 0.01-0.65) and hair (RR = 0.27-0.63) colors as protective factors. Eleven SNPs were significantly associated with BCC, three of which were also associated with M-inv. Gene-environment analysis yielded 201 SNP-environment interactions across 90 genes (FDR-adjusted q < 0.05). SNPs from the FANCA gene showed interactions with at least one clinical factor in all cancer groups, of which three (rs9926296, rs3743860, rs2376883) showed interaction with nearly every factor in BCC and M-inv. CONCLUSIONS: We identified novel risk factors for keratinocyte carcinomas and melanoma, highlighted the prognostic value of several FANCA alleles among individuals with a history of sunlamp use and childhood sunburns, and demonstrated the importance of combining genetic and clinical data in disease risk stratification. IMPACT: This study revealed genome-wide associations with important implications for understanding skin cancer risk in the context of the rapidly-evolving field of precision medicine. Major individual factors (including sex, hair and skin color, and sun protection use) were significant mediators for all skin cancers, interacting with >200 SNPs across four skin cancer types.

Comparison of the Basal Cell Carcinoma (BCC) Tumour Microenvironment to Other Solid Malignancies.

Basal cell carcinoma (BCC) is the most common form of skin cancer, contributing to nearly a third of new cancer cases in Western countries. Most BCCs are considered low risk "routine" lesions that can either be excised through surgery or treated with chemotherapeutic agents. However, around 1-2% of BCC cases are locally aggressive, present a high risk of metastasis, and often develop chemoresistance, termed advanced BCC. There currently exists no animal model or cell line that can recapitulate advanced BCC, let alone intermediate-risk and high-risk early BCC. We previously found that aggressive BCC tumours presented a Th2 cytokine inflammation profile, mesenchymal stem cell properties, and macrophage-induced tumoral inflammation. In this study, we aimed to identify potential BCC "relatives" among solid-organ malignancies who present similar immune cell proportions in their microenvironment compositions. Using immune cell type deconvolution by CIBERSORTx, and cell type enrichment by xCell, we determined three cancers with the most similar tumour microenvironments as compared to BCC. Specifically, chromophobe renal cell carcinoma, sarcoma, and skin cutaneous melanoma presented significance in multiple cell types, namely in CD4+ T lymphocytes, gammadelta T lymphocytes, and NK cell populations. Consequently, further literature analysis was conducted to understand similarities between BCC and its "relatives", as well as investigating novel treatment targets. By identifying cancers most like BCC, we hope to propose prospective druggable pathways, as well as to gain insight on developing a reliable animal or cell line model to represent advanced BCC.

Ectopically Expressed Meiosis-Specific Cancer Testis Antigen HORMAD1 Promotes Genomic Instability in Squamous Cell Carcinomas.

Genomic instability is a prominent hallmark of cancer, however the mechanisms that drive and sustain this process remain elusive. Research demonstrates that numerous cancers with increased levels of genomic instability ectopically express meiosis-specific genes and undergo meiomitosis, the clash of mitotic and meiotic processes. These meiotic genes may represent novel therapeutic targets for the treatment of cancer. We studied the relationship between the expression of the meiosis protein HORMAD1 and genomic instability in squamous cell carcinomas (SCCs). First, we assessed markers of DNA damage and genomic instability following knockdown and overexpression of HORMAD1 in different cell lines representing SCCs and epithelial cancers. shRNA-mediated depletion of HORMAD1 expression resulted in increased genomic instability, DNA damage, increased sensitivity to etoposide, and decreased expression of DNA damage response/repair genes. Conversely, overexpression of HORMAD1 exhibited protective effects leading to decreased DNA damage, enhanced survival and decreased sensitivity to etoposide. Furthermore, we identified a meiotic molecular pathway that regulates HORMAD1 expression by targeting the upstream meiosis transcription factor STRA8. Our results highlight a specific relationship between HORMAD1 and genomic instability in SCCs, suggesting that selectively inhibiting HORMAD1, possibly, through STRA8 signaling, may provide a new paradigm of treatment options for HORMAD1-expressing SCCs.

Serotonin transporter 5-HTTLPR genotype moderates the effects of childhood adversity on posttraumatic stress disorder risk: a replication study.

We reported that the 5-HTTLPR polymorphism in the promoter region of the serotonin transporter gene (SLC6A4) moderates the effect of childhood adversity on posttraumatic stress disorder (PTSD) risk [Xie et al. (2009); Arch Gen Psychiatry 66 (11): 1201-1209]. In the present study, we considered 5,178 subjects (a group with generally high substance dependence comorbidity, as for our previous study) using similar methodology to replicate our previous results. We used logistic regression analyses to explore the interaction effect of 5-HTTLPR genotype and childhood adversity on PTSD risk. We found that, as reported in our previous study, in individuals with childhood adversity, the presence of one or two copies of the S allele of 5-HTTLPR increased the risk to develop PTSD. This gene-environment interaction effect was present in European Americans (EAs), but not in African Americans (AAs; EAs, OR = 1.49, 95% CI = 1.07-2.08, P = 0.019; AAs, OR = 0.90, 95% CI = 0.60-1.35, P = 0.62). The statistical power to detect this interaction effect was increased when data were combined with those from our previous study [Xie et al. (2009); Arch Gen Psychiatry 66 (11): 1201-1209]. The findings reported here replicate those from our previous work, adding to a growing body of research demonstrating that the 5-HTTLPR genotype moderates risk for anxiety and depression phenotypes in the context of stress and adverse events.

Analysis of multiple basal cell carcinomas (BCCs) arising in one individual highlights genetic tumor heterogeneity and identifies novel driver mutations.

Basal cell carcinoma (BCC) is the most common human cancer, especially in individuals with light skin phototypes (i.e., Fitzpatrick I-II skin type). Many affected develop multiple BCCs during their lifetime. It is not uncommon to observe elderly patients with >5 BCCs. In this study, we explored whether for patients diagnosed with multiple BCCs, analyzing the genomic mutations in one tumor could be sufficient to derive meaningful molecular/genetic conclusions regarding the other BCC tumors. Following the Genome Analysis Toolkit (GATK) best practices we have completed the study of 6 BCCs that occurred in an 83-year-old Caucasian male due to sun exposure. We have analyzed exome sequencing data of each BCC tumor and matched normal skin samples. We identified that BCCs from the same patient shared some of the key driver mutations, but they also displayed significant intertumoral heterogeneity. This finding may in part explain the different clinical progression/evolution of BCCs observed in the same patient. This work also highlights the value of characterizing multiple BCCs in one individual to identify patient-specific genetic events with a potential link to other malignancies and implications for personalized medicine.

A genome-wide screen for variants influencing certolizumab pegol response in a moderate to severe rheumatoid arthritis population.

Certolizumab pegol (CZP) is a PEGylated Fc-free tumor necrosis factor (TNF) inhibitor antibody approved for use in the treatment of rheumatoid arthritis (RA), Crohn's disease, psoriatic arthritis, axial spondyloarthritis and psoriasis. In a clinical trial of patients with severe RA, CZP improved disease symptoms in approximately half of patients. However, variability in CZP efficacy remains a problem for clinicians, thus, the aim of this study was to identify genetic variants predictive of CZP response. We performed a genome-wide association study (GWAS) of 302 RA patients treated with CZP in the REALISTIC trial to identify common single nucleotide polymorphisms (SNPs) associated with treatment response. Whole-exome sequencing was also performed for 74 CZP extreme responders and non-responders within the same population, as well as 1546 population controls. No common SNPs or rare functional variants were significantly associated with CZP response, though a non-significant enrichment in the RA-implicated KCNK5 gene was observed. Two SNPs near spondin-1 and semaphorin-4G approached genome-wide significance. The results of the current study did not provide an unambiguous predictor of CZP response.

The transcriptional landscape analysis of basal cell carcinomas reveals novel signalling pathways and actionable targets.

Basal cell carcinoma (BCC) is the most common skin cancer and human malignancy. Although most BCCs are easily managed, some are aggressive locally, require Mohs micrographic surgery, or can even metastasize. In the latter, resistance to Sonic Hedgehog inhibitors may occur. Despite their frequent occurrence in clinical practice, their transcriptional landscape remains poorly understood. By analyzing BCC RNA sequencing data according to clinically important features (all BCCs versus normal skin, high-risk versus low-risk BCCs based solely on histopathological subtypes with aggressive features, advanced versus non-advanced BCCs, and vismodegib-resistant versus vismodegib-sensitive tumors), we have identified novel differentially regulated genes and new targetable pathways implicated in BCC tumorigenesis. Pathways as diverse as IL-17, TLR, Akt/PI3K, cadherins, integrins, PDGF, and Wnt/ß-catenin are promising therapeutic avenues for local and systemic agents in managing this common malignancy, including through drug re-purposing of existing medications. We experimentally validated several of these targets as biomarkers in our patient-derived cohort of primary BCC tumors.

In silico analyses of the tumor microenvironment highlight tumoral inflammation, a Th2 cytokine shift and a mesenchymal stem cell-like phenotype in advanced in basal cell carcinomas.

Basal Cell Carcinoma (BCC) represents the most common form of all cancers. BCC is characteristically surrounded by a fibromyxoid stroma. Previous studies have suggested a shift towards a Th2 response, an increase in T regulatory lymphocytes and the presence of cancer-associated fibroblasts in the BCC tumor microenvironment. In this study, we aimed to further characterize the BCC tumor microenvironment in detail by analyzing BCC RNA-Sequencing data and correlating it with clinically-relevant features via in silico RNA deconvolution. Using immune cell type deconvolution by CIBERSORT, we have identified a brisk lymphocytic infiltration, and more abundant macrophages in BCC tumors compared to normal skin. Using cell type enrichment by xCell, we confirmed the observed immune infiltration in BCC tumors and compared them to normal skin. We observed a shift towards Th2 immunity in advanced and vismodegib-resistant tumors. Tumoral inflammation induced by macrophage activity was associated with advanced BCCs, while lymphocytic infiltration was most significant in non-advanced tumors, likely related to an adaptive anti-tumoral response. In advanced and vismodegib-resistant BCCs, mesenchymal stem cell-like properties were observed. Particularly in vismodegib-resistant BCCs, fibroblasts and adipocytes were found at high number, which ultimately may contribute to the decreased drug delivery to the tumor. In conclusion, this study has revealed notable BCC tumor microenvironment findings associated with important clinical features. Microenvironment-altering agents may be used locally for "routine" BCCs and systematically for advanced or resistant BCCs.