A concerted neuron-astrocyte program declines in ageing and schizophrenia.

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA sequencing to analyse the prefrontal cortex of 191 human donors aged 22-97 years, including healthy individuals and people with schizophrenia. Latent-factor analysis of these data revealed that, in people whose cortical neurons more strongly expressed genes encoding synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the synaptic neuron and astrocyte program (SNAP). In schizophrenia and ageing-two conditions that involve declines in cognitive flexibility and plasticity1,2-cells divested from SNAP: astrocytes, glutamatergic (excitatory) neurons and GABAergic (inhibitory) neurons all showed reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy people of similar age, may underlie many aspects of normal human interindividual differences and may be an important point of convergence for multiple kinds of pathophysiology.

Genetic drivers and cellular selection of female mosaic X chromosome loss.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Rare coding variants in ten genes confer substantial risk for schizophrenia.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Monogenic and polygenic inheritance become instruments for clonal selection.

Clonally expanded blood cells that contain somatic mutations (clonal haematopoiesis) are commonly acquired with age and increase the risk of blood cancer1-9. The blood clones identified so far contain diverse large-scale mosaic chromosomal alterations (deletions, duplications and copy-neutral loss of heterozygosity (CN-LOH)) on all chromosomes1,2,5,6,9, but the sources of selective advantage that drive the expansion of most clones remain unknown. Here, to identify genes, mutations and biological processes that give selective advantage to mutant clones, we analysed genotyping data from the blood-derived DNA of 482,789 participants from the UK Biobank10. We identified 19,632 autosomal mosaic chromosomal alterations and analysed these for relationships to inherited genetic variation. We found 52 inherited, rare, large-effect coding or splice variants in 7 genes that were associated with greatly increased vulnerability to clonal haematopoiesis with specific acquired CN-LOH mutations. Acquired mutations systematically replaced the inherited risk alleles (at MPL) or duplicated them to the homologous chromosome (at FH, NBN, MRE11, ATM, SH2B3 and TM2D3). Three of the genes (MRE11, NBN and ATM) encode components of the MRN-ATM pathway, which limits cell division after DNA damage and telomere attrition11-13; another two (MPL and SH2B3) encode proteins that regulate the self-renewal of stem cells14-16. In addition, we found that CN-LOH mutations across the genome tended to cause chromosomal segments with alleles that promote the expansion of haematopoietic cells to replace their homologous (allelic) counterparts, increasing polygenic drive for blood-cell proliferation traits. Readily acquired mutations that replace chromosomal segments with their homologous counterparts seem to interact with pervasive inherited variation to create a challenge for lifelong cytopoiesis.

Serum biomarkers are altered in UK Biobank participants with mosaic chromosomal alterations.

Age-related clonal expansion of cells harbouring mosaic chromosomal alterations (mCAs) is one manifestation of clonal haematopoiesis. Identifying factors that influence the generation and promotion of clonal expansion of mCAs are key to investigate the role of mCAs in health and disease. Herein, we report on widely measured serum biomarkers and their possible association with mCAs, which could provide new insights into molecular alterations that promote acquisition and clonal expansion. We performed a cross-sectional investigation of the association of 32 widely measured serum biomarkers with autosomal mCAs, mosaic loss of the Y chromosome, and mosaic loss of the X chromosome in 436 784 cancer-free participants from the UK Biobank. mCAs were associated with a range of commonly measured serum biomarkers such as lipid levels, circulating sex hormones, blood sugar homeostasis, inflammation and immune function, vitamins and minerals, kidney function, and liver function. Biomarker levels in participants with mCAs were estimated to differ by up to 5% relative to mCA-free participants, and individuals with higher cell fraction mCAs had greater deviation in mean biomarker values. Polygenic scores associated with sex hormone binding globulin, vitamin D, and total cholesterol were also associated with mCAs. Overall, we observed commonly used clinical serum biomarkers related to disease risk are associated with mCAs, suggesting mechanisms involved in these diseases could be related to mCA proliferation and clonal expansion.

BCFtools/liftover: an accurate and comprehensive tool to convert genetic variants across genome assemblies.

MOTIVATION: Many genetics studies report results tied to genomic coordinates of a legacy genome assembly. However, as assemblies are updated and improved, researchers are faced with either realigning raw sequence data using the updated coordinate system or converting legacy datasets to the updated coordinate system to be able to combine results with newer datasets. Currently available tools to perform the conversion of genetic variants have numerous shortcomings, including poor support for indels and multi-allelic variants, that lead to a higher rate of variants being dropped or incorrectly converted. As a result, many researchers continue to work with and publish using legacy genomic coordinates. RESULTS: Here we present BCFtools/liftover, a tool to convert genomic coordinates across genome assemblies for variants encoded in the variant call format with improved support for indels represented by different reference alleles across genome assemblies and full support for multi-allelic variants. It further supports variant annotation fields updates whenever the reference allele changes across genome assemblies. The tool has the lowest rate of variants being dropped with an order of magnitude less indels dropped or incorrectly converted and is an order of magnitude faster than other tools typically used for the same task. It is particularly suited for converting variant callsets from large cohorts to novel telomere-to-telomere assemblies as well as summary statistics from genome-wide association studies tied to legacy genome assemblies. AVAILABILITY AND IMPLEMENTATION: The tool is written in C and freely available under the MIT open source license as a BCFtools plugin available at http://github.com/freeseek/score.

Genetic predisposition to mosaic Y chromosome loss in blood.

Mosaic loss of chromosome Y (LOY) in circulating white blood cells is the most common form of clonal mosaicism1-5, yet our knowledge of the causes and consequences of this is limited. Here, using a computational approach, we estimate that 20% of the male population represented in the UK Biobank study (n = 205,011) has detectable LOY. We identify 156 autosomal genetic determinants of LOY, which we replicate in 757,114 men of European and Japanese ancestry. These loci highlight genes that are involved in cell-cycle regulation and cancer susceptibility, as well as somatic drivers of tumour growth and targets of cancer therapy. We demonstrate that genetic susceptibility to LOY is associated with non-haematological effects on health in both men and women, which supports the hypothesis that clonal haematopoiesis is a biomarker of genomic instability in other tissues. Single-cell RNA sequencing identifies dysregulated expression of autosomal genes in leukocytes with LOY and provides insights into why clonal expansion of these cells may occur. Collectively, these data highlight the value of studying clonal mosaicism to uncover fundamental mechanisms that underlie cancer and other ageing-related diseases.

A Phenome-Wide Association Study of genes associated with COVID-19 severity reveals shared genetics with complex diseases in the Million Veteran Program.

The study aims to determine the shared genetic architecture between COVID-19 severity with existing medical conditions using electronic health record (EHR) data. We conducted a Phenome-Wide Association Study (PheWAS) of genetic variants associated with critical illness (n = 35) or hospitalization (n = 42) due to severe COVID-19 using genome-wide association summary data from the Host Genetics Initiative. PheWAS analysis was performed using genotype-phenotype data from the Veterans Affairs Million Veteran Program (MVP). Phenotypes were defined by International Classification of Diseases (ICD) codes mapped to clinically relevant groups using published PheWAS methods. Among 658,582 Veterans, variants associated with severe COVID-19 were tested for association across 1,559 phenotypes. Variants at the ABO locus (rs495828, rs505922) associated with the largest number of phenotypes (nrs495828 = 53 and nrs505922 = 59); strongest association with venous embolism, odds ratio (ORrs495828 1.33 (p = 1.32 x 10-199), and thrombosis ORrs505922 1.33, p = 2.2 x10-265. Among 67 respiratory conditions tested, 11 had significant associations including MUC5B locus (rs35705950) with increased risk of idiopathic fibrosing alveolitis OR 2.83, p = 4.12 × 10-191; CRHR1 (rs61667602) associated with reduced risk of pulmonary fibrosis, OR 0.84, p = 2.26× 10-12. The TYK2 locus (rs11085727) associated with reduced risk for autoimmune conditions, e.g., psoriasis OR 0.88, p = 6.48 x10-23, lupus OR 0.84, p = 3.97 x 10-06. PheWAS stratified by ancestry demonstrated differences in genotype-phenotype associations. LMNA (rs581342) associated with neutropenia OR 1.29 p = 4.1 x 10-13 among Veterans of African and Hispanic ancestry but not European. Overall, we observed a shared genetic architecture between COVID-19 severity and conditions related to underlying risk factors for severe and poor COVID-19 outcomes. Differing associations between genotype-phenotype across ancestries may inform heterogenous outcomes observed with COVID-19. Divergent associations between risk for severe COVID-19 with autoimmune inflammatory conditions both respiratory and non-respiratory highlights the shared pathways and fine balance of immune host response and autoimmunity and caution required when considering treatment targets.

An extra X chromosome among adult women in the Million Veteran Program: A more benign perspective of trisomy X.

Despite affecting in 1 in every 1000 females, remarkably little is known about trisomy X syndrome (47,XXX), especially among older adults who are undiagnosed. In this study, we aimed to determine the prevalence of 47,XXX among females enrolled in the Million Veterans Program (MVP; mean age 50.2 ± 13.6 years), and compare broad health outcomes between females with 47,XXX and 46,XX matched controls. We identified 61 females with an additional X chromosome, corresponding to a prevalence of 103 per 100,000 females; 27.9% had been clinically diagnosed. Females with 47,XXX had taller stature (+6.1 cm, p < 0.001), greater rate of outpatient encounters (p = 0.026), higher odds of kidney disease (odds ratio [OR] = 12.3; 95% confidence interval [CI] 2.9-51.8), glaucoma (OR = 5.1; 95% CI 1.5-13.9), and congestive heart failure (OR = 5.6; 95% CI 1.4-24.2), and were more likely to be unemployed (p = 0.008) with lower annual income (p = 0.021) when compared with 46,XX controls of the same age and genetic ancestry. However, there were no differences in the rates of other encounter types, Charlson Comorbidity Index, all other medical and psychological diagnoses, military service history or quality of life metrics. In conclusion, in this aging and predominately undiagnosed sample, 47,XXX conferred few differences when compared with matched controls, offering a more reassuring perspective to the trisomy X literature.

A Neanderthal haplotype introgressed into the human genome confers protection against membranous nephropathy.

Class 2 HLA and PLA2R1 alleles are exceptionally strong genetic risk factors for membranous nephropathy (MN), leading, through an unknown mechanism, to a targeted autoimmune response. Introgressed archaic haplotypes (introduced from an archaic human genome into the modern human genome) might influence phenotypes through gene dysregulation. Here, we investigated the genomic region surrounding the PLA2R1 gene. We reconstructed the phylogeny of Neanderthal and modern haplotypes in this region and calculated the probability of the observed clustering being the result of introgression or common descent. We imputed variants for the participants in our previous genome-wide association study and we compared the distribution of Neanderthal variants between MN cases and controls. The region associated with the lead MN risk locus in the PLA2R1 gene was confirmed and showed that, within a 507 kb region enriched in introgressed sequence, a stringently defined 105 kb haplotype, intersecting the coding regions for PLA2R1 and ITGB6, is inherited from Neanderthals. Thus, introgressed Neanderthal haplotypes overlapping PLA2R1 are differentially represented in MN cases and controls, with enrichment In controls suggesting a protective effect.

The genomic landscape across 474 surgically accessible epileptogenic human brain lesions.

Understanding the exact molecular mechanisms involved in the aetiology of epileptogenic pathologies with or without tumour activity is essential for improving treatment of drug-resistant focal epilepsy. Here, we characterize the landscape of somatic genetic variants in resected brain specimens from 474 individuals with drug-resistant focal epilepsy using deep whole-exome sequencing (>350×) and whole-genome genotyping. Across the exome, we observe a greater number of somatic single-nucleotide variants in low-grade epilepsy-associated tumours (7.92 ± 5.65 single-nucleotide variants) than in brain tissue from malformations of cortical development (6.11 ± 4 single-nucleotide variants) or hippocampal sclerosis (5.1 ± 3.04 single-nucleotide variants). Tumour tissues also had the largest number of likely pathogenic variant carrying cells. low-grade epilepsy-associated tumours had the highest proportion of samples with one or more somatic copy-number variants (24.7%), followed by malformations of cortical development (5.4%) and hippocampal sclerosis (4.1%). Recurring somatic whole chromosome duplications affecting Chromosome 7 (16.8%), chromosome 5 (10.9%), and chromosome 20 (9.9%) were observed among low-grade epilepsy-associated tumours. For germline variant-associated malformations of cortical development genes such as TSC2, DEPDC5 and PTEN, germline single-nucleotide variants were frequently identified within large loss of heterozygosity regions, supporting the recently proposed 'second hit' disease mechanism in these genes. We detect somatic variants in 12 established lesional epilepsy genes and demonstrate exome-wide statistical support for three of these in the aetiology of low-grade epilepsy-associated tumours (e.g. BRAF) and malformations of cortical development (e.g. SLC35A2 and MTOR). We also identify novel significant associations for PTPN11 with low-grade epilepsy-associated tumours and NRAS Q61 mutated protein with a complex malformation of cortical development characterized by polymicrogyria and nodular heterotopia. The variants identified in NRAS are known from cancer studies to lead to hyperactivation of NRAS, which can be targeted pharmacologically. We identify large recurrent 1q21-q44 duplication including AKT3 in association with focal cortical dysplasia type 2a with hyaline astrocytic inclusions, another rare and possibly under-recognized brain lesion. The clinical-genetic analyses showed that the numbers of somatic single-nucleotide variant across the exome and the fraction of affected cells were positively correlated with the age at seizure onset and surgery in individuals with low-grade epilepsy-associated tumours. In summary, our comprehensive genetic screen sheds light on the genome-scale landscape of genetic variants in epileptic brain lesions, informs the design of gene panels for clinical diagnostic screening and guides future directions for clinical implementation of epilepsy surgery genetics.

Insights into clonal haematopoiesis from 8,342 mosaic chromosomal alterations.

The selective pressures that shape clonal evolution in healthy individuals are largely unknown. Here we investigate 8,342 mosaic chromosomal alterations, from 50 kb to 249 Mb long, that we uncovered in blood-derived DNA from 151,202 UK Biobank participants using phase-based computational techniques (estimated false discovery rate, 6-9%). We found six loci at which inherited variants associated strongly with the acquisition of deletions or loss of heterozygosity in cis. At three such loci (MPL, TM2D3-TARSL2, and FRA10B), we identified a likely causal variant that acted with high penetrance (5-50%). Inherited alleles at one locus appeared to affect the probability of somatic mutation, and at three other loci to be objects of positive or negative clonal selection. Several specific mosaic chromosomal alterations were strongly associated with future haematological malignancies. Our results reveal a multitude of paths towards clonal expansions with a wide range of effects on human health.

Chromosomal instability drives metastasis through a cytosolic DNA response.

Chromosomal instability is a hallmark of cancer that results from ongoing errors in chromosome segregation during mitosis. Although chromosomal instability is a major driver of tumour evolution, its role in metastasis has not been established. Here we show that chromosomal instability promotes metastasis by sustaining a tumour cell-autonomous response to cytosolic DNA. Errors in chromosome segregation create a preponderance of micronuclei whose rupture spills genomic DNA into the cytosol. This leads to the activation of the cGAS-STING (cyclic GMP-AMP synthase-stimulator of interferon genes) cytosolic DNA-sensing pathway and downstream noncanonical NF-κB signalling. Genetic suppression of chromosomal instability markedly delays metastasis even in highly aneuploid tumour models, whereas continuous chromosome segregation errors promote cellular invasion and metastasis in a STING-dependent manner. By subverting lethal epithelial responses to cytosolic DNA, chromosomally unstable tumour cells co-opt chronic activation of innate immune pathways to spread to distant organs.

Analysing an allelic series of rare missense variants of CACNA1I in a Swedish schizophrenia cohort.

CACNA1I is implicated in the susceptibility to schizophrenia by large-scale genetic association studies of single nucleotide polymorphisms. However, the channelopathy of CACNA1I in schizophrenia is unknown. CACNA1I encodes CaV3.3, a neuronal voltage-gated calcium channel that underlies a subtype of T-type current that is important for neuronal excitability in the thalamic reticular nucleus and other regions of the brain. Here, we present an extensive functional characterization of 57 naturally occurring rare and common missense variants of CACNA1I derived from a Swedish schizophrenia cohort of more than 10 000 individuals. Our analysis of this allelic series of coding CACNA1I variants revealed that reduced CaV3.3 channel current density was the dominant phenotype associated with rare CACNA1I coding alleles derived from control subjects, whereas rare CACNA1I alleles from schizophrenia patients encoded CaV3.3 channels with altered responses to voltages. CACNA1I variants associated with altered current density primarily impact the ionic channel pore and those associated with altered responses to voltage impact the voltage-sensing domain. CaV3.3 variants associated with altered voltage dependence of the CaV3.3 channel and those associated with peak current density deficits were significantly segregated across affected and unaffected groups (Fisher's exact test, P = 0.034). Our results, together with recent data from the SCHEMA (Schizophrenia Exome Sequencing Meta-Analysis) cohort, suggest that reduced CaV3.3 function may protect against schizophrenia risk in rare cases. We subsequently modelled the effect of the biophysical properties of CaV3.3 channel variants on thalamic reticular nucleus excitability and found that compared with common variants, ultrarare CaV3.3-coding variants derived from control subjects significantly decreased thalamic reticular nucleus excitability (P = 0.011). When all rare variants were analysed, there was a non-significant trend between variants that reduced thalamic reticular nucleus excitability and variants that either had no effect or increased thalamic reticular nucleus excitability across disease status. Taken together, the results of our functional analysis of an allelic series of >50 CACNA1I variants in a schizophrenia cohort reveal that loss of function of CaV3.3 is a molecular phenotype associated with reduced disease risk burden, and our approach may serve as a template strategy for channelopathies in polygenic disorders.

Human pluripotent stem cells recurrently acquire and expand dominant negative P53 mutations.

Human pluripotent stem cells (hPS cells) can self-renew indefinitely, making them an attractive source for regenerative therapies. This expansion potential has been linked with the acquisition of large copy number variants that provide mutated cells with a growth advantage in culture. The nature, extent and functional effects of other acquired genome sequence mutations in cultured hPS cells are not known. Here we sequence the protein-coding genes (exomes) of 140 independent human embryonic stem cell (hES cell) lines, including 26 lines prepared for potential clinical use. We then apply computational strategies for identifying mutations present in a subset of cells in each hES cell line. Although such mosaic mutations were generally rare, we identified five unrelated hES cell lines that carried six mutations in the TP53 gene that encodes the tumour suppressor P53. The TP53 mutations we observed are dominant negative and are the mutations most commonly seen in human cancers. We found that the TP53 mutant allelic fraction increased with passage number under standard culture conditions, suggesting that the P53 mutations confer selective advantage. We then mined published RNA sequencing data from 117 hPS cell lines, and observed another nine TP53 mutations, all resulting in coding changes in the DNA-binding domain of P53. In three lines, the allelic fraction exceeded 50%, suggesting additional selective advantage resulting from the loss of heterozygosity at the TP53 locus. As the acquisition and expansion of cancer-associated mutations in hPS cells may go unnoticed during most applications, we suggest that careful genetic characterization of hPS cells and their differentiated derivatives be carried out before clinical use.

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder.

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

Schizophrenia risk from complex variation of complement component 4.

Schizophrenia is a heritable brain illness with unknown pathogenic mechanisms. Schizophrenia's strongest genetic association at a population level involves variation in the major histocompatibility complex (MHC) locus, but the genes and molecular mechanisms accounting for this have been challenging to identify. Here we show that this association arises in part from many structurally diverse alleles of the complement component 4 (C4) genes. We found that these alleles generated widely varying levels of C4A and C4B expression in the brain, with each common C4 allele associating with schizophrenia in proportion to its tendency to generate greater expression of C4A. Human C4 protein localized to neuronal synapses, dendrites, axons, and cell bodies. In mice, C4 mediated synapse elimination during postnatal development. These results implicate excessive complement activity in the development of schizophrenia and may help explain the reduced numbers of synapses in the brains of individuals with schizophrenia.

Haplotype Sharing Provides Insights into Fine-Scale Population History and Disease in Finland.

Finland provides unique opportunities to investigate population and medical genomics because of its adoption of unified national electronic health records, detailed historical and birth records, and serial population bottlenecks. We assembled a comprehensive view of recent population history (≤100 generations), the timespan during which most rare-disease-causing alleles arose, by comparing pairwise haplotype sharing from 43,254 Finns to that of 16,060 Swedes, Estonians, Russians, and Hungarians from geographically and linguistically adjacent countries with different population histories. We find much more extensive sharing in Finns, with at least one ≥ 5 cM tract on average between pairs of unrelated individuals. By coupling haplotype sharing with fine-scale birth records from more than 25,000 individuals, we find that although haplotype sharing broadly decays with geographical distance, there are pockets of excess haplotype sharing; individuals from northeast Finland typically share several-fold more of their genome in identity-by-descent segments than individuals from southwest regions. We estimate recent effective population-size changes through time across regions of Finland, and we find that there was more continuous gene flow as Finns migrated from southwest to northeast between the early- and late-settlement regions than was dichotomously described previously. Lastly, we show that haplotype sharing is locally enriched by an order of magnitude among pairs of individuals sharing rare alleles and especially among pairs sharing rare disease-causing variants. Our work provides a general framework for using haplotype sharing to reconstruct an integrative view of recent population history and gain insight into the evolutionary origins of rare variants contributing to disease.

The Genetic Landscape of Diamond-Blackfan Anemia.

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 7 out of 1,000,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole-exome sequencing (WES). We identified relevant rare and predicted damaging mutations for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and located in 1 of 19 previously reported ribosomal protein (RP)-encoding genes. Using exon coverage estimates, we identified and validated 31 deletions in RP genes. We also observed an enrichment for extended splice site mutations and validated their diverse effects using RNA sequencing in cell lines obtained from individuals with DBA. Leveraging the size of our cohort, we observed robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. We further identified rare mutations in seven previously unreported RP genes that may cause DBA, as well as several distinct disorders that appear to phenocopy DBA, including nine individuals with biallelic CECR1 mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain >5% of DBA-affected case subjects. Overall, this report should inform not only clinical practice for DBA-affected individuals, but also the design and analysis of rare variant studies for heterogeneous Mendelian disorders.