It has recently been shown that ancestors of New Guineans and Bougainville Islanders have inherited a proportion of their ancestry from Denisovans, an archaic hominin group from Siberia. However, only a sparse sampling of populations from Southeast Asia and Oceania were analyzed. Here, we quantify Denisova admixture in 33 additional populations from Asia and Oceania. Aboriginal Australians, Near Oceanians, Polynesians, Fijians, east Indonesians, and Mamanwa (a "Negrito" group from the Philippines) have all inherited genetic material from Denisovans, but mainland East Asians, western Indonesians, Jehai (a Negrito group from Malaysia), and Onge (a Negrito group from the Andaman Islands) have not. These results indicate that Denisova gene flow occurred into the common ancestors of New Guineans, Australians, and Mamanwa but not into the ancestors of the Jehai and Onge and suggest that relatives of present-day East Asians were not in Southeast Asia when the Denisova gene flow occurred. Our finding that descendants of the earliest inhabitants of Southeast Asia do not all harbor Denisova admixture is inconsistent with a history in which the Denisova interbreeding occurred in mainland Asia and then spread over Southeast Asia, leading to all its earliest modern human inhabitants. Instead, the data can be most parsimoniously explained if the Denisova gene flow occurred in Southeast Asia itself. Thus, archaic Denisovans must have lived over an extraordinarily broad geographic and ecological range, from Siberia to tropical Asia.
Whole-genome sequencing across multiple samples in a population provides an unprecedented opportunity for comprehensively characterizing the polymorphic variants in the population. Although the 1000 Genomes Project (1KGP) has offered brief insights into the value of population-level sequencing, the low coverage has compromised the ability to confidently detect rare and low-frequency variants. In addition, the composition of populations in the 1KGP is not complete, despite the fact that the study design has been extended to more than 2,500 samples from more than 20 population groups. The Malays are one of the Austronesian groups predominantly present in Southeast Asia and Oceania, and the Singapore Sequencing Malay Project (SSMP) aims to perform deep whole-genome sequencing of 100 healthy Malays. By sequencing at a minimum of 30× coverage, we have illustrated the higher sensitivity at detecting low-frequency and rare variants and the ability to investigate the presence of hotspots of functional mutations. Compared to the low-pass sequencing in the 1KGP, the deeper coverage allows more functional variants to be identified for each person. A comparison of the fidelity of genotype imputation of Malays indicated that a population-specific reference panel, such as the SSMP, outperforms a cosmopolitan panel with larger number of individuals for common SNPs. For lower-frequency (<5%) markers, a larger number of individuals might have to be whole-genome sequenced so that the accuracy currently afforded by the 1KGP can be achieved. The SSMP data are expected to be the benchmark for evaluating the value of deep population-level sequencing versus low-pass sequencing, especially in populations that are poorly represented in population-genetics studies.
Tibetan high-altitude adaptation (HAA) has been studied extensively, and many candidate genes have been reported. Subsequent efforts targeting HAA functional variants, however, have not been that successful (e.g., no functional variant has been suggested for the top candidate HAA gene, EPAS1). With WinXPCNVer, a method developed in this study, we detected in microarray data a Tibetan-enriched deletion (TED) carried by 90% of Tibetans; 50% were homozygous for the deletion, whereas only 3% carried the TED and 0% carried the homozygous deletion in 2,792 worldwide samples (p < 10(-15)). We employed long PCR and Sanger sequencing technologies to determine the exact copy number and breakpoints of the TED in 70 additional Tibetan and 182 diverse samples. The TED had identical boundaries (chr2: 46,694,276-46,697,683; hg19) and was 80 kb downstream of EPAS1. Notably, the TED was in strong linkage disequilibrium (LD; r(2) = 0.8) with EPAS1 variants associated with reduced blood concentrations of hemoglobin. It was also in complete LD with the 5-SNP motif, which was suspected to be introgressed from Denisovans, but the deletion itself was absent from the Denisovan sequence. Correspondingly, we detected that footprints of positive selection for the TED occurred 12,803 (95% confidence interval = 12,075-14,725) years ago. We further whole-genome deep sequenced (>60×) seven Tibetans and verified the TED but failed to identify any other copy-number variations with comparable patterns, giving this TED top priority for further study. We speculate that the specific patterns of the TED resulted from its own functionality in HAA of Tibetans or LD with a functional variant of EPAS1.
The jungle habitat of the Temuan aborigines harbors a variety of infectious diseases, the most notable being malaria. Our study of 15 genetic systems in the Temuan revealed substantial polymorphism and within-population genetic diversity. The polymorphisms for Hb beta, G6PD, and El are of interest in regard to genetic adaptation to malaria. Among the polymorphisms investigated we conclude that G6PD deficiency and elliptocytosis are likely to have malaria-resistant effects as evidenced by their low association with malarial parasitemia or their higher frequency in adults than in children. These findings suggest that the malarial habitat of the Temuans is livable in the long range sense for them because of the cluster of malaria-resistant alleles in their gene pool (G6PD)-, El, and possibly, but not tested here because of its low frequency, Hb beta E). The same condition probably holds for the Semai, the nearest aborigine neighbors of the Temuan (although the Semai have not been tested for malarial parasitemia and for these polymorphisms simultaneously), since the Semai have substantial Hb betaE, G6PD-, and El. The Temuan have a cultural identity system of rituals, beliefs, and certain aspects of language which effectively isolates them genetically from Malays and other nonaborigines. This system hinders the dilution of the malaria-resistant alleles of the Temuan gene pool with the malaria-susceptible alleles of the nonaborigine gene pools.
De novo variants (DNVs) cause many genetic diseases. When DNVs are examined in the whole coding regions of genes in next-generation sequencing analyses, pathogenic DNVs often cluster in a specific region. One such region is the last exon and the last 50 bp of the penultimate exon, where truncating DNVs cause escape from nonsense-mediated mRNA decay [NMD(-) region]. Such variants can have dominant-negative or gain-of-function effects. Here, we first developed a resource of rates of truncating DNVs in NMD(-) regions under the null model of DNVs. Utilizing this resource, we performed enrichment analysis of truncating DNVs in NMD(-) regions in 346 developmental and epileptic encephalopathy (DEE) trios. We observed statistically significant enrichment of truncating DNVs in semaphorin 6B (SEMA6B) (p value: 2.8 × 10-8; exome-wide threshold: 2.5 × 10-6). The initial analysis of the 346 individuals and additional screening of 1,406 and 4,293 independent individuals affected by DEE and developmental disorders collectively identified four truncating DNVs in the SEMA6B NMD(-) region in five individuals who came from unrelated families (p value: 1.9 × 10-13) and consistently showed progressive myoclonic epilepsy. RNA analysis of lymphoblastoid cells established from an affected individual showed that the mutant allele escaped NMD, indicating stable production of the truncated protein. Importantly, heterozygous truncating variants in the NMD(+) region of SEMA6B are observed in general populations, and SEMA6B is most likely loss-of-function tolerant. Zebrafish expressing truncating variants in the NMD(-) region of SEMA6B orthologs displayed defective development of brain neurons and enhanced pentylenetetrazole-induced seizure behavior. In summary, we show that truncating DNVs in the final exon of SEMA6B cause progressive myoclonic epilepsy.
The identification of disease alleles underlying human autoinflammatory diseases can provide important insights into the mechanisms that maintain neutrophil homeostasis. Here, we focused our attention on generalized pustular psoriasis (GPP), a potentially life-threatening disorder presenting with cutaneous and systemic neutrophilia. Following the whole-exome sequencing of 19 unrelated affected individuals, we identified a subject harboring a homozygous splice-site mutation (c.2031-2A>C) in MPO. This encodes myeloperoxidase, an essential component of neutrophil azurophil granules. MPO screening in conditions phenotypically related to GPP uncovered further disease alleles in one subject with acral pustular psoriasis (c.2031-2A>C;c.2031-2A>C) and in two individuals with acute generalized exanthematous pustulosis (c.1705C>T;c.2031-2A>C and c.1552_1565del;c.1552_1565del). A subsequent analysis of UK Biobank data demonstrated that the c.2031-2A>C and c.1705C>T (p.Arg569Trp) disease alleles were also associated with increased neutrophil abundance in the general population (p = 5.1 × 10-6 and p = 3.6 × 10-5, respectively). The same applied to three further deleterious variants that had been genotyped in the cohort, with two alleles (c.995C>T [p.Ala332Val] and c.752T>C [p.Met251Thr]) yielding p values < 10-10. Finally, treatment of healthy neutrophils with an MPO inhibitor (4-Aminobenzoic acid hydrazide) increased cell viability and delayed apoptosis, highlighting a mechanism whereby MPO mutations affect granulocyte numbers. These findings identify MPO as a genetic determinant of pustular skin disease and neutrophil abundance. Given the recent interest in the development of MPO antagonists for the treatment of neurodegenerative disease, our results also suggest that the pro-inflammatory effects of these agents should be closely monitored.
We ascertained a multi-generation Malaysian family with Joubert syndrome (JS). The presence of asymptomatic obligate carrier females suggested an X-linked recessive inheritance pattern. Affected males presented with mental retardation accompanied by postaxial polydactyly and retinitis pigmentosa. Brain MRIs showed the presence of a "molar tooth sign," which classifies this syndrome as classic JS with retinal involvement. Linkage analysis showed linkage to Xpter-Xp22.2 and a maximum LOD score of 2.06 for marker DXS8022. Mutation analysis revealed a frameshift mutation, p.K948NfsX8, in exon 21 of OFD1. In an isolated male with JS, a second frameshift mutation, p.E923KfsX3, in the same exon was identified. OFD1 has previously been associated with oral-facial-digital type 1 (OFD1) syndrome, a male-lethal X-linked dominant condition, and with X-linked recessive Simpson-Golabi-Behmel syndrome type 2 (SGBS2). In a yeast two-hybrid screen of a retinal cDNA library, we identified OFD1 as an interacting partner of the LCA5-encoded ciliary protein lebercilin. We show that X-linked recessive mutations in OFD1 reduce, but do not eliminate, the interaction with lebercilin, whereas X-linked dominant OFD1 mutations completely abolish binding to lebercilin. In addition, recessive mutations in OFD1 did not affect the pericentriolar localization of the recombinant protein in hTERT-RPE1 cells, whereas this localization was lost for dominant mutations. These findings offer a molecular explanation for the phenotypic spectrum observed for OFD1 mutations; this spectrum now includes OFD1 syndrome, SGBS2, and JS.
Previous research has shown that polygenic risk scores (PRSs) can be used to stratify women according to their risk of developing primary invasive breast cancer. This study aimed to evaluate the association between a recently validated PRS of 313 germline variants (PRS313) and contralateral breast cancer (CBC) risk. We included 56,068 women of European ancestry diagnosed with first invasive breast cancer from 1990 onward with follow-up from the Breast Cancer Association Consortium. Metachronous CBC risk (N = 1,027) according to the distribution of PRS313 was quantified using Cox regression analyses. We assessed PRS313 interaction with age at first diagnosis, family history, morphology, ER status, PR status, and HER2 status, and (neo)adjuvant therapy. In studies of Asian women, with limited follow-up, CBC risk associated with PRS313 was assessed using logistic regression for 340 women with CBC compared with 12,133 women with unilateral breast cancer. Higher PRS313 was associated with increased CBC risk: hazard ratio per standard deviation (SD) = 1.25 (95%CI = 1.18-1.33) for Europeans, and an OR per SD = 1.15 (95%CI = 1.02-1.29) for Asians. The absolute lifetime risks of CBC, accounting for death as competing risk, were 12.4% for European women at the 10th percentile and 20.5% at the 90th percentile of PRS313. We found no evidence of confounding by or interaction with individual characteristics, characteristics of the primary tumor, or treatment. The C-index for the PRS313 alone was 0.563 (95%CI = 0.547-0.586). In conclusion, PRS313 is an independent factor associated with CBC risk and can be incorporated into CBC risk prediction models to help improve stratification and optimize surveillance and treatment strategies.
Whereas large-scale statistical analyses can robustly identify disease-gene relationships, they do not accurately capture genotype-phenotype correlations or disease mechanisms. We use multiple lines of independent evidence to show that different variant types in a single gene, SATB1, cause clinically overlapping but distinct neurodevelopmental disorders. Clinical evaluation of 42 individuals carrying SATB1 variants identified overt genotype-phenotype relationships, associated with different pathophysiological mechanisms, established by functional assays. Missense variants in the CUT1 and CUT2 DNA-binding domains result in stronger chromatin binding, increased transcriptional repression, and a severe phenotype. In contrast, variants predicted to result in haploinsufficiency are associated with a milder clinical presentation. A similarly mild phenotype is observed for individuals with premature protein truncating variants that escape nonsense-mediated decay, which are transcriptionally active but mislocalized in the cell. Our results suggest that in-depth mutation-specific genotype-phenotype studies are essential to capture full disease complexity and to explain phenotypic variability.
Progressive myoclonus epilepsies (PMEs) comprise a group of clinically and genetically heterogeneous rare diseases. Over 70% of PME cases can now be molecularly solved. Known PME genes encode a variety of proteins, many involved in lysosomal and endosomal function. We performed whole-exome sequencing (WES) in 84 (78 unrelated) unsolved PME-affected individuals, with or without additional family members, to discover novel causes. We identified likely disease-causing variants in 24 out of 78 (31%) unrelated individuals, despite previous genetic analyses. The diagnostic yield was significantly higher for individuals studied as trios or families (14/28) versus singletons (10/50) (OR = 3.9, p value = 0.01, Fisher's exact test). The 24 likely solved cases of PME involved 18 genes. First, we found and functionally validated five heterozygous variants in NUS1 and DHDDS and a homozygous variant in ALG10, with no previous disease associations. All three genes are involved in dolichol-dependent protein glycosylation, a pathway not previously implicated in PME. Second, we independently validate SEMA6B as a dominant PME gene in two unrelated individuals. Third, in five families, we identified variants in established PME genes; three with intronic or copy-number changes (CLN6, GBA, NEU1) and two very rare causes (ASAH1, CERS1). Fourth, we found a group of genes usually associated with developmental and epileptic encephalopathies, but here, remarkably, presenting as PME, with or without prior developmental delay. Our systematic analysis of these cases suggests that the small residuum of unsolved cases will most likely be a collection of very rare, genetically heterogeneous etiologies.
We describe four families with affected siblings showing unique clinical features: early-onset (before 1 year of age) progressive diffuse brain atrophy with regression, postnatal microcephaly, postnatal growth retardation, muscle weakness/atrophy, and respiratory failure. By whole-exome sequencing, we identified biallelic TBCD mutations in eight affected individuals from the four families. TBCD encodes TBCD (tubulin folding co-factor D), which is one of five tubulin-specific chaperones playing a pivotal role in microtubule assembly in all cells. A total of seven mutations were found: five missense mutations, one nonsense, and one splice site mutation resulting in a frameshift. In vitro cell experiments revealed the impaired binding between most mutant TBCD proteins and ARL2, TBCE, and β-tubulin. The in vivo experiments using olfactory projection neurons in Drosophila melanogaster indicated that the TBCD mutations caused loss of function. The wide range of clinical severity seen in this neurodegenerative encephalopathy may result from the residual function of mutant TBCD proteins. Furthermore, the autopsied brain from one deceased individual showed characteristic neurodegenerative findings: cactus and somatic sprout formations in the residual Purkinje cells in the cerebellum, which are also seen in some diseases associated with mitochondrial impairment. Defects of microtubule formation caused by TBCD mutations may underlie the pathomechanism of this neurodegenerative encephalopathy.
The ABO, MN and Rh blood types, and the Hp, Tf, and Gm [Gm (a), Gm (x), Gm(b), and Gm-like] factors were determined for 128 unrelated Indians (parents of families, 63 with two parents tested and two with one parent tested), and 90 unrelated Chinese (parents of 46 families, 44 with two parents tested and two with one parent tested), and for the offspring from these families. The frequencies of the several blood types are presented. They were done primarily to aid in paternity testing. They compare favorably with the findings of previous studies. The allele Hp1 is rare in the Indian population (.09) and relatively infrequent in the Chinese (.29). Unfortunately, the data shed no light on the problem of the inheritance of the phenotype Hp O. Only Tf C was found among the Indians. About four per cent of the Chinese were heterozygous for Tf CD,, all other were Tf CC. The Indians have a high frequency of Gm(a) and of Gm (x), and a low frequency of Gm (b). They appear to have alleles Gma, Gmax, and Gmb in the following frequencies: .535, .234(5), and .230(5), respectively. Three families appear to have a GMxb allele, providing the offspring are not extra-marital. The Chinese appear to have the alleles Gm^ab, Gm^a, and Gm^ax in the following frequencies: .741, .231, and .028, respectively.
Genetic risk prediction has several potential applications in medical research and clinical practice and could be used, for example, to stratify a heterogeneous population of patients by their predicted genetic risk. However, for polygenic traits, such as psychiatric disorders, the accuracy of risk prediction is low. Here we use a multivariate linear mixed model and apply multi-trait genomic best linear unbiased prediction for genetic risk prediction. This method exploits correlations between disorders and simultaneously evaluates individual risk for each disorder. We show that the multivariate approach significantly increases the prediction accuracy for schizophrenia, bipolar disorder, and major depressive disorder in the discovery as well as in independent validation datasets. By grouping SNPs based on genome annotation and fitting multiple random effects, we show that the prediction accuracy could be further improved. The gain in prediction accuracy of the multivariate approach is equivalent to an increase in sample size of 34% for schizophrenia, 68% for bipolar disorder, and 76% for major depressive disorders using single trait models. Because our approach can be readily applied to any number of GWAS datasets of correlated traits, it is a flexible and powerful tool to maximize prediction accuracy. With current sample size, risk predictors are not useful in a clinical setting but already are a valuable research tool, for example in experimental designs comparing cases with high and low polygenic risk.
RAC1 is a widely studied Rho GTPase, a class of molecules that modulate numerous cellular functions essential for normal development. RAC1 is highly conserved across species and is under strict mutational constraint. We report seven individuals with distinct de novo missense RAC1 mutations and varying degrees of developmental delay, brain malformations, and additional phenotypes. Four individuals, each harboring one of c.53G>A (p.Cys18Tyr), c.116A>G (p.Asn39Ser), c.218C>T (p.Pro73Leu), and c.470G>A (p.Cys157Tyr) variants, were microcephalic, with head circumferences between -2.5 to -5 SD. In contrast, two individuals with c.151G>A (p.Val51Met) and c.151G>C (p.Val51Leu) alleles were macrocephalic with head circumferences of +4.16 and +4.5 SD. One individual harboring a c.190T>G (p.Tyr64Asp) allele had head circumference in the normal range. Collectively, we observed an extraordinary spread of ∼10 SD of head circumferences orchestrated by distinct mutations in the same gene. In silico modeling, mouse fibroblasts spreading assays, and in vivo overexpression assays using zebrafish as a surrogate model demonstrated that the p.Cys18Tyr and p.Asn39Ser RAC1 variants function as dominant-negative alleles and result in microcephaly, reduced neuronal proliferation, and cerebellar abnormalities in vivo. Conversely, the p.Tyr64Asp substitution is constitutively active. The remaining mutations are probably weakly dominant negative or their effects are context dependent. These findings highlight the importance of RAC1 in neuronal development. Along with TRIO and HACE1, a sub-category of rare developmental disorders is emerging with RAC1 as the central player. We show that ultra-rare disorders caused by private, non-recurrent missense mutations that result in varying phenotypes are challenging to dissect, but can be delineated through focused international collaboration.
By combining data from 160,500 individuals with breast cancer and 226,196 controls of Asian and European ancestry, we conducted genome- and transcriptome-wide association studies of breast cancer. We identified 222 genetic risk loci and 137 genes that were associated with breast cancer risk at a p
In polygenic score (PGS) analysis, the coefficient of determination (R2) is a key statistic to evaluate efficacy. R2 is the proportion of phenotypic variance explained by the PGS, calculated in a cohort that is independent of the genome-wide association study (GWAS) that provided estimates of allelic effect sizes. The SNP-based heritability (hSNP2, the proportion of total phenotypic variances attributable to all common SNPs) is the theoretical upper limit of the out-of-sample prediction R2. However, in real data analyses R2 has been reported to exceed hSNP2, which occurs in parallel with the observation that hSNP2 estimates tend to decline as the number of cohorts being meta-analyzed increases. Here, we quantify why and when these observations are expected. Using theory and simulation, we show that if heterogeneities in cohort-specific hSNP2 exist, or if genetic correlations between cohorts are less than one, hSNP2 estimates can decrease as the number of cohorts being meta-analyzed increases. We derive conditions when the out-of-sample prediction R2 will be greater than hSNP2 and show the validity of our derivations with real data from a binary trait (major depression) and a continuous trait (educational attainment). Our research calls for a better approach to integrating information from multiple cohorts to address issues of between-cohort heterogeneity.
Leucine zipper-EF-hand containing transmembrane protein 1 (LETM1) encodes an inner mitochondrial membrane protein with an osmoregulatory function controlling mitochondrial volume and ion homeostasis. The putative association of LETM1 with a human disease was initially suggested in Wolf-Hirschhorn syndrome, a disorder that results from de novo monoallelic deletion of chromosome 4p16.3, a region encompassing LETM1. Utilizing exome sequencing and international gene-matching efforts, we have identified 18 affected individuals from 11 unrelated families harboring ultra-rare bi-allelic missense and loss-of-function LETM1 variants and clinical presentations highly suggestive of mitochondrial disease. These manifested as a spectrum of predominantly infantile-onset (14/18, 78%) and variably progressive neurological, metabolic, and dysmorphic symptoms, plus multiple organ dysfunction associated with neurodegeneration. The common features included respiratory chain complex deficiencies (100%), global developmental delay (94%), optic atrophy (83%), sensorineural hearing loss (78%), and cerebellar ataxia (78%) followed by epilepsy (67%), spasticity (53%), and myopathy (50%). Other features included bilateral cataracts (42%), cardiomyopathy (36%), and diabetes (27%). To better understand the pathogenic mechanism of the identified LETM1 variants, we performed biochemical and morphological studies on mitochondrial K+/H+ exchange activity, proteins, and shape in proband-derived fibroblasts and muscles and in Saccharomyces cerevisiae, which is an important model organism for mitochondrial osmotic regulation. Our results demonstrate that bi-allelic LETM1 variants are associated with defective mitochondrial K+ efflux, swollen mitochondrial matrix structures, and loss of important mitochondrial oxidative phosphorylation protein components, thus highlighting the implication of perturbed mitochondrial osmoregulation caused by LETM1 variants in neurological and mitochondrial pathologies.
Pathogenic variants in multiple genes on the X chromosome have been implicated in syndromic and non-syndromic intellectual disability disorders. ZFX on Xp22.11 encodes a transcription factor that has been linked to diverse processes including oncogenesis and development, but germline variants have not been characterized in association with disease. Here, we present clinical and molecular characterization of 18 individuals with germline ZFX variants. Exome or genome sequencing revealed 11 variants in 18 subjects (14 males and 4 females) from 16 unrelated families. Four missense variants were identified in 11 subjects, with seven truncation variants in the remaining individuals. Clinical findings included developmental delay/intellectual disability, behavioral abnormalities, hypotonia, and congenital anomalies. Overlapping and recurrent facial features were identified in all subjects, including thickening and medial broadening of eyebrows, variations in the shape of the face, external eye abnormalities, smooth and/or long philtrum, and ear abnormalities. Hyperparathyroidism was found in four families with missense variants, and enrichment of different tumor types was observed. In molecular studies, DNA-binding domain variants elicited differential expression of a small set of target genes relative to wild-type ZFX in cultured cells, suggesting a gain or loss of transcriptional activity. Additionally, a zebrafish model of ZFX loss displayed an altered behavioral phenotype, providing additional evidence for the functional significance of ZFX. Our clinical and experimental data support that variants in ZFX are associated with an X-linked intellectual disability syndrome characterized by a recurrent facial gestalt, neurocognitive and behavioral abnormalities, and an increased risk for congenital anomalies and hyperparathyroidism.