Small Pedigrees Research Articles

Versatile Identification of Copy Number Variants with Canvas.

Versatile and efficient variant calling tools are needed to analyze large-scale sequencing datasets. In particular, identification of copy number changes remains a challenging task due to their complexity, susceptibility to sequencing biases, variation in coverage data and dependence on genome-wide sample properties, such as tumor polyploidy, polyclonality in cancer samples, or frequency of de novo variation in germline genomes of pedigrees. The frequent need of core sequencing facilities to process samples from both normal and tumor sources favors multipurpose variant calling tools with functionality to process these diverse sets within a single software framework. This not only simplifies the overall bioinformatics workflow but also streamlines maintenance by shortening the software update cycle and requiring only limited staff training. Here we introduce Canvas, a tool for identification of copy number changes from diverse sequencing experiments including whole-genome matched tumor-normal, small pedigree, and single-sample normal resequencing, as well as whole-exome matched and unmatched tumor-normal studies. In addition to variant calling, Canvas infers genome-wide parameters such as cancer ploidy, purity, and heterogeneity. It provides fast and easy-to-run workflows that can scale to thousands of samples and can be easily incorporated into variant calling pipelines.

A pipeline combining multiple strategies for prioritizing heterozygous variants for the identification of candidate genes in exome datasets

BackgroundThe identification of disease-causing variants in autosomal dominant diseases using exome-sequencing data remains a difficult task in small pedigrees. We combined several strategies to improve filtering and prioritizing of heterozygous variants using exome-sequencing datasets in familial Meniere disease: an in-house Pathogenic Variant (PAVAR) score, the Variant Annotation Analysis and Search Tool (VAAST-Phevor), Exomiser-v2, CADD, and FATHMM. We also validated the method by a benchmarking procedure including causal mutations in synthetic exome datasets.ResultsPAVAR and VAAST were able to select the same sets of candidate variants independently of the studied disease. In contrast, Exomiser V2 and VAAST-Phevor had a variable correlation depending on the phenotypic information available for the disease on each family. Nevertheless, all the selected diseases ranked a limited number of concordant variants in the top 10 ranking, using the three systems or other combined algorithm such as CADD or FATHMM.Benchmarking analyses confirmed that the combination of systems with different approaches improves the prediction of candidate variants compared with the use of a single method. The overall efficiency of combined tools ranges between 68 and 71% in the top 10 ranked variants.ConclusionsOur pipeline prioritizes a short list of heterozygous variants in exome datasets based on the top 10 concordant variants combining multiple systems.

Copy number variants in the sheep genome detected using multiple approaches.

BackgroundCopy number variants (CNVs) are a type of polymorphism found to underlie phenotypic variation, both in humans and livestock. Most surveys of CNV in livestock have been conducted in the cattle genome, and often utilise only a single approach for the detection of copy number differences. Here we performed a study of CNV in sheep, using multiple methods to identify and characterise copy number changes. Comprehensive information from small pedigrees (trios) was collected using multiple platforms (array CGH, SNP chip and whole genome sequence data), with these data then analysed via multiple approaches to identify and verify CNVs.ResultsIn total, 3,488 autosomal CNV regions (CNVRs) were identified in this study, which substantially builds on an initial survey of the sheep genome that identified 135 CNVRs. The average length of the identified CNVRs was 19 kb (range of 1 kb to 3.6 Mb), with shorter CNVRs being more frequent than longer CNVRs. The total length of all CNVRs was 67.6Mbps, which equates to 2.7 % of the sheep autosomes. For individuals this value ranged from 0.24 to 0.55 %, and the majority of CNVRs were identified in single animals. Rather than being uniformly distributed throughout the genome, CNVRs tended to be clustered. Application of three independent approaches for CNVR detection facilitated a comparison of validation rates. CNVs identified on the Roche-NimbleGen 2.1M CGH array generally had low validation rates with lower density arrays, while whole genome sequence data had the highest validation rate (>60 %).ConclusionsThis study represents the first comprehensive survey of the distribution, prevalence and characteristics of CNVR in sheep. Multiple approaches were used to detect CNV regions and it appears that the best method for verifying CNVR on a large scale involves using a combination of detection methodologies. The characteristics of the 3,488 autosomal CNV regions identified in this study are comparable to other CNV regions reported in the literature and provide a valuable and sizeable addition to the small subset of published sheep CNVs.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2754-7) contains supplementary material, which is available to authorized users.

The genetics of amyotrophic lateral sclerosis: current insights.

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disorder that results in loss of the upper and lower motor neurons from motor cortex, brainstem, and spinal cord. While the majority of cases are sporadic, approximately 10% show familial inheritance. ALS is usually inherited in an autosomal dominant manner, although autosomal recessive and X-linked inheritance do occur. To date, 24 of the genes at 26 loci have been identified; these include loci linked to ALS and to frontotemporal dementia-ALS, where family pedigrees contain individuals with frontotemporal dementia with/without ALS. The most commonly established genetic causes of familial ALS (FALS) to date are the presence of a hexanucleotide repeat expansion in the C9ORF72 gene (39.3% FALS) and mutation of SOD1, TARDBP, and FUS, with frequencies of 12%–23.5%, 5%, and 4.1%, respectively. However, with the increasing use of next-generation sequencing of small family pedigrees, this has led to an increasing number of genes being associated with ALS. This review provides a comprehensive review on the genetics of ALS and an update of the pathogenic mechanisms associated with these genes. Commonly implicated pathways have been established, including RNA processing, the protein degradation pathways of autophagy and ubiquitin–proteasome system, as well as protein trafficking and cytoskeletal function. Elucidating the role genetics plays in both FALS and sporadic ALS is essential for understanding the subsequent cellular dysregulation that leads to motor neuron loss, in order to develop future effective therapeutic strategies.

A Pipeline for Classifying Relationships Using Dense SNP/SNV Data and Putative Pedigree Information.

When genome-wide association studies (GWAS) or sequencing studies are performed on family-based datasets, the genotype data can be used to check the structure of putative pedigrees. Even in datasets of putatively unrelated people, close relationships can often be detected using dense single-nucleotide polymorphism/variant (SNP/SNV) data. A number of methods for finding relationships using dense genetic data exist, but they all have certain limitations, including that they typically use average genetic sharing, which is only a subset of the available information. Here, we present a set of approaches for classifying relationships in GWAS datasets or large-scale sequencing datasets. We first propose an empirical method for detecting identity by descent segments in close relative pairs using un-phased dense SNP data and demonstrate how that information can assist in building a relationship classifier. We then develop a strategy to take advantage of putative pedigree information to enhance classification accuracy. Our methods are tested and illustrated with two datasets from two distinct populations. Finally, we propose classification pipelines for checking and identifying relationships in datasets containing a large number of small pedigrees.

Short communication: The effect of genotyping cows to improve the reliability of genomic predictions for selection candidates

In this study we investigate the potential of enlarging the reference population for genomic prediction in dairy cattle by routinely genotyping a random sample of the first-crop daughters of every AI bull in the breeding program. We analyzed small nuclear pedigrees, each consisting of a genotyped selection candidate and 3 generations of genotyped male ancestors. Genotypes were taken from the genomic routine evaluation of Fleckvieh cattle in Germany and Austria. The phenotypic information of a daughter of any one male in each of these pedigrees was either considered to be part of the daughter yield deviation of the corresponding sire, or was assumed to be an individually observed genotyped daughter of this sire. Daughter genotypes in this case were simulated from phased haplotypes of their sires and random maternal gametes drawn from a haplotype library. We measured the gain from genotyping daughters as the increase in model-based theoretical reliability of the genomic prediction for a putative selection candidate. We expressed the improvements as a marginal increase, corresponding to an increase in reliability at a reliability baseline level of zero, to simplify comparisons. Results were encouraging with 2 to 40% of marginal reliability increase for selection candidates depending on the assumed heritability of the trait and the number of daughters modeled to be genotyped in the design.

Next generation sequencing in a large cohort of patients presenting with neuromuscular disease before or at birth.

BackgroundFetal akinesia/hypokinesia, arthrogryposis and severe congenital myopathies are heterogeneous conditions usually presenting before or at birth. Although numerous causative genes have been identified for each of these disease groups, in many cases a specific genetic diagnosis remains elusive. Due to the emergence of next generation sequencing, virtually the entire coding region of an individual’s DNA can now be analysed through “whole” exome sequencing, enabling almost all known and novel disease genes to be investigated for disorders such as these.MethodsGenomic DNA samples from 45 patients with fetal akinesia/hypokinesia, arthrogryposis or severe congenital myopathies from 38 unrelated families were subjected to next generation sequencing. Clinical features and diagnoses for each patient were supplied by referring clinicians. Genomic DNA was used for either whole exome sequencing or a custom-designed neuromuscular sub-exomic supercapture array containing 277 genes responsible for various neuromuscular diseases. Candidate disease-causing variants were investigated and confirmed using Sanger sequencing. Some of the cases within this cohort study have been published previously as separate studies.ResultsA conclusive genetic diagnosis was achieved for 18 of the 38 families. Within this cohort, mutations were found in eight previously known neuromuscular disease genes (CHRND, CHNRG, ECEL1, GBE1, MTM1, MYH3, NEB and RYR1) and four novel neuromuscular disease genes were identified and have been published as separate reports (GPR126, KLHL40, KLHL41 and SPEG). In addition, novel mutations were identified in CHRND, KLHL40, NEB and RYR1. Autosomal dominant, autosomal recessive, X-linked, and de novo modes of inheritance were observed.ConclusionsBy using next generation sequencing on a cohort of 38 unrelated families with fetal akinesia/hypokinesia, arthrogryposis, or severe congenital myopathy we therefore obtained a genetic diagnosis for 47 % of families. This study highlights the power and capacity of next generation sequencing (i) to determine the aetiology of genetically heterogeneous neuromuscular diseases, (ii) to identify novel disease genes in small pedigrees or isolated cases and (iii) to refine the interplay between genetic diagnosis and clinical evaluation and management.Electronic supplementary materialThe online version of this article (doi:10.1186/s13023-015-0364-0) contains supplementary material, which is available to authorized users.

PBAP: a pipeline for file processing and quality control of pedigree data with dense genetic markers.

Huge genetic datasets with dense marker panels are now common. With the availability of sequence data and recognition of importance of rare variants, smaller studies based on pedigrees are again also common. Pedigree-based samples often start with a dense marker panel, a subset of which may be used for linkage analysis to reduce computational burden and to limit linkage disequilibrium between single-nucleotide polymorphisms (SNPs). Programs attempting to select markers for linkage panels exist but lack flexibility. We developed a pedigree-based analysis pipeline (PBAP) suite of programs geared towards SNPs and sequence data. PBAP performs quality control, marker selection and file preparation. PBAP sets up files for MORGAN, which can handle analyses for small and large pedigrees, typically human, and results can be used with other programs and for downstream analyses. We evaluate and illustrate its features with two real datasets. PBAP scripts may be downloaded from http://faculty.washington.edu/wijsman/software.shtml. wijsman@uw.edu. Supplementary data are available at Bioinformatics online.

A rare variant in MCF2L identified using exclusion linkage in a pedigree with premature atherosclerosis

Cardiovascular disease (CVD) is a major cause of death in Western societies. CVD risk is largely genetically determined. The molecular pathology is, however, not elucidated in a large number of families suffering from CVD. We applied exclusion linkage analysis and next-generation sequencing to elucidate the molecular defect underlying premature CVD in a small pedigree, comprising two generations of which six members suffered from premature CVD. A total of three variants showed co-segregation with the disease status in the family. Two of these variants were excluded from further analysis based on the prevalence in replication cohorts, whereas a non-synonymous variant in MCF.2 Cell Line Derived Transforming Sequence-like protein (MCF2L, c.2066A>G; p.(Asp689Gly); NM_001112732.1), located in the DH domain, was only present in the studied family. MCF2L is a guanine exchange factor that potentially links pathways that signal through Rac1 and RhoA. Indeed, in HeLa cells, MCF2L689Gly failed to activate Rac1 as well as RhoA, resulting in impaired stress fiber formation. Moreover, MCF2L protein was found in human atherosclerotic lesions but not in healthy tissue segments. In conclusion, a rare functional variant in MCF2L, leading to impaired DH function, was identified in a small pedigree with premature CVD. The presence of MCF2L in human atherosclerotic plaque specimen lends further support to its potential role in atherosclerosis.

Examining KIT, EDNRB, MITF, PMEL17, MLANA Genes to Determine the Genetic Basis of Equine Deafness

There appears to be a correlation between congenital deafness and lack of pigmentation in a small pedigree of Spanish Mustangs. At this time, no hereditary explanation has been presented for this pattern as it appears in horses. Based on the previously determined correlation of melanocyte function and/or migration to deafness in humans, we have conducted a candidate gene study to attempt to link the white coat pigmentation, blue eyes and cochleosaccular deafness within this pedigree. Candidate genes were chosen on the basis of correlation to similar disease states in humans, dogs and cats, including KIT, EDNRB, MITF, and MLANA. KIT, a receptor tyrosine kinase in the pigment production pathway, has been implicated in multiple mammalian deafness studies to date. A mutation (KI16 + 1037A) in the 5′ splice junction of exon 17 is known to cause sabino coloration, but was determined to not be linked to deafness. Alternatively, a heterozygous glycine to serine non‐synonymous mutation was identified within exon 6 of MITF, a transcription factor necessary for proper melanocyte differentiation, but doesn't appear to correlate to deafness. In addition, no mutations within exons were found in neither EDNRB (endothelin receptor type B) nor MLANA, also known as melan A or MART1 (melanoma antigen recognized by t‐cells). Further analysis of additional candidate genes, including PMEL17 and SOX10, and the inclusion of intronic sequence analysis for all candidate genes, will be necessary to determine the cause of cochleosaccular deafness, through mutations in these genes or others.

Gene network and canonical pathway analysis in canine myxomatous mitral valve disease: A microarray study

Myxomatous mitral valve disease (MMVD) is the single most common acquired heart disease of the dog and is particularly common in small pedigree breed dogs such as the Cavalier King Charles spaniel (CKCS). There are limited data on the mitral valve transcriptome and the aim of this study was to use the microarray technology in conjunction with bioinformatics platforms to analyse transcript changes in MMVD in CKCS compared to normal dogs (non-CKCS). Differentially expressed genes (n = 5397) were identified using cut-off settings of fold change, false discovery rate (FDR) and P < 0.05. In total, 4002 genes were annotated to a specific transcript in the Affymetrix canine database, and after further filtering, 591 annotated canine genes were identified: 322 (55%) were up-regulated and 269 (45%) were down-regulated. Canine microRNAs (cfa-miR; n = 59) were also identified.Gene ontology and network analysis platforms identified between six and 10 significantly different biological function clusters from which the following were selected as relevant to MMVD: inflammation, cell movement, cardiovascular development, extracellular matrix organisation and epithelial-to-mesenchymal (EMT) transition. Ingenuity Pathway Analysis identified three canonical pathways relevant to MMVD: caveolar-mediated endocytosis, remodelling of epithelial adherens junctions, and endothelin-1 signalling. Considering the biological relevance to MMVD, the gene families of importance with significant difference between groups included collagens, ADAMTS peptidases, proteoglycans, matrix metalloproteinases (MMPs) and their inhibitors, basement membrane components, cathepsin S, integrins, tight junction cell adhesion proteins, cadherins, other matrix-associated proteins, and members of the serotonin (5-HT)/transforming growth factor -β signalling pathway.

A novel missense mutation in ADAMTS10 in Norwegian Elkhound primary glaucoma.

Primary glaucoma is one of the most common causes of irreversible blindness both in humans and in dogs. Glaucoma is an optic neuropathy affecting the retinal ganglion cells and optic nerve, and elevated intraocular pressure is commonly associated with the disease. Glaucoma is broadly classified into primary open angle (POAG), primary closed angle (PCAG) and primary congenital glaucoma (PCG). Human glaucomas are genetically heterogeneous and multiple loci have been identified. Glaucoma affects several dog breeds but only three loci and one gene have been implicated so far. We have investigated the genetics of primary glaucoma in the Norwegian Elkhound (NE). We established a small pedigree around the affected NEs collected from Finland, US and UK and performed a genome-wide association study with 9 cases and 8 controls to map the glaucoma gene to 750 kb region on canine chromosome 20 (praw = 4.93×10−6, pgenome = 0.025). The associated region contains a previously identified glaucoma gene, ADAMTS10, which was subjected to mutation screening in the coding regions. A fully segregating missense mutation (p.A387T) in exon 9 was found in 14 cases and 572 unaffected NEs (pFisher = 3.5×10−27) with a high carrier frequency (25.3%). The mutation interrupts a highly conserved residue in the metalloprotease domain of ADAMTS10, likely affecting its functional capacity. Our study identifies the genetic cause of primary glaucoma in NEs and enables the development of a genetic test for breeding purposes. This study establishes also a new spontaneous canine model for glaucoma research to study the ADAMTS10 biology in optical neuropathy.

Novel Family-control Analysis Can Highly Prioritize Sequence Variants in Familial Cardiomyopathy

Detecting causative genes for cardiomyopathy among hundreds of candidate variants is still challenging, especially in small pedigrees where conventional linkage analysis is underpowered or even uninformative. To overcome this problem, we developed a novel approach, combining family data with population controls. Not knowing the true disease-causing variant, the idea is that a set of variants surrounding a true disease locus should have different frequencies in affected family members and unaffected controls.

First case of bovine psoroptic mange in Scotland in over 30 years

Cattle scab diagnosed in an imported calfDandy-Walker malformation in small pedigree Suffolk flockCorynebacterium renale causing renal disease in pregnant ewesAtrophic rhinitis in a previously unaffected pig herdInfectious bronchitis diagnosed in...

Recessive mutations in CAKUT and VACTERL association

Understanding the complex genetic makeup underlying congenital anomalies of the kidney and urinary tract (CAKUT) is of primary importance to improve diagnosis, stratify risk for later-onset complications, and develop therapeutic strategies. Saisawat et al. used homozygosity mapping coupled with next-generation sequencing to identify recessive mutations in TRAP1 in families with isolated CAKUT and with VACTERL association. This study points to a novel player in kidney development, possibly affecting apoptosis and endoplasmic reticulum stress signaling.

Epidemiology and Pathophysiology of Mitral Valve Prolapse

Mitral valve (MV) prolapse (MVP) is a common disorder, afflicting 2% to 3% of the general population.1,2 It is characterized by typical fibromyxomatous changes in the mitral leaflet tissue with superior displacement of 1 or both leaflets into the left atrium.3,4 With a prevalence of 2% to 3%, MVP is expected to affect ≈7.8 million individuals in the United States and >176 million people worldwide. MVP can be associated with significant mitral regurgitation (MR), bacterial endocarditis, congestive heart failure, and even sudden death.5–7 MVP is a clinical entity that is not fully understood, despite being known for more than a century. A “midsystolic click” was first described in 1887 by Cuffer and Barbillon.8 In 1963, Barlow and Pocock9 demonstrated the presence of MR by angiography in patients with the “click-murmur” syndrome. Criley et al10 subsequently coined the term MVP. MVP may be familial or sporadic. Despite being the most common cause of isolated MR requiring surgical repair,11 little is known about the genetic mechanisms underlying the pathogenesis and progression of MVP. Studies on the heritable features of MVP have been limited by the analysis of relatively small pedigrees and by self-referral and selection biases, including a preponderance of data from hospital-based cohorts.12,13 Nonetheless, the majority of data favor an autosomal-dominant pattern of inheritance in a large proportion of individuals with MVP.12,13 Despite the variability in clinical features, familial MVP might be considered a prevalent mendelian cardiac abnormality in humans. Although filamin-A has been identified as causing an X-linked form of MVP,14 the causative genes for the more common form of autosomal-dominant MVP have yet to be defined. In this review, we summarize our current knowledge of the diagnosis, epidemiology, prognosis, …

DEPDC5 mutations in families presenting as autosomal dominant nocturnal frontal lobe epilepsy.

To study the prevalence of DEPDC5 mutations in a series of 30 small European families with a phenotype compatible with autosomal dominant nocturnal frontal lobe epilepsy (ADNFLE). Thirty unrelated families referred with ADNFLE were recruited in France, Italy, Germany, Belgium, and Norway. Whole-exome sequencing was performed in 10 probands and direct sequencing of the DEPDC5 coding sequence in 20 probands. Testing for nonsense-mediated messenger RNA decay (NMD) was performed in lymphoblastic cells. Exome sequencing revealed a splice acceptor mutation (c.2355-2A>G) in DEPDC5 in the proband of a German family. In addition, 3 nonsense DEPDC5 mutations (p.Arg487*, p.Arg1087*, and p.Trp1369*) were detected in the probands of 2 French and one Belgian family. The nonsense mutations p.Arg487* and p.Arg1087* were targeted by NMD, leading to the degradation of the mutated transcripts. At the clinical level, 78% of the patients with DEPDC5 mutations were drug resistant. DEPDC5 loss-of-function mutations were found in 13% of the families with a presentation of ADNFLE. The rate of drug resistance was high in patients with DEPDC5 mutations. Small ADNFLE pedigrees with DEPDC5 mutations might actually represent a part of the broader familial focal epilepsy with variable foci phenotype.

Genetic analysis of high bone mass cases from the BARCOS cohort of Spanish postmenopausal women.

The aims of the study were to establish the prevalence of high bone mass (HBM) in a cohort of Spanish postmenopausal women (BARCOS) and to assess the contribution of LRP5 and DKK1 mutations and of common bone mineral density (BMD) variants to a HBM phenotype. Furthermore, we describe the expression of several osteoblast-specific and Wnt-pathway genes in primary osteoblasts from two HBM cases. A 0.6% of individuals (10/1600) displayed Z-scores in the HBM range (sum Z-score >4). While no mutation in the relevant exons of LRP5 was detected, a rare missense change in DKK1 was found (p.Y74F), which cosegregated with the phenotype in a small pedigree. Fifty-five BMD SNPs from Estrada et al. [NatGenet 44:491-501,2012] were genotyped in the HBM cases to obtain risk scores for each individual. In this small group of samples, Z-scores were found inversely related to risk scores, suggestive of a polygenic etiology. There was a single exception, which may be explained by a rare penetrant genetic variant, counterbalancing the additive effect of the risk alleles. The expression analysis in primary osteoblasts from two HBM cases and five controls suggested that IL6R, DLX3, TWIST1 and PPARG are negatively related to Z-score. One HBM case presented with high levels of RUNX2, while the other displayed very low SOX6. In conclusion, we provide evidence of lack of LRP5 mutations and of a putative HBM-causing mutation in DKK1. Additionally, we present SNP genotyping and expression results that suggest additive effects of several genes for HBM.

Using VAAST to Identify Disease-Associated Variants in Next-Generation Sequencing Data.

The VAAST pipeline is specifically designed to identify disease-associated alleles in next-generation sequencing data. In the protocols presented in this paper, we outline the best practices for variant prioritization using VAAST. Examples and test data are provided for case-control, small pedigree, and large pedigree analyses. These protocols will teach users the fundamentals of VAAST, VAAST 2.0, and pVAAST analyses.

Exome sequencing and in vitro studies identified podocalyxin as a candidate gene for focal and segmental glomerulosclerosis

Our understanding of focal and segmental glomerulosclerosis (FSGS) has advanced significantly from the studies of rare, monogenic forms of the disease. These studies have demonstrated the critical roles of multiple aspects of podocyte function in maintaining glomerular function. A substantial body of research has suggested that the integral membrane protein podocalyxin (PODXL) is required for proper function of podocytes, possibly by preserving the patency of the slit diaphragm by negative charge-based repulsion. Exome sequencing of affected cousins from an autosomal dominant pedigree with FSGS identified a co-segregating private variant, PODXL p.L442R, affecting the transmembrane region of the protein. Of the remaining 11 shared gene variants, two segregated with disease but their gene products were not detected in the glomerulus. In comparison to wild type, this disease-segregating PODXL variant facilitated dimerization. By contrast, this change does not alter protein stability, extracellular domain glycosylation, cell surface expression, global subcellular localization, or interaction with its intracellular binding partner ezrin. Thus, a variant form of PODXL remains the most likely candidate causing FSGS in one family with autosomal dominant inheritance, but its full effect on protein function remains unknown. Our work highlights the challenge faced in the clinical interpretation of whole exome data for small pedigrees with autosomal dominant diseases.