Genomic variants in exons and introns: identifying the splicing spoilers

SpliceTransformer predicts tissue-specific splicing linked to human diseases

We present SpliceTransformer (SpTransformer), a deep-learning framework that predicts tissue-specific RNA splicing alterations linked to human diseases based on genomic sequence. SpTransformer outperforms all previous methods on splicing prediction. Application to approximately 1.3 million genetic variants in the ClinVar database reveals that splicing alterations account for 60% of intronic and synonymous pathogenic mutations, and occur at different frequencies across tissue types. Importantly, tissue-specific splicing alterations match their clinical manifestations independent of gene expression variation. We validate the enrichment in three brain disease datasets involving over 164,000 individuals. Additionally, we identify single nucleotide variations that cause brain-specific splicing alterations, and find disease-associated genes harboring these single nucleotide variations with distinct expression patterns involved in diverse biological processes. Finally, SpTransformer analysis of whole exon sequencing data from blood samples of patients with diabetic nephropathy predicts kidney-specific RNA splicing alterations with 83% accuracy, demonstrating the potential to infer disease-causing tissue-specific splicing events. SpTransformer provides a powerful tool to guide biological and clinical interpretations of human diseases.

Association of PPARGC1A gene polymorphism and mtDNA methylation with coal-burning fluorosis: a case-control study

BACKGROUND

Coal-burning fluorosis is a chronic poisoning resulting from the prolonged use of locally available high-fluoride coal for heating and cooking. Prolonged fluoride exposure has been demonstrated to decrease PPARGC1A levels. Therefore, this case-control aims to evaluate the genetic association of PPARGC1A gene polymorphisms and methylation of the mitochondrial D-loop region with coal-burning fluorosis.

RESULT

The results showed that the TT genotype at rs13131226 and the AA genotype at rs1873532 increased the risk of coal-burning fluorosis (OR = 1.84, P = 0.004; OR = 1.97, P = 0.007), the CT and CC genotypes at rs7665116 decreased the risk of coal-burning fluorosis (OR = 0.54, P = 0.003). The TT genotype at the rs2970847 site and the AA genotype at the rs2970870 site increase the risk of developing skeletal fluorosis (OR = 4.12, P = 0.003; OR = 2.22, P = 0.011). Haplotype AG constructed by rs3736265-rs1873532 increased the risk of the prevalence of coal-burning fluorosis (OR = 1.465, P = 0.005); CG decreased the risk of the prevalence of coal-burning fluorosis (OR = 0.726, P = 0.020). Haplotype CGGT constructed by rs6821591-rs768695-rs3736265-rs2970847 increased the risk of the prevalence of skeletal fluorosis (OR = 1.558, P = 0.027). A 1% increase in CpG_4 methylation levels in the mtDNA D-loop region is associated with a 2.3% increase in the risk of coal-burning fluorosis. Additionally. There was a significant interaction between rs13131226 and rs1873532; CpG_4 and CpG_8.9; rs13131224,rs6821591 and rs7665116 were observed in the occurrence of fluorosis in the Guizhou population (χ2 = 16.917, P < 0.001; χ2 = 21.198, P < 0.001; χ2 = 36.078, P < 0.001).

CONCLUSION

PPARGC1A polymorphisms rs13131226 and rs1873532 and the mitochondrial DNA D-loop methylation site CpG_4 have been associated with an increased risk of fluorosis, conversely polymorphism rs7665116 was associated with a decreased risk of fluorosis. Polymorphisms rs2970870 were associated with increased risk of skeletal fluorosis, and polymorphism rs2970847 was associated with decreased risk of skeletal fluorosis. These SNPs and CpG can be used as potential targets to assess fluorosis risk.

Oligonucleotide therapies for nonalcoholic steatohepatitis

Nonalcoholic steatohepatitis (NASH) represents a severe disease subtype of nonalcoholic fatty liver disease (NAFLD) that is thought to be highly associated with systemic metabolic abnormalities. It is characterized by a series of substantial liver damage, including hepatocellular steatosis, inflammation, and fibrosis. The end stage of NASH, in some cases, may result in cirrhosis and hepatocellular carcinoma (HCC). Nowadays a large number of investigations are actively under way to test various therapeutic strategies, including emerging oligonucleotide drugs (e.g., antisense oligonucleotide, small interfering RNA, microRNA, mimic/inhibitor RNA, and small activating RNA) that have shown high potential in treating this fatal liver disease. This article systematically reviews the pathogenesis of NASH/NAFLD, the promising druggable targets proven by current studies in chemical compounds or biological drug development, and the feasibility and limitations of oligonucleotide-based therapeutic approaches under clinical or pre-clinical studies.

Klumpy: A tool to evaluate the integrity of long-read genome assemblies and illusive sequence motifs

The improvement and decreasing costs of third-generation sequencing technologies has widened the scope of biological questions researchers can address with de novo genome assemblies. With the increasing number of reference genomes, validating their integrity with minimal overhead is vital for establishing confident results in their applications. Here, we present Klumpy, a tool for detecting and visualizing both misassembled regions in a genome assembly and genetic elements (e.g. genes) of interest in a set of sequences. By leveraging the initial raw reads in combination with their respective genome assembly, we illustrate Klumpy's utility by investigating antifreeze glycoprotein (afgp) loci across two icefishes, by searching for a reported absent gene in the northern snakehead fish, and by scanning the reference genomes of a mudskipper and bumblebee for misassembled regions. In the two former cases, we were able to provide support for the noncanonical placement of an afgp locus in the icefishes and locate the missing snakehead gene. Furthermore, our genome scans were able identify an unmappable locus in the mudskipper reference genome and identify a putative repetitive element shared among several species of bees.

MC1R and age heteroclassification of face phenotypes in the Rio Grande do Sul population

BACKGROUND

Forensic DNA phenotyping (FDP) consists of the use of methodologies for predicting externally visible characteristics (EVCs) from the genetic material of biological samples found in crime scenes and has proven to be a promising tool in aiding human identification in police activities. Currently, methods based on multiplex assays and statistical models of prediction of EVCs related to hair, skin, and iris pigmentation using panels of SNP and INDEL biomarkers have already been developed and validated by the forensic scientific community. As well as traces of pigmentation, an individual's perceived age (PA) can also be considered an EVC and its estimation in unknown individuals can be useful for the progress of investigations. Liu and colleagues (2016) were pioneers in evidencing that, in addition to lifestyle and environmental factors, the presence of SNP and INDEL variants in the MC1R gene - which encodes a transmembrane receptor responsible for regulating melanin production - seems to contribute to an individual's PA. The group highlighted the association between these MC1R gene polymorphisms and the PA in the European population, where carriers of risk haplotypes appeared to be up to 2 years older in comparison to their chronological age (CA).

PURPOSE

Understanding that genotype-phenotype relationships cannot be extrapolated between different population groups, this study aimed to test this hypothesis and verify the applicability of this variant panel in the Rio Grande do Sul admixed population.

METHODS

Based on genomic data from a sample of 261 volunteers representative of gaucho population and using a multiple linear regression (MLR) model, our group was able to verify a significant association among nine intronic variants in loci adjacent to MC1R (e.g., AFG3L1P, TUBB3, FANCA) and facial age appearance, whose PA was defined after age heteroclassification of standard frontal face images through 11 assessors.

RESULTS

Different from that observed in European populations, our results show that the presence of effect alleles (R) of the selected variants in our sample influenced both younger and older face phenotypes. The influence of each variant on PA is expressed as β values.

CONCLUSIONS

There are important molecular mechanisms behind the effects of MC1R locus on PA, and the genomic background of each population seems to be crucial to determine this influence.

Splice-Switching Antisense Oligonucleotides Correct Phenylalanine Hydroxylase Exon 11 Skipping Defects and Rescue Enzyme Activity in Phenylketonuria

The PAH gene encodes the hepatic enzyme phenylalanine hydroxylase (PAH), and its deficiency, known as phenylketonuria (PKU), leads to neurotoxic high levels of phenylalanine. PAH exon 11 is weakly defined, and several missense and intronic variants identified in patients affect the splicing process. Recently, we identified a novel intron 11 splicing regulatory element where U1snRNP binds, participating in exon 11 definition. In this work, we describe the implementation of an antisense strategy targeting intron 11 sequences to correct the effect of PAH mis-splicing variants. We used an in vitro assay with minigenes and identified splice-switching antisense oligonucleotides (SSOs) that correct the exon skipping defect of PAH variants c.1199+17G>A, c.1199+20G>C, c.1144T>C, and c.1066-3C>T. To examine the functional rescue induced by the SSOs, we generated a hepatoma cell model with variant c.1199+17G>A using CRISPR/Cas9. The edited cell line reproduces the exon 11 skipping pattern observed from minigenes, leading to reduced PAH protein levels and activity. SSO transfection results in an increase in exon 11 inclusion and corrects PAH deficiency. Our results provide proof of concept of the potential therapeutic use of a single SSO for different exonic and intronic splicing variants causing PAH exon 11 skipping in PKU.

Pleiotropic effects of different exonic nucleotide changes at the same position contribute to hemophilia B phenotypic variation

BACKGROUND

The disease-causing effects of genetic variations often depend on their location within a gene. Exonic changes generally lead to alterations in protein production, secretion, activity, or clearance. However, owing to the overlap between proteins and splicing codes, missense variants can also affect messenger RNA splicing, thus adding a layer of complexity and influencing disease phenotypes.

OBJECTIVES

To extensively characterize a panel of 13 exonic variants in the F9 gene occurring at 6 different factor IX positions and associated with varying severities of hemophilia B (HB).

METHODS

Computational predictions, splicing analysis, and recombinant factor IX assays were exploited to characterize F9 variants.

RESULTS

We demonstrated that 5 (38%) of 13 selected F9 exonic variants have pleiotropic effects. Although bioinformatic approaches accurately classified effects, extensive experimental assays were required to elucidate and deepen the molecular mechanisms underlying the pleiotropic effects. Importantly, their characterization was instrumental in developing tailored RNA therapeutics based on engineered U7 small nuclear RNA to mask cryptic splice sites and compensatory U1 small nuclear RNA to enhance exon definition.

CONCLUSION

Overall, albeit a multitool bioinformatic approach suggested the molecular effects of multiple HB variants, the deep investigation of molecular mechanisms revealed insights into the HB phenotype-genotype relationship, enabling accurate classification of HB variants. Importantly, knowledge of molecular mechanisms allowed the development of tailored RNA therapeutics, which can also be translated to other genetic diseases.

MAOB expression correlates with a favourable prognosis in prostate cancer, and its genetic variants are associated with the metastasis of the disease

Monoamine oxidase B (MAOB), a neurotransmitter-degrading enzyme, was reported to reveal conflicting roles in various cancers. However, the functional role of MAOB and impacts of its genetic variants on prostate cancer (PCa) is unknown. Herein, we genotyped four loci of MAOB single-nucleotide polymorphisms (SNPs), including rs1799836 (A/G), rs3027452 (G/A), rs6651806 (A/C) and rs6324 (G/A) in 702 PCa Taiwanese patients. We discovered that PCa patients carrying the MAOB rs6324 A-allele exhibited an increased risk of having a high initial prostate-specific antigen (iPSA) level (>10 ng/mL). Additionally, patients with the rs3027452 A-allele had a higher risk of developing distal metastasis, particularly in the subpopulation with high iPSA levels. In a subpopulation without postoperative biochemical recurrence, patients carrying the rs1799836 G-allele had a higher risk of developing lymph node metastasis and recurrence compared to those carrying the A-allele. Furthermore, genotype screening in PCa cell lines revealed that cells carrying the rs1799836 G-allele expressed lower MAOB levels than those carrying the A-allele. Functionally, overexpression and knockdown of MAOB in PCa cells respectively suppressed and enhanced cell motility and proliferation. In clinical observations, correlations of lower MAOB expression levels with higher Gleason scores, advanced clinical T stages, tumour metastasis, and poorer prognosis in PCa patients were noted. Our findings suggest that MAOB may act as a suppressor of PCa progression, and the rs3027452 and rs1799836 genetic variants of MAOB are linked to PCa metastasis within the Taiwanese population.

Expression of the GCG gene and secretion of active glucagon-like peptide-1 varies along the length of intestinal tract in horses

BACKGROUND

Active glucagon-like peptide-1 (aGLP-1) has been implicated in the pathogenesis of equine insulin dysregulation (ID), but its role is unclear. Cleavage of proglucagon (coded by the GCG gene) produces aGLP-1 in enteral L cells.

OBJECTIVES

The aim in vivo was to examine the sequence of the exons of GCG in horses with and without ID, where aGLP-1 was higher in the group with ID. The aims in vitro were to identify and quantify the expression of GCG in the equine intestine (as a marker of L cells) and determine intestinal secretion of aGLP-1.

STUDY DESIGN

Genomic studies were case-control studies. Expression and secretion studies in vitro were cross-sectional.

METHODS

The GCG gene sequence of the exons was determined using a hybridisation capture protocol. Expression and quantification of GCG in samples of stomach duodenum, jejunum, ileum, caecum and ascending and descending colon was achieved with droplet digital PCR. For secretory studies tissue explants were incubated with 12 mM glucose and aGLP-1 secretion was measured with an ELISA.

RESULTS

Although the median [IQR] post-prandial aGLP-1 concentrations were higher (p = 0.03) in animals with ID (10.2 [8.79-15.5]), compared with healthy animals (8.47 [6.12-11.7]), there was 100% pairwise identity of the exons of the GCG sequence for the cohort. The mRNA concentrations of GCG and secretion of aGLP-1 differed (p < 0.001) throughout the intestine.

MAIN LIMITATIONS

Only the exons of the GCG gene were sequenced and breeds were not compared. The horses used for the study in vitro were not assessed for ID and different horses were used for the small, and large, intestinal studies.

CONCLUSIONS

Differences in post-prandial aGLP-1 concentration were not due to a variant in the exons of the GCG gene sequence in this cohort. Both the large and small intestine are sites of GLP-1 secretion.

Alternative Splicing Landscape of Head and Neck Squamous Cell Carcinoma

Head and neck malignancies are a significant global health concern, with head and neck squamous cell carcinoma (HNSCC) being the sixth most common cancer worldwide accounting for > 90% of cases. In recent years, there has been growing recognition of the potential role of alternative splicing (AS) in the etiology of cancer. Increasing evidence suggests that AS is associated with various aspects of cancer progression, including tumor occurrence, invasion, metastasis, and drug resistance. Additionally, AS is involved in shaping the tumor microenvironment, which plays a crucial role in tumor development and response to therapy. AS can influence the expression of factors involved in angiogenesis, immune response, and extracellular matrix remodeling, all of which contribute to the formation of a supportive microenvironment for tumor growth. Exploring the mechanism of AS events in HNSCC could provide insights into the development and progression of this cancer, as well as its interaction with the tumor microenvironment. Understanding how AS contributes to the molecular changes in HNSCC cells and influences the tumor microenvironment could lead to the identification of new therapeutic targets. Targeted chemotherapy and immunotherapy strategies tailored to the specific AS patterns in HNSCC could potentially improve treatment outcomes and reduce side effects. This review explores the concept, types, processes, and technological advancements of AS, focusing on its role in the initiation, progression, treatment, and prognosis of HNSCC.

The role of cell adhesion molecule IgSF9b at the inhibitory synapse and psychiatric disease

Understanding perturbations in synaptic function between health and disease states is crucial to the treatment of neuropsychiatric illness. While genome-wide association studies have identified several genetic loci implicated in synaptic dysfunction in disorders such as autism and schizophrenia, many have not been rigorously characterized. Here, we highlight immunoglobulin superfamily member 9b (IgSF9b), a cell adhesion molecule thought to localize exclusively to inhibitory synapses in the brain. While both pre-clinical and clinical studies suggest its association with psychiatric diseases, our understanding of IgSF9b in synaptic maintenance, neural circuits, and behavioral phenotypes remains rudimentary. Moreover, these functions wield undiscovered influences on neurodevelopment. This review evaluates current literature and publicly available gene expression databases to explore the implications of IgSF9b dysfunction in rodents and humans. Through a focused analysis of one high-risk gene locus, we identify areas requiring further investigation and unearth clues related to broader mechanisms contributing to the synaptic etiology of psychiatric disorders.

The effect of insulin-like growth factor 1 receptor gene single nucleotide polymorphism on growth and milk production traits in two native Indian tropical goat breeds

Insulin like growth factor1(IGF-1) is an essential growth factor that mediates the growth-promoting functions of pituitary growth hormone. Insulin like growth factor 1 receptor (IGF1R) is a tyrosine kinase receptor that mediates the actions of IGF1. Therefore, IGF1R is a candidate gene for examining SNPs linked with growth and production traits. The objective of this study was to detect the c.546 + 179170A > T transversion in intron 2 of the gene encoding IGF1R in two goat breeds, Attappady Black and Malabari, and examine the association of this polymorphism with growth and milk production. For the identification of the SNP, the T-ARMS-PCR was utilized. All three genotypes were present in the two investigated breeds. The polymorphism was found to be significantly (p < 0.05) linked with growth traits. At birth, 3 and 6 months of age, Attappady goats with the AT genotype had significantly (p < 0.05) higher body weights than those with the AA and TT genotypes. Malabari goats with the AT genotype had significantly (p < 0.05) higher body weights at birth and at 3 months of age. The genotypes of the IGF1R gene had no effect on total or peak milk production. Therefore, this SNP could be used as a molecular marker in selection of meat-producing goat breeds.HIGHLIGHTSc.546 + 179170A > T IGF1R transversion was detected using T-ARMS-PCR in two indigenous goat breeds from Kerala.Attappady Black and Malabari goat breeds of Kerala possessed all the three genotypesIn these breeds, there was a significant correlation between this SNP, c.546 + 179170A > T IGF1R, and body weight.In both the breeds, there was no association between this SNP and milk production traits.

Construction of a risk scoring system using clinical factors and RYR2 polymorphisms for bleeding complications in patients on direct oral anticoagulants

Introduction: Bleeding is one of the most undesirable complications of direct oral anticoagulants (DOACs). While the ryanodine receptor (RYR2) has been related to cardiac diseases, research on bleeding complications is lacking. This study aimed to elucidate the association between RYR2 and bleeding risk to develop the risk scoring system in patients treated with DOACs. Methods: This study was a retrospective analysis of prospectively collected samples. We selected ten SNPs within the RYR2 gene, and two models were constructed (Model I: demographic factors only, Model II: demographic and genetic factors) in multivariable analysis. Independent risk factors for bleeding were used to develop a risk scoring system. Results: A total of 447 patients were included, and 49 experienced either major bleeding or clinically relevant non-major bleeding. In Model I, patients using rivaroxaban and experiencing anemia exhibited an increased bleeding risk after adjusting for covariates. Upon incorporating genetic factors into Model I, a significant association with bleeding was also observed in cases of overdosing on DOACs and in patients with a creatinine clearance (CrCl) < 30 mL/min, in addition to rivaroxaban and anemia (Model II). Among genetic factors, RYR2 rs12594 GG, rs17682073 AA, rs3766871 GG, and rs6678625 T alleles were associated with bleeding complications. The area under the receiver operating characteristic curve (AUROC) of Model I was 0.670, whereas that of Model II increased to 0.803, demonstrating better performance with the inclusion of genetic factors. Using the significant variables in Model II, a risk scoring system was constructed. The predicted bleeding risks for scores of 0, 1-2, 3-4, 5-6, 7-8, and 9-10 points were 0%, 1.2%, 4.6%, 15.7%, 41.7%, and 73.3%, respectively. Conclusion: This study revealed an association between RYR2 and bleeding complications among patients taking DOACs and established a risk scoring system to support individualized DOAC treatment for these patients.

Splicing defects in rare diseases: transcriptomics and machine learning strategies towards genetic diagnosis

Genomic variants affecting pre-messenger RNA splicing and its regulation are known to underlie many rare genetic diseases. However, common workflows for genetic diagnosis and clinical variant interpretation frequently overlook splice-altering variants. To better serve patient populations and advance biomedical knowledge, it has become increasingly important to develop and refine approaches for detecting and interpreting pathogenic splicing variants. In this review, we will summarize a few recent developments and challenges in using RNA sequencing technologies for rare disease investigation. Moreover, we will discuss how recent computational splicing prediction tools have emerged as complementary approaches for revealing disease-causing variants underlying splicing defects. We speculate that continuous improvements to sequencing technologies and predictive modeling will not only expand our understanding of splicing regulation but also bring us closer to filling the diagnostic gap for rare disease patients.

Novel COL4A3 synonymous mutation causes Alport syndrome coexistent with immunoglobulin A nephropathy in a woman: A case report

BACKGROUND

Alport syndrome (AS) is an inherited disease of the glomerular basement membrane caused by mutations in genes encoding α3, α4, or α5 chains of type IV collagen. It manifests with hematuria or proteinuria, which is often accompanied by hearing impairments and ocular abnormalities. Histopathologically, AS shows mesangial proliferation and sometimes incidental immunoglobulin A (IgA) deposition. Hematuria or proteinuria is also a common presentation in patients with IgA nephropathy that makes it difficult to differentially diagnose AS and IgA nephropathy solely based on these clinical and pathological features.

CASE SUMMARY

Herein, we present the case of a 59-year-old female patient who was admitted to our hospital with persistent microscopic hematuria and occasional proteinuria that had lasted for > 2 years. This patient had a familial history of renal disease and was diagnosed with autosomal dominant AS (ADAS) and IgA nephropathy based on the findings of renal biopsy as well as genetic testing performed using whole-exome sequencing, which suggested that the patient carried a novel heterozygous variation (c.888G>A:p.Gln296Gln) in the COL4A3 gene that enriches the mutation spectrum of ADAS. The proband received an angiotensin receptor blocker therapy after a definitive diagnosis was established. After one year of therapy, a significant reduction in proteinuria was observed. The number of microscopic red blood cells per high-power field decreased to one-quarter of the baseline levels. Renal function also maintained well during the follow-up.

CONCLUSION

Our case highlights the significance of performing kidney biopsy and genetic testing in the diagnosis of AS and familial IgA nephropathy.

Novel Gene Polymorphisms for Stable Warfarin Dose in a Korean Population: Genome-Wide Association Study

Warfarin has a narrow therapeutic window and high intra- and inter-individual variability. Considering that many published papers on genotype-guided dosing are derived from European populations, the aim of this study was to investigate novel genetic variants associated with the variability of stable warfarin dose in the Korean population with cardiac valve replacement, using the GWAS approach. This retrospective cohort study was performed from January 1982 to December 2020 at the Severance Cardiovascular Hospital of Yonsei University College of Medicine. GWAS was performed to identify associations between genotypes and the warfarin maintenance dose, by comparing the allele frequency of genetic variants between individuals. Then, the extent of genetic and non-genetic factors on the dose variability was determined by multivariable regression analysis. The study enrolled 214 participants, and the most robust signal cluster was detected on chromosome 16 around VKORC1. Followed by VKORC1, three novel variants (NKX2-6 rs310279, FRAS1 rs4386623, and FAM201A rs1890109) showed an association with stable warfarin dose requirement in univariate analysis. The algorithm was constructed by using multivariable analysis that includes genetic and non-genetic factors, and it could explain 58.5% of the variations in stable warfarin doses. In this variability, VKORC1 rs9934438 and FRAS1 rs4386623 accounted for 33.0% and 9.9%, respectively. This GWAS analysis identified the fact that three novel variants (NKX2-6 rs310279, FRAS1 rs4386623, and FAM201A rs1890109) were associated with stable warfarin doses. Additional research is necessary to validate the results and establish personalized treatment strategies for the Korean population.

Interactive Associations between PPARγ and PPARGC1A and Bisphosphonate-Related Osteonecrosis of the Jaw in Patients with Osteoporosis

Bisphosphonate-related osteonecrosis of the jaw (BRONJ) is a rare but severe adverse effect that can occur as a result of bisphosphonate treatment. This study aimed to examine the relationship between PPARγ and PPARGC1A polymorphisms and the BRONJ development in female osteoporosis patients undergoing bisphosphonate treatment. We prospectively conducted this nested case-control study at the Ewha Womans University Mokdong Hospital between 2014 and 2018. We assessed five single-nucleotide polymorphisms (SNPs) of PPARγ and six SNPs of PPARGC1A and performed a multivariable logistic regression analysis to determine the independent risk factors for developing BRONJ. There were a total of 123 patients included in this study and 56 patients (45.5%) developed BRONJ. In the univariate analysis, PPARGC1A rs2946385 and rs10020457 polymorphisms were significantly associated with BRONJ (p = 0.034, p = 0.020, respectively), although the results were not statistically significant in the multivariable analysis. Patients with the combined genotypes of GG in both PPARγ rs1151999 and PPARGC1A rs2946385 showed a 3.03-fold higher risk of BRONJ compared to individuals with other genotype combinations after adjusting for confounders (95% confidence interval (CI): 1.01-9.11). Old age (≥70 years) and duration of bisphosphonate use (≥60 months) increased the risk of BRONJ. The area under the receiver operating characteristic curve for the predicted probability was 0.78 (95% CI: 0.69-0.87, p < 0.001), demonstrating a satisfactory level of discriminatory power. Our study elucidated that PPARγ and PPARGC1A polymorphisms were interactively associated with BRONJ development. These results have potential implications for tailoring personalized treatments for females undergoing bisphosphonate therapy for osteoporosis.

Effects of the T-type calcium channel CaV3.2 R1584P mutation on absence seizure susceptibility in GAERS and NEC congenic rats models

RATIONALE

Low-voltage-activated or T-type Ca²⁺ channels play a key role in the generation of seizures in absence epilepsy. We have described a homozygous, gain of function substitution mutation (R1584P) in the Ca_V3.2 T-type Ca²⁺ channel gene (Cacna1h) in the Genetic Absence Epilepsy Rats from Strasbourg (GAERS). The non-epileptic control (NEC) rats, derived from the same original Wistar strains as GAERS but selectively in-breed not to express seizures, are null for the R1584P mutation. To study the effects of this mutation in rats who otherwise have a GAERS or NEC genetic background, we bred congenic GAERS-Cacna1hNEC (GAERS null for R1584P mutation) and congenic NEC-Cacna1hGAERS (NEC homozygous for R1584P mutation) and evaluated the seizure and behavioral phenotype of these strains in comparison to the original GAERS and NEC strains.

METHODS

To evaluate seizure expression in the congenic strains, EEG electrodes were implanted in NEC, GAERS, GAERS^-Cacna1hNEC without the R1584P mutation, and NEC^{-Cacna1hGAERS} with the R1584P mutation rats. In the first study, continuous EEG recordings were acquired from week 4 (when seizures begin to develop in GAERS) to week 14 of age (when GAERS display hundreds of seizures per day). In the second study, the seizure and behavioral phenotype of GAERS and NEC^{-Cacna1hGAERS} strains were evaluated during young age (6 weeks of age) and adulthood (16 weeks of age) of GAERS, NEC, GAERS^-Cacna1hNEC and NEC^{-Cacna1hGAERS}. The Open field test (OFT) and sucrose preference test (SPT) were performed to evaluate anxiety-like and depressive-like behavior, respectively. This was followed by EEG recordings at 18 weeks of age to quantify the seizures, and spike-wave discharge (SWD) cycle frequency. At the end of the study, the whole thalamus was collected for T-type calcium channel mRNA expression analysis.

RESULTS

GAERS had a significantly shorter latency to first seizures and an increased number of seizures per day compared to GAERS^-Cacna1hNEC. On the other hand, the presence of the R1584P mutation in the NEC^{-Cacna1hGAERS} was not enough to generate spontaneous seizures in their seizure-resistant background. 6 and 16-week-old GAERS and GAERS^-Cacna1hNEC rats showed anxiety-like behavior in the OFT, in contrast to NEC and NEC^{-Cacna1hGAERS}. Results from the SPT showed that the GAERS developed depressive-like in the SPT compared to GAERS^-Cacna1hNEC, NEC, and NEC^{-Cacna1hGAERS}. Analysis of the EEG at 18 weeks of age showed that the GAERS had an increased number of seizures per day, increased total seizure duration and a higher cycle frequency of SWD relative to GAERS^-Cacna1hNEC. However, the average seizure duration was not significantly different between strains. Quantitative real-time PCR showed that the T-type Ca²⁺ channel isoform Ca_V3.2 channel expression was significantly increased in GAERS compared to NEC, GAERS^-Cacna1hNEC and NEC^{-Cacna1hGAERS}. The presence of the R1584P mutation increased the total ratio of Ca_V3.2 + 25/-25 splice variants in GAERS and NEC^{-Cacna1hGAERS} compared to NEC and GAERS^-Cacna1hNEC.

DISCUSSION

The data from this study demonstrate that the R1584P mutation in isolation on a seizure-resistant NEC genetic background was insufficient to generate absence seizures, and that a GAERS genetic background can cause seizures even without the mutation. However, the study provides evidence that the R1584P mutation acts as a modulator of seizures development and expression, and depressive-like behavior in the SPT, but not the anxiety phenotype of the GAERS model of absence epilepsy.

Intronic position +9 and -9 are potentially splicing sites boundary from intronic variants analysis of whole exome sequencing data

Whole exome sequencing (WES) can also detect some intronic variants, which may affect splicing and gene expression, but how to use these intronic variants, and the characteristics about them has not been reported. This study aims to reveal the characteristics of intronic variant in WES data, to further improve the clinical diagnostic value of WES. A total of 269 WES data was analyzed, 688,778 raw variants were called, among these 367,469 intronic variants were in intronic regions flanking exons which was upstream/downstream region of the exon (default is 200 bps). Contrary to expectation, the number of intronic variants with quality control (QC) passed was the lowest at the +2 and -2 positions but not at the +1 and -1 positions. The plausible explanation was that the former had the worst effect on trans-splicing, whereas the latter did not completely abolish splicing. And surprisingly, the number of intronic variants that passed QC was the highest at the +9 and -9 positions, indicating a potential splicing site boundary. The proportion of variants which could not pass QC filtering (false variants) in the intronic regions flanking exons generally accord with "S"-shaped curve. At +5 and -5 positions, the number of variants predicted damaging by software was most. This was also the position at which many pathogenic variants had been reported in recent years. Our study revealed the characteristics of intronic variant in WES data for the first time, we found the +9 and -9 positions might be a potentially splicing sites boundary and +5 and -5 positions were potentially important sites affecting splicing or gene expression, the +2 and -2 positions seem more important splicing site than +1 and -1 positions, and we found variants in intronic regions flanking exons over ± 50 bps may be unreliable. This result can help researchers find more useful variants and demonstrate that WES data is valuable for intronic variants analysis.

Effect of OR51E1 Single Nucleotide Polymorphisms on Glioma Susceptibility in the Chinese Han Population

BACKGROUND

Glioma is one of the common primary intracranial tumors, which is heterogeneous among individuals with a low cure rate. Our study aimed to investigate the association between single nucleotide polymorphisms (SNPs) of the OR51E1 gene and glioma susceptibility in the Chinese Han population.

METHODS

A total of six SNPs on OR51E1 in 1,026 subjects (526 cases and 500 controls) were genotyped by MassARRAY iPLEX GOLD assay. The association between these SNPs and glioma susceptibility was analyzed using logistic regression, and odds ratios (ORs) and 95% confidence intervals (CIs) were also calculated. The multifactor dimensionality reduction (MDR) method was applied to detect "SNP-SNP" interactions.

RESULTS

In the overall sample, polymorphisms rs10768148, rs7102992, and rs10500608 were identified to be associated with glioma risk. In the stratified analysis based on gender, only polymorphism rs10768148 was observed to be associated with the risk of glioma. In the age-stratified analysis, rs7102992, rs74052483, and rs10500609 contributed to the risk of glioma in subjects aged > 40 years. And polymorphisms rs10768148 and rs7102992 were associated with the risk of glioma in subjects aged ≤ 40 years and subjects with astrocytoma. In addition, a strong synergistic relationship between rs74052483 and rs10768148, and a strong redundant relationship between rs7102992 and rs10768148 were identified in the study.

CONCLUSIONS

This study demonstrated the association of OR51E1 polymorphisms with glioma susceptibility, providing a basis for assessing glioma risk-associated variants in the Chinese Han population.

Platelet Serotonin (5-HT) Concentration, Platelet Monoamine Oxidase B (MAO-B) Activity and HTR2A, HTR2C, and MAOB Gene Polymorphisms in Asthma

The complex role of the serotonin system in respiratory function and inflammatory diseases such as asthma is unclear. Our study investigated platelet serotonin (5-HT) levels and platelet monoamine oxidase B (MAO-B) activity, as well as associations with HTR2A (rs6314; rs6313), HTR2C (rs3813929; rs518147), and MAOB (rs1799836; rs6651806) gene polymorphisms in 120 healthy individuals and 120 asthma patients of different severity and phenotypes. Platelet 5-HT concentration was significantly lower, while platelet MAO-B activity was considerably higher in asthma patients; however, they did not differ between patients with different asthma severity or phenotypes. Only the healthy subjects, but not the asthma patients, carrying the MAOB rs1799836 TT genotype had significantly lower platelet MAO-B activity than the C allele carriers. No significant differences in the frequency of the genotypes, alleles, or haplotypes for any of the investigated HTR2A, HTR2C and MAOB gene polymorphisms have been observed between asthma patients and healthy subjects or between patients with various asthma phenotypes. However, the carriers of the HTR2C rs518147 CC genotype or C allele were significantly less frequent in severe asthma patients than in the G allele carriers. Further studies are necessary to elucidate the involvement of the serotonergic system in asthma pathophysiology.

ASAS-EGB: A statistical framework for estimating allele-specific alternative splicing events using transcriptome data

Ubiquitous Alternative splicing (AS) in eukaryotes greatly increases the biodiversity of proteins and involves in disease and cancer. Allele-specific AS studies can facilitate the identification of cis-acting elements because both alleles share the same cellular environment. Due to the limited information provided on the exons defined by AS events, we propose a statistical framework and algorithm ASAS-EGB for ASAS analysis using the gene transcriptome. The framework obtains exclusively compatible sets of gene isoforms supporting each event isoform, and uses both phased and non-phased SNPs within the exons on the gene isoforms for inference. Using this strategy, we have demonstrated ASAS-EGB can yield better ASAS inferential performance than using event isoforms. ASAS-EGB supports both single-end and paired-end RNA-seq data, and we have proved its robustness using RNA-seq replicates of individual NA12878. ASAS-EGB builds Bayesian models for ASAS analysis, and the MCMC method is used to solve the problem. With more detailed annotations for individual genomes and transcriptomes appearing in the future, the algorithm proposed by the paper can provide better support for these data to reveal the regulatory mechanisms of individual genomes.

Relation between single nucleotide polymorphisms in circadian clock relevant genes and cholesterol metabolism

Single nucleotide polymorphisms (SNPs) in circadian clock relevant genes are associated with several metabolic health variables, but little is known about their associations with human cholesterol metabolism. Therefore, this study examined associations between SNPs in ARNTL, ARNTL2, CLOCK, CRY1, CRY2, PER2, and PER3 with the intestinal cholesterol absorption markers campesterol and sitosterol, the endogenous cholesterol synthesis marker lathosterol, and total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C) concentrations in 456 healthy individuals from Western European descent. One SNP in ARNTL2 (rs1037924) showed a significant association with lathosterol. Several SNPs in ARNTL (rs4146388, rs58901760, rs6486121), ARNTL2 (rs73075788), CLOCK (rs13113518, rs35115774, rs6832769), and CRY1 (rs2078074) were significantly associated with intestinal cholesterol absorption. Genetic variants in CRY2, PER2, and PER3 were not significantly associated with intestinal cholesterol absorption or endogenous cholesterol synthesis. None of the SNPs were associated with TC or LDL-C, except for one SNP in PER2 (rs11894491) with serum LDL-C concentrations. The findings suggest that various SNPs in ARNTL, ARNTL2, CLOCK and CRY1 play a role in intestinal cholesterol absorption and endogenous cholesterol synthesis, which was not reflected in TC and LDL-C concentrations. The significant associations between SNPs and intestinal cholesterol absorption and endogenous cholesterol synthesis should be validated in other cohorts.

Optimized Detection of Unknown MET Exon 14 Skipping Mutations in Routine Testing for Patients With Non-Small-Cell Lung Cancer

PURPOSE

MET exon 14 (METex14) skipping is an actionable biomarker in non-small-cell lung cancer. However, MET variants are highly complex and diverse, and not all variants lead to exon 14 skipping. Assessing the skipping effect of unknown variants is still a key issue in molecular diagnosis.

MATERIALS AND METHODS

We retrospectively collected MET variants around exon 14 from 4,233 patients with non-small-cell lung cancer who underwent next-generation sequencing testing using DNA, as well as two published data sets.

RESULTS

Among the 4,233 patients, 44 unique variants including 29 novel variants (65.9%) were discovered from 53 patients. Notably, 31 samples (58.5%) failed RNA verification. Using RNA verification, nine novel skipping variants and five nonskipping variants were confirmed. We further used SpliceAI with the delta score cutoff of 0.315 to aid the classification of novel variants (sensitivity = 98.88% and specificity = 100%). When applied to the reported variants, we also found three wrongly classified nonskipping variants. Finally, an optimized knowledge-based interpretation procedure for clinical routine was built according to the mutation type and location, and five more skipping mutations from the 13 unknown variants were determined, which improved the population determination rate to 0.92%.

CONCLUSION

This study discovered more METex14 skipping variants and optimized an innovative approach that could be adapted for the interpretation of infrequent or novel METex14 variants timely without experimental validation.

Antisense oligonucleotide is a promising intervention for liver diseases

As the body's critical metabolic organ, the liver plays an essential role in maintaining proper body homeostasis. However, as people's living standards have improved and the number of unhealthy lifestyles has increased, the liver has become overburdened. These have made liver disease one of the leading causes of death worldwide. Under the influence of adverse factors, liver disease progresses from simple steatosis to hepatitis, to liver fibrosis, and finally to cirrhosis and cancer, followed by increased mortality. Until now, there has been a lack of accepted effective treatments for liver disease. Based on current research, antisense oligonucleotide (ASO), as an alternative intervention for liver diseases, is expected to be an effective treatment due to its high efficiency, low toxicity, low dosage, strong specificity, and additional positive characteristics. In this review, we will first introduce the design, modification, delivery, and the mechanisms of ASO, and then summarize the application of ASO in liver disease treatment, including in non-alcoholic fatty liver disease (NAFLD), hepatitis, liver fibrosis, and liver cancer. Finally, we discuss challenges and perspectives on the transfer of ASO drugs into clinical use. This review provides a current and comprehensive understanding of the integrative and systematic functions of ASO for its use in liver disease.

FOXM1 Variant Contributes to Gefitinib Resistance via Activating Wnt/β-Catenin Signal Pathway in Patients with Non-Small Cell Lung Cancer

PURPOSE

Although gefitinib prolonged the progression-free survival (PFS) of patients with non-small cell lung cancer (NSCLC), unpredictable resistance limited its clinical efficacy. Novel predictive biomarkers with explicit mechanisms are urgently needed.

EXPERIMENTAL DESIGN

A total of 282 patients with NSCLC with gefitinib treatment were randomly assigned in a 7:3 ratio to exploratory (n = 192) and validation (n = 90) cohorts. The candidate polymorphisms were selected with Haploview4.2 in Hapmap and genotyped by a MassARRAY system, and the feature variables were identified through Randomforest Survival analysis. Tanswell and clonogenic assays, base editing and cell-derived tumor xenograft model were performed to uncover the underlying mechanism.

RESULTS

We found that the germline missense polymorphism rs3742076 (A>G, S628P), located in transactivation domain of FOXM1, was associated with PFS in exploratory (median PFS: GG vs. GA&AA, 9.20 vs. 13.37 months, P = 0.00039, HR = 2.399) and validation (median PFS: GG vs. GA&AA, 8.13 vs. 13.80 months, P = 0.048, HR = 2.628) cohorts. We elucidated that rs3742076_G conferred resistance to gefitinib by increasing protein stability of FOXM1 and facilitating an aggressive phenotype in vitro and in vivo through activating wnt/β-catenin signaling pathway. Meanwhile, FOXM1 level was highly associated with prognosis in patients with EGFR-mutant NSCLC. Mechanistically, FOXM1 rs3742076_G upregulated wnt/β-catenin activity by directly binding to β-catenin in cytoplasm and promoting transcription of β-catenin in nucleus. Remarkably, inhibition of β-catenin markedly reversed rs3742076_G-induced gefitinib resistance and aggressive phenotypes.

CONCLUSIONS

These findings characterized rs3742076_G as a gain-of-function polymorphism in mediating gefitinib resistance and tumor aggressiveness, and highlighted the variant as a predictive biomarker in guiding gefitinib treatment.

A Study of Combined Genotype Effects of SHCBP1 on Wool Quality Traits in Chinese Merino

SHCBP1 (Shc SH2-domain binding protein 1) is a member of the Src and collagen homolog (Shc) protein family and is closely associated with multiple signaling pathways that play important roles during hair follicle induction, morphogenesis, and cycling. The purpose of this study was to investigate SHCBP1 gene expression, polymorphisms, and the association between SHCBP1 and wool quality traits in Chinese Merino sheep. The SHCBP1 gene was shown, by qPCR, to be ubiquitously expressed in sheep tissues and differentially expressed in the adult skin of Chinese Merino and Suffolk sheep. Four SNPs (termed SHCBP1SNPs 1-4) were identified by Sanger sequencing and were located in exon 2, intron 9, intron 12, and exon 13 of the sheep SHCBP1 gene, respectively. SHCBP1SNPs 3 and 4 (rs411176240 and rs160910635) were significantly associated with wool crimp (P < 0.05). The combined polymorphism (SHCBP1SNP3-SHCBP1SNP4) was significantly associated with wool crimp (P < 0.05). Bioinformatics analysis showed that the SNPs associated with wool crimp (SHCBP1SNPs 3 and 4) might affect the pre-mRNA splicing by creating binding sites for serine-arginine-rich proteins and that SHCBP1SNP4 might alter the SHCBP1 mRNA and protein secondary structure. Our results suggest that SHCBP1 influences wool crimp and SHCBP1SNPs 3 and 4 might be useful markers for marker-assisted selection and sheep breeding.

Association of the SLC47A1 Gene Variant With Responses to Metformin Monotherapy in Drug-naive Patients With Type 2 Diabetes

CONTEXT

Although metformin is the first-line treatment for type 2 diabetes, the blood sugar-lowering effect of metformin varies among populations. SLC47A1 plays an important role in metformin pharmacokinetics and pharmacodynamics.

OBJECTIVE

We performed a systematic review and meta-analysis to investigate the association between SLC47A1 rs2289669 (G > A) and the metformin response in drug-naive patients with type 2 diabetes.

METHODS

Studies published until January 27, 2022, were retrieved from Cochrane CENTRAL, Embase, PubMed, and Web of Science. Two reviewers independently screened titles, abstracts, and full-text articles. Studies conducted in newly diagnosed or drug-naive patients with type 2 diabetes who received metformin monotherapy were included. A total of 6 studies involving 953 patients were included in this meta-analysis. We extracted the study characteristics and changes in glycated hemoglobin (HbA1c) levels before and after treatment according to the SLC47A1 rs2289669 genotype. Changes in HbA1c levels were analyzed using mean differences (MDs) and 95% CIs. SLC47A1 rs2289669 was associated with changes in HbA1c levels (A carrier vs GG; MD = -0.55; 95% CI, -0.91 to - 0.20; I² = 63%). The sensitivity analysis yielded similar results to the main analysis (MD range, -0.64 to -0.37). When comparing all 3 genotypes, there were significant differences in HbA1c level changes between AA vs GG and GA vs GG, but not in GA vs AA.

CONCLUSION

This meta-analysis showed that SLC47A1 rs2289669 is associated with the glycemic response to metformin in drug-naive patients with type 2 diabetes.

A novel nonsense variant in ARID1B causing simultaneous RNA decay and exon skipping is associated with Coffin-Siris syndrome

Coffin-Siris syndrome (CSS) is a congenital disorder that is characterized by an absent/hypoplastic fifth distal phalanx, psychomotor developmental delay, and coarse facial features. One of the causative genes, ARID1B (AT-rich interactive domain-containing protein 1B), encodes components of the BAF chromatin remodeling complexes. Here, we report a case of a 3-year 8-month-old male with a novel nonsense variant (NM_001374820.1:c.4282C > T, p.(Gln1428*)) in the ARID1B gene, which was identified with whole-exome sequencing. He showed clinical symptoms of cleft soft palate, distinctive facial features (flat nasal bridge, thick eyebrows, and long eyelashes), right cryptorchidism, and hypertrichosis that partially overlapped with CSS. One of the most characteristic features of CSS is absent/hypoplastic fifth distal phalanx. He showed no obvious clinical finding in the lengths of his fingers or in the formation of his fingernails. However, radiographic analyses of the metacarpophalangeal bones revealed shortening of all the distal phalanges and fifth middle phalanges, suggesting brachydactyly. We performed mRNA analyses and revealed that both nonsense-mediated decay and nonsense-associated altered splicing were simultaneously caused by the c.4282C > T nonsense variant. The proband's clinical manifestations fit the previously reported criteria of disease for CSS or intellectual disability with ARID1B variant. Altogether, we suggest that c.4282C > T is a pathogenic variant that causes this clinical phenotype.

Deciphering associations between three RNA splicing-related genetic variants and lung cancer risk

Limited efforts have been made in assessing the effect of genome-wide profiling of RNA splicing-related variation on lung cancer risk. In the present study, we first identified RNA splicing-related genetic variants linked to lung cancer in a genome-wide profiling analysis and then conducted a two-stage (discovery and replication) association study in populations of European ancestry. Discovery and validation were conducted sequentially with a total of 29,266 cases and 56,450 controls from both the Transdisciplinary Research in Cancer of the Lung and the International Lung Cancer Consortium as well as the OncoArray database. For those variants identified as significant in the two datasets, we further performed stratified analyses by smoking status and histological type and investigated their effects on gene expression and potential regulatory mechanisms. We identified three genetic variants significantly associated with lung cancer risk: rs329118 in JADE2 (P = 8.80E-09), rs2285521 in GGA2 (P = 4.43E-08), and rs198459 in MYRF (P = 1.60E-06). The combined effects of all three SNPs were more evident in lung squamous cell carcinomas (P = 1.81E-08, P = 6.21E-08, and P = 7.93E-04, respectively) than in lung adenocarcinomas and in ever smokers (P = 9.80E-05, P = 2.70E-04, and P = 2.90E-05, respectively) than in never smokers. Gene expression quantitative trait analysis suggested a role for the SNPs in regulating transcriptional expression of the corresponding target genes. In conclusion, we report that three RNA splicing-related genetic variants contribute to lung cancer susceptibility in European populations. However, additional validation is needed, and specific splicing mechanisms of the target genes underlying the observed associations also warrants further exploration.

Genetic screening in patients of Meige syndrome and blepharospasm

OBJECTIVE

Meige syndrome (MS) is cranial dystonia, including bilateral eyelid spasms (blepharospasm; BSP) and involuntary movements of the jaw muscles (oromandibular dystonia; OMD). Up to now, the pathogenic genes of MS and BSP are still unclear.

METHODS

We performed Sanger sequencing of GNAL, TOR1A, TOR2A, THAP1, and REEP4 exons on 78 patients, including 53 BSP and 25 MS and 96 healthy controls.

RESULTS

c.845G > C[R282P] of TOR1A, c.629delC[p.Gly210AlafsTer60] of TOR2A, c.1322A > G[N441S] of GNAL, c.446G > A[R149Q], and c.649C > T[R217C] of REEP4 were identified and predicated as deleterious probably damaging variants. Three potential alterations of splicing variants of TOR1A and TOR2A were identified in patients. The frequencies of TOR1A rs1435566780 and THAP1 rs545930392 were higher in patients than in controls.

CONCLUSIONS

TOR1A rs1435566780 (c.*16G > C(G > A)) and THAP1 rs545930392 (c.192G > A[K64K]) may contribute to the etiology of MS and BSP. Other identified rare mutations predicted as deleterious probably damaging need further confirmation. Larger MS and BSP cohorts and functional studies will need to be performed further to elucidate the association between these genes and the diseases.

SERPING1 Variants and C1-INH Biological Function: A Close Relationship With C1-INH-HAE

Hereditary angioedema with C1 Inhibitor deficiency (C1-INH-HAE) is caused by a constellation of variants of the SERPING1 gene (n = 809; 1,494 pedigrees), accounting for 86.8% of HAE families, showing a pronounced mutagenic liability of SERPING1 and pertaining to 5.6% de novo variants. C1-INH is the major control serpin of the kallikrein–kinin system (KKS). In addition, C1-INH controls complement C1 and plasminogen activation, both systems contributing to inflammation. Recognizing the failed control of C1s protease or KKS provides the diagnosis of C1-INH-HAE. SERPING1 variants usually behave in an autosomal-dominant character with an incomplete penetrance and a low prevalence. A great majority of variants (809/893; 90.5%) that were introduced into online database have been considered as pathogenic/likely pathogenic. Haploinsufficiency is a common feature in C1-INH-HAE where a dominant-negative variant product impacts the wild-type allele and renders it inactive. Small (36.2%) and large (8.3%) deletions/duplications are common, with exon 4 as the most affected one. Point substitutions with missense variants (32.2%) are of interest for the serpin structure–function relationship. Canonical splice sites can be affected by variants within introns and exons also (14.3%). For noncanonical sequences, exon skipping has been confirmed by splicing analyses of patients' blood-derived RNAs (n = 25). Exonic variants (n = 6) can affect exon splicing. Rare deep-intron variants (n = 6), putatively acting as pseudo-exon activating mutations, have been characterized as pathogenic. Some variants have been characterized as benign/likely benign/of uncertain significance (n = 74). This category includes some homozygous (n = 10) or compound heterozygous variants (n = 11). They are presenting with minor allele frequency (MAF) below 0.00002 (i.e., lower than C1-INH-HAE frequency), and may be quantitatively unable to cause haploinsufficiency. Rare benign variants could contribute as disease modifiers. Gonadal mosaicism in C1-INH-HAE is rare and must be distinguished from a de novo variant. Situations with paternal or maternal disomy have been recorded (n = 3). Genotypes must be interpreted with biological investigation fitting with C1-INH expression and typing. Any SERPING1 variant reminiscent of the dysfunctional phenotype of serpin with multimerization or latency should be identified as serpinopathy.

The Contribution of COL4A5 Splicing Variants to the Pathogenesis of X-Linked Alport Syndrome

X-linked Alport syndrome (XLAS) is caused by pathogenic variants in COL4A5 and is characterized by progressive kidney disease, hearing loss, and ocular abnormalities. Recent advances in genetic analysis and further understanding of genotype-phenotype correlations in affected male patients raises the importance of detecting splicing variants in COL4A5. Aberrant splicing of COL4A5 is caused not only by canonical splice site variants but also non-canonical splice site variants such as deep intronic changes or even substitutions in exons. Patients with splicing variants account for ~15% of all cases in XLAS. In addition, it has been shown that there is a significant difference in kidney survival depending on the aberrant splicing patterns of transcripts- in particular in-frame or out-of-frame nucleotide changes in transcripts. Therefore, cDNA analysis of patient mRNA is necessary to determine the impact of splice site variants and to confirm a diagnosis of XLAS and to predict the kidney prognosis. However, it is usually difficult to amplify COL4A5 transcripts extracted from peripheral blood leukocytes. For these cases, in vitro minigene assays or RNA sequence extracted from urine derived cells can confirm aberrant splicing patterns. Moreover, controlling aberrant splicing by nucleic acids or small molecular compounds in genetic diseases are attracting attention as a potential therapeutic strategy. Here, we review the frequency of splicing variants in COL4A5, the latest diagnostic strategies, and the prospects for new therapeutic approaches.

Genomic architecture of phenotypic extremes in a wild cervid

Identifying the genes underlying fitness-related traits such as body size and male ornamentation can provide tools for conservation and management and are often subject to various selective pressures. Here we performed high-depth whole genome re-sequencing of pools of individuals representing the phenotypic extremes for antler and body size in white-tailed deer (Odocoileus virginianus). Samples were selected from a tissue repository containing phenotypic data for 4,466 male white-tailed deer from Anticosti Island, Quebec, with four pools representing the extreme phenotypes for antler and body size after controlling for age. Our results revealed a largely homogenous population but detected highly divergent windows between pools for both traits, with the mean allele frequency difference of 14% for and 13% for antler and body SNPs in outlier windows, respectively. Genes in outlier antler windows were enriched for pathways associated with cell death and protein metabolism and some of the most differentiated windows included genes associated with oncogenic pathways and reproduction, processes consistent with antler evolution and growth. Genes associated with body size were more nuanced, suggestive of a highly complex trait. Overall, this study revealed the complex genomic make-up of both antler morphology and body size in free-ranging white-tailed deer and identified target loci for additional analyses.

Biomedical data, computational methods and tools for evaluating disease-disease associations

In recent decades, exploring potential relationships between diseases has been an active research field. With the rapid accumulation of disease-related biomedical data, a lot of computational methods and tools/platforms have been developed to reveal intrinsic relationship between diseases, which can provide useful insights to the study of complex diseases, e.g. understanding molecular mechanisms of diseases and discovering new treatment of diseases. Human complex diseases involve both external phenotypic abnormalities and complex internal molecular mechanisms in organisms. Computational methods with different types of biomedical data from phenotype to genotype can evaluate disease-disease associations at different levels, providing a comprehensive perspective for understanding diseases. In this review, available biomedical data and databases for evaluating disease-disease associations are first summarized. Then, existing computational methods for disease-disease associations are reviewed and classified into five groups in terms of the usages of biomedical data, including disease semantic-based, phenotype-based, function-based, representation learning-based and text mining-based methods. Further, we summarize software tools/platforms for computation and analysis of disease-disease associations. Finally, we give a discussion and summary on the research of disease-disease associations. This review provides a systematic overview for current disease association research, which could promote the development and applications of computational methods and tools/platforms for disease-disease associations.

Genetic Factors of Renin–Angiotensin System Associated with Major Bleeding for Patients Treated with Direct Oral Anticoagulants

The purpose of this study was to identify the renin–angiotensin system (RAS)-related genetic factors associated with bleeding and develop the bleeding risk scoring system in patients receiving direct oral anticoagulants (DOACs). This study was a retrospective analysis of prospectively collected samples from June 2018 to May 2020. To investigate the associations between RAS-related genetic factors and major bleeding, we selected 16 single nucleotide polymorphisms (SNPs) from five genes (namely, AGT, REN, ACE, AGTR1, and AGTR2). Multivariable logistic regression analysis was employed to investigate the independent risk factors for bleeding and to develop a risk scoring system. A total of 172 patients were included in the analysis, including 33 major bleeding cases. Both old age (≥65 years) and moderate to severe renal impairment (CrCl < 50 mL/min) increased the risk of bleeding in the multivariable analysis. Among RAS-related polymorphisms, patients carrying TT genotype of rs5050 and A allele of rs4353 experienced a 3.6-fold (95% CI: 1.4–9.3) and 3.1-fold (95% CI: 1.1–9.3) increase in bleeding, respectively. The bleeding risk increased exponentially with a higher score; the risks were 0%, 2.8%, 16.9%, 32.7%, and 75% in patients with 0, 1, 2, 3, and 4 points, respectively. Although this study is limited to a retrospective study design, this is the first study to suggest RAS-related genetic markers and risk scoring systems, including both clinical and genetic factors, for major bleeding in patients receiving DOAC treatment.

Alternative Splicing in Human Biology and Disease

Alternative pre-mRNA splicing allows for the production of multiple mRNAs from an individual gene, which not only expands the protein-coding potential of the genome but also enables complex mechanisms for the post-transcriptional control of gene expression. Regulation of alternative splicing entails a combinatorial interplay between an abundance of trans-acting splicing factors, cis-acting regulatory sequence elements and their concerted effects on the core splicing machinery. Given the extent and biological significance of alternative splicing in humans, it is not surprising that aberrant splicing patterns can cause or contribute to a wide range of diseases. In this introductory chapter, we outline the mechanisms that govern alternative pre-mRNA splicing and its regulation and discuss how dysregulated splicing contributes to human diseases affecting the motor system and the brain.

Vulnerable exons, like ACADM exon 5, are highly dependent on maintaining a correct balance between splicing enhancers and silencers

It is now widely accepted that aberrant splicing of constitutive exons is often caused by mutations affecting cis-acting splicing regulatory elements, but there is a misconception that all exons have an equal dependency on splicing regulatory elements and thus a similar susceptibility to aberrant splicing. We investigated exonic mutations in ACADM exon 5 to experimentally examine their effect on splicing and found that 7 out of 11 tested mutations affected exon inclusion, demonstrating that this constitutive exon is particularly vulnerable to exonic splicing mutations. Employing ACADM exon 5 and 6 as models, we demonstrate that the balance between splicing enhancers and silencers, flanking intron length, and flanking splice site strength are important factors that determine exon definition and splicing efficiency of the exon in question. Our study shows that two constitutive exons in ACADM have different inherent vulnerabilities to exonic splicing mutations. This suggests that in silico prediction of potential pathogenic effects on splicing from exonic mutations may be improved by also considering the inherent vulnerability of the exon. Moreover, we show that single nucleotide polymorphism that affect either of two different exonic splicing silencers, located far apart in exon 5, all protect against both immediately flanking and more distant exonic splicing mutations.

Alternative splicing and liver disease

Alternative splicing produces complex and dynamic changes in the protein isoforms that are necessary for the proper biological functioning of the metabolic pathways involved in liver development and hepatocyte homeostasis. Changes in the physiological state of alternatively spliced forms are increasingly linked to liver pathologies. This may occur when the expression or function of the set of proteins controlling the alternative splicing processes are altered by external effectors such as oxidative stress and other environmental variations. Studies addressing these modifications reveal a complex interplay between the expression levels of different proteins that regulate the alternative splicing process as well as the changes in alternative splicing. This interplay results in a cascade of different protein isoforms that correlate with the progression of non-alcoholic fatty liver disease, hepatocellular carcinoma, and alcoholic liver disease. However, research on the detailed molecular mechanism underlying the production of these isoforms is needed. It is imperative to identify the physiological processes affected by the differentially spliced isoforms and confirm their role on the onset and maintenance of the pathology. This is required to design potential therapeutic approaches targeting the key splicing changes to revert the pathological condition as well as identify prognostic markers. In this review, we describe the complexity of the splicing process through an example to encourage researchers to go down this path. Subsequently, rather than a catalog of splicing events we have hand-picked and discuss a few selected studies of specific liver pathologies and suggested ways to focus research on these areas.

C Deletion at the re74650330 Locus of the SLC39A8 Gene (rs74650330) Increases the Risk of Coronary Artery Disease in Individuals with Low-Density Lipoprotein Cholesterol Levels

Background: Genetic variants of the SLC39A8 gene are associated with several cardiovascular disease risk factors, including body mass index, systolic blood pressure (SBP), diastolic blood pressure (DBP), N-terminal pro-B-type natriuretic peptide (NT-proBNP) and high-density lipoprotein cholesterol (HDL-C) levels. The present study aimed to investigate the association between the SLC39A8 SNPs rs13107325 and rs74650330 and CAD in the Han population in Jiangsu (China). Methods: Genotyping of these SNPs was performed in 258 patients with CAD and 170 healthy controls using the base-quenched probe technique. The association between the alleles of the rs74650330 locus and blood lipid and glucose profiles was investigated. Receiver operating characteristic (ROC) curve analysis was used to quantify the optimal thresholds for lipid and FBG levels and the risk factors for CAD were estimated by logistic regression analysis. Results: The rs13107325 polymorphism was not found in the 428 Chinese individuals enrolled in the current study. For rs74650330, individuals harboring the C allele had significantly higher HDL levels than those without this allele in the control group (p = 0.039), while the opposite was true for low-density lipoprotein cholesterol (LDL-C) levels (p = 0.046). Further analysis indicated that when LDL-C levels were lower than 2.365 mmol/L, subjects with C/del and del/del had a 7.293-fold increased risk of CAD compared with that of controls without the mutation (odds ratio: 7.293; 95% confidence interval: 0.953-55.79). Conclusions: The susceptibility of SLC39A8 polymorphisms to CAD were studied and revealed a possible role for the deletion variant of rs74650330 in increasing the risk of CAD among the Chinese Han population.

mRNA Analysis of Frameshift Mutations with Stop Codon in the Last Exon: The Case of Hemoglobins Campania [α1 cod95 (-C)] and Sciacca [α1 cod109 (-C)]

An insertion or deletion of a nucleotide (nt) in the penultimate or the last exon can result in a frameshift and premature termination codon (PTC), giving rise to an unstable protein variant, showing a dominant phenotype. We described two α-globin mutants created by the deletion of a nucleotide in the penultimate or the last exon of the α1-globin gene: the Hb Campania or α1 cod95 (−C), causing a frameshift resulting in a PTC at codon 102, and the Hb Sciacca or α1 cod109 (−C), causing a frameshift and formation of a PTC at codon 133. The carriers showed α-thalassemia alterations (mild microcytosis with normal Hb A2) and lacked hemoglobin variants. The 3D model indicated the α-chain variants’ instability, due to the severe structural alterations with impairment of the chaperone alpha-hemoglobin stabilizing protein (AHSP) interaction. The qualitative and semiquantitative analyses of the α1mRNA from the reticulocytes of carriers highlighted a reduction in the variant cDNAs that constituted 34% (Hb Campania) and 15% (Hb Sciacca) of the total α1-globin cDNA, respectively. We developed a workflow for the in silico analysis of mechanisms triggering no-go decay, and its results suggested that the reduction in the variant mRNA was likely due to no-go decay caused by the presence of a rare triplet, and, in the case of Hb Sciacca, also by the mRNA’s secondary structure variation. It would be interesting to correlate the phenotype with the quantity of other frameshift mRNA variants, but very few data concerning α- and β-globin variants are available.

Knockin mouse model of the human CFL2 p.A35T mutation results in a unique splicing defect and severe myopathy phenotype

Cofilin-2 is an actin-binding protein that is predominantly expressed in skeletal and cardiac muscles and belongs to the AC group of proteins, which includes cofilin-1 and destrin. In humans, cofilin-2 (CFL2) mutations have been associated with congenital myopathies that include nemaline and myofibrillar myopathy. To understand the pathogenicity of the human CFL2 mutation, p.A35T, that first linked cofilin-2 with the human disease, we created a knock-in mouse model. The Cfl2A35T/A35T (KI) mice were indistinguishable from their wild-type littermates at birth, but they rapidly worsened and died by postnatal day 9. The phenotypic, histopathologic and molecular findings mimicked the constitutive Cfl2-knockout (KO) mice described previously, including sarcomeric disruption and actin accumulations in skeletal muscles and negligible amounts of cofilin-2 protein. In addition, KI mice demonstrated a marked reduction in Cfl2 mRNA levels in various tissues including skeletal muscles. Further investigation revealed evidence of alternative splicing with the presence of two alternate transcripts of smaller size. These alternate transcripts were expressed at very low levels in the wild-type mice and were significantly upregulated in the mutant mice, indicating that pre-translational splicing defects may be a critical component of the disease mechanism associated with the mutation. Evidence of reduced expression of the full-length CFL2 transcript was also observed in the muscle biopsy sample of the patient with p.A35T mutation.

Dissection of pleiotropic effects of variants in and adjacent to F8 exon 19 and rescue of mRNA splicing and protein function

The pathogenic significance of nucleotide variants commonly relies on nucleotide position within the gene, with exonic changes generally attributed to quantitative or qualitative alteration of protein biosynthesis, secretion, activity, or clearance. However, these changes may exert pleiotropic effects on both protein biology and mRNA splicing due to the overlapping of the amino acid and splicing codes, thus shaping the disease phenotypes. Here, we focused on hemophilia A, in which the definition of F8 variants' causative role and association to bleeding phenotypes is crucial for proper classification, genetic counseling, and management of affected individuals. We extensively characterized a large panel of hemophilia A-causing variants (n = 30) within F8 exon 19 by combining and comparing in silico and recombinant expression analyses. We identified exonic variants with pleiotropic effects and dissected the altered protein features of all missense changes. Importantly, results from multiple prediction algorithms provided qualitative results, while recombinant assays allowed us to correctly infer the likely phenotype severity for 90% of variants. Molecular characterization of pathogenic variants was also instrumental for the development of tailored correction approaches to rescue splicing affecting variants or missense changes impairing protein folding. A single engineered U1snRNA rescued mRNA splicing of nine different variants and the use of a chaperone-like drug resulted in improved factor VIII protein secretion for four missense variants. Overall, dissection of the molecular mechanisms of a large panel of HA variants allowed precise classification of HA-affected individuals and favored the development of personalized therapeutic approaches.

The USDA-ARS Ag100Pest Initiative: High-Quality Genome Assemblies for Agricultural Pest Arthropod Research

Simple Summary

High-quality genome assemblies are essential tools for modern biological research. In the past, creating genome assemblies was prohibitively expensive and time-consuming for most non-model insect species due to, in part, the technical challenge of isolating the necessary quantity and quality of DNA from many species. Sequencing methods have now improved such that many insect genomes can be sequenced and assembled at scale. We created the Ag100Pest Initiative to propel agricultural research forward by assembling reference-quality genomes of important arthropod pest species. Here, we describe the Ag100Pest Initiative’s processes and experimental procedures. We show that the Ag100Pest Initiative will greatly expand the diversity of publicly available arthropod genome assemblies. We also demonstrate the high quality of preliminary contig assemblies. We share arthropod-specific technical details and insights that we have gained during the project. The methods and preliminary results presented herein should help other researchers attain similarly high-quality assemblies, effectively changing the landscape of insect genomics.

Abstract

The phylum Arthropoda includes species crucial for ecosystem stability, soil health, crop production, and others that present obstacles to crop and animal agriculture. The United States Department of Agriculture’s Agricultural Research Service initiated the Ag100Pest Initiative to generate reference genome assemblies of arthropods that are (or may become) pests to agricultural production and global food security. We describe the project goals, process, status, and future. The first three years of the project were focused on species selection, specimen collection, and the construction of lab and bioinformatics pipelines for the efficient production of assemblies at scale. Contig-level assemblies of 47 species are presented, all of which were generated from single specimens. Lessons learned and optimizations leading to the current pipeline are discussed. The project name implies a target of 100 species, but the efficiencies gained during the project have supported an expansion of the original goal and a total of 158 species are currently in the pipeline. We anticipate that the processes described in the paper will help other arthropod research groups or other consortia considering genome assembly at scale.

usDSM: a novel method for deleterious synonymous mutation prediction using undersampling scheme

Although synonymous mutations do not alter the encoded amino acids, they may impact protein function by interfering with the regulation of RNA splicing or altering transcript splicing. New progress on next-generation sequencing technologies has put the exploration of synonymous mutations at the forefront of precision medicine. Several approaches have been proposed for predicting the deleterious synonymous mutations specifically, but their performance is limited by imbalance of the positive and negative samples. In this study, we firstly expanded the number of samples greatly from various data sources and compared six undersampling strategies to solve the problem of the imbalanced datasets. The results suggested that cluster centroid is the most effective scheme. Secondly, we presented a computational model, undersampling scheme based method for deleterious synonymous mutation (usDSM) prediction, using 14-dimensional biology features and random forest classifier to detect the deleterious synonymous mutation. The results on the test datasets indicated that the proposed usDSM model can attain superior performance in comparison with other state-of-the-art machine learning methods. Lastly, we found that the deep learning model did not play a substantial role in deleterious synonymous mutation prediction through a lot of experiments, although it achieves superior results in other fields. In conclusion, we hope our work will contribute to the future development of computational methods for a more accurate prediction of the deleterious effect of human synonymous mutation. The web server of usDSM is freely accessible at http://usdsm.xialab.info/.

A verified genomic reference sample for assessing performance of cancer panels detecting small variants of low allele frequency

BACKGROUND

Oncopanel genomic testing, which identifies important somatic variants, is increasingly common in medical practice and especially in clinical trials. Currently, there is a paucity of reliable genomic reference samples having a suitably large number of pre-identified variants for properly assessing oncopanel assay analytical quality and performance. The FDA-led Sequencing and Quality Control Phase 2 (SEQC2) consortium analyze ten diverse cancer cell lines individually and their pool, termed Sample A, to develop a reference sample with suitably large numbers of coding positions with known (variant) positives and negatives for properly evaluating oncopanel analytical performance.

RESULTS

In reference Sample A, we identify more than 40,000 variants down to 1% allele frequency with more than 25,000 variants having less than 20% allele frequency with 1653 variants in COSMIC-related genes. This is 5-100× more than existing commercially available samples. We also identify an unprecedented number of negative positions in coding regions, allowing statistical rigor in assessing limit-of-detection, sensitivity, and precision. Over 300 loci are randomly selected and independently verified via droplet digital PCR with 100% concordance. Agilent normal reference Sample B can be admixed with Sample A to create new samples with a similar number of known variants at much lower allele frequency than what exists in Sample A natively, including known variants having allele frequency of 0.02%, a range suitable for assessing liquid biopsy panels.

CONCLUSION

These new reference samples and their admixtures provide superior capability for performing oncopanel quality control, analytical accuracy, and validation for small to large oncopanels and liquid biopsy assays.

GALNT2 Gene Variant rs4846914 Is Associated with Insulin and Insulin Resistance Depending on BMI in PCOS Patients: a Case-Control Study

GALNT2 is a GalNAc transferase that regulates insulin signaling, lipogenesis, and serum lipid fractions. The objective of this study was to investigate the association of GALNT2 rs2144300 and rs4846914 single nucleotide polymorphisms (SNPs) with the risk of polycystic ovary syndrome (PCOS) and related traits. The two SNPs were genotyped in 616 PCOS patients and 482 control subjects. Genetic associations with related traits were also analyzed. The genotype distributions of the two SNPs in PCOS patients were similar to those of normal controls. However, significant differences were noted across the three groups of genotypes with respect to the examined variables. In the PCOS group, subjects with genotype AA at the rs4846914 SNP exhibited an increased fasting serum insulin and homeostasis model insulin resistance (HOMA-IR) index compared with that of corresponding GG or GA genotype carriers (all P < 0.05). When PCOS patients were further separated into obese and non-obese subgroups, the genotype-related effects on insulin and HOMA-IR were more obvious, and variations in BMI and FSH levels were exclusively observed in obese PCOS subjects (all P < 0.05). In addition, fasting plasma glucose levels were affected by the genotypes of the rs2144300 SNP in normal control women (P < 0.05). rs4846914 and rs2144300 polymorphisms in the GALNT2 gene are associated with insulin and HOMA-IR, BMI, and FSH levels in obese PCOS patients and glucose levels in normal control women, respectively, but not with PCOS. GALNT2 rs4846914 AA carrier status may be associated with insulin resistance and related traits in obese patients.

Quantitative Evaluation of CFTR Pre-mRNA Splicing Dependent on the (TG)mTn Poly-Variant Tract

Genetic analysis in cystic fibrosis (CF) is a difficult task. Within the many causes of variability and uncertainty, a major determinant is poor knowledge of the functional effect of most DNA variants of the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene. In turn, knowledge of the effect of a CFTR variant has dramatic diagnostic, prognostic and, in the era of CF precision medicine, also therapeutic consequences. One of the most challenging CFTR variants is the (TG)mTn haplotype, which has variable functional effect and controversial clinical consequences. The exact quantification of the anomalous splicing of CFTR exon 10 (in the HGVS name; exon 9 in the legacy name) and, consequently, of the residual wild-type functional CFTR mRNA, should be mandatory in clinical assessment of patients with potentially pathological haplotype of this tract. Here, we present a real time-based assay for the quantification of the proportion of exon 10+/exon 10− CFTR mRNA, starting from nasal brushing. Our assay proved rapid, economic and easy to perform. Specific primers used for this assay are either disclosed or commercially available, allowing any laboratory to easily perform it. A simplified analysis of the data is provided, facilitating the interpretation of the results. This method helps to enhance the comprehension of the genotype–phenotype relationship in CF and CFTR-related disorders (CFTR-RD), crucial for the diagnosis, prognosis and personalized therapy of CF.