Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing.

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

Double-Stranded Break Repair in Mammalian Cells and Precise Genome Editing.

In mammalian cells, double-strand breaks (DSBs) are repaired predominantly by error-prone non-homologous end joining (NHEJ), but less prevalently by error-free template-dependent homologous recombination (HR). DSB repair pathway selection is the bedrock for genome editing. NHEJ results in random mutations when repairing DSB, while HR induces high-fidelity sequence-specific variations, but with an undesirable low efficiency. In this review, we first discuss the latest insights into the action mode of NHEJ and HR in a panoramic view. We then propose the future direction of genome editing by virtue of these advancements. We suggest that by switching NHEJ to HR, full fidelity genome editing and robust gene knock-in could be enabled. We also envision that RNA molecules could be repurposed by RNA-templated DSB repair to mediate precise genetic editing.

Myostatin increases the expression of matrix metalloproteinase genes to promote preadipocytes differentiation in pigs.

ABSTACTMyostatin (MSTN) resulted in reduced backfat thickness in MSTN-knockout (MSTN-KO) pigs, whereas the underlying mechanism remains elusive. In this study, RNA sequencing (RNA-seq) was used to screen differentially expressed genes (DEGs) in porcine fat tissues. We identified 285 DEGs, including 4 adipocyte differentiation-related genes (ADRGs). Matrix Metalloproteinase-2/7 (MMP-2/7), fibronectin (FN), and laminin (LN) were differentially expressed in MSTN-KO pigs compared with wild-type (WT) pigs. To investigate the molecular mechanism, we treated the preadipocytes with siRNA and recombinant MSTN protein. The results indicated that MSTN increased the expression of MMP-2/7/9 and promoted the preadipocyte differentiation. To further validate the effect of MSTN on MMP-2/7/9 expression, we treated MSTN-KO PK15 cells with recombinant MSTN protein and detected the expression of MMP-2/7/9. The data showed that MSTN increases the expression of MMP-2/7/9 in PK15. This study revealed that MSTN promoted preadipocyte differentiation and provided the basis for the mechanism of fatty deposition in pigs.

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.

The effects of some heavy metals on some fish species.

The current research was designed to investigate the adverse effects of dense ores in different organs exposed to these metals polluted water on three edible fish species (Wallago attu, Catla catla and Tilapia nilotica), sampled from River Sutlej at Head Islam (Hasilpur), Pakistan. The assimilation of different elements Ni, Fe, Cd, Cr and Pb in body parts of normal fish took from the fish farm and also made a comparison among understudied fishes. Overall absorption of ores in fishes was ranked as Fe > Ni > Cr > Cd > Pb. A remarkable variation in the deposition of elements in different organs of three fish species and was analyzed statistically (p < 0.05). High levels of ore assimilation in fish organs pointed out potential health risks for the fish, other aquatic organisms and human beings too.

Insilico modeling and analysis of small molecules binding to the PHLPP1 protein by molecular dynamics simulation.

The PHLPP (Pleckstrin homology domain leucine-rich repeat protein phosphatases) is a newly discovered group of genes which includes PHLPP1 and PHLPP2 and plays an integral part in several cellular processes like apoptosis, cell signaling cell survival, and cell proliferation etc. Both the activation and deactivation of these genes can have vital role in several ailments like heart diseases, circadian rhythm and most importantly the cancer, hence encouraging the growth of novel therapeutic elements. To give new directions into the development of PHLPP1- targeting drugs, the interaction mechanism between PHLPP1 and five important ligands 4IP, B39, 635, ATP and GTA were investigated through docking and Molecular Dynamics Simulation. It is also noteworthy to be mentioned here that there is no previous crystal structure of PHLPP1 available. The in-silico results can provide potential base for advancements in development of new therapeutic elements targeting different diseases, mainly cancer. In this study, we employed homology modeling technique to develop a high-quality structure model of PHLPP1. The PHLPP1 model was then used in docking interaction analysis and Molecular Dynamics Simulation, to study binding pockets and interactions of PHLPP1 ligands and finding actively contributing residues in binding pocket. In final step, Free Energy Estimation was performed to observe ligand binding's quantitative characteristics.

Isolation and genetic characterization of virulent strains of avian paramyxovirus-1 from multiple avian species in Azad Jammu and Kashmir 2017-2018.

Despite intensive vaccination, endemicity of Avian paramyxoviruses-1 (APMV-1) is a significant problem in developing countries in Africa, Middle East, and Asia. Given the importance of APMV-1 in poultry and multiple non-poultry avian species, it is important to continue surveillance programs, routine monitoring and characterization of field isolates in the region where viruses are endemic. The purpose of this study was to pathotyped and genetically characterized 21 APMV-1s isolated from multiple avian species reared in different regions of Azad Jammu and Kashmir (AJK). Phylogenetic analysis based on complete fusion (F) gene sequences showed that 17 APMV-1 isolates obtained from commercial poultry and backyard birds belonged to sub-genotype VIIi. Though, one pigeon-origin APMV-1 isolate was clustered in sub-genotype VIg and three in recently designated new sub-genotype VIm of genotype VI. The pigeon-origin isolates had the following two motifs 113-RKKR↓F-117 and 113-RQRR↓F-117, while all other isolates had the polybasic amino acid sequence 113-RQKR↓F-117 at the F-cleavage site, which is characteristic of virulent APMV-1 strains. These results are consistent with the five viruses that had intracerebral pathogenicity indices (ICPIs) of between 1.50 and 1.73, corresponding to a velogenic pathotype. The APMV-1s isolated from commercial poultry and backyard birds in this study showed low nucleotide distance (0.3-0.9%) and genetically closely related (> 97%) to viruses repeatedly isolated (2011-2017) from multiple avian species in other states of Pakistan. Strengthened surveillance programs in both commercial poultry and backyard flocks are needed to better assess the commercial-backyard bird interface and form a basis for evidence-based measures to limit and prevent APMV-1 transmission.

Isolation and Characterization of Avian Influenza H9N2 Viruses from Different Avian Species in Pakistan 2016-17.

Avian influenza (AI) virus (AIV) subtype H9N2 continues to cause significant outbreaks among commercial and backyard poultry in Pakistan. Despite this, the characterization of H9N2 viruses in avian hosts other than chickens in Pakistan has not been thoroughly investigated. In this study, 12 low pathogenicity avian influenza viruses subtype H9N2 were isolated from peacocks (n = 4), ducks (n = 4), pheasants (n = 2), geese (n = 1), and black swans (n = 1) in Pakistan during 2016 and were characterized on the basis of the hemagglutinin (HA) and neuraminidase genes. All of the viruses possessed an amino acid substitution Q226L in the receptor-binding site of the HA protein, which is known to contribute to increased viral replication and virulence in mammals. In addition, phylogenetic studies showed that these H9N2 AIVs belonged to the Middle East B genetic group of sublineage G1 and were very similar to viruses isolated from an outbreak in chickens in Pakistan in 2017. This demonstrates an epidemiologic link between poultry and other avian species, which is a fact to consider in future H9N2 disease management programs.

Nota de investigación ­ Aplicación in ovo y al primer día de edad de una vacuna viva contra el virus de la enfermedad infecciosa de la bolsa en pollos de engorde comerciales. La enfermedad infecciosa de la bolsa (IBD) es una enfermedad económicamente importante de pollos jóvenes causada por un Avibirnavirus, el virus de la enfermedad infecciosa de la bolsa (IBDV). El virus causal es altamente resistente a los ambientes avícolas y la vacunación es la medida más efectiva para el control de esta enfermedad. Sin embargo, la presunta neutralización de las cepas altamente atenuadas por los anticuerpos maternos y la supuesta virulencia de las cepas parcialmente atenuadas ha limitado la implementación de cepas de vacunas vivas convencionales contra la enfermedad de Gumboro en pollitos antes o después de la eclosión. Sin embargo, los datos preliminares han generado dudas sobre la validez de este dogma prevaleciente. Para analizar la posible aplicación de una cepa de vacuna viva intermedia plus del virus de Gumboro, cepa MB-1, en embriones de pollo con inmunidad materna y a pollos de un día de edad, se realizaron cuatro ensayos de campo a gran escala en distintos lugares del mundo. Los cuatro ensayos midieron la seguridad relativa, los parámetros de inmunización para Gumboro y los rendimientos de producción de MB-1 frente a las vacunas de Gumboro con complejos inmunes establecidas en una variedad de sistemas comerciales de engorde. El desempeño general de la salud y la producción en los cuatro ensayos ha sido similar o mejor en los grupos con la cepa MB-1. Los resultados desafían la noción prevaleciente de que las cepas vivas de del virus de la enfermedad infecciosa de la bolsa pueden neutralizarse o sobrepasar la inmunidad materna e inducir daño permanente a la respuesta inmune de los pollos de engorde jóvenes. Se ha observado un fenómeno de replicación retardada después de la administración parenteral de la cepa de vacuna viva contra el virus de la enfermedad de Gumboro, sin embargo, este mecanismo de replicación retardada aún no se ha dilucidado. Los resultados de este estudio justifican una mayor investigación de las cepas de vacunas vivas convencionales contra la enfermedad de Gumboro como una alternativa para la inmunización activa de pollos de engorde antes y después de la eclosión.

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.

Review - MicroRNAs: A new paradigm towards mechanistic insight of diseases.

In 1993 miRNAs were discovered during a research on Caenorhabditis elegans conducted by Victor Ambros and Gary Ruvkun. The gene lin-4 that played important role in development in C. elgans was observed not encoding any protein but a very small RNA molecule of just 22 nucleotides. Main objective of this review is to highlight the significance of miRNAs in regulating the expression of many genes, which are either directly or indirectly involved in many diseases. One of the major causes of illness and death in developed countries of the world is cardiovascular disease. Some of the miRNAs have certain role to play in heart that are not specified for heart. So miRNAs have been found to be in other tissues like fibroblasts, endothelial cells and smooth muscle cells that are part of physiological study of cardiovascular system. Adult heart has limited capacity of regeneration therefore lost cardiomyocytes due to myocardial ischemia or infarction can result in low performance of heart. miRNAs have been shown to play a role in apoptotic regulation of cardiomyocytes in vivo. Many studies have shown that miR146a and 155 are up regulated in peripheral blood mononuclear cells, synovial fibroblasts, synovial fluid and Th-17 cells from rheumatoid arthritis patients as compared to healthy persons. Several types of miRNAs are playing important roles in type 1 diabetes mellitus including miR-375 and miR-375 with intolerance to glucose and decreased beta cells account due to impaired proliferation. Up regulation of miR-125a in WAT of type 2 Diabetes mellitus have been observed. miRNAs have proved to be the important regulators of cytokines and growth factor expression. Thus, suggested as a good biomarker and target of therapy. miRNA profiling techniques have revealed the role of miRNAs in Multiple sclerosis.

Pathotyping and genetic characterization of avian avulavirus-1 from domestic and wild waterfowl, geese and black swans in Pakistan, 2014 to 2017.

Twenty-nine avian avulavirus-1 viruses (AAvV-1s) from healthy domestic and wild ducks, geese and black swans collected in Pakistan between 2014-2017 have been pathotyped and genetically characterized. A phylogenetic analysis revealed that 21 of the isolates belonged to sub-genotype VIIi, whereas eight isolates were highly similar to vaccine-like viruses of genotype II. In addition to confirming the continued presence of sub-genotype VIIi AAvV-1s in Pakistan, this study identifies the probable spill-over of vaccine-like viruses from vaccinated poultry to wild and domestic waterfowl and, as such, has important implications for the control and management of Newcastle disease in Pakistan.

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.

Osteopontin gene polymorphism association with milk traits and its expression analysis in milk of riverine buffalo.

Osteopontin gene is regarded as a plausible candidate in mammary gland differentiation and development, expressed by variety of cells, tissues, and biological fluids including milk. The current study was performed in two phases. In the first phase, Osteopontin gene polymorphisms were identified and associated with milk composition such as ash, milk fat, SNF, lactose, and protein. In the second phase, milk samples from five healthy mastitis-free Nili Ravi buffaloes were analyzed for expression of Osteopontin gene at transition (day 15), mid (day 90), and end (day 250) stage of their second lactation. Briefly, blood samples were collected from Nili Ravi buffalo to isolate the genomic DNA, specific primers were designed for PCR amplification. The amplified PCR products were sequenced bi-directionally. Six polymorphisms were identified in the coding region and four in the intronic region of the gene. The results showed that SNP g.38329758 T > C causing substitution of valine to alanine (V127A) was associated with high milk protein. For mRNA expression analysis, somatic cells were separated from milk samples for RNA isolation. Analysis of differential gene expression data has permitted us to illustrate the expression pattern of osteopontin gene in lactating buffalo. The Osteopontin gene was found to be transcribed among all three lactation stages, but expression was observed with the highest value (fold change) in peak lactation and remained elevated till the end of lactation. Identified gene marker may be helpful for the prediction of superior animal for selection. The presented study also gave an insight into the genetic screening and lactation biology of riverine buffalo, offering direction for future research in lactating buffalo.

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.

Polymorphism analysis of prion protein gene in 11 Pakistani goat breeds.

The association between caprine PrP gene polymorphisms and its susceptibility to scrapie has been investigated in current years. As the ORF of the PrP gene is extremely erratic in different breeds of goats, we studied the PrP gene polymorphisms in 80 goats which belong to 11 Pakistani indigenous goat breeds from all provinces of Pakistan. A total of 6 distinct polymorphic sites (one novel) with amino acid substitutions were identified in the PrP gene which includes 126 (A -> G), 304 (G -> T), 379 (A -> G), 414 (C -> T), 428 (A -> G) and 718 (C -> T). The locus c.428 was found highly polymorphic in all breeds as compare to other loci. On the basis of these PrP variants NJ phylogenetic tree was constructed through MEGA6.1 which showed that all goat breeds along with domestic sheep and Mauflon sheep appeared as in one clade and sharing its most recent common ancestors (MRCA) with deer species while Protein analysis has shown that these polymorphisms can lead to varied primary, secondary and tertiary structure of protein. Based on these polymorphic variants, genetic distance, multidimensional scaling plot and principal component analyses revealed the clear picture regarding greater number of substitutions in cattle PrP regions as compared to the small ruminant species. In particular these findings may pinpoint the fundamental control over the scrapie in Capra hircus on genetic basis.

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.

Identification of functional consequence of a novel selection signature in CYP11b1 gene for milk fat content in Bubalus bubalis.

Genomic selection for traits of economic importance is an emerging approach carrying tremendous potentials. Many of polygenic traits as milk fat, protein and yield have been characterize at genomic level and important selection signatures have been identified. Cytochrome P450 enzymes are potential loci for affecting many of dairy capabilities. Present study was conducted for genomic dissection of CYP11b1 gene in riverine buffaloes and seven genetic variations were identified. Out of these, one novel polymorphism (p.A313T) was found well associated with milk fat %age. AB genotyped buffaloes were found to have higher milk fat %age (8.9%) for this loci. p.A313T was further validated at larger data set by restriction digestion using CviAII enzyme. Functional consequences of this locus were also predicted by studying three dimensional structure of CYP11b1 protein. For this purpose, 3D protein model was predicted by homology modeling, secondary structural attributes were determined, signal peptide was predicted and a transmembrane helix was also identified. One of polymorphism (p.Y205L) was found in the vicinity of functionally significant F-G loop region, which is the part of protein gets attached to the inner mitochondrial membrane. But this variation could not be associated and needs further investigation. p.A30V, a popular selection marker in cattle, was found in buffaloes as well but could not be associated and might need further confirmation on larger data set. Results of this study illustrate the impending potential of this gene in determining dairy capabilities of buffaloes and might have a role in selection of superior dairy buffaloes.

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.

Sympatric species distribution, genetic diversity and population structure of Haemonchus isolates from domestic ruminants in Pakistan.

Haemonchus species are major gastro-intestinal parasites affecting ruminants across the world. The present study aimed to assess the sympatric species distribution, genetic diversity, population structure and frequency of ß-tubulin isotype 1 alleles associated with benzimidazole resistance. Internal transcribed spacer 2 (ITS2) sequences revealed three sympatric species of Haemonchus, H. contortus, H. placei and H. longistipes with 12 distinct genotypes circulating among ruminant hosts in Pakistan. High genetic variability was observed in Pakistani Haemonchus isolates at nicotine amide dehydrogenase subunit 4 (ND4) and cytochrome oxidase subunit 1 (COI) gene loci. Intra-population diversity parameters were higher in H. contortus isolates than H. placei. Phylogenetic analysis of ND4 and COI sequences did not reveal clustering of haplotypes originating from a particular host indicating high rate of gene flow among Haemonchus parasites infecting sheep, goat and cattle in Pakistan. ND4 and COI haplotypes from Pakistan were compared to sequences of Haemonchus isolates from 11 countries to elucidate the population structure. Multidimensional scaling (MDS) plot of pairwise FST derived from 531 ND4 haplotypes revealed clustering together of H. contortus from Pakistan, China, Malaysia and Italy while the isolates from Yemen and United States were found to be genetically distinct. With respect to H. placei, isolates from Pakistan were found to be genetically differentiated from isolates of other countries. The tests for selective neutrality revealed negative D statistics and did not reveal significant deviations in Pakistani Haemonchus populations while significant deviation (P < 0.05) was observed in Brazilian and Chinese H. contortus populations. Median Joining (MJ) network of ND4 haplotypes revealed Yemenese H. contortus being closer to H. placei cluster. ß-tubulin isotype 1 genotyping revealed 7.86% frequency of Y allele associated with benzimidazole resistance at F200Y locus in Pakistani Haemonchus isolates.