A sensitive and affordable multiplex RT-qPCR assay for SARS-CoV-2 detection.

With the ongoing COVID-19 (Coronavirus Disease 2019) pandemic, caused by the novel coronavirus SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2), there is a need for sensitive, specific, and affordable diagnostic tests to identify infected individuals, not all of whom are symptomatic. The most sensitive test involves the detection of viral RNA using RT-qPCR (quantitative reverse transcription PCR), with many commercial kits now available for this purpose. However, these are expensive, and supply of such kits in sufficient numbers cannot always be guaranteed. We therefore developed a multiplex assay using well-established SARS-CoV-2 targets alongside a human cellular control (RPP30) and a viral spike-in control (Phocine Herpes Virus 1 [PhHV-1]), which monitor sample quality and nucleic acid extraction efficiency, respectively. Here, we establish that this test performs as well as widely used commercial assays, but at substantially reduced cost. Furthermore, we demonstrate >1,000-fold variability in material routinely collected by combined nose and throat swabbing and establish a statistically significant correlation between the detected level of human and SARS-CoV-2 nucleic acids. The inclusion of the human control probe in our assay therefore provides a quantitative measure of sample quality that could help reduce false-negative rates. We demonstrate the feasibility of establishing a robust RT-qPCR assay at approximately 10% of the cost of equivalent commercial assays, which could benefit low-resource environments and make high-volume testing affordable.

Targeted next-generation sequencing makes new molecular diagnoses and expands genotype-phenotype relationship in Ehlers-Danlos syndrome.

PURPOSE: Ehlers-Danlos syndrome (EDS) comprises a group of overlapping hereditary disorders of connective tissue with significant morbidity and mortality, including major vascular complications. We sought to identify the diagnostic utility of a next-generation sequencing (NGS) panel in a mixed EDS cohort. METHODS: We developed and applied PCR-based NGS assays for targeted, unbiased sequencing of 12 collagen and aortopathy genes to a cohort of 177 unrelated EDS patients. Variants were scored blind to previous genetic testing and then compared with results of previous Sanger sequencing. RESULTS: Twenty-eight pathogenic variants in COL5A1/2, COL3A1, FBN1, and COL1A1 and four likely pathogenic variants in COL1A1, TGFBR1/2, and SMAD3 were identified by the NGS assays. These included all previously detected single-nucleotide and other short pathogenic variants in these genes, and seven newly detected pathogenic or likely pathogenic variants leading to clinically significant diagnostic revisions. Twenty-two variants of uncertain significance were identified, seven of which were in aortopathy genes and required clinical follow-up. CONCLUSION: Unbiased NGS-based sequencing made new molecular diagnoses outside the expected EDS genotype-phenotype relationship and identified previously undetected clinically actionable variants in aortopathy susceptibility genes. These data may be of value in guiding future clinical pathways for genetic diagnosis in EDS.Genet Med 18 11, 1119-1127.

Inferring mechanisms of copy number change from haplotype structures at the human DEFA1A3 locus.

BACKGROUND: The determination of structural haplotypes at copy number variable regions can indicate the mechanisms responsible for changes in copy number, as well as explain the relationship between gene copy number and expression. However, obtaining spatial information at regions displaying extensive copy number variation, such as the DEFA1A3 locus, is complex, because of the difficulty in the phasing and assembly of these regions. The DEFA1A3 locus is intriguing in that it falls within a region of high linkage disequilibrium, despite its high variability in copy number (n = 3-16); hence, the mechanisms responsible for changes in copy number at this locus are unclear. RESULTS: In this study, a region flanking the DEFA1A3 locus was sequenced across 120 independent haplotypes with European ancestry, identifying five common classes of DEFA1A3 haplotype. Assigning DEFA1A3 class to haplotypes within the 1000 Genomes project highlights a significant difference in DEFA1A3 class frequencies between populations with different ancestry. The features of each DEFA1A3 class, for example, the associated DEFA1A3 copy numbers, were initially assessed in a European cohort (n = 599) and replicated in the 1000 Genomes samples, showing within-class similarity, but between-class and between-population differences in the features of the DEFA1A3 locus. Emulsion haplotype fusion-PCR was used to generate 61 structural haplotypes at the DEFA1A3 locus, showing a high within-class similarity in structure. CONCLUSIONS: Structural haplotypes across the DEFA1A3 locus indicate that intra-allelic rearrangement is the predominant mechanism responsible for changes in DEFA1A3 copy number, explaining the conservation of linkage disequilibrium across the locus. The identification of common structural haplotypes at the DEFA1A3 locus could aid studies into how DEFA1A3 copy number influences expression, which is currently unclear.

Whole genome sequencing enhances molecular diagnosis of primary ciliary dyskinesia.

BACKGROUND: Primary ciliary dyskinesia (PCD) is a genetic disorder affecting motile cilia. Most cases are inherited recessively, due to variants in >50 genes that result in abnormal or absent motile cilia. This leads to chronic upper and lower airway disease, subfertility, and laterality defects. Given overlapping clinical features and genetic heterogeneity, diagnosis can be difficult and often occurs late. Of those tested an estimated 30% of genetically screened PCD patients still lack a molecular diagnosis. A molecular diagnosis allows for appropriate clinical management including prediction of phenotypic features correlated to genotype. Here, we aimed to identify how readily a genetic diagnosis could be made using whole genome sequencing (WGS) to facilitate identification of pathogenic variants in known genes as well as novel PCD candidate genes. METHODS: WGS was used to screen for pathogenic variants in eight patients with PCD. RESULTS: 7/8 cases had homozygous or biallelic variants in DNAH5, DNAAF4 or DNAH11 classified as pathogenic or likely pathogenic. Three identified variants were deletions, ranging from 3 to 13 kb, for which WGS identified precise breakpoints, permitting confirmation by Sanger sequencing. WGS yielded identification of a de novo variant in a novel PCD gene TUBB4B. CONCLUSION: Here, WGS uplifted genetic diagnosis of PCD by identifying structural variants and novel modes of inheritance in new candidate genes. WGS could be an important component of the PCD diagnostic toolkit, increasing molecular diagnostic yield from current (70%) levels, and enhancing our understanding of fundamental biology of motile cilia and variants in the noncoding genome.

Ciliopathy patient variants reveal organelle-specific functions for TUBB4B in axonemal microtubules.

Tubulin, one of the most abundant cytoskeletal building blocks, has numerous isotypes in metazoans encoded by different conserved genes. Whether these distinct isotypes form cell type- and context-specific microtubule structures is poorly understood. Based on a cohort of 12 patients with primary ciliary dyskinesia as well as mouse mutants, we identified and characterized variants in the TUBB4B isotype that specifically perturbed centriole and cilium biogenesis. Distinct TUBB4B variants differentially affected microtubule dynamics and cilia formation in a dominant-negative manner. Structure-function studies revealed that different TUBB4B variants disrupted distinct tubulin interfaces, thereby enabling stratification of patients into three classes of ciliopathic diseases. These findings show that specific tubulin isotypes have distinct and nonredundant subcellular functions and establish a link between tubulinopathies and ciliopathies.

Accurate measurement of gene copy number for human alpha-defensin DEFA1A3.

BACKGROUND: Multi-allelic copy number variants include examples of extensive variation between individuals in the copy number of important genes, most notably genes involved in immune function. The definition of this variation, and analysis of its impact on function, has been hampered by the technical difficulty of large-scale but accurate typing of genomic copy number. The copy-variable alpha-defensin locus DEFA1A3 on human chromosome 8 commonly varies between 4 and 10 copies per diploid genome, and presents considerable challenges for accurate high-throughput typing. RESULTS: In this study, we developed two paralogue ratio tests and three allelic ratio measurements that, in combination, provide an accurate and scalable method for measurement of DEFA1A3 gene number. We combined information from different measurements in a maximum-likelihood framework which suggests that most samples can be assigned to an integer copy number with high confidence, and applied it to typing 589 unrelated European DNA samples. Typing the members of three-generation pedigrees provided further reassurance that correct integer copy numbers had been assigned. Our results have allowed us to discover that the SNP rs4300027 is strongly associated with DEFA1A3 gene copy number in European samples. CONCLUSIONS: We have developed an accurate and robust method for measurement of DEFA1A3 copy number. Interrogation of rs4300027 and associated SNPs in Genome-Wide Association Study SNP data provides no evidence that alpha-defensin copy number is a strong risk factor for phenotypes such as Crohn's disease, type I diabetes, HIV progression and multiple sclerosis.

Absolute measurement of the tissue origins of cell-free DNA in the healthy state and following paracetamol overdose.

BACKGROUND: Despite the emergence of cell-free DNA (cfDNA) as a clinical biomarker in cancer, the tissue origins of cfDNA in healthy individuals have to date been inferred only by indirect and relative measurement methods, such as tissue-specific methylation and nucleosomal profiling. METHODS: We performed the first direct, absolute measurement of the tissue origins of cfDNA, using tissue-specific knockout mouse strains, in both healthy mice and following paracetamol (APAP) overdose. We then investigated the utility of total cfDNA and the percentage of liver-specific cfDNA as clinical biomarkers in patients presenting with APAP overdose. RESULTS: Analysis of cfDNA from healthy tissue-specific knockout mice showed that cfDNA originates predominantly from white and red blood cell lineages, with minor contribution from hepatocytes, and no detectable contribution from skeletal and cardiac muscle. Following APAP overdose in mice, total plasma cfDNA and the percentage fraction originating from hepatocytes increased by ~ 100 and ~ 19-fold respectively. Total cfDNA increased by an average of more than 236-fold in clinical samples from APAP overdose patients with biochemical evidence of liver injury, and 18-fold in patients without biochemically apparent liver injury. Measurement of liver-specific cfDNA, using droplet digital PCR and methylation analysis, revealed that the contribution of liver to cfDNA was increased by an average of 175-fold in APAP overdose patients with biochemically apparent liver injury compared to healthy subjects, but was not increased in overdose patients with normal liver function tests. CONCLUSIONS: We present a novel method for measurement of the tissue origins of cfDNA in healthy and disease states and demonstrate the potential of cfDNA as a clinical biomarker in APAP overdose.

Genetic epidemiology of motor neuron disease-associated variants in the Scottish population.

Genetic understanding of motor neuron disease (MND) has evolved greatly in the past 10 years, including the recent identification of association between MND and variants in TBK1 and NEK1. Our aim was to determine the frequency of pathogenic variants in known MND genes and to assess whether variants in TBK1 and NEK1 contribute to the burden of MND in the Scottish population. SOD1, TARDBP, OPTN, TBK1, and NEK1 were sequenced in 441 cases and 400 controls. In addition to 44 cases known to carry a C9orf72 hexanucleotide repeat expansion, we identified 31 cases and 2 controls that carried a loss-of-function or pathogenic variant. Loss-of-function variants were found in TBK1 in 3 cases and no controls and, separately, in NEK1 in 3 cases and no controls. This study provides an accurate description of the genetic epidemiology of MND in Scotland and provides support for the contribution of both TBK1 and NEK1 to MND susceptibility in the Scottish population.