Whole exome sequencing identifies a novel splice-site mutation in ADAMTS17 in an Indian family with Weill-Marchesani syndrome

Purpose: Weill-Marchesani syndrome (WMS) is a rare connective tissue disorder, characterized by short stature, micro-spherophakic lens, and stubby hands and feet (brachydactyly). WMS is caused by mutations in the FBN1, ADAMTS10, and LTBP2 genes. Mutations in the LTBP2 and ADAMTS17 genes cause a WMS-like syndrome, in which the affected individuals show major features of WMS but do not display brachydactyly and joint stiffness. The main purpose of our study was to determine the genetic cause of WMS in an Indian family. Methods: Whole exome sequencing (WES) was used to identify the genetic cause of WMS in the family. The cosegregation of the mutation was determined with Sanger sequencing. Reverse transcription (RT)-PCR analysis was used to assess the effect of a splice-site mutation on splicing of the ADAMTS17 transcript. Results: The WES analysis identified a homozygous novel splice-site mutation c.873+1G>T in a known WMS-like syndrome gene, ADAMTS17, in the family. RT-PCR analysis in the patient showed that exon 5 was skipped, which resulted in the deletion of 28 amino acids in the ADAMTS17 protein. Conclusions: The mutation in the WMS-like syndrome gene ADAMTS17 also causes WMS in an Indian family. The present study will be helpful in genetic diagnosis of this family and increases the number of mutations of this gene to six.

A short in-frame deletion in NTRK1 tyrosine kinase domain caused by a novel splice site mutation in a patient with congenital insensitivity to pain with anhidrosis

Background: Congenital insensitivity to pain with anhidrosis (CIPA) is a rare autosomal recessive genetic disease characterized by the lack of reaction to noxious stimuli and anhidrosis. It is caused by mutations in the NTRK1 gene, which encodes the high affinity tyrosine kinase receptor I for Neurotrophic Growth Factor (NGF). -- Case Presentation: We present the case of a female patient diagnosed with CIPA at the age of 8 months. The patient is currently 6 years old and her psychomotor development conforms to her age (RMN, SPECT and psychological study are in the range of normality). PCR amplification of DNA, followed by direct sequencing, was used to investigate the presence of NTRK1 gene mutations. Reverse transcriptase (RT)-PCR amplification of RNA, followed by cloning and sequencing of isolated RT-PCR products was used to characterize the effect of the mutations on NTRK1 mRNA splicing. The clinical diagnosis of CIPA was confirmed by the detection of two splice-site mutations in NTRK1, revealing that the patient was a compound heterozygote at this gene. One of these alterations, c.574+1G > A, is located at the splice donor site of intron 5. We also found a second mutation, c.2206-2 A > G, not previously reported in the literature, which is located at the splice acceptor site of intron 16. Each parent was confirmed to be a carrier for one of the mutations by DNA sequencing analysis. It has been proposed that the c.574+1G > A mutation would cause exon 5 skipping during NTRK1 mRNA splicing. We could confirm this prediction and, more importantly, we provide evidence that the novel c.2206-2A > G mutation also disrupts normal NTRK1 splicing, leading to the use of an alternative splice acceptor site within exon 17. As a consequence, this mutation would result in the production of a mutant NTRK1 protein with a seven aminoacid in-frame deletion in its tyrosine kinase domain. --Conclusions: We present the first description of a CIPA-associated NTRK1 mutation causing a short interstitial deletion in the tyrosine kinase domain of the receptor. The possible phenotypical implications of this mutation are discussed.

Dominant Splice Site Mutations in PIK3R1 Cause Hyper IgM Syndrome, Lymphadenopathy and Short Stature.

The purpose of this research was to use next generation sequencing to identify mutations in patients with primary immunodeficiency diseases whose pathogenic gene mutations had not been identified. Remarkably, four unrelated patients were found by next generation sequencing to have the same heterozygous mutation in an essential donor splice site of PIK3R1 (NM_181523.2:c.1425 + 1G > A) found in three prior reports. All four had the Hyper IgM syndrome, lymphadenopathy and short stature, and one also had SHORT syndrome. They were investigated with in vitro immune studies, RT-PCR, and immunoblotting studies of the mutation's effect on mTOR pathway signaling. All patients had very low percentages of memory B cells and class-switched memory B cells and reduced numbers of naïve CD4+ and CD8+ T cells. RT-PCR confirmed the presence of both an abnormal 273 base-pair (bp) size and a normal 399 bp size band in the patient and only the normal band was present in the parents. Following anti-CD40 stimulation, patient's EBV-B cells displayed higher levels of S6 phosphorylation (mTOR complex 1 dependent event), Akt phosphorylation at serine 473 (mTOR complex 2 dependent event), and Akt phosphorylation at threonine 308 (PI3K/PDK1 dependent event) than controls, suggesting elevated mTOR signaling downstream of CD40. These observations suggest that amino acids 435-474 in PIK3R1 are important for its stability and also its ability to restrain PI3K activity. Deletion of Exon 11 leads to constitutive activation of PI3K signaling. This is the first report of this mutation and immunologic abnormalities in SHORT syndrome.

Autosomal dominant Alport syndrome caused by a COL4A3 splice site mutation

Alport syndrome (AS) is a clinically and genetically heterogeneous renal disorder, predominantly affecting the type IV collagen alpha 3/alpha 4/alpha 5 network of the glomerular basement membrane (GBM). AS can be caused by mutations in any of the three genes encoding these type IV collagen chains. The majority of AS families (85%) are X-linked (XL-AS) involving mutations in the COL4A5 gene. Mutations in the COL4A3 and COL4A4 genes cause autosomal recessive AS (AR-AS), accounting for approximately 14% of the cases. Recently, autosomal dominant AS (AD-AS) was linked to the COL4A3/COL4A4 locus in a large family.

Growth hormone (GH)-releasing hormone increases the expression of the dominant-negative GH isoform in cases of isolated GH deficiency due to GH splice-site mutations

An autosomal dominant form of isolated GH deficiency (IGHD II) can result from heterozygous splice site mutations that weaken recognition of exon 3 leading to aberrant splicing of GH-1 transcripts and production of a dominant-negative 17.5-kDa GH isoform. Previous studies suggested that the extent of missplicing varies with different mutations and the level of GH expression and/or secretion. To study this, wt-hGH and/or different hGH-splice site mutants (GH-IVS+2, GH-IVS+6, GH-ISE+28) were transfected in rat pituitary cells expressing human GHRH receptor (GC-GHRHR). Upon GHRH stimulation, GC-GHRHR cells coexpressing wt-hGH and each of the mutants displayed reduced hGH secretion and intracellular GH content when compared with cells expressing only wt-hGH, confirming the dominant-negative effect of 17.5-kDa isoform on the secretion of 22-kDa GH. Furthermore, increased amount of 17.5-kDa isoform produced after GHRH stimulation in cells expressing GH-splice site mutants reduced production of endogenous rat GH, which was not observed after GHRH-induced increase in wt-hGH. In conclusion, our results support the hypothesis that after GHRH stimulation, the severity of IGHD II depends on the position of splice site mutation leading to the production of increasing amounts of 17.5-kDa protein, which reduces the storage and secretion of wt-GH in the most severely affected cases. Due to the absence of GH and IGF-I-negative feedback in IGHD II, a chronic up-regulation of GHRH would lead to an increased stimulatory drive to somatotrophs to produce more 17.5-kDa GH from the severest mutant alleles, thereby accelerating autodestruction of somatotrophs in a vicious cycle.

Isolated autosomal dominant growth hormone deficiency: stimulating mutant GH-1 gene expression drives GH-1 splice-site selection, cell proliferation, and apoptosis

The majority of mutations that cause isolated GH deficiency type II (IGHD II) affect splicing of GH-1 transcripts and produce a dominant-negative GH isoform lacking exon 3 resulting in a 17.5-kDa isoform, which further leads to disruption of the GH secretory pathway. A clinical variability in the severity of the IGHD II phenotype depending on the GH-1 gene alteration has been reported, and in vitro and transgenic animal data suggest that the onset and severity of the phenotype relates to the proportion of 17.5-kDa produced. The removal of GH in IGHD creates a positive feedback loop driving more GH expression, which may itself increase 17.5-kDa isoform productions from alternate splice sites in the mutated GH-1 allele. In this study, we aimed to test this idea by comparing the impact of stimulated expression by glucocorticoids on the production of different GH isoforms from wild-type (wt) and mutant GH-1 genes, relying on the glucocorticoid regulatory element within intron 1 in the GH-1 gene. AtT-20 cells were transfected with wt-GH or mutated GH-1 variants (5'IVS-3 + 2-bp T->C; 5'IVS-3 + 6 bp T->C; ISEm1: IVS-3 + 28 G->A) known to cause clinical IGHD II of varying severity. Cells were stimulated with 1 and 10 mum dexamethasone (DEX) for 24 h, after which the relative amounts of GH-1 splice variants were determined by semiquantitative and quantitative (TaqMan) RT-PCR. In the absence of DEX, only around 1% wt-GH-1 transcripts were the 17.5-kDa isoform, whereas the three mutant GH-1 variants produced 29, 39, and 78% of the 17.5-kDa isoform. DEX stimulated total GH-1 gene transcription from all constructs. Notably, however, DEX increased the amount of 17.5-kDa GH isoform relative to the 22- and 20-kDa isoforms produced from the mutated GH-1 variants, but not from wt-GH-1. This DEX-induced enhancement of 17.5-kDa GH isoform production, up to 100% in the most severe case, was completely blocked by the addition of RU486. In other studies, we measured cell proliferation rates, annexin V staining, and DNA fragmentation in cells transfected with the same GH-1 constructs. The results showed that that the 5'IVS-3 + 2-bp GH-1 gene mutation had a more severe impact on those measures than the splice site mutations within 5'IVS-3 + 6 bp or ISE +28, in line with the clinical severity observed with these mutations. Our findings that the proportion of 17.5-kDa produced from mutant GH-1 alleles increases with increased drive for gene expression may help to explain the variable onset progression, and severity observed in IGHD II.

Sensitivity of splice sites to antisense oligonucleotides in vivo.

A series of HeLa cell lines which stably express beta-globin pre-mRNAs carrying point mutations at nt 654, 705, or 745 of intron 2 has been developed. The mutations generate aberrant 5' splice sites and activate a common 3' cryptic splice site upstream leading to aberrantly spliced beta-globin mRNA. Antisense oligonucleotides, which in vivo blocked aberrant splice sites and restored correct splicing of the pre-mRNA, revealed major differences in the sensitivity of these sites to antisense probes. Although the targeted pre-mRNAs differed only by single point mutations, the effective concentrations of the oligonucleotides required for correction of splicing varied up to 750-fold. The differences among the aberrant 5' splice sites affected sensitivity of both the 5' and 3' splice sites; in particular, sensitivity of both splice sites was severely reduced by modification of the aberrant 5' splice sites to the consensus sequence. These results suggest large differences in splicing of very similar pre-mRNAs in vivo. They also indicate that antisense oligonucleotides may provide useful tools for studying the interactions of splicing machinery with pre-mRNA.

Role of adenine functional groups in the recognition of the 3′-splice-site AG during the second step of pre-mRNA splicing

The AG dinucleotide at the 3′ splice sites of metazoan nuclear pre-mRNAs plays a critical role in catalytic step II of the splicing reaction. Previous studies have shown that replacement of the guanine by adenine in the AG (AG → GG) inhibits this step. We find that the second step was even more severely inhibited by cytosine (AG → CG) or uracil (AG → UG) substitutions at this position. By contrast, a relatively moderate inhibition was observed with a hypoxanthine substitution (AG → HG). When adenine was replaced by a purine base (AG → PG) or by 7-deazaadenine (AG → c7AG), little effect on the second step was observed, suggesting that the 6-NH2 and N7 groups do not play a critical role in adenine recognition. Finally, replacement of adenine by 2-aminopurine (AG → 2-APG) had no effect on the second step. Taken together, our results suggest that the N1 group of adenine functions as an essential determinant in adenine recognition during the second step of pre-mRNA splicing.

Conserved stem–loop structures in the HIV-1 RNA region containing the A3 3′ splice site and its cis-regulatory element: possible involvement in RNA splicing

The HIV-1 transcript is alternatively spliced to over 30 different mRNAs. Whether RNA secondary structure can influence HIV-1 RNA alternative splicing has not previously been examined. Here we have determined the secondary structure of the HIV-1/BRU RNA segment, containing the alternative A3, A4a, A4b, A4c and A5 3′ splice sites. Site A3, required for tat mRNA production, is contained in the terminal loop of a stem–loop structure (SLS2), which is highly conserved in HIV-1 and related SIVcpz strains. The exon splicing silencer (ESS2) acting on site A3 is located in a long irregular stem–loop structure (SLS3). Two SLS3 domains were protected by nuclear components under splicing condition assays. One contains the A4c branch points and a putative SR protein binding site. The other one is adjacent to ESS2. Unexpectedly, only the 3′ A residue of ESS2 was protected. The suboptimal A3 polypyrimidine tract (PPT) is base paired. Using site-directed mutagenesis and transfection of a mini-HIV-1 cDNA into HeLa cells, we found that, in a wild-type PPT context, a mutation of the A3 downstream sequence that reinforced SLS2 stability decreased site A3 utilization. This was not the case with an optimized PPT. Hence, sequence and secondary structure of the PPT may cooperate in limiting site A3 utilization.

GeneSplicer: a new computational method for splice site prediction

GeneSplicer is a new, flexible system for detecting splice sites in the genomic DNA of various eukaryotes. The system has been tested successfully using DNA from two reference organisms: the model plant Arabidopsis thaliana and human. It was compared to six programs representing the leading splice site detectors for each of these species: NetPlantGene, NetGene2, HSPL, NNSplice, GENIO and SpliceView. In each case GeneSplicer performed comparably to the best alternative, in terms of both accuracy and computational efficiency.

Modulation of 5' splice site choice in pre-messenger RNA by two distinct steps.

Ser/Arg-rich proteins (SR proteins) are essential splicing factors that commit pre-messenger RNAs to splicing and also modulate 5' splice site choice in the presence or absence of functional U1 small nuclear ribonucleoproteins (snRNPs). Here, we perturbed the U1 snRNP in HeLa cell nuclear extract by detaching the U1-specific A protein using a 2'-O-methyl oligonucleotide (L2) complementary to its binding site in U1 RNA. In this extract, the standard adenovirus substrate is spliced normally, but excess amounts of SR proteins do not exclusively switch splicing from the normal 5' splice site to a proximal site (site 125 within the adenovirus intron), suggesting that modulation of 5' splice site choice exerted by SR proteins requires integrity of the U1 snRNP. The observation that splicing does not necessarily follow U1 binding indicates that interactions between the U1 snRNP and components assembled on the 3' splice site via SR proteins may also be critical for 5' splice site selection. Accordingly, we found that SR proteins promote the binding of the U2 snRNP to the branch site and stabilize the complex formed on a 3'-half substrate in the presence or absence of functional U1 snRNPs. A novel U2/U6/3'-half substrate crosslink was also detected and promoted by SR proteins. Our results suggest that SR proteins in collaboration with the U1 snRNP function in two distinct steps to modulate 5' splice site selection.

Lethal chondrodysplasia in a family of Holstein cattle is associated with a de novo splice site variant of COL2A1

Background Lethal chondrodysplasia (bulldog syndrome) is a well-known congenital syndrome in cattle and occurs sporadically in many breeds. In 2015, it was noticed that about 12 % of the offspring of the phenotypically normal Danish Holstein sire VH Cadiz Captivo showed chondrodysplasia resembling previously reported bulldog calves. Pedigree analysis of affected calves did not display obvious inbreeding to a common ancestor, suggesting the causative allele was not a rare recessive. The normal phenotype of the sire suggested a dominant inheritance with incomplete penetrance or a mosaic mutation. Results Three malformed calves were examined by necropsy, histopathology, radiology, and computed tomography scanning. These calves were morphologically similar and displayed severe disproportionate dwarfism and reduced body weight. The syndrome was characterized by shortening and compression of the body due to reduced length of the spine and the long bones of the limbs. The vicerocranium had severe dysplasia and palatoschisis. The bones had small irregular diaphyses and enlarged epiphyses consisting only of chondroid tissue. The sire and a total of four affected half-sib offspring and their dams were genotyped with the BovineHD SNP array to map the defect in the genome. Significant genetic linkage was obtained for several regions of the bovine genome including chromosome 5 where whole genome sequencing of an affected calf revealed a COL2A1 point mutation (g.32473300 G > A). This private sequence variant was predicted to affect splicing as it altered the conserved splice donor sequence GT at the 5’-end of COL2A1 intron 36, which was changed to AT. All five available cases carried the mutant allele in heterozygous state and all five dams were homozygous wild type. The sire VH Cadiz Captivo was shown to be a gonadal and somatic mosaic as assessed by the presence of the mutant allele at levels of about 5 % in peripheral blood and 15 % in semen. Conclusions The phenotypic and genetic findings are comparable to a previously reported COL2A1 missense mutation underlying lethal chondrodysplasia in the offspring of a mosaic French Holstein sire (Igale Masc). The identified independent spontaneous splice site variant in COL2A1 most likely caused chondrodysplasia and must have occurred during the early foetal development of the sire. This study provides a first example of a dominant COL2A1 splice site variant as candidate causal mutation of a severe lethal chondrodysplasia phenotype. Germline mosaicism is a relatively frequent mechanism in the origin of genetic disorders and explains the prevalence of a certain fraction of affected offspring. Paternal dominant de novo mutations are a risk in cattle breeding, especially because the ratio of defective offspring may be very high and be associated with significant animal welfare problems.