A homemade MLPA assay detects known CTNS mutations and identifies a novel deletion in a previously unresolved cystinosis family

Infantile nephropatic cystinosis is a rare, recessive, and genetically homogeneous disorder impairing renal function. It is caused by mutations in CTNS. Several large copy number aberrations have been identified but, for the majority of these, heterozygous patients and carriers can not easily be identified. We therefore developed a multiplex ligation-dependent probe amplification assay targeting eight of the twelve CTNS exons. We show that this assay is valid in detecting known deletions in both the homozygous and heterozygous state. The application to a family previously found mutation-negative by conventional screening revealed a novel large deletion which, as the first of its kind, does not involve the coding region. We conclude that our assay represents a valid tool for further completing the CTNS mutation spectrum and for simplified carrier testing in cystinosis families harboring copy number mutations. More generally, our study exemplifies the use of synthetic, homemade MLPA probesets as cheap, efficient, and rapidly available screening tools for small genes and/or very rare diseases.     

Mosaicism for FMR1 gene full mutation and intermediate allele in a female foetus: A postzygotic retraction event

Fragile X syndrome is caused by the expansion of an unstable CGG repeat in the 5′UTR of FMR1 gene. The occurrence of mosaicism is not uncommon, especially in male patients, whereas in females it is not so often reported. Here we report a female foetus that was subject to prenatal diagnosis, because of her mother being a premutation carrier. The foetus was identified as being a mosaic for an intermediate allele and a full mutation of FMR1 gene, in the presence of a normal allele. The mosaic status was confirmed in three different tissues of the foetus – amniotic fluid, skin biopsy and blood – the last two obtained after pregnancy termination. Karyotype analysis and X-chromosome STR markers analysis do not support the mosaicism as inheritance of both maternal alleles. Oligonucleotide array-CGH excluded an imbalance that could contain the primer binding site with a different repeat size. The obtained results give compelling evidence for a postzygotic expansion mechanism where the foetus mosaic pattern originated from expansion of the mother's premutation into a full mutation and consequent regression to an intermediate allele in a proportion of cells. These events occurred in early embryogenesis before the commitment of cells into the different tissues, as the three tested tissues of the foetus have the same mosaic pattern. The couple has a son with Fragile X mental retardation syndrome and choose to terminate this pregnancy after genetic counselling.     

Assessment of phylogenetic inter-relationships in the genus Cymbidium (Orchidaceae) based on internal transcribed spacer region of rDNA

Sequence data obtained from nrITS region were used to assess phylogenetic inter-relationships and infrageneric classification of ten Cymbidium species collected from north-east India. The final aligned data matrix of combined ITS 1, 5.8S and ITS 2 yielded 684 characters. The ITS 1 and ITS 2 regions showed variable sequence lengths and G + C content (%). The 5.8S region was found to be more conserved (98.71%) followed by ITS 1 (86.12%) and ITS 2 (69.40%). ITS 2 recorded highest percentage of parsimony informative sites (7.46%), high sequence divergence with indels (24.63%), high number of transitions and transversions. ITS sequence data determined the phylogeny of Asiatic Cymbidiums with high bootstrap values. All three proposed subgenera could be distinguished clearly by all four (MP, ML, NJ, and BI) phylogenetic methods. This study validates the utility of ITS rDNA region as a reliable indicator of phylogenetic relationships, especially ITS 2 as probable DNA barcode at higher levels and can serve as an additional approach for identification of broader range of plant taxa especially orchids.     

Trans-channel interactions in batrachotoxin-modified rat skeletal muscle sodium channels: kinetic analysis of mutual inhibition between mu-conotoxin GIIIA derivatives and amine blockers

R13X derivatives of μ-conotoxin GIIIA bind externally to single sodium channels and block current incompletely with mean “blocked” durations of several seconds. We studied interactions between two classes of blockers (μ-conotoxins and amines) by steady state, kinetic analysis of block of BTX-modified Na channels in planar bilayers. The amines cause all-or-none block at a site internal to the selectivity filter. TPrA and DEA block single Na channels with very different kinetics. TPrA induces discrete, all-or-none, blocked events (mean blocked durations, ∼100 ms), whereas DEA produces a concentration-dependent reduction of the apparent single channel amplitude (“fast” block). These distinct modes of action allow simultaneous evaluation of block by TPrA and DEA, showing a classical, competitive interaction between them. The apparent affinity of TPrA decreases with increasing [DEA], based on a decrease in the association rate for TPrA. When an R13X μ-conotoxin derivative and one of the amines are applied simultaneously on opposite sides of the membrane, a mutually inhibitory interaction is observed. Dissociation constants, at +50 mV, for TPrA (∼4 mM) and DEA (∼30 mM) increase by ∼20%-50% when R13E (nominal net charge, +4) or R13Q (+5) is bound. Analysis of the slow blocking kinetics for the two toxin derivatives showed comparable decreases in affinity of the μ-conotoxins in the presence of an amine. Although this mutual inhibition seems to be qualitatively consistent with an electrostatic interaction across the selectivity filter, quantitative considerations raise questions about the mechanistic details of the interaction.     

Structural and phylogenetic comparison of napsin genes: The duplication,loss of function and human-specific pseudogenization of napsin B

Aspartic proteinases form a widely distributed protein superfamily, including cathepsin D, cathepsin E, pepsins, renin, BACE and napsin. Human napsin genes are located on human chromosome 19q13, which comprises napsin A and napsin B. Napsin B has been annotated as a pseudogene because it lacks an in-frame stop codon; its nascent chains are cotranslationally degraded. Until recently, there have been no studies concerning the molecular evolution of the napsin protein family in the human genome. In the present study, we investigated the evolution and gene organization of the napsin protein family. Napsin B orthologs are primarily distributed in primates, while napsin A orthologs are the only napsin genes in other species. The corresponding regions of napsin B in the available sequences from primate species contain an in-frame stop codon at a position equivalent to that of human napsin A. In addition, a rare single-nucleotide polymorphism (SNP) that creates a proper stop codon in human napsin B was identified using HapMap populations. Recombinant protein expression and three-dimensional comparative modeling revealed that napsin B exhibits residual activity toward synthetic aspartic protease substrates compared with napsin A, presumably through a napsin B-specific Arg287 residue. Thus, napsin B was duplicated from napsin A during the early stages of primate evolution, and the subsequent loss of napsin B function during primate evolution reflected ongoing human-specific napsin B pseudogenization.     

SHOX gene defects and selected dysmorphic signs in patients of idiopathic short stature and Léri-Weill dyschondrosteosis

The aim of the study was to analyze frequency of SHOX gene defects and selected dysmorphic signs in patients of both idiopathic short stature (ISS) and Léri-Weill dyschondrosteosis (LWD), all derived from the Czech population.Overall, 98 subjects were analyzed in the study. Inclusion criteria were the presence of short stature (− 2.0 SD), in combination with at least one of the selected dysmorphic signs for the ISS + group; and the presence of Madelung deformity, without positive karyotyping for the LWD + group. Each proband was analyzed by use of P018 MLPA kit, which covers SHOX and its regulatory sequences. Additionally, mutational analysis was done of the coding portions of the SHOX.Both extent and breakpoint localizations in the deletions/duplications found were quite variable. Some PAR1 rearrangements were detected, without obvious phenotypic association. In the ISS + group, MLPA analysis detected four PAR1 deletions associated with a SHOX gene defect, PAR1 duplication with an ambiguous effect, and two SHOX mutations (13.7%). In the LWD + group, MLPA analysis detected nine deletions in PAR1 region, with a deleterious effect on SHOX, first reported case of isolated SHOX enhancer duplication, and SHOX mutation (68.8%). In both ISS + and LWD + groups were positivity associated with a disproportionately short stature; in the ISS + group, in combination with muscular hypertrophy.It seems that small PAR1 rearrangements might be quite frequent in the population. Our study suggests disproportionateness, especially in combination with muscular hypertrophy, as relevant indicators of ISS to be the effect of SHOX defect.     

Polymorphisms and phenotypic analysis of cytochrome P450 2D6 in the Tibetan population

The cytochrome P450 2D6 enzyme (CYP2D6) metabolizes about 25% of prescribed drugs in the endoplasmic reticulum, and genetic polymorphisms in CYP2D6 can greatly affect its activity and lead to differences among individuals in drug efficacy and adverse drug reactions. To investigate genetic polymorphisms in CYP2D6 among Tibetan Chinese, we directly sequenced the whole gene in 96 unrelated, healthy Tibetans from The Tibet Autonomous Region of China and screened for genetic variants in the promoter, exons, introns, and 3′UTR. We detected fifty-one genetic polymorphisms in CYP2D6, and 16 of them are novel. The allele frequencies of CYP2D6*1, *2, *5, *10, *41, and *49 were 0.25, 0.43, 0.02, 0.29, 0.02, and 0.01, respectively. The frequency of CYP2D6*10, a putative poor-metabolizer allele, was lower in our sample population compared with that in the Han Chinese population (p < 0.001). In addition, haplotype analysis allowed 15 CYP2D6 haplotypes to be classified into three groups. In conclusion, our results provide basic information about CPY2D6 alleles in Tibetans and suggest that the enzymatic activities of CYP2D6 may differ among the diverse ethnic populations of China. Our results provide a basis for safer drug administration and better therapeutic treatment among Tibetans.     

The evolution and functional divergence of the beta-carotene oxygenase gene family in teleost fish—Exemplified by Atlantic salmon

In mammals, two carotenoid cleaving oxygenases are known; beta-carotene 15,15′-monooxygenase (BCMO1) and beta-carotene 9′,10′-oxygenase (BCO2). BCMO1 is a key enzyme in vitamin A synthesis by symmetrically cleaving beta-carotene into 2 molecules of all-trans-retinal, while BCO2 is responsible for asymmetric cleavage of a broader range of carotenoids. Here, we show that the Atlantic salmon beta-carotene oxygenase (bco) gene family contains 5 members, three bco2 and two bcmo1 paralogs. Using public sequence databases, multiple bco genes were also found in several additional teleost species. Phylogenetic analysis indicates that bco2a and bco2b originate from the teleost fish specific genome duplication (FSGD or 3R), while the third and more distant paralog, bco2 like, might stem from a prior duplication event in the teleost lineage. The two bcmo1 paralogs (bcmo1 and bcmo1 like) appear to be the result of an ancient duplication event that took place before the divergence of ray-finned (Actinopterygii) and lobe-finned fish (Sarcopterygii), with subsequent nonfunctionalization and loss of one Sarcopterygii paralog. Gene expression analysis of the bcmo1 and bco2 paralogs in Atlantic salmon reveals regulatory divergence with tissue specific expression profiles, suggesting that the beta-carotene oxygenase subtypes have evolved functional divergences. We suggest that teleost fish have evolved and maintained an extended repertoire of beta-carotene oxygenases compared to the investigated Sarcopterygii species, and hypothesize that the main driver behind this functional divergence is the exposure to a diverse set of carotenoids in the aquatic environment.     

Three novel homozygous mutations in the GNPTG gene that cause mucolipidosis type III gamma

Background

Mucolipidosis type III gamma (MLIII gamma) is an autosomal recessive disease caused by a mutation in the GNPTG gene, which encodes the γ subunit of the N-acetylglucosamine-1-phosphotransferase (GlcNAc-1-phosphotransferase). This protein plays a key role in the transport of lysosomal hydrolases to the lysosome.

Methods

Three Chinese children with typical skeletal abnormalities of MLIII were identified, who were from unrelated consanguineous families. After obtaining informed consent, genomic DNA was isolated from the patients and their parents. Direct sequencing of the GNPTG and GNPTAB genes was performed using standard PCR reactions.

Results

The three probands showed clinical features typical of MLIII gamma, such as joint stiffness and vertebral scoliosis without coarsened facial features. Mutation analysis of the GNPTG gene showed that three novel mutations were identified, two in exon seven [c.425G>A (p.Cys142Val)] and [c.515dupC (p.His172Profs27X)], and one in exon eight [c.609+1G>C]. Their parents were determined to be heterozygous carriers when compared to the reference sequence in GenBank on NCBI.

Conclusions

Mutation of the GNPTG gene is the cause of MLIII gamma in our patients. Our findings expand the mutation spectrum of the GNPTG gene and extend the knowledge of the phenotype–genotype correlation of the disease.     

An RNA electrophoretic mobility shift and mutational analysis of rnp-4f 5′-UTR intron splicing regulatory proteins in Drosophila reveals a novel new role for a dADAR protein isoform

Alternative splicing greatly enhances the diversity of proteins encoded by eukaryotic genomes, and is also important in gene expression control. In contrast to the great depth of knowledge as to molecular mechanisms in the splicing pathway itself, relatively little is known about the regulatory events behind this process. The 5′-UTR and 3′-UTR in pre-mRNAs play a variety of roles in controlling eukaryotic gene expression, including translational modulation, and nearly 4000 of the roughly 14,000 protein coding genes in Drosophila contain introns of unknown functional significance in their 5′-UTR. Here we report the results of an RNA electrophoretic mobility shift analysis of Drosophila rnp-4f 5′-UTR intron 0 splicing regulatory proteins. The pre-mRNA potential regulatory element consists of an evolutionarily-conserved 177-nt stem-loop arising from pairing of intron 0 with part of adjacent exon 2. Incubation of in vitro transcribed probe with embryo protein extract is shown to result in two shifted RNA–protein bands, and protein extract from a dADAR null mutant fly line results in only one shifted band. A mutated stem-loop in which the conserved exon 2 primary sequence is changed but secondary structure maintained by introducing compensatory base changes results in diminished band shifts. To test the hypothesis that dADAR plays a role in intron splicing regulation in vivo, levels of unspliced rnp-4f mRNA in dADAR mutant were compared to wild-type via real-time qRT-PCR. The results show that during embryogenesis unspliced rnp-4f mRNA levels fall by up to 85% in the mutant, in support of the hypothesis. Taken together, these results demonstrate a novel role for dADAR protein in rnp-4f 5′-UTR alternative intron splicing regulation which is consistent with a previously proposed model.     

Characterization of Oncorhynchus mykiss 5-hydroxyisourate hydrolase/transthyretin superfamily: Evolutionary and functional analyses

Teleosts have highly diverged genomes that resulted from whole genome duplication, which leads to an extensive diversity of paralogous genes. Transthyretin (TTR), an extracellular thyroid hormone (TH) binding protein, is thought to have evolved from an ancestral 5-hydroxyisourate hydrolase (HIUHase) by gene duplication at some stage of chordate evolution. To characterize the functions of proteins that arose from duplicated genes in teleosts, we investigated the phylogenetic relationship of teleost HIUHase and TTR aa sequences, the expression levels of Oncorhynchus mykiss HIUHase and TTR mRNA in various tissues and the biological activities of the O. mykiss re-HIUHase and re-TTR. Phylogenetic analysis of the teleost aa sequences revealed the presence of two HIUHase subfamilies, HIUHase 1 (which has an N-terminal peroxisomal targeting signal-2 [PTS2]) and HIUHase 2 (which does not have an N-terminal PTS2), and one TTR family. The tissue distributions of HIUHase 1 and TTR mRNA were similar in juvenile O. mykiss and the mRNA levels were highest in the liver. The O. mykiss re-HIUHase and re-TTR proteins were both 40–50 kDa homotetramers consisting of 14–15 kDa subunits, with 30% identity. HIUHase had 5-hydroxyisourate (5-HIU) hydrolysis activity with Zn^2 + sensitivity, whereas TTR had ligand binding activity with a preference for THs and several environmental chemicals, such as halogenated phenols. Our results suggest that O. mykiss HIUHase and TTR have diverged from a common ancestral HIHUase with no functional complementation.     

Defective pre-mRNA splicing in PKD1 due to presumed missense and synonymous mutations causing autosomal dominant polycystic disease

Autosomal dominant polycystic kidney disease is the most common human monogenic disorder and is caused by mutations in the PKD1 or PKD2 genes. Most patients with the disease present mutations in PKD1, and a considerable number of these alterations are single base substitutions within the coding sequence that are usually predicted to lead to missense or synonymous mutations. There is growing evidence that some of these mutations can be detrimental by affecting the pre-mRNA splicing process. The aim of our study was to test PKD1 mutations, described as missense or synonymous in the literature or databases, for their effects on exon inclusion. Bioinformatics tools were used to select mutations with a potential effect on pre-mRNA splicing. Mutations were experimentally tested using minigene assays. Exons and adjacent intronic sequences were PCR-amplified and cloned in the splicing reporter minigene, and selected mutations were introduced by site-directed mutagenesis. Minigenes were transfected into kidney derived cell lines. RNA from cultured cells was analyzed by RT-PCR and DNA sequencing. Analysis of thirty-three PKD1 exonic mutations revealed three mutations that induce splicing defects. The substitution c.11156G > A, previously predicted as missense mutation p.R3719Q, abolished the donor splice site of intron 38 and resulted in the incorporation of exon 38 with 117 bp of intron 38 and skipping of exon 39. Two synonymous variants, c.327A > T (p.G109G) and c.11257C > A (p.R3753R), generated strong donor splice sites within exons 3 and 39 respectively, resulting in incorporation of incomplete exons. These three nucleotide substitutions represent the first PKD1 exonic mutations that induce aberrant mRNAs. Our results strengthen the importance to evaluate the consequences of presumed missense and synonymous mutations at the mRNA level.