Causal Variants Research Articles

A novel pathogenic variant in the KCTD7 gene in a patient with neuronal ceroid lipofuscinosis (CLN14): a case report and review of the literature

BackgroundNeuronal ceroid lipofuscinosis (NCL) is a heterogeneous group of 13 rare, progressive neurodegenerative diseases of the brain and retina. CLN14 is a very rare subtype of NCL caused by pathogenic variants in the KCTD7 gene. Only four cases of this subtype have been reported in the literature.Case presentationA nine-month-old, previously healthy male who was firstborn to first-cousin parents presented with progressive psychomotor regression, dysmorphic facial features, myoclonus, and vision loss. Neurological examination showed generalized hypotonia and brisk reflexes. He continued to deteriorate until age 18 months, when he developed his first generalized tonic-clonic seizure. An ophthalmological examination showed a hypopigmented fundus and slight temporal disc pallor. Brain MRI showed mild generalized brain atrophy and white matter disease. EEG revealed a severely abnormal trace marked by generalized, high amplitude, sharply contoured, polymorphic delta slowing intermixed with theta slowing and some alpha activity, with disorganized and scattered spikes and sharp waves. The patient continued to have uncontrolled seizures and further psychomotor regression until he died of status epilepticus and pneumonia at the age of 44 months. WES identified a novel homozygous variant c.413T > C, p.(Leu138Pro) in the KCTD7 gene, causing an amino acid transition from leucine to proline at position 138. Both parents were carriers of the same variant.ConclusionsWe present the fifth known case of CLN14 in the literature and report the clinical course and a novel underlying likely causative variant in the KCTD7 gene. The improving accessibility and affordability of genetic testing will likely uncover more NCL cases and further expand the disease’s genotypic and phenotypic spectrum.

Гено-фенотипические особенности муковисцидоза у российских детей из Чеченской и Карачаево-Черкесской Республик

Background: Cystic fibrosis (CF) is an autosomal recessive disease that occurs with a frequency of 1:1,500 to 1:5,000 newborns, according to the World Health Organization. In the Russian Federation, the average incidence of the disease is 1:10,000 newborns. The medical and social significance of this disease is associated with the early disability of patients, the need for long-term treatment and constant follow-up, as well as the heterogeneity of phenotypic manifestations, which in turn requires early diagnosis using molecular genetic methods. The aim of the study: To study the clinical and molecular genetic characteristics of children with CF from the Chechen and Karachay-Cherkess Republics. Materials and methods: The study included 237 patients with a confirmed diagnosis of CF. Molecular genetic diagnostics was performed using the method of mass parallel sequencing using a hybridization target panel that includes the entire CFTR gene. Sequencing was performed on the MiSeq platform. All causal nucleotide variants were validated using Sanger sequencing. Results: The most common amino acid variants of the CFTR gene in patients with CF from the Chechen Republic are: p.Y515* (94 alleles/78.3%), p.E92K (18 alleles/15%). The presence of variants p.Y515* and p.E92K in the heterozygous state in the genotype is more often a favorable prognostic factor for the absence of pancreatic lesion. The main distinctive clinical feature of these patients is the onset of the disease with manifestations of pseudo-Bartter syndrome. The major amino acid variant of the CFTR gene in patients with CF from the Karachay-Cherkess Republic is: p.W1282* (28 alleles/70%). The number of patients with pseudo-Bartter syndrome was 54.5%. In 90% of cases, in patients who are homo- and heterozygous according to p.W1282*, a severe degree of pancreatic insufficiency was noted. Conclusion: Patients with cystic fibrosis from the Chechen and Karachay-Cherkess Republics are a difficult category of patients, despite neonatal screening, and cystic fibrosis therapy developed and implemented. This is partly due to the peculiarities of the phenotype and the unique distribution of alleles and genotypes of the CFTR gene. The high number of homozygous variants of the CFTR gene in patients from the Chechen and Karachay-Cherkess Republics is probably due to monoethnic marital assortativity.

A Splice Site Variant in ADAMTS3 Is the Likely Causal Variant for Pulmonary Hypoplasia with Anasarca in Persian/Persian-Cross Sheep.

Pulmonary hypoplasia with anasarca, or hydrops fetalis, is characterized by stillbirth, diffuse oedema, and generalized lymph node hypoplasia. The enlarged fetus frequently causes dystocia. The disease has been reported in cattle and sheep as an inherited condition with a recessive mode of inheritance. This is the first report of the disease in Persian/Persian-cross sheep in Australia. Affected fetuses were reported from three flocks, and a total of eleven affected, eleven obligate carrier, and 188 related Persian/Persian-cross animals were available for analysis, as well as unrelated control animals. SNP genotyping revealed a region of homozygosity in affected animals on ovine chromosome six, which contained the functional candidate gene ADAMTS3. Whole genome sequencing of two affected fetuses and one obligate carrier ewe revealed a single nucleotide deletion, ENSOARG00000013204:g.87124344delC, located 3 bp downstream from a donor splice site region in the ADAMTS3 gene. Sanger sequencing of cDNA containing this variant further revealed that it is likely to introduce an early splice site in exon 14, resulting in a loss of 6 amino acids at the junction of exon 14 and intron 14/15. A genotyping assay was developed, and the ENSOARG00000013204:g.87124344delC segregated with disease in 209 animals, allowing for effective identification of carrier animals.

Identification of candidate causal variants and target genes at 41 breast cancer risk loci through differential allelic expression analysis

Understanding breast cancer genetic risk relies on identifying causal variants and candidate target genes in risk loci identified by genome-wide association studies (GWAS), which remains challenging. Since most loci fall in active gene regulatory regions, we developed a novel approach facilitated by pinpointing the variants with greater regulatory potential in the disease’s tissue of origin. Through genome-wide differential allelic expression (DAE) analysis, using microarray data from 64 normal breast tissue samples, we mapped the variants associated with DAE (daeQTLs). Then, we intersected these with GWAS data to reveal candidate risk regulatory variants and analysed their cis-acting regulatory potential. Finally, we validated our approach by extensive functional analysis of the 5q14.1 breast cancer risk locus. We observed widespread gene expression regulation by cis-acting variants in breast tissue, with 65% of coding and noncoding expressed genes displaying DAE (daeGenes). We identified over 54 K daeQTLs for 6761 (26%) daeGenes, including 385 daeGenes harbouring variants previously associated with BC risk. We found 1431 daeQTLs mapped to 93 different loci in strong linkage disequilibrium with risk-associated variants (risk-daeQTLs), suggesting a link between risk-causing variants and cis-regulation. There were 122 risk-daeQTL with stronger cis-acting potential in active regulatory regions with protein binding evidence. These variants mapped to 41 risk loci, of which 29 had no previous report of target genes and were candidates for regulating the expression levels of 65 genes. As validation, we identified and functionally characterised five candidate causal variants at the 5q14.1 risk locus targeting the ATG10 and ATP6AP1L genes, likely acting via modulation of alternative transcription and transcription factor binding. Our study demonstrates the power of DAE analysis and daeQTL mapping to identify causal regulatory variants and target genes at breast cancer risk loci, including those with complex regulatory landscapes. It additionally provides a genome-wide resource of variants associated with DAE for future functional studies.

Phenotype-driven genomics enhance diagnosis in children with unresolved neuromuscular diseases.

Establishing a molecular diagnosis remains challenging in half of individuals with childhood-onset neuromuscular diseases (NMDs) despite exome sequencing. This study evaluates the diagnostic utility of combining genomic approaches in undiagnosed NMD patients. We performed deep phenotyping of 58 individuals with unsolved childhood-onset NMDs that have previously undergone inconclusive exome studies. Genomic approaches included trio genome sequencing and RNASeq. Genetic diagnoses were reached in 23 out of 58 individuals (40%). Twenty-one individuals carried causal single nucleotide variants (SNVs) or small insertions and deletions, while 2 carried pathogenic structural variants (SVs). Genomic sequencing identified pathogenic variants in coding regions or at the splice site in 17 out of 21 resolved cases, while RNA sequencing was additionally required for the diagnosis of 4 cases. Reasons for previous diagnostic failures included low coverage in exonic regions harboring the second pathogenic variant and involvement of genes that were not yet linked to human diseases at the time of the first NGS analysis. In summary, our systematic genetic analysis, integrating deep phenotyping, trio genome sequencing and RNASeq, proved effective in diagnosing unsolved childhood-onset NMDs. This approach holds promise for similar cohorts, offering potential improvements in diagnostic rates and clinical management of individuals with NMDs.

Proteome-wide mendelian randomization identifies novel therapeutic targets for chronic kidney disease.

There is an urgent need to pinpoint novel targets for drug discovery in the context of chronic kidney disease (CKD), and the proteome represents a significant pool of potential therapeutic targets. To address this, we performed proteome-wide analyses using Mendelian randomization (MR) and colocalization techniques to uncover potential targets for CKD. We extracted summary-level data from the ARIC study, focusing on 7213 European American (EA) individuals and 4657 plasma proteins. To broaden our analysis, we incorporated genetic association data from Icelandic cohorts, thereby enhancing our investigation into the correlations with chronic kidney disease (CKD), creatinine-based estimated glomerular filtration rate (eGFRcrea), and estimated glomerular filtration rate (eGFR). We utilized genetic association data from the GWAS Catalog, including CKD (765,348, 625,219 European ancestry and 140,129 non-European ancestry), eGFRcrea (1,004,040, European ancestry), and eGFR (567,460, European ancestry). Employing MR analysis, we estimated the associations between proteins and CKD risk. Additionally, we conducted colocalization analysis to evaluate the existence of shared causal variants between the identified proteins and CKD. We detected notable correlations between levels predicted based on genetics of three circulating proteins and CKD, eGFRcrea, and eGFR. Notably, our colocalization analysis provided robust evidence supporting these associations. Specifically, genetically predicted levels of Transcription elongation factor A protein 2 (TCEA2) and Neuregulin-4 (NRG4) exhibited an inverse relationship with CKD risk, while Glucokinase regulatory protein (GCKR) showed an increased risk of CKD. Furthermore, our colocalization analysis also supported the associations of TCEA2, NRG4, and GCKR with the risk of eGFRcrea and eGFR.

EDA Variants Are Responsible for Approximately 90% of Deciduous Tooth Agenesis.

Deciduous tooth agenesis is a severe craniofacial developmental defect because it affects masticatory function from infancy and may result in delayed growth and development. Here, we aimed to identify the crucial pathogenic genes and clinical features of patients with deciduous tooth agenesis. We recruited 84 patients with severe deciduous tooth agenesis. Whole-exome and Sanger sequencing were used to identify the causative variants. Phenotype-genotype correlation analysis was conducted. We identified 54 different variants in 8 genes in 84 patients, including EDA (73, 86.9%), PAX9 (2, 2.4%), LRP6 (2, 2.4%), MSX1 (2, 2.4%), BMP4 (1, 1.2%), WNT10A (1, 1.2%), PITX2 (1, 1.2%), and EDARADD (1, 1.2%). Variants in ectodysplasin A (EDA) accounted for 86.9% of patients with deciduous tooth agenesis. Patients with the EDA variants had an average of 15.4 missing deciduous teeth. Mandibular deciduous central incisors had the highest missing rate (100%), followed by maxillary deciduous lateral incisors (98.8%) and mandibular deciduous lateral incisors (97.7%). Our results indicated that EDA gene variants are major pathogenic factors for deciduous tooth agenesis, and EDA is specifically required for deciduous tooth development. The results provide guidance for clinical diagnosis and genetic counseling of deciduous tooth agenesis.

Characterization of non-coding variants associated with transcription factor binding through ATAC-seq-defined footprint QTLs in liver.

Non-coding variants discovered by genome-wide association studies (GWAS) are enriched in regulatory elements harboring transcription factor (TF) binding motifs, strongly suggesting a connection between disease association and the disruption of cis-regulatory sequences. Occupancy of a TF inside a region of open chromatin can be detected in ATAC-seq where bound TFs block the transposase Tn5, leaving a pattern of relatively depleted Tn5 insertions known as a "footprint". Here, we sought to identify variants associated with TF-binding, or "footprint quantitative trait loci" (fpQTLs) in ATAC-seq data generated from 170 human liver samples. We used computational tools to scan the ATAC-seq reads to quantify TF binding likelihood as "footprint scores" at variants derived from whole genome sequencing generated in the same samples. We tested for association between genotype and footprint score and observed 693 fpQTLs associated with footprint-inferred TF binding (FDR < 5%). Given that Tn5 insertion sites are measured with base-pair resolution, we show that fpQTLs can aid GWAS and QTL fine-mapping by precisely pinpointing TF activity within broad trait-associated loci where the underlying causal variant is unknown. Liver fpQTLs were strongly enriched across ChIP-seq peaks, liver expression QTLs (eQTLs), and liver-related GWAS loci, and their inferred effect on TF binding was concordant with their effect on underlying sequence motifs in 80% of cases. We conclude that fpQTLs can reveal causal GWAS variants, define the role of TF binding site disruption in disease and provide functional insights into non-coding variants, ultimately informing novel treatments for common diseases.

Transcriptomic Insights into the Atrial Fibrillation Susceptibility Locus near the MYOZ1 and SYNPO2L Genes

Genome-wide association studies have identified a locus on chromosome 10q22, where many co-inherited single nucleotide polymorphisms (SNPs) are associated with atrial fibrillation (AF). This study seeks to identify the impact of this locus on gene expression at the transcript isoform level in human left atria and to gain insight into potential causal variants. Bulk RNA sequencing was analyzed to identify myozenin 1 (MYOZ1) and synaptopodin 2-like (SYNPO2L) transcript isoforms and the association of common SNPs in this region with transcript isoform expression levels. Chromatin marks were used to suggest candidate regulatory SNPs in this region. Protein amino acid changes were examined for predicted functional consequences. Transfection of MYOZ1 and two SYNPO2L isoforms were performed to localize their encoded proteins in cardiomyocytes derived from stem cells. We identified one MYOZ1 transcript isoform and four SYNPO2L transcript isoforms, two of which encode proteins, while the other two encode long noncoding RNAs (lncRNAs). The risk allele of the strongest AF susceptibility SNP on chromosome 10q22 is associated with decreased MYOZ1 expression and increased expression of the two SNYPO2L lncRNA isoforms. There are many SNPs co-inherited with the top AF-associated SNP due to linkage disequilibrium (LD), including rs11000728, which we propose as the MYOZ1 regulatory SNP, confirmed by reporter gene transfection. In addition, this LD block includes three missense SNPs in the SYNPO2L gene, with the minor protective haplotype predicted to be detrimental to protein function. MYOZ1 and both protein isoforms of SYNPO2L were localized to the sarcomere. This is a complex locus with the potential for several SNPs in a haplotype to alter AF susceptibility by opposing effects on MYOZ1 and SYNPO2L lncRNA expression, along with effects on SYNPO2L protein function.

Causality of genetically determined serum metabolites on lower back pain or/and sciatica: a comprehensive Mendelian randomized study.

There is an urgent need to confirm biomarkers reflecting the pathogenesis and targeted drugs of lower back pain or/and sciatica in clinical practice. This study aimed to conduct a two sample bidirectional Mendelian randomization (MR) analysis to explore the causal link between 486 serum metabolites and lower back pain or/and sciatica. All data come from two public shared databases of European ancestry and single nucleotide polymorphisms (SNPs) for lower back pain or/and sciatica acted as instrumental variables. The traditional inverse variance weighting (IVW) method, weighted-median method, MR-Egger methodand other methods were used to estimate causality. The horizontal pleiotropy, heterogeneities were also verified through the MR-Egger intercept test, Cochran's Q test, MR-PRESSO test and the leave-one-out sensitivity analysis. Reverse MR analysis was employed to evaluate the direct impact of metabolites on lower back pain or/and sciatica. Additionally, we conducted the colocalization analysis to reflect the causality deeply. Furthermore, metabolic pathway analysis was performed. 28 metabolites (18 known metabolites, 1 identified metabolites and 9 unknown metabolites) relevant to the risk of sciatica or/and lower back pain after using genetic variants as probes at PIVW < 0.05 were identifed. Among them, 8 serum metabolites decreased risk of sciatica or/and lower back pain significantly (P < 0.05), and 14 serum metabolites increased risk of sciatica or/and lower back pain significantly (P < 0.05). No reverse causal association was found between 28 metabolites and sciatica or/and lower back pain. Colocalization analysis results showed that the associations between sciatica or/and lower back pain and the 28 identified metabolites were not due to shared causal variant sites. Moreover, pathway enrichment analysis identifed 11 signifcant metabolic pathways, which are mainly involved in the pathological mechanism of sciatica or/and lower back pain (P < 0.05). There was no horizontal pleiotropy or heterogeneity in the other analyses. Our analyses provided robust evidence of causal associations between blood metabolites on sciatica or/and lower back pain. However, the underlying mechanisms remain to be further investigated.

Whole-genome sequencing reveals rare variants associated with gout in Taiwanese males.

To identify rare variants (RVs) of gout, we sequenced the whole genomes of 321 male gout patients and combined these with those of 64 male gout patients and 682 normal controls at Taiwan Biobank. We performed ACAT-O to identify 682 significant RVs (p < 3.8 × 10-8) clustered on chromosomes 1, 7, 10, 16, and 18. To prioritize causal variants effectively, we sifted them by Combined Annotation-Dependent Depletion score >10 or |effect size| ≥ 1.5 for those without CADD scores. In particular, to the best of our knowledge, we identified the rare variants rs559954634, rs186763678, and 13-85340782-G-A for the first time to be associated with gout in Taiwanese males. Importantly, the RV rs559954634 positively affects gout, and its neighboring gene NPHS2 is involved in serum urate and expressed in kidney tissues. The kidneys play a major role in regulating uric acid levels. This suggests that rs559954634 may be involved in gout. Furthermore, rs186763678 is in the intron of NFIA that interacts with SLC2A9, which has the most significant effect on serum urate. Note that gene-gene interaction NFIA-SLC2A9 is significantly associated with serum urate in the Italian MICROS population and a Croatian population. Moreover, 13-85340782-G-A significantly affects gout susceptibility (odds ratio 6.0; P = 0.038). The >1% carrier frequencies of these potentially pathogenic (protective) RVs in cases (controls) suggest the revealed associations may be true; these RVs deserve further studies for the mechanism. Finally, multivariate logistic regression analysis shows that the rare variants rs559954634 and 13-85340782-G-A jointly are significantly associated with gout susceptibility.

PSAP-Genomic-Regions: A Method Leveraging Population Data to Prioritize Coding and Non-Coding Variants in Whole Genome Sequencing for Rare Disease Diagnosis.

The introduction of Next-Generation Sequencing technologies in the clinics has improved rare disease diagnosis. Nonetheless, for very heterogeneous or very rare diseases, more than half of cases still lack molecular diagnosis. Novel strategies are needed to prioritize variants within a single individual. The Population Sampling Probability (PSAP) method was developed to meet this aim but only for coding variants in exome data. Here, we propose an extension of the PSAP method to the non-coding genome called PSAP-genomic-regions. In this extension, instead of considering genes as testing units (PSAP-genes strategy), we use genomic regions defined over the whole genome that pinpoint potential functional constraints. We conceived an evaluation protocol for our method using artificially generated disease exomes and genomes, by inserting coding and non-coding pathogenic ClinVar variants in large data sets of exomes and genomes from the general population. PSAP-genomic-regions significantly improves the ranking of these variants compared to using a pathogenicity score alone. Using PSAP-genomic-regions, more than 50% of non-coding ClinVar variants were among the top 10 variants of the genome. On real sequencing data from six patients with Cerebral Small Vessel Disease and nine patients with male infertility, all causal variants were ranked in the top 100 variants with PSAP-genomic-regions. By revisiting the testing units used in the PSAP method to include non-coding variants, we have developed PSAP-genomic-regions, an efficient whole-genome prioritization tool which offers promising results for the diagnosis of unresolved rare diseases.

Genetic Insights into Externalizing and Internalizing Traits through Integration of the Research Domain Criteria and Hierarchical Taxonomy of Psychopathology.

There is considerable comorbidity between externalizing (EXT) and internalizing (INT) psychopathology. Understanding the shared genetic underpinnings of these spectra is crucial for advancing knowledge of their biological bases and potential health impacts, and for informing empirical models like the Research Domain Criteria (RDoC) and Hierarchical Taxonomy of Psychopathology (HiTOP). We conducted a multivariate genome-wide association study (GWAS) of EXT and INT psychopathology by applying genomic structural equation modeling to summary statistics from 16 EXT and INT traits in European-ancestry individuals (n = 16,400 to 1,074,629). Downstream analyses explored associations across RDoC units of analysis (i.e., genes, molecules, cells, circuits, physiology, and behaviors). The GWAS identified 409 lead single nucleotide polymorphisms (SNPs) for EXT, 85 for INT, and 256 for EXT+INT (i.e., shared) traits. Bivariate causal mixture models estimated that nearly all EXT and INT causal variants overlapped, despite a genetic correlation of 0.37 (SE = 0.02). Drug repurposing analyses identified potential therapeutic targets, including perturbagens affecting dopamine and serotonin pathways. EXT genes had enriched expression in GABAergic, cortical, and hippocampal neurons, while INT genes were more narrowly linked to GABAergic neurons. EXT+INT liability was associated with reduced grey matter volumes in the amygdala and subcallosal cortex. These findings reveal both genetic overlap and distinct molecular and neurobiological pathways underlying EXT and INT psychopathology. By integrating genomic insights with the RDoC and HiTOP frameworks, this study advances our understanding of the mechanisms driving these dimensions of psychopathology.

Low vitamin C status and hypermobility‐related disorders in patients with bleeding disorder of unknown cause

AbstractIntroductionBleeding disorder of unknown cause (BDUC) is a challenging diagnosis that predominantly affects women. Previous investigations into connective tissue disorders (CTD) and vitamin C have not been conducted.AimTo examine the association between hypermobility‐related disorders, vitamin C status and BDUC.MethodsPatients were selected following laboratory and genetic screening that yielded negative results for known hemostasis disorders. Sixty patients with BDUC and an ISTH BAT score ≥ 10 underwent clinically examination for skin hyperextensibility and for hypermobility assessed by Beighton score. Vitamin C was analyzed by high‐performance liquid chromatography. Genetic screening for causal variants in 42 CTD genes was performed.ResultsThe majority of patients were female (56/60). Median ISTH BAT score was 13 (range 10–23). Beighton score was positive in 29/60 patients compared to 1/20 healthy controls (HC) (p &lt; .001). Hyperextensive skin was observed in (18/60) patients, and none (0/20) of the HC (p = .0041). Ten patients met the clinical diagnostic criteria for hypermobile Ehlers–Danlos syndrome (hEDS), and one patient was diagnosed with Noonan syndrome. Genetic screening excluded various subtypes of EDS with known genetic backgrounds. Average vitamin C level was adequate, but lower than in HC (55.9 vs. 70.4 μmol/L; p = .001). Suboptimal, or low vitamin C were identified in 19/60 compared to 1/20 HC (p = .018).ConclusionOur study demonstrates that BDUC is frequently associated with hypermobility disorders and low vitamin C status. Our results could pave the way for a randomized study of vitamin C supplementation in patients with BDUC.

Genome-wide association study of hypersensitivity skin reactions induced by nonionic iodinated contrast media.

In Taiwan, nonionic iodinated contrast media (ICMs) are commonly used but can occasionally cause severe side effects. The infrequency of these adverse events, coupled with the complexities in establishing direct causality, poses significant challenges for genetic research. : To investigate the genetic factors associated with skin reactions mediated by nonionic ICMs on a genome-wide scale. A hospital-based cohort from the China Medical University Hospital biobank was utilized to conduct a comprehensive genome-wide association study (GWAS) using PLINK v1.9. The study incorporated two distinct cohorts: one based on adverse drug reaction (ADR) reports, capturing immediate reactions, and the other based on self-reports, which primarily reflected delayed reactions. Known loci were determined by the GWAS catalog. Fine mapping was conducted by FINEMAP to predict causal variants. Pathway enrichment analysis was performed by clusterProfiler to reveal the biological function of the identified genetic signatures. The ADR-based cohort included 120 cases and 3640 controls. GWAS identified 6 candidate risk loci, namely rs150515068, rs6847491, rs192044153, rs191908641, rs376660317, and rs368821335. The self-report-based cohort, consisting of 275 cases and 8338 controls, revealed 36 additional candidate risk loci. Fine mapping further identified 4 causal variants within each cohort. Pathway analysis showed that immediate HSR-related genes are linked to growth hormone response and signaling, while non-immediate HSR genes are involved in neurotransmission. This study offers new perspectives on the genetic foundation of nonionic ICM-induced skin reactions within the Taiwanese population, suggesting that the genes contributing to immediate and non-immediate HSRs might have different functional roles.

A rare transcript homozygous variants in CLRN1(USH3A) causes Usher syndrome type 3 in a Chinese family

BackgroundUsher syndrome type 3 (USH3) is an autosomal recessive inherited disorder caused by pathogenic variants in the CLRN1 gene.ObjectTo evaluate the genotype-phenotype correlation of Usher syndrome type 3 (USH3) in a deaf-blind Chinese family of 3 generations with 2 patients.MethodsWe collected blood samples and clinical data from all of the pedigree family members. Genomic DNA was isolated from peripheral leukocytes using standard method. Targeted next generation sequencing and Sanger sequencing were performed to find the pathogenic variants in this family. Digital PCR and plasmid overexpression assay were used to verify the pathogenicity of variant sites in different transcripts.ResultsAll patients developed bilateral sensorineural hearing loss (SHL), progressive vision loss and nyctalopia. NGS of genes for Usher syndrome, deafness and retinal dystrophy identified a locus mutation in CLRN1 that caused completely different amino acid changes in different transcripts[CLRN1:c.474T > A(P.Cys158Ter) at NM_001256819.2 or c.302T > A(p.Val101Asp) at NM_174878.3], and plasmid overexpression experiments confirmed that the c.474T > A(P.Cys158Ter, NM_001256819.2) was a pathogenic variant which has never been associated with Usher syndrome in China, and the transcript of this mutation was not the version commonly found worldwide.ConclusionsThe CLRN1c.474T > A(NM_001256819.2) mutation is the causative variant in the Chinese family with USH3. The pathogenicity of different transcripts should be particularly considered in pathogenicity analysis.

Exploring Genetic Drug Targets in Acne Vulgaris: A Comprehensive Proteome-Wide Mendelian Randomization Study.

Acne vulgaris presents a substantial clinical challenge due to its complex pathophysiology and significant impact on quality of life. Identification of novel therapeutic targets for acne using genetic tools can guide the development of more effective treatments. Utilizing a dataset comprising 35 559 Icelandic individuals, we performed proteomic analyses to quantify 4709 circulating proteins. We integrated these data with acne-specific genome-wide association studies (GWAS) encompassing 34 422 acne patients and 364 991 controls. Mendelian randomization (MR) analyses employed the TwoSampleMR tool and Summary-data-based Mendelian Randomization (SMR) to estimate the causal effects of identified proteins on acne risk. Colocalization analyses assessed the likelihood of shared genetic etiology between protein levels and acne using the "coloc" R package. Our proteome-wide MR analysis initially identified 128 proteins potentially associated with acne risk. Following multiple testing corrections using the Benjamini-Hochberg method, fatty acid synthase (FASN) and tissue inhibitor of metalloproteinases 4 (TIMP4) remained significantly associated with acne risk. FASN exhibited a protective effect against acne (OR = 0.768, 95% CI: 0.676-0.872, p = 4.685E-05), while TIMP4 was associated with an increased risk (OR = 1.169, 95% CI: 1.103-1.241, p = 1.956E-07). Colocalization analysis supported a shared genetic basis for these protein-acne associations, with posterior probabilities indicating strong evidence of shared causal variants. Our findings highlight the utility of integrative genomic approaches in identifying potential therapeutic targets for acne. FASN and TIMP4, in particular, demonstrate strong potential as targets for therapeutic intervention, pending further validation through clinical research. These results offer a foundation for targeted acne treatment development, aligning with personalized medicine principles.

SLC39A8.p.(Ala391Thr) is associated with poorer cognitive ability: a cross-sectional study of schizophrenia and the general UK population.

The missense SNP NC_000004.12:g.102267552C>T (SLC39A8.p.(Ala391Thr), rs13107325) in SLC39A8, which encodes a zinc transporter, has been linked to schizophrenia and is the likely causal variant for one of the genome-wide association loci associated with the disorder. We tested whether the schizophrenia-risk allele at p.(Ala391Thr) was associated with schizophrenia-related phenotypes, including positive, negative, and disorganised symptoms, cognitive ability, educational attainment, and age of psychosis onset, within three schizophrenia cohorts (combined N=1,232) and, with equivalent phenotypes, in a sample of population controls (UK Biobank, N=355,069). We used regression analyses controlling for age, sex, and population stratification. Within the schizophrenia cohorts, after correction for multiple testing, p.(Ala391Thr) was not significantly associated with any schizophrenia-related phenotypes. In the unaffected participants from the UK Biobank, the schizophrenia-risk allele at p.(Ala391Thr) was associated with significantly poorer cognitive ability and fluid intelligence, a lower probability of obtaining GCSEs or a degree-level qualification, and fewer years in education. There was no association between p.(Ala391Thr) and self-reported psychotic experiences in this cohort. The schizophrenia-risk allele was associated with poorer cognitive ability, but not psychotic experiences, in a volunteer sample drawn from of the general population. To determine whether p.(Ala391Thr) is associated with cognitive phenotypes in people with schizophrenia, and to understand the role of p.(Ala391Thr) in the aetiology of cognitive impairment in schizophrenia, larger independent samples are required.

Age-dependent genetic architectures of chicken body weight explored by multidimensional GWAS and molQTL analyses

Chicken body weight (BW) is a critical trait in breeding. Although genetic variants associated with BW have been investigated by genome-wide association studies (GWAS), the contributions of causal variants and their molecular mechanisms remain largely unclear in chickens. In this study, we construct a comprehensive genetic atlas of chicken BW by an integrative analysis of 30 age points and 5 quantitative trait loci (QTL) across 27 tissues. We find that chicken growth is a cumulative non-linear process, which can be divided into three distinct stages. Our GWAS analysis reveals that BW-related genetic variations show ordered patterns in these three stages. Genetic variations in chromosome 1 may regulate the overall growth process, likely by modulating the hypothalamus-specific expression of SLC25A30 and retina-specific expression of NEK3. Moreover, genetic variations in chromosome 4 and chromosome 27 may play dominant roles in regulating BW during stage Ⅱ (8–22 weeks) and stage Ⅲ (23–72 weeks), respectively. In summary, our study presents a comprehensive genetic atlas regulating developmental stage-specific changes in chicken BW, thus providing important resources for genomic selection in breeding programs.

Human hypofunctional NCF1 variants promote pulmonary fibrosis in the bleomycin-induced mouse model and patients with systemic sclerosis via expansion of SPP1+ monocytes-derived macrophages

ObjectiveWe assessed the role of a systemic lupus erythematosus causal hypofunctional variant, neutrophil cytosolic factor 1 (NCF1)-p.Arg90His (p.R90H) substitution, in systemic sclerosis (SSc).MethodsAssociation of NCF1-H90 with SSc was performed in...