Discovery and genotyping of novel sequence insertions in many sequenced individuals.

MOTIVATION: Despite recent advances in algorithms design to characterize structural variation using high-throughput short read sequencing (HTS) data, characterization of novel sequence insertions longer than the average read length remains a challenging task. This is mainly due to both computational difficulties and the complexities imposed by genomic repeats in generating reliable assemblies to accurately detect both the sequence content and the exact location of such insertions. Additionally, de novo genome assembly algorithms typically require a very high depth of coverage, which may be a limiting factor for most genome studies. Therefore, characterization of novel sequence insertions is not a routine part of most sequencing projects. RESULT: Here, we present Pamir, a new algorithm to efficiently and accurately discover and genotype novel sequence insertions using either single or multiple genome sequencing datasets. Pamir is able to detect breakpoint locations of the insertions and calculate their zygosity (i.e. heterozygous versus homozygous) by analyzing multiple sequence signatures, matching one-end-anchored sequences to small-scale de novo assemblies of unmapped reads, and conducting strand-aware local assembly. We test the efficacy of Pamir on both simulated and real data, and demonstrate its potential use in accurate and routine identification of novel sequence insertions in genome projects. AVAILABILITY AND IMPLEMENTATION: Pamir is available at https://github.com/vpc-ccg/pamir . CONTACT: fhach@{sfu.ca, prostatecentre.com } or calkan@cs.bilkent.edu.tr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comparative transcriptome profiling approach to glean virulence and immunomodulation-related genes of Fasciola hepatica.

BACKGROUND: Fasciola hepatica causes chronic liver disease, fasciolosis, leading to significant losses in the livestock economy and concerns for human health in many countries. The identification of F. hepatica genes involved in the parasite's virulence through modulation of host immune system is utmost important to comprehend evasion mechanisms of the parasite and develop more effective strategies against fasciolosis. In this study, to identify the parasite's putative virulence genes which are associated with host immunomodulation, we explored whole transcriptome of an adult F. hepatica using current transcriptome profiling approaches integrated with detailed in silico analyses. In brief, the comparison of the parasite transcripts with the specialised public databases containing sequence data of non-parasitic organisms (Dugesiidae species and Caenorhabditis elegans) or of numerous pathogens and investigation of the sequences in terms of nucleotide evolution (directional selection) and cytokine signaling relation were conducted. RESULTS: NGS of the whole transcriptome resulted in 19,534,766 sequence reads, yielding a total of 40,260 transcripts (N50 = 522 bp). A number of the parasite transcripts (n = 1,671) were predicted to be virulence-related on the basis of the exclusive homology with the pathogen-associated data, positive selection or relationship with cytokine signaling. Of these, a group of the virulence-related genes (n = 62), not previously described, were found likely to be associated with immunomodulation based on in silico functional categorisation, showing significant sequence similarities with various immune receptors (i.e. MHC I class, TGF-ß receptor, toll/interleukin-1 receptor, T-cell receptor, TNF receptor, and IL-18 receptor accessory protein), cytokines (i.e. TGF-ß, interleukin-4/interleukin-13 and TNF-α), cluster of differentiations (e.g. CD48 and CD147) or molecules associated with other immunomodulatory mechanisms (such as regulation of macrophage activation). Some of the genes (n = 5) appeared to be under positive selection (Ka/Ks > 1), imitating proteins associated with cytokine signaling (through sequence homologies with thrombospondin type 1, toll/interleukin-1 receptor, TGF-ß receptor and CD147). CONCLUSIONS: With a comparative transcriptome profiling approach, we have identified a number of potential immunomodulator genes of F. hepatica (n = 62), which are firstly described here, could be employed for the development of better strategies (including RNAi) in the battle against both zoonotically and economically important disease, fasciolosis.

Generating a detailed protein profile of Fasciola hepatica during the chronic stage of infection in cattle.

Fasciola hepatica is a trematode helminth causing a damaging disease, fasciolosis, in ruminants and humans. Comprehensive proteomic studies broaden our knowledge of the parasite's protein profile, and provide new insights into the development of more effective strategies to deal with fasciolosis. The objective of this study was to generate a comprehensive profile of F. hepatica proteins expressed during the chronic stage of infection in cattle by building on previous efforts in this area. The approach included an improved sample preparation procedure for surface and internal layers of the parasite, the application of nano-UPLC-ESI-qTOF-MS (nano-ultra-performance LC and ESI quadrupole TOF MS) integrated with different acquisition methods and in silico database search against various protein databases and a transcript database including a new assembly of publically available EST. Of a total of 776 identified proteins, 206 and 332 were specific to the surface and internal layers of the parasite, respectively. Furthermore, 238 proteins were common to both layers, with comparative differences of 172 proteins detected. Specific proteins not previously identified in F. hepatica, but shown to be immunomodulatory or potential drug targets for other parasites, are discussed.

Whole genome sequencing of Turkish genomes reveals functional private alleles and impact of genetic interactions with Europe, Asia and Africa.

BACKGROUND: Turkey is a crossroads of major population movements throughout history and has been a hotspot of cultural interactions. Several studies have investigated the complex population history of Turkey through a limited set of genetic markers. However, to date, there have been no studies to assess the genetic variation at the whole genome level using whole genome sequencing. Here, we present whole genome sequences of 16 Turkish individuals resequenced at high coverage (32×-48×). RESULTS: We show that the genetic variation of the contemporary Turkish population clusters with South European populations, as expected, but also shows signatures of relatively recent contribution from ancestral East Asian populations. In addition, we document a significant enrichment of non-synonymous private alleles, consistent with recent observations in European populations. A number of variants associated with skin color and total cholesterol levels show frequency differentiation between the Turkish populations and European populations. Furthermore, we have analyzed the 17q21.31 inversion polymorphism region (MAPT locus) and found increased allele frequency of 31.25% for H1/H2 inversion polymorphism when compared to European populations that show about 25% of allele frequency. CONCLUSION: This study provides the first map of common genetic variation from 16 western Asian individuals and thus helps fill an important geographical gap in analyzing natural human variation and human migration. Our data will help develop population-specific experimental designs for studies investigating disease associations and demographic history in Turkey.

Robustness of Massively Parallel Sequencing Platforms.

The improvements in high throughput sequencing technologies (HTS) made clinical sequencing projects such as ClinSeq and Genomics England feasible. Although there are significant improvements in accuracy and reproducibility of HTS based analyses, the usability of these types of data for diagnostic and prognostic applications necessitates a near perfect data generation. To assess the usability of a widely used HTS platform for accurate and reproducible clinical applications in terms of robustness, we generated whole genome shotgun (WGS) sequence data from the genomes of two human individuals in two different genome sequencing centers. After analyzing the data to characterize SNPs and indels using the same tools (BWA, SAMtools, and GATK), we observed significant number of discrepancies in the call sets. As expected, the most of the disagreements between the call sets were found within genomic regions containing common repeats and segmental duplications, albeit only a small fraction of the discordant variants were within the exons and other functionally relevant regions such as promoters. We conclude that although HTS platforms are sufficiently powerful for providing data for first-pass clinical tests, the variant predictions still need to be confirmed using orthogonal methods before using in clinical applications.

Activating mutations of STAT5B and STAT3 in lymphomas derived from γδ-T or NK cells.

Lymphomas arising from NK or γδ-T cells are very aggressive diseases and little is known regarding their pathogenesis. Here we report frequent activating mutations of STAT3 and STAT5B in NK/T-cell lymphomas (n=51), γδ-T-cell lymphomas (n=43) and their cell lines (n=9) through next generation and/or Sanger sequencing. STAT5B N642H is particularly frequent in all forms of γδ-T-cell lymphomas. STAT3 and STAT5B mutations are associated with increased phosphorylated protein and a growth advantage to transduced cell lines or normal NK cells. Growth-promoting activity of the mutants can be partially inhibited by a JAK1/2 inhibitor. Molecular modelling and surface plasmon resonance measurements of the N642H mutant indicate a marked increase in binding affinity of the phosphotyrosine-Y699 with the mutant histidine. This is associated with the prolonged persistence of the mutant phosphoSTAT5B and marked increase of binding to target sites. Our findings suggest that JAK-STAT pathway inhibition may represent a therapeutic strategy.

The distinct genetic pattern of ALS in Turkey and novel mutations.

The frequency of amyotrophic lateral sclerosis (ALS) mutations has been extensively investigated in several populations; however, a systematic analysis in Turkish cases has not been reported so far. In this study, we screened 477 ALS patients for mutations, including 116 familial ALS patients from 82 families and 361 sporadic ALS (sALS) cases. Patients were genotyped for C9orf72 (18.3%), SOD1 (12.2%), FUS (5%), TARDBP (3.7%), and UBQLN2 (2.4%) gene mutations, which together account for approximately 40% of familial ALS in Turkey. No SOD1 mutations were detected in sALS patients; however, C9orf72 (3.1%) and UBQLN2 (0.6%) explained 3.7% of sALS in the population. Exome sequencing revealed mutations in OPTN, SPG11, DJ1, PLEKHG5, SYNE1, TRPM7, and SQSTM1 genes, many of them novel. The spectrum of mutations reflect both the distinct genetic background and the heterogeneous nature of the Turkish ALS population.