Prediction of human mRNA donor and acceptor sites from the DNA sequence

Functional Analysis of A Novel Splicing Mutation in The Mutase Gene of Two Unrelated Pedigrees

Objective

Methylmalonic acidura (MMA) is a rare autosomal recessive inborn error of metabolism. In this study we present a novel nucleotide change in the mutase (MUT) gene of two unrelated Iranian pedigrees and introduce the methods used for its functional analysis.

Materials and Methods

Two probands with definite diagnosis of MMA and a common novel variant in the MUT were included in a descriptive study. Bioinformatic prediction of the splicing variant was done with different prediction servers. Reverse transcriptionpolymerase chain reaction (RT-PCR) was done for splicing analysis and the products were analyzed by sequencing.

Results

The included index patients showed elevated levels of propionylcarnitine (C3). Urine organic acid analysis confirmed the diagnosis of MMA, and screening for mutations in the MUT revealed a novel C to G variation at the 3´ splice acceptor site in intron 12. In silico analysis suggested the change as a mutation in a conserved sequence. The splicing analysis showed that the C to G nucleotide change at position -3 in the acceptor splice site can lead to retention of the intron 12 sequence.

Conclusion

This is the first report of a mutation at the position -3 in the MUT intron 12 (c.2125-3C>G). The results suggest that the identified variation can be associated with the typical clinical manifestations of MMA.

Reproductive neuropeptides that stimulate spawning in the Sydney Rock Oyster (Saccostrea glomerata)

The Sydney Rock Oyster, Saccostrea glomerata, is a socioeconomically important species in Australia, yet little is known about the molecular mechanism that regulates its reproduction. To address this gap, we have performed a combination of high throughput transcriptomic and peptidomic analysis, to identify genes and neuropeptides that are expressed in the key regulatory tissues of S. glomerata; the visceral ganglia and gonads. Neuropeptides are known to encompass a diverse class of peptide messengers that play functional roles in many aspects of an animal's life, including reproduction. Approximately 28 neuropeptide genes were identified, primarily within the visceral ganglia transcriptome, that encode precursor proteins containing numerous neuropeptides; some were confirmed through mass spectral peptidomics analysis of the visceral ganglia. Of those, 28 bioactive neuropeptides were synthesized, and then tested for their capacity to induce gonad development and spawning in S. glomerata. Egg laying hormone, gonadotropin-releasing hormone, APGWamide, buccalin, CCAP and LFRFamide were neuropeptides found to trigger spawning in ripe animals. Additional testing of APGWa and buccalin demonstrated their capacity to advance conditioning and gonadal maturation. In summary, our analysis of S. glomerata has identified neuropeptides that can influence the reproductive cycle of this species, specifically by accelerating gonadal maturation and triggering spawning. Other molluscan neuropeptides identified in this study will enable further research into understanding the neuroendocrinology of oysters, which may benefit their cultivation.

Identification of a female spawn-associated Kazal-type inhibitor from the tropical abalone Haliotis asinina

Abalone (Haliotis) undergoes a period of reproductive maturation, followed by the synchronous release of gametes, called broadcast spawning. Field and laboratory studies have shown that the tropical species Haliotis asinina undergoes a two-week spawning cycle, thus providing an excellent opportunity to investigate the presence of endogenous spawning-associated peptides. In female H. asinina, we have isolated a peptide (5145 Da) whose relative abundance in hemolymph increases substantially just prior to spawning and is still detected using reverse-phase high-performance liquid chromatography chromatograms up to 1-day post-spawn. We have isolated this peptide from female hemolymph as well as samples prepared from the gravid female gonad, and demonstrated through comparative sequence analysis that it contains features characteristic of Kazal-type proteinase inhibitors (KPIs). Has-KPI is expressed specifically within the gonad of adult females. A recombinant Has-KPI was generated using a yeast expression system. The recombinant Has-KPI does not induce premature spawning of female H. asinina when administered intramuscularly. However it displays homomeric aggregations and interaction with at least one mollusc-type neuropeptide (LRDFVamide), suggesting a role for it in regulating neuropeptide endocrine communication. This research provides new understanding of a peptide that can regulate reproductive processes in female abalone, which has the potential to lead to the development of greater control over abalone spawning. The findings also highlight the need to further explore abalone reproduction to clearly define a role for novel spawning-associated peptide in sexual maturation and spawning. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.

A novel MED12 mutation: Evidence for a fourth phenotype

Mutations of the MED12 gene have been reported mainly in males with FG (Opitz-Kaveggia), Lujan-Fryns, or X-linked Ohdo syndromes. Recently, a different phenotype characterized by minor anomalies, severe intellectual disability (ID), and absent language was reported in female and male patients belonging to the same family and carrying a frameshift MED12 mutation (c.5898dupC). Here, we report on two brothers and their niece affected by severe and mild ID, respectively, where whole exome sequencing combined with variant analysis within a panel of ID-related genes, disclosed a novel c.2312T>C (p.Ile771Thr) MED12 mutation. This variant, which has not been reported as a polymorphism, was not present in a third unaffected brother, and was predicted to be deleterious by five bioinformatic databases. This finding together with the phenotypic analogies shared with the carriers of c.5898dupC mutation suggests the existence of a fourth MED12-related disorder, characterized by severe ID, absent or deficient language and, milder, clinical manifestation in heterozygotes. We have reviewed the literature on MED12 heterozygotes, their clinical manifestations, and discuss the possible biological causes of this condition. © 2016 Wiley Periodicals, Inc.

CADD score has limited clinical validity for the identification of pathogenic variants in noncoding regions in a hereditary cancer panel

PURPOSE

Several in silico tools have been shown to have reasonable research sensitivity and specificity for classifying sequence variants in coding regions. The recently-developed Combined Annotation Dependent Depletion (CADD) method generates predictive scores for single nucleotide variants (SNVs) in all areas of the genome, including non-coding regions. We sought to determine the clinical validity of non-coding variant CADD scores.

METHODS

We evaluated 12,391 unique SNVs in 624 patient samples submitted for germline mutation testing in a cancer-related gene panel. We compared the distributions of CADD scores of rare SNVs, common SNVs in our patient population, and the null distribution of all possible SNVs stratifying by genomic region.

RESULTS

The median CADD scores of intronic and nonsynonymous variants were significantly different between rare and common SNVs (p<0.0001). Despite these different distributions, no individual variants could be identified as plausibly causative among rare intronic variants with the highest scores. The ROC AUC for non-coding variants is modest, and the positive predictive value of CADD for intronic variants in panel testing was found to be 0.088.

CONCLUSION

Focused in-silico scoring systems with much higher predictive value will be necessary for clinical genomic applications.

Expression and interaction of the mango ethylene receptor MiETR1 and different receptor versions of MiERS1

Different versions of the mango ethylene receptor MiERS1 were identified and the analysis indicates that, in addition to MiERS1, two short versions of this receptor (MiERS1m, MiERS1s), representing truncated proteins with central deletions of functional domains, are present in mango. The short receptor versions reveal a different expression pattern compared to MiERS1, and they are highly variably transcribed. With transient expression assays using fluorescent fusion proteins, the localisation and the interaction of the receptors were determined in leaf cells of the tobacco model. MiERS1, MiETR1, and the short MiERS1 receptor versions are anchored in the endoplasmic reticulum (ER) membrane and co-localise with each other and with an ER-marker. Furthermore, ectopic expression of the mango receptors appears to induce a re-organisation of the ER resulting in accumulation of ER bodies. Interaction assays suggest that both short MiERS1 receptor versions can bind to proteins located in the ER. Bi-molecular fluorescence complementation (BiFC) assays indicate, that MiERS1m may dimerise with itself and can also interact with MiERS1, but not with MiETR1. Further, it as found that MiETR1 can interact with MiERS1. Interaction of MiERS1s with the other ethylene receptors could not be detected, although it was located in the ER membrane system.

Frequency of rare mutations and common genetic variations in severe hypertriglyceridemia in the general population of Spain

Background

Hypertriglyceridemia (HTG) is a common complex metabolic trait that results of the accumulation of relatively common genetic variants in combination with other modifier genes and environmental factors resulting in increased plasma triglyceride (TG) levels. The majority of severe primary hypertriglyceridemias is diagnosed in adulthood and their molecular bases have not been fully defined yet. The prevalence of HTG is highly variable among populations, possibly caused by differences in environmental factors and genetic background. However, the prevalence of very high TG and the frequency of rare mutations causing HTG in a whole non-selected population have not been previously studied.

Methods

The total of 23,310 subjects over 18 years from a primary care-district in a middle-class area of Zaragoza (Spain) with TG >500 mg/dL were selected to establish HTG prevalence. Those affected of primary HTG were considered for further genetic analisys. The promoters, coding regions and exon-intron boundaries of LPL, LMF1, APOC2, APOA5, APOE and GPIHBP1 genes were sequenced. The frequency of rare variants identified was studied in 90 controls.

Results

One hundred ninety-four subjects (1.04 %) had HTG and 90 subjects (46.4 %) met the inclusion criteria for primary HTG. In this subgroup, nine patients (12.3 %) were carriers of 7 rare variants in LPL, LMF1, APOA5, GPIHBP1 or APOE genes. Three of these mutations are described for the first time in this work. The presence of a rare pathogenic mutation did not confer a differential phenotype or a higher family history of HTG.

Conclusion

The prevalence of rare mutations in candidate genes in subjects with primary HTG is low. The low frequency of rare mutations, the absence of a more severe phenotype or the dominant transmission of the HTG would not suggest the use of genetic analysis in the clinical practice in this population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12944-016-0251-2) contains supplementary material, which is available to authorized users.

Single nucleotide polymorphism analysis reveals heterogeneity within a seedling tree population of a polyembryonic mango cultivar

Within a polyembryonic mango seedling tree population, the genetic background of individuals should be identical because vigorous plants for cultivation are expected to develop from nucellar embryos representing maternal clones. Due to the fact that the mango cultivar 'Hôi' is assigned to the polyembryonic ecotype, an intra-cultivar variability of ethylene receptor genes was unexpected. Ethylene receptors in plants are conserved, but the number of receptors or receptor isoforms is variable regarding different plant species. However, it is shown here that the ethylene receptor MiETR1 is present in various isoforms within the mango cultivar 'Hôi'. The investigation of single nucleotide polymorphisms revealed that different MiETR1 isoforms can not be discriminated simply by individual single nucleotide exchanges but by the specific arrangement of single nucleotide polymorphisms at certain positions in the exons of MiETR1. Furthermore, an MiETR1 isoform devoid of introns in the genomic sequence was identified. The investigation demonstrates some limitations of high resolution melting and ScreenClust analysis and points out the necessity of sequencing to identify individual isoforms and to determine the variability within the tree population.

GNAO1 encephalopathy: further delineation of a severe neurodevelopmental syndrome affecting females

Background

De novo heterozygous mutations in the GNAO1 gene, encoding the Gα o subunit of G-proteins, are the cause of a severe neurodevelopmental disorder, featuring early infantile seizures, profound cognitive dysfunction and, occasionally, movement disorder (early infantile epileptic encephalopathy-17).

Methods

We report a further case of this association in a 20 month-old Spanish girl with neonatal-onset refractory seizures, progressive microcephaly, oral-lingual dyskinesia and nearly absent psychomotor development. We performed whole-exome sequencing, a computational structural analysis of the novel gene variant identified and reviewed the previously reported cases.

Results

Trio whole-exome-sequencing uncovered a de novo p.Leu199Pro GNAO1 mutation. Computational structural analysis indicates this novel variant adversely affects the stability of the G-protein heterotrimeric complex as a whole. Of note, our patient showed a sustained seizure reduction while on a ketogenic diet.

Conclusions

With this observation, a total of twelve patients with GNAO1 encephalopathy have been reported. Oral-lingual dyskinesia and responsiveness of seizures to ketogenic diet are novel features. The distorted sex ratio (12/12 females) of the condition remains unexplained; a differential gender effect of the disruption of G-protein- mediated signal transduction on the developing brain can be hypothesized.

Electronic supplementary material

The online version of this article (doi:10.1186/s13023-016-0416-0) contains supplementary material, which is available to authorized users.

Correlating Photoreceptor Mosaic Structure to Clinical Findings in Stargardt Disease

Purpose

To demonstrate a method for correlating photoreceptor mosaic structure with optical coherence tomography (OCT) and microperimetry findings in patients with Stargardt disease.

Methods

A total of 14 patients with clinically diagnosed Stargardt disease were imaged using confocal and split-detection adaptive optics scanning light ophthalmoscopy. Cone photoreceptors were identified manually in a band along the temporal meridian. Resulting values were compared to a normative database (n = 9) to generate cone density deviation (CDD) maps. Manual measurement of outer nuclear layer plus Henle fiber layer (ONL+HFL) thickness was performed, in addition to determination of the presence of ellipsoid zone (EZ) and interdigitation zone (IZ) bands on OCT. These results, along with microperimetry data, were overlaid with the CDD maps.

Results

Wide variation in foveal structure and CDD maps was seen within this small group. Disruption of ONL+HFL and/or IZ band was seen in all patients, with EZ band preservation in regions with low cone density in 38% of locations analyzed. Normality of retinal lamellar structure on OCT corresponded with cone density and visual function at 50/78 locations analyzed. Outer retinal tubulations containing photoreceptor-like structures were observed in 3 patients.

Conclusions

The use of CDD color-coded maps enables direct comparison of cone mosaic local density with other measures of retinal structure and function. Larger normative datasets and improved tools for automation of image alignment are needed.

Translational Relevance

The approach described facilitates comparison of complex multimodal data sets from patients with inherited retinal degeneration, and can be expanded to incorporate other structural imaging or functional testing.

MYBPH acts as modifier of cardiac hypertrophy in hypertrophic cardiomyopathy (HCM) patients

Left ventricular hypertrophy is a risk factor for cardiovascular morbidity and mortality. Hypertrophic cardiomyopathy (HCM) is considered a model disease to study causal molecular factors underlying isolated cardiac hypertrophy. However, HCM manifests with various clinical symptoms, even in families bearing the same genetic defects, suggesting that additional factors contribute to hypertrophy. The gene encoding the cardiac myosin binding protein C (cMYBPC) is one of the most frequently implicated genes in HCM. Recently another myosin binding protein, myosin binding protein H (MYBPH) was shown to function in concert with cMYBPC in regulating cardiomyocyte contraction. Given the similarity in sequence, structure and the critical role MYBPH plays in sarcomere contraction, we proposed that MYBPH may be involved in HCM pathogenesis. Family-based genetic association analysis was employed to investigate the contribution of MYBPH in modifying hypertrophy. Seven single nucleotide polymorphisms and haplotypes in MYBPH were investigated for hypertrophy modifying effects in 388 individuals (27 families), in which three unique South African HCM-causing founder mutations (p.R403W and pA797T in β-myosin heavy chain gene (MYH7) and p.R92W in the cardiac troponin T gene (TNNT2)) segregate. We observed a significant association between rs2250509 and hypertrophy traits in the p.A797T MYH7 mutation group. Additionally, haplotype GGTACTT significantly affected hypertrophy within the same mutation group. MYBPH was for the first time assessed as a candidate hypertrophy modifying gene. We identified a novel association between MYBPH and hypertrophy traits in HCM patients carrying the p.A797T MYH7 mutation, suggesting that variation in MYBPH can modulate the severity of hypertrophy in HCM.

CELF4 Variant and Anthracycline-Related Cardiomyopathy: A Children's Oncology Group Genome-Wide Association Study

PURPOSE

Interindividual variability in the dose-dependent association between anthracyclines and cardiomyopathy suggests that genetic susceptibility could play a role. The current study uses an agnostic approach to identify genetic variants that could modify cardiomyopathy risk.

METHODS

A genome-wide association study was conducted in childhood cancer survivors with and without cardiomyopathy (cases and controls, respectively). Single-nucleotide polymorphisms (SNPs) that surpassed a prespecified threshold for statistical significance were independently replicated. The possible mechanistic significance of validated SNP(s) was sought by using healthy heart samples.

RESULTS

No SNP was marginally associated with cardiomyopathy. However, SNP rs1786814 on the CELF4 gene passed the significance cutoff for gene-environment interaction (Pge = 1.14 × 10(-5)). Multivariable analyses adjusted for age at cancer diagnosis, sex, anthracycline dose, and chest radiation revealed that, among patients with the A allele, cardiomyopathy was infrequent and not dose related. However, among those exposed to greater than 300 mg/m(2) of anthracyclines, the rs1786814 GG genotype conferred a 10.2-fold (95% CI, 3.8- to 27.3-fold; P < .001) increased risk of cardiomyopathy compared with those who had GA/AA genotypes and anthracycline exposure of 300 mg/m(2) or less. This gene-environment interaction was successfully replicated in an independent set of anthracycline-related cardiomyopathy. CUG-BP and ETR-3-like factor proteins control developmentally regulated splicing of TNNT2, the gene that encodes for cardiac troponin T (cTnT), a biomarker of myocardial injury. Coexistence of more than one cTnT variant results in a temporally split myofilament response to calcium, which causes decreased contractility. Analysis of TNNT2 splicing variants in healthy human hearts suggested an association between the rs1786814 GG genotype and coexistence of more than one TNNT2 splicing variant (90.5% GG v 41.7% GA/AA; P = .005).

CONCLUSION

We report a modifying effect of a polymorphism of CELF4 (rs1786814) on the dose-dependent association between anthracyclines and cardiomyopathy, which possibly occurs through a pathway that involves the expression of abnormally spliced TNNT2 variants.

Adipokinetic hormone receptor gene identification and its role in triacylglycerol metabolism in the blood-sucking insect Rhodnius prolixus

Adipokinetic hormone (AKH) has been associated with the control of energy metabolism in a large number of arthropod species due to its role on the stimulation of lipid, carbohydrate and amino acid mobilization/release. In the insect Rhodnius prolixus, a vector of Chagas' disease, triacylglycerol (TAG) stores must be mobilized to sustain the metabolic requirements during moments of exercise or starvation. Besides the recent identification of the R. prolixus AKH peptide, other components required for the AKH signaling cascade and its mode of action remain uncharacterized in this insect. In the present study, we identified and investigated the expression profile of the gene encoding the AKH receptor of R. prolixus (RhoprAkhr). This gene is highly conserved in comparison to other sequences already described and its transcript is abundant in the fat body and the flight muscle of the kissing bug. Moreover, RhoprAkhr expression is induced in the fat body at moments of increased TAG mobilization; the knockdown of this gene resulted in TAG accumulation both in fat body and flight muscle after starvation. The inhibition of Rhopr-AKHR transcription as well as the treatment of insects with the peptide Rhopr-AKH in its synthetic form altered the transcript levels of two genes involved in lipid metabolism, the acyl-CoA-binding protein-1 (RhoprAcbp1) and the mitochondrial glycerol-3-phosphate acyltransferase-1 (RhoprGpat1). These results indicate that the AKH receptor is regulated at transcriptional level and is required for TAG mobilization under starvation. In addition to the classical view of AKH as a direct regulator of enzymatic activity, we propose here that AKH signaling may account for the regulation of nutrient metabolism by affecting the expression profile of target genes.

Trafficking protein particle complex 6A delta (TRAPPC6AΔ) is an extracellular plaque-forming protein in the brain

Tumor suppressor WWOX is involved in the progression of cancer and neurodegeneration. Here, we examined whether protein aggregation occurs in the brain of nondemented, middle-aged humans and whether this is associated with WWOX downregulation. We isolated an N-terminal internal deletion isoform, TPC6AΔ, derived from alternative splicing of the TRAPPC6A (TPC6A) gene transcript. TPC6AΔ proteins are present as aggregates or plaques in the extracellular matrix of the brain such as in the cortex. Filter retardation assays revealed that aggregate formation of TPC6AΔ occurs preceding Aβ generation in the hippocampi of middle-aged postmortem normal humans. In a Wwox gene knockout mouse model, we showed the plaques of pT181-Tau and TPC6AΔ in the cortex and hippocampus in 3-week-old mice, suggesting a role of WWOX in limiting TPC6AΔ aggregation. To support this hypothesis, in vitro analysis revealed that TGF-β1 induces dissociation of the ectopic complex of TPC6AΔ and WWOX in cells, and then TPC6AΔ undergoes Ser35 phosphorylation-dependent polymerization and induces caspase 3 activation and Aβ production. Similarly, knockdown of WWOX by siRNA resulted in dramatic aggregation of TPC6AΔ. Together, when WWOX is downregulated, TPC6AΔ is phosphorylated at Ser35 and becomes aggregated for causing caspase activation that leads to Tau aggregation and Aβ formation.

Molecular genetic analysis of PKHD1 by next-generation sequencing in Czech families with autosomal recessive polycystic kidney disease

Background

Autosomal recessive polycystic kidney disease (ARPKD) is an early-onset form of polycystic kidney disease that often leads to devastating outcomes for patients. ARPKD is caused by mutations in the PKHD1 gene, an extensive gene that encodes for the ciliary protein fibrocystin/polyductin. Next-generation sequencing is presently the best option for molecular diagnosis of ARPKD. Our aim was to set up the first study of ARPKD patients from the Czech Republic, to determine the composition of their mutations and genotype-phenotype correlations, along with establishment of next-generation sequencing of the PKHD1 gene that could be used for the diagnosis of ARPKD patients.

Methods

Mutational analysis of the PKHD1 gene was performed in 24 families using the amplicon-based next-generation sequencing (NGS) technique. In patients without 2 causal mutations identified by NGS, subsequent MLPA analysis of the PKHD1 gene was carried out.

Results

Two underlying mutations were detected in 54 % of families (n = 13), one mutation in 13 % of families (n = 3), and in 33 % of families (n = 8) no mutation could be detected. Overall, seventeen different mutations (5 novel) were detected, including deletion of one exon. The detection rate in our study reached 60 % in the entire cohort of patients; but 90 % in the group of patients who fulfilled all clinical criteria of ARPKD, and 42 % in the group of patients with unknown kidney pathology. The most frequent mutation was T36M, accounting for nearly 21 % of all identified mutations.

Conclusions

Next-generation sequencing of the PKHD1 gene is a very useful method of molecular diagnosis in patients with a full clinical picture of ARPKD, and it has a high detection rate. Furthermore, its relatively low costs and rapidity allow the molecular genetic analysis of patients without the full clinical criteria of ARPKD, who might also have mutations in the PKHD1 gene.

Development of viable TAP-tagged dengue virus for investigation of host-virus interactions in viral replication

Dengue virus (DENV) is a mosquito-borne flavivirus responsible for life-threatening dengue haemorrhagic fever (DHF) and dengue shock syndrome (DSS). The viral replication machinery containing the core non-structural protein 5 (NS5) is implicated in severe dengue symptoms but molecular details remain obscure. To date, studies seeking to catalogue and characterize interaction networks between viral NS5 and host proteins have been limited to the yeast two-hybrid system, computational prediction and co-immunoprecipitation (IP) of ectopically expressed NS5. However, these traditional approaches do not reproduce a natural course of infection in which a number of DENV NS proteins colocalize and tightly associate during the replication process. Here, we demonstrate the development of a recombinant DENV that harbours a TAP tag in NS5 to study host-virus interactions in vivo. We show that our engineered DENV was infective in several human cell lines and that the tags were stable over multiple viral passages, suggesting negligible structural and functional disturbance of NS5. We further provide proof-of-concept for the use of rationally tagged virus by revealing a high confidence NS5 interaction network in human hepatic cells. Our analysis uncovered previously unrecognized hnRNP complexes and several low-abundance fatty acid metabolism genes, which have been implicated in the viral life cycle. This study sets a new standard for investigation of host-flavivirus interactions.

Identification of novel PKD1 and PKD2 mutations in a Chinese population with autosomal dominant polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) is one of the most frequently inherited renal diseases caused by mutations in PKD1 and PKD2. We performed mutational analyses of PKD genes in 49 unrelated patients using direct PCR-sequencing and multiplex ligation-dependent probe amplification (MLPA) for PKD1 and PKD2. RT-PCR analysis was also performed in a family with a novel PKD2 splicing mutation. Disease-causing mutations were identified in 44 (89.8%) of the patients: 42 (95.5%) of the patients showed mutations in PKD1, and 2 (4.5%) showed mutations in PKD2. Ten nonsense, 17 frameshift, 4 splicing and one in-frame mutation were found in 32 of the patients. Large rearrangements were found in 3 patients, and missense mutations were found in 9 patients. Approximately 61.4% (27/44) of the mutations are first reported with a known mutation rate of 38.6%. RNA analysis of a novel PKD2 mutation (c.595_595 + 14delGGTAAGAGCGCGCGA) suggested monoallelic expression of the wild-type allele. Furthermore, patients with PKD1-truncating mutations reached end-stage renal disease (ESRD) earlier than patients with non-truncating mutations (47 ± 3.522 years vs. 59 ± 11.687 years, P = 0.016). The mutation screening of PKD genes in Chinese ADPKD patients will enrich our mutation database and significantly contribute to improve genetic counselling for ADPKD patients.

From Whole Gene Deletion to Point Mutations of EP300-Positive Rubinstein-Taybi Patients: New Insights into the Mutational Spectrum and Peculiar Clinical Hallmarks

Rubinstein-Taybi syndrome (RSTS) is a rare congenital neurodevelopmental disorder characterized by growth deficiency, skeletal abnormalities, dysmorphic features, and intellectual disability. Causative mutations in CREBBP and EP300 genes have been identified in ∼55% and ∼8% of affected individuals. To date, only 28 EP300 alterations in 29 RSTS clinically described patients have been reported. EP300 analysis of 22 CREBBP-negative RSTS patients from our cohort led us to identify six novel mutations: a 376-kb deletion depleting EP300 gene; an exons 17-19 deletion (c.(3141+1_3142-1)_(3590+1_3591-1)del/p.(Ile1047Serfs*30)); two stop mutations, (c.3829A>T/p.(Lys1277*) and c.4585C>T/p.(Arg1529*)); a splicing mutation (c.1878-12A>G/p.(Ala627Glnfs*11)), and a duplication (c.4640dupA/p.(Asn1547Lysfs*3)). All EP300-mutated individuals show a mild RSTS phenotype and peculiar findings including maternal gestosis, skin manifestation, especially nevi or keloids, back malformations, and a behavior predisposing to anxiety. Furthermore, the patient carrying the complete EP300 deletion does not show a markedly severe clinical picture, even if a more composite phenotype was noticed. By characterizing six novel EP300-mutated patients, this study provides further insights into the EP300-specific clinical presentation and expands the mutational repertoire including the first case of a whole gene deletion. These new data will enhance EP300-mutated cases identification highlighting distinctive features and will improve the clinical practice allowing a better genotype-phenotype correlation.

Genetic characterization of interleukins (IL-1α, IL-1β, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18) with relevant biological roles in lagomorphs

ILs, as essential innate immune modulators, are involved in an array of biological processes. In the European rabbit (Oryctolagus cuniculus) IL-1α, IL-1β, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18 have been implicated in inflammatory processes and in the immune response against rabbit hemorrhagic disease virus and myxoma virus infections. In this study we characterized these ILs in six Lagomorpha species (European rabbit, pygmy rabbit, two cottontail rabbit species, European brown hare and American pika). Overall, these ILs are conserved between lagomorphs, including in their exon/intron structure. Most differences were observed between leporids and American pika. Indeed, when comparing both, some relevant differences were observed in American pika, such as the location of the stop codon in IL-1α and IL-2, the existence of a different transcript in IL8 and the number of cysteine residues in IL-1β. Changes at N-glycosylation motifs were also detected in IL-1, IL-10, IL-12B and IL-15. IL-1α is the protein that presents the highest evolutionary distances, which is in contrast to IL-12A where the distances between lagomorphs are the lowest. For all these ILs, sequences of human and European rabbit are more closely related than between human and mouse or European rabbit and mouse.

Survey of Programs Used to Detect Alternative Splicing Isoforms from Deep Sequencing Data In Silico

Next-generation sequencing techniques have been rapidly emerging. However, the massive sequencing reads hide a great deal of unknown important information. Advances have enabled researchers to discover alternative splicing (AS) sites and isoforms using computational approaches instead of molecular experiments. Given the importance of AS for gene expression and protein diversity in eukaryotes, detecting alternative splicing and isoforms represents a hot topic in systems biology and epigenetics research. The computational methods applied to AS prediction have improved since the emergence of next-generation sequencing. In this study, we introduce state-of-the-art research on AS and then compare the research methods and software tools available for AS based on next-generation sequencing reads. Finally, we discuss the prospects of computational methods related to AS.

[Identification of variants in LMF1 gene associated with primary hypertriglyceridemia]

The majority of severe primary hypertriglyceridemia (HTG) are diagnosed in adults, and their molecular bases have not yet been fully defined. The promoter, coding regions and intron-exon boundaries of LMF1 were sequenced in 112 patients with severe primary hipertrigliceridemia (defined as TG above 500mg/dl). Five patients (4.46%) were carriers of four rare variants in the LMF1 gene associated with HTG, which participate in lipoprotein lipase (LpL) function. Also, we have identified two common variants, c.194-28 T>G and c.729+18C>G that were associated with HTG, with a different allelic frequency to that observed in the general population. A bioinformatic analysis of all found variants was conducted, defining the following as potentially harmful: p.Arg364Gln, p.Arg451Trp, p.Pro562Arg and p.Leu85Leu. Our results suggest that LMF1 mutations are involved in a substantial proportion of cases with severe HTG, putting together the moderate-aggressive effect of rare mutations with polymorphisms classically associated with this disease.

A New Homozygous IGF1R Variant Defines a Clinically Recognizable Incomplete Dominant form of SHORT Syndrome

Here, we describe a child, born from consanguineous parents, with clinical features of SHORT syndrome, high IGF1 levels, developmental delay, CNS defects, and marked progeroid appearance. By exome sequencing, we identified a new homozygous c.2201G>T missense mutation in the IGF1R gene. Proband's parents and other relatives, all heterozygous carriers of the mutation, presented with milder phenotype including high IGFI levels, short stature, and type 2 diabetes. Functional studies using patient's cell lines showed a lower IGF1R expression that leads to the alteration of IGF1R-mediated PI3K/AKT/mTOR downstream pathways, including autophagy. This study defines a clinically recognizable incomplete dominant form of SHORT syndrome, and provides relevant insights into the pathophysiological and phenotypical consequences of IGF1R mutations.

Architecture and biological applications of artificial neural networks: a tuberculosis perspective

Advancement of science and technology has prompted researchers to develop new intelligent systems that can solve a variety of problems such as pattern recognition, prediction, and optimization. The ability of the human brain to learn in a fashion that tolerates noise and error has attracted many researchers and provided the starting point for the development of artificial neural networks: the intelligent systems. Intelligent systems can acclimatize to the environment or data and can maximize the chances of success or improve the efficiency of a search. Due to massive parallelism with large numbers of interconnected processers and their ability to learn from the data, neural networks can solve a variety of challenging computational problems. Neural networks have the ability to derive meaning from complicated and imprecise data; they are used in detecting patterns, and trends that are too complex for humans, or other computer systems. Solutions to the toughest problems will not be found through one narrow specialization; therefore we need to combine interdisciplinary approaches to discover the solutions to a variety of problems. Many researchers in different disciplines such as medicine, bioinformatics, molecular biology, and pharmacology have successfully applied artificial neural networks. This chapter helps the reader in understanding the basics of artificial neural networks, their applications, and methodology; it also outlines the network learning process and architecture. We present a brief outline of the application of neural networks to medical diagnosis, drug discovery, gene identification, and protein structure prediction. We conclude with a summary of the results from our study on tuberculosis data using neural networks, in diagnosing active tuberculosis, and predicting chronic vs. infiltrative forms of tuberculosis.

CIC inactivating mutations identify aggressive subset of 1p19q codeleted gliomas

OBJECTIVE

CIC gene is frequently mutated in oligodendroglial tumors with 1p19q codeletion. However, clinical and biological impact remain poorly understood.

METHODS

We sequenced the CIC gene on 127 oligodendroglial tumors (109 with the 1p19q codeletion) and analyzed patients' outcome. We compared magnetic resonance imaging, transcriptomic profile, and CIC protein expression of CIC wild-type (WT) and mutant gliomas. We compared the level of expression of CIC target genes on Hs683-IDH1(R132H) cells transfected with lentivirus encoding mutant and WT CIC.

RESULTS

We found 63 mutations affecting 60 of 127 patients, virtually all 1p19q codeleted and IDH mutated (59 of 60). In the 1p19q codeleted gliomas, CIC mutations were associated with a poorer outcome by uni- (p = 0.001) and multivariate analysis (p < 0.016). CIC mutation prognostic impact was validated on the TCGA cohort. CIC mutant grade II codeleted gliomas spontaneously grew faster than WTs. Transcriptomic analysis revealed an enrichment of proliferative pathways and oligodendrocyte precursor cell gene expression profile in CIC mutant gliomas, with upregulation of normally CIC repressed genes ETV1, ETV4, ETV5, and CCND1. Various missense mutations resulted in CIC protein expression loss. Moreover, a truncating CIC mutation resulted in a defect of nuclear targeting of CIC protein to the nucleus in a human glioma cell line expressing IDH1(R132H) and overexpression of CCND1 and other new target genes of CIC, such as DUSP4 and SPRED1.

INTERPRETATION

CIC mutations result in protein inactivation with upregulation of CIC target genes, activation of proliferative pathways, inhibition of differentiation, and poorer outcome in patients with a 1p19q codeleted glioma.

Two novel UBR1 gene mutations ın a patient with Johanson Blizzard Syndrome: A mild phenotype without mental retardation

Johanson-Blizzard Syndrome (JBS) (MIM #243800) is a rare autosomal recessive genetic disorder characterized by exocrine pancreatic insufficiency, abnormal facial appearance and varying degrees of mental retardation. Mutations in UBR1 gene (MIM *605981) are considered to be responsible for the syndrome. Here, we report a 3 year-old mentally normal JBS girl. The patient presented with exocrine pancreatic insufficiency as well as failure-to-thrive. On dysmorphological examination, she was noted to have an abnormal hair pattern with frontal upsweep and alae nasi hypoplasia. With these findings, JBS diagnosis was established clinically. Molecular analysis of the UBR1 gene revealed two inherited novel mutations; one coming from each parent. These novel mutations were c. 1280T>G and c. 2432+5G>C, and they were found to be disease causing via in-silico analysis. In conclusion, for patients with longstanding exocrine pancreatic insufficiency, it should be considered as being symptomatic of a far broader picture. To omit connection with rare genetic diseases, such as Johanson-Blizzard Syndrome, a detailed dysmorphological examination ought to be performed.

Mutations in a TGF-β ligand, TGFB3, cause syndromic aortic aneurysms and dissections

Background

Aneurysms affecting the aorta are a common condition associated with high mortality as a result of aortic dissection or rupture. Investigations of the pathogenic mechanisms involved in syndromic types of thoracic aortic aneurysms, such as Marfan and Loeys-Dietz syndromes, have revealed an important contribution of disturbed transforming growth factor (TGF)-β signaling.

Objectives

This study sought to discover a novel gene causing syndromic aortic aneurysms in order to unravel the underlying pathogenesis.

Methods

We combined genome-wide linkage analysis, exome sequencing, and candidate gene Sanger sequencing in a total of 470 index cases with thoracic aortic aneurysms. Extensive cardiological examination, including physical examination, electrocardiography, and transthoracic echocardiography was performed. In adults, imaging of the entire aorta using computed tomography or magnetic resonance imaging was done.

Results

Here, we report on 43 patients from 11 families with syndromic presentations of aortic aneurysms caused by TGFB3 mutations. We demonstrate that TGFB3 mutations are associated with significant cardiovascular involvement, including thoracic/abdominal aortic aneurysm and dissection, and mitral valve disease. Other systemic features overlap clinically with Loeys-Dietz, Shprintzen-Goldberg, and Marfan syndromes, including cleft palate, bifid uvula, skeletal overgrowth, cervical spine instability and clubfoot deformity. In line with previous observations in aortic wall tissues of patients with mutations in effectors of TGF-β signaling (TGFBR1/2, SMAD3, and TGFB2), we confirm a paradoxical up-regulation of both canonical and noncanonical TGF-β signaling in association with up-regulation of the expression of TGF-β ligands.

Conclusions

Our findings emphasize the broad clinical variability associated with TGFB3 mutations and highlight the importance of early recognition of the disease because of high cardiovascular risk.

RNA-Binding Proteins: Splicing Factors and Disease

Pre-mRNA splicing is mediated by interactions of the Core Spliceosome and an array of accessory RNA binding proteins with cis-sequence elements. Splicing is a major regulatory component in higher eukaryotes. Disruptions in splicing are a major contributor to human disease. One in three hereditary disease alleles are believed to cause aberrant splicing. Hereditary disease alleles can alter splicing by disrupting a splicing element, creating a toxic RNA, or affecting splicing factors. One of the challenges of medical genetics is identifying causal variants from the thousands of possibilities discovered in a clinical sequencing experiment. Here we review the basic biochemistry of splicing, the mechanisms of splicing mutations, the methods for identifying splicing mutants, and the potential of therapeutic interventions.

Mitochondrial DNA Polymerase γ Mutations and Their Implications in mtDNA Alterations in Colorectal Cancer

Mitochondrial DNA was found to be highly mutated in colorectal cancer cells. One of the key molecules involved in the maintenance of the mitochondrial genome is the nuclear-encoded polymerase gamma. The aim of our study was to determine if there is a link between polymorphisms within the polymerase gamma gene (POLG) and somatic mutations within the mitochondrial genome in cancer cells. We investigated POLG sequence variability in 50 colorectal cancer patients whose complete mitochondrial genome sequences were determined. Relative mtDNA copy number was also determined. We identified 251 sequence variants in the POLG gene. Most of them were germline-specific (∼92%). Twenty-one somatic changes in POLG were found in 10 colorectal cancer patients. We have found no association between the occurrence of mtDNA somatic mutations and the somatically occurring variants in POLG. MtDNA content was reduced in patients carrying somatic variants in POLG or germline nucleotide variants located in the region encoding the POLG polymerase domain, but the difference did not reach statistical significance. Our findings suggest that somatic mtDNA mutations occurring in colorectal cancer are not a consequence of somatic mutations in POLG. Nevertheless, POLG nucleotide variants may lead to a decrease in mtDNA content, and consequently result in mitochondrial dysfunction.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Analysis of HGD Gene Mutations in Patients with Alkaptonuria from the United Kingdom: Identification of Novel Mutations

Alkaptonuria (AKU) is a rare autosomal recessive disorder with incidence ranging from 1:100,000 to 1:250,000. The disorder is caused by a deficiency of the enzyme homogentisate 1,2-dioxygenase (HGD), which results from defects in the HGD gene. This enzyme converts homogentisic acid to maleylacetoacetate and has a major role in the catabolism of phenylalanine and tyrosine. To elucidate the mutation spectrum of the HGD gene in patients with alkaptonuria from 42 patients attending the National Alkaptonuria Centre, 14 exons of the HGD gene and the intron-exon boundaries were analysed by PCR-based sequencing. A total of 34 sequence variants was observed, confirming the genetic heterogeneity of AKU. Of these mutations, 26 were missense substitutions and four splice site mutations. There were two deletions and one duplication giving rise to frame shifts and one substitution abolishing the translation termination codon (no stop). Nine of the mutations were previously unreported novel variants. Using computational approaches based on the 3D structure, these novel mutations are predicted to affect the activity of the protein complex through destabilisation of the individual protomer structure or through disruption of protomer-protomer interactions.

Detection of novel polyomaviruses in fruit bats in Indonesia

Bats are an important natural reservoir for a variety of viral pathogens, including polyomaviruses (PyVs). The aims of this study were: (i) to determine which PyVs are present in bats in Indonesia and (ii) to analyze the evolutionary relationships between bat PyVs and other known PyVs. Using broad-spectrum polymerase chain reaction (PCR)-based assays, we screened PyV DNA isolated from spleen samples from 82 wild fruit bats captured in Indonesia. Fragments of the PyV genome were detected in 10 of the 82 spleen samples screened, and eight full-length viral genome sequences were obtained using an inverse PCR method. A phylogenetic analysis of eight whole viral genome sequences showed that BatPyVs form two distinct genetic clusters within the proposed genus Orthopolyomavirus that are genetically different from previously described BatPyVs. Interestingly, one group of BatPyVs is genetically related to the primate PyVs, including human PyV9 and trichodysplasia spinulosa-associated PyV. This study has identified the presence of novel PyVs in fruit bats in Indonesia and provides genetic information about these BatPyVs.

OSBPL2 encodes a protein of inner and outer hair cell stereocilia and is mutated in autosomal dominant hearing loss (DFNA67)

Background

Early-onset hearing loss is mostly of genetic origin. The complexity of the hearing process is reflected by its extensive genetic heterogeneity, with probably many causative genes remaining to be identified. Here, we aimed at identifying the genetic basis for autosomal dominant non-syndromic hearing loss (ADNSHL) in a large German family.

Methods

A panel of 66 known deafness genes was analyzed for mutations by next-generation sequencing (NGS) in the index patient. We then conducted genome-wide linkage analysis, and whole-exome sequencing was carried out with samples of two patients. Expression of Osbpl2 in the mouse cochlea was determined by immunohistochemistry. Because Osbpl2 has been proposed as a target of miR-96, we investigated homozygous Mir96 mutant mice for its upregulation.

Results

Onset of hearing loss in the investigated ADNSHL family is in childhood, initially affecting the high frequencies and progressing to profound deafness in adulthood. However, there is considerable intrafamilial variability. We mapped a novel ADNSHL locus, DFNA67, to chromosome 20q13.2-q13.33, and subsequently identified a co-segregating heterozygous frameshift mutation, c.141_142delTG (p.Arg50Alafs*103), in OSBPL2, encoding a protein known to interact with the DFNA1 protein, DIAPH1. In mice, Osbpl2 was prominently expressed in stereocilia of cochlear outer and inner hair cells. We found no significant Osbpl2 upregulation at the mRNA level in homozygous Mir96 mutant mice.

Conclusion

The function of OSBPL2 in the hearing process remains to be determined. Our study and the recent description of another frameshift mutation in a Chinese ADNSHL family identify OSBPL2 as a novel gene for progressive deafness.

Electronic supplementary material

The online version of this article (doi:10.1186/s13023-015-0238-5) contains supplementary material, which is available to authorized users.

Genetic variation and RNA binding proteins: tools and techniques to detect functional polymorphisms

At its most fundamental level the goal of genetics is to connect genotype to phenotype. This question is asked at a basic level evaluating the role of genes and pathways in genetic model organism. Increasingly, this question is being asked in the clinic. Genomes of individuals and populations are being sequenced and compared. The challenge often comes at the stage of analysis. The variant positions are analyzed with the hope of understanding human disease. However after a genome or exome has been sequenced, the researcher is often deluged with hundreds of potentially relevant variations. Traditionally, amino-acid changing mutations were considered the tractable class of disease-causing mutations; however, mutations that disrupt noncoding elements are the subject of growing interest. These noncoding changes are a major avenue of disease (e.g., one in three hereditary disease alleles are predicted to affect splicing). Here, we review some current practices of medical genetics, the basic theory behind biochemical binding and functional assays, and then explore technical advances in how variations that alter RNA protein recognition events are detected and studied. These advances are advances in scale-high-throughput implementations of traditional biochemical assays that are feasible to perform in any molecular biology laboratory. This chapter utilizes a case study approach to illustrate some methods for analyzing polymorphisms. The first characterizes a functional intronic SNP that deletes a high affinity PTB site using traditional low-throughput biochemical and functional assays. From here we demonstrate the utility of high-throughput splicing and spliceosome assembly assays for screening large sets of SNPs and disease alleles for allelic differences in gene expression. Finally we perform three pilot drug screens with small molecules (G418, tetracycline, and valproic acid) that illustrate how compounds that rescue specific instances of differential pre-mRNA processing can be discovered.

In silico prediction of splice-altering single nucleotide variants in the human genome

In silico tools have been developed to predict variants that may have an impact on pre-mRNA splicing. The major limitation of the application of these tools to basic research and clinical practice is the difficulty in interpreting the output. Most tools only predict potential splice sites given a DNA sequence without measuring splicing signal changes caused by a variant. Another limitation is the lack of large-scale evaluation studies of these tools. We compared eight in silico tools on 2959 single nucleotide variants within splicing consensus regions (scSNVs) using receiver operating characteristic analysis. The Position Weight Matrix model and MaxEntScan outperformed other methods. Two ensemble learning methods, adaptive boosting and random forests, were used to construct models that take advantage of individual methods. Both models further improved prediction, with outputs of directly interpretable prediction scores. We applied our ensemble scores to scSNVs from the Catalogue of Somatic Mutations in Cancer database. Analysis showed that predicted splice-altering scSNVs are enriched in recurrent scSNVs and known cancer genes. We pre-computed our ensemble scores for all potential scSNVs across the human genome, providing a whole genome level resource for identifying splice-altering scSNVs discovered from large-scale sequencing studies.

SNPlice: variants that modulate Intron retention from RNA-sequencing data

RATIONALE

The growing recognition of the importance of splicing, together with rapidly accumulating RNA-sequencing data, demand robust high-throughput approaches, which efficiently analyze experimentally derived whole-transcriptome splice profiles.

RESULTS

We have developed a computational approach, called SNPlice, for identifying cis-acting, splice-modulating variants from RNA-seq datasets. SNPlice mines RNA-seq datasets to find reads that span single-nucleotide variant (SNV) loci and nearby splice junctions, assessing the co-occurrence of variants and molecules that remain unspliced at nearby exon-intron boundaries. Hence, SNPlice highlights variants preferentially occurring on intron-containing molecules, possibly resulting from altered splicing. To illustrate co-occurrence of variant nucleotide and exon-intron boundary, allele-specific sequencing was used. SNPlice results are generally consistent with splice-prediction tools, but also indicate splice-modulating elements missed by other algorithms. SNPlice can be applied to identify variants that correlate with unexpected splicing events, and to measure the splice-modulating potential of canonical splice-site SNVs.

AVAILABILITY AND IMPLEMENTATION

SNPlice is freely available for download from https://code.google.com/p/snplice/ as a self-contained binary package for 64-bit Linux computers and as python source-code.

CONTACT

pmudvari@gwu.edu or horvatha@gwu.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of six novel mutations in ZEB1 and description of the associated phenotypes in patients with posterior polymorphous corneal dystrophy 3

Posterior polymorphous corneal dystrophy 3 (PPCD3) is a rare autosomal dominant disorder caused by mutations in ZEB1. To date all identified disease-causing variants were unique to the studied families, except for c.1576dup. We have detected six novel ZEB1 mutations; c.1749_1750del; p.(Pro584*) and c.1717_1718del; p.(Val573Phefs*12) in two Czech families, c.1176dup; p.(Ala393Serfs*19), c.1100C>A; p.(Ser367*), c.627del; p.(Phe209Leufs*11) in three British families and a splice site mutation, c.685-2A>G, in a patient of Sri Lankan origin. An additional British proband had the c.1576dup; p.(Val526Glyfs*3) mutation previously reported in other populations. Clinical findings were variable and included bilateral congenital corneal opacity in one proband, development of opacity before the age of 2 years in another individual and bilateral iris flocculi in yet another subject. The majority of eyes examined by corneal topography (10 out of 16) had an abnormally steep cornea (flat keratometry 46.5-52.7 diopters, steep keratometry 48.1-54.0 diopters). One proband underwent surgery for cryptorchidism. Our study further demonstrates that PPCD3 can present as corneal edema in early childhood, and that an abnormally steep keratometry is a common feature of this condition. As cryptorchidism has been previously observed in two other PPCD3 cases, its association with the disease warrants further investigation.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Genome-wide identification and analysis of basic helix-loop-helix domains in dog, Canis lupus familiaris

The basic helix-loop-helix (bHLH) domain is a highly conserved amino acid motif that defines a group of DNA-binding transcription factors. bHLH proteins play essential regulatory roles in a variety of biological processes in animal, plant, and fungus. The domestic dog, Canis lupus familiaris, is a good model organism for genetic, physiological, and behavioral studies. In this study, we identified 115 putative bHLH genes in the dog genome. Based on a phylogenetic analysis, 51, 26, 14, 4, 12, and 4 dog bHLH genes were assigned to six separate groups (A-F); four bHLH genes were categorized as ''orphans''. Within-group evolutionary relationships inferred from the phylogenetic analysis were consistent with positional conservation, other conserved domains flanking the bHLH motif, and highly conserved intron/exon patterns in other vertebrates. Our analytical results confirmed the GenBank annotations of 89 dog bHLH proteins and provided information that could be used to update the annotations of the remaining 26 dog bHLH proteins. These data will provide good references for further studies on the structures and regulatory functions of bHLH proteins in the growth and development of dogs, which may help in understanding the mechanisms that underlie the physical and behavioral differences between dogs and wolves.

The ORF012 gene of Marek's disease virus type 1 produces a spliced transcript and encodes a novel nuclear phosphoprotein essential for virus growth

Marek's disease virus (MDV), an alphaherpesvirus, is the causative agent of a lethal disease in chickens characterized by generalized nerve inflammation and rapid lymphoma development. The extensive colinearity of the MDV genome with those of related herpesviruses has eased functional characterization of many MDV genes. However, MDV carries a number of unique open reading frames (ORFs) that have not yet been investigated regarding their coding potentials and the functions of their products. Among these unique ORFs are two putative ORFs, ORF011 and ORF012, which are found at the extreme left end of the MDV unique long region. Using reverse transcriptase PCR, we showed that ORF011 and ORF012 are not individual genes but form a single gene through mRNA splicing of a small intron, resulting in the novel ORF012. We generated an ORF012-null virus using an infectious clone of MDV strain RB-1B. The deletion virus had a marked growth defect in vitro and could not be passaged in cultured cells, suggesting an essential role for the ORF012 product in virus replication. Further studies revealed that protein 012 (p012) localized to the nucleus in transfected and infected cells, and we identified by site-directed mutagenesis and green fluorescent protein (GFP) reporter fusion assays a nuclear localization signal (NLS) that was mapped to a 23-amino-acid sequence at the protein's C terminus. Nuclear export was blocked using leptomycin B, suggesting a potential role for p012 as a nuclear/cytoplasmic shuttling protein. Finally, p012 is phosphorylated at multiple residues, a modification that could possibly regulate its subcellular distribution.

IMPORTANCE

Marek's disease virus (MDV) causes a devastating oncogenic disease in chickens with high morbidity and mortality. The costs for disease prevention reach several billion dollars annually. The functional investigation of MDV genes is necessary to understand its complex replication cycle, which eventually could help us to interfere with MDV and herpesviral pathogenesis. We have identified a previously unidentified phosphoprotein encoded by MDV ORF012. We were able to show experimentally that predicted splicing of the gene based on bioinformatics data does indeed occur during replication. The newly identified p012 is essential for MDV replication and localizes to the nucleus due to the presence of a transferable nuclear localization signal at its C terminus. Our results also imply that p012 could constitute a nucleocytoplasmic shuttle protein, a feature that could prove interesting and important.

Rapid identification of a natural knockout allele of ARMADILLO REPEAT-CONTAINING KINESIN1 that causes root hair branching by mapping-by-sequencing

In Arabidopsis (Arabidopsis thaliana), branched root hairs are an indicator of defects in root hair tip growth. Among 62 accessions, one accession (Heiligkreuztal2 [HKT2.4]) displayed branched root hairs, suggesting that this accession carries a mutation in a gene of importance for tip growth. We determined 200- to 300-kb mapping intervals using a mapping-by-sequencing approach of F2 pools from crossings of HKT2.4 with three different accessions. The intersection of these mapping intervals was 80 kb in size featuring not more than 36 HKT2.4-specific single nucleotide polymorphisms, only two of which changed the coding potential of genes. Among them, we identified the causative single nucleotide polymorphism changing a splicing site in ARMADILLO REPEAT-CONTAINING KINESIN1. The applied strategies have the potential to complement statistical methods in high-throughput phenotyping studies using different natural accessions to identify causative genes for distinct phenotypes represented by only one or a few accessions.

Neuropeptides encoded by the genomes of the Akoya pearl oyster Pinctata fucata and Pacific oyster Crassostrea gigas: a bioinformatic and peptidomic survey

BACKGROUND

Oysters impart significant socio-ecological benefits from primary production of food supply, to estuarine ecosystems via reduction of water column nutrients, plankton and seston biomass. Little though is known at the molecular level of what genes are responsible for how oysters reproduce, filter nutrients, survive stressful physiological events and form reef communities. Neuropeptides represent a diverse class of chemical messengers, instrumental in orchestrating these complex physiological events in other species.

RESULTS

By a combination of in silico data mining and peptide analysis of ganglia, 74 putative neuropeptide genes were identified from genome and transcriptome databases of the Akoya pearl oyster, Pinctata fucata and the Pacific oyster, Crassostrea gigas, encoding precursors for over 300 predicted bioactive peptide products, including three newly identified neuropeptide precursors PFGx8amide, RxIamide and Wx3Yamide. Our findings also include a gene for the gonadotropin-releasing hormone (GnRH) and two egg-laying hormones (ELH) which were identified from both oysters. Multiple sequence alignments and phylogenetic analysis supports similar global organization of these mature peptides. Computer-based peptide modeling of the molecular tertiary structures of ELH highlights the structural homologies within ELH family, which may facilitate ELH activity leading to the release of gametes.

CONCLUSION

Our analysis demonstrates that oysters possess conserved molluscan neuropeptide domains and overall precursor organization whilst highlighting many previously unrecognized bivalve idiosyncrasies. This genomic analysis provides a solid foundation from which further studies aimed at the functional characterization of these molluscan neuropeptides can be conducted to further stimulate advances in understanding the ecology and cultivation of oysters.

GALNS mutations in Indian patients with mucopolysaccharidosis IVA

Mucopolysaccharidosis IV A (Morquio syndrome A, MPS IVA) is a lysosomal storage disease caused by the deficiency of N-acetylgalactosamine-6-sulfatase (GALNS). The mutation spectrum in this condition is yet to be determined in Indians. We aimed to analyze the mutations in the GALNS gene in Asian Indians with MPS IVA. All the exons and the adjacent intronic regions of the gene were amplified and sequenced in sixty-eight unrelated Indian families. We identified 136 mutant alleles comprising of 40 different mutations. We report twenty-two novel mutations that comprise of seventeen missense (p.Asn32Thr, p.Leu36Arg, p.Pro52Leu, p.Pro77Ser, p.Cys79Arg, p.His142Pro, p.Tyr191Asp, p.Asn204Thr, p.Gly188Ser, p.Phe216Ser, p.Trp230Cys, p.Ala291Ser, p.Gly317Arg, p.His329Pro, p.Arg386Ser, p.Glu450Gly, p.Cys501Ser), three splice-site variants (c.120+1G>C, c.1003-3C>G, c.1139+1G>A), one nonsense mutation (p.Gln414*) and one frameshift mutation (p.Pro420Leufs*440). Eighteen mutations have been reported earlier. Among these p.Ser287Leu (8.82%), p.Phe216Ser (7.35%), p.Asn32Thr (6.61%) and p.Ala291Ser (5.88%) were the most frequent mutations in Indian patients but were rare in the mutational profiles reported in other populations. These results indicate that the Indian patients may have a distinct mutation spectrum compared to those of other populations. Mutant alleles in exon 1, 7 and 8 accounted for 44.8% of the mutations, and sequencing of these exons initially may be a cost-effective approach in Asian Indian patients. This is the largest study on molecular analysis of patients with MPS IVA reported in the literature, and the first report from India.

Gammaretroviral pol sequences act in cis to direct polysome loading and NXF1/NXT-dependent protein production by gag-encoded RNA

Background

All retroviruses synthesize essential proteins via alternatively spliced mRNAs. Retrovirus genera, though, exploit different mechanisms to coordinate the synthesis of proteins from alternatively spliced mRNAs. The best studied of these retroviral, post-transcriptional effectors are the trans-acting Rev protein of lentiviruses and the cis-acting constitutive transport element (CTE) of the betaretrovirus Mason-Pfizer monkey virus (MPMV). How members of the gammaretrovirus genus translate protein from unspliced RNA has not been elucidated.

Results

The mechanism by which two gammaretroviruses, XMRV and MLV, synthesize the Gag polyprotein (Pr65^Gag) from full-length, unspliced mRNA was investigated here. The yield of Pr65^Gag from a gag–only expression plasmid was found to be at least 30-fold less than that from an otherwise isogenic gag-pol expression plasmid. A frameshift mutation disrupting the pol open reading frame within the gag-pol expression plasmid did not decrease Pr65^Gag production and 398 silent nucleotide changes engineered into gag rendered Pr65^Gag synthesis pol-independent. These results are consistent with pol-encoded RNA acting in cis to promote Pr65^Gag translation. Two independently-acting pol fragments were identified by screening 17 pol deletion mutations. To determine the mechanism by which pol promoted Pr65^Gag synthesis, gag RNA in total and cytoplasmic fractions was quantitated by northern blot and by RT-PCR. The pol sequences caused, maximally, three-fold increase in total or cytoplasmic gag mRNA. Instead, pol sequences increased gag mRNA association with polyribosomes ~100-fold, a magnitude sufficient to explain the increase in Pr65^Gag translation efficiency. The MPMV CTE, an NXF1-binding element, substituted for pol in promoting Pr65^Gag synthesis. A pol RNA stem-loop resembling the CTE promoted Pr65^Gag synthesis. Over-expression of NXF1 and NXT, host factors that bind to the MPMV CTE, synergized with pol to promote gammaretroviral gag RNA loading onto polysomes and to increase Pr65^Gag synthesis. Conversely, Gag polyprotein synthesis was decreased by NXF1 knockdown. Finally, overexpression of SRp20, a shuttling protein that binds to NXF1 and promotes NXF1 binding to RNA, also increased gag RNA loading onto polysomes and increased Pr65^Gag synthesis.

Conclusion

These experiments demonstrate that gammaretroviral pol sequences act in cis to recruit NXF1 and SRp20 to promote polysome loading of gag RNA and, thereby license the synthesis of Pr65^Gag from unspliced mRNA.

NAMPT and NAPRT1: novel polymorphisms and distribution of variants between normal tissues and tumor samples

Nicotinamide phosphoribosyltransferase (NAMPT) and nicotinate phosphoribosyltransferase domain containing 1 (NAPRT1) are the main human NAD salvage enzymes. NAD regulates energy metabolism and cell signaling, and the enzymes that control NAD availability are linked to pathologies such as cancer and neurodegeneration. Here, we have screened normal and tumor samples from different tissues and populations of origin for mutations in human NAMPT and NAPRT1, and evaluated their potential pathogenicity. We have identified several novel polymorphisms and showed that NAPRT1 has a greater genetic diversity than NAMPT, where any alteration can have a greater functional impact. Some variants presented different frequencies between normal and tumor samples that were most likely related to their population of origin. The novel mutations described that affect protein structure or expression levels can be functionally relevant and should be considered in a disease context. Particularly, mutations that decrease NAPRT1 expression can predict the usefulness of Nicotinic Acid in tumor treatments with NAMPT inhibitors.

Comprehensive screening for mutations associated with colorectal cancer in unselected cases reveals penetrant and nonpenetrant mutations

Germline mutation testing in patients with colorectal cancer (CRC) is offered only to a subset of patients with a clinical presentation or tumor histology suggestive of familial CRC syndromes, probably underestimating familial CRC predisposition. The aim of our study was to determine whether unbiased screening of newly diagnosed CRC cases with next generation sequencing (NGS) increases the overall detection rate of germline mutations. We analyzed 152 consecutive CRC patients for germline mutations in 18 CRC-associated genes using NGS. All patients were also evaluated for Bethesda criteria and all tumors were investigated for microsatellite instability, immunohistochemistry for mismatch repair proteins and the BRAF*V600E somatic mutation. NGS based sequencing identified 27 variants in 9 genes in 23 out of 152 patients studied (18%). Three of them were already reported as pathogenic and 12 were class 3 germline variants with an uncertain prediction of pathogenicity. Only 1 of these patients fulfilled Bethesda criteria and had a microsatellite instable tumor and an MLH1 germline mutation. The others would have been missed with current approaches: 2 with a MSH6 premature termination mutation and 12 uncertain, potentially pathogenic class 3 variants in APC, MLH1, MSH2, MSH6, MSH3 and MLH3. The higher NGS mutation detection rate compared with current testing strategies based on clinicopathological criteria is probably due to the large genetic heterogeneity and overlapping clinical presentation of the various CRC syndromes. It can also identify apparently nonpenetrant germline mutations complicating the clinical management of the patients and their families.

Analysis of POFUT1 gene mutation in a Chinese family with Dowling-Degos disease

Dowling-Degos disease (DDD) is an autosomal dominant genodermatosis characterized by reticular pigmented anomaly mainly affecting flexures. Though KRT5 has been identified to be the causal gene of DDD, the heterogeneity of this disease was displayed: for example, POFUT1 and POGLUT1 were recently identified and confirmed to be additional pathogenic genes of DDD. To identify other DDD causative genes, we performed genome-wide linkage and exome sequencing analyses in a multiplex Chinese DDD family, in which the KRT5 mutation was absent. Only a novel 1-bp deletion (c.246+5delG) in POFUT1 was found. No other novel mutation or this deletion was detected in POFUT1 in a second DDD family and a sporadic DDD case by Sanger Sequencing. The result shows the genetic-heterogeneity and complexity of DDD and will contribute to the further understanding of DDD genotype/phenotype correlations and to the pathogenesis of this disease.

Targeted next-generation sequencing reveals novel USH2A mutations associated with diverse disease phenotypes: implications for clinical and molecular diagnosis

USH2A mutations have been implicated in the disease etiology of several inherited diseases, including Usher syndrome type 2 (USH2), nonsyndromic retinitis pigmentosa (RP), and nonsyndromic deafness. The complex genetic and phenotypic spectrums relevant to USH2A defects make it difficult to manage patients with such mutations. In the present study, we aim to determine the genetic etiology and to characterize the correlated clinical phenotypes for three Chinese pedigrees with nonsyndromic RP, one with RP sine pigmento (RPSP), and one with USH2. Family histories and clinical details for all included patients were reviewed. Ophthalmic examinations included best corrected visual acuities, visual field measurements, funduscopy, and electroretinography. Targeted next-generation sequencing (NGS) was applied using two sequence capture arrays to reveal the disease causative mutations for each family. Genotype-phenotype correlations were also annotated. Seven USH2A mutations, including four missense substitutions (p.P2762A, p.G3320C, p.R3719H, and p.G4763R), two splice site variants (c.8223+1G>A and c.8559-2T>C), and a nonsense mutation (p.Y3745*), were identified as disease causative in the five investigated families, of which three reported to have consanguineous marriage. Among all seven mutations, six were novel, and one was recurrent. Two homozygous missense mutations (p.P2762A and p.G3320C) were found in one individual family suggesting a potential double hit effect. Significant phenotypic divergences were revealed among the five families. Three families of the five families were affected with early, moderated, or late onset RP, one with RPSP, and the other one with USH2. Our study expands the genotypic and phenotypic variability relevant to USH2A mutations, which would help with a clear insight into the complex genetic and phenotypic spectrums relevant to USH2A defects, and is complementary for a better management of patients with such mutations. We have also demonstrated that a targeted NGS approach is a valuable tool for the genetic diagnosis of USH2 and RP.