Expression Changes Confirm Genomic Variants Predicted to Result in Allele-Specific, Alternative mRNA Splicing

Deciphering the role of alternative splicing in neoplastic diseases for immune-oncological therapies

Alternative splicing (AS) is an important molecular biological mechanism regulated by complex mechanisms involving a plethora of cis and trans-acting elements. Furthermore, AS is tissue specific and altered in various pathologies, including infectious, inflammatory, and neoplastic diseases. Recently developed immuno-oncological therapies include monoclonal antibodies (mAbs) and chimeric antigen receptor (CAR) T cells targeting, among others, immune checkpoint (ICP) molecules. Despite therapeutic successes have been demonstrated, only a limited number of patients showed long-term benefit from these therapies with tumor entity-related differential response rates were observed. Interestingly, splice variants of common immunotherapeutic targets generated by AS are able to completely escape and/or reduce the efficacy of mAb- and/or CAR-based tumor immunotherapies. Therefore, the analyses of splicing patterns of targeted molecules in tumor specimens prior to therapy might help correct stratification, thereby increasing therapy success by antibody panel selection and antibody dosages. In addition, the expression of certain splicing factors has been linked with the patients' outcome, thereby highlighting their putative prognostic potential. Outstanding questions are addressed to translate the findings into clinical application. This review article provides an overview of the role of AS in (tumor) diseases, its molecular mechanisms, clinical relevance, and therapy response.

An elevated level of interleukin-17A in a Senegalese malaria cohort is associated with rs8193038 IL-17A genetic variant

Malaria infection is a multifactorial disease partly modulated by host immuno-genetic factors. Recent evidence has demonstrated the importance of Interleukin-17 family proinflammatory cytokines and their genetic variants in host immunity. However, limited knowledge exists about their role in parasitic infections such as malaria. We aimed to investigate IL-17A serum levels in patients with severe and uncomplicated malaria and gene polymorphism's influence on the IL-17A serum levels. In this research, 125 severe (SM) and uncomplicated (UM) malaria patients and 48 free malaria controls were enrolled. IL-17A serum levels were measured with ELISA. PCR and DNA sequencing were used to assess host genetic polymorphisms in IL-17A. We performed a multivariate regression to estimate the impact of human IL-17A variants on IL-17A serum levels and malaria outcomes. Elevated serum IL-17A levels accompanied by increased parasitemia were found in SM patients compared to UM and controls (P < 0.0001). Also, the IL-17A levels were lower in SM patients who were deceased than in those who survived. In addition, the minor allele frequencies (MAF) of two IL-17A polymorphisms (rs3819024 and rs3748067) were more prevalent in SM patients than UM patients, indicating an essential role in SM. Interestingly, the heterozygous rs8193038 AG genotype was significantly associated with higher levels of IL-17A than the homozygous wild type (AA). According to our results, it can be concluded that the IL-17A gene rs8193038 polymorphism significantly affects IL-17A gene expression. Our results fill a gap in the implication of IL-17A gene polymorphisms on the cytokine level in a malaria cohort. IL-17A gene polymorphisms also may influence cytokine production in response to Plasmodium infections and may contribute to the hyperinflammatory responses during severe malaria outcomes.

L-type calcium channels and neuropsychiatric diseases: Insights into genetic risk variant-associated genomic regulation and impact on brain development

Recent human genetic studies have linked a variety of genetic variants in the CACNA1C and CACNA1D genes to neuropsychiatric and neurodevelopmental disorders. This is not surprising given the work from multiple laboratories using cell and animal models that have established that Cav1.2 and Cav1.3 L-type calcium channels (LTCCs), encoded by CACNA1C and CACNA1D, respectively, play a key role in various neuronal processes that are essential for normal brain development, connectivity, and experience-dependent plasticity. Of the multiple genetic aberrations reported, genome-wide association studies (GWASs) have identified multiple single nucleotide polymorphisms (SNPs) in CACNA1C and CACNA1D that are present within introns, in accordance with the growing body of literature establishing that large numbers of SNPs associated with complex diseases, including neuropsychiatric disorders, are present within non-coding regions. How these intronic SNPs affect gene expression has remained a question. Here, we review recent studies that are beginning to shed light on how neuropsychiatric-linked non-coding genetic variants can impact gene expression via regulation at the genomic and chromatin levels. We additionally review recent studies that are uncovering how altered calcium signaling through LTCCs impact some of the neuronal developmental processes, such as neurogenesis, neuron migration, and neuron differentiation. Together, the described changes in genomic regulation and disruptions in neurodevelopment provide possible mechanisms by which genetic variants of LTCC genes contribute to neuropsychiatric and neurodevelopmental disorders.

Case Report: A prenatal diagnosis of osteogenesis imperfecta in a patient with a novel pathogenic variant in COL1A2

Osteogenesis imperfecta is considered a rare genetic condition which is characterized by bone fragility. In 85% of cases, it is caused by mutations in COL1A1 and COL1A2 genes which are essential to produce type I collagen. We report the case of a female neonate delivered to a 27-year-old women at San Bartolomé Teaching Hospital with a family history of clavicle fracture. A prenatal control with ultrasound was performed to the mother at 29 weeks. A fetus with altered morphology and multiple fractures was found. Therefore, a prenatal diagnosis of osteogenesis imperfecta was performed. The neonate was born with a respiratory distress syndrome and an acyanotic congenital heart disease. Therefore, she remained in NICU until her death. We highlight the importance of prenatal diagnosis, genetic counseling and a multidisciplinary evaluation in this type of pathologies and report a new probably pathogenic variant in the COL1A2 gene detected by exomic sequencing in amniotic fluid.

Identification of a Novel Non-Canonical Splice-Site Variant in ABCD1

Cerebral adrenoleukodystrophy (CALD) is a fatal genetic disease characterized by rapid, devastating neurological decline, with a narrow curative treatment window in the early stage. Non-canonical splice-site (NCSS) variants can easily be missed during genomic DNA analyses, and only a few of them in ABCD1 have been explored. Here, we studied a Chinese patient with clinical features similar to those of early-stage CALD but with a negative molecular diagnosis and a sibling who had presumably died of CALD. Trio-based whole-exome sequencing (trio-WES) and RNA sequencing (RNA-Seq) revealed a novel hemizygote NCSS variant c.901-25_901-9 del in ABCD1 intron 1, resulting in a complex splicing pattern. The in vitro minigene assay revealed that the c.901-25_901-9 del construct contained two aberrant transcripts that caused skipping of exon 2 and a small 48-bp deletion on left of the same exon. We identified a novel NCSS variant, that extends the spectrum of the known ABCD1 variants, and demonstrated the pathogenicity of this gene variant. Our findings highlight the importance of combining RNA-Seq and WES techniques for prompt diagnosis of leukodystrophy with NCSS variants.

Congenital myopathy as a new phenotype caused by two undescribed variants in ASCC1 gene

We present a patient with congenital myopathy and an inborn epiphysiolysis of the ulna. Whole-exome sequencing analysis revealed two novel mutations in Activation Signal Cointegrator Complex 1 (ASCC1) gene in a compound heterozygous state-a splicing variant c.395-2A>G and a deletion of the first two coding exons. Homozygous and compound heterozygous LoF variants in ASCC1 gene lead to a severe phenotype of spinal muscular atrophy with congenital bone fractures 2 (SMABF2). All patients described to date presented with a severe muscular hypotony, inborn fractures, and passed away shortly after birth while our proband had moderate hypotony, no fractures, but epiphysiolysis and he was 3.5 years old at the time of examination. To explain the phenotype of our patient, we performed an RNA analysis of all family members. We discovered that the c.395-2A>G variant results in two aberrant mRNA isoforms. We also validated the deletion of two exons in ASCC1 gene that lead to the increased expression of this truncated transcript by 1.8 times. To investigate the possible impact of this deletion on the phenotype we predicted a new Kozak sequence in exon 4 that could lead to the formation of a truncated protein with shortened KH domain and a full RNA ligase-like domain. We suggest that this unexpectedly different phenotype of the proband with ASCC1-related disorder could be explained by the presence of the truncated protein with an increased expression.

Late Onset Subacute Profound Biotinidase Deficiency Caused by a Novel Homozygous Variant c.466-3T>G in the BTD Gene

Biotinidase deficiency (BD) is an autosomal recessive disorder caused by bi-allelic mutation in the BTD gene. Clinical manifestations in BD mainly depends on residual biotinidase enzyme activity, although there are some exceptions. Broadly BD disorders are classified as profound BD and partial BD. Further profound BD can be early onset, late onset, and sometimes may be asymptomatic. Clinically late-onset profound BD can present with spectrum of manifestations ranging from single organ to multiple organ involvement, typically affecting function of brain, eye, ear, and skin. Here, a first-born child to consanguineous parents with late-onset profound BD presenting with hyperventilation secondary to lactic acidosis, hypotonia, evolving spasticity, and abnormal neuroimaging findings caused by novel homozygous variant, c.466-3T>G in the BTD gene is reported.

Dissecting the phenotypic variability of osteogenesis imperfecta

Osteogenesis imperfecta (OI) is a heterogeneous family of collagen type I-related diseases characterized by bone fragility. OI is most commonly caused by single-nucleotide substitutions that replace glycine residues or exon splicing defects in the COL1A1 and COL1A2 genes that encode the α1(I) and α2(I) collagen chains. Mutant collagen is partially retained intracellularly, impairing cell homeostasis. Upon secretion, it assembles in disorganized fibrils, altering mineralization. OI is characterized by a wide range of clinical outcomes, even in the presence of identical sequence variants. Given the heterotrimeric nature of collagen I, its amino acid composition and the peculiarity of its folding, several causes may underlie the phenotypic variability of OI. A deep analysis of entries regarding glycine and splice site collagen substitution of the largest publicly available patient database reveals a higher risk of lethal phenotype for carriers of variants in α1(I) than in α2(I) chain. However, splice site variants are predominantly associated with lethal phenotype when they occur in COL1A2. In addition, lethality is increased when mutations occur in regions of importance for extracellular matrix interactions. Both extracellular and intracellular determinants of OI clinical severity are discussed in light of the findings from in vitro and in vivo OI models. Combined with meticulous tracking of clinical cases via a publicly available database, the available OI animal models have proven to be a unique tool to shed light on new modulators of phenotype determination for this rare heterogeneous disease.

Genetic association study of intron variants in the forkhead box protein P3 gene in Chinese patients diagnosed with cervical cancer

The aim of this study was to investigate the effects of forkhead box protein P3 (FOXP3) intron single nucleotide variants (SNVs) in high‐risk human papilloma virus (HR‐HPV) infection and cervical cancer (CC) malignant lesions. We performed FOXP3 genotyping in 350 patients with CC and 350 healthy controls using the ImLDR multiple single nucleotide polymorphism genotyping technology. The heterozygous mutation TC in rs2294021 decreased the risk of HR‐HPV infection and CC malignant lesions (TC vs. TT: OR = 0.71, 95% CI = 0.51–0.99); the dominant model TC+CC and allele C in rs2294021 decreased the risk of CC malignant lesions (TC+CC vs. TT: OR = 0.69, 95% CI = 0.50–0.95; C vs. T: OR = 0.78, 95% CI = 0.63–0.97). The heterozygous mutation GA, dominant model GA+AA and allele A in rs3761549 also decreased the risk of HR‐HPV infection and CC malignant lesions (GA vs. GG: OR = 0.70, 95% CI = 0.51–0.96; GA+AA vs. GG: OR = 0.69, 95% CI = 0.51–0.94; A vs. G: OR = 0.75, 95% CI = 0.58–0.96). Patients with CC and HR‐HPV infection carrying rs2294021 TC and rs3761549 GA had lower expression of FOXP3 protein. Haplotype analysis revealed that T‐C‐A decreased the risk of HR‐HPV infection. Furthermore, we found a significant association between immune cells infiltration and prognosis in patients with CC. Our findings demonstrated that rs2294021 and rs3761549 variants may protect against HR‐HPV and CC malignant lesions by downregulating FOXP3 and that FOXP3 was associated with immune cells infiltration, which affected the prognosis of CC.

A hypomorphic mutation in the mouse Csn1s1 gene generated by CRISPR/Cas9 pronuclear microinjection

Caseins are major milk proteins that have an evolutionarily conserved role in nutrition. Sequence variations in the casein genes affect milk composition in livestock species. Regulatory elements of the casein genes could be used to direct the expression of desired transgenes into the milk of transgenic animals. Dozens of casein alleles have been identified for goats, cows, sheep, camels and horses, and these sequence variants are associated with altered gene expression and milk protein content. Most of the known mutations affecting casein genes’ expression are located in the promoter and 3’-untranslated regions. We performed pronuclear microinjections with Cas9 mRNA and sgRNA against the first coding exon of the mouse Csn1s1 gene to introduce random mutations in the α-casein (Csn1s1) signal peptide sequence at the beginning of the mouse gene. Sanger sequencing of the founder mice identified 40 mutations. As expected, mutations clustered around the sgRNA cut site (3 bp from PAM). Most of the mutations represented small deletions (1–10 bp), but we detected several larger deletions as well (100–300 bp). Functionally most mutations led to gene knockout due to a frameshift or a start codon loss. Some of the mutations represented in-frame indels in the first coding exon. Of these, we describe a novel hypomorphic Csn1s1 (Csn1s1c.4-5insTCC) allele. We measured Csn1s1 protein levels and confirmed that the mutation has a negative effect on milk composition, which shows a 50 % reduction in gene expression and a 40–80 % decrease in Csn1s1 protein amount, compared to the wild-type allele. We assumed that mutation affected transcript stability or splicing by an unknown mechanism. This mutation can potentially serve as a genetic marker for low Csn1s1 expression.

The Association Between Single-Nucleotide Polymorphisms of Co-Stimulatory Genes Within Non-HLA Region and the Prognosis of Leukemia Patients With Hematopoietic Stem Cell Transplantation

To avoid graft rejection, the hematopoietic stem cells with matched classical human leukocyte antigen (HLA) alleles are the primary choice for clinical allogeneic transplantation. However, even if the fully HLA-matched hematopoietic stem cells are used for transplantation, some patients still have poor prognosis after hematopoietic stem cell transplantation (HSCT), suggesting that the HLA system was not the only determinant of the outcomes of HSCT. In this study, we investigated whether the single-nucleotide polymorphisms (SNPs) of the co-stimulatory genes within non-HLA regions were related to the outcomes of HSCT. The genomic DNAs of 163 patients who had acute leukemia and received HSCT and their respective donors were collected for analysis. Thirty-four SNPs located in the four co-stimulatory genes including cytotoxic T-lymphocyte associated protein 4 (CTLA4), CD28, tumor necrosis factor ligand superfamily 4 (TNFSF4), and programmed cell death protein 1 (PDCD1) were selected to explore their relationship with the adverse outcomes after transplantation, including mortality, cytomegalovirus infection, graft-versus-host disease, and relapse. Our results revealed that nine SNPs in the CTLA4 gene, five SNPs in the PDCD1 gene, two SNPs in the TNFSF4 gene, and four SNPs in the CD28 gene were significantly associated with the occurrence of adverse outcomes post-HSCT. These SNPs may play important roles in immune response to allografts post-HSCT and can be the targets for developing strategy to identify appropriate donors.

(Re-)activity in the caregiving situation: Genetic diversity within Oxytocin-Vasopressin Pathway is associated with salivary oxytocin and vasopressin concentrations in response to contact with a crying infant-simulator

Oxytocin (OT) and vasopressin (AVP) hormones as well as their receptors (OXTR and AVPR1a) have been deemed crucial for caregiving and sensitive responsiveness to infant cues. However, previous research on genetic polymorphisms and OT and AVP levels in the context of caregiving were sparse and have brought contradictory findings. The aim of this reported observational study was to examine the impact of genetic variations within genes related to OT and AVP signaling pathway on hormones levels' changes in response to the caregiving situation. A total of 221 adult intimate couples (110 childless, non-pregnant and 111 expectant couples) participated in three 10 min sessions, during which they were taking care of a crying life-like simulator. 30 min prior to the first session salivary samples to analyze basal OT and AVP, and polymorphisms in OXTR, AVPR1a and CD38 genes were collected. Subsequent OT and AVP levels were measured 15 min after each session. The two most frequently studied OXTR SNPs (rs53576 and rs2254298) had no or a minor impact on higher OT levels, which were linked to rs1042778, rs13316193, rs2228485, rs2268490, rs4686302 genotypes. AVP levels were affected by rs1042778, rs13316193, rs4686302 and rs237887. OT levels varied depending on the OT (rs2770378, rs4813625), CD38 (rs379686), and 5-HTR2A (rs6314) genotype. OT and AVP levels were also associated with rs6314 (5-HTR2A). AVP levels were linked to ESR1 (rs1884051) and SIM1 (rs3734354) variations. Shorter variants of RS3 and RS1 were associated with lower levels of AVP. In conclusion, analyzed polymorphisms were associated with both the level and changes in OT and AVP hormone levels in the standardized situation of caregiving reactions to infant crying.

Estimating partial-body ionizing radiation exposure by automated cytogenetic biodosimetry

PURPOSE

Inhomogeneous exposures to ionizing radiation can be detected and quantified with the dicentric chromosome assay (DCA) of metaphase cells. Complete automation of interpretation of the DCA for whole-body irradiation has significantly improved throughput without compromising accuracy, however, low levels of residual false positive dicentric chromosomes (DCs) have confounded its application for partial-body exposure determination.

MATERIALS AND METHODS

We describe a method of estimating and correcting for false positive DCs in digitally processed images of metaphase cells. Nearly all DCs detected in unirradiated calibration samples are introduced by digital image processing. DC frequencies of irradiated calibration samples and those exposed to unknown radiation levels are corrected subtracting this false positive fraction from each. In partial-body exposures, the fraction of cells exposed, and radiation dose can be quantified after applying this modification of the contaminated Poisson method.

RESULTS

Dose estimates of three partially irradiated samples diverged 0.2-2.5 Gy from physical doses and irradiated cell fractions deviated by 2.3%-15.8% from the known levels. Synthetic partial-body samples comprised of unirradiated and 3 Gy samples from 4 laboratories were correctly discriminated as inhomogeneous by multiple criteria. Root mean squared errors of these dose estimates ranged from 0.52 to 1.14 Gy² and from 8.1 to 33.3%² for the fraction of cells irradiated.

CONCLUSIONS

Automated DCA can differentiate whole- from partial-body radiation exposures and provides timely quantification of estimated whole-body equivalent dose.