A clinically validated whole genome pipeline for structural variant detection and analysis

Expansion of the prenatal phenotype of Baraitser-Winter syndrome: Presentation of two cases of multiple congenital anomaly syndrome

Baraitser-Winter cerebrofrontofacial syndrome (BWCFF) is a variable multiple congenital anomaly condition, typically presenting postnatally with neurocognitive delays, distinctive facial features, cortical brain malformations, and in some, a variety of additional congenital malformations. However, only a few cases have reported the prenatal presentation of this syndrome. Here, we report two cases of BWCFF and their associated prenatal findings. One case presented with non-immune hydrops fetalis and a horseshoe kidney and was found to have a de novo heterozygous variant in ACTB (c.158A>G). The second case presented with gastroschisis, bilateral cleft lip and palate, and oligohydramnios, and was found to harbor a different de novo variant in ACTB (c.826G>A). Limited reports exist describing prenatally identified anomalies that include fetal growth restriction, increased nuchal fold, bilateral hydronephrosis, rocker bottom foot, talipes, cystic hygroma, omphalocele, and hydrops fetalis. In addition, only three of these cases have included detailed prenatal imaging findings. The two prenatal cases presented here demonstrate an expansion of the prenatal phenotype of BWCFF to include gastroschisis, lymphatic involvement, and oligohydramnios, which should each warrant consideration of this diagnosis in the setting of additional anomalies.

Long read sequencing on its way to the routine diagnostics of genetic diseases

The clinical application of technological progress in the identification of DNA alterations has always led to improvements of diagnostic yields in genetic medicine. At chromosome side, from cytogenetic techniques evaluating number and gross structural defects to genomic microarrays detecting cryptic copy number variants, and at molecular level, from Sanger method studying the nucleotide sequence of single genes to the high-throughput next-generation sequencing (NGS) technologies, resolution and sensitivity progressively increased expanding considerably the range of detectable DNA anomalies and alongside of Mendelian disorders with known genetic causes. However, particular genomic regions (i.e., repetitive and GC-rich sequences) are inefficiently analyzed by standard genetic tests, still relying on laborious, time-consuming and low-sensitive approaches (i.e., southern-blot for repeat expansion or long-PCR for genes with highly homologous pseudogenes), accounting for at least part of the patients with undiagnosed genetic disorders. Third generation sequencing, generating long reads with improved mappability, is more suitable for the detection of structural alterations and defects in hardly accessible genomic regions. Although recently implemented and not yet clinically available, long read sequencing (LRS) technologies have already shown their potential in genetic medicine research that might greatly impact on diagnostic yield and reporting times, through their translation to clinical settings. The main investigated LRS application concerns the identification of structural variants and repeat expansions, probably because techniques for their detection have not evolved as rapidly as those dedicated to single nucleotide variants (SNV) identification: gold standard analyses are karyotyping and microarrays for balanced and unbalanced chromosome rearrangements, respectively, and southern blot and repeat-primed PCR for the amplification and sizing of expanded alleles, impaired by limited resolution and sensitivity that have not been significantly improved by the advent of NGS. Nevertheless, more recently, with the increased accuracy provided by the latest product releases, LRS has been tested also for SNV detection, especially in genes with highly homologous pseudogenes and for haplotype reconstruction to assess the parental origin of alleles with de novo pathogenic variants. We provide a review of relevant recent scientific papers exploring LRS potential in the diagnosis of genetic diseases and its potential future applications in routine genetic testing.

Reanalysis of Trio Whole-Genome Sequencing Data Doubles the Yield in Autism Spectrum Disorder: De Novo Variants Present in Half

Autism spectrum disorder (ASD) is a common condition with lifelong implications. The last decade has seen dramatic improvements in DNA sequencing and related bioinformatics and databases. We analyzed the raw DNA sequencing files on the Variantyx® bioinformatics platform for the last 50 ASD patients evaluated with trio whole-genome sequencing (trio-WGS). "Qualified" variants were defined as coding, rare, and evolutionarily conserved. Primary Diagnostic Variants (PDV), additionally, were present in genes directly linked to ASD and matched clinical correlation. A PDV was identified in 34/50 (68%) of cases, including 25 (50%) cases with heterozygous de novo and 10 (20%) with inherited variants. De novo variants in genes directly associated with ASD were far more likely to be Qualifying than non-Qualifying versus a control group of genes (p = 0.0002), validating that most are indeed disease related. Sequence reanalysis increased diagnostic yield from 28% to 68%, mostly through inclusion of de novo PDVs in genes not yet reported as ASD associated. Thirty-three subjects (66%) had treatment recommendation(s) based on DNA analyses. Our results demonstrate a high yield of trio-WGS for revealing molecular diagnoses in ASD, which is greatly enhanced by reanalyzing DNA sequencing files. In contrast to previous reports, de novo variants dominate the findings, mostly representing novel conditions. This has implications to the cause and rising prevalence of autism.

ONT long-read WGS for variant discovery and orthogonal confirmation of short read WGS derived genetic variants in clinical genetic testing

Technological advances in Next-Generation Sequencing dramatically increased clinical efficiency of genetic testing, allowing detection of a wide variety of variants, from single nucleotide events to large structural aberrations. Whole Genome Sequencing (WGS) has allowed exploration of areas of the genome that might not have been targeted by other approaches, such as intergenic regions. A single technique detecting all genetic variants at once is intended to expedite the diagnostic process while making it more comprehensive and efficient. Nevertheless, there are still several shortcomings that cannot be effectively addressed by short read sequencing, such as determination of the precise size of short tandem repeat (STR) expansions, phasing of potentially compound recessive variants, resolution of some structural variants and exact determination of their boundaries, etc. Therefore, in some cases variants can only be tentatively detected by short reads sequencing and require orthogonal confirmation, particularly for clinical reporting purposes. Moreover, certain regulatory authorities, for example, New York state CLIA, require orthogonal confirmation of every reportable variant. Such orthogonal confirmations often involve numerous different techniques, not necessarily available in the same laboratory and not always performed in an expedited manner, thus negating the advantages of "one-technique-for-all" approach, and making the process lengthy, prone to logistical and analytical faults, and financially inefficient. Fortunately, those weak spots of short read sequencing can be compensated by long read technology that have comparable or better detection of some types of variants while lacking the mentioned above limitations of short read sequencing. At Variantyx we have developed an integrated clinical genetic testing approach, augmenting short read WGS-based variant detection with Oxford Nanopore Technologies (ONT) long read sequencing, providing simultaneous orthogonal confirmation of all types of variants with the additional benefit of improved identification of exact size and position of the detected aberrations. The validation study of this augmented test has demonstrated that Oxford Nanopore Technologies sequencing can efficiently verify multiple types of reportable variants, thus ensuring highly reliable detection and a quick turnaround time for WGS-based clinical genetic testing.

Familial Early-Onset Alzheimer's Caused by Novel Genetic Variant and APP Duplication: A Cross-Sectional Study

BACKGROUND

The clinical characteristics of symptomatic and asymptomatic carriers of early- onset autosomal dominant Alzheimer's (EOADAD) due to a yet-undescribed chromosomal rearrangement may add to the available body of knowledge about Alzheimer's disease and may enlighten novel and modifier genes. We report the clinical and genetic characteristics of asymptomatic and symptomatic individuals carrying a novel APP duplication rearrangement.

METHODS

Individuals belonging to a seven-generation pedigree with familial cognitive decline or intracerebral hemorrhages were recruited. Participants underwent medical, neurological, and neuropsychological evaluations. The genetic analysis included chromosomal microarray, Karyotype, fluorescence in situ hybridization, and whole genome sequencing.

RESULTS

Of 68 individuals, six females presented with dementia, and four males presented with intracerebral hemorrhage. Of these, nine were found to carry Chromosome 21 copy number gain (chr21:27,224,097-27,871,284, GRCh37/hg19) including the APP locus (APP-dup). In seven, Chromosome 5 copy number gain (Chr5: 24,786,234-29,446,070, GRCh37/hg19) (Chr5-CNG) cosegregated with the APP-dup. Both duplications co-localized to chromosome 18q21.1 and segregated in 25 pre-symptomatic carriers. Compared to non-carriers, asymptomatic carriers manifested cognitive decline in their mid-thirties. A third of the affected individuals carried a diagnosis of a dis-immune condition.

CONCLUSION

APP extra dosage, even in isolation and when located outside chromosome 21, is pathogenic. The clinical presentation of APP duplication varies and may be gender specific, i.e., ICH in males and cognitive-behavioral deterioration in females. The association with immune disorders is presently unclear but may prove relevant. The implication of Chr5-CNG co-segregation and the surrounding chromosome 18 genetic sequence needs further clarification.

Retrotransposition disrupting EBP in a girl and her mother with X-linked dominant chondrodysplasia punctata

X-linked dominant chondrodysplasia punctata (CDPX2) is a rare congenital disorder caused by pathogenic variants in EBP on Xp11.23. We encountered a girl and her mother with CDPX2-compatible phenotypes including punctiform calcification in the neonatal period of the girl, and asymmetric limb shortening and ichthyosis following the Blaschko lines in both subjects. Although Sanger direct sequencing failed to reveal a disease-causing variant in EBP, whole genome sequencing (WGS) followed by Manta analysis identified a ~ 4.5 kb insertion at EBP exon 2 of both subjects. The insertion was associated with the hallmarks of retrotransposition such as an antisense poly(A) tail, a target site duplication, and a consensus endonuclease cleavage site, and the inserted sequence harbored full-length SVA_F1 element with 5'- and 3'-transductions containing the Alu sequence. The results imply the relevance of retrotransposition to the human genetic diseases and the usefulness of WGS in the identification of retrotransposition.

Characterization of Genomic Variation from Lotus (Nelumbo Adans.) Mutants with Wide and Narrow Tepals

Compared with rose, chrysanthemum, and water lily, the absence of short-wide and long-narrow tepals of ornamental lotus (Nelumbo Adans.) limits the commercial value of flowers. In this study, the genomes of two groups of lotus mutants with wide-short and narrow-long tepals were resequenced to uncover the genomic variation and candidate genes associated with tepal shape. In group NL (short for N. lutea, containing two mutants and one control of N. lutea), 716,656 single nucleotide polymorphisms (SNPs) and 221,688 insertion-deletion mutations (Indels) were obtained, while 639,953 SNPs and 134,6118 Indels were obtained in group WSH (short for ‘Weishan Hong’, containing one mutant and two controls of N. nucifera ‘Weishan Hong’). Only a small proportion of these SNPs and Indels was mapped to exonic regions of genome: 1.92% and 0.47%, respectively, in the NL group, and 1.66% and 0.48%, respectively, in the WSH group. Gene Ontology (GO) analysis showed that out of 4890 (NL group) and 1272 (WSH group) annotated variant genes, 125 and 62 genes were enriched (Q < 0.05), respectively. Additionally, in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, 104 genes (NL group) and 35 genes (WSH group) were selected (p < 0.05). Finally, there were 306 candidate genes that were sieved to determine the development of tepal shape in lotus plants. It will be an essential reference for future identification of tepal-shaped control genes in lotus plants. This is the first comprehensive report of genomic variation controlling tepal shape in lotus, and the mutants in this study are promising materials for breeding novel lotus cultivars with special tepals.

Deciphering complex genome rearrangements in C. elegans using short-read whole genome sequencing

Genomic rearrangements cause congenital disorders, cancer, and complex diseases in human. Yet, they are still understudied in rare diseases because their detection is challenging, despite the advent of whole genome sequencing (WGS) technologies. Short-read (srWGS) and long-read WGS approaches are regularly compared, and the latter is commonly recommended in studies focusing on genomic rearrangements. However, srWGS is currently the most economical, accurate, and widely supported technology. In Caenorhabditis elegans (C. elegans), such variants, induced by various mutagenesis processes, have been used for decades to balance large genomic regions by preventing chromosomal crossover events and allowing the maintenance of lethal mutations. Interestingly, those chromosomal rearrangements have rarely been characterized on a molecular level. To evaluate the ability of srWGS to detect various types of complex genomic rearrangements, we sequenced three balancer strains using short-read Illumina technology. As we experimentally validated the breakpoints uncovered by srWGS, we showed that, by combining several types of analyses, srWGS enables the detection of a reciprocal translocation (eT1), a free duplication (sDp3), a large deletion (sC4), and chromoanagenesis events. Thus, applying srWGS to decipher real complex genomic rearrangements in model organisms may help designing efficient bioinformatics pipelines with systematic detection of complex rearrangements in human genomes.

Genome sequencing data analysis for rare disease gene discovery

Rare diseases occur in a smaller proportion of the general population, which is variedly defined as less than 200 000 individuals (US) or in less than 1 in 2000 individuals (Europe). Although rare, they collectively make up to approximately 7000 different disorders, with majority having a genetic origin, and affect roughly 300 million people globally. Most of the patients and their families undergo a long and frustrating diagnostic odyssey. However, advances in the field of genomics have started to facilitate the process of diagnosis, though it is hindered by the difficulty in genome data analysis and interpretation. A major impediment in diagnosis is in the understanding of the diverse approaches, tools and datasets available for variant prioritization, the most important step in the analysis of millions of variants to select a few potential variants. Here we present a review of the latest methodological developments and spectrum of tools available for rare disease genetic variant discovery and recommend appropriate data interpretation methods for variant prioritization. We have categorized the resources based on various steps of the variant interpretation workflow, starting from data processing, variant calling, annotation, filtration and finally prioritization, with a special emphasis on the last two steps. The methods discussed here pertain to elucidating the genetic basis of disease in individual patient cases via trio- or family-based analysis of the genome data. We advocate the use of a combination of tools and datasets and to follow multiple iterative approaches to elucidate the potential causative variant.

The Temple Grandin Genome: Comprehensive Analysis in a Scientist with High-Functioning Autism

Autism spectrum disorder (ASD) is a heterogeneous condition with a complex genetic etiology. The objective of this study is to identify the complex genetic factors that underlie the ASD phenotype and other clinical features of Professor Temple Grandin, an animal scientist and woman with high-functioning ASD. Identifying the underlying genetic cause for ASD can impact medical management, personalize services and treatment, and uncover other medical risks that are associated with the genetic diagnosis. Prof. Grandin underwent chromosomal microarray analysis, whole exome sequencing, and whole genome sequencing, as well as a comprehensive clinical and family history intake. The raw data were analyzed in order to identify possible genotype-phenotype correlations. Genetic testing identified variants in three genes (SHANK2, ALX1, and RELN) that are candidate risk factors for ASD. We identified variants in MEFV and WNT10A, reported to be disease-associated in previous studies, which are likely to contribute to some of her additional clinical features. Moreover, candidate variants in genes encoding metabolic enzymes and transporters were identified, some of which suggest potential therapies. This case report describes the genomic findings in Prof. Grandin and it serves as an example to discuss state-of-the-art clinical diagnostics for individuals with ASD, as well as the medical, logistical, and economic hurdles that are involved in clinical genetic testing for an individual on the autism spectrum.

Development of a new expanded next-generation sequencing panel for genetic diseases involved in dyslipidemia

The aim of this study was to provide an efficient tool: reliable, able to increase the molecular diagnosis performance, to facilitate the detection of copy number variants (CNV), to assess genetic risk scores (wGRS) and to offer the opportunity to explore candidate genes. Custom SeqCap EZ libraries, NextSeq500 sequencing and a homemade pipeline enable the analysis of 311 dyslipidemia-related genes. In the training group (48 DNA from patients with a well-established molecular diagnosis), this next-generation sequencing (NGS) workflow showed an analytical sensitivity >99% (n = 532 variants) without any false negative including a partial deletion of one exon. In the prospective group, from 25 DNA from patients without prior molecular analyses, 18 rare variants were identified in the first intention panel genes, allowing the diagnosis of monogenic dyslipidemia in 11 patients. In six other patients, the analysis of minor genes and wGRS determination provided a hypothesis to explain the dyslipidemia. Remaining data from the whole NGS workflow identified four patients with potentially deleterious variants. This NGS process gives a major opportunity to accede to an enhanced understanding of the genetic of dyslipidemia by simultaneous assessment of multiple genetic determinants.

Detection of non-targeted transgenes by whole-genome resequencing for gene-doping control

Gene doping has raised concerns in human and equestrian sports and the horseracing industry. There are two possible types of gene doping in the sports and racing industry: (1) administration of a gene-doping substance to postnatal animals and (2) generation of genetically engineered animals by modifying eggs. In this study, we aimed to identify genetically engineered animals by whole-genome resequencing (WGR) for gene-doping control. Transgenic cell lines, in which the erythropoietin gene (EPO) cDNA form was inserted into the genome of horse fibroblasts, were constructed as a model of genetically modified horse. Genome-wide screening of non-targeted transgenes was performed to find structural variation using DELLY based on split-read and paired-end algorithms and Control-FREEC based on read-depth algorithm. We detected the EPO transgene as an intron deletion in the WGR data by the split-read algorithm of DELLY. In addition, single-nucleotide polymorphisms and insertions/deletions artificially introduced in the EPO transgene were identified by WGR. Therefore, genome-wide screening using WGR can contribute to gene-doping control even if the targets are unknown. This is the first study to detect transgenes as intron deletions for gene-doping detection.

A family history of SCID and unrevealing WES: An approach to management and guidance of patients

Severe Combined Immunodeficiency (SCID) is a genetically heterogeneous group of disorders characterized by severe T cell lymphopenia and defective T and B cell function. Without prompt diagnosis and early intervention, patients with SCID typically die from infection within the first year of life. Advances in molecular genetics have led to rapid and efficient diagnosis of SCID cases, particularly when paired with newborn screening. However, some cases remain unsolved, and this is of particular relevance to families that plan to have more children. Here we report a patient who died from complications of SCID in whom whole exome sequencing failed to reveal a candidate variant. We describe how Sanger sequencing of parents was used to study the genomic regions that were poorly covered by WES, and how immune phenotyping results were used in the setting of genetic counseling.

Long-read sequencing for rare human genetic diseases

During the past decade, the search for pathogenic mutations in rare human genetic diseases has involved huge efforts to sequence coding regions, or the entire genome, using massively parallel short-read sequencers. However, the approximate current diagnostic rate is <50% using these approaches, and there remain many rare genetic diseases with unknown cause. There may be many reasons for this, but one plausible explanation is that the responsible mutations are in regions of the genome that are difficult to sequence using conventional technologies (e.g., tandem-repeat expansion or complex chromosomal structural aberrations). Despite the drawbacks of high cost and a shortage of standard analytical methods, several studies have analyzed pathogenic changes in the genome using long-read sequencers. The results of these studies provide hope that further application of long-read sequencers to identify the causative mutations in unsolved genetic diseases may expand our understanding of the human genome and diseases. Such approaches may also be applied to molecular diagnosis and therapeutic strategies for patients with genetic diseases in the future.