At a glance: the largest Niemann-Pick type C1 cohort with 602 patients diagnosed over 15 years.

Niemann-Pick type C1 disease (NPC1 [OMIM 257220]) is a rare and severe autosomal recessive disorder, characterized by a multitude of neurovisceral clinical manifestations and a fatal outcome with no effective treatment to date. Aiming to gain insights into the genetic aspects of the disease, clinical, genetic, and biomarker PPCS data from 602 patients referred from 47 countries and diagnosed with NPC1 in our laboratory were analyzed. Patients' clinical data were dissected using Human Phenotype Ontology (HPO) terms, and genotype-phenotype analysis was performed. The median age at diagnosis was 10.6 years (range 0-64.5 years), with 287 unique pathogenic/likely pathogenic (P/LP) variants identified, expanding NPC1 allelic heterogeneity. Importantly, 73 P/LP variants were previously unpublished. The most frequent variants detected were: c.3019C > G, p.(P1007A), c.3104C > T, p.(A1035V), and c.2861C > T, p.(S954L). Loss of function (LoF) variants were significantly associated with earlier age at diagnosis, highly increased biomarker levels, and a visceral phenotype (abnormal abdomen and liver morphology). On the other hand, the variants p.(P1007A) and p.(S954L) were significantly associated with later age at diagnosis (p < 0.001) and mildly elevated biomarker levels (p ≤ 0.002), consistent with the juvenile/adult form of NPC1. In addition, p.(I1061T), p.(S954L), and p.(A1035V) were associated with abnormality of eye movements (vertical supranuclear gaze palsy, p ≤ 0.05). We describe the largest and most heterogenous cohort of NPC1 patients published to date. Our results suggest that besides its utility in variant classification, the biomarker PPCS might serve to indicate disease severity/progression. In addition, we establish new genotype-phenotype relationships for "frequent" NPC1 variants.

Clinical and genetic characterization of a cohort of 97 CLN6 patients tested at a single center.

BACKGROUND: Ceroid lipofuscinoses neuronal 6 (CLN6) disease belongs to the neuronal ceroid lipofuscinoses (NCLs), complex and genetically heterogeneous disorders with wide geographical and phenotypic variation. The first clinical signs usually appear between 18 months and 8 years, but examples of later-onset have also been reported. Common manifestations include ataxia, seizures, vision impairment, and developmental regression. Because these are shared by other neurological diseases, identification of CLN6 genetic variants is imperative for early diagnosis. RESULTS: We present one of the largest cohorts to date of genetically diagnosed CLN6 patients screened at a single center. In total 97 subjects, originating from 20 countries were screened between 2010 and 2020. They comprised 86 late-infantile, eight juvenile, and three adult-onset cases (two patients with Kufs disease type A, and one with teenage progressive myoclonic epilepsy). The male to female ratio was 1.06: 1.00. The age at referral was between six months and 33 years. The time from disease onset to referral ranged from less than 1 month to 8.3 years. The clinical phenotype consisted of a combination of symptoms, as reported before. We characterized a total of 45 distinct variants defining 45 distinct genotypes. Twenty-four were novel variants, some with distinct geographic associations. Remarkably, c.257A > G (p.H86R) was present in five out of 23 unrelated Egyptian individuals but in no patients from other countries. The most common genotype was homozygosity for the c.794_796del in-frame deletion. It was present in about one-third of CLN6 patients (28 unrelated cases, and 2 familial cases), all with late-infantile onset. Variants with a high likelihood of causing loss of CLN6 function were found in 21% of cases and made up 33% of all distinct variants. Forty-four percent of variants were classified as pathogenic or likely pathogenic. CONCLUSIONS: Our study significantly expands the number of published clinical cases and the mutational spectrum of disease-associated CLN6 variants, especially for the Middle Eastern and North African regions. We confirm previous observations regarding the most prevalent symptoms and recommend including CLN6 in the genetic diagnosis of patients presenting with early-onset abnormalities of the nervous system, musculoskeletal system, and eye.

Combining exome/genome sequencing with data repository analysis reveals novel gene-disease associations for a wide range of genetic disorders.

PURPOSE: Within this study, we aimed to discover novel gene-disease associations in patients with no genetic diagnosis after exome/genome sequencing (ES/GS). METHODS: We followed two approaches: (1) a patient-centered approach, which after routine diagnostic analysis systematically interrogates variants in genes not yet associated to human diseases; and (2) a gene variant centered approach. For the latter, we focused on de novo variants in patients that presented with neurodevelopmental delay (NDD) and/or intellectual disability (ID), which are the most common reasons for genetic testing referrals. Gene-disease association was assessed using our data repository that combines ES/GS data and Human Phenotype Ontology terms from over 33,000 patients. RESULTS: We propose six novel gene-disease associations based on 38 patients with variants in the BLOC1S1, IPO8, MMP15, PLK1, RAP1GDS1, and ZNF699 genes. Furthermore, our results support causality of 31 additional candidate genes that had little published evidence and no registered OMIM phenotype (56 patients). The phenotypes included syndromic/nonsyndromic NDD/ID, oral-facial-digital syndrome, cardiomyopathies, malformation syndrome, short stature, skeletal dysplasia, and ciliary dyskinesia. CONCLUSION: Our results demonstrate the value of data repositories which combine clinical and genetic data for discovering and confirming gene-disease associations. Genetic laboratories should be encouraged to pursue such analyses for the benefit of undiagnosed patients and their families.

Successful application of genome sequencing in a diagnostic setting: 1007 index cases from a clinically heterogeneous cohort.

Despite clear technical superiority of genome sequencing (GS) over other diagnostic methods such as exome sequencing (ES), few studies are available regarding the advantages of its clinical application. We analyzed 1007 consecutive index cases for whom GS was performed in a diagnostic setting over a 2-year period. We reported pathogenic and likely pathogenic (P/LP) variants that explain the patients' phenotype in 212 of the 1007 cases (21.1%). In 245 additional cases (24.3%), a variant of unknown significance (VUS) related to the phenotype was reported. We especially investigated patients which had had ES with no genetic diagnosis (n = 358). For this group, GS diagnostic yield was 14.5% (52 patients with P/LP out of 358). GS should be especially indicated for ES-negative cases since up to 29.6% of them  could benefit from GS testing (14.5% with P/LP, n = 52 and 15.1% with VUS, n = 54). Genetic diagnoses in most of the ES-negative/GS-positive cases were determined by technical superiority of GS, i.e., access to noncoding regions and more uniform coverage. Importantly, we reported 79 noncoding variants, of which, 41 variants were classified as P/LP. Interpretation of noncoding variants remains challenging, and in many cases, complementary methods based on direct enzyme assessment, biomarker testing and RNA analysis are needed for variant classification and diagnosis. We present the largest cohort of patients with GS performed in a clinical setting to date. The results of this study should direct the decision for GS as standard second-line, or even first-line stand-alone test.

EIF2AK2 Missense Variants Associated with Early Onset Generalized Dystonia.

OBJECTIVE: The study was undertaken to identify a monogenic cause of early onset, generalized dystonia. METHODS: Methods consisted of genome-wide linkage analysis, exome and Sanger sequencing, clinical neurological examination, brain magnetic resonance imaging, and protein expression studies in skin fibroblasts from patients. RESULTS: We identified a heterozygous variant, c.388G>A, p.Gly130Arg, in the eukaryotic translation initiation factor 2 alpha kinase 2 (EIF2AK2) gene, segregating with early onset isolated generalized dystonia in 5 patients of a Taiwanese family. EIF2AK2 sequencing in 191 unrelated patients with unexplained dystonia yielded 2 unrelated Caucasian patients with an identical heterozygous c.388G>A, p.Gly130Arg variant, occurring de novo in one case, another patient carrying a different heterozygous variant, c.413G>C, p.Gly138Ala, and one last patient, born from consanguineous parents, carrying a third, homozygous variant c.95A>C, p.Asn32Thr. These 3 missense variants are absent from gnomAD, and are located in functional domains of the encoded protein. In 3 patients, additional neurological manifestations were present, including intellectual disability and spasticity. EIF2AK2 encodes a kinase (protein kinase R [PKR]) that phosphorylates eukaryotic translation initiation factor 2 alpha (eIF2α), which orchestrates the cellular stress response. Our expression studies showed abnormally enhanced activation of the cellular stress response, monitored by PKR-mediated phosphorylation of eIF2α, in fibroblasts from patients with EIF2AK2 variants. Intriguingly, PKR can also be regulated by PRKRA (protein interferon-inducible double-stranded RNA-dependent protein kinase activator A), the product of another gene causing monogenic dystonia. INTERPRETATION: We identified EIF2AK2 variants implicated in early onset generalized dystonia, which can be dominantly or recessively inherited, or occur de novo. Our findings provide direct evidence for a key role of a dysfunctional eIF2α pathway in the pathogenesis of dystonia. ANN NEUROL 2021;89:485-497.

Rapid Large-Scale COVID-19 Testing During Shortages.

The Coronavirus disease 2019 (COVID-19) pandemic caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has resulted in economic and social lockdowns in most countries all over the globe. Early identification of infected individuals is regarded as one of the most important prerequisites for fighting the pandemic and for returning to a 'New Normal'. Large-scale testing is therefore crucial, but is facing several challenges including shortage of sample collection tools and of molecular biological reagents, and the need for safe electronic communication of medical reports. We present the successful establishment of a holistic SARS-CoV-2 testing platform that covers proband registration, sample collection and shipment, sample testing, and report issuing. The RT-PCR-based virus detection, being central to the platform, was extensively validated: sensitivity and specificity were defined as 96.8% and 100%, respectively; intra-run and inter-run precision were <3%. A novel type of sample swab and an in-house-developed RNA extraction system were shown to perform as good as commercially available products. The resulting flexibility guarantees independence from the current bottlenecks in SARS-CoV-2 testing. Based on our technology, we offered testing at local, national, and global levels. In the present study, we report the results from approx. 18,000 SARS-CoV-2 tests in almost 10,000 individuals from a low-frequency SARS-CoV-2 pandemic area in a homogenous geographical region in north-eastern Germany for a period of 10 weeks (21 March to 31 May 2020). Among the probands, five SARS-CoV-2 positive cases were identified. Comparative analysis of corresponding virus genomes revealed a diverse origin from three of the five currently recognized SARS-CoV-2 phylogenetic clades. Our study exemplifies how preventive SARS-CoV-2 testing can be set up in a rapid and flexible manner. The application of our test has enabled a safe maintenance/resume of critical local infrastructure, e.g., nursing homes where more than 5000 elderlies and caretakers got tested. The strategy outlined by the present study may serve as a blueprint for the implementation of large-scale preventive SARS-CoV-2 testing elsewhere.

Novel clinical and genetic insight into CXorf56-associated intellectual disability.

Intellectual disability (ID) is one of most frequent reasons for genetic consultation. The complex molecular anatomy of ID ranges from complete chromosomal imbalances to single nucleotide variant changes occurring de novo, with thousands of genes identified. This extreme genetic heterogeneity challenges the molecular diagnosis, which mostly requires a genomic approach. CXorf56 is largely uncharacterized and was recently proposed as a candidate ID gene based on findings in a single Dutch family. Here, we describe nine cases (six males and three females) from three unrelated families. Exome sequencing and combined database analyses, identified family-specific CXorf56 variants (NM_022101.3:c.498_503del, p.(Glu167_Glu168del) and c.303_304delCTinsACCC, p.(Phe101Leufs*20)) that segregated with the ID phenotype. These variants are presumably leading to loss-of-function, which is the proposed disease mechanism. Clinically, CXorf56-related disease is a slowly progressive neurological disorder. The phenotype is more severe in hemizygote males, but might also manifests in heterozygote females, which showed skewed X-inactivation patterns in blood. Male patients might present previously unreported neurological features such as epilepsy, abnormal gait, tremor, and clonus, which extends the clinical spectrum of the disorder. In conclusion, we confirm the causative role of variants in CXorf56 for an X-linked form of intellectual disability with additional neurological features. The gene should be considered for molecular diagnostics of patients with ID, specifically when family history is suggestive of X-linked inheritance. Further work is needed to understand the role of this gene in neurodevelopment and intellectual disability.

Loss of UGP2 in brain leads to a severe epileptic encephalopathy, emphasizing that bi-allelic isoform-specific start-loss mutations of essential genes can cause genetic diseases.

Developmental and/or epileptic encephalopathies (DEEs) are a group of devastating genetic disorders, resulting in early-onset, therapy-resistant seizures and developmental delay. Here we report on 22 individuals from 15 families presenting with a severe form of intractable epilepsy, severe developmental delay, progressive microcephaly, visual disturbance and similar minor dysmorphisms. Whole exome sequencing identified a recurrent, homozygous variant (chr2:64083454A > G) in the essential UDP-glucose pyrophosphorylase (UGP2) gene in all probands. This rare variant results in a tolerable Met12Val missense change of the longer UGP2 protein isoform but causes a disruption of the start codon of the shorter isoform, which is predominant in brain. We show that the absence of the shorter isoform leads to a reduction of functional UGP2 enzyme in neural stem cells, leading to altered glycogen metabolism, upregulated unfolded protein response and premature neuronal differentiation, as modeled during pluripotent stem cell differentiation in vitro. In contrast, the complete lack of all UGP2 isoforms leads to differentiation defects in multiple lineages in human cells. Reduced expression of Ugp2a/Ugp2b in vivo in zebrafish mimics visual disturbance and mutant animals show a behavioral phenotype. Our study identifies a recurrent start codon mutation in UGP2 as a cause of a novel autosomal recessive DEE syndrome. Importantly, it also shows that isoform-specific start-loss mutations causing expression loss of a tissue-relevant isoform of an essential protein can cause a genetic disease, even when an organism-wide protein absence is incompatible with life. We provide additional examples where a similar disease mechanism applies.

Novel NAXE variants as a cause for neurometabolic disorder: implications for treatment.

Neurometabolic disorders are often inherited and complex disorders that result from abnormalities of enzymes important for development and function of the nervous system. Recently, biallelic mutations in NAXE (APOA1BP) were found in patients with an infantile, lethal, neurometabolic disease. Here, exome sequencing was performed in two affected sisters and their healthy parents. The best candidate, NAXE, was tested for replication in exome sequencing data from 4351 patients with neurodevelopmental disorders. Quantitative RT-PCR, western blot and form factor analysis were performed to assess NAXE expression, protein levels and to analyze mitochondrial morphology in fibroblasts. Vitamin B3 was administered to one patient. Compound heterozygous missense (c.757G>A: p.Gly253Ser) and splicing (c.665-1G>A) variants in NAXE were identified in both affected sisters. In contrast to the previously reported patients with biallelic NAXE variants, our patients showed a milder phenotype with disease onset in early adulthood with psychosis, cognitive impairment, seizures, cerebellar ataxia and spasticity. The symptoms fluctuated. Additional screening of NAXE identified three novel homozygous missense variants (p.Lys245Gln, p.Asp218Asn, p.Ile214Val) in three patients with overlapping phenotype (fluctuating disease course, respiratory insufficiency, movement disorder). Lastly, patients with the c.665-1G>A splicing variant showed a significant reduction of NAXE expression compared to control fibroblasts and undetectable NAXE protein levels compared to control fibroblasts. Based on the metabolic pathway, vitamin B3 and coenzyme Q treatment was introduced in one patient in addition to antiepileptic treatment. This combination and avoidance of triggers was associated with continuous motor and cognitive improvement. The NAXE variants identified in this study suggest a loss-of-function mechanism leading to an insufficient NAD(P)HX repair system. Importantly, symptoms of patients with NAXE variants may improve with vitamin B3/coenzyme Q administration.

Novel pathogenic variants and multiple molecular diagnoses in neurodevelopmental disorders.

BACKGROUND: Rare denovo variants represent a significant cause of neurodevelopmental delay and intellectual disability (ID). METHODS: Exome sequencing was performed on 4351 patients with global developmental delay, seizures, microcephaly, macrocephaly, motor delay, delayed speech and language development, or ID according to Human Phenotype Ontology (HPO) terms. All patients had previously undergone whole exome sequencing as part of diagnostic genetic testing with a focus on variants in genes implicated in neurodevelopmental disorders up to January 2017. This resulted in a genetic diagnosis in 1336 of the patients. In this study, we specifically searched for variants in 14 recently implicated novel neurodevelopmental disorder (NDD) genes. RESULTS: We identified 65 rare, protein-changing variants in 11 of these 14 novel candidate genes. Fourteen variants in CDK13, CHD4, KCNQ3, KMT5B, TCF20, and ZBTB18 were scored pathogenic or likely pathogenic. Of note, two of these patients had a previously identified cause of their disease, and thus, multiple molecular diagnoses were made including pathogenic/likely pathogenic variants in FOXG1 and CDK13 or in TMEM237 and KMT5B. CONCLUSIONS: Looking for pathogenic variants in newly identified NDD genes enabled us to provide a molecular diagnosis to 14 patients and their close relatives and caregivers. This underlines the relevance of re-evaluation of existing exome data on a regular basis to improve the diagnostic yield and serve the needs of our patients.

De Novo Variants in TAOK1 Cause Neurodevelopmental Disorders.

De novo variants represent a significant cause of neurodevelopmental delay and intellectual disability. A genetic basis can be identified in only half of individuals who have neurodevelopmental disorders (NDDs); this indicates that additional causes need to be elucidated. We compared the frequency of de novo variants in patient-parent trios with (n = 2,030) versus without (n = 2,755) NDDs. We identified de novo variants in TAOK1 (thousand and one [TAO] amino acid kinase 1), which encodes the serine/threonine-protein kinase TAO1, in three individuals with NDDs but not in persons who did not have NDDs. Through further screening and the use of GeneMatcher, five additional individuals with NDDs were found to have de novo variants. All eight variants were absent from gnomAD (Genome Aggregation Database). The variant carriers shared a non-specific phenotype of developmental delay, and six individuals had additional muscular hypotonia. We established a fibroblast line of one mutation carrier, and we demonstrated that reduced mRNA levels of TAOK1 could be increased upon cycloheximide treatment. These results indicate nonsense-mediated mRNA decay. Further, there was neither detectable phosphorylated TAO1 kinase nor phosphorylated tau in these cells, and mitochondrial morphology was altered. Knockdown of the ortholog gene Tao1 (Tao, CG14217) in Drosophila resulted in delayed early development. The majority of the Tao1-knockdown flies did not survive beyond the third instar larval stage. When compared to control flies, Tao1 knockdown flies revealed changed morphology of the ventral nerve cord and the neuromuscular junctions as well as a decreased number of endings (boutons). Furthermore, mitochondria in mutant flies showed altered distribution and decreased size in axons of motor neurons. Thus, we provide compelling evidence that de novo variants in TAOK1 cause NDDs.

Biallelic variants in the transcription factor PAX7 are a new genetic cause of myopathy.

PURPOSE: Skeletal muscle growth and regeneration rely on muscle stem cells, called satellite cells. Specific transcription factors, particularly PAX7, are key regulators of the function of these cells. Knockout of this factor in mice leads to poor postnatal survival; however, the consequences of a lack of PAX7 in humans have not been established. METHODS: Here, we study five individuals with myopathy of variable severity from four unrelated consanguineous couples. Exome sequencing identified pathogenic variants in the PAX7 gene. Clinical examination, laboratory tests, and muscle biopsies were performed to characterize the disease. RESULTS: The disease was characterized by hypotonia, ptosis, muscular atrophy, scoliosis, and mildly dysmorphic facial features. The disease spectrum ranged from mild to severe and appears to be progressive. Muscle biopsies showed the presence of atrophic fibers and fibroadipose tissue replacement, with the absence of myofiber necrosis. A lack of PAX7 expression was associated with satellite cell pool exhaustion; however, the presence of residual myoblasts together with regenerating myofibers suggest that a population of PAX7-independent myogenic cells partially contributes to muscle regeneration. CONCLUSION: These findings show that biallelic variants in the master transcription factor PAX7 cause a new type of myopathy that specifically affects satellite cell survival.

Development of an evidence-based algorithm that optimizes sensitivity and specificity in ES-based diagnostics of a clinically heterogeneous patient population.

PURPOSE: Next-generation sequencing (NGS) is rapidly replacing Sanger sequencing in genetic diagnostics. Sensitivity and specificity of NGS approaches are not well-defined, but can be estimated from applying NGS and Sanger sequencing in parallel. Utilizing this strategy, we aimed at optimizing exome sequencing (ES)-based diagnostics of a clinically diverse patient population. METHODS: Consecutive DNA samples from unrelated patients with suspected genetic disease were exome-sequenced; comparatively nonstringent criteria were applied in variant calling. One thousand forty-eight variants in genes compatible with the clinical diagnosis were followed up by Sanger sequencing. Based on a set of variant-specific features, predictors for true positives and true negatives were developed. RESULTS: Sanger sequencing confirmed 81.9% of ES-derived variants. Calls from the lower end of stringency accounted for the majority of the false positives, but also contained ~5% of the true positives. A predictor incorporating three variant-specific features classified 91.7% of variants with 100% specificity and 99.75% sensitivity. Confirmation status of the remaining variants (8.3%) was not predictable. CONCLUSIONS: Criteria for variant calling in ES-based diagnostics impact on specificity and sensitivity. Confirmatory sequencing for a proportion of variants, therefore, remains a necessity. Our study exemplifies how these variants can be defined on an empirical basis.

Clinical exome sequencing: results from 2819 samples reflecting 1000 families.

We report our results of 1000 diagnostic WES cases based on 2819 sequenced samples from 54 countries with a wide phenotypic spectrum. Clinical information given by the requesting physicians was translated to HPO terms. WES processes were performed according to standardized settings. We identified the underlying pathogenic or likely pathogenic variants in 307 families (30.7%). In further 253 families (25.3%) a variant of unknown significance, possibly explaining the clinical symptoms of the index patient was identified. WES enabled timely diagnosing of genetic diseases, validation of causality of specific genetic disorders of PTPN23, KCTD3, SCN3A, PPOX, FRMPD4, and SCN1B, and setting dual diagnoses by detecting two causative variants in distinct genes in the same patient. We observed a better diagnostic yield in consanguineous families, in severe and in syndromic phenotypes. Our results suggest that WES has a better yield in patients that present with several symptoms, rather than an isolated abnormality. We also validate the clinical benefit of WES as an effective diagnostic tool, particularly in nonspecific or heterogeneous phenotypes. We recommend WES as a first-line diagnostic in all cases without a clear differential diagnosis, to facilitate personal medical care.

Validation of a semiconductor next-generation sequencing assay for the clinical genetic screening of CFTR.

Genetic testing for cystic fibrosis and CFTR-related disorders mostly relies on laborious molecular tools that use Sanger sequencing to scan for mutations in the CFTR gene. We have explored a more efficient genetic screening strategy based on next-generation sequencing (NGS) of the CFTR gene. We validated this approach in a cohort of 177 patients with previously known CFTR mutations and polymorphisms. Genomic DNA was amplified using the Ion AmpliSeq™ CFTR panel. The DNA libraries were pooled, barcoded, and sequenced using an Ion Torrent PGM sequencer. The combination of different robust bioinformatics tools allowed us to detect previously known pathogenic mutations and polymorphisms in the 177 samples, without detecting spurious pathogenic calls. In summary, the assay achieves a sensitivity of 94.45% (95% CI: 92% to 96.9%), with a specificity of detecting nonvariant sites from the CFTR reference sequence of 100% (95% CI: 100% to 100%), a positive predictive value of 100% (95% CI: 100% to 100%), and a negative predictive value of 99.99% (95% CI: 99.99% to 100%). In addition, we describe the observed allelic frequencies of 94 unique definitely and likely pathogenic, uncertain, and neutral CFTR variants, some of them not previously annotated in the public databases. Strikingly, a seven exon spanning deletion as well as several more technically challenging variants such as pathogenic poly-thymidine-guanine and poly-thymidine (poly-TG-T) tracts were also detected. Targeted NGS is ready to substitute classical molecular methods to perform genetic testing on the CFTR gene.

EcmPred: prediction of extracellular matrix proteins based on random forest with maximum relevance minimum redundancy feature selection.

The extracellular matrix (ECM) is a major component of tissues of multicellular organisms. It consists of secreted macromolecules, mainly polysaccharides and glycoproteins. Malfunctions of ECM proteins lead to severe disorders such as marfan syndrome, osteogenesis imperfecta, numerous chondrodysplasias, and skin diseases. In this work, we report a random forest approach, EcmPred, for the prediction of ECM proteins from protein sequences. EcmPred was trained on a dataset containing 300 ECM and 300 non-ECM and tested on a dataset containing 145 ECM and 4187 non-ECM proteins. EcmPred achieved 83% accuracy on the training and 77% on the test dataset. EcmPred predicted 15 out of 20 experimentally verified ECM proteins. By scanning the entire human proteome, we predicted novel ECM proteins validated with gene ontology and InterPro. The dataset and standalone version of the EcmPred software is available at http://www.inb.uni-luebeck.de/tools-demos/Extracellular_matrix_proteins/EcmPred.

RSARF: prediction of residue solvent accessibility from protein sequence using random forest method.

Prediction of protein structure from its amino acid sequence is still a challenging problem. The complete physicochemical understanding of protein folding is essential for the accurate structure prediction. Knowledge of residue solvent accessibility gives useful insights into protein structure prediction and function prediction. In this work, we propose a random forest method, RSARF, to predict residue accessible surface area from protein sequence information. The training and testing was performed using 120 proteins containing 22006 residues. For each residue, buried and exposed state was computed using five thresholds (0%, 5%, 10%, 25%, and 50%). The prediction accuracy for 0%, 5%, 10%, 25%, and 50% thresholds are 72.9%, 78.25%, 78.12%, 77.57% and 72.07% respectively. Further, comparison of RSARF with other methods using a benchmark dataset containing 20 proteins shows that our approach is useful for prediction of residue solvent accessibility from protein sequence without using structural information. The RSARF program, datasets and supplementary data are available at http://caps.ncbs.res.in/download/pugal/RSARF/.