The use of pharmacogenomics in clinical practice is becoming standard of care. However, due to the complex genetic makeup of pharmacogenes, not all genetic variation is currently accounted for.... Show moreThe use of pharmacogenomics in clinical practice is becoming standard of care. However, due to the complex genetic makeup of pharmacogenes, not all genetic variation is currently accounted for. Here, we show the utility of long-read sequencing to resolve complex pharmacogenes by analyzing a well-characterised sample. This data consists of long reads that were processed to resolve phased haploblocks. 73% of pharmacogenes were fully covered in one phased haploblock, including 9/15 genes that are 100% complex. Variant calling accuracy in the pharmacogenes was high, with 99.8% recall and 100% precision for SNVs and 98.7% precision and 98.0% recall for Indels. For the majority of gene-drug interactions in the DPWG and CPIC guidelines, the associated genes could be fully resolved (62% and 63% respectively). Together, these findings suggest that long-read sequencing data offers promising opportunities in elucidating complex pharmacogenes and haplotype phasing while maintaining accurate variant calling. Show less
Lee, M. van der; Allard, W.G.; Vossen, R.H.A.M.; Baak-Pablo, R.F.; Menafra, R.; Deiman, B.A.L.M.; ... ; Anvar, S.Y. 2021
Pharmacogenomics is a key component of personalized medicine that promises safer and more effective drug treatment by individualizing drug choice and dose based on genetic profiles. In clinical... Show morePharmacogenomics is a key component of personalized medicine that promises safer and more effective drug treatment by individualizing drug choice and dose based on genetic profiles. In clinical practice, genetic biomarkers are used to categorize patients into *-alleles to predict CYP450 enzyme activity and adjust drug dosages accordingly. However, this approach leaves a large part of variability in drug response unexplained. Here, we present a proof-of-concept approach that uses continuous-scale (instead of categorical) assignments to predict enzyme activity. We used full CYP2D6 gene sequences obtained with long-read amplicon-based sequencing and cytochrome P450 (CYP) 2D6-mediated tamoxifen metabolism data from a prospective study of 561 patients with breast cancer to train a neural network. The model explained 79% of interindividual variability in CYP2D6 activity compared to 54% with the conventional *-allele approach, assigned enzyme activities to known alleles with previously reported effects, and predicted the activity of previously uncharacterized combinations of variants. The results were replicated in an independent cohort of tamoxifen-treated patients (model R-2 adjusted = 0.66 versus *-allele R-2 adjusted = 0.35) and a cohort of patients treated with the CYP2D6 substrate venlafaxine (model R2 adjusted = 0.64 versus *-allele R-2 adjusted = 0.55). Human embryonic kidney cells were used to confirm the effect of five genetic variants on metabolism of the CYP2D6 substrate bufuralol in vitro. These results demonstrate the advantage of a continuous scale and a completely phased genotype for prediction of CYP2D6 enzyme activity and could potentially enable more accurate prediction of individual drug response. Show less
During the transition from sexual to asexual reproduction, a suite of reproduction-related sexual traits become superfluous, and may be selected against if costly. Female functional virginity... Show moreDuring the transition from sexual to asexual reproduction, a suite of reproduction-related sexual traits become superfluous, and may be selected against if costly. Female functional virginity refers to asexual females resisting to mate or not fertilizing eggs after mating. These traits appear to be among the first that evolve during transitions from sexual to asexual reproduction. The genetic basis of female functional virginity remains elusive. Previously, we reported that female functional virginity segregates as expected for a single recessive locus in the asexual parasitoid wasp Asobara japonica. Here, we investigate the genetic basis of this trait by quantitative trait loci (QTL) mapping and candidate gene analyses. Consistent with the segregation of phenotypes, we found a single QTL of large effect, spanning over 4.23 Mb and comprising at least 131 protein-coding genes, of which 15 featured sex-biased expression in the related sexual species Asobara tabida. Two of the 15 sex-biased genes were previously identified to differ between related sexual and asexual population/species: CD151 antigen and nuclear pore complex protein Nup50. A third gene, hormone receptor 4, is involved in steroid hormone mediated mating behaviour. Overall, our results are consistent with a single locus, or a cluster of closely linked loci, underlying rapid evolution of female functional virginity in the transition to asexuality. Once this variant, causing rejection to mate, has swept through a population, the flanking region does not get smaller owing to lack of recombination in asexuals. Show less
Jong, A. de; Dirven, R.J.; Boender, J.; Atiq, F.; Anvar, S.Y.; Leebeek, F.W.G.; ... ; Eikenboom, J. 2020
Von Willebrand disease (VWD) is the most common inherited bleeding disorder and is mainly caused by dominant-negative mutations in the multimeric protein von Willebrand factor (VWF). These... Show moreVon Willebrand disease (VWD) is the most common inherited bleeding disorder and is mainly caused by dominant-negative mutations in the multimeric protein von Willebrand factor (VWF). These mutations may either result in quantitative or qualitative defects in VWF. VWF is an endothelial protein that is secreted to the circulation upon endothelial activation. Once secreted, VWF multimers bind platelets and chaperone coagulation factor VIII in the circulation. Treatment of VWD focuses on increasing VWF plasma levels, but production and secretion of mutant VWF remain uninterrupted. Presence of circulating mutant VWF might, however, still affect normal hemostasis or functionalities of VWF beyond hemostasis. We hypothesized that inhibition of the production of mutant VWF improves the function of VWF overall and ameliorates VWD phenotypes. We previously proposed the use of allele-specific small-interfering RNAs (siRNAs) that target frequentVWFsingle nucleotide polymorphisms to inhibit mutantVWF. The aim of this study is to prove the functionality of these allele-specific siRNAs in endothelial colony-forming cells (ECFCs). We isolated ECFCs from a VWD type 2A patient with an intracellular multimerization defect, reduced VWF collagen binding, and a defective processing of proVWF to VWF. After transfection of an allele-specific siRNA that specifically inhibited expression of mutant VWF, we showed amelioration of the laboratory phenotype, with normalization of the VWF collagen binding, improvement in VWF multimers, and enhanced VWF processing. Altogether, we prove that allele-specific inhibition of the production of mutant VWF by siRNAs is a promising therapeutic strategy to improve VWD phenotypes. Show less
Tzovaras, B.G.; Segers, F.H.I.D.; Bicker, A.; Grande, F. dal; Otte, J.; Anvar, S.Y.; ... ; Ebersberger, I. 2020
Lichens are valuable models in symbiosis research and promising sources of biosynthetic genes for biotechnological applications. Most lichenized fungi grow slowly, resist aposymbiotic cultivation,... Show moreLichens are valuable models in symbiosis research and promising sources of biosynthetic genes for biotechnological applications. Most lichenized fungi grow slowly, resist aposymbiotic cultivation, and are poor candidates for experimentation. Obtaining contiguous, high-quality genomes for such symbiotic communities is technically challenging. Here, we present the first assembly of a lichen holo-genome from metagenomic whole-genome shotgun data comprising both PacBio long reads and Illumina short reads. The nuclear genomes of the two primary components of the lichen symbiosis-the fungus Umbilicaria pustulata (33 Mb) and the green alga Trebouxia sp. (53 Mb)-were assembled at contiguities comparable to single-species assemblies. The analysis of the read coverage pattern revealed a relative abundance of fungal to algal nuclei of similar to 20:1. Gap-free, circular sequences for all organellar genomes were obtained. The bacterial community is dominated by Acidobacteriaceae and encompasses strains closely related to bacteria isolated from other lichens. Gene set analyses showed no evidence of horizontal gene transfer from algae or bacteria into the fungal genome. Our data suggest a lineage-specific loss of a putative gibberellin-20-oxidase in the fungus, a gene fusion in the fungal mitochondrion, and a relocation of an algal chloroplast gene to the algal nucleus. Major technical obstacles during reconstruction of the holo-genome were coverage differences among individual genomes surpassing three orders of magnitude. Moreover, we show that GC-rich inverted repeats paired with nonrandom sequencing error in PacBio data can result in missing gene predictions. This likely poses a general problem for genome assemblies based on long reads. Show less
Lee, M. van der; Allard, W.G.; Bollen, S.; Santen, G.W.E.; Ruivenkamp, C.A.L.; Hoffer, M.J.V.; ... ; Swen, J.J. 2019
For similar to 80 drugs, widely recognized pharmacogenetics dosing guidelines are available. However, the use of these guidelines in clinical practice remains limited as only a fraction of patients... Show moreFor similar to 80 drugs, widely recognized pharmacogenetics dosing guidelines are available. However, the use of these guidelines in clinical practice remains limited as only a fraction of patients is subjected to pharmacogenetic screening. We investigated the feasibility of repurposing whole exome sequencing (WES) data for a panel of 42 variants in 11 pharmacogenes to provide a pharmacogenomic profile. Existing diagnostic WES-data from child-parent trios totaling 1,583 individuals were used. Results were successfully extracted for 39 variants. No information could be extracted for three variants, located in CYP2C19, UGT1A1, and CYP3A5, and for CYP2D6 copy number. At least one actionable phenotype was present in 86% of the individuals. Haplotype phasing proved relevant for CYP2B6 assignments as 1.5% of the phenotypes were corrected after phasing. In conclusion, repurposing WES-data can yield meaningful pharmacogenetic profiles for 7 of 11 important pharmacogenes, which can be used to guide drug treatment. Show less
Background: Theory predicts that dependency within host-endosymbiont interactions results in endosymbiont genome size reduction. Unexpectedly, the largest Wolbachia genome was found in the obligate... Show moreBackground: Theory predicts that dependency within host-endosymbiont interactions results in endosymbiont genome size reduction. Unexpectedly, the largest Wolbachia genome was found in the obligate, parthenogenesis-associated wFol. In this study, we investigate possible processes underlying this genome expansion by comparing a re-annotated wFol genome to other Wolbachia genomes. In addition, we also search for candidate genes related to parthenogenesis induction (PI).Results: Within wFol, we found five phage WO regions representing 25.4% of the complete genome, few pseudogenized genes, and an expansion of DNA-repair genes in comparison to other Wolbachia. These signs of genome conservation were mirrored in the wFol host, the springtail F. candida, which also had an expanded DNA-repair gene family and many horizontally transferred genes. Across all Wolbachia genomes, there was a strong correlation between gene numbers of Wolbachia strains and their hosts. In order to identify genes with a potential link to PI, we assembled the genome of an additional PI strain, wLcla. Comparisons between four PI Wolbachia, including wFol and wLcla, and fourteen non-PI Wolbachia yielded a small set of potential candidate genes for further investigation.Conclusions: The strong similarities in genome content of wFol and its host, as well as the correlation between host and Wolbachia gene numbers suggest that there may be some form of convergent evolution between endosymbiont and host genomes. If such convergent evolution would be strong enough to overcome the evolutionary forces causing genome reduction, it would enable expanded genomes within long-term obligate endosymbionts. Show less
Weerts, M.J.A.; Timmermans, E.C.; Stolpe, A. van de; Vossen, R.H.A.M.; Anvar, S.Y.; Foekens, J.A.; ... ; Martens, J.W.M. 2018
BACKGROUND\nThe common carp (Cyprinus carpio) is the oldest, most domesticated and one of the most cultured fish species for food consumption. Besides its economic importance, the common carp is... Show moreBACKGROUND\nThe common carp (Cyprinus carpio) is the oldest, most domesticated and one of the most cultured fish species for food consumption. Besides its economic importance, the common carp is also highly suitable for comparative physiological and disease studies in combination with the animal model zebrafish (Danio rerio). They are genetically closely related but offer complementary benefits for fundamental research, with the large body mass of common carp presenting possibilities for obtaining sufficient cell material for advanced transcriptome and proteome studies.\nRESULTS\nHere we have used 19 different tissues from an F1 hybrid strain of the common carp to perform transcriptome analyses using RNA-Seq. For a subset of the tissues we also have performed deep proteomic studies. As a reference, we updated the European common carp genome assembly using low coverage Pacific Biosciences sequencing to permit high-quality gene annotation. These annotated gene lists were linked to zebrafish homologs, enabling direct comparisons with published datasets. Using clustering, we have identified sets of genes that are potential selective markers for various types of tissues. In addition, we provide a script for a schematic anatomical viewer for visualizing organ-specific expression data.\nCONCLUSIONS\nThe identified transcriptome and proteome data for carp tissues represent a useful resource for further translational studies of tissue-specific markers for this economically important fish species that can lead to new markers for organ development. The similarity to zebrafish expression patterns confirms the value of common carp as a resource for studying tissue-specific expression in cyprinid fish. The availability of the annotated gene set of common carp will enable further research with both applied and fundamental purposes. Show less