Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy.

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

Three generation families: Analysis of de novo variants in autism.

De novo variants (DNVs) analysis has proven to be a powerful approach to gene discovery in Autism Spectrum Disorder (ASD), which has not yet been shown in a Brazilian ASD cohort. The relevance of inherited rare variants has also been suggested, particularly in oligogenic models. We hypothesized that three-generation analyses of DNVs could provide new insights into the relevance of de novo and inherited variants across generations. To accomplish this goal, we performed whole-exome sequencing of 33 septet families composed of probands, parents, and grandparents (n = 231 individuals) and compared DNV rates (DNVr) between generations and those from two control cohorts. The DNVr in the probands (DNVr = 1.16) was marginally higher than in parents (DNVr = 0.60; p = 0.054), and in controls (DNVr = 0.68; p = 0.035, congenital heart disorder and DNVr = 0.70; p = 0.047, unaffected ASD siblings from Simons Simplex Collection). Moreover, most of the DNVs were found to have paternal origin in both generations (84.6%). Finally, we observed that 40% (6/15) of the DNVs in parents transmitted for probands are in ASD or ASD candidate genes, representing recently emerged risk variants to ASD in their families and suggest ZNF536, MSL2 and HDAC9 as ASD candidate genes. We did not observe an enrichment of risk variants nor sex bias of transmitted variants in the three generations, that can be due to sample size. These results further reinforce the relevance of de novo variants in ASD.

Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

Trio genome sequencing for developmental delay and pediatric heart conditions: A comparative microcost analysis.

PURPOSE: Genome sequencing (GS) can aid clinical management of multiple pediatric conditions. Insurers require accurate cost information to inform funding and implementation decisions. The objective was to compare the laboratory workflows and microcosts of trio GS testing in children with developmental delay (DD) and in children with cardiac conditions. METHODS: Cost items related to each step in trio GS (child and 2 parents) for both populations were identified and measured. Program costs over 5 years were estimated. Probabilistic and deterministic analyses were conducted. RESULTS: The mean cost per trio GS was CAD$6634.11 (95% CI = 6352.29-6913.40) for DD and CAD$8053.10 (95% CI = 7699.30-8558.10) for cardiac conditions. The 5-year program cost was CAD$28.11 million (95% CI = 26.91-29.29) for DD and CAD$5.63 million (95% CI = 5.38-5.98) for cardiac conditions. Supplies constituted the largest cost component for both populations. The higher cost per sample for the population with cardiac conditions was due to the inclusion of pharmacogenomics, higher bioinformatics labor costs, and a more labor intensive case review. CONCLUSION: This analysis indicated important variation in trio GS workflow and costs between pediatric populations in a single institution. Enhanced understanding of the clinical utility and costs of GS can inform harmonization and implementation decision-making.

Expanding the search for genetic biomarkers of Parkinson's disease into the living brain.

Altered gene expression related to Parkinson's Disease (PD) has not been described in the living brain, yet this information may support novel discovery pertinent to disease pathophysiology and treatment. This study compared the transcriptome in brain biopsies obtained from living PD and Control patients. To evaluate the novelty of this data, a comprehensive literature review also compared differentially expressed gene (DEGs) identified in the current study with those reported in PD cadaveric brain and peripheral tissues. RNA was extracted from rapidly cryopreserved frontal lobe specimens collected from PD and Control patients undergoing neurosurgical procedures. RNA sequencing (RNA-Seq) was performed and validated using quantitative polymerase chain reaction. DEG data was assessed using bioinformatics and subsequently included within a comparative analysis of PD RNA-Seq studies. 370 DEGs identified in living brain specimens reflected diverse gene groups and included key members of trophic signaling, apoptosis, inflammation and cell metabolism pathways. The comprehensive literature review yielded 7 RNA-Seq datasets generated from blood, skin and cadaveric brain but none from a living brain source. From the current dataset, 123 DEGs were identified only within the living brain and 267 DEGs were either newly found or had distinct directional change in living brain relative to other tissues. This is the first known study to analyze the transcriptome in brain tissue from living PD and Control patients. The data produced using these methods offer a unique, unexplored resource with potential to advance insight into the genetic associations of PD.

The Cardiac Genome Clinic: implementing genome sequencing in pediatric heart disease.

PURPOSE: This study investigated the diagnostic utility of nontargeted genomic testing in patients with pediatric heart disease. METHODS: We analyzed genome sequencing data of 111 families with cardiac lesions for rare, disease-associated variation. RESULTS: In 14 families (12.6%), we identified causative variants: seven were de novo (ANKRD11, KMT2D, NR2F2, POGZ, PTPN11, PURA, SALL1) and six were inherited from parents with no or subclinical heart phenotypes (FLT4, DNAH9, MYH11, NEXMIF, NIPBL, PTPN11). Outcome of the testing was associated with the presence of extracardiac features (p = 0.02), but not a positive family history for cardiac lesions (p = 0.67). We also report novel plausible gene-disease associations for tetralogy of Fallot/pulmonary stenosis (CDC42BPA, FGD5), hypoplastic left or right heart (SMARCC1, TLN2, TRPM4, VASP), congenitally corrected transposition of the great arteries (UBXN10), and early-onset cardiomyopathy (TPCN1). The identified candidate genes have critical functions in heart development, such as angiogenesis, mechanotransduction, regulation of heart size, chromatin remodeling, or ciliogenesis. CONCLUSION: This data set demonstrates the diagnostic and scientific value of genome sequencing in pediatric heart disease, anticipating its role as a first-tier diagnostic test. The genetic heterogeneity will necessitate large-scale genomic initiatives for delineating novel gene-disease associations.

The CF Canada-Sick Kids Program in individual CF therapy: A resource for the advancement of personalized medicine in CF.

BACKGROUND: Therapies targeting certain CFTR mutants have been approved, yet variations in clinical response highlight the need for in-vitro and genetic tools that predict patient-specific clinical outcomes. Toward this goal, the CF Canada-Sick Kids Program in Individual CF Therapy (CFIT) is generating a "first of its kind", comprehensive resource containing patient-specific cell cultures and data from 100 CF individuals that will enable modeling of therapeutic responses. METHODS: The CFIT program is generating: 1) nasal cells from drug naïve patients suitable for culture and the study of drug responses in vitro, 2) matched gene expression data obtained by sequencing the RNA from the primary nasal tissue, 3) whole genome sequencing of blood derived DNA from each of the 100 participants, 4) induced pluripotent stem cells (iPSCs) generated from each participant's blood sample, 5) CRISPR-edited isogenic control iPSC lines and 6) prospective clinical data from patients treated with CF modulators. RESULTS: To date, we have recruited 57 of 100 individuals to CFIT, most of whom are homozygous for F508del (to assess in-vitro: in-vivo correlations with respect to ORKAMBI response) or heterozygous for F508del and a minimal function mutation. In addition, several donors are homozygous for rare nonsense and missense mutations. Nasal epithelial cell cultures and matched iPSC lines are available for many of these donors. CONCLUSIONS: This accessible resource will enable development of tools that predict individual outcomes to current and emerging modulators targeting F508del-CFTR and facilitate therapy discovery for rare CF causing mutations.

The Personal Genome Project Canada: findings from whole genome sequences of the inaugural 56 participants.

BACKGROUND: The Personal Genome Project Canada is a comprehensive public data resource that integrates whole genome sequencing data and health information. We describe genomic variation identified in the initial recruitment cohort of 56 volunteers. METHODS: Volunteers were screened for eligibility and provided informed consent for open data sharing. Using blood DNA, we performed whole genome sequencing and identified all possible classes of DNA variants. A genetic counsellor explained the implication of the results to each participant. RESULTS: Whole genome sequencing of the first 56 participants identified 207 662 805 sequence variants and 27 494 copy number variations. We analyzed a prioritized disease-associated data set (n = 1606 variants) according to standardized guidelines, and interpreted 19 variants in 14 participants (25%) as having obvious health implications. Six of these variants (e.g., in BRCA1 or mosaic loss of an X chromosome) were pathogenic or likely pathogenic. Seven were risk factors for cancer, cardiovascular or neurobehavioural conditions. Four other variants - associated with cancer, cardiac or neurodegenerative phenotypes - remained of uncertain significance because of discrepancies among databases. We also identified a large structural chromosome aberration and a likely pathogenic mitochondrial variant. There were 172 recessive disease alleles (e.g., 5 individuals carried mutations for cystic fibrosis). Pharmacogenomics analyses revealed another 3.9 potentially relevant genotypes per individual. INTERPRETATION: Our analyses identified a spectrum of genetic variants with potential health impact in 25% of participants. When also considering recessive alleles and variants with potential pharmacologic relevance, all 56 participants had medically relevant findings. Although access is mostly limited to research, whole genome sequencing can provide specific and novel information with the potential of major impact for health care.

A Comprehensive Workflow for Read Depth-Based Identification of Copy-Number Variation from Whole-Genome Sequence Data.

A remaining hurdle to whole-genome sequencing (WGS) becoming a first-tier genetic test has been accurate detection of copy-number variations (CNVs). Here, we used several datasets to empirically develop a detailed workflow for identifying germline CNVs >1 kb from short-read WGS data using read depth-based algorithms. Our workflow is comprehensive in that it addresses all stages of the CNV-detection process, including DNA library preparation, sequencing, quality control, reference mapping, and computational CNV identification. We used our workflow to detect rare, genic CNVs in individuals with autism spectrum disorder (ASD), and 120/120 such CNVs tested using orthogonal methods were successfully confirmed. We also identified 71 putative genic de novo CNVs in this cohort, which had a confirmation rate of 70%; the remainder were incorrectly identified as de novo due to false positives in the proband (7%) or parental false negatives (23%). In individuals with an ASD diagnosis in which both microarray and WGS experiments were performed, our workflow detected all clinically relevant CNVs identified by microarrays, as well as additional potentially pathogenic CNVs < 20 kb. Thus, CNVs of clinical relevance can be discovered from WGS with a detection rate exceeding microarrays, positioning WGS as a single assay for genetic variation detection.

A homozygous mutation in the stem II domain of RNU4ATAC causes typical Roifman syndrome.

Roifman syndrome (OMIM# 616651) is a complex syndrome encompassing skeletal dysplasia, immunodeficiency, retinal dystrophy and developmental delay, and is caused by compound heterozygous mutations involving the Stem II region and one of the other domains of the RNU4ATAC gene. This small nuclear RNA gene is essential for minor intron splicing. The Canadian Centre for Primary Immunodeficiency Registry and Repository were used to derive patient information as well as tissues. Utilising RNA sequencing methodologies, we analysed samples from patients with Roifman syndrome and assessed intron retention. We demonstrate that a homozygous mutation in Stem II is sufficient to cause the full spectrum of features associated with typical Roifman syndrome. Further, we demonstrate the same pattern of aberration in minor intron retention as found in cases with compound heterozygous mutations.

A microcosting and cost-consequence analysis of clinical genomic testing strategies in autism spectrum disorder.

PurposeWhole-exome (WES) and whole-genome sequencing (WGS) increase the diagnostic yield in autism spectrum disorder (ASD) compared to chromosomal microarray (CMA), but there have been no comprehensive cost analyses. The objective was to perform such an assessment of CMA, WES, and WGS and compare the incremental cost per additional positive finding in hypothetical testing scenarios.MethodsFive-year patient and program costs were estimated from an institutional perspective. WES and WGS estimates were based on HiSeq 2500 with an additional WGS estimate for HiSeq X platforms. Parameter uncertainty was assessed with probabilistic and deterministic sensitivity analysis.ResultsThe cost per ASD sample was CAD$1,655 (95% CI: 1,611; 1,699) for WES, CAD$2,851 (95% CI: 2,750; 2,956) for WGS on HiSeq X, and CAD$5,519 (95% CI: 5,244; 5,785) on HiSeq 2500, compared to CAD$744 (95% CI 714, 773) for CMA. The incremental cost was over CAD$25,000 per additional positive finding if CMA was replaced by newer technology.ConclusionWhile costs for WES and WGS remain high, future reductions in material and equipment costs, and increased understanding of newly discovered variants and variants of unknown significance will lead to improved value.

Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder.

We are performing whole-genome sequencing of families with autism spectrum disorder (ASD) to build a resource (MSSNG) for subcategorizing the phenotypes and underlying genetic factors involved. Here we report sequencing of 5,205 samples from families with ASD, accompanied by clinical information, creating a database accessible on a cloud platform and through a controlled-access internet portal. We found an average of 73.8 de novo single nucleotide variants and 12.6 de novo insertions and deletions or copy number variations per ASD subject. We identified 18 new candidate ASD-risk genes and found that participants bearing mutations in susceptibility genes had significantly lower adaptive ability (P = 6 × 10-4). In 294 of 2,620 (11.2%) of ASD cases, a molecular basis could be determined and 7.2% of these carried copy number variations and/or chromosomal abnormalities, emphasizing the importance of detecting all forms of genetic variation as diagnostic and therapeutic targets in ASD.

De Novo Genome and Transcriptome Assembly of the Canadian Beaver (Castor canadensis).

The Canadian beaver (Castor canadensis) is the largest indigenous rodent in North America. We report a draft annotated assembly of the beaver genome, the first for a large rodent and the first mammalian genome assembled directly from uncorrected and moderate coverage (< 30 ×) long reads generated by single-molecule sequencing. The genome size is 2.7 Gb estimated by k-mer analysis. We assembled the beaver genome using the new Canu assembler optimized for noisy reads. The resulting assembly was refined using Pilon supported by short reads (80 ×) and checked for accuracy by congruency against an independent short read assembly. We scaffolded the assembly using the exon-gene models derived from 9805 full-length open reading frames (FL-ORFs) constructed from the beaver leukocyte and muscle transcriptomes. The final assembly comprised 22,515 contigs with an N50 of 278,680 bp and an N50-scaffold of 317,558 bp. Maximum contig and scaffold lengths were 3.3 and 4.2 Mb, respectively, with a combined scaffold length representing 92% of the estimated genome size. The completeness and accuracy of the scaffold assembly was demonstrated by the precise exon placement for 91.1% of the 9805 assembled FL-ORFs and 83.1% of the BUSCO (Benchmarking Universal Single-Copy Orthologs) gene set used to assess the quality of genome assemblies. Well-represented were genes involved in dentition and enamel deposition, defining characteristics of rodents with which the beaver is well-endowed. The study provides insights for genome assembly and an important genomics resource for Castoridae and rodent evolutionary biology.

The genetic diversity of Epstein-Barr virus in the setting of transplantation relative to non-transplant settings: A feasibility study.

This study examines EBV strains from transplant patients and patients with IM by sequencing major EBV genes. We also used NGS to detect EBV DNA within total genomic DNA, and to evaluate its genetic variation. Sanger sequencing of major EBV genes was used to compare SNVs from samples taken from transplant patients vs. patients with IM. We sequenced EBV DNA from a healthy EBV-seropositive individual on a HiSeq 2000 instrument. Data were mapped to the EBV reference genomes (AG876 and B95-8). The number of EBNA2 SNVs was higher than for EBNA1 and the other genes sequenced within comparable reference coordinates. For EBNA2, there was a median of 15 SNV among transplant samples compared with 10 among IM samples (p = 0.036). EBNA1 showed little variation between samples. For NGS, we identified 640 and 892 variants at an unadjusted p value of 5 × 10(-8) for AG876 and B95-8 genomes, respectively. We used complementary sequence strategies to examine EBV genetic diversity and its application to transplantation. The results provide the framework for further characterization of EBV strains and related outcomes after organ transplantation.

Compound heterozygous mutations in the noncoding RNU4ATAC cause Roifman Syndrome by disrupting minor intron splicing.

Roifman Syndrome is a rare congenital disorder characterized by growth retardation, cognitive delay, spondyloepiphyseal dysplasia and antibody deficiency. Here we utilize whole-genome sequencing of Roifman Syndrome patients to reveal compound heterozygous rare variants that disrupt highly conserved positions of the RNU4ATAC small nuclear RNA gene, a minor spliceosome component that is essential for minor intron splicing. Targeted sequencing confirms allele segregation in six cases from four unrelated families. RNU4ATAC rare variants have been recently reported to cause microcephalic osteodysplastic primordial dwarfism, type I (MOPD1), whose phenotype is distinct from Roifman Syndrome. Strikingly, all six of the Roifman Syndrome cases have one variant that overlaps MOPD1-implicated structural elements, while the other variant overlaps a highly conserved structural element not previously implicated in disease. RNA-seq analysis confirms extensive and specific defects of minor intron splicing. Available allele frequency data suggest that recessive genetic disorders caused by RNU4ATAC rare variants may be more prevalent than previously reported.

Whole-exome analysis of foetal autopsy tissue reveals a frameshift mutation in OBSL1, consistent with a diagnosis of 3-M Syndrome.

BACKGROUND: We report a consanguineous couple that has experienced three consecutive pregnancy losses following the foetal ultrasound finding of short limbs. Post-termination examination revealed no skeletal dysplasia, but some subtle proximal limb shortening in two foetuses, and a spectrum of mildly dysmorphic features. Karyotype was normal in all three foetuses (46, XX) and comparative genomic hybridization microarray analysis detected no pathogenic copy number variants. RESULTS: Whole-exome sequencing and genome-wide homozygosity mapping revealed a previously reported frameshift mutation in the OBSL1 gene (c.1273insA p.T425nfsX40), consistent with a diagnosis of 3-M Syndrome 2 (OMIM #612921), which had not been anticipated from the clinical findings. CONCLUSIONS: Our study provides novel insight into the early clinical manifestations of this form of 3-M syndrome, and demonstrates the utility of whole exome sequencing as a tool for prenatal diagnosis in particular when there is a family history suggestive of a recurrent set of clinical symptoms.

Fatal combined immunodeficiency associated with heterozygous mutation in STAT1.

BACKGROUND: Mutations in the gene for the signal transducer and activator of transcription 1, STAT1, have been shown to be associated with death at an early age due to overwhelming viral infection (complete STAT1 deficiency) or, more commonly, selective deficiencies to mycobacterial or fungal infection (typically heterozygous STAT1 mutations). OBJECTIVES: To define the molecular basis of progressive combined immunodeficiency in a group of patients with fatal infections. METHODS: We studied a group of unrelated patients who displayed an unusual progressive form of combined immunodeficiency. Whole exome sequencing assisted in confirming a common genetic defect in this group, which consisted of a heterozygous mutation of the STAT1 gene. STAT1 protein level as well as function was assessed, and a detailed evaluation of the immune system, including analysis of thymus tissue, was performed. RESULTS: Patients were found to carry de novo heterozygous mutations in STAT1 encoding T385A, I294T, or C284R amino acid substitutions. STAT1 expression appeared significantly decreased as a result of these changes but not completely absent, with diminished signaling responses. This group display progressive loss in lymphocyte number and function accompanied by increasing autoimmune features as well as severe, fatal infections. CONCLUSIONS: These findings show that some heterozygous aberrations of STAT1 can be associated with progressive combined immunodeficiency, quite distinct from the limited susceptibilities to infection previously reported for heterozygous STAT1 mutations. These mutations were not inherited, rather, arose de novo in each case. Accompanied by significant patient mortality, this finding suggests that this class of STAT1 mutation is ultimately fatal due to overwhelming infection.

Adult siblings with homozygous G6PC3 mutations expand our understanding of the severe congenital neutropenia type 4 (SCN4) phenotype.

BACKGROUND: Severe congenital neutropenia type 4 (SCN4) is an autosomal recessive disorder caused by mutations in the third subunit of the enzyme glucose-6-phosphatase (G6PC3). Its core features are congenital neutropenia and a prominent venous skin pattern, and affected individuals have variable birth defects. Oculocutaneous albinism type 4 (OCA4) is caused by autosomal recessive mutations in SLC45A2. METHODS: We report a sister and brother from Newfoundland, Canada with complex phenotypes. The sister was previously reported by Cullinane et al., 2011. We performed homozygosity mapping, next generation sequencing and conventional Sanger sequencing to identify mutations that cause the phenotype in this family. We have also summarized clinical data from 49 previously reported SCN4 cases with overlapping phenotypes and interpret the medical histories of these siblings in the context of the literature. RESULTS: The siblings' phenotype is due in part to a homozygous mutation in G6PC3, [c.829C > T, p.Gln277X]. Their ages are 38 and 37 years respectively and they are the oldest SCN4 patients published to date. Both presented with congenital neutropenia and later developed Crohn disease. We suggest that the latter is a previously unrecognized SCN4 manifestation and that not all affected individuals have an intellectual disability. The sister also has a homozygous mutation in SLC45A2, which explains her severe oculocutaneous hypopigmentation. Her brother carried one SLC45A2 mutation and was diagnosed with "partial OCA" in childhood. CONCLUSIONS: This family highlights that apparently novel syndromes can in fact be caused by two known autosomal recessive disorders.

The enigmatic monotypic crab plover Dromas ardeola is closely related to pratincoles and coursers (Aves, Charadriiformes, Glareolidae).

The phylogenetic placement of the monotypic crab plover Dromasardeola (Aves, Charadriiformes) remains controversial. Phylogenetic analysis of anatomical and behavioral traits using phenetic and cladistic methods of tree inference have resulted in conflicting tree topologies, suggesting a close association of Dromas to members of different suborders and lineages within Charadriiformes. Here, we revisited the issue by applying Bayesian and parsimony methods of tree inference to 2,012 anatomical and 5,183 molecular characters to a set of 22 shorebird genera (including Turnix). Our results suggest that Bayesian analysis of anatomical characters does not resolve the phylogenetic relationship of shorebirds with strong statistical support. In contrast, Bayesian and parsimony tree inference from molecular data provided much stronger support for the phylogenetic relationships within shorebirds, and support a sister relationship of Dromas to Glareolidae (pratincoles and coursers), in agreement with previously published DNA-DNA hybridization studies.

The enigmatic monotypic crab plover Dromas ardeola is closely related to pratincoles and coursers (Aves, Charadriiformes, Glareolidae)

The phylogenetic placement of the monotypic crab plover Dromas ardeola (Aves, Charadriiformes) remains controversial. Phylogenetic analysis of anatomical and behavioral traits using phenetic and cladistic methods of tree inference have resulted in conflicting tree topologies, suggesting a close association of Dromas to members of different suborders and lineages within Charadriiformes. Here, we revisited the issue by applying Bayesian and parsimony methods of tree inference to 2,012 anatomical and 5,183 molecular characters to a set of 22 shorebird genera (including Turnix). Our results suggest that Bayesian analysis of anatomical characters does not resolve the phylogenetic relationship of shorebirds with strong statistical support. In contrast, Bayesian and parsimony tree inference from molecular data provided much stronger support for the phylogenetic relationships within shorebirds, and support a sister relationship of Dromas to Glareolidae (pratincoles and coursers), in agreement with previously published DNA-DNA hybridization studies.