A novel 11 base pair deletion in KMT2C resulting in Kleefstra syndrome 2

BACKGROUND

Haploinsufficiency of the Lysine Methyltransferase 2C (KMT2C) gene results in the autosomal dominant disorder, Kleefstra syndrome 2. It is an extremely rare neurodevelopmental condition, with 14 previous reports describing varied clinical manifestations including dysmorphic features, delayed psychomotor development and delayed growth.

METHODS

Here, we describe a female with global developmental delay, attention deficit disorder, dyspraxia, short stature and subtle non-specific dysmorphic features. To identify causative mutations, whole exome sequencing was performed on the proband and her younger brother with discrete clinical presentation.

RESULTS

Whole exome sequencing identified a novel de novo heterozygous 11 bp deletion in KMT2C (c.1759_1769del), resulting in a frameshift mutation and early termination of the protein (p.Gln587SerfsTer7). This variant is the second-most N-terminal reported mutation, located 4171 amino acids upstream of the critical enzymatically active SET domain (required for chromatin modification and histone methylation).

CONCLUSION

The majority of the other reported mutations are frameshift mutations upstream of the SET domain and are predicted to result in protein truncation. It is thought that truncation of the SET domain, results functionally in an inability to modify chromatin through histone methylation. This report expands the clinical and genetic characterisation of Kleefstra syndrome 2.

Proof of concept for multiplex amplicon sequencing for mutation identification using the MinION nanopore sequencer

Rapid, cost-effective identification of genetic variants in small candidate genomic regions remains a challenge, particularly for less well equipped or lower throughput laboratories. The application of Oxford Nanopore Technologies’ MinION sequencer has the potential to fulfil this requirement. We demonstrate a proof of concept for a multiplexing assay that pools PCR amplicons for MinION sequencing to enable sequencing of multiple templates from multiple individuals, which could be applied to gene-targeted diagnostics. A combined strategy of barcoding and sample pooling was developed for simultaneous multiplex MinION sequencing of 100 PCR amplicons. The amplicons are family-specific, spanning a total of 30 loci in DNA isolated from 82 human neurodevelopmental cases and family members. The target regions were chosen for further interrogation because a potentially disease-causative variant had been identified in affected individuals following Illumina exome sequencing. The pooled MinION sequences were deconvoluted by aligning to custom references using the minimap2 aligner software. Our multiplexing approach produced an interpretable and expected sequence from 29 of the 30 targeted genetic loci. The sequence variant which was not correctly resolved in the MinION sequence was adjacent to a five nucleotide homopolymer. It is already known that homopolymers present a resolution problem with the MinION approach. Interestingly despite equimolar quantities of PCR amplicon pooled for sequencing, significant variation in the depth of coverage (127×–19,626×; mean = 8321×, std err = 452.99) was observed. We observed independent relationships between depth of coverage and target length, and depth of coverage and GC content. These relationships demonstrate biases of the MinION sequencer for longer templates and those with lower GC content. We demonstrate an efficient approach for variant discovery or confirmation from short DNA templates using the MinION sequencing device. With less than 130 × depth of coverage required for accurate genotyping, the methodology described here allows for rapid highly multiplexed targeted sequencing of large numbers of samples in a minimally equipped laboratory with a potential cost as much 200 × less than that from Sanger sequencing.

Contaminating DNA in human saliva alters the detection of variants from whole genome sequencing

Cells obtained from human saliva are commonly used as an alternative DNA source when blood is difficult or less convenient to collect. Although DNA extracted from saliva is considered to be of comparable quality to that derived from blood, recent studies have shown that non-human contaminating DNA derived from saliva can confound whole genome sequencing results. The most concerning complication is that non-human reads align to the human reference genome using standard methodology, which can critically affect the resulting variant genotypes identified in a genome. We identified clusters of anomalous variants in saliva DNA derived reads which aligned in an atypical manner. These reads had only short regions of identity to the human reference sequence, flanked by soft clipped sequence. Sequence comparisons of atypically aligning reads from eight human saliva-derived samples to RefSeq genomes revealed the majority to be of bacterial origin (63.46%). To partition the non-human reads during the alignment step, a decoy of the most prevalent bacterial genome sequences was designed and utilised. This reduced the number of atypically aligning reads when trialled on the eight saliva-derived samples by 44% and most importantly prevented the associated anomalous genotype calls. Saliva derived DNA is often contaminated by DNA from other species. This can lead to non-human reads aligning to the human reference genome using current alignment best-practices, impacting variant identification. This problem can be diminished by using a bacterial decoy in the alignment process.

Identifying, understanding, and correcting technical artifacts on the sex chromosomes in next-generation sequencing data

Background

Mammalian X and Y chromosomes share a common evolutionary origin and retain regions of high sequence similarity. Similar sequence content can confound the mapping of short next-generation sequencing reads to a reference genome. It is therefore possible that the presence of both sex chromosomes in a reference genome can cause technical artifacts in genomic data and affect downstream analyses and applications. Understanding this problem is critical for medical genomics and population genomic inference.

Results

Here, we characterize how sequence homology can affect analyses on the sex chromosomes and present XYalign, a new tool that (1) facilitates the inference of sex chromosome complement from next-generation sequencing data; (2) corrects erroneous read mapping on the sex chromosomes; and (3) tabulates and visualizes important metrics for quality control such as mapping quality, sequencing depth, and allele balance. We find that sequence homology affects read mapping on the sex chromosomes and this has downstream effects on variant calling. However, we show that XYalign can correct mismapping, resulting in more accurate variant calling. We also show how metrics output by XYalign can be used to identify XX and XY individuals across diverse sequencing experiments, including low- and high-coverage whole-genome sequencing, and exome sequencing. Finally, we discuss how the flexibility of the XYalign framework can be leveraged for other uses including the identification of aneuploidy on the autosomes. XYalign is available open source under the GNU General Public License (version 3).

Conclusions

Sex chromsome sequence homology causes the mismapping of short reads, which in turn affects downstream analyses. XYalign provides a reproducible framework to correct mismapping and improve variant calling on the sex chromsomes.

RBV: Read balance validator, a tool for prioritising copy number variations in germline conditions

The popularisation and decreased cost of genome resequencing has resulted in an increased use in molecular diagnostics. While there are a number of established and high quality bioinfomatic tools for identifying small genetic variants including single nucleotide variants and indels, currently there is no established standard for the detection of copy number variants (CNVs) from sequence data. The requirement for CNV detection from high throughput sequencing has resulted in the development of a large number of software packages. These tools typically utilise the sequence data characteristics: read depth, split reads, read pairs, and assembly-based techniques. However, the additional source of information from read balance (defined as relative proportion of reads of each allele at each position) has been underutilised in the existing applications. Here we present Read Balance Validator (RBV), a bioinformatic tool that uses read balance for prioritisation and validation of putative CNVs. The software simultaneously interrogates nominated regions for the presence of deletions or multiplications, and can differentiate larger CNVs from diploid regions. Additionally, the utility of RBV to test for inheritance of CNVs is demonstrated in this report. RBV is a CNV validation and prioritisation bioinformatic tool for both genome and exome sequencing available as a python package from https://github.com/whitneywhitford/RBV.

Evaluation of the performance of copy number variant prediction tools for the detection of deletions from whole genome sequencing data

BACKGROUND

Whole genome sequencing (WGS) has increased in popularity and decreased in cost over the past decade, rendering this approach as a viable and sensitive method for variant detection. In addition to its utility for single nucleotide variant detection, WGS data has the potential to detect Copy Number Variants (CNV) to fine resolution. Many CNV detection software packages have been developed exploiting four main types of data: read pair, split read, read depth, and assembly based methods. The aim of this study was to evaluate the efficiency of each of these main approaches in detecting germline deletions.

METHODS

WGS data and high confidence deletion calls for the individual NA12878 from the Genome in a Bottle consortium were the benchmark dataset. The performance of BreakDancer, CNVnator, Delly, FermiKit, and Pindel was assessed by comparing the accuracy and sensitivity of each software package in detecting deletions exceeding 1 kb.

RESULTS

There was considerable variability in the outputs of the different WGS CNV detection programs. The best performance was seen from BreakDancer and Delly, with 92.6% and 96.7% sensitivity, respectively and 34.5% and 68.5% false discovery rate (FDR), respectively. In comparison, Pindel, CNVnator, and FermiKit were less effective with sensitivities of 69.1%, 66.0%, and 15.8%, respectively and FDR of 91.3%, 69.0%, and 31.7%, respectively. Concordance across software packages was poor, with only 27 of the total 612 benchmark deletions identified by all five methodologies.

CONCLUSIONS

The WGS based CNV detection tools evaluated show disparate performance in identifying deletions ≥1 kb, particularly those utilising different input data characteristics. Software that exploits read pair based data had the highest sensitivity, namely BreakDancer and Delly. BreakDancer also had the second lowest false discovery rate. Therefore, in this analysis read pair methods (BreakDancer in particular) were the best performing approaches for the identification of deletions ≥1 kb, balancing accuracy and sensitivity. There is potential for improvement in the detection algorithms, particularly for reducing FDR. This analysis has validated the utility of WGS based CNV detection software to reliably identify deletions, and these findings will be of use when choosing appropriate software for deletion detection, in both research and diagnostic medicine.

Compound heterozygous SLC19A3 mutations further refine the critical promoter region for biotin-thiamine-responsive basal ganglia disease

Mutations in the gene SLC19A3 result in thiamine metabolism dysfunction syndrome 2, also known as biotin-thiamine-responsive basal ganglia disease (BTBGD). This neurometabolic disease typically presents in early childhood with progressive neurodegeneration, including confusion, seizures, and dysphagia, advancing to coma and death. Treatment is possible via supplement of biotin and/or thiamine, with early treatment resulting in significant lifelong improvements. Here we report two siblings who received a refined diagnosis of BTBGD following whole-genome sequencing. Both children inherited compound heterozygous mutations from unaffected parents; a missense single-nucleotide variant (p.G23V) in the first transmembrane domain of the protein, and a 4808-bp deletion in exon 1 encompassing the 5′ UTR and minimal promoter region. This deletion is the smallest promoter deletion reported to date, further defining the minimal promoter region of SLC19A3. Unfortunately, one of the siblings died prior to diagnosis, but the other is showing significant improvement after commencement of therapy. This case demonstrates the power of whole-genome sequencing for the identification of structural variants and subsequent diagnosis of rare neurodevelopmental disorders.

Compound Heterozygous Inheritance of Mutations in Coenzyme Q8A Results in Autosomal Recessive Cerebellar Ataxia and Coenzyme Q10 Deficiency in a Female Sib-Pair

Autosomal recessive ataxias are characterised by a fundamental loss in coordination of gait with associated atrophy of the cerebellum. There is significant clinical and genetic heterogeneity amongst inherited ataxias; however, an early molecular diagnosis is essential with low-risk treatments available for some of these conditions. We describe two female siblings who presented early in life with unsteady gait and cerebellar atrophy. Whole exome sequencing revealed compound heterozygous inheritance of two pathogenic mutations (p.Leu277Pro, c.1506+1G>A) in the coenzyme Q8A gene (COQ8A), a gene central to biosynthesis of coenzyme Q (CoQ). The paternally derived p.Leu277Pro mutation is predicted to disrupt a conserved motif in the substrate-binding pocket of the protein, resulting in inhibition of CoQ₁₀ production. The maternal c.1506+1G>A mutation destroys a canonical splice donor site in exon 12 affecting transcript processing and subsequent protein translation. Mutations in this gene can result in primary coenzyme Q₁₀ deficiency type 4, which is characterized by childhood onset of cerebellar ataxia and exercise intolerance, both of which were observed in this sib-pair. Muscle biopsies revealed unequivocally low levels of CoQ_10, and the siblings were subsequently established on a therapeutic dose of CoQ₁₀ with distinct clinical evidence of improvement after 1 year of treatment. This case emphasises the importance of an early and accurate molecular diagnosis for suspected inherited ataxias, particularly given the availability of approved treatments for some subtypes.

A Label-Free, Sensitive, Real-Time, Semiquantitative Electrochemical Measurement Method for DNA Polymerase Amplification (ePCR)

Oligonucleotide hybridization to a complementary sequence that is covalently attached to an electrochemically active conducting polymer (ECP) coating the working electrode of an electrochemical cell causes an increase in reaction impedance for the ferro-ferricyanide redox couple. We demonstrate the use of this effect to measure, in real time, the progress of DNA polymerase chain reaction (PCR) amplification of a minor component of a DNA extract. The forward primer is attached to the ECP. The solution contains other PCR components and the redox couple. Each cycle of amplification gives an easily measurable impedance increase. Target concentration can be estimated by cycle count to reach a threshold impedance. As proof of principle, we demonstrate an electrochemical real-time quantitative PCR (e-PCR) measurement in the total DNA extracted from chicken blood of an 844 base pair region of the mitochondrial Cytochrome c oxidase gene, present at ∼1 ppm of total DNA. We show that the detection and semiquantitation of as few as 2 copies/μL of target can be achieved within less than 10 PCR cycles.

The use of iodixanol as a density gradient material for separating human sperm from semen

Iodixanol, a new nonionic density gradient material with relatively low osmolality and high density, was evaluated to determine its suitability for the separation of human sperm from semen for their subsequent therapeutic use. Using a three-layer iodixanol gradient (approximately 1.17/1.15/1.05 g/mL), sperm were centrifuged at 1000g for 30 min and collected from the 1.05/1.15 interface. Using this method, a mean of 78% of the motile and 99% of the morphologically normal sperm originally present in the semen were recovered at the interface. There was no significant increase in the percentage of motile or morphologically normal sperm in the final preparation compared to the original semen. Sperm survived iodixanol density gradient centrifugation well, showing only modest declines in motility (18%) and velocity (35%) during a subsequent 24-h incubation period. When iodixanol was compared to Percoll density gradient centrifugation with semen from the same ejaculate, there was no significant difference between methods with regard to sperm yield, enrichment of motile or morphologically normal sperm in the preparation or sperm survival following separation. Iodixanol provides a suitable, nontoxic alternative to Percoll for the preparation of human sperm for therapeutic use.

Phosphorylatable and epitope-tagged human erythropoietins: utility and purification of native baculovirus-derived forms

The hematopoietic glycopeptide erythropoietin (EPO) is a prime regulator of red cell production in mammals, yet the precise nature of its interaction with specific cell surface receptors is poorly understood. Towards defining domains of EPO that are involved in receptor activation, we have developed (i) conditions for the expression of recombinant human EPO (rhEPO) at high levels in SF9 cells using modified 2- and 5-liter stirred reactors, (ii) a two-step procedure for the purification of this EPO without denaturation, and (iii) forms of EPO tagged with either a hemagglutinin influenza virus epitope or a consensus sequence for in vitro phosphorylation. Compared to EPO expressed in mammalian cells, rhEPO from SF9 cells in N-glycosylated with simple, neutral oligosaccharides of limited size, yet as purified presently using nondenaturing procedures, possesses exceptionally high in vitro activity (> or = 500,000 U/mg). Thus, this form of EPO should prove advantageous for direct physicochemical analyses. Regarding epitope-tagged and phosphorylatable EPOs, forms modified at the amino terminus (Ala1) fully retained receptor binding and in vitro biological activities. In contrast, forms modified at the carboxy terminus (Cys161) were inactive and did not compete for receptor binding, indicating that integrity of this domain is essential for receptor recognition. For active amino-terminal-modified forms, the specific binding of MAb 12CA5 to native HAI-EPO and the utility of 32P-labeled PHOS-EPO in receptor binding and internalization studies also were demonstrated. The development of these unique, highly active forms of human EPO should advance studies of essential interactions between this cytokine and its cell surface receptor.

Persistent reovirus infection of CHO cells resulting in virus resistance

We obtained a persistently infected line of Chinese hamster ovary cells by selection for resistance to reovirus infection. The cells were persistently infected by a population of viruses that were (i) cytopathic for parental chinese hamster ovary cells and (ii) similar to wild-type reovirus in molecular characteristics. The growth rate, plating efficiency, and morphology of the cells were altered. A large majority of the cells in the population were infected. There was no detectable interferon present in the medium. The cells were relatively resistant to a wide range of viruses.