Search

ABSTRACT Worldwide, one in five men aged over 50 years will experience osteoporosis or a clinical bone fracture, with a greater fracture‐related mortality rate than women. However, the genetic etiology of osteoporosis in men is still poorly understood. We aimed to identify the genetic variants and candidate genes associated with extremely low or high BMD for a better understanding of the biology underlying low bone density that may point to potential therapeutic targets for increasing bone mass. Subjects from the Osteoporotic Fractures in Men Study (MrOS) cohort were evaluated by age and BMI‐adjusted total hip BMD. Those with BMD values 3 SDs away from the mean were selected and the remaining individuals whose adjusted BMD ranked at the highest or lowest 100 were included. Men with the lowest adjusted BMD (N = 98) and highest adjusted BMD (N = 110) were chosen for exome sequencing. Controls (N = 82) were men of Northern and Western European descent from the US Utah population of the 1000 Genomes Project. Fisher's exact test was performed to compare low‐ or high‐BMD subjects with controls for single‐gene associations. Additionally, sets of candidate genes causative of heritable disorders of connective tissue, including osteogenesis imperfecta (OI) and Ehlers‐Danlos syndrome (EDS), were grouped for multigene and mutation burden analyses. No single‐gene associations with rare variants were found for either the low BMD group (33 genes) or high BMD group (18 genes). In the group of OI genes, we detected a significant threefold increased accumulation of rare variants in low‐BMD subjects compared with controls (p = 0.009). Additionally, genes associated with EDS had a twofold increased frequency in low‐BMD subjects compared with controls (p = 0.03). These findings reveal a rare variant burden in OI and EDS disease genes at low BMD, which suggests a potential gene‐panel approach to screen for multivariant associations in larger cohorts. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

Abstract Background Multiplexed Assays of Variant Effects (MAVEs) can test all possible single variants in a gene of interest. The resulting saturation-style functional data may help resolve variant classification disparities between populations, especially for Variants of Uncertain Significance (VUS). Methods We analyzed clinical significance classifications in 213,663 individuals of European-like genetic ancestry versus 206,975 individuals of non-European-like genetic ancestry from All of Us and the Genome Aggregation Database. Then, we incorporated clinically calibrated MAVE data into the Clinical Genome Resource’s Variant Curation Expert Panel rules to automate VUS reclassification for BRCA1, TP53, and PTEN. Results Using two orthogonal statistical approaches, we show a higher prevalence (p ≤ 5.95e − 06) of VUS in individuals of non-European-like genetic ancestry across all medical specialties assessed in all three databases. Further, in the non-European-like genetic ancestry group, higher rates of Benign or Likely Benign and variants with no clinical designation (p ≤ 2.5e − 05) were found across many medical specialties, whereas Pathogenic or Likely Pathogenic assignments were increased in individuals of European-like genetic ancestry (p ≤ 2.5e − 05). Using MAVE data, we reclassified VUS in individuals of non-European-like genetic ancestry at a significantly higher rate in comparison to reclassified VUS from European-like genetic ancestry (p = 9.1e − 03) effectively compensating for the VUS disparity. Further, essential code analysis showed equitable impact of MAVE evidence codes but inequitable impact of allele frequency (p = 7.47e − 06) and computational predictor (p = 6.92e − 05) evidence codes for individuals of non-European-like genetic ancestry. Conclusions Generation of saturation-style MAVE data should be a priority to reduce VUS disparities and produce equitable training data for future computational predictors.

We identified a homozygous missense alteration (c.75C>A, p.D25E) in CLCC1, encoding a presumptive intracellular chloride channel highly expressed in the retina, associated with autosomal recessive retinitis pigmentosa (arRP) in eight consanguineous families of Pakistani descent. The p.D25E alteration decreased CLCC1 channel function accompanied by accumulation of mutant protein in granules within the ER lumen, while siRNA knockdown of CLCC1 mRNA induced apoptosis in cultured ARPE-19 cells. TALEN KO in zebrafish was lethal 11 days post fertilization. The depressed electroretinogram (ERG) cone response and cone spectral sensitivity of 5 dpf KO zebrafish and reduced eye size, retinal thickness, and expression of rod and cone opsins could be rescued by injection of wild type CLCC1 mRNA. Clcc1+/- KO mice showed decreased ERGs and photoreceptor number. Together these results strongly suggest that intracellular chloride transport by CLCC1 is a critical process in maintaining retinal integrity, and CLCC1 is crucial for survival and function of retinal cells.

Dominant mutations in CACNA1A, encoding the α-1A subunit of the neuronal P/Q type voltage-dependent Ca2+ channel, can cause diverse neurological phenotypes. Rare cases of markedly severe early onset developmental delay and congenital ataxia can be due to de novo CACNA1A missense alleles, with variants affecting the S4 transmembrane segments of the channel, some of which are reported to be loss-of-function. Exome sequencing in five individuals with severe early onset ataxia identified one novel variant (p.R1673P), in a girl with global developmental delay and progressive cerebellar atrophy, and a recurrent, de novo p.R1664Q variant, in four individuals with global developmental delay, hypotonia, and ophthalmologic abnormalities. Given the severity of these phenotypes we explored their functional impact in Drosophila. We previously generated null and partial loss-of-function alleles of cac, the homolog of CACNA1A in Drosophila. Here, we created transgenic wild type and mutant genomic rescue constructs with the two noted conserved point mutations. The p.R1673P mutant failed to rescue cac lethality, displayed a gain-of-function phenotype in electroretinograms (ERG) recorded from mutant clones, and evolved a neurodegenerative phenotype in aging flies, based on ERGs and transmission electron microscopy. In contrast, the p.R1664Q variant exhibited loss of function and failed to develop a neurodegenerative phenotype. Hence, the novel R1673P allele produces neurodegenerative phenotypes in flies and human, likely due to a toxic gain of function.

Abstract Background X-linked acrogigantism (X-LAG; MIM: 300942) is a severe form of pituitary gigantism caused by chromosome Xq26.3 duplications involving GPR101. X-LAG-associated duplications disrupt the integrity of the topologically associating domain (TAD) containing GPR101 and lead to the formation of a neo-TAD that drives pituitary GPR101 misexpression and gigantism. As X-LAG is fully penetrant and heritable, duplications involving GPR101 identified on prenatal screening studies, like amniocentesis, can pose an interpretation challenge for medical geneticists and raise important concerns for patients and families. Therefore, providing robust information on the functional genomic impact of such duplications has important research and clinical value with respect to gene regulation and triplosensitivity traits. Methods We employed 4C/HiC-seq as a clinical tool to determine the functional impact of incidentally discovered GPR101 duplications on TAD integrity in three families. After defining duplications and breakpoints around GPR101 by clinical-grade and high-density aCGH, we constructed 4C/HiC chromatin contact maps for our study population and compared them with normal and active (X-LAG) controls. Results We showed that duplications involving GPR101 that preserved the centromeric invariant TAD boundary did not generate a pathogenic neo-TAD and that ectopic enhancers were not adopted. This allowed us to discount presumptive/suspected X-LAG diagnoses and GPR101 misexpression, obviating the need for intensive clinical follow-up. Conclusions This study highlights the importance of TAD boundaries and chromatin interactions in determining the functional impact of copy number variants and provides proof-of-concept for using 4C/HiC-seq as a clinical tool to acquire crucial information for genetic counseling and to support clinical decision-making in cases of suspected TADopathies.

Abstract Developmental and epileptic encephalopathies (DEEs) feature altered brain development, developmental delay and seizures, with seizures exacerbating developmental delay. Here we identify a cohort with biallelic variants in DENND5A, encoding a membrane trafficking protein, and develop animal models with phenotypes like the human syndrome. We demonstrate that DENND5A interacts with Pals1/MUPP1, components of the Crumbs apical polarity complex required for symmetrical division of neural progenitor cells. Human induced pluripotent stem cells lacking DENND5A fail to undergo symmetric cell division with an inherent propensity to differentiate into neurons. These phenotypes result from misalignment of the mitotic spindle in apical neural progenitors. Cells lacking DENND5A orient away from the proliferative apical domain surrounding the ventricles, biasing daughter cells towards a more fate-committed state, ultimately shortening the period of neurogenesis. This study provides a mechanism for DENND5A-related DEE that may be generalizable to other developmental conditions and provides variant-specific clinical information for physicians and families.

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.

Many loci in the human genome harbor complex genomic structures that can result in susceptibility to genomic rearrangements leading to various genomic disorders. Nephronophthisis 1 (NPHP1, MIM# 256100) is an autosomal recessive disorder that can be caused by defects of NPHP1; the gene maps within the human 2q13 region where low copy repeats (LCRs) are abundant. Loss of function of NPHP1 is responsible for approximately 85% of the NPHP1 cases-about 80% of such individuals carry a large recurrent homozygous NPHP1 deletion that occurs via nonallelic homologous recombination (NAHR) between two flanking directly oriented ~45 kb LCRs. Published data revealed a non-pathogenic inversion polymorphism involving the NPHP1 gene flanked by two inverted ~358 kb LCRs. Using optical mapping and array-comparative genomic hybridization, we identified three potential novel structural variant (SV) haplotypes at the NPHP1 locus that may protect a haploid genome from the NPHP1 deletion. Inter-species comparative genomic analyses among primate genomes revealed massive genomic changes during evolution. The aggregated data suggest that dynamic genomic rearrangements occurred historically within the NPHP1 locus and generated SV haplotypes observed in the human population today, which may confer differential susceptibility to genomic instability and the NPHP1 deletion within a personal genome. Our study documents diverse SV haplotypes at a complex LCR-laden human genomic region. Comparative analyses provide a model for how this complex region arose during primate evolution, and studies among humans suggest that intra-species polymorphism may potentially modulate an individual's susceptibility to acquiring disease-associated alleles.

Abstract Background We previously described the KINSSHIP syndrome, an autosomal dominant disorder associated with intellectual disability (ID), mesomelic dysplasia and horseshoe kidney, caused by de novo variants in the degron of AFF3. Mouse knock-ins and overexpression in zebrafish provided evidence for a dominant-negative mode of action, wherein an increased level of AFF3 resulted in pathological effects. Methods Evolutionary constraints suggest that other modes-of-inheritance could be at play. We challenged this hypothesis by screening ID cohorts for individuals with predicted-to-be damaging variants in AFF3. We used both animal and cellular models to assess the deleteriousness of the identified variants. Results We identified an individual with a KINSSHIP-like phenotype carrying a de novo partial duplication of AFF3 further strengthening the hypothesis that an increased level of AFF3 is pathological. We also detected seventeen individuals displaying a milder syndrome with either heterozygous Loss-of-Function (LoF) or biallelic missense variants in AFF3. Consistent with semi-dominance, we discovered three patients with homozygous LoF and one compound heterozygote for a LoF and a missense variant, who presented more severe phenotypes than their heterozygous parents. Matching zebrafish knockdowns exhibit neurological defects that could be rescued by expressing human AFF3 mRNA, confirming their association with the ablation of aff3. Conversely, some of the human AFF3 mRNAs carrying missense variants identified in affected individuals did not rescue these phenotypes. Overexpression of mutated AFF3 mRNAs in zebrafish embryos produced a significant increase of abnormal larvae compared to wild-type overexpression further demonstrating deleteriousness. To further assess the effect of AFF3 variation, we profiled the transcriptome of fibroblasts from affected individuals and engineered isogenic cells harboring + / + , KINSSHIP/KINSSHIP, LoF/ + , LoF/LoF or KINSSHIP/LoF AFF3 genotypes. The expression of more than a third of the AFF3 bound loci is modified in either the KINSSHIP/KINSSHIP or the LoF/LoF lines. While the same pathways are affected, only about one third of the differentially expressed genes are common to the homozygote datasets, indicating that AFF3 LoF and KINSSHIP variants largely modulate transcriptomes differently, e.g. the DNA repair pathway displayed opposite modulation. Conclusions Our results and the high pleiotropy shown by variation at this locus suggest that minute changes in AFF3 function are deleterious.

Abstract Background Multilocus pathogenic variants (MPVs) are genetic changes that affect multiple gene loci or regions of the genome, collectively leading to multiple molecular diagnoses. MPVs may also contribute to intrafamilial phenotypic variability between affected individuals within a nuclear family. In this study, we aim to gain further insights into the influence of MPVs on a disease manifestation in individual research subjects and explore the complexities of the human genome within a familial context. Methods We conducted a systematic reanalysis of exome sequencing data and runs of homozygosity (ROH) regions of 47 sibling pairs previously diagnosed with various neurodevelopmental disorders (NDD). Results We found siblings with MPVs driven by long ROH regions in 8.5% of families (4/47). The patients with MPVs exhibited significantly higher FROH values (p-value = 1.4e-2) and larger total ROH length (p-value = 1.8e-2). Long ROH regions mainly contribute to this pattern; the siblings with MPVs have a larger total size of long ROH regions than their siblings in all families (p-value = 6.9e-3). Whereas the short ROH regions in the siblings with MPVs are lower in total size compared to their sibling pairs with single locus pathogenic variants (p-value = 0.029), and there are no statistically significant differences in medium ROH regions between sibling pairs (p-value = 0.52). Conclusion This study sheds light on the significance of considering MPVs in families with affected sibling pairs and the role of ROH as an adjuvant tool in explaining clinical variability within families. Identifying individuals carrying MPVs may have implications for disease management, identification of possible disease risks to different family members, genetic counseling and exploring personalized treatment approaches.

Abstract Background NODAL signaling plays a critical role in embryonic patterning and heart development in vertebrates. Genetic variants resulting in perturbations of the TGF-β/NODAL signaling pathway have reproducibly been shown to cause laterality defects in humans. To further explore this association and improve genetic diagnosis, the study aims to identify and characterize a broader range of NODAL variants in a large number of individuals with laterality defects. Methods We re-analyzed a cohort of 321 proband-only exomes of individuals with clinically diagnosed laterality congenital heart disease (CHD) using family-based, rare variant genomic analyses. To this cohort we added 12 affected subjects with known NODAL variants and CHD from institutional research and clinical cohorts to investigate an allelic series. For those with candidate contributory variants, variant allele confirmation and segregation analysis were studied by Sanger sequencing in available family members. Array comparative genomic hybridization and droplet digital PCR were utilized for copy number variants (CNV) validation and characterization. We performed Human Phenotype Ontology (HPO)-based quantitative phenotypic analyses to dissect allele-specific phenotypic differences. Results Missense, nonsense, splice site, indels, and/or structural variants of NODAL were identified as potential causes of heterotaxy and other laterality defects in 33 CHD cases. We describe a recurrent complex indel variant for which the nucleic acid secondary structure predictions implicate secondary structure mutagenesis as a possible mechanism for formation. We identified two CNV deletion alleles spanning NODAL in two unrelated CHD cases. Furthermore, 17 CHD individuals were found (16/17 with known Hispanic ancestry) to have the c.778G > A:p.G260R NODAL missense variant which we propose reclassification from variant of uncertain significance (VUS) to likely pathogenic. Quantitative HPO-based analyses of the observed clinical phenotype for all cases with p.G260R variation, including heterozygous, homozygous, and compound heterozygous cases, reveal clustering of individuals with biallelic variation. This finding provides evidence for a genotypic-phenotypic correlation and an allele-specific gene dosage model. Conclusion Our data further support a role for rare deleterious variants in NODAL as a cause for sporadic human laterality defects, expand the repertoire of observed anatomical complexity of potential cardiovascular anomalies, and implicate an allele specific gene dosage model.

Inverted repeats (IRs) can facilitate structural variation as crucibles of genomic rearrangement. Complex duplication-inverted triplication-duplication (DUP-TRP/INV-DUP) rearrangements that contain breakpoint junctions within IRs have been recently associated with both MECP2 duplication syndrome (MIM#300260) and Pelizaeus-Merzbacher disease (PMD, MIM#312080). We investigated 17 unrelated PMD subjects with copy number gains at the PLP1 locus including triplication and quadruplication of specific genomic intervals-16/17 were found to have a DUP-TRP/INV-DUP rearrangement product. An IR distal to PLP1 facilitates DUP-TRP/INV-DUP formation as well as an inversion structural variation found frequently amongst normal individuals. We show that a homology-or homeology-driven replicative mechanism of DNA repair can apparently mediate template switches within stretches of microhomology. Moreover, we provide evidence that quadruplication and potentially higher order amplification of a genomic interval can occur in a manner consistent with rolling circle amplification as predicted by the microhomology-mediated break induced replication (MMBIR) model.

Megacystis-microcolon-intestinal hypoperistalsis syndrome (MMIHS) is a rare disorder of enteric smooth muscle function affecting the intestine and bladder. Patients with this severe phenotype are dependent on total parenteral nutrition and urinary catheterization. The cause of this syndrome has remained a mystery since Berdon's initial description in 1976. No genes have been clearly linked to MMIHS. We used whole-exome sequencing for gene discovery followed by targeted Sanger sequencing in a cohort of patients with MMIHS and intestinal pseudo-obstruction. We identified heterozygous ACTG2 missense variants in 15 unrelated subjects, ten being apparent de novo mutations. Ten unique variants were detected, of which six affected CpG dinucleotides and resulted in missense mutations at arginine residues, perhaps related to biased usage of CpG containing codons within actin genes. We also found some of the same heterozygous mutations that we observed as apparent de novo mutations in MMIHS segregating in families with intestinal pseudo-obstruction, suggesting that ACTG2 is responsible for a spectrum of smooth muscle disease. ACTG2 encodes γ2 enteric actin and is the first gene to be clearly associated with MMIHS, suggesting an important role for contractile proteins in enteric smooth muscle disease.

Normal development of the genitourinary (GU) tract is a complex process that frequently goes awry. In male children the most frequent congenital GU anomalies are cryptorchidism (1-4%), hypospadias (1%) and micropenis (0.35%). Bladder exstrophy and epispadias complex (BEEC) (1∶47000) occurs less frequently but significantly impacts patients' lives. Array comparative genomic hybridization (aCGH) identified seven individuals with overlapping deletions in the 2p15 region (66.0 kb-5.6 Mb). Six of these patients have GU defects, while the remaining patient has no GU defect. These deletions encompass the transcription factor OTX1. Subjects 2-7 had large de novo CNVs (2.39-6.31 Mb) and exhibited features similar to those associated with the 2p15p16.1 and 2p15p14 microdeletion syndromes, including developmental delay, short stature, and variable GU defects. Subject-1 with BEEC had the smallest deletion (66 kb), which deleted only one copy of OTX1. Otx1-null mice have seizures, prepubescent transient growth retardation and gonadal defects. Two subjects have short stature, two have seizures, and six have GU defects, mainly affecting the external genitalia. The presence of GU defects in six patients in our cohort and eight of thirteen patients reported with deletions within 2p14p16.1 (two with deletion of OTX1) suggest that genes in 2p15 are important for GU development. Genitalia defects in these patients could result from the effect of OTX1 on pituitary hormone secretion or on the regulation of SHH signaling, which is crucial for development of the bladder and genitalia.

IntroductionSystemic lupus erythematosus is a multi-faceted autoimmune disorder of complex etiology. Pre-pubertal onset of pediatric systemic lupus erythematosus (pSLE) is uncommon and should raise suspicion for a genetic driver of disease. Autosomal recessive p40phox deficiency is a rare immunologic disorder characterized by defective but not abolished NADPH oxidase activity with residual production of reactive oxygen species (ROS) by phagocytic cells.Case presentationWe report the case of a now 18-year-old female with pSLE onset at 7 years of age. She presented with recurrent fever and malar rash. Aspects of her immune dysregulation over time have included typical pSLE features including production of autoantibodies, hematologic manifestations, and hypocomplementemia, as well as chronic suppurative skin lesions and recurrent infections. Genetic analysis revealed biallelic pathogenic variants in NCF4 resulting in p40phox deficiency. Comprehensive NADPH oxidase activity studies confirmed significantly decreased production of reactive oxygen species, confirming the cellular phenotype seen in p40phox deficient patients.ConclusionsHere, we present a patient with pSLE harboring biallelic variants in NCF4. Our patient represents a unique clinical presentation of severe onset autoimmunity in the setting of a rare inborn error of immunity affecting NADPH oxidase activity. This case underscores the need to consider genetic causes of pSLE in cases of pre-pubertal onset and atypical disease.

Abstract CELSR3 codes for a planar cell polarity protein. We describe twelve affected individuals from eleven independent families with bi-allelic variants in CELSR3. Affected individuals presented with an overlapping phenotypic spectrum comprising central nervous system (CNS) anomalies (7/12), combined CNS anomalies and congenital anomalies of the kidneys and urinary tract (CAKUT) (3/12) and CAKUT only (2/12). Computational simulation of the 3D protein structure suggests the position of the identified variants to be implicated in penetrance and phenotype expression. CELSR3 immunolocalization in human embryonic urinary tract and transient suppression and rescue experiments of Celsr3 in fluorescent zebrafish reporter lines further support an embryonic role of CELSR3 in CNS and urinary tract formation.

Genomic disorders are often caused by recurrent copy number variations (CNVs), with nonallelic homologous recombination (NAHR) as the underlying mechanism. Recently, several microhomology-mediated repair mechanisms--such as microhomology-mediated end-joining (MMEJ), fork stalling and template switching (FoSTeS), microhomology-mediated break-induced replication (MMBIR), serial replication slippage (SRS), and break-induced SRS (BISRS)--were described in the etiology of non-recurrent CNVs in human disease. In addition, their formation may be stimulated by genomic architectural features. It is, however, largely unexplored to what extent these mechanisms contribute to rare, locus-specific pathogenic CNVs. Here, fine-mapping of 42 microdeletions of the FOXL2 locus, encompassing FOXL2 (32) or its regulatory domain (10), serves as a model for rare, locus-specific CNVs implicated in genetic disease. These deletions lead to blepharophimosis syndrome (BPES), a developmental condition affecting the eyelids and the ovary. For breakpoint mapping we used targeted array-based comparative genomic hybridization (aCGH), quantitative PCR (qPCR), long-range PCR, and Sanger sequencing of the junction products. Microhomology, ranging from 1 bp to 66 bp, was found in 91.7% of 24 characterized breakpoint junctions, being significantly enriched in comparison with a random control sample. Our results show that microhomology-mediated repair mechanisms underlie at least 50% of these microdeletions. Moreover, genomic architectural features, like sequence motifs, non-B DNA conformations, and repetitive elements, were found in all breakpoint regions. In conclusion, the majority of these microdeletions result from microhomology-mediated mechanisms like MMEJ, FoSTeS, MMBIR, SRS, or BISRS. Moreover, we hypothesize that the genomic architecture might drive their formation by increasing the susceptibility for DNA breakage or promote replication fork stalling. Finally, our locus-centered study, elucidating the etiology of a large set of rare microdeletions involved in a monogenic disorder, can serve as a model for other clustered, non-recurrent microdeletions in genetic disease.

Curation and interpretation of copy number variants identified by genome-wide testing is challenged by the large number of events harbored in each personal genome. Conventional determination of phenotypic relevance relies on patterns of higher frequency in affected individuals versus controls; however, an increasing amount of ascertained variation is rare or private to clans. Consequently, frequency data have less utility to resolve pathogenic from benign. One solution is disease-specific algorithms that leverage gene knowledge together with variant frequency to aid prioritization. We used large-scale resources including Gene Ontology, protein-protein interactions and other annotation systems together with a broad set of 83 genes with known associations to epilepsy to construct a pathogenicity score for the phenotype. We evaluated the score for all annotated human genes and applied Bayesian methods to combine the derived pathogenicity score with frequency information from our diagnostic laboratory. Analysis determined Bayes factors and posterior distributions for each gene. We applied our method to subjects with abnormal chromosomal microarray results and confirmed epilepsy diagnoses gathered by electronic medical record review. Genes deleted in our subjects with epilepsy had significantly higher pathogenicity scores and Bayes factors compared to subjects referred for non-neurologic indications. We also applied our scores to identify a recently validated epilepsy gene in a complex genomic region and to reveal candidate genes for epilepsy. We propose a potential use in clinical decision support for our results in the context of genome-wide screening. Our approach demonstrates the utility of integrative data in medical genomics.

The hotspots of structural polymorphisms and structural mutability in the human genome remain to be explained mechanistically. We examine associations of structural mutability with germline DNA methylation and with non-allelic homologous recombination (NAHR) mediated by low-copy repeats (LCRs). Combined evidence from four human sperm methylome maps, human genome evolution, structural polymorphisms in the human population, and previous genomic and disease studies consistently points to a strong association of germline hypomethylation and genomic instability. Specifically, methylation deserts, the ~1% fraction of the human genome with the lowest methylation in the germline, show a tenfold enrichment for structural rearrangements that occurred in the human genome since the branching of chimpanzee and are highly enriched for fast-evolving loci that regulate tissue-specific gene expression. Analysis of copy number variants (CNVs) from 400 human samples identified using a custom-designed array comparative genomic hybridization (aCGH) chip, combined with publicly available structural variation data, indicates that association of structural mutability with germline hypomethylation is comparable in magnitude to the association of structural mutability with LCR-mediated NAHR. Moreover, rare CNVs occurring in the genomes of individuals diagnosed with schizophrenia, bipolar disorder, and developmental delay and de novo CNVs occurring in those diagnosed with autism are significantly more concentrated within hypomethylated regions. These findings suggest a new connection between the epigenome, selective mutability, evolution, and human disease.

The functional contribution of CNV to human biology and disease pathophysiology has undergone limited exploration. Recent observations in humans indicate a tentative link between CNV and weight regulation. Smith-Magenis syndrome (SMS), manifesting obesity and hypercholesterolemia, results from a deletion CNV at 17p11.2, but is sometimes due to haploinsufficiency of a single gene, RAI1. The reciprocal duplication in 17p11.2 causes Potocki-Lupski syndrome (PTLS). We previously constructed mouse strains with a deletion, Df(11)17, or duplication, Dp(11)17, of the mouse genomic interval syntenic to the SMS/PTLS region. We demonstrate that Dp(11)17 is obesity-opposing; it conveys a highly penetrant, strain-independent phenotype of reduced weight, leaner body composition, lower TC/LDL, and increased insulin sensitivity that is not due to alteration in food intake or activity level. When fed with a high-fat diet, Dp(11)17/+ mice display much less weight gain and metabolic change than WT mice, demonstrating that the Dp(11)17 CNV protects against metabolic syndrome. Reciprocally, Df(11)17/+ mice with the deletion CNV have increased weight, higher fat content, decreased HDL, and reduced insulin sensitivity, manifesting a bona fide metabolic syndrome. These observations in the deficiency animal model are supported by human data from 76 SMS subjects. Further, studies on knockout/transgenic mice showed that the metabolic consequences of Dp(11)17 and Df(11)17 CNVs are not only due to dosage alterations of Rai1, the predominant dosage-sensitive gene for SMS and likely also PTLS. Our experiments in chromosome-engineered mouse CNV models for human genomic disorders demonstrate that a CNV can be causative for weight/metabolic phenotypes. Furthermore, we explored the biology underlying the contribution of CNV to the physiology of weight control and energy metabolism. The high penetrance, strain independence, and resistance to dietary influences associated with the CNVs in this study are features distinct from most SNP-associated metabolic traits and further highlight the potential importance of CNV in the etiology of both obesity and MetS as well as in the protection from these traits.

A novel microduplication syndrome involving various-sized contiguous duplications in 17p13.3 has recently been described, suggesting that increased copy number of genes in 17p13.3, particularly PAFAH1B1, is associated with clinical features including facial dysmorphism, developmental delay, and autism spectrum disorder. We have previously shown that patient-derived cell lines from individuals with haploinsufficiency of RPA1, a gene within 17p13.3, exhibit an impaired ATR-dependent DNA damage response (DDR). Here, we show that cell lines from patients with duplications specifically incorporating RPA1 exhibit a different although characteristic spectrum of DDR defects including abnormal S phase distribution, attenuated DNA double strand break (DSB)-induced RAD51 chromatin retention, elevated genomic instability, and increased sensitivity to DNA damaging agents. Using controlled conditional over-expression of RPA1 in a human model cell system, we also see attenuated DSB-induced RAD51 chromatin retention. Furthermore, we find that transient over-expression of RPA1 can impact on homologous recombination (HR) pathways following DSB formation, favouring engagement in aberrant forms of recombination and repair. Our data identifies unanticipated defects in the DDR associated with duplications in 17p13.3 in humans involving modest RPA1 over-expression.

Myelin protein zero (MPZ) is a critical structural component of myelin in the peripheral nervous system. The MPZ gene is regulated, in part, by the transcription factors SOX10 and EGR2. Mutations in MPZ, SOX10, and EGR2 have been implicated in demyelinating peripheral neuropathies, suggesting that components of this transcriptional network are candidates for harboring disease-causing mutations (or otherwise functional variants) that affect MPZ expression.We utilized a combination of multi-species sequence comparisons, transcription factor-binding site predictions, targeted human DNA re-sequencing, and in vitro and in vivo enhancer assays to study human non-coding MPZ variants.Our efforts revealed a variant within the first intron of MPZ that resides within a previously described SOX10 binding site is associated with decreased enhancer activity, and alters binding of nuclear proteins. Additionally, the genomic segment harboring this variant directs tissue-relevant reporter gene expression in zebrafish.This is the first reported MPZ variant within a cis-acting transcriptional regulatory element. While we were unable to implicate this variant in disease onset, our data suggests that similar non-coding sequences should be screened for mutations in patients with neurological disease. Furthermore, our multi-faceted approach for examining the functional significance of non-coding variants can be readily generalized to study other loci important for myelin structure and function.

A large fraction of genome variation between individuals is comprised of submicroscopic copy number variation of genomic DNA segments. We assessed the relative contribution of structural changes and gene dosage alterations on phenotypic outcomes with mouse models of Smith-Magenis and Potocki-Lupski syndromes. We phenotyped mice with 1n (Deletion/+), 2n (+/+), 3n (Duplication/+), and balanced 2n compound heterozygous (Deletion/Duplication) copies of the same region. Parallel to the observations made in humans, such variation in gene copy number was sufficient to generate phenotypic consequences: in a number of cases diametrically opposing phenotypes were associated with gain versus loss of gene content. Surprisingly, some neurobehavioral traits were not rescued by restoration of the normal gene copy number. Transcriptome profiling showed that a highly significant propensity of transcriptional changes map to the engineered interval in the five assessed tissues. A statistically significant overrepresentation of the genes mapping to the entire length of the engineered chromosome was also found in the top-ranked differentially expressed genes in the mice containing rearranged chromosomes, regardless of the nature of the rearrangement, an observation robust across different cell lineages of the central nervous system. Our data indicate that a structural change at a given position of the human genome may affect not only locus and adjacent gene expression but also "genome regulation." Furthermore, structural change can cause the same perturbation in particular pathways regardless of gene dosage. Thus, the presence of a genomic structural change, as well as gene dosage imbalance, contributes to the ultimate phenotype.

An accurate and precisely annotated genome assembly is a fundamental requirement for functional genomic analysis. Here, the complete DNA sequence and gene annotation of mouse Chromosome 11 was used to test the efficacy of large-scale sequencing for mutation identification. We re-sequenced the 14,000 annotated exons and boundaries from over 900 genes in 41 recessive mutant mouse lines that were isolated in an N-ethyl-N-nitrosourea (ENU) mutation screen targeted to mouse Chromosome 11. Fifty-nine sequence variants were identified in 55 genes from 31 mutant lines. 39% of the lesions lie in coding sequences and create primarily missense mutations. The other 61% lie in noncoding regions, many of them in highly conserved sequences. A lesion in the perinatal lethal line l11Jus13 alters a consensus splice site of nucleoredoxin (Nxn), inserting 10 amino acids into the resulting protein. We conclude that point mutations can be accurately and sensitively recovered by large-scale sequencing, and that conserved noncoding regions should be included for disease mutation identification. Only seven of the candidate genes we report have been previously targeted by mutation in mice or rats, showing that despite ongoing efforts to functionally annotate genes in the mammalian genome, an enormous gap remains between phenotype and function. Our data show that the classical positional mapping approach of disease mutation identification can be extended to large target regions using high-throughput sequencing.

Chromosome structural changes with nonrecurrent endpoints associated with genomic disorders offer windows into the mechanism of origin of copy number variation (CNV). A recent report of nonrecurrent duplications associated with Pelizaeus-Merzbacher disease identified three distinctive characteristics. First, the majority of events can be seen to be complex, showing discontinuous duplications mixed with deletions, inverted duplications, and triplications. Second, junctions at endpoints show microhomology of 2-5 base pairs (bp). Third, endpoints occur near pre-existing low copy repeats (LCRs). Using these observations and evidence from DNA repair in other organisms, we derive a model of microhomology-mediated break-induced replication (MMBIR) for the origin of CNV and, ultimately, of LCRs. We propose that breakage of replication forks in stressed cells that are deficient in homologous recombination induces an aberrant repair process with features of break-induced replication (BIR). Under these circumstances, single-strand 3' tails from broken replication forks will anneal with microhomology on any single-stranded DNA nearby, priming low-processivity polymerization with multiple template switches generating complex rearrangements, and eventual re-establishment of processive replication.

The modern synthetic view of human evolution proposes that the fixation of novel mutations is driven by the balance among selective advantage, selective disadvantage, and genetic drift. When considering the global architecture of the human genome, the same model can be applied to understanding the rapid acquisition and proliferation of exogenous DNA. To explore the evolutionary forces that might have morphed human genome architecture, we investigated the origin, composition, and functional potential of numts (nuclear mitochondrial pseudogenes), partial copies of the mitochondrial genome found abundantly in chromosomal DNA. Our data indicate that these elements are unlikely to be advantageous, since they possess no gross positional, transcriptional, or translational features that might indicate beneficial functionality subsequent to integration. Using sequence analysis and fossil dating, we also show a probable burst of integration of numts in the primate lineage that centers on the prosimian-anthropoid split, mimics closely the temporal distribution of Alu and processed pseudogene acquisition, and coincides with the major climatic change at the Paleocene-Eocene boundary. We therefore propose a model according to which the gross architecture and repeat distribution of the human genome can be largely accounted for by a population bottleneck early in the anthropoid lineage and subsequent effectively neutral fixation of repetitive DNA, rather than positive selection or unusual insertion pressures.

BACKGROUND: Array Comparative Genomic Hybridization (a-CGH) is a powerful molecular cytogenetic tool to detect genomic imbalances and study disease mechanism and pathogenesis. We report our experience with the clinical implementation of this high resolution human genome analysis, referred to as Chromosomal Microarray Analysis (CMA). METHODS AND FINDINGS: CMA was performed clinically on 2513 postnatal samples from patients referred with a variety of clinical phenotypes. The initial 775 samples were studied using CMA array version 4 and the remaining 1738 samples were analyzed with CMA version 5 containing expanded genomic coverage. Overall, CMA identified clinically relevant genomic imbalances in 8.5% of patients: 7.6% using V4 and 8.9% using V5. Among 117 cases referred for additional investigation of a known cytogenetically detectable rearrangement, CMA identified the majority (92.5%) of the genomic imbalances. Importantly, abnormal CMA findings were observed in 5.2% of patients (98/1872) with normal karyotypes/FISH results, and V5, with expanded genomic coverage, enabled a higher detection rate in this category than V4. For cases without cytogenetic results available, 8.0% (42/524) abnormal CMA results were detected; again, V5 demonstrated an increased ability to detect abnormality. Improved diagnostic potential of CMA is illustrated by 90 cases identified with 51 cryptic microdeletions and 39 predicted apparent reciprocal microduplications in 13 specific chromosomal regions associated with 11 known genomic disorders. In addition, CMA identified copy number variations (CNVs) of uncertain significance in 262 probands; however, parental studies usually facilitated clinical interpretation. Of these, 217 were interpreted as familial variants and 11 were determined to be de novo; the remaining 34 await parental studies to resolve the clinical significance. CONCLUSIONS: This large set of clinical results demonstrates the significantly improved sensitivity of CMA for the detection of clinically relevant genomic imbalances and highlights the need for comprehensive genetic counseling to facilitate accurate clinical correlation and interpretation.

Rearrangements of our genome can be responsible for inherited as well as sporadic traits. The analyses of chromosome breakpoints in the proximal short arm of Chromosome 17 (17p) reveal nonallelic homologous recombination (NAHR) as a major mechanism for recurrent rearrangements whereas nonhomologous end-joining (NHEJ) can be responsible for many of the nonrecurrent rearrangements. Genome architectural features consisting of low-copy repeats (LCRs), or segmental duplications, can stimulate and mediate NAHR, and there are hotspots for the crossovers within the LCRs. Rearrangements introduce variation into our genome for selection to act upon and as such serve an evolutionary function analogous to base pair changes. Genomic rearrangements may cause Mendelian diseases, produce complex traits such as behaviors, or represent benign polymorphic changes. The mechanisms by which rearrangements convey phenotypes are diverse and include gene dosage, gene interruption, generation of a fusion gene, position effects, unmasking of recessive coding region mutations (single nucleotide polymorphisms, SNPs, in coding DNA) or other functional SNPs, and perhaps by effects on transvection.

Purpose: The variome of the Turkish (TK) population, a population with a considerable history of admixture and consanguinity, has not been deeply investigated for insights on the genomic architecture of disease. Methods: We generated and analyzed a database of variants derived from exome sequencing data of 773 TK unrelated, clinically affected individuals with various suspected Mendelian disease traits and 643 unaffected relatives. Results: Using uniform manifold approximation and projection, we showed that the TK genomes are more similar to those of Europeans and consist of 2 main subpopulations: clusters 1 and 2 (N = 235 and 1181, respectively), which differ in admixture proportion and variome (https://turkishvariomedb.shinyapps.io/tvdb/). Furthermore, the higher inbreeding coefficient values observed in the TK affected compared with unaffected individuals correlated with a larger median span of long-sized (>2.64 Mb) runs of homozygosity (ROH) regions (P value = 2.09e-18). We show that long-sized ROHs are more likely to be formed on recently configured haplotypes enriched for rare homozygous deleterious variants in the TK affected compared with TK unaffected individuals (P value = 3.35e-11). Analysis of genotype-phenotype correlations reveals that genes with rare homozygous deleterious variants in long-sized ROHs provide the most comprehensive set of molecular diagnoses for the observed disease traits with a systematic quantitative analysis of Human Phenotype Ontology terms. Conclusion: Our findings support the notion that novel rare variants on newly configured haplotypes arising within the recent past generations of a family or clan contribute significantly to recessive disease traits in the TK population.

Abstract Craniofacial microsomia (CFM; also known as Goldenhar syndrome), is a craniofacial developmental disorder of variable expressivity and severity with a recognizable set of abnormalities. These birth defects are associated with structures derived from the first and second pharyngeal arches, can occur unilaterally and include ear dysplasia, microtia, preauricular tags and pits, facial asymmetry and other malformations. The inheritance pattern is controversial, and the molecular etiology of this syndrome is largely unknown. A total of 670 patients belonging to unrelated pedigrees with European and Chinese ancestry with CFM, are investigated. We identify 18 likely pathogenic variants in 21 probands (3.1%) in FOXI3. Biochemical experiments on transcriptional activity and subcellular localization of the likely pathogenic FOXI3 variants, and knock-in mouse studies strongly support the involvement of FOXI3 in CFM. Our findings indicate autosomal dominant inheritance with reduced penetrance, and/or autosomal recessive inheritance. The phenotypic expression of the FOXI3 variants is variable. The penetrance of the likely pathogenic variants in the seemingly dominant form is reduced, since a considerable number of such variants in affected individuals were inherited from non-affected parents. Here we provide suggestive evidence that common variation in the FOXI3 allele in trans with the pathogenic variant could modify the phenotypic severity and accounts for the incomplete penetrance.

Summary: Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

NK cell deficiencies (NKD) are a type of primary immune deficiency in which the major immunologic abnormality affects NK cell number, maturity, or function. Since NK cells contribute to immune defense against virally infected cells, patients with NKD experience higher susceptibility to chronic, recurrent, and fatal viral infections. An individual with recurrent viral infections and mild hypogammaglobulinemia was identified to have an X-linked damaging variant in the transcription factor gene ELF4. The variant does not decrease expression but disrupts ELF4 protein interactions and DNA binding, reducing transcriptional activation of target genes and selectively impairing ELF4 function. Corroborating previous murine models of ELF4 deficiency (Elf4–/–) and using a knockdown human NK cell line, we determined that ELF4 is necessary for normal NK cell development, terminal maturation, and function. Through characterization of the NK cells of the proband, expression of the proband’s variant in Elf4–/– mouse hematopoietic precursor cells, and a human in vitro NK cell maturation model, we established this ELF4 variant as a potentially novel cause of NKD.

Human NK cell deficiency (NKD) is a primary immunodeficiency in which the main clinically relevant immunological defect involves missing or dysfunctional NK cells. Here, we describe a familial NKD case in which 2 siblings had a substantive NKD and neutropenia in the absence of other immune system abnormalities. Exome sequencing identified compound heterozygous variants in Go-Ichi-Ni-San (GINS) complex subunit 4 (GINS4, also known as SLD5), an essential component of the human replicative helicase, which we demonstrate to have a damaging impact upon the expression and assembly of the GINS complex. Cells derived from affected individuals and a GINS4-knockdown cell line demonstrate delayed cell cycle progression, without signs of improper DNA synthesis or increased replication stress. By modeling partial GINS4 depletion in differentiating NK cells in vitro, we demonstrate the causal relationship between the genotype and the NK cell phenotype, as well as a cell-intrinsic defect in NK cell development. Thus, biallelic partial loss-of-function mutations in GINS4 define a potentially novel disease-causing gene underlying NKD with neutropenia. Together with the previously described mutations in other helicase genes causing NKD, and with the mild defects observed in other human cells, these variants underscore the importance of this pathway in NK cell biology.