Canfam GSD: De novo chromosome-length genome assembly of the German Shepherd Dog (Canis lupus familiaris) using a combination of long reads, optical mapping, and Hi-C

Background: The German Shepherd Dog (GSD) is one of the most common breeds on earth and has been bred for its utility and intelligence. It is often first choice for police and military work, as well as protection, disability assistance, and search-and-rescue. Yet, GSDs are well known to be susceptible to a range of genetic diseases that can interfere with their training. Such diseases are of particular concern when they occur later in life, and fully trained animals are not able to continue their duties. Findings: Here, we provide the draft genome sequence of a healthy German Shepherd female as a reference for future disease and evolutionary studies. We generated this improved canid reference genome (CanFam GSD) utilizing a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. The GSD assembly is ∼80 times as contiguous as the current canid reference genome (20.9 vs 0.267 Mb contig N50), containing far fewer gaps (306 vs 23,876) and fewer scaffolds (429 vs 3,310) than the current canid reference genome CanFamv3.1. Two chromosomes (4 and 35) are assembled into single scaffolds with no gaps. BUSCO analyses of the genome assembly results show that 93.0% of the conserved single-copy genes are complete in the GSD assembly compared with 92.2% for CanFam v3.1. Homology-based gene annotation increases this value to ∼99%. Detailed examination of the evolutionarily important pancreatic amylase region reveals that there are most likely 7 copies of the gene, indicative of a duplication of 4 ancestral copies and the disruption of 1 copy. Conclusions: GSD genome assembly and annotation were produced with major improvement in completeness, continuity, and quality over the existing canid reference. This resource will enable further research related to canine diseases, the evolutionary relationships of canids, and other aspects of canid biology

The Australian dingo is an early offshoot of modern breed dogs

Dogs are uniquely associated with human dispersal and bring transformational insight into the domestication process. Dingoes represent an intriguing case within canine evolution being geographically isolated for thousands of years. Here, we present a high-quality de novo assembly of a pure dingo (CanFam_DDS). We identified large chromosomal differences relative to the current dog reference (CanFam3.1) and confirmed no expanded pancreatic amylase gene as found in breed dogs. Phylogenetic analyses using variant pairwise matrices show that the dingo is distinct from five breed dogs with 100% bootstrap support when using Greenland wolf as the outgroup. Functionally, we observe differences in methylation patterns between the dingo and German shepherd dog genomes and differences in serum biochemistry and microbiome makeup. Our results suggest that distinct demographic and environmental conditions have shaped the dingo genome. In contrast, artificial human selection has likely shaped the genomes of domestic breed dogs after divergence from the dingo

Whole genome sequencing for the genetic diagnosis of heterogenous dystonia phenotypes

Introduction: Dystonia is a clinically and genetically heterogeneous disorder and a genetic cause is often difficult to elucidate. This is the first study to use whole genome sequencing (WGS) to investigate dystonia in a large sample of affected individuals. Methods: WGS was performed on 111 probands with heterogenous dystonia phenotypes. We performed analysis for coding and non-coding variants, copy number variants (CNVs), and structural variants (SVs). We assessed for an association between dystonia and 10 known dystonia risk variants. Results: A genetic diagnosis was obtained for 11.7% (13/111) of individuals. We found that a genetic diagnosis was more likely in those with an earlier age at onset, younger age at testing, and a combined dystonia phenotype. We identified pathogenic/likely-pathogenic variants in ADCY5 (n = 1), ATM (n = 1), GNAL (n = 2), GLB1 (n = 1), KMT2B (n = 2), PRKN (n = 2), PRRT2 (n = 1), SGCE (n = 2), and THAP1 (n = 1). CNVs were detected in 3 individuals. We found an association between the known risk variant ARSG rs11655081 and dystonia (p = 0.003). Conclusion: A genetic diagnosis was found in 11.7% of individuals with dystonia. The diagnostic yield was higher in those with an earlier age of onset, younger age at testing, and a combined dystonia phenotype. WGS may be particularly relevant for dystonia given that it allows for the detection of CNVs, which accounted for 23% of the genetically diagnosed cases. © 2019 The Author

Epigenetic profiling of heterochromatic satellite DNA

Zakrzewski F, Weisshaar B, Fuchs J, et al. Epigenetic profiling of heterochromatic satellite DNA. Chromosoma. 2011;120(4):409-422.Sugar beet (Beta vulgaris) chromosomes consist of large heterochromatic blocks in pericentromeric, centromeric, and intercalary regions comprised of two different highly abundant DNA satellite families. To investigate DNA methylation at single base resolution at heterochromatic regions, we applied a method for strand-specific bisulfite sequencing of more than 1,000 satellite monomers followed by statistical analyses. As a result, we uncovered diversity in the distribution of different methylation patterns in both satellite families. Heavily methylated CG and CHG (H=A, T, or C) sites occur more frequently in intercalary heterochromatin, while CHH sites, with the exception of CAA, are only sparsely methylated, in both intercalary and pericentromeric/centromeric heterochromatin. We show that the difference in DNA methylation intensity is correlated to unequal distribution of heterochromatic histone H3 methylation marks. While clusters of H3K9me2 were absent from pericentromeric heterochromatin and restricted only to intercalary heterochromatic regions, H3K9me1 and H3K27me1 were observed in all types of heterochromatin. By sequencing of a small RNA library consisting of 6.76 million small RNAs, we identified small interfering RNAs (siRNAs) of 24 nucleotides in size which originated from both strands of the satellite DNAs. We hypothesize an involvement of these siRNAs in the regulation of DNA and histone methylation for maintaining heterochromatin

Epigenetic profiling of heterochromatic satellite DNA

Zakrzewski F, Weisshaar B, Fuchs J, et al. Epigenetic profiling of heterochromatic satellite DNA. Chromosoma. 2011;120(4):409-422.Sugar beet (Beta vulgaris) chromosomes consist of large heterochromatic blocks in pericentromeric, centromeric, and intercalary regions comprised of two different highly abundant DNA satellite families. To investigate DNA methylation at single base resolution at heterochromatic regions, we applied a method for strand-specific bisulfite sequencing of more than 1,000 satellite monomers followed by statistical analyses. As a result, we uncovered diversity in the distribution of different methylation patterns in both satellite families. Heavily methylated CG and CHG (H=A, T, or C) sites occur more frequently in intercalary heterochromatin, while CHH sites, with the exception of CAA, are only sparsely methylated, in both intercalary and pericentromeric/centromeric heterochromatin. We show that the difference in DNA methylation intensity is correlated to unequal distribution of heterochromatic histone H3 methylation marks. While clusters of H3K9me2 were absent from pericentromeric heterochromatin and restricted only to intercalary heterochromatic regions, H3K9me1 and H3K27me1 were observed in all types of heterochromatin. By sequencing of a small RNA library consisting of 6.76 million small RNAs, we identified small interfering RNAs (siRNAs) of 24 nucleotides in size which originated from both strands of the satellite DNAs. We hypothesize an involvement of these siRNAs in the regulation of DNA and histone methylation for maintaining heterochromatin

Application of Genome Sequencing from Blood to Diagnose Mitochondrial Diseases

Mitochondrial diseases can be caused by pathogenic variants in nuclear or mitochondrial DNA-encoded genes that often lead to multisystemic symptoms and can have any mode of inheritance. Using a single test, Genome Sequencing (GS) can effectively identify variants in both genomes, but it has not yet been universally used as a first-line approach to diagnosing mitochondrial diseases due to related costs and challenges in data analysis. In this article, we report three patients with mitochondrial disease molecularly diagnosed through GS performed on DNA extracted from blood to demonstrate different diagnostic advantages of this technology, including the detection of a low-level heteroplasmic pathogenic variant, an intragenic nuclear DNA deletion, and a large mtDNA deletion. Current technical improvements and cost reductions are likely to lead to an expanded routine diagnostic usage of GS and of the complementary “Omic” technologies in mitochondrial diseases

Whole Genome Sequencing, Focused Assays and Functional Studies Increasing Understanding in Cryptic Inherited Retinal Dystrophies

The inherited retinal dystrophies (IRDs) are a clinically and genetically complex group of disorders primarily affecting the rod and cone photoreceptors or other retinal neuronal layers, with emerging therapies heralding the need for accurate molecular diagnosis. Targeted capture and panel-based strategies examining the partial or full exome deliver molecular diagnoses in many IRD families tested. However, approximately one in three families remain unsolved and unable to obtain personalised recurrence risk or access to new clinical trials or therapy. In this study, we investigated whole genome sequencing (WGS), focused assays and functional studies to assist with unsolved IRD cases and facilitate integration of these approaches to a broad molecular diagnostic clinical service. The WGS approach identified variants not covered or underinvestigated by targeted capture panel-based clinical testing strategies in six families. This included structural variants, with notable benefit of the WGS approach in repetitive regions demonstrated by a family with a hybrid gene and hemizygous missense variant involving the opsin genes, OPN1LW and OPN1MW. There was also benefit in investigation of the repetitive GC-rich ORF15 region of RPGR. Further molecular investigations were facilitated by focused assays in these regions. Deep intronic variants were identified in IQCB1 and ABCA4, with functional RNA based studies of the IQCB1 variant revealing activation of a cryptic splice acceptor site. While targeted capture panel-based methods are successful in achieving an efficient molecular diagnosis in a proportion of cases, this study highlights the additional benefit and clinical value that may be derived from WGS, focused assays and functional genomics in the highly heterogeneous IRDs

The diagnostic utility of genome sequencing in a pediatric cohort with suspected mitochondrial disease

Purpose: The utility of genome sequencing (GS) in the diagnosis of suspected pediatric mitochondrial disease (MD) was investigated. Methods: An Australian cohort of 40 pediatric patients with clinical features suggestive of MD were classified using the modified Nijmegen mitochondrial disease severity scoring into definite (17), probable (17), and possible (6) MD groups. Trio GS was performed using DNA extracted from patient and parent blood. Data were analyzed for single-nucleotide variants, indels, mitochondrial DNA variants, and structural variants. Results: A definitive MD gene molecular diagnosis was made in 15 cases and a likely MD molecular diagnosis in a further five cases. Causative mitochondrial DNA (mtDNA) variants were identified in four of these cases. Three potential novel MD genes were identified. In seven cases, causative variants were identified in known disease genes with no previous evidence of causing a primary MD. Diagnostic rates were higher in patients classified as having definite MD. Conclusion: GS efficiently identifies variants in MD genes of both nuclear and mitochondrial origin. A likely molecular diagnosis was identified in 67% of cases and a definitive molecular diagnosis achieved in 55% of cases. This study highlights the value of GS for a phenotypically and genetically heterogeneous disorder like MD

De Novo Variants Disrupting the HX Repeat Motif of ATN1 Cause a Recognizable Non-Progressive Neurocognitive Syndrome

Polyglutamine expansions in the transcriptional co-repressor Atrophin-1, encoded by ATN1, cause the neurodegenerative condition dentatorubral-pallidoluysian atrophy (DRPLA) via a proposed novel toxic gain of function. We present detailed phenotypic information on eight unrelated individuals who have de novo missense and insertion variants within a conserved 16-amino-acid "HX repeat" motif of ATN1. Each of the affected individuals has severe cognitive impairment and hypotonia, a recognizable facial gestalt, and variable congenital anomalies. However, they lack the progressive symptoms typical of DRPLA neurodegeneration. To distinguish this subset of affected individuals from the DRPLA diagnosis, we suggest using the term CHEDDA (congenital hypotonia, epilepsy, developmental delay, digit abnormalities) to classify the condition. CHEDDA-related variants alter the particular structural features of the HX repeat motif, suggesting that CHEDDA results from perturbation of the structural and functional integrity of the HX repeat. We found several non-homologous human genes containing similar motifs of eight to 10 HX repeat sequences, including RERE, where disruptive variants in this motif have also been linked to a separate condition that causes neurocognitive and congenital anomalies. These findings suggest that perturbation of the HX motif might explain other Mendelian human conditions