A multi-ancestry GWAS of Fuchs corneal dystrophy highlights the contributions of laminins, collagen, and endothelial cell regulation.

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular etiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program followed by multi-ancestry meta-analysis with the previous largest FECD GWAS, for a total of 3970 cases and 333,794 controls. We confirm the previous four loci, and identify eight novel loci: SSBP3, THSD7A, LAMB1, PIDD1, RORA, HS3ST3B1, LAMA5, and COL18A1. We further confirm the TCF4 locus in GWAS for admixed African and Hispanic/Latino ancestries and show an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling, validated through homology, suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and colocalization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Autoimmune alleles at the major histocompatibility locus modify melanoma susceptibility.

Autoimmunity and cancer represent two different aspects of immune dysfunction. Autoimmunity is characterized by breakdowns in immune self-tolerance, while impaired immune surveillance can allow for tumorigenesis. The class I major histocompatibility complex (MHC-I), which displays derivatives of the cellular peptidome for immune surveillance by CD8+ T cells, serves as a common genetic link between these conditions. As melanoma-specific CD8+ T cells have been shown to target melanocyte-specific peptide antigens more often than melanoma-specific antigens, we investigated whether vitiligo- and psoriasis-predisposing MHC-I alleles conferred a melanoma-protective effect. In individuals with cutaneous melanoma from both The Cancer Genome Atlas (n = 451) and an independent validation set (n = 586), MHC-I autoimmune-allele carrier status was significantly associated with a later age of melanoma diagnosis. Furthermore, MHC-I autoimmune-allele carriers were significantly associated with decreased risk of developing melanoma in the Million Veteran Program (OR = 0.962, p = 0.024). Existing melanoma polygenic risk scores (PRSs) did not predict autoimmune-allele carrier status, suggesting these alleles provide orthogonal risk-relevant information. Mechanisms of autoimmune protection were neither associated with improved melanoma-driver mutation association nor improved gene-level conserved antigen presentation relative to common alleles. However, autoimmune alleles showed higher affinity relative to common alleles for particular windows of melanocyte-conserved antigens and loss of heterozygosity of autoimmune alleles caused the greatest reduction in presentation for several conserved antigens across individuals with loss of HLA alleles. Overall, this study presents evidence that MHC-I autoimmune-risk alleles modulate melanoma risk unaccounted for by current PRSs.

Multi-ancestry GWAS of Fuchs corneal dystrophy highlights roles of laminins, collagen, and endothelial cell regulation.

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular pathophysiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program (MVP) and meta-analyzed with the previous largest FECD GWAS, finding twelve significant loci (eight novel). We further confirmed the TCF4 locus in admixed African and Hispanic/Latino ancestries, and found an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and co-localization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Pharmacogenetic allele variant frequencies: An analysis of the VA's Million Veteran Program (MVP) as a representation of the diversity in US population.

We present allele frequencies of pharmacogenomics relevant variants across multiple ancestry in a sample representative of the US population. We analyzed 658,582 individuals with genotype data and extracted pharmacogenomics relevant single nucleotide variant (SNV) alleles, human leukocyte antigens (HLA) 4-digit alleles and an important copy number variant (CNV), the full deletion/duplication of CYP2D6. We compiled distinct allele frequency tables for European, African American, Hispanic, and Asian ancestry individuals. In addition, we compiled allele frequencies based on local ancestry reconstruction in the African-American (2-way deconvolution) and Hispanic (3-way deconvolution) cohorts.

Congenital X-linked neutropenia with myelodysplasia and somatic tetraploidy due to a germline mutation in SEPT6.

Septins play key roles in mammalian cell division and cytokinesis but have not previously been implicated in a germline human disorder. A male infant with severe neutropenia and progressive dysmyelopoiesis with tetraploid myeloid precursors was identified. No known genetic etiologies for neutropenia or bone marrow failure were found. However, next-generation sequencing of germline samples from the patient revealed a novel, de novo germline stop-loss mutation in the X-linked gene SEPT6 that resulted in reduced SEPT6 staining in bone marrow granulocyte precursors and megakaryocytes. Patient skin fibroblast-derived induced pluripotent stem cells (iPSCs) produced reduced myeloid colonies, particularly of the granulocyte lineage. CRISPR/Cas9 knock-in of the patient's mutation or complete knock-out of SEPT6 was not tolerated in non-patient-derived iPSCs or human myeloid cell lines, but SEPT6 knock-out was successful in an erythroid cell line and resulting clones revealed a propensity to multinucleation. In silico analysis predicts that the mutated protein hinders the dimerization of SEPT6 coiled-coils in both parallel and antiparallel arrangements, which could in turn impair filament formation. These data demonstrate a critical role for SEPT6 in chromosomal segregation in myeloid progenitors that can account for the unusual predisposition to aneuploidy and dysmyelopoiesis.

Homozygous deletions implicate non-coding epigenetic marks in Autism spectrum disorder.

More than 98% of the human genome is made up of non-coding DNA, but techniques to ascertain its contribution to human disease have lagged far behind our understanding of protein coding variations. Autism spectrum disorder (ASD) has been mostly associated with coding variations via de novo single nucleotide variants (SNVs), recessive/homozygous SNVs, or de novo copy number variants (CNVs); however, most ASD cases continue to lack a genetic diagnosis. We analyzed 187 consanguineous ASD families for biallelic CNVs. Recessive deletions were significantly enriched in affected individuals relative to their unaffected siblings (17% versus 4%, p < 0.001). Only a small subset of biallelic deletions were predicted to result in coding exon disruption. In contrast, biallelic deletions in individuals with ASD were enriched for overlap with regulatory regions, with 23/28 CNVs disrupting histone peaks in ENCODE (p < 0.009). Overlap with regulatory regions was further demonstrated by comparisons to the 127-epigenome dataset released by the Roadmap Epigenomics project, with enrichment for enhancers found in primary brain tissue and neuronal progenitor cells. Our results suggest a novel noncoding mechanism of ASD, describe a powerful method to identify important noncoding regions in the human genome, and emphasize the potential significance of gene activation and regulation in cognitive and social function.

Congenital macrothrombocytopenia with focal myelofibrosis due to mutations in human G6b-B is rescued in humanized mice.

Unlike primary myelofibrosis (PMF) in adults, myelofibrosis in children is rare. Congenital (inherited) forms of myelofibrosis (cMF) have been described, but the underlying genetic mechanisms remain elusive. Here we describe 4 families with autosomal recessive inherited macrothrombocytopenia with focal myelofibrosis due to germ line loss-of-function mutations in the megakaryocyte-specific immunoreceptor tyrosine-based inhibitory motif (ITIM)-containing receptor G6b-B (G6b, C6orf25, or MPIG6B). Patients presented with a mild-to-moderate bleeding diathesis, macrothrombocytopenia, anemia, leukocytosis and atypical megakaryocytes associated with a distinctive, focal, perimegakaryocytic pattern of bone marrow fibrosis. In addition to identifying the responsible gene, the description of G6b-B as the mutated protein potentially implicates aberrant G6b-B megakaryocytic signaling and activation in the pathogenesis of myelofibrosis. Targeted insertion of human G6b in mice rescued the knockout phenotype and a copy number effect of human G6b-B expression was observed. Homozygous knockin mice expressed 25% of human G6b-B and exhibited a marginal reduction in platelet count and mild alterations in platelet function; these phenotypes were more severe in heterozygous mice that expressed only 12% of human G6b-B. This study establishes G6b-B as a critical regulator of platelet homeostasis in humans and mice. In addition, the humanized G6b mouse will provide an invaluable tool for further investigating the physiological functions of human G6b-B as well as testing the efficacy of drugs targeting this receptor.

Biallelic mutations in human DCC cause developmental split-brain syndrome.

Motor, sensory, and integrative activities of the brain are coordinated by a series of midline-bridging neuronal commissures whose development is tightly regulated. Here we report a new human syndrome in which these commissures are widely disrupted, thus causing clinical manifestations of horizontal gaze palsy, scoliosis, and intellectual disability. Affected individuals were found to possess biallelic loss-of-function mutations in the gene encoding the axon-guidance receptor 'deleted in colorectal carcinoma' (DCC), which has been implicated in congenital mirror movements when it is mutated in the heterozygous state but whose biallelic loss-of-function human phenotype has not been reported. Structural MRI and diffusion tractography demonstrated broad disorganization of white-matter tracts throughout the human central nervous system (CNS), including loss of all commissural tracts at multiple levels of the neuraxis. Combined with data from animal models, these findings show that DCC is a master regulator of midline crossing and development of white-matter projections throughout the human CNS.

Ringed sideroblasts in ß-thalassemia.

Symptomatic ß-thalassemia is one of the globally most common inherited disorders. The initial clinical presentation is variable. Although common hematological analyses are typically sufficient to diagnose the disease, sometimes the diagnosis can be more challenging. We describe a series of patients with ß-thalassemia whose diagnosis was delayed, required bone marrow examination in one affected member of each family, and revealed ringed sideroblasts, highlighting the association of this morphological finding with these disorders. Thus, in the absence of characteristic congenital sideroblastic mutations or causes of acquired sideroblastic anemia, the presence of ringed sideroblasts should raise the suspicion of ß-thalassemia.

Hyperammonemia as a Presenting Feature in Two Siblings with FBXL4 Variants.

Early-onset mitochondrial encephalomyopathy is a rare disorder that presents in the neonatal period with lactic acidosis, hypotonia, and developmental delay. Sequence variants in the nuclear-encoded gene FBXL4 have been previously demonstrated to be a cause of early-onset mitochondrial encephalomyopathy in several unrelated families. We have identified a pair of siblings with mutations in FBXL4 who each presented in the neonatal period with hyperammonemia, low plasma levels of aspartate, low urine levels of tricarboxylic acid cycle intermediates suggesting a defect in anaplerosis, and cerebellar hypoplasia in addition to lactic acidosis and other classic signs of mitochondrial encephalomyopathy. After initial clinical stabilization, both subjects continued to have episodic exacerbations characterized by lactic acidosis and hyperammonemia. Previously reported cases of FBXL4 mutations are reviewed and compared with these affected siblings. These two new cases add to the spectrum of disease caused by mutations in FBLX4 and suggest possible benefit from anaplerotic therapies.

A recurring mutation in the respiratory complex 1 protein NDUFB11 is responsible for a novel form of X-linked sideroblastic anemia.

The congenital sideroblastic anemias (CSAs) are a heterogeneous group of inherited blood disorders characterized by pathological mitochondrial iron deposition in erythroid precursors. Each known cause has been attributed to a mutation in a protein associated with heme biosynthesis, iron-sulfur cluster biogenesis, mitochondrial translation, or a component of the mitochondrial respiratory chain. Here, we describe a recurring mutation, c.276_278del, p.F93del, in NDUFB11, a mitochondrial respiratory complex I-associated protein encoded on the X chromosome, in 5 males with a variably syndromic, normocytic CSA. The p.F93del mutation results in respiratory insufficiency and loss of complex I stability and activity in patient-derived fibroblasts. Targeted introduction of this allele into K562 erythroleukemia cells results in a proliferation defect with minimal effect on erythroid differentiation potential, suggesting the mechanism of anemia in this disorder.

Mutations in the substrate binding glycine-rich loop of the mitochondrial processing peptidase-α protein (PMPCA) cause a severe mitochondrial disease.

We describe a large Lebanese family with two affected members, a young female proband and her male cousin, who had multisystem involvement including profound global developmental delay, severe hypotonia and weakness, respiratory insufficiency, blindness, and lactic acidemia-findings consistent with an underlying mitochondrial disorder. Whole-exome sequencing was performed on DNA from the proband and both parents. The proband and her cousin carried compound heterozygous mutations in the PMPCA gene that encodes for α-mitochondrial processing peptidase (α-MPP), a protein likely involved in the processing of mitochondrial proteins. The variants were located close to and postulated to affect the substrate binding glycine-rich loop of the α-MPP protein. Functional assays including immunofluorescence and western blot analysis on patient's fibroblasts revealed that these variants reduced α-MPP levels and impaired frataxin production and processing. We further determined that those defects could be rescued through the expression of exogenous wild-type PMPCA cDNA. Our findings link defective α-MPP protein to a severe mitochondrial disease.

Evolution of GOUNDRY, a cryptic subgroup of Anopheles gambiae s.l., and its impact on susceptibility to Plasmodium infection.

The recent discovery of a previously unknown genetic subgroup of Anopheles gambiae sensu lato underscores our incomplete understanding of complexities of vector population demographics in Anopheles. This subgroup, named GOUNDRY, does not rest indoors as adults and is highly susceptible to Plasmodium infection in the laboratory. Initial description of GOUNDRY suggested it differed from other known Anopheles taxa in surprising and sometimes contradictory ways, raising a number of questions about its age, population size and relationship to known subgroups. To address these questions, we sequenced the complete genomes of 12 wild-caught GOUNDRY specimens and compared these genomes to a panel of Anopheles genomes. We show that GOUNDRY is most closely related to Anopheles coluzzii, and the timing of cladogenesis is not recent, substantially predating the advent of agriculture. We find a large region of the X chromosome that has swept to fixation in GOUNDRY within the last 100 years, which may be an inversion that serves as a partial barrier to contemporary gene flow. Interestingly, we show that GOUNDRY has a history of inbreeding that is significantly associated with susceptibility to Plasmodium infection in the laboratory. Our results illuminate the genomic evolution of one of probably several cryptic, ecologically specialized subgroups of Anopheles and provide a potent example of how vector population dynamics may complicate efforts to control or eradicate malaria.

Genetic Structure of a Local Population of the Anopheles gambiae Complex in Burkina Faso.

Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. Population heterogeneities for ecological and behavioral attributes expand and stabilize malaria transmission over space and time, and populations may change in response to vector control, urbanization and other factors. There is a need for approaches to comprehensively describe the structure and characteristics of a sympatric local mosquito population, because incomplete knowledge of vector population composition may hinder control efforts. To this end, we used a genome-wide custom SNP typing array to analyze a population collection from a single geographic region in West Africa. The combination of sample depth (n = 456) and marker density (n = 1536) unambiguously resolved population subgroups, which were also compared for their relative susceptibility to natural genotypes of Plasmodium falciparum malaria. The population subgroups display fluctuating patterns of differentiation or sharing across the genome. Analysis of linkage disequilibrium identified 19 new candidate genes for association with underlying population divergence between sister taxa, A. coluzzii (M-form) and A. gambiae (S-form).

LARS2 Variants Associated with Hydrops, Lactic Acidosis, Sideroblastic Anemia, and Multisystem Failure.

Pathogenic variants in mitochondrial aminoacyl-tRNA synthetases result in a broad range of mitochondrial respiratory chain disorders despite their shared role in mitochondrial protein synthesis. LARS2 encodes the mitochondrial leucyl-tRNA synthetase, which attaches leucine to its cognate tRNA. Sequence variants in LARS2 have previously been associated with Perrault syndrome, characterized by premature ovarian failure and hearing loss (OMIM #615300). In this study, we report variants in LARS2 that are associated with a severe multisystem metabolic disorder. The proband was born prematurely with severe lactic acidosis, hydrops, and sideroblastic anemia. She had multisystem complications with hyaline membrane disease, impaired cardiac function, a coagulopathy, pulmonary hypertension, and progressive renal disease and succumbed at 5 days of age. Whole exome sequencing of patient DNA revealed compound heterozygous variants in LARS2 (c.1289C>T; p.Ala430Val and c.1565C>A; p.Thr522Asn). The c.1565C>A (p.Thr522Asn) LARS2 variant has previously been associated with Perrault syndrome and both identified variants are predicted to be damaging (SIFT, PolyPhen). Muscle and liver samples from the proband did not display marked mitochondrial respiratory chain enzyme deficiency. Immunoblotting of patient muscle and liver showed LARS2 levels were reduced in liver and complex I protein levels were reduced in patient muscle and liver. Aminoacylation assays revealed p.Ala430Val LARS2 had an 18-fold loss of catalytic efficiency and p.Thr522Asn a 9-fold loss compared to wild-type LARS2. We suggest that the identified LARS2 variants are responsible for the severe multisystem clinical phenotype seen in this baby and that mutations in LARS2 can result in variable phenotypes.

Congenital sideroblastic anemia due to mutations in the mitochondrial HSP70 homologue HSPA9.

The congenital sideroblastic anemias (CSAs) are relatively uncommon diseases characterized by defects in mitochondrial heme synthesis, iron-sulfur (Fe-S) cluster biogenesis, or protein synthesis. Here we demonstrate that mutations in HSPA9, a mitochondrial HSP70 homolog located in the chromosome 5q deletion syndrome 5q33 critical deletion interval and involved in mitochondrial Fe-S biogenesis, result in CSA inherited as an autosomal recessive trait. In a fraction of patients with just 1 severe loss-of-function allele, expression of the clinical phenotype is associated with a common coding single nucleotide polymorphism in trans that correlates with reduced messenger RNA expression and results in a pseudodominant pattern of inheritance.

Association mapping by pooled sequencing identifies TOLL 11 as a protective factor against Plasmodium falciparum in Anopheles gambiae.

BACKGROUND: The genome-wide association study (GWAS) techniques that have been used for genetic mapping in other organisms have not been successfully applied to mosquitoes, which have genetic characteristics of high nucleotide diversity, low linkage disequilibrium, and complex population stratification that render population-based GWAS essentially unfeasible at realistic sample size and marker density. METHODS: We designed a novel mapping strategy for the mosquito system that combines the power of linkage mapping with the resolution afforded by genetic association. We established founder colonies from West Africa, controlled for diversity, linkage disequilibrium and population stratification. Colonies were challenged by feeding on the infectious stage of the human malaria parasite, Plasmodium falciparum, mosquitoes were phenotyped for parasite load, and DNA pools for phenotypically similar mosquitoes were Illumina sequenced. Phenotype-genotype mapping was carried out in two stages, coarse and fine. RESULTS: In the first mapping stage, pooled sequences were analysed genome-wide for intervals displaying relativereduction in diversity between phenotype pools, and candidate genomic loci were identified for influence upon parasite infection levels. In the second mapping stage, focused genotyping of SNPs from the first mapping stage was carried out in unpooled individual mosquitoes and replicates. The second stage confirmed significant SNPs in a locus encoding two Toll-family proteins. RNAi-mediated gene silencing and infection challenge revealed that TOLL 11 protects mosquitoes against P. falciparum infection. CONCLUSIONS: We present an efficient and cost-effective method for genetic mapping using natural variation segregating in defined recent Anopheles founder colonies, and demonstrate its applicability for mapping in a complex non-model genome. This approach is a practical and preferred alternative to population-based GWAS for first-pass mapping of phenotypes in Anopheles. This design should facilitate mapping of other traits involved in physiology, epidemiology, and behaviour.

The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: genetic correlation and functional testing.

BACKGROUND: Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. It is known that a large haplotype shared between An. gambiae and Anopheles coluzzii by introgression carries point mutations of the voltage-gated sodium channel gene para, including the L1014F kdr mutation associated with insensitivity to pyrethroid insecticides. Carriage of L1014F kdr is also correlated with higher susceptibility to infection with Plasmodium falciparum. However, the genetic mechanism and causative gene(s) underlying the parasite susceptibility phenotype are not known. METHODS: Mosquitoes from the wild Burkina Faso population were challenged by feeding on natural P. falciparum gametocytes. Oocyst infection phenotypes were determined and were tested for association with SNP genotypes. Candidate genes in the detected locus were prioritized and RNAi-mediated gene silencing was used to functionally test for gene effects on P. falciparum susceptibility. RESULTS: A genetic locus, Pfin6, was identified that influences infection levels of P. falciparum in mosquitoes. The locus segregates as a ~3 Mb haplotype carrying 65 predicted genes including the para gene. The haplotype carrying the kdr allele of para is linked to increased parasite infection prevalence, but many single nucleotide polymorphisms on the haplotype are also equally linked to the infection phenotype. Candidate genes in the haplotype were prioritized and functionally tested. Silencing of para did not influence P. falciparum infection, while silencing of a predicted immune gene, serine protease ClipC9, allowed development of significantly increased parasite numbers. CONCLUSIONS: Genetic variation influencing Plasmodium infection in wild Anopheles is linked to a natural ~3 megabase haplotype on chromosome 2L that carries the kdr allele of the para gene. Evidence suggests that para gene function does not directly influence parasite susceptibility, and the association of kdr with infection may be due to tight linkage of kdr with other gene(s) on the haplotype. Further work will be required to determine if ClipC9 influences the outcome of P. falciparum infection in nature, as well as to confirm the absence of a direct influence by para.

Mutations in PYCR2, Encoding Pyrroline-5-Carboxylate Reductase 2, Cause Microcephaly and Hypomyelination.

Despite recent advances in understanding the genetic bases of microcephaly, a large number of cases of microcephaly remain unexplained, suggesting that many microcephaly syndromes and associated genes have yet to be identified. Here, we report mutations in PYCR2, which encodes an enzyme in the proline biosynthesis pathway, as the cause of a unique syndrome characterized by postnatal microcephaly, hypomyelination, and reduced cerebral white-matter volume. Linkage mapping and whole-exome sequencing identified homozygous mutations (c.355C>T [p.Arg119Cys] and c.751C>T [p.Arg251Cys]) in PYCR2 in the affected individuals of two consanguineous families. A lymphoblastoid cell line from one affected individual showed a strong reduction in the amount of PYCR2. When mutant cDNAs were transfected into HEK293FT cells, both variant proteins retained normal mitochondrial localization but had lower amounts than the wild-type protein, suggesting that the variant proteins were less stable. A PYCR2-deficient HEK293FT cell line generated by genome editing with the clustered regularly interspaced short palindromic repeat (CRISPR)-Cas9 system showed that PYCR2 loss of function led to decreased mitochondrial membrane potential and increased susceptibility to apoptosis under oxidative stress. Morpholino-based knockdown of a zebrafish PYCR2 ortholog, pycr1b, recapitulated the human microcephaly phenotype, which was rescued by wild-type human PYCR2 mRNA, but not by mutant mRNAs, further supporting the pathogenicity of the identified variants. Hypomyelination and the absence of lax, wrinkly skin distinguishes this condition from that caused by previously reported mutations in the gene encoding PYCR2's isozyme, PYCR1, suggesting a unique and indispensable role for PYCR2 in the human CNS during development.