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 3^rd‐ and 4^th‐generation sequencing methods. This study will help people to adopt the best suited sequencing platform for their research work, clinical or diagnostic activities.

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.

Single Nucleotide Polymorphisms (SNPs) are the most common candidate mutations in human beings that play a vital role in the genetic basis of certain diseases. Previous studies revealed that Solute Carrier Family 26 Member 4 (SLC26A4) being an essential gene of the multi-faceted transporter family SLC26 facilitates reflexive movement of Iodide into follicular lumen through apical membrane of thyrocyte. SLC26A4 gene encodes Pendred protein, a membrane glycoprotein, highly hydrophobic in nature, present at the apical membrane of thyrocyte functioning as transporter of iodide for thyroid cells. A minor genetic variation in SLC26A4 can cause Pendred syndrome, a syndrome associated with thyroid glands and deafness. In this study, we performed in-silico analysis of 674 missense SNPs of SLC26A4 using different computational platforms. The bunch of tools including SNPNEXUS, SNAP-2, PhD-SNP, SNPs&GO, I-Mutant, ConSurf, and ModPred were used to predict 23 highly confident damaging and disease causing nsSNPs (G209V, G197R, L458P, S427P, Q101P, W472R, N392Y, V359E, R409C, Q235R, R409P, G139V, G497S, H723R, D87G, Y127H, F667C, G334A, G95R, S427C, R291W, Q383H and E384G) that could potentially alter the SLC26A4 gene. Moreover, protein structure prediction, protein-ligand docking and Molecular Dynamics simulation were performed to confirm the impact of two evident alterations (Y127H and G334A) on the protein structure and function.

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.pk M. T. Pervaiz is a co-corresponding author. # Authors have an equal contribution.

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 photoreceptor 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.

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.