BMPR2 as a Novel Predisposition Gene for Hereditary Colorectal Polyposis.

BACKGROUND & AIMS: Colorectal cancer (CRC) is one of the most prevalent tumors worldwide, with incidence quickly increasing (particularly in the context of early-onset cases), despite important prevention efforts, mainly in the form of population-wide screening programs. Although many cases present a clear familial component, the current list of hereditary CRC genes leaves a considerable proportion of the cases unexplained. METHODS: In this work, we used whole-exome sequencing approaches on 19 unrelated patients with unexplained colonic polyposis to identify candidate CRC predisposition genes. The candidate genes were then validated in an additional series of 365 patients. CRISPR-Cas9 models were used to validate BMPR2 as a potential candidate for CRC risk. RESULTS: We found 8 individuals carrying 6 different variants in the BMPR2 gene (approximately 2% of our cohort of patients with unexplained colonic polyposis). CRISPR-Cas9 models of 3 of these variants showed that the p.(Asn442Thrfs∗32) truncating variant completely abrogated BMP pathway function in a similar way to the BMPR2 knockout. Missense variants p.(Asn565Ser), p.(Ser967Pro) had varying effects on cell proliferation levels, with the former impairing cell control inhibition via noncanonical pathways. CONCLUSIONS: Collectively, these results support loss-of-function BMPR2 variants as candidates to be involved in CRC germline predisposition.

Design, synthesis and validation of a new Crimped Head-Piece for DNA-Encoded libraries generation.

Codification of DNA Encoded Libraries (DELs) is critical for successful ligand identification of molecules that bind a protein of interest (POI). There are different encoding strategies that permit, for instance, the customization of a DEL for testing single or dual pharmacophores (single strand DNA) or for producing and screening large diversity libraries of small molecules (double strand DNA). Both approaches challenges, either from the synthetic and encoding point of view, or from the selection methodology to be utilized for the screening. The Head-Piece contains the DNA sequence that is attached to a chemical compound, allowing the encoding of each molecule with a unique DNA tag. Designing the Head-Piece for a DNA-encoded library involves careful consideration of several key aspects including DNA barcode identity, sequence length and attachment chemistry. Here we describe a double stranded DNA versatile Head-Piece that can be used for the generation of single or dual pharmacophore libraries, but also shows other advanced DEL functionalities, stability and enlarged encoding capacity.

CSVS, a crowdsourcing database of the Spanish population genetic variability.

The knowledge of the genetic variability of the local population is of utmost importance in personalized medicine and has been revealed as a critical factor for the discovery of new disease variants. Here, we present the Collaborative Spanish Variability Server (CSVS), which currently contains more than 2000 genomes and exomes of unrelated Spanish individuals. This database has been generated in a collaborative crowdsourcing effort collecting sequencing data produced by local genomic projects and for other purposes. Sequences have been grouped by ICD10 upper categories. A web interface allows querying the database removing one or more ICD10 categories. In this way, aggregated counts of allele frequencies of the pseudo-control Spanish population can be obtained for diseases belonging to the category removed. Interestingly, in addition to pseudo-control studies, some population studies can be made, as, for example, prevalence of pharmacogenomic variants, etc. In addition, this genomic data has been used to define the first Spanish Genome Reference Panel (SGRP1.0) for imputation. This is the first local repository of variability entirely produced by a crowdsourcing effort and constitutes an example for future initiatives to characterize local variability worldwide. CSVS is also part of the GA4GH Beacon network. CSVS can be accessed at: http://csvs.babelomics.org/.

Normalization of DNA encoded library affinity selection results driven by high throughput sequencing and HPLC purification.

The output of an affinity selection screening results in a huge amount of valuable data that, after conducting the appropriate analysis, lead to the correct identification of the compounds enriched in the target of interest. The approach chosen to perform these analyses has become a key step in the development of a successful DNA Encoded Library platform. In this paper, we describe the combination of High Performance Liquid Chromatography purification during the library production with the Next Generation Sequencing analysis of the libraries to assess the yield of the chemical reactions prior to the affinity selection. This process allows us, apart from achieving higher quality libraries, to enable a normalization analysis of the affinity selection output, thus minimizing the bias induced by the chemical yield of each reaction as a misleading factor within the analysis and subsequent compound short-listing for off-DNA synthesis.

Identification of putative second genetic hits in schizophrenia carriers of high-risk copy number variants and resequencing in additional samples.

Copy number variants (CNVs) conferring risk of schizophrenia present incomplete penetrance, suggesting the existence of second genetic hits. Identification of second hits may help to find genes with rare variants of susceptibility to schizophrenia. The aim of this work was to search for second hits of moderate/high risk in schizophrenia carriers of risk CNVs and resequencing of the relevant genes in additional samples. To this end, ten patients with risk CNVs at cytobands 15q11.2, 15q11.2-13.1, 16p11.2, or 16p13.11, were subjected to whole-exome sequencing. Rare single nucleotide variants, defined as those absent from main public databases, were classified according to bioinformatic prediction of pathogenicity by CADD scores. The average number of rare predicted pathogenic variants per sample was 13.6 (SD 2.01). Two genes, BFAR and SYNJ1, presented rare predicted pathogenic variants in more than one sample. Follow-up resequencing of these genes in 432 additional cases and 432 controls identified a significant excess of rare predicted pathogenic variants in case samples at SYNJ1. Taking into account its function in clathrin-mediated synaptic vesicle endocytosis at presynaptic terminals, our results suggest an impairment of this process in schizophrenia.

Postmortem genetic testing should be recommended in sudden cardiac death cases due to thoracic aortic dissection.

BACKGROUND: Acute thoracic aortic dissections and ruptures, the main life-threatening complications of the corresponding aneurysms, are an important cause of sudden cardiac death. Despite the usefulness of the molecular diagnosis of these conditions in the clinical setting, the corresponding forensic field remains largely unexplored. The main goal of this study was to explore and validate a new massive parallel sequencing candidate gene​ assay as a diagnostic tool for acute thoracic aortic dissection autopsy cases. MATERIALS AND METHODS: Massive parallel sequencing of 22 thoracic aortic disease candidate genes performed in 17 cases of thoracic aortic dissection using AmpliSeq and Ion Proton technologies. Genetic variants were filtered by location, type, and frequency at the Exome Aggregation Consortium and an internal database and further classified based on the American College of Medical Genetics and Genomics (ACMG) recommendations published in 2015. All prioritized results were confirmed by traditional sequencing. RESULTS: From the total of 10 potentially pathogenic genetic variants identified in 7 out of the 17 initial samples, 2 of them were further classified as pathogenic, 2 as likely pathogenic, 1 as possibly benign, and the remaining 5 as variants of uncertain significance, reaching a molecular autopsy yield of 23%, approximately. CONCLUSIONS: This massive parallel sequencing candidate gene approach proved useful for the molecular autopsy of aortic dissection sudden cardiac death cases and should therefore be progressively incorporated into the forensic field, being especially beneficial for the anticipated diagnosis and risk stratification of any other family member at risk of developing the same condition.

BigBWA: approaching the Burrows-Wheeler aligner to Big Data technologies.

UNLABELLED: BigBWA is a new tool that uses the Big Data technology Hadoop to boost the performance of the Burrows-Wheeler aligner (BWA). Important reductions in the execution times were observed when using this tool. In addition, BigBWA is fault tolerant and it does not require any modification of the original BWA source code. AVAILABILITY AND IMPLEMENTATION: BigBWA is available at the project GitHub repository: https://github.com/citiususc/BigBWA.

A novel stop mutation in the vascular endothelial growth factor-C gene (VEGFC) results in Milroy-like disease.

BACKGROUND: Milroy and Milroy-like disease are rare disorders characterised by congenital lymphoedema caused by dysfunctional lymphatic vessel formation. Loss of extracellular response mediated by vascular endothelial growth factor receptor 3 (VEGFR-3) is associated with Milroy disease, and VEGFR-3 gene is mutated in around 70% of the cases diagnosed. The only genetic alteration known to be associated with Milroy-like disease was recently identified in a family with a frameshift mutation in vascular endothelial growth factor C (VEGFC) gene, which encodes a VEGFR3 ligand. METHODS AND RESULTS: We report a newborn patient with an external phenotype consistent with Milroy disease and a truncating mutation (p.R210X) in the VEGFC gene detected by exome sequence analysis. Subsequent analysis, by lymphoscintigraphic scan, performed for research purposes, allowed us to correct the diagnosis, confirming patient's disease as Milroy-like. The mutation segregates with the phenotype in the family according to a dominant model with full penetrance. CONCLUSIONS: The clinical presentation, similar to Milroy disease, indicates an overlapping of the external phenotype of both diseases, suggesting that genetic analysis of VEGFC would be useful in diagnosing patients that present with Milroy features but have no mutation in VEGFR-3. Establishing a well-defined genetic pattern would help with differential diagnosis.

Adaptive selection of an incretin gene in Eurasian populations.

Diversities in human physiology have been partially shaped by adaptation to natural environments and changing cultures. Recent genomic analyses have revealed single nucleotide polymorphisms (SNPs) that are associated with adaptations in immune responses, obvious changes in human body forms, or adaptations to extreme climates in select human populations. Here, we report that the human GIP locus was differentially selected among human populations based on the analysis of a nonsynonymous SNP (rs2291725). Comparative and functional analyses showed that the human GIP gene encodes a cryptic glucose-dependent insulinotropic polypeptide (GIP) isoform (GIP55S or GIP55G) that encompasses the SNP and is resistant to serum degradation relative to the known mature GIP peptide. Importantly, we found that GIP55G, which is encoded by the derived allele, exhibits a higher bioactivity compared with GIP55S, which is derived from the ancestral allele. Haplotype structure analysis suggests that the derived allele at rs2291725 arose to dominance in East Asians ∼8100 yr ago due to positive selection. The combined results suggested that rs2291725 represents a functional mutation and may contribute to the population genetics observation. Given that GIP signaling plays a critical role in homeostasis regulation at both the enteroinsular and enteroadipocyte axes, our study highlights the importance of understanding adaptations in energy-balance regulation in the face of the emerging diabetes and obesity epidemics.

Next generation sequencing challenges in the analysis of cardiac sudden death due to arrhythmogenic disorders.

Inherited arrhythmogenic disorders is a relatively common cause of cardiac sudden death in young people. Diagnosis has been difficult so far due to the genetic heterogeneity of the disease. Next generation sequencing (NGS) is offering a new scenario for diagnosis. The purpose of our study was to validate NGS for the analysis of twenty-eight genes known to be associated with inherited arrhythmogenic disorders and therefore with sudden cardiac death. SureSelect hybridization was used to enrich DNA from 53 samples, prior to be sequenced with the SOLID™ System of Life Technologies. Depth of coverage, consistency of coverage across samples, and location of variants identified were assessed. All the samples showed a depth of coverage over 200×, except one of them discarded because of its coverage below 30×. Average percentage of target bp covered at least 20× was 96.45%. In the remaining samples, following a prioritization process 46 possible variants in 31 samples were found, of which 45 were confirmed by Sanger sequencing. After filtering variants according to their minor allele frequency in the Exome Sequencing Project 27 putative pathogenic variants in 20 samples remained. With the use of in silico tools, 13 variants in 11 samples were classified as likely pathogenic. In conclusion, NGS allowed us to accurately detect arrhythmogenic disease causing mutations in a fast and cost-efficient manner that is suitable for daily clinical and forensic practice of genetic testing of this type of disorders.

Clinical and Genetic Analysis of Patients With TK2 Deficiency.

Objectives: Thymidine kinase 2 deficiency (TK2d) is a rare autosomal recessive disorder that stems from a perturbation of the mitochondrial DNA maintenance. Nucleoside treatment has recently shown promise as a disease-modifying therapy. TK2d was initially associated with rapidly progressive fatal myopathy in children featuring mitochondrial DNA depletion. Subsequently, less severe variants of the disease were described, with onset of symptoms during adolescence or adulthood and associated with the presence of multiple mtDNA deletions. These less severe phenotypes have been reported in only 15% of the approximately 120 patients described worldwide. However, some reports suggest that these juvenile and adult-onset presentations may be more common. The objective of this study was to describe the clinical phenotype in a sample of patients from Spain. Methods: This study includes 53 patients harboring biallelic TK2 pathogenic variants, compiling data retrospectively from 7 Spanish centers. We analyzed allele frequency, investigated the most recent common ancestor of core haplotypes, and used the Runs of Homozygosity approach to investigate variant coalescence. Results: Symptom onset distribution revealed that 32 patients (60%) experienced symptoms beyond 12 years of age. Approximately 30% of patients died of respiratory insufficiency, while 56% of surviving patients needed mechanical ventilation. Genetic analysis identified 16 distinct variants in TK2. Two variants, p.Lys202del and p.Thr108Met, exhibited significantly higher prevalence in the Spanish population than that reported in gnomAD database (86-fold and 13-fold, respectively). These variants are estimated to have originated approximately 16.8 generations ago for p.Thr108Met and 95.2 generations ago for p.Lys202del within the Spanish population, with the increase in frequency attributed to various forms of inbreeding. In late-onset cases, 46.9% carried the p.Lys202del variant. Discussion: The higher frequency of TK2d in Spain can be partially attributed to the increased prevalence of 2 variants and consanguinity. Notably, in 60% of the cohort, the disease was late-onset, emphasizing the potential underdiagnosis of this subgroup of patients in other regions. Raising awareness of this potentially treatable disorder is of utmost importance because early interventions can significantly affect the quality of life and survival of affected individuals.

Detection of the Copy Number Variants of Genes in Patients with Familial Cardiac Diseases by Massively Parallel Sequencing.

INTRODUCTION: The implication of copy number variations in familial heart disease is known, although in-depth knowledge is lacking; hence, more studies are needed to further our understanding. Massively parallel sequencing, thanks to its recent surge in use, is emerging as a valid tool for the detection of this type of variant, through the use of appropriate software. METHODS: We conducted a study with 182 patients diagnosed with mendelian cardiovascular diseases who underwent sequencing using a cardiac gene panel and then a specific calling process for copy number variations (CNVs) with ExomeDepth software, which provides us with a Bayes factor (BF), a score of the probability that a CNV detected is true. RESULTS: After a rigorous CNV prioritization process, we confirmed the variants obtained by MLPA or SNP-based array, finding three real CNVs in five individuals in the MYH11, FBN1 and PDMI7 genes. CONCLUSION: The confirmed CNVs present in all cases BF values > 60, thus establishing a threshold to consider real CNVs in the calling process carried out by ExomeDepth on our gene panel.

Genome-wide identification of hypoxia-inducible factor binding sites and target genes by a probabilistic model integrating transcription-profiling data and in silico binding site prediction.

The transcriptional response driven by Hypoxia-inducible factor (HIF) is central to the adaptation to oxygen restriction. Hence, the complete identification of HIF targets is essential for understanding the cellular responses to hypoxia. Herein we describe a computational strategy based on the combination of phylogenetic footprinting and transcription profiling meta-analysis for the identification of HIF-target genes. Comparison of the resulting candidates with published HIF1a genome-wide chromatin immunoprecipitation indicates a high sensitivity (78%) and specificity (97.8%). To validate our strategy, we performed HIF1a chromatin immunoprecipitation on a set of putative targets. Our results confirm the robustness of the computational strategy in predicting HIF-binding sites and reveal several novel HIF targets, including RE1-silencing transcription factor co-repressor (RCOR2). In addition, mapping of described polymorphisms to the predicted HIF-binding sites identified several single-nucleotide polymorphisms (SNPs) that could alter HIF binding. As a proof of principle, we demonstrate that SNP rs17004038, mapping to a functional hypoxia response element in the macrophage migration inhibitory factor (MIF) locus, prevents induction of this gene by hypoxia. Altogether, our results show that the proposed strategy is a powerful tool for the identification of HIF direct targets that expands our knowledge of the cellular adaptation to hypoxia and provides cues on the inter-individual variation in this response.

Trio-based exome sequencing reveals a high rate of the de novo variants in intellectual disability.

Intellectual disability (ID), a neurodevelopmental disorder affecting 1-3% of the general population, is characterized by limitations in both intellectual function and adaptive skills. The high number of conditions associated with ID underlines its heterogeneous origin and reveals the difficulty of obtaining a rapid and accurate genetic diagnosis. However, the Next Generation Sequencing, and the whole exome sequencing (WES) in particular, has boosted the diagnosis rate associated with ID. In this study, WES performed on 244 trios of patients clinically diagnosed with isolated or syndromic ID and their respective unaffected parents has allowed the identification of the underlying genetic basis of ID in 64 patients, yielding a diagnosis rate of 25.2%. Our results suggest that trio-based WES facilitates ID's genetic diagnosis, particularly in patients who have been extensively waiting for a definitive molecular diagnosis. Moreover, genotypic information from parents provided by trio-based WES enabled the detection of a high percentage (61.5%) of de novo variants inside our cohort. Establishing a quick genetic diagnosis of ID would allow early intervention and better clinical management, thus improving the quality of life of these patients and their families.

ENGINES: exploring single nucleotide variation in entire human genomes.

BACKGROUND: Next generation ultra-sequencing technologies are starting to produce extensive quantities of data from entire human genome or exome sequences, and therefore new software is needed to present and analyse this vast amount of information. The 1000 Genomes project has recently released raw data for 629 complete genomes representing several human populations through their Phase I interim analysis and, although there are certain public tools available that allow exploration of these genomes, to date there is no tool that permits comprehensive population analysis of the variation catalogued by such data. DESCRIPTION: We have developed a genetic variant site explorer able to retrieve data for Single Nucleotide Variation (SNVs), population by population, from entire genomes without compromising future scalability and agility. ENGINES (ENtire Genome INterface for Exploring SNVs) uses data from the 1000 Genomes Phase I to demonstrate its capacity to handle large amounts of genetic variation (>7.3 billion genotypes and 28 million SNVs), as well as deriving summary statistics of interest for medical and population genetics applications. The whole dataset is pre-processed and summarized into a data mart accessible through a web interface. The query system allows the combination and comparison of each available population sample, while searching by rs-number list, chromosome region, or genes of interest. Frequency and FST filters are available to further refine queries, while results can be visually compared with other large-scale Single Nucleotide Polymorphism (SNP) repositories such as HapMap or Perlegen. CONCLUSIONS: ENGINES is capable of accessing large-scale variation data repositories in a fast and comprehensive manner. It allows quick browsing of whole genome variation, while providing statistical information for each variant site such as allele frequency, heterozygosity or FST values for genetic differentiation. Access to the data mart generating scripts and to the web interface is granted from http://spsmart.cesga.es/engines.php.

Call for participation in the neurogenetics consortium within the Human Variome Project.

The rate of DNA variation discovery has accelerated the need to collate, store and interpret the data in a standardised coherent way and is becoming a critical step in maximising the impact of discovery on the understanding and treatment of human disease. This particularly applies to the field of neurology as neurological function is impaired in many human disorders. Furthermore, the field of neurogenetics has been proven to show remarkably complex genotype-to-phenotype relationships. To facilitate the collection of DNA sequence variation pertaining to neurogenetic disorders, we have initiated the "Neurogenetics Consortium" under the umbrella of the Human Variome Project. The Consortium's founding group consisted of basic researchers, clinicians, informaticians and database creators. This report outlines the strategic aims established at the preliminary meetings of the Neurogenetics Consortium and calls for the involvement of the wider neurogenetic community in enabling the development of this important resource.

The impact of modern migrations on present-day multi-ethnic Argentina as recorded on the mitochondrial DNA genome.

BACKGROUND: The genetic background of Argentineans is a mosaic of different continental ancestries. From colonial to present times, the genetic contribution of Europeans and sub-Saharan Africans has superposed to or replaced the indigenous genetic 'stratum'. A sample of 384 individuals representing different Argentinean provinces was collected and genotyped for the first and the second mitochondrial DNA (mtDNA) hypervariable regions, and selectively genotyped for mtDNA SNPs. This data was analyzed together with additional 440 profiles from rural and urban populations plus 304 from Native American Argentineans, all available from the literature. A worldwide database was used for phylogeographic inferences, inter-population comparisons, and admixture analysis. Samples identified as belonging to hg (hg) H2a5 were sequenced for the entire mtDNA genome. RESULTS: Phylogenetic and admixture analyses indicate that only half of the Native American component in urban Argentineans might be attributed to the legacy of extinct ancestral Argentineans and that the Spanish genetic contribution is slightly higher than the Italian one. Entire H2a5 genomes linked these Argentinean mtDNAs to the Basque Country and improved the phylogeny of this Basque autochthonous clade. The fingerprint of African slaves in urban Argentinean mtDNAs was low and it can be phylogeographically attributed predominantly to western African. The European component is significantly more prevalent in the Buenos Aires province, the main gate of entrance for Atlantic immigration to Argentina, while the Native American component is larger in North and South Argentina. AMOVA, Principal Component Analysis and hgs/haplotype patterns in Argentina revealed an important level of genetic sub-structure in the country. CONCLUSIONS: Studies aimed to compare mtDNA frequency profiles from different Argentinean geographical regions (e.g., forensic and case-control studies) should take into account the important genetic heterogeneity of the country in order to prevent false positive claims of association in disease studies or inadequate evaluation of forensic evidence.

Longitudinal analysis on parasite diversity in honeybee colonies: new taxa, high frequency of mixed infections and seasonal patterns of variation.

To evaluate the influence that parasites have on the losses of Apis mellifera it is essential to monitor their presence in the colonies over time. Here we analysed the occurrence of nosematids, trypanosomatids and neogregarines in five homogeneous colonies for up to 21 months until they collapsed. The study, which combined the use of several molecular markers with the application of a massive parallel sequencing technology, provided valuable insights into the epidemiology of these parasites: (I) it enabled the detection of parasite species rarely reported in honeybees (Nosema thomsoni, Crithidia bombi, Crithidia acanthocephali) and the identification of two novel taxa; (II) it revealed the existence of a high rate of co-infections (80% of the samples harboured more than one parasite species); (III) it uncovered an identical pattern of seasonal variation for nosematids and trypanosomatids, that was different from that of neogregarines; (IV) it showed that there were no significant differences in the fraction of positive samples, nor in the levels of species diversity, between interior and exterior bees; and (V) it unveiled that the variation in the number of parasite species was not directly linked with the failure of the colonies.

Broadening the Applicability of a Custom Multi-Platform Panel of Microhaplotypes: Bio-Geographical Ancestry Inference and Expanded Reference Data.

The development of microhaplotype (MH) panels for massively parallel sequencing (MPS) platforms is gaining increasing relevance for forensic analysis. Here, we expand the applicability of a 102 autosomal and 11 X-chromosome panel of MHs, previously validated with both MiSeq and Ion S5 MPS platforms and designed for identification purposes. We have broadened reference population data for identification purposes, including data from 240 HGDP-CEPH individuals of native populations from North Africa, the Middle East, Oceania and America. Using the enhanced population data, the panel was evaluated as a marker set for bio-geographical ancestry (BGA) inference, providing a clear differentiation of the five main continental groups of Africa, Europe, East Asia, Native America, and Oceania. An informative degree of differentiation was also achieved for the population variation encompassing North Africa, Middle East, Europe, South Asia, and East Asia. In addition, we explored the potential for individual BGA inference from simple mixed DNA, by simulation of mixed profiles followed by deconvolution of mixture components.

Viability of in-house datamarting approaches for population genetics analysis of SNP genotypes.

BACKGROUND: Databases containing very large amounts of SNP (Single Nucleotide Polymorphism) data are now freely available for researchers interested in medical and/or population genetics applications. While many of these SNP repositories have implemented data retrieval tools for general-purpose mining, these alone cannot cover the broad spectrum of needs of most medical and population genetics studies. RESULTS: To address this limitation, we have built in-house customized data marts from the raw data provided by the largest public databases. In particular, for population genetics analysis based on genotypes we have built a set of data processing scripts that deal with raw data coming from the major SNP variation databases (e.g. HapMap, Perlegen), stripping them into single genotypes and then grouping them into populations, then merged with additional complementary descriptive information extracted from dbSNP. This allows not only in-house standardization and normalization of the genotyping data retrieved from different repositories, but also the calculation of statistical indices from simple allele frequency estimates to more elaborate genetic differentiation tests within populations, together with the ability to combine population samples from different databases. CONCLUSION: The present study demonstrates the viability of implementing scripts for handling extensive datasets of SNP genotypes with low computational costs, dealing with certain complex issues that arise from the divergent nature and configuration of the most popular SNP repositories. The information contained in these databases can also be enriched with additional information obtained from other complementary databases, in order to build a dedicated data mart. Updating the data structure is straightforward, as well as permitting easy implementation of new external data and the computation of supplementary statistical indices of interest.