The existence of pleiotropy in disorders with multi-organ involvement can suggest therapeutic targets that could ameliorate overall disease severity. Here we assessed pleiotropy of modifier genes in ...cystic fibrosis (CF). CF, caused by mutations in the cystic fibrosis transmembrane conductance regulator (
CFTR)
gene, affects the lungs, liver, pancreas and intestines. However, modifier genes contribute to variable disease severity across affected organs, even in individuals with the same
CFTR
genotype. We sought to determine whether
SLC26A9
,
SLC9A3
and
SLC6A14,
that contribute to meconium ileus in CF, are pleiotropic for other early-affecting CF co-morbidities. In the Canadian CF population, we assessed evidence for pleiotropic effects on (1) pediatric lung disease severity (
n
= 815), (2) age at first acquisition of
Pseudomonas aeruginosa
(
P. aeruginosa
) (
n
= 730), and (3) prenatal pancreatic damage measured by immunoreactive trypsinogen (
n
= 126). A multiple-phenotype analytic strategy assessed evidence for pleiotropy in the presence of phenotypic correlation. We required the same alleles to be associated with detrimental effects.
SLC26A9
was pleiotropic for meconium ileus and pancreatic damage (
p
= 0.002 at rs7512462),
SLC9A3
for meconium ileus and lung disease (
p
= 1.5 × 10
−6
at rs17563161), and
SLC6A14
for meconium ileus and both lung disease and age at first
P. aeruginosa
infection (
p
= 0.0002 and
p
= 0.006 at rs3788766, respectively). The meconium ileus risk alleles in
SLC26A9
,
SLC9A3
and
SLC6A14
are pleiotropic, increasing risk for other early CF co-morbidities. Furthermore, co-morbidities affecting the same organ tended to associate with the same genes. The existence of pleiotropy within this single disorder suggests that complementary therapeutic strategies to augment solute transport will benefit multiple CF-associated tissues.
ATPase family AAA-domain containing protein 3A (ATAD3A) is a nuclear-encoded mitochondrial membrane protein implicated in mitochondrial dynamics, nucleoid organization, protein translation, cell ...growth, and cholesterol metabolism. We identified a recurrent de novo ATAD3A c.1582C>T (p.Arg528Trp) variant by whole-exome sequencing (WES) in five unrelated individuals with a core phenotype of global developmental delay, hypotonia, optic atrophy, axonal neuropathy, and hypertrophic cardiomyopathy. We also describe two families with biallelic variants in ATAD3A, including a homozygous variant in two siblings, and biallelic ATAD3A deletions mediated by nonallelic homologous recombination (NAHR) between ATAD3A and gene family members ATAD3B and ATAD3C. Tissue-specific overexpression of borR534W, the Drosophila mutation homologous to the human c.1582C>T (p.Arg528Trp) variant, resulted in a dramatic decrease in mitochondrial content, aberrant mitochondrial morphology, and increased autophagy. Homozygous null bor larvae showed a significant decrease of mitochondria, while overexpression of borWT resulted in larger, elongated mitochondria. Finally, fibroblasts of an affected individual exhibited increased mitophagy. We conclude that the p.Arg528Trp variant functions through a dominant-negative mechanism that results in small mitochondria that trigger mitophagy, resulting in a reduction in mitochondrial content. ATAD3A variation represents an additional link between mitochondrial dynamics and recognizable neurological syndromes, as seen with MFN2, OPA1, DNM1L, and STAT2 mutations.
Given the rarity of most single-gene Mendelian disorders, concerted efforts of data exchange between clinical and scientific communities are critical to optimize molecular diagnosis and novel disease ...gene discovery.
We designed and implemented protocols for the study of cases for which a plausible molecular diagnosis was not achieved in a clinical genomics diagnostic laboratory (i.e. unsolved clinical exomes). Such cases were recruited to a research laboratory for further analyses, in order to potentially: (1) accelerate novel disease gene discovery; (2) increase the molecular diagnostic yield of whole exome sequencing (WES); and (3) gain insight into the genetic mechanisms of disease. Pilot project data included 74 families, consisting mostly of parent-offspring trios. Analyses performed on a research basis employed both WES from additional family members and complementary bioinformatics approaches and protocols.
Analysis of all possible modes of Mendelian inheritance, focusing on both single nucleotide variants (SNV) and copy number variant (CNV) alleles, yielded a likely contributory variant in 36% (27/74) of cases. If one includes candidate genes with variants identified within a single family, a potential contributory variant was identified in a total of ~51% (38/74) of cases enrolled in this pilot study. The molecular diagnosis was achieved in 30/63 trios (47.6%). Besides this, the analysis workflow yielded evidence for pathogenic variants in disease-associated genes in 4/6 singleton cases (66.6%), 1/1 multiplex family involving three affected siblings, and 3/4 (75%) quartet families. Both the analytical pipeline and the collaborative efforts between the diagnostic and research laboratories provided insights that allowed recent disease gene discoveries (PURA, TANGO2, EMC1, GNB5, ATAD3A, and MIPEP) and increased the number of novel genes, defined in this study as genes identified in more than one family (DHX30 and EBF3).
An efficient genomics pipeline in which clinical sequencing in a diagnostic laboratory is followed by the detailed reanalysis of unsolved cases in a research environment, supplemented with WES data from additional family members, and subject to adjuvant bioinformatics analyses including relaxed variant filtering parameters in informatics pipelines, can enhance the molecular diagnostic yield and provide mechanistic insights into Mendelian disorders. Implementing these approaches requires collaborative clinical molecular diagnostic and research efforts.
The underlying genetic etiology of rhabdomyolysis remains elusive in a significant fraction of individuals presenting with recurrent metabolic crises and muscle weakness. Using exome sequencing, we ...identified bi-allelic mutations in TANGO2 encoding transport and Golgi organization 2 homolog (Drosophila) in 12 subjects with episodic rhabdomyolysis, hypoglycemia, hyperammonemia, and susceptibility to life-threatening cardiac tachyarrhythmias. A recurrent homozygous c.460G>A (p.Gly154Arg) mutation was found in four unrelated individuals of Hispanic/Latino origin, and a homozygous ∼34 kb deletion affecting exons 3–9 was observed in two families of European ancestry. One individual of mixed Hispanic/European descent was found to be compound heterozygous for c.460G>A (p.Gly154Arg) and the deletion of exons 3–9. Additionally, a homozygous exons 4–6 deletion was identified in a consanguineous Middle Eastern Arab family. No homozygotes have been reported for these changes in control databases. Fibroblasts derived from a subject with the recurrent c.460G>A (p.Gly154Arg) mutation showed evidence of increased endoplasmic reticulum stress and a reduction in Golgi volume density in comparison to control. Our results show that the c.460G>A (p.Gly154Arg) mutation and the exons 3–9 heterozygous deletion in TANGO2 are recurrent pathogenic alleles present in the Latino/Hispanic and European populations, respectively, causing considerable morbidity in the homozygotes in these populations.
Improved resolution of molecular diagnostic technologies enabled detection of smaller sized exonic level copy-number variants (CNVs). The contribution of CNVs to autosomal recessive (AR) conditions ...may be better recognized using a large clinical cohort.
We retrospectively investigated the CNVs’ contribution to AR conditions in cases subjected to chromosomal microarray analysis (CMA, N=~70,000) and/or clinical exome sequencing (ES, N=~12,000) at Baylor Genetics; most had pediatric onset neurodevelopmental disorders.
CNVs contributed to biallelic variations in 87 cases, including 81 singletons and three affected sibling pairs. Seventy cases had CNVs affecting both alleles, and 17 had a CNV and a single-nucleotide variant (SNV)/indel in trans. In total, 94.3% of AR-CNVs affected one gene; among these 41.4% were single-exon and 35.0% were multiexon partial-gene events. Sixty-nine percent of homozygous AR-CNVs were embedded in homozygous genomic intervals. Five cases had large deletions unmasking an SNV/indel on the intact allele for a recessive condition, resulting in multiple molecular diagnoses.
AR-CNVs are often smaller in size, transmitted through generations, and underrecognized due to limitations in clinical CNV detection methods. Our findings from a large clinical cohort emphasized integrated CNV and SNV/indel analyses for precise clinical and molecular diagnosis especially in the context of genomic disorders.
Using a multistage genetic association approach comprising 7,480 affected individuals and 7,779 controls, we identified markers in chromosomal region 8q24 associated with colorectal cancer. In stage ...1, we genotyped 99,632 SNPs in 1,257 affected individuals and 1,336 controls from Ontario. In stages 2-4, we performed serial replication studies using 4,024 affected individuals and 4,042 controls from Seattle, Newfoundland and Scotland. We identified one locus on chromosome 8q24 and another on 9p24 having combined odds ratios (OR) for stages 1-4 of 1.18 (trend; P = 1.41 × 10−8) and 1.14 (trend; P = 1.32 × 10−5), respectively. Additional analyses in 2,199 affected individuals and 2,401 controls from France and Europe supported the association at the 8q24 locus (OR = 1.16, trend; 95% confidence interval (c.i.): 1.07-1.26; P = 5.05 × 10−4). A summary across all seven studies at the 8q24 locus was highly significant (OR = 1.17, c.i.: 1.12-1.23; P = 3.16 × 10−11). This locus has also been implicated in prostate cancer.
Metastasis is responsible for 90% of human cancer mortality, yet it remains a challenge to model human cancer metastasis in vivo. Here we describe mouse models of high-grade serous ovarian cancer, ...also known as high-grade serous carcinoma (HGSC), the most common and deadliest human ovarian cancer type. Mice genetically engineered to harbor Dicer1 and Pten inactivation and mutant p53 robustly replicate the peritoneal metastases of human HGSC with complete penetrance. Arising from the fallopian tube, tumors spread to the ovary and metastasize throughout the pelvic and peritoneal cavities, invariably inducing hemorrhagic ascites. Widespread and abundant peritoneal metastases ultimately cause mouse deaths (100%). Besides the phenotypic and histopathological similarities, mouse HGSCs also display marked chromosomal instability, impaired DNA repair, and chemosensitivity. Faithfully recapitulating the clinical metastases as well as molecular and genomic features of human HGSC, this murine model will be valuable for elucidating the mechanisms underlying the development and progression of metastatic ovarian cancer and also for evaluating potential therapies.
Sufficiently powered case-control studies with next-generation sequence (NGS) data remain prohibitively expensive for many investigators. If feasible, a more efficient strategy would be to include ...publicly available sequenced controls. However, these studies can be confounded by differences in sequencing platform; alignment, single nucleotide polymorphism and variant calling algorithms; read depth; and selection thresholds. Assuming one can match cases and controls on the basis of ethnicity and other potential confounding factors, and one has access to the aligned reads in both groups, we investigate the effect of systematic differences in read depth and selection threshold when comparing allele frequencies between cases and controls. We propose a novel likelihood-based method, the robust variance score (RVS), that substitutes genotype calls by their expected values given observed sequence data.
We show theoretically that the RVS eliminates read depth bias in the estimation of minor allele frequency. We also demonstrate that, using simulated and real NGS data, the RVS method controls Type I error and has comparable power to the 'gold standard' analysis with the true underlying genotypes for both common and rare variants.
An RVS R script and instructions can be found at strug.research.sickkids.ca, and at https://github.com/strug-lab/RVS.
lisa.strug@utoronto.ca
Supplementary data are available at Bioinformatics online.