Shen Gu

Primary immunodeficiency diseases: Genomic approaches delineate heterogeneous Mendelian disorders

Background: Primary immunodeficiency diseases (PIDDs) are clinically and genetically heterogeneous disorders thus far associated with mutations in more than 300 genes. The clinical phenotypes derived from distinct genotypes can overlap. Genetic etiology can be a prognostic indicator of disease severity and can influence treatment decisions. Objective: We sought to investigate the ability of whole‐exome screening methods to detect disease‐causing variants in patients with PIDDs. Methods: Patients with PIDDs from 278 families from 22 countries were investigated by using whole‐exome sequencing. Computational copy number variant (CNV) prediction pipelines and an exome‐tiling chromosomal microarray were also applied to identify intragenic CNVs. Analytic approaches initially focused on 475 known or candidate PIDD genes but were nonexclusive and further tailored based on clinical data, family history, and immunophenotyping. Results: A likely molecular diagnosis was achieved in 110 (40%) unrelated probands. Clinical diagnosis was revised in about half (60/110) and management was directly altered in nearly a quarter (26/110) of families based on molecular findings. Twelve PIDD‐causing CNVs were detected, including 7 smaller than 30 Kb that would not have been detected with conventional diagnostic CNV arrays. Conclusion: This high‐throughput genomic approach enabled detection of disease‐related variants in unexpected genes; permitted detection of low‐grade constitutional, somatic, and revertant mosaicism; and provided evidence of a mutational burden in mixed PIDD immunophenotypes.

Recurrent De Novo and Biallelic Variation of ATAD3A, Encoding a Mitochondrial Membrane Protein, Results in Distinct Neurological Syndromes

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.

Alu-mediated diverse and complex pathogenic copy-number variants within human chromosome 17 at p13.3

Alu repetitive elements are known to be major contributors to genome instability by generating Alu-mediated copy-number variants (CNVs). Most of the reported Alu-mediated CNVs are simple deletions and duplications, and the mechanism underlying Alu-Alu-mediated rearrangement has been attributed to non-allelic homologous recombination (NAHR). Chromosome 17 at the p13.3 genomic region lacks extensive low-copy repeat architecture; however, it is highly enriched for Alu repetitive elements, with a fraction of 30% of total sequence annotated in the human reference genome, compared with the 10% genome-wide and 18% on chromosome 17. We conducted mechanistic studies of the 17p13.3 CNVs by performing high-density oligonucleotide array comparative genomic hybridization, specifically interrogating the 17p13.3 region with ∼150 bp per probe density; CNV breakpoint junctions were mapped to nucleotide resolution by polymerase chain reaction and Sanger sequencing. Studied rearrangements include 5 interstitial deletions, 14 tandem duplications, 7 terminal deletions and 13 complex genomic rearrangements (CGRs). Within the 17p13.3 region, Alu-Alu-mediated rearrangements were identified in 80% of the interstitial deletions, 46% of the tandem duplications and 50% of the CGRs, indicating that this mechanism was a major contributor for formation of breakpoint junctions. Our studies suggest that Alu repetitive elements facilitate formation of non-recurrent CNVs, CGRs and other structural aberrations of chromosome 17 at p13.3. The common observation of Alu-mediated rearrangement in CGRs and breakpoint junction sequences analysis further demonstrates that this type of mechanism is unlikely attributed to NAHR, but rather may be due to a recombination-coupled DNA replicative repair process.

Whole-exome sequencing in the molecular diagnosis of individuals with congenital anomalies of the kidney and urinary tract and identification of a new causative gene

Purpose:To investigate the utility of whole-exome sequencing (WES) to define a molecular diagnosis for patients clinically diagnosed with congenital anomalies of kidney and urinary tract (CAKUT).Methods:WES was performed in 62 families with CAKUT. WES data were analyzed for single-nucleotide variants (SNVs) in 35 known CAKUT genes, putatively deleterious sequence changes in new candidate genes, and potentially disease-associated copy-number variants (CNVs).Results:In approximately 5% of families, pathogenic SNVs were identified in PAX2, HNF1B, and EYA1. Observed phenotypes in these families expand the current understanding about the role of these genes in CAKUT. Four pathogenic CNVs were also identified using two CNV detection tools. In addition, we found one deleterious de novo SNV in FOXP1 among the 62 families with CAKUT. The clinical database of the Baylor Miraca Genetics laboratory was queried and seven additional unrelated individuals with novel de novo SNVs in FOXP1 were identified. Six of these eight individuals with FOXP1 SNVs have syndromic urinary tract defects, implicating this gene in urinary tract development.Conclusion:We conclude that WES can be used to identify molecular etiology (SNVs, CNVs) in a subset of individuals with CAKUT. WES can also help identify novel CAKUT genes.Genet Med 19 4, 412–420.

Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort

Abstract We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor–Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5% (82% for single-exonic and 88% for multi-exonic calls) and precision of 78% (53% single-exonic and 96% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17–50% of pathogenic CNVs in different disease cohorts where 7.1–11% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

Nonrecurrent 17p11.2p12 Rearrangement Events that Result in Two Concomitant Genomic Disorders: The PMP22-RAI1 Contiguous Gene Duplication Syndrome

The genomic duplication associated with Potocki-Lupski syndrome (PTLS) maps in close proximity to the duplication associated with Charcot-Marie-Tooth disease type 1A (CMT1A). PTLS is characterized by hypotonia, failure to thrive, reduced body weight, intellectual disability, and autistic features. CMT1A is a common autosomal dominant distal symmetric peripheral polyneuropathy. The key dosage-sensitive genes RAI1 and PMP22 are respectively associated with PTLS and CMT1A. Recurrent duplications accounting for the majority of subjects with these conditions are mediated by nonallelic homologous recombination between distinct low-copy repeat (LCR) substrates. The LCRs flanking a contiguous genomic interval encompassing both RAI1 and PMP22 do not share extensive homology; thus, duplications encompassing both loci are rare and potentially generated by a different mutational mechanism. We characterized genomic rearrangements that simultaneously duplicate PMP22 and RAI1, including nine potential complex genomic rearrangements, in 23 subjects by high-resolution array comparative genomic hybridization and breakpoint junction sequencing. Insertions and microhomologies were found at the breakpoint junctions, suggesting potential replicative mechanisms for rearrangement formation. At the breakpoint junctions of these nonrecurrent rearrangements, enrichment of repetitive DNA sequences was observed, indicating that they might predispose to genomic instability and rearrangement. Clinical evaluation revealed blended PTLS and CMT1A phenotypes with a potential earlier onset of neuropathy. Moreover, additional clinical findings might be observed due to the extra duplicated material included in the rearrangements. Our genomic analysis suggests replicative mechanisms as a predominant mechanism underlying PMP22-RAI1 contiguous gene duplications and provides further evidence supporting the role of complex genomic architecture in genomic instability.

Mechanisms for Complex Chromosomal Insertions

Chromosomal insertions are genomic rearrangements with a chromosome segment inserted into a non-homologous chromosome or a non-adjacent locus on the same chromosome or the other homologue, constituting ~2% of nonrecurrent copy-number gains. Little is known about the molecular mechanisms of their formation. We identified 16 individuals with complex insertions among 56,000 individuals tested at Baylor Genetics using clinical array comparative genomic hybridization (aCGH) and fluorescence in situ hybridization (FISH). Custom high-density aCGH was performed on 10 individuals with available DNA, and breakpoint junctions were fine-mapped at nucleotide resolution by long-range PCR and DNA sequencing in 6 individuals to glean insights into potential mechanisms of formation. We observed microhomologies and templated insertions at the breakpoint junctions, resembling the breakpoint junction signatures found in complex genomic rearrangements generated by replication-based mechanism(s) with iterative template switches. In addition, we analyzed 5 families with apparently balanced insertion in one parent detected by FISH analysis and found that 3 parents had additional small copy-number variants (CNVs) at one or both sides of the inserting fragments as well as at the inserted sites. We propose that replicative repair can result in interchromosomal complex insertions generated through chromothripsis-like chromoanasynthesis involving two or three chromosomes, and cause a significant fraction of apparently balanced insertions harboring small flanking CNVs.

Mechanisms for the Generation of Two Quadruplications Associated with Split-Hand Malformation

Germline copy‐number variants (CNVs) involving quadruplications are rare and the mechanisms generating them are largely unknown. Previously, we reported a 20‐week gestation fetus with split‐hand malformation; clinical microarray detected two maternally inherited triplications separated by a copy‐number neutral region at 17p13.3, involving BHLHA9 and part of YWHAE. Here, we describe an 18‐month‐old male sibling of the previously described fetus with split‐hand malformation. Custom high‐density microarray and digital droplet PCR revealed the copy‐number gains were actually quadruplications in the mother, the fetus, and her later born son. This quadruplication‐normal‐quadruplication pattern was shown to be expanded from the triplication‐normal‐triplication CNV at the same loci in the maternal grandmother. We mapped two breakpoint junctions and demonstrated that both are mediated by Alu repetitive elements and identical in these four individuals. We propose a three‐step process combining Alu‐mediated replicative‐repair‐based mechanism(s) and intergenerational, intrachromosomal nonallelic homologous recombination to generate the quadruplications in this family.

Hutterite-type cataract maps to chromosome 6p21.32-p21.31, cosegregates with a homozygous mutation in LEMD2, and is associated with sudden cardiac death

Juvenile‐onset cataracts are known among the Hutterites of North America. Despite being identified over 30 years ago, this autosomal recessive condition has not been mapped, and the disease gene is unknown.

Predicting human genes susceptible to genomic instability associated with Alu/Alu-mediated rearrangements

Alu elements, the short interspersed element numbering more than 1 million copies per human genome, can mediate the formation of copy number variants (CNVs) between substrate pairs. These Alu/Alu-mediated rearrangements (AAMRs) can result in pathogenic variants that cause diseases. To investigate the impact of AAMR on gene variation and human health, we first characterized Alus that are involved in mediating CNVs (CNV-Alus) and observed that these Alus tend to be evolutionarily younger. We then computationally generated, with the assistance of a supercomputer, a test data set consisting of 78 million Alu pairs and predicted ∼18% of them are potentially susceptible to AAMR. We further determined the relative risk of AAMR in 12,074 OMIM genes using the count of predicted CNV-Alu pairs and experimentally validated the predictions with 89 samples selected by correlating predicted hotspots with a database of CNVs identified by clinical chromosomal microarrays (CMAs) on the genomes of approximately 54,000 subjects. We fine-mapped 47 duplications, 40 deletions, and two complex rearrangements and examined a total of 52 breakpoint junctions of simple CNVs. Overall, 94% of the candidate breakpoints were at least partially Alu mediated. We successfully predicted all (100%) of Alu pairs that mediated deletions (n = 21) and achieved an 87% positive predictive value overall when including AAMR-generated deletions and duplications. We provided a tool, AluAluCNVpredictor, for assessing AAMR hotspots and their role in human disease. These results demonstrate the utility of our predictive model and provide insights into the genomic features and molecular mechanisms underlying AAMR.