Potential Causal Variants Research Articles

ATP2B4 regulatory genetic variants are associated with mild malaria

BackgroundGenome-wide association studies have identified ATP2B4 as a severe malaria resistance gene. Recently, 8 potential causal regulatory variants have been shown to be associated with severe malaria.MethodsGenotyping of rs10900585, rs11240734, rs1541252, rs1541253, rs1541254, rs1541255, rs10751450, rs10751451 and rs10751452 was performed in 154 unrelated individuals (79 controls and 75 mild malaria patients). rs10751450, rs10751451 and rs10751452 were genotyped by Taqman assays, whereas the fragment of the ATP2B4 gene containing the remaining SNPs was sequenced. Logistic regression analysis was used to assess the association between the SNPs and mild malaria.ResultsThe results showed that mild malaria was associated with rs10900585, rs11240734, rs1541252, rs1541253, rs1541254, rs1541255, rs10751450, rs10751451 and rs10751452. The homozygous genotypes for the major alleles were associated with an increased risk of mild malaria. Furthermore, the haplotype containing the major alleles and that containing the minor alleles were the most frequent haplotypes. Individuals with the major haplotypes had a significantly higher risk of mild malaria compared to the carriers of the minor allele haplotype.ConclusionsATP2B4 polymorphisms that have been associated with severe malaria are also associated with mild malaria.

Combining multi-omics approaches to prioritize the variant-regulated functional long non-coding RNAs in autism spectrum disorder.

Rising evidence has indicated that long non-coding RNA (lncRNA) may play an essential role in the development of autism spectrum disorder (ASD). However, identifying the lncRNAs associated with ASD and the risk loci on them remains a major challenge. This study aims to identify potential causative variants and explore the related mechanisms. By leveraging differential expression analysis, WGCNA analysis and cis-expression quantitative analysis, our study mined functional SNPs with the regulated long non-coding RNA genes in brain tissues. We recruited 611 ASD children and 645 healthy children in the case-control study. Total 68 different expressed lncRNAs were validated by calculating the brain tissue-specific expression using RNA-seq data. By the WGCNA method, 9 functional lncRNAs classified as e-lncRNA were found to interact with 976 ASD risk genes. Furthermore, we mined functional SNPs regulated long non-coding RNAs in brain tissues. We analyzed the association between candidate SNPs and ASD risks in Chinese children, which showed BDNF-AS rs1565228 allele G to C reduced the risk of ASD (OR = 0.81, 95%CI: 0.66-0.98). Further bioinformatics analysis showed that the variant rs1565228 C>G with the low binding affinity of transcription factor SRF caused the decreased expression of lncRNA BDNF-AS. Our study revealed that rs2295412 in the non-coding sequence of the lncRNA gene region was significantly associated with the risk of ASD. These findings suggested that the SNPs in the non-coding region of lncRNA may play important roles in the genetic susceptibility of ASD, which may facilitate the early screening of ASD.

Meta-Analysis and Multivariate GWAS Analyses in 80,950 Individuals of African Ancestry Identify Novel Variants Associated with Blood Pressure Traits.

High blood pressure (HBP) has been implicated as a major risk factor for cardiovascular diseases in several populations, including individuals of African ancestry. Despite the elevated burden of HBP-induced cardiovascular diseases in Africa and other populations of African descent, limited genetic studies have been carried out to explore the genetic mechanism driving this phenomenon. We performed genome-wide association univariate and multivariate analyses of both systolic (SBP) and diastolic blood pressure (DBP) traits in 80,950 individuals of African ancestry. We used summary statistics data from six independent cohorts, including the African Partnership for Chronic Disease Research (APCDR), the UK Biobank, and the Million Veteran Program (MVP). FUMA was used to annotate, prioritize, visualize, and interpret our findings to gain a better understanding of the molecular mechanism(s) underlying the genetics of BP traits. Finally, we undertook a Bayesian fine-mapping analysis to identify potential causal variants. Our meta-analysis identified 10 independent variants associated with SBP and 9 with DBP traits. Whilst our multivariate GWAS method identified 21 independent signals, 18 of these SNPs have been previously identified. SBP was linked to gene sets involved in biological processes such as synapse assembly and cell-cell adhesion via plasma membrane adhesion. Of the 19 independent SNPs identified in the BP meta-analysis, only 11 variants had posterior probability (PP) of > 50%, including one novel variant: rs562545 (MOBP, PP = 77%). To facilitate further research and fine-mapping of high-risk loci/variants in highly susceptible groups for cardiovascular disease and other related traits, large-scale genomic datasets are needed. Our findings highlight the importance of including ancestrally diverse populations in large GWASs and the need for diversity in genetic research.

Sequencing-based fine-mapping and in silico functional characterization of the 10q24.32 arsenic metabolism efficiency locus across multiple arsenic-exposed populations.

Inorganic arsenic is highly toxic and carcinogenic to humans. Exposed individuals vary in their ability to metabolize arsenic, and variability in arsenic metabolism efficiency (AME) is associated with risks of arsenic-related toxicities. Inherited genetic variation in the 10q24.32 region, near the arsenic methyltransferase (AS3MT) gene, is associated with urine-based measures of AME in multiple arsenic-exposed populations. To identify potential causal variants in this region, we applied fine mapping approaches to targeted sequencing data generated for exposed individuals from Bangladeshi, American Indian, and European American populations (n = 2,357, 557, and 648 respectively). We identified three independent association signals for Bangladeshis, two for American Indians, and one for European Americans. The size of the confidence sets for each signal varied from 4 to 85 variants. There was one signal shared across all three populations, represented by the same SNP in American Indians and European Americans (rs191177668) and in strong linkage disequilibrium (LD) with a lead SNP in Bangladesh (rs145537350). Beyond this shared signal, differences in LD patterns, minor allele frequency (MAF) (e.g., rs12573221 ~13% in Bangladesh ~0.2% among American Indians), and/or heterogeneity in effect sizes across populations likely contributed to the apparent population specificity of the additional identified signals. One of our potential causal variants influences AS3MT expression and nearby DNA methylation in numerous GTEx tissue types (with rs4919690 as a likely causal variant). Several SNPs in our confidence sets overlap transcription factor binding sites and cis-regulatory elements (from ENCODE). Taken together, our analyses reveal multiple potential causal variants in the 10q24.32 region influencing AME, including a variant shared across populations, and elucidate potential biological mechanisms underlying the impact of genetic variation on AME.

Two Distinct Deleterious Causative Variants in a Family with Multiple Cancer-Affected Patients.

Only 5 to 10% of cancers are hereditary, but they are particularly important since they can be passed down from generation to generation, and family members are at elevated risk. Although screening methods are one of the essential strategies for dealing with hereditary cancers, they do not have high specificity and sensitivity. The emergence of whole-exome sequencing (WES) causes a significant increase in the diagnostic rate of cancer-causing variants in at-risk families. We performed WES on the proband's DNA sample from an Iranian family with multiple cancer-affected members to identify potential causative variants. Multiple in silico tools were used to evaluate the candidate variants' pathogenicity and their effects on the protein's structure, function, and stability. Moreover, the candidate variants were co-segregated in the family with Sanger sequencing. The WES data analysis identified two pathogenic variants (CHEK2: NM_007194.4: c.538C>T, p.Arg180Cys and MLH1: NM_000249.4, c.844G>A, p.Ala282Thr). Sanger sequencing data showed each of the variants was incompletely segregated with phenotype, but both of them explained the patient's phenotype together. Also, the structural analysis demonstrated that due to the variant (c.538C>T), a salt bridge between arginine 180 and glutamic acid 149 was lost. Indeed, several protein stability tools described both variants as destabilizing. Herein, we interestingly identify two distinct deleterious causative variants (CHEK2: NM_007194.4: c.538C>T, p.Arg180Cys and MLH1: NM_000249.4, c.844G>A, p.Ala282Thr) in a family with several cancer-affected members. Furthermore, this study's findings established the utility of WES in the genetic diagnostics of cancer.

GWAS of Genetic Resilience to Age‐Related Risk for Alzheimer’s Disease in Admixed African Ancestry Individuals

AbstractBackgroundWe sought to identify genetic variants that confer resilience to age‐related Alzheimer’s disease (AD) risk in subjects of admixed African ancestry.MethodsParticipants were ages 60+, of African ancestry (> = 25%), and diagnosed as cases or controls. Genetic data were available from SNP arrays imputed to TOPMed, whole genome sequencing (WGS), or whole exome sequencing (WES). Genome‐wide association studies (GWAS) were performed per data type, split by Hispanic/non‐Hispanic participants (Table‐1), followed by meta‐analysis (Plink v2.0; GWAMA v2.2.2). GWAS performed multiple linear regression on an AD‐age score that models resilience to age‐related risk for AD (Le Guen & Belloy et al. 2021; Figure‐1). Models adjusted for sex, APOE*4/APOE*2 dosage, the first five genetic principal components (PC‐AiR; GENESIS; R v3.6), and array/sequencing center. Brain amyloid Positron Emission Tomography (PET) qualitative reads together with SNP array data imputed to TOPMed were also available for African American individuals in A4. These data were independent from the primary GWAS and were used to evaluate significant GWAS variants for their association with binarized amyloid outcome (negative/positive), using a consistent model that additionally adjusted for age‐at‐scan.ResultsWe found a novel protective genome‐wide significant intergenic variant (rs77450754) ∼200kb downstream of ATXN8OS/KLHL1 (Figure‐2‐3). It has a MAF of 13% in African ancestry samples, but is rare (0.1%) in Europeans and (<2%) Latino/Admixed Americans in gnomAD v.3. The latter may explain why effects appeared less pronounced in the considered Hispanic samples (since linkage to potential causal nearby variants may differ) (Figure‐3). The variant also displayed a protective, non‐significant, association with brain amyloid in African ancestry samples (Table‐2). Interestingly, KLHL1, also known as MRP2, is a member of the superfamily of ATP‐binding cassette (ABC) transporters, and has potential functions related to neurite outgrowth and efflux transport at the blood‐brain barrier, supporting its potential relevance to AD pathogenesis.ConclusionsOur results provide new insights into genetic resilience to AD across aging and emphasize the importance of including ancestrally‐diverse populations in genetic studies. Our results also highlight the need for the field to compose ancestrally‐diverse multi‐omics and endophenotype data to support further validation analyses of genetic risk loci prioritized from GWAS.

Extensive germline-somatic interplay contributes to prostate cancer progression through HNF1B co-option of TMPRSS2-ERG

Genome-wide association studies have identified 270 loci conferring risk for prostate cancer (PCa), yet the underlying biology and clinical impact remain to be investigated. Here we observe an enrichment of transcription factor genes including HNF1B within PCa risk-associated regions. While focused on the 17q12/HNF1B locus, we find a strong eQTL for HNF1B and multiple potential causal variants involved in the regulation of HNF1B expression in PCa. An unbiased genome-wide co-expression analysis reveals PCa-specific somatic TMPRSS2-ERG fusion as a transcriptional mediator of this locus and the HNF1B eQTL signal is ERG fusion status dependent. We investigate the role of HNF1B and find its involvement in several pathways related to cell cycle progression and PCa severity. Furthermore, HNF1B interacts with TMPRSS2-ERG to co-occupy large proportion of genomic regions with a remarkable enrichment of additional PCa risk alleles. We finally show that HNF1B co-opts ERG fusion to mediate mechanistic and biological effects of the PCa risk-associated locus 17p13.3/VPS53/FAM57A/GEMIN4. Taken together, we report an extensive germline-somatic interaction between TMPRSS2-ERG fusion and genetic variations underpinning PCa risk association and progression.

Identification of nine novel variants across PAX3, SOX10, EDNRB, and MITF genes in Waardenburg syndrome with next-generation sequencing.

Waardenburg syndrome (WS) is a hereditary, genetically heterogeneous disorder characterized by variable presentations of sensorineural hearing impairment and pigmentation anomalies. This study aimed to investigate the clinical features of WS in detail and determine the genetic causes of patients with clinically suspected WS. A total of 24 patients from 21 Han-Taiwanese families were enrolled and underwent comprehensive physical and audiological examinations. We applied targeted next-generation sequencing (NGS) to investigate the potential causative variants in these patients and further validated the candidate variants through Sanger sequencing. We identified 19 causative variants of WS in our cohort. Of these variants, nine were novel and discovered in PAX3, SOX10, EDNRB, and MITF genes, including missense, nonsense, deletion, and splice site variants. Several patients presented with skeletal deformities, hypotonia, megacolon, and neurological disorders that were rarely seen in WS. This study revealed highly phenotypic variability in Taiwanese WS patients and demonstrated that targeted NGS allowed us to clarify the genetic diagnosis and extend the genetic variant spectrum of WS.

Construction of a transposase accessible chromatin landscape reveals chromatin state of repeat elements and potential causal variant for complex traits in pigs

BackgroundA comprehensive landscape of chromatin states for multiple mammalian tissues is essential for elucidating the molecular mechanism underlying regulatory variants on complex traits. However, the genome-wide chromatin accessibility has been only reported in limited tissue types in pigs.ResultsHere we report a genome-wide landscape of chromatin accessibility of 20 tissues in two female pigs at ages of 6 months using ATAC-seq, and identified 557,273 merged peaks, which greatly expanded the pig regulatory element repository. We revealed tissue-specific regulatory elements which were associated with tissue-relevant biological functions. We identified both positive and negative significant correlations between the regulatory elements and gene transcripts, which showed distinct distributions in terms of their strength and distances from corresponding genes. We investigated the presence of transposable elements (TEs) in open chromatin regions across all tissues, these included identifications of porcine endogenous retroviruses (PERVs) exhibiting high accessibility in liver and homology of porcine specific virus sequences to universally accessible transposable elements. Furthermore, we prioritized a potential causal variant for polyunsaturated fatty acid in the muscle.ConclusionsOur data provides a novel multi-tissues accessible chromatin landscape that serve as an important resource for interpreting regulatory sequences in tissue-specific and conserved biological functions, as well as regulatory variants of loci associated with complex traits in pigs.

Whole-Exome Sequencing identified Olfactory Receptor genes as a key contributor to extreme obesity with progression to nonalcoholic steatohepatitis in Mexican patients: Olfactory receptor genes in obese NASH patients

Obesity is a global health problem that triggers fat liver accumulation. The prevalence of obesity and the risk of non-alcoholic steatohepatitis (NASH) among young obese Mexican is high. Furthermore, genetic predisposition is a key factor in weight gain and disrupts metabolism. Herein, we used Whole-Exome Sequencing to identify potential causal variants and the biological processes that lead to obesity with progression to NASH among Mexican patients. Whole-Exome Sequencing was performed in nine obese patients with NASH diagnosis with a BMI ≥30 kg/m2 and one control (BMI=24.2 kg/m2) by using the Ion S5TM platform. Genetic variants were determined by Ion Reporter software. Enriched GO biological set genes were identified by the WebGestalt tool. Genetic variants within ≥2 obese NASH patients and having scores of SIFT 0.0-0.05 and Polyphen 0.85-1.0 were categorized as pathogenic. A total of 1359 variants with a probable pathogenic effect were determined in obese patients with NASH diagnosis. After several filtering steps, the most frequent pathogenic variants found were rs25640-HSD17B4, rs8105737-OR1I1, rs998544-OR5R1, and rs4916685, rs10037067, and rs2366926 in ADGRV1. Notably, the primary biological processes affected by these pathogenic variants were the sensory perception and detection of chemical stimulus pathways in which the olfactory receptor gene family was the most enriched. Variants in the olfactory receptor genes were highly enriched in Mexican obese patients that progress to NASH and could be potential targets of association studies.

Genetic Spectrum and Clinical Heterogeneity of Chinese Frontotemporal Dementia Patients: Data from PUMCH Dementia Cohort.

Background:There are relatively few data on the genetic spectrum of Chinese frontotemporal dementia (FTD) population.Objective:With the dementia cohort of Peking Union Medical College Hospital, we aim to illustrate the genetic spectrum of FTD patients, as well as the phenotypic heterogeneity of FTD-gene variant carriers.Methods:204 unrelated, clinically diagnosed FTD patients of Chinese ancestry were enrolled. All the participants received demographic survey, history inquiry, physical examination, cognitive assessment, blood biochemical test, brain CT/MRI, and gene sequencing.Results:56.4% (115/204) participants were clinically diagnosed with behavioral variant of FTD, 20.6% (42/204) with nonfluent/agrammatic variant primary progressive aphasia (PPA), 20.1% (41/204) with semantic variant PPA, and 2.9% (6/204) with mixed variant PPA. 11.8% (24/204) subjects harbored the potential causative variants in FTD-related genes, including the MAPT (n = 7), TBK1 (n = 7), GRN (n = 2), TBK1+GRN (n = 1), VCP (n = 1), TARDBP (n = 1), UBQLN2 (n = 1), SQSTM1 (n = 1), DCTN1 (n = 1), HNRNPA1 (n = 1), and C9orf72 GGGGCC repeats (n = 1). The TBK1 T31fs, T457fs, K622fs, c.359-1G>A, the VCP P188T, and the GRN P50fs, P439fs were novel pathogenic/likely pathogenic variants. The TBK1 carriers showed a later disease onset and a higher incidence of parietal atrophy relative to the MAPTcarriers.Conclusion:There is genetic and clinical heterogeneity among Chinese FTD population. The TBK1 has a high mutation frequency in Chinese FTD patients.

A scalable Bayesian functional GWAS method accounting for multivariate quantitative functional annotations with applications for studying Alzheimer disease.

SummaryExisting methods for integrating functional annotations in genome-wide association studies (GWASs) to fine-map and prioritize potential causal variants are limited to using non-overlapped categorical annotations or limited by the computation burden of modeling genome-wide variants. To overcome these limitations, we propose a scalable Bayesian functional GWAS method to account for multivariate quantitative functional annotations (BFGWAS_QUANT), accompanied by a scalable computation algorithm enabling joint modeling of genome-wide variants. Simulation studies validated the performance of BFGWAS_QUANT for accurately quantifying annotation enrichment and improving GWAS power. Applying BFGWAS_QUANT to study five Alzheimer disease (AD)-related phenotypes using individual-level GWAS data (n = ∼1,000), we found that histone modification annotations have higher enrichment than expression quantitative trait locus (eQTL) annotations for all considered phenotypes, with the highest enrichment in H3K27me3 (polycomb regression). We also found that cis-eQTLs in microglia had higher enrichment than eQTLs of bulk brain frontal cortex tissue for all considered phenotypes. A similar enrichment pattern was also identified using the International Genomics of Alzheimer’s Project (IGAP) summary-level GWAS data of AD (n = ∼54,000). The strongest known APOE E4 risk allele was identified for all five phenotypes, and the APOE locus was validated using the IGAP data. BFGWAS_QUANT fine-mapped 32 significant variants from 1,073 genome-wide significant variants in the IGAP data. We also demonstrated that the polygenic risk scores (PRSs) using effect size estimates by BFGWAS_QUANT had a similar prediction accuracy as other methods assuming a sparse causal model. Overall, BFGWAS_QUANT is a useful GWAS tool for quantifying annotation enrichment and prioritizing potential causal variants.

Targeted Next Generation Sequencing Revealed Novel Variants in the PKD1 and PKD2 Genes of Iranian Patients with Autosomal Dominant Polycystic Kidney Disease.

Autosomal dominant polycystic kidney disease (ADPKD), one of the common inherited disorders in humans, is characterized by the development and enlargement of renal cysts, often leading to end-stage renal disease (ESRD). In this study, Iranian ADPKD families were subjected to high-throughput DNA sequencing to find potential causative variants facilitating the way toward risk assessment and targeted therapy. Our protocol was based on the targeted next generation sequencing (NGS) panel previously developed in our center comprising 12 genes involved in PKD. This panel has been applied to investigate the genetic causes of 32 patients with a clinical suspicion of ADPKD. We identified a total of 31 variants for 32 individuals, two of which were each detected in two individuals. Twenty-seven out of 31 detected variants were interpreted as pathogenic/likely pathogenic and the remaining 4 of uncertain significance with a molecular diagnostic success rate of 87.5%. Among these variants, 25 PKD1/2 pathogenic/likely pathogenic variants were detected in 32 index patients (78.1%), and variants of uncertain significance in four individuals (12.5% in PKD1/2). The majority of variants was identified in PKD1 (74.2%). Autosomal recessive PKD was identified in one patient, indicating the similarities between recessive and dominant PKD. In concordance with earlier studies, this biallelic PKD1 variant, p.Arg3277Cys, leads to rapidly progressive and severe disease with very early-onset ADPKD. Our findings suggest that targeted gene panel sequencing is expected to be the method of choice to improve diagnostic and prognostic accuracy in PKD patients with heterogeneity in genetic background.

Abstract P3109: Identification Of Drug Transporter Genomic Variants And Inhibitors That Protect Against Doxorubicin Induced Cardiotoxicity

Multiple pharmacogenomic studies have identified the synonymous genomic variant rs7853758 (G > A, L461L) and the intronic variant rs885004 in SLC28A3 (solute carrier family 28 member 3) as statistically associated with a lower incidence of anthracycline-induced cardiotoxicity. However, the true causal variant(s), the cardioprotective mechanism of this locus, the role of SLC28A3 and other solute carrier (SLC) transporters in anthracycline-induced cardiotoxicity, and the suitability of SLC transporters as targets for cardioprotective drugs has not been investigated. To answer these questions, six well-phenotyped, doxorubicin-treated pediatric patients from the original association study cohort were recruited again, and human-induced pluripotent stem cell-derived cardiomyocytes were generated. Patient-specific doxorubicin-induced cardiotoxicity (DIC) was then characterized using assays of cell viability, activated caspase 3/7, and doxorubicin uptake. We demonstrate that patient-derived cardiomyocytes recapitulate the cardioprotective effect of the SLC28A3 locus and that SLC28A3 expression influences the severity of DIC. The role of SLC28A3 in DIC was then confirmed using overexpression and knockout of SLC28A3 in isogenic human-induced pluripotent stem cell-derived cardiomyocytes. To identify the potential causal variant, Nanopore-based fine-mapping of the SLC28A3 locus was then completed after SLC28A3 resequencing. An extended in silico haplotype and functional analysis followed by base editing, identified a novel cardioprotective single nucleotide polymorphism, rs11140490, in the SLC28A3 locus. Next, using high-throughput drug screening of 1200 drugs in patient-derived cardiomyocytes and whole organism validation in mice, we identify the SLC competitive inhibitor desipramine as protective against DIC. Finally, using eight cancer cell lines, we demonstrate that desipramine does not impede doxorubicin cytotoxicity.

Investigation of Genetic Causes in Patients with Congenital Heart Disease in Qatar: Findings from the Sidra Cardiac Registry.

Congenital heart disease (CHD) is one of the most common forms of birth defects worldwide, with a prevalence of 1–2% in newborns. CHD is a multifactorial disease partially caused by genetic defects, including chromosomal abnormalities and single gene mutations. Here, we describe the Sidra Cardiac Registry, which includes 52 families and a total of 178 individuals, and investigate the genetic etiology of CHD in Qatar. We reviewed the results of genetic tests conducted in patients as part of their clinical evaluation, including chromosomal testing. We also performed whole exome sequencing (WES) to identify potential causative variants. Sixteen patients with CHD had chromosomal abnormalities that explained their complex CHD phenotype, including six patients with trisomy 21. Moreover, using exome analysis, we identified potential CHD variants in 24 patients, revealing 65 potential variants in 56 genes. Four variants were classified as pathogenic/likely pathogenic based on the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) classification; these variants were detected in four patients. This study sheds light on several potential genetic variants contributing to the development of CHD. Additional functional studies are needed to better understand the role of the identified variants in the pathogenesis of CHD.

Functional significance of the rare rs35667974 IFIH1 gene polymorphism, associated with multiple autoimmune diseases, using a structural biological approach

Autoimmune diseases, which affect approximately 5% of human population, are a range of diseases in which the immune response to self-antigens results in damage or dysfunction of tissues. Recent genome wide association studies (GWAS) have successfully identified novel autoimmune disease-associated loci, with many of them shared by multiple disease-associated pathways but much of the genetics and pathophysiological mechanisms remain still obscure. Considering that most of the potential causal variants are still unknown, many studies showed that the missense variant rs35667974 at interferon-induced with helicase C domain 1 (IFIH1) gene is protective for type 1 diabetes (T1D), psoriasis (PS) and psoriatic arthritis (PsA). Recently, this variant was found to be also associated with ankylosing spondylitis (AS), Crohn’s disease (CD) and ulcerative colitis (UC). The IFIH1 gene encodes a cytoplasmic RNA helicase otherwise known as melanoma differentiation-associated 5 (MDA5) that recognizes viral RNA and is involved in innate immunity through recognition of viral RNA. In the present study we sought to investigate the association of the rare rs35667974 variant of IFIH1 gene, which resides in exon 14 and changes a conserved isoleucine at position #923 to valine, in the development of various autoimmune diseases and give a reason for the selectivity affecting different autoimmune diseases. Evolutionary studies and three-dimensional (3 D) homology modelling were employed on the MDA5 protein product, through its association with dsRNA, recognition factor controlling cytokine and chemokine signalling, to investigate the protective role of the MDA5 variant for certain autoimmune diseases.

Exploring Potential Causal Genes for Uterine Leiomyomas: A Summary Data-Based Mendelian Randomization and FUMA Analysis

Objective: To explore potential causal genetic variants and genes underlying the pathogenesis of uterine leiomyomas (ULs). Methods: We conducted the summary data-based Mendelian randomization (SMR) analyses and performed functional mapping and annotation using FUMA to examine genetic variants and genes that are potentially involved in the pathogenies of ULs. Both analyses used summarized data of a recent genome-wide association study (GWAS) on ULs, which has a total sample size of 244,324 (20,406 cases and 223,918 controls). We performed separate SMR analysis using CAGE and GTEx eQTL data. Results: Using the CAGE eQTL data, our SMR analysis identified 13 probes tagging 10 unique genes that were pleiotropically/potentially causally associated with ULs, with the top three probes being ILMN_1675156 (tagging CDC42, PSMR = 8.03 × 10−9), ILMN_1705330 (tagging CDC42, PSMR = 1.02 × 10−7) and ILMN_2343048 (tagging ABCB9, PSMR = 9.37 × 10−7). Using GTEx eQTL data, our SMR analysis did not identify any significant genes after correction for multiple testing. FUMA analysis identified 106 independent SNPs, 24 genomic loci and 137 genes that are potentially involved in the pathogenesis of ULs, seven of which were also identified by the SMR analysis. Conclusions: We identified many genetic variants, genes, and genomic loci that are potentially involved in the pathogenesis of ULs. More studies are needed to explore the exact underlying mechanisms in the etiology of ULs.

A unique Smith-Magenis patient with a de novo intragenic deletion on the maternally inherited overexpressed RAI1 allele

RAI1 is a dosage-sensitive gene whose decreased or increased expression by recurrent and non-recurrent 17p11.2 deletions or duplications causes Smith-Magenis (SMS) or Potocki-Lupski syndromes (PTLS), respectively. Here we report on a 21-year-old female patient showing SMS phenotype who was found to carry a 3.4 kb de novo intragenic RAI1 deletion. Interestingly, a significant increase in RAI1 transcript levels was identified in the patient’s, brother’s and mother’s peripheral blood cells. Allele-specific dosage analysis revealed that the patient’s maternally inherited overexpressed RAI1 allele harbors the intragenic deletion, confirming the SMS diagnosis due to the presence of a single wild-type RAI1 functional allele. The mother and brother do not present any PTLS neurologic/behavioral clinical features. Extensive sequencing of RAI1 promoter and predicted regulatory regions showed no potential causative variants accounting for gene overexpression. However, the mother and both children share a novel private missense variant in RAI1 exon 3, currently classified as a VUS (uncertain significance), though predicted by two bioinformatic tools to disrupt the binding site of one specific transcription factor. The reported familial case, the second showing RAI1 overexpression in the absence of RAI1 duplication, may help to understand the regulation of RAI1 dosage sensitivity although its phenotypic effect remains to be determined.

Identification of a Novel Functional Non-synonymous Single Nucleotide Polymorphism in Frizzled Class Receptor 6 Gene for Involvement in Depressive Symptoms.

Although numerous susceptibility loci for depression have been identified in recent years, their biological function and molecular mechanism remain largely unknown. By using an exome-wide association study for depressive symptoms assessed by the Center for Epidemiological Studies Depression (CES-D) score, we discovered a novel missense single nucleotide polymorphism (SNP), rs61753730 (Q152E), located in the fourth exon of the frizzled class receptor 6 gene (FZD6), which is a potential causal variant and is significantly associated with the CES-D score. Computer-based in silico analysis revealed that the protein configuration and stability, as well as the secondary structure of FZD6 differed greatly between the wild-type (WT) and Q152E mutant. We further found that rs61753730 significantly affected the luciferase activity and expression of FZD6 in an allele-specific way. Finally, we generated Fzd6-knockin (Fzd6-KI) mice with rs61753730 mutation using the CRISPR/Cas9 genome editing system and found that these mice presented greater immobility in the forced swimming test, less preference for sucrose in the sucrose preference test, as well as decreased center entries, center time, and distance traveled in the open filed test compared with WT mice after exposed to chronic social defeat stress. These results indicate the involvement of rs61753730 in depression. Taken together, our findings demonstrate that SNP rs61753730 is a novel functional variant and plays an important role in depressive symptoms.

Flashfm-ivis: interactive visualization for fine-mapping of multiple quantitative traits.

Summary flashfm-ivis provides a suite of interactive visualization plots to view potential causal genetic variants that underlie associations that are shared or distinct between multiple quantitative traits and compares results between single- and multi-trait fine-mapping. Unique features include network diagrams that show joint effects between variants for each trait and regional association plots that integrate fine-mapping results, all with user-controlled zoom features for an interactive exploration of potential causal variants across traits.Availability and implementation flashfm-ivis is an open-source software under the MIT license. It is available as an interactive web-based tool (http://shiny.mrc-bsu.cam.ac.uk/apps/flashfm-ivis/) and as an R package. Code and documentation are available at https://github.com/fz-cambridge/flashfm-ivis and https://zenodo.org/record/6376244#.YjnarC-l2X0. Additional features can be downloaded as standalone R libraries to encourage reuse.Supplementary information Supplementary information are available at Bioinformatics online.