An Integrated Framework for Analysis and Prediction of Impact of Single Nucleotide Polymorphism Associated with Human Diseases.

Single nucleotide polymorphisms are most common type of genetic variation in human genome. Analyzing genetic variants can help us better understand the genetic basis of diseases and develop predictive models which are useful to identify individuals who are at increased risk for certain diseases. Several SNP analysis tools have already been developed. For running these tools, the user needs to collect data from various databases. Secondly, often researchers have to use multiple variant analysis tools for cross validating their results and increase confidence in their findings. Extracting data from multiple databases and running multiple tools at a time, increases complexity and time required for analysis. There are some web-based tools that integrate multiple genetic variant databases and provide variant annotations for a few tools. These approaches have some limitations such as retrieving annotation information, filtering common pathogenic variants. The proposed web-based tool, namely IPSNP: An Integrated Platform for Predicting Impact of SNPs is written in Django which is a python-based framework. It uses RESTful API of MyVariant.info to extract annotation information of variants associated with a given gene, rsID, HGVS format variants specified in a VCF file for 29 tools. The results are in the form of a CSV file of predictions (1) derived from the consensus decision, (2) a file having annotations for the variants associated with the given gene, (3) a file showing variants declared as pathogenic commonly by the selected tools, and (4) a CSV file containing chromosome coordinates based on GRCh37 and GRCh38 genome assemblies, rsIDs and proteomic data, so that users may use tools of their choice and avoiding manual parameter collection for each tool. IPSNP is a valuable resource for researchers and clinicians and it can help to save time and effort in discovering the novel disease-associated variants and the development of personalized treatments.

A Missense Pathogenic Variant in a Conserved Region of CNTNAP2 Is Associated with Obesity, Seizures, and Language Impairment in a Pakistani Family.

Introduction: In a consanguineous family, seven siblings born in three sibships showed a syndromic disorder characterized by obesity, seizures, and language impairment phenotypes, which appeared at early age or developed during early childhood. Methods: By whole-exome sequencing and subsequent Sanger sequencing, a novel homozygous missense variant (c.3371 T>A [p.Ile1124Asn]) in exon 20 of the CNTNAP2 gene was identified. Results: The pathogenic variant in this family is located within one of the laminin G-like 4 domains of CASPR2 and may cause loss of hydrophobic interactions of CASPR2 with its partner proteins. Single nucleotide and copy number variants in this gene have previously been related to Gilles de la Tourette syndrome, cortical dysplasia-focal epilepsy syndrome, schizophrenia, Pitt-Hopkins syndrome, and autism spectrum, attention deficit hyperactivity, and obsessive compulsive disorders. Yet, few studies described patients with CNTNAP2 variants showing diet-induced obesity. Conclusion: This report expands the phenotypic spectrum of this rare syndrome and provides deeper insights by documenting the clinical features and genetic findings of the patients.

A Comprehensive Review of Performance of Next-Generation Sequencing Platforms.

Background: Next-generation sequencing methods have been developed and proposed to investigate any query in genomics or clinical activity involving DNA. Technical advancement in these sequencing methods has enhanced sequencing volume to several billion nucleotides within a very short time and low cost. During the last few years, the usage of the latest DNA sequencing platforms in a large number of research projects helped to improve the sequencing methods and technologies, thus enabling a wide variety of research/review publications and applications of sequencing technologies. Objective: The proposed study is aimed at highlighting the most fast and accurate NGS instruments developed by various companies by comparing output per hour, quality of the reads, maximum read length, reads per run, and their applications in various domains. This will help research institutions and biological/clinical laboratories to choose the sequencing instrument best suited to their environment. The end users will have a general overview about the history of the sequencing technologies, latest developments, and improvements made in the sequencing technologies till now. Results: The proposed study, based on previous studies and manufacturers' descriptions, highlighted that in terms of output per hour, Nanopore PromethION outperformed all sequencers. BGI was on the second position, and Illumina was on the third position. Conclusion: The proposed study investigated various sequencing instruments and highlighted that, overall, Nanopore PromethION is the fastest sequencing approach. BGI and Nanopore can beat Illumina, which is currently the most popular sequencing company. With respect to quality, Ion Torrent NGS instruments are on the top of the list, Illumina is on the second position, and BGI DNB is on the third position. Secondly, memory- and time-saving algorithms and databases need to be developed to analyze data produced by the 3rd- and 4th-generation sequencing methods. This study will help people to adopt the best suited sequencing platform for their research work, clinical or diagnostic activities.

Structural and Functional Impact of Damaging Nonsynonymous Single Nucleotide Polymorphisms (nsSNPs) on Human VPS35 Protein Using Computational Approaches.

Parkinson's disease is the second most common progressive neurodegenerative movement disorder. Mutations in retromer complex subunit and VPS35 represent the second most common cause of late-onset familial Parkinson's disease. The mutation in VPS35 can disrupt the normal protein functions resulting in Parkinson's disease. The aim of this study was the identification of deleterious missense Single Nucleotide Polymorphisms (nsSNPs) and their structural and functional impact on the VPS35 protein. In this study, several insilico tools were used to identify deleterious and disease-associated nsSNPs. 3D structure of VPS35 protein was constructed through MODELLER 9.2, normalized using FOLDX, and evaluated through RAMPAGE and ERRAT whereas, FOLDX was used for mutagenesis. 25 ligands were obtained from literature and docked using PyRx 0.8 software. Based on the binding affinity, five ligands i.e., PG4, MSE, GOL, EDO, and CAF were further analyzed. Molecular Dynamic simulation analysis was performed using GROMACS 5.1.4, where temperature, pressure, density, RMSD, RMSF, Rg, and SASA graphs were analyzed. The results showed that the mutations Y67H, R524W, and D620N had a structural and functional impact on the VPS35 protein. The current findings will help in appropriate drug design against the disease caused by these mutations in a large population using in-vitro study.

Compound Homozygous Rare Mutations in PLCE1 and HPS1 Genes Associated with Autosomal Recessive Retinitis Pigmentosa in Pakistani Families.

Background: Retinitis pigmentosa (RP) belongs to pigmentary retinopathies, a generic name for all retinal dystrophies with a major phenotypical and genotypical variation, characterized by progressive reduction of photo-receptor functionality of the rod and cone. Global prevalence of RP is ~ 1/4000 and it can be inherited as autosomal dominant (adRP), autosomal recessive (arRP) or X- linked (xlRP). We designed this study to identify causative mutations in Pakistani families affected with arRP. Methods: In 2019, we recruited two unrelated Pakistani consanguineous families affected with progressive vision loss and night blindness from Punjab region. Clinical diagnosis confirmed the; bone spicule pigmentation of the retina, and an altered electroretinogram (EGR) response. Proband and healthy individual from each family were subjected for whole-exome sequencing (WES). Various computational tools were used to analyze the Next Generation Sequencing (NGS) data and to predict the pathogenicity of the identified mutations. Results: WES data analysis highlighted two missense homozygous variants at position c.T1405A (p.S469T) in PLCE1 and c.T11C (p.V4A) in HPS1 genes in proband of both families. Healthy individuals of two families were tested negative for p.S469T and p.V4A mutations. The variant analysis study including molecular dynamic simulations predicted mutations as disease causing. Conclusion: Compound effect of mutations in rarely linked PLCE1 and HPS1 genes could also cause RP. This study highlights the potential application of WES for a rapid and precise molecular diagnosis for heterogeneous genetic diseases such as RP.

In Silico Analysis of Hepatitis B Virus Genotype D Subgenotype D1 Circulating in Pakistan, China, and India.

The focus of this study was the computational analysis of hepatitis B virus (HBV) genotype D subgenotype D1 in Pakistan, China, and India. In total, 54 complete genome sequences of HBV genotype D subgenotype D1 were downloaded from National Center for Biotechnology Information (NCBI). Of these, 6 complete genome sequences were from Pakistan, 14 were from China, and 34 were from India. Sequence alignment showed less than 4% divergence in these sequences. C and X genes showed divergence of less than 3%. Comparison over the S gene showed more than 97% similarity among the nucleotide sequences of genotype D subgenotype D1. The identity and similarity matrix of 54 nucleotide sequences of HBV genotype D subgenotype D1 from Pakistan, China, and India revealed more than 93% identity and 93% similarity. Phylogenetic analysis highlighted that complete genome isolates of HBV circulating in Pakistan had the closest evolutionary relationship with its neighboring countries China and India. China's (HQ833466) and Pakistan's (AB583680.1) isolates shared the same ancestor. Gene structure analysis showed that "P" gene exons were the longest, about three-fourth of the genome size, whereas gene "S" had the second longest coding regions with 2 exons and 1 intron. However, "C" and "X" genes had 1 smallest exon. X proteins had proven role in spreading of the HBV infection diseases. For HBx analysis, 1 X protein sequence of HBV genotype D subgenotype D1 belonging to each country was obtained. Homology models of the 3 X proteins generated using SWISS-MODEL revealed GMQE (Global Model Quality Estimation) = 0.1. Global and local quality estimate scores including Z-scores for Qualitative Model Energy Analysis (QMEAN) C-beta, all-atom, solvation, and torsion energy scores were similar indicating good quality, accuracy, and reliability of the predicted models. Three-dimensional (3D) visualization showed similar structures and Ramachandran plots showed a high percentage of protein residues into the favorable region for X protein models.

SAliBASE: A Database of Simulated Protein Alignments.

Simulated alignments are alternatives to manually constructed multiple sequence alignments for evaluating performance of multiple sequence alignment tools. The importance of simulated sequences is recognized because their true evolutionary history is known, which is very helpful for reconstructing accurate phylogenetic trees and alignments. However, generating simulated alignments require expertise to use bioinformatics tools and consume several hours for reconstructing even a few hundreds of simulated sequences. It becomes a tedious job for an end user who needs a few datasets of variety of simulated sequences. Currently, there is no databank available which may help researchers to download simulated sequences/alignments for their study. Major focus of our study was to develop a database of simulated protein sequences (SAliBASE) based on different varying parameters such as insertion rate, deletion rate, sequence length, number of sequences, and indel size. Each dataset has corresponding alignment as well. This repository is very useful for evaluating multiple alignment methods.

A Comprehensive Study of De Novo Genome Assemblers: Current Challenges and Future Prospective.

BACKGROUND: Current advancements in next-generation sequencing technology have made possible to sequence whole genome but assembling a large number of short sequence reads is still a big challenge. In this article, we present the comparative study of seven assemblers, namely, ABySS, Velvet, Edena, SGA, Ray, SSAKE, and Perga, using prokaryotic and eukaryotic paired-end as well as single-end data sets from Illumina platform. RESULTS: Results showed that in case of single-end data sets, Velvet and ABySS outperformed in all the seven assemblers with comparatively low assembling time and high genome fraction. Velvet consumed the least amount of memory than any other assembler. In case of paired-end data sets, Velvet consumed least amount of time and produced high genome fraction after ABySS and Ray. In terms of low memory usage, SGA and Edena outperformed in all the assemblers. Ray also showed good genome fraction; however, extremely high assembling time consumed by the Ray might make it prohibitively slow on larger data sets of single and paired-end data. CONCLUSIONS: Our comparison study will provide assistance to the scientists for selecting the suitable assembler according to their data sets and will also assist the developers to upgrade or develop a new assembler for de novo assembling.

MPD3: a useful medicinal plants database for drug designing.

Medicinal plants are the main natural pools for the discovery and development of new drugs. In the modern era of computer-aided drug designing (CADD), there is need of prompt efforts to design and construct useful database management system that allows proper data storage, retrieval and management with user-friendly interface. An inclusive database having information about classification, activity and ready-to-dock library of medicinal plant's phytochemicals is therefore required to assist the researchers in the field of CADD. The present work was designed to merge activities of phytochemicals from medicinal plants, their targets and literature references into a single comprehensive database named as Medicinal Plants Database for Drug Designing (MPD3). The newly designed online and downloadable MPD3 contains information about more than 5000 phytochemicals from around 1000 medicinal plants with 80 different activities, more than 900 literature references and 200 plus targets. The designed database is deemed to be very useful for the researchers who are engaged in medicinal plants research, CADD and drug discovery/development with ease of operation and increased efficiency. The designed MPD3 is a comprehensive database which provides most of the information related to the medicinal plants at a single platform. MPD3 is freely available at: http://bioinform.info .

A Comparative Study of Land Cover Classification by Using Multispectral and Texture Data.

The main objective of this study is to find out the importance of machine vision approach for the classification of five types of land cover data such as bare land, desert rangeland, green pasture, fertile cultivated land, and Sutlej river land. A novel spectra-statistical framework is designed to classify the subjective land cover data types accurately. Multispectral data of these land covers were acquired by using a handheld device named multispectral radiometer in the form of five spectral bands (blue, green, red, near infrared, and shortwave infrared) while texture data were acquired with a digital camera by the transformation of acquired images into 229 texture features for each image. The most discriminant 30 features of each image were obtained by integrating the three statistical features selection techniques such as Fisher, Probability of Error plus Average Correlation, and Mutual Information (F + PA + MI). Selected texture data clustering was verified by nonlinear discriminant analysis while linear discriminant analysis approach was applied for multispectral data. For classification, the texture and multispectral data were deployed to artificial neural network (ANN: n-class). By implementing a cross validation method (80-20), we received an accuracy of 91.332% for texture data and 96.40% for multispectral data, respectively.

The accuracy of protein structure alignment servers

Background: Protein structural alignment is one of the most fundamental and crucial areas of research in the domain of computational structural biology. Comparison of a protein structure with known structures helps to classify it as a new or belonging to a known group of proteins. This, in turn, is useful to determine the function of protein, its evolutionary relationship with other protein molecules and grasping principles underlying protein architecture and folding. Results: A large number of protein structure alignment methods are available. Each protein structure alignment tool has its own strengths and weaknesses that need to be highlighted. We compared and presented results ofsix most popular and publically available servers for protein structure comparison. These web-based servers were compared with the respect to functionality (features provided by these servers) and accuracy (how well the structural comparison is performed). The CATH was used as a reference. The results showed that overall CE was top performer. DALI and PhyreStorm showed similar results whereas PDBeFold showed the lowest performance. In case of few secondary structural elements, CE, DALI and PhyreStorm gave 100% success rate. Conclusion: Overall none of the structural alignment servers showed 100% success rate. Studies of overall performance, effect of mainly alpha and effect of mainly beta showed consistent performance. CE, DALI, FatCat and PhyreStorm showed more than 90% success rate.

IVisTMSA: Interactive Visual Tools for Multiple Sequence Alignments.

IVisTMSA is a software package of seven graphical tools for multiple sequence alignments. MSApad is an editing and analysis tool. It can load 409% more data than Jalview, STRAP, CINEMA, and Base-by-Base. MSA comparator allows the user to visualize consistent and inconsistent regions of reference and test alignments of more than 21-MB size in less than 12 seconds. MSA comparator is 5,200% efficient and more than 40% efficient as compared to BALiBASE c program and FastSP, respectively. MSA reconstruction tool provides graphical user interfaces for four popular aligners and allows the user to load several sequence files at a time. FASTA generator converts seven formats of alignments of unlimited size into FASTA format in a few seconds. MSA ID calculator calculates identity matrix of more than 11,000 sequences with a sequence length of 2,696 base pairs in less than 100 seconds. Tree and Distance Matrix calculation tools generate phylogenetic tree and distance matrix, respectively, using neighbor joining% identity and BLOSUM 62 matrix.

Evaluating the accuracy and efficiency of multiple sequence alignment methods.

A comparison of 10 most popular Multiple Sequence Alignment (MSA) tools, namely, MUSCLE, MAFFT(L-INS-i), MAFFT (FFT-NS-2), T-Coffee, ProbCons, SATe, Clustal Omega, Kalign, Multalin, and Dialign-TX is presented. We also focused on the significance of some implementations embedded in algorithm of each tool. Based on 10 simulated trees of different number of taxa generated by R, 400 known alignments and sequence files were constructed using indel-Seq-Gen. A total of 4000 test alignments were generated to study the effect of sequence length, indel size, deletion rate, and insertion rate. Results showed that alignment quality was highly dependent on the number of deletions and insertions in the sequences and that the sequence length and indel size had a weaker effect. Overall, ProbCons was consistently on the top of list of the evaluated MSA tools. SATe, being little less accurate, was 529.10% faster than ProbCons and 236.72% faster than MAFFT(L-INS-i). Among other tools, Kalign and MUSCLE achieved the highest sum of pairs. We also considered BALiBASE benchmark datasets and the results relative to BAliBASE- and indel-Seq-Gen-generated alignments were consistent in the most cases.