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.

Abstract PR10: FBXO11 is recurrently mutated in Burkitt lymphoma and its inactivation accelerates lymphomagenesis in Eμ-myc mice

Introduction: BCL6 is a key oncogene in lymphoma pathogenesis. Malignant lymphoid cells exploit several mechanisms to deregulate the expression of BCL6, including chromosomal translocations, somatic mutations in the promoter regulatory regions, or inactivation of the pathway that controls its degradation. FBXO11 was recently identified as a major ubiquitin ligase involved in the degradation of BCL6 and was found to be frequently inactivated by mutations or deletion in diffuse large B cell lymphoma (DLBCL). Thus, FBXO11 acts as an oncosuppressor in DLBCL by promoting the accumulation of BCL6. Given the prominent role of FBXO11 in regulating BCL6 stability, we searched for FBXO11 mutations in BCL6-positive lymphomas, other than DLBCL, and we investigated its role in lymphoma development in vivo.

Methods: We sequenced the entire FBXO11 coding sequence by classical Sanger sequencing in 100 cases of follicular lymphoma (FL), 36 cases of Burkitt lymphoma (BL), 8 BL cell lines and 8 anaplastic large cell lymphoma cell lines, all BCL6-positive lymphomas. Moreover, we sequenced 50 cases of marginal zone B cell lymphoma (MZL), which show variable expression of BCL6.

We functionally validated the FBXO11 mutations by developing mutant constructs and testing their ability to induce BCL6 and SNAIL degradation. We then applied the CRISPR/Cas9 system to disrupt the endogenous FBXO11 gene in BL cells and evaluated its effect on BCL6 stability.

To dissect the in vivo role of FBXO11 in lymphomagenesis we generated conditional FBXO11 knock-out (KO) mice (FBXO11fl/fl) to delete protein expression in CD19-positive B cells. To investigate whether FBXO11 inactivation cooperates with c-myc in lymphomagenesis we crossed CD19/Cre-FBXO11fl/fl mice with Eμ-myc transgenic mice.

Results: We identified FBXO11 mutations in BL cases and cell lines (10/44, 22.7%), one case of FL (1/100) and one case of MZL (1/50). In FL and MZL, the mutational analysis of tissue collected by microdissection showed that the mutation was specifically acquired by the high-grade component. The frequency of FBXO11 mutation in BL was further validated on an independent cohort of 51 BL cases (6/51, 11.7%), obtained from published data of a previous whole-exome sequencing study (Love et al. Nat Genet 2012).

BL mutations were mostly missense mutations located in the functional CASH domains and also splice-site mutations that were never described before and that induced alternative splicing, as confirmed on the mRNA extracted from the tumor samples. All mutations produced a mutant FBXO11 with impaired ability to induce BCL6 and SNAIL degradation. CRISPR/Cas9 mediated KO of FBXO11 in BL cells resulted in an almost complete stabilization of BCL6, thus suggesting that FBXO11 is the main, if not unique, ubiquitin ligase that controls BCL6 stability in BL.

Finally, we observed an acceleration of lymphoma development in the CD19/Cre-FBXO11fl/fl mice crossed with Eμ-myc transgenic mice. The lymphomas showed histologic features of high-grade disease with a more mature B-cell phenotype than the Eμ-myc tumors alone.

Conclusions: Overall our results demonstrate that FBXO11 is frequently mutated in BL. All mutants identified impair the ability of FBXO11 to regulate the degradation of BCL6. Together with our previous findings (Duan et al. Nature 2012), this study shows that FBXO11 is mostly mutated in aggressive lymphomas such as DLBCL and BL, and suggests that FBXO11 mutations could contribute to their aggressiveness. In fact, we show that FBXO11 inactivation cooperates with c-myc in accelerating lymphomagenesis. We have established a novel murine lymphoma model that resembles more closely the human BL, providing a novel promising tool for the study of lymphomagenesis and, potentially, preclinical testing of therapeutic approaches with BCL6 and/or c-myc inhibitors.

This abstract is also being presented as Poster 32.

Citation Format: Chiara Pighi, Mara Compagno, Taek-Chin Cheong, Teresa Poggio, Qi Wang, Fernanda Langellotto, Anoop Sendamarai, Kyriacos Markianos, Paola Francia di Celle, Alberto Zamò, Roberto Chiarle. FBXO11 is recurrently mutated in Burkitt lymphoma and its inactivation accelerates lymphomagenesis in Eμ-myc mice [abstract]. In: Proceedings of the Second AACR Conference on Hematologic Malignancies: Translating Discoveries to Novel Therapies; May 6-9, 2017; Boston, MA. Philadelphia (PA): AACR; Clin Cancer Res 2017;23(24_Suppl):Abstract nr PR10.

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.

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.

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.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2009-z) contains supplementary material, which is available to authorized users.

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

Background

Members of the Anophelesgambiae species complex are primary vectors of human malaria in Africa. It is known that a large haplotype shared between An. gambiae and Anophelescoluzzii 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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0924-8) contains supplementary material, which is available to authorized users.

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.

Loss of PCLO function underlies pontocerebellar hypoplasia type III

OBJECTIVE

To identify the genetic cause of pontocerebellar hypoplasia type III (PCH3).

METHODS

We studied the original reported pedigree of PCH3 and performed genetic analysis including genome-wide single nucleotide polymorphism genotyping, linkage analysis, whole-exome sequencing, and Sanger sequencing. Human fetal brain RNA sequencing data were then analyzed for the identified candidate gene.

RESULTS

The affected individuals presented with severe global developmental delay and seizures starting in the first year of life. Brain MRI of an affected individual showed diffuse atrophy of the cerebrum, cerebellum, and brainstem. Genome-wide single nucleotide polymorphism analysis confirmed the linkage to chromosome 7q we previously reported, and showed no other genomic areas of linkage. Whole-exome sequencing of 2 affected individuals identified a shared homozygous, nonsense variant in the PCLO (piccolo) gene. This variant segregated with the disease phenotype in the pedigree was rare in the population and was predicted to eliminate the PDZ and C2 domains in the C-terminus of the protein. RNA sequencing data of human fetal brain showed that PCLO was moderately expressed in the developing cerebral cortex.

CONCLUSIONS

Here, we show that a homozygous, nonsense PCLO mutation underlies the autosomal recessive neurodegenerative disorder, PCH3. PCLO is a component of the presynaptic cytoskeletal matrix, and is thought to be involved in regulation of presynaptic proteins and synaptic vesicles. Our findings suggest that PCLO is crucial for the development and survival of a wide range of neuronal types in the human brain.

Expanding the phenotype associated with the NEFL mutation: neuromuscular disease in a family with overlapping myopathic and neurogenic findings

IMPORTANCE

Newer sequencing technologies in combination with traditional gene mapping techniques, such as linkage analysis, can help identify the genetic basis of disease for patients with rare disorders of uncertain etiology. This approach may expand the phenotypic spectrum of disease associated with those genetic mutations.

OBJECTIVE

To elucidate the molecular cause of a neuromuscular disease among a family in which 4 members, a mother and her 3 sons, were affected.

DESIGN, SETTING, AND PARTICIPANTS

Two of 4 affected members manifested nemaline myopathy, a common subtype of congenital myopathy, while the other 2 had a nonspecific myopathy. Single-nucleotide polymorphism-based linkage analysis was performed on DNA samples from the 4 affected family members, and whole-genome sequencing was performed in the proband. Real-time quantitative reverse transcription-polymerase chain reaction, immunofluorescence, and Western blot analysis were performed on muscle biopsy specimens.

MAIN OUTCOMES AND MEASURES

Whole-genome sequencing and linkage analysis identified a variant in a gene that explains the phenotype.

RESULTS

We identified a novel neurofilament light polypeptide (NEFL) nonsense mutation in all affected members. NEFL mutations have been previously linked to Charcot-Marie-Tooth disease in humans. This led us to reevaluate the diagnosis, and we recognized that several of the findings, especially those related to the muscle biopsy specimens and electromyography, were consistent with a neurogenic disease.

CONCLUSIONS AND RELEVANCE

NEFL mutations are known to cause Charcot-Marie-Tooth disease in humans and motor neuron disease in mice. We report the identification of an NEFL mutation in a family clinically manifesting congenital myopathy. We also describe potential overlap between myopathic and neurogenic findings in this family. These findings expand the phenotypic spectrum of diseases associated with NEFL mutations. This study is an example of the power of genomic approaches to identify potentially pathogenic mutations in unsuspected genes responsible for heterogeneous neuromuscular diseases.

X-linked sideroblastic anemia due to ALAS2 intron 1 enhancer element GATA-binding site mutations

X-linked sideroblastic anemia (XLSA) is the most common form of congenital sideroblastic anemia. In affected males, it is uniformly associated with partial loss-of-function missense mutations in the erythroid-specific heme biosynthesis protein 5-aminolevulinate synthase 2 (ALAS2). Here, we report five families with XLSA owing to mutations in a GATA transcription factor binding site located in a transcriptional enhancer element in intron 1 of the ALAS2 gene. As such, this study defines a new class of mutations that should be evaluated in patients undergoing genetic testing for a suspected diagnosis of XLSA.

Comparative RNA editing in autistic and neurotypical cerebella

Adenosine-to-inosine (A-to-I) RNA editing is a neurodevelopmentally regulated epigenetic modification shown to modulate complex behavior in animals. Little is known about human A-to-I editing, but it is thought to constitute one of many molecular mechanisms connecting environmental stimuli and behavioral outputs. Thus, comprehensive exploration of A-to-I RNA editing in human brains may shed light on gene-environment interactions underlying complex behavior in health and disease. Synaptic function is a main target of A-to-I editing, which can selectively recode key amino acids in synaptic genes, directly altering synaptic strength and duration in response to environmental signals. Here, we performed a high-resolution survey of synaptic A-to-I RNA editing in a human population, and examined how it varies in autism, a neurodevelopmental disorder in which synaptic abnormalities are a common finding. Using ultra-deep (>1000 × ) sequencing, we quantified the levels of A-to-I editing of 10 synaptic genes in postmortem cerebella from 14 neurotypical and 11 autistic individuals. A high dynamic range of editing levels was detected across individuals and editing sites, from 99.6% to below detection limits. In most sites, the extreme ends of the population editing distributions were individuals with autism. Editing was correlated with isoform usage, clusters of correlated sites were identified, and differential editing patterns examined. Finally, a dysfunctional form of the editing enzyme adenosine deaminase acting on RNA B1 was found more commonly in postmortem cerebella from individuals with autism. These results provide a population-level, high-resolution view of A-to-I RNA editing in human cerebella and suggest that A-to-I editing of synaptic genes may be informative for assessing the epigenetic risk for autism.

Recessive truncating titin gene, TTN, mutations presenting as centronuclear myopathy

OBJECTIVE

To identify causative genes for centronuclear myopathies (CNM), a heterogeneous group of rare inherited muscle disorders that often present in infancy or early life with weakness and hypotonia, using next-generation sequencing of whole exomes and genomes.

METHODS

Whole-exome or -genome sequencing was performed in a cohort of 29 unrelated patients with clinicopathologic diagnoses of CNM or related myopathy depleted for cases with mutations of MTM1, DNM2, and BIN1. Immunofluorescence analyses on muscle biopsies, splicing assays, and gel electrophoresis of patient muscle proteins were performed to determine the molecular consequences of mutations of interest.

RESULTS

Autosomal recessive compound heterozygous truncating mutations of the titin gene, TTN, were identified in 5 individuals. Biochemical analyses demonstrated increased titin degradation and truncated titin proteins in patient muscles, establishing the impact of the mutations.

CONCLUSIONS

Our study identifies truncating TTN mutations as a cause of congenital myopathy that is reported as CNM. Unlike the classic CNM genes that are all involved in excitation-contraction coupling at the triad, TTN encodes the giant sarcomeric protein titin, which forms a myofibrillar backbone for the components of the contractile machinery. This study expands the phenotypic spectrum associated with TTN mutations and indicates that TTN mutation analysis should be considered in cases of possible CNM without mutations in the classic CNM genes.

A novel rat model of hereditary hemochromatosis due to a mutation in transferrin receptor 2

Sporadic iron overload in rats has been reported, but whether it is due to genetic or environmental causes is unknown. In the current study, phenotypic analysis of Hsd:HHCL Wistar rats revealed a low incidence of histologically detected liver iron overload. Here we characterized the pathophysiology of the iron overload and showed that the phenotype is heritable and due to a mutation in a single gene. We identified a single male rat among the 132 screened animals that exhibited predominantly periportal, hepatocellular iron accumulation. This rat expressed low RNA levels of the iron regulatory hormone hepcidin and low protein levels of transferrin receptor 2 (Tfr2), a membrane protein essential for hepcidin expression in humans and mice and mutated in forms of hereditary hemochromatosis. Sequencing of Tfr2 in the iron-overloaded rat revealed a novel Ala679Gly polymorphism in a highly conserved residue. Quantitative trait locus mapping indicated that this polymorphism correlated strongly with serum iron and transferrin saturations in male rats. Expression of the Gly679 variant in tissue culture cell lines revealed decreased steady-state levels of Tfr2. Characterization of iron metabolism in the progeny of polymorphic rats suggested that homozygosity for the Ala679Gly allele leads to a hemochromatosis phenotype. However, we currently cannot exclude the possibility that a polymorphism or mutation in the noncoding region of Tfr2 contributes to the iron-overload phenotype. Hsd:HHCL rats are the first genetic rat model of hereditary hemochromatosis and may prove useful for understanding the molecular mechanisms underlying the regulation of iron metabolism.

Using whole-exome sequencing to identify inherited causes of autism

Despite significant heritability of autism spectrum disorders (ASDs), their extreme genetic heterogeneity has proven challenging for gene discovery. Studies of primarily simplex families have implicated de novo copy number changes and point mutations, but are not optimally designed to identify inherited risk alleles. We apply whole-exome sequencing (WES) to ASD families enriched for inherited causes due to consanguinity and find familial ASD associated with biallelic mutations in disease genes (AMT, PEX7, SYNE1, VPS13B, PAH, and POMGNT1). At least some of these genes show biallelic mutations in nonconsanguineous families as well. These mutations are often only partially disabling or present atypically, with patients lacking diagnostic features of the Mendelian disorders with which these genes are classically associated. Our study shows the utility of WES for identifying specific genetic conditions not clinically suspected and the importance of partial loss of gene function in ASDs.

A cryptic subgroup of Anopheles gambiae is highly susceptible to human malaria parasites

Population subgroups of the African malaria vector Anopheles gambiae have not been comprehensively characterized owing to the lack of unbiased sampling methods. In the arid savanna zone of West Africa, where potential oviposition sites are scarce, widespread collection from larval pools in the peridomestic human habitat yielded a comprehensive genetic survey of local A. gambiae population subgroups, independent of adult resting behavior and ecological preference. A previously unknown subgroup of exophilic A. gambiae is sympatric with the known endophilic A. gambiae in this region. The exophilic subgroup is abundant, lacks differentiation into M and S molecular forms, and is highly susceptible to infection with wild Plasmodium falciparum. These findings might have implications for the epidemiology of malaria transmission and control.

Lack of association of the serotonin transporter polymorphism with the sudden infant death syndrome in the San Diego Dataset

Dysfunction of medullary serotonin (5-HT)-mediated respiratory and autonomic function is postulated to underlie the pathogenesis of the majority of sudden infant death syndrome (SIDS) cases. Several studies have reported an increased frequency of the LL genotype and L allele of the 5-HT transporter (5-HTT) gene promoter polymorphism (5-HTTLPR), which is associated with increased transcriptional activity and 5-HT transport in vitro, in SIDS cases compared with controls. These findings raise the possibility that this polymorphism contributes to or exacerbates existing medullary 5-HT dysfunction in SIDS. In this study, we tested the hypothesis that the frequency of LL genotype and L allele are higher in 179 SIDS cases compared with 139 controls of multiple ethnicities in the San Diego SIDS Dataset. We observed no significant association of genotype or allele with SIDS cases either in the total cohort or on stratification for ethnicity. These observations do not support previous findings that the L allele and/or LL genotype of the 5-HTTLPR are associated with SIDS.

Serotonin-related FEV gene variant in the sudden infant death syndrome is a common polymorphism in the African-American population

An important subset of the sudden infant death syndrome (SIDS) is associated with multiple serotonergic (5-HT) abnormalities in regions of the medulla oblongata. The mouse ortholog of the fifth Ewing variant gene (FEV) is critical for 5-HT neuronal development. A putatively rare intronic variant [IVS2-191_190insA, here referred to as c.128-(191_192)dupA] has been reported as a SIDS-associated mutation in an African-American population. We tested this association in an independent dataset: 137 autopsied cases (78 SIDS, 59 controls) and an additional 296 control DNA samples from Coriell Cell Repositories. In addition to the c.128-(191_192)dupA variant, we observed an associated single-base deletion [c.128-(301-306)delG] in a subset of the samples. Neither of the two FEV variants showed significant association with SIDS in either the African-American subgroup or the overall cohort. Although we found a significant association of c.128-(191_192)dupA with SIDS when San Diego Hispanic SIDS cases were compared with San Diego Hispanic controls plus Mexican controls (p = 0.04), this became nonsignificant after multiple testing correction. Among Coriell controls, 33 of 99 (33%) African-American and 0 of 197 (0%) of the remaining controls carry the polymorphism (c.128-(191_192)dupA). The polymorphism seems to be a common, likely nonpathogenic, variant in the African-American population.

Identifying autism loci and genes by tracing recent shared ancestry

To find inherited causes of autism-spectrum disorders, we studied families in which parents share ancestors, enhancing the role of inherited factors. We mapped several loci, some containing large, inherited, homozygous deletions that are likely mutations. The largest deletions implicated genes, including PCDH10 (protocadherin 10) and DIA1 (deleted in autism1, or c3orf58), whose level of expression changes in response to neuronal activity, a marker of genes involved in synaptic changes that underlie learning. A subset of genes, including NHE9 (Na+/H+ exchanger 9), showed additional potential mutations in patients with unrelated parents. Our findings highlight the utility of "homozygosity mapping" in heterogeneous disorders like autism but also suggest that defective regulation of gene expression after neural activity may be a mechanism common to seemingly diverse autism mutations.

Mutations in TMPRSS6 cause iron-refractory iron deficiency anemia (IRIDA)

Iron deficiency is usually attributed to chronic blood loss or inadequate dietary intake. Here, we show that iron deficiency anemia refractory to oral iron therapy can be caused by germline mutations in TMPRSS6, which encodes a type II transmembrane serine protease produced by the liver that regulates the expression of the systemic iron regulatory hormone hepcidin. These findings demonstrate that TMPRSS6 is essential for normal systemic iron homeostasis in humans.

A major genetic locus controlling natural Plasmodium falciparum infection is shared by East and West African Anopheles gambiae

Background

Genetic linkage mapping identified a region of chromosome 2L in the Anopheles gambiae genome that exerts major control over natural infection by Plasmodium falciparum. This 2L Plasmodium-resistance interval was mapped in mosquitoes from a natural population in Mali, West Africa, and controls the numbers of P. falciparum oocysts that develop on the vector midgut. An important question is whether genetic variation with respect to Plasmodium-resistance exists across Africa, and if so whether the same or multiple geographically distinct resistance mechanisms are responsible for the trait.

Methods

To identify P falciparum resistance loci in pedigrees generated and infected in Kenya, East Africa, 28 microsatellite loci were typed across the mosquito genome. Genetic linkage mapping was used to detect significant linkage between genotype and numbers of midgut oocysts surviving to 7–8 days post-infection.

Results

A major malaria-control locus was identified on chromosome 2L in East African mosquitoes, in the same apparent position originally identified from the West African population. Presence of this resistance locus explains 75% of parasite free mosquitoes. The Kenyan resistance locus is named EA_Pfin1 (East Africa_ Plasmodium falciparum Infection Intensity).

Conclusion

Detection of a malaria-control locus at the same chromosomal location in both East and West African mosquitoes indicates that, to the level of genetic resolution of the analysis, the same mechanism of Plasmodium-resistance, or a mechanism controlled by the same genomic region, is found across Africa, and thus probably operates in A. gambiae throughout its entire range.

Natural Malaria Infection in Anopheles gambiae Is Regulated by a Single Genomic Control Region

We surveyed an Anopheles gambiae population in a West African malaria transmission zone for naturally occurring genetic loci that control mosquito infection with the human malaria parasite, Plasmodium falciparum. The strongest Plasmodium resistance loci cluster in a small region of chromosome 2L and each locus explains at least 89% of parasite-free mosquitoes in independent pedigrees. Together, the clustered loci form a genomic Plasmodium-resistance island that explains most of the genetic variation for malaria parasite infection of mosquitoes in nature. Among the candidate genes in this chromosome region, RNA interference knockdown assays confirm a role in Plasmodium resistance for Anopheles Plasmodium-responsive leucine-rich repeat 1 (APL1), encoding a leucine-rich repeat protein that is similar to molecules involved in natural pathogen resistance mechanisms in plants and mammals.

Finding prostate cancer susceptibility genes

Prostate cancer is a heterogeneous disease with multiple loci contributing to susceptibility. Traditionally, genome-wide scans using high-risk families have utilized stratification by number of affected individuals, family history of other cancers, or family age at diagnosis to improve genetic homogeneity. In addition to locus heterogeneity, for later onset diseases such as prostate cancer, a major limitation to mapping efforts is that key parental DNA samples are rarely available. The lack of available samples from upper generations reduces inheritance information, and as a result, the standard 10-cM genome scan does not provide full power to detect linkage. To increase the ability to find disease-associated loci, much denser genome-wide scans must be undertaken in multiple ethnic groups. In addition, new ways of defining homogenous subsets of families need to be developed.

Sequence-based linkage analysis

The rapid decrease in the cost of DNA sequencing will enable its use for novel applications. Here, we investigate the use of DNA sequencing for simultaneous discovery and genotyping of polymorphisms in family linkage studies. In the proposed approach, short contiguous segments of genomic DNA, regularly spaced across the genome, are resequenced in each pedigree member, and all sequence polymorphisms discovered within a pedigree are used as genetic markers. We use computer simulations consistent with observed human sequence diversity to show that segments of 500-1,000 base pairs, spaced at intervals of 1-2 Mb across the genome, provide linkage information that equals or exceeds that of traditional marker-based approaches. We validate these results experimentally by implementing the sequence-based linkage approach for chromosome 19 in CEPH pedigrees.

A 3.9-centimorgan-resolution human single-nucleotide polymorphism linkage map and screening set

Recent advances in technologies for high-throughout single-nucleotide polymorphism (SNP)-based genotyping have improved efficiency and cost so that it is now becoming reasonable to consider the use of SNPs for genomewide linkage analysis. However, a suitable screening set of SNPs and a corresponding linkage map have yet to be described. The SNP maps described here fill this void and provide a resource for fast genome scanning for disease genes. We have evaluated 6,297 SNPs in a diversity panel composed of European Americans, African Americans, and Asians. The markers were assessed for assay robustness, suitable allele frequencies, and informativeness of multi-SNP clusters. Individuals from 56 Centre d'Etude du Polymorphisme Humain pedigrees, with >770 potentially informative meioses altogether, were genotyped with a subset of 2,988 SNPs, for map construction. Extensive genotyping-error analysis was performed, and the resulting SNP linkage map has an average map resolution of 3.9 cM, with map positions containing either a single SNP or several tightly linked SNPs. The order of markers on this map compares favorably with several other linkage and physical maps. We compared map distances between the SNP linkage map and the interpolated SNP linkage map constructed by the deCode Genetics group. We also evaluated cM/Mb distance ratios in females and males, along each chromosome, showing broadly defined regions of increased and decreased rates of recombination. Evaluations indicate that this SNP screening set is more informative than the Marshfield Clinic's commonly used microsatellite-based screening set.

Genetic loci affecting resistance to human malaria parasites in a West African mosquito vector population

Successful propagation of the malaria parasite Plasmodium falciparum within a susceptible mosquito vector is a prerequisite for the transmission of malaria. A field-based genetic analysis of the major human malaria vector, Anopheles gambiae, has revealed natural factors that reduce the transmission of P. falciparum. Differences in P. falciparum oocyst numbers between mosquito isofemale families fed on the same infected blood indicated a large genetic component affecting resistance to the parasite, and genome-wide scanning in pedigrees of wild mosquitoes detected segregating resistance alleles. The apparently high natural frequency of resistance alleles suggests that malaria parasites (or a similar pathogen) exert a significant selective pressure on vector populations.

A joint analysis of asthma affection status and IgE levels in multiple data sets collected for asthma

We present a joint linkage analysis of eight data sets collected for asthma. Three of the data sets are full genome scans, while the remaining five concentrate on a 40-cM region on chromosome 5. We perform the analysis using one qualitative and one quantitative phenotype: asthma status and IgE level. Considering all data sets simultaneously, we do not find evidence for linkage to asthma affection status beyond the level expected to occur by chance twice per genome scan. In contrast, we observe significant linkage to IgE level on chromosome 6.

Lambda(0) polarization in 800-GeV/cpp --> p(f)(Lambda0K+)

We report results from a study of Lambda(0) polarization in the exclusive reaction pp --> p(f)(Lambda0K+) at 800-GeV/c. We observe a dependence of the polarization on the Lambda0K+ invariant mass with large (+71%) positive polarization at small mass (1.63-GeV/c(2)) and large (-43%) negative polarization at large mass (2.75-GeV/c(2)). This observation confirms the result of the CERN ISR R608 experiment and extends the range over which the effect is observed. The strong dependence of the polarization on the Lambda0K+ invariant mass suggests that the origin of the polarization is closely related to the production dynamics of the diffractively produced Lambda0K+ system.

Diffractively produced charm final states in 800-GeV/c pp collisions

We report the first observation of diffractively produced open charm in 800-GeV/c pp collisions of the type pp-->pD*X. We measure cross sections of sigma(diff)(D*+) = (0.185+/-0.044+/-0.054) (mu)b and sigma(diff)(D(*-)) = (0.174+/-0.034+/-0.029) (mu)b. Our measurements are based on 4.3x10(9) events recorded by FNAL E690 in the fixed-target run of 1991. We compare our results with previous fixed-target charm experiments.

Efficient multipoint linkage analysis through reduction of inheritance space

Computational constraints currently limit exact multipoint linkage analysis to pedigrees of moderate size. We introduce new algorithms that allow analysis of larger pedigrees by reducing the time and memory requirements of the computation. We use the observed pedigree genotypes to reduce the number of inheritance patterns that need to be considered. The algorithms are implemented in a new version (version 2.1) of the software package GENEHUNTER. Performance gains depend on marker heterozygosity and on the number of pedigree members available for genotyping, but typically are 10-1,000-fold, compared with the performance of the previous release (version 2.0). As a result, families with up to 30 bits of inheritance information have been analyzed, and further increases in family size are feasible. In addition to computation of linkage statistics and haplotype determination, GENEHUNTER can also perform single-locus and multilocus transmission/disequilibrium tests. We describe and implement a set of permutation tests that allow determination of empirical significance levels in the presence of linkage disequilibrium among marker loci.