American journal of human genetics最新文献

Whereas 16p11.2 BP4-5 copy-number variants (CNVs) represent one of the most pleiotropic etiologies of genomic syndromes in both clinical and population cohorts, the mechanisms leading to such pleiotropy remain understudied. Identifying 73 deletion and 89 duplication carrier individuals among unrelated White British UK Biobank participants, we performed a phenome-wide association study (PheWAS) between the region's copy number and 117 complex traits and diseases, mimicking four dosage models. Forty-six phenotypes (39%) were affected by 16p11.2 BP4-5 CNVs, with the deletion-only, mirror, U-shape, and duplication-only models being the best fit for 30, 10, 4, and 2 phenotypes, respectively, aligning with the stronger deleteriousness of the deletion. Upon individually adjusting CNV effects for either body mass index (BMI), height, or educational attainment (EA), we found that sixteen testable deletion-driven associations-primarily with cardiovascular and metabolic traits-were BMI dependent, with EA playing a more subtle role and no association depending on height. Bidirectional Mendelian randomization supported that 13 out of these 16 associations were secondary consequences of the CNV's impact on BMI. For the 23 traits that remained significantly associated upon individual adjustment for mediators, matched-control analyses found that 10 phenotypes, including musculoskeletal traits, liver enzymes, fluid intelligence, platelet count, and pneumonia and acute kidney injury risk, remained associated under strict Bonferroni correction, with 10 additional nominally significant associations. These results paint a complex picture of 16p11.2 BP4-5's pleiotropic pattern that involves direct effects on multiple physiological systems and indirect co-morbidities consequential to the CNV's impact on BMI and EA, acting through trait-specific dosage mechanisms.

Microtubule affinity-regulating kinase 2 (MARK2) contributes to establishing neuronal polarity and developing dendritic spines. Although large-scale sequencing studies have associated MARK2 variants with autism spectrum disorder (ASD), the clinical features and variant spectrum in affected individuals with MARK2 variants, early developmental phenotypes in mutant human neurons, and the pathogenic mechanism underlying effects on neuronal development have remained unclear. Here, we report 31 individuals with MARK2 variants and presenting with ASD, other neurodevelopmental disorders, and distinctive facial features. Loss-of-function (LoF) variants predominate (81%) in affected individuals, while computational analysis and in vitro expression assay of missense variants supported the effect of MARK2 loss. Using proband-derived and CRISPR-engineered isogenic induced pluripotent stem cells (iPSCs), we show that MARK2 loss leads to early neuronal developmental and functional deficits, including anomalous polarity and dis-organization in neural rosettes, as well as imbalanced proliferation and differentiation in neural progenitor cells (NPCs). Mark2^+/- mice showed abnormal cortical formation and partition and ASD-like behavior. Through the use of RNA sequencing (RNA-seq) and lithium treatment, we link MARK2 loss to downregulation of the WNT/β-catenin signaling pathway and identify lithium as a potential drug for treating MARK2-associated ASD.

The ClinGen Hereditary Breast, Ovarian, and Pancreatic Cancer (HBOP) Variant Curation Expert Panel (VCEP) is composed of internationally recognized experts in clinical genetics, molecular biology, and variant interpretation. This VCEP made specifications for the American College of Medical Genetics and Association for Molecular Pathology (ACMG/AMP) guidelines for the ataxia telangiectasia mutated (ATM) gene according to the ClinGen protocol. These gene-specific rules for ATM were modified from the ACMG/AMP guidelines and were tested against 33 ATM variants of various types and classifications in a pilot curation phase. The pilot revealed a majority agreement between the HBOP VCEP classifications and the ClinVar-deposited classifications. Six pilot variants had conflicting interpretations in ClinVar, and re-evaluation with the VCEP's ATM-specific rules resulted in four that were classified as benign, one as likely pathogenic, and one as a variant of uncertain significance (VUS) by the VCEP, improving the certainty of interpretations in the public domain. Overall, 28 of the 33 pilot variants were not VUS, leading to an 85% classification rate. The ClinGen-approved, modified rules demonstrated value for improved interpretation of variants in ATM.

As more patients receive genome-wide sequencing, the number of individuals diagnosed with multiple monogenic conditions is increasing. We sought to investigate the relative phenotypic contribution of dual diagnoses using both manual curation and computational approaches. First, we computed 1,003,236 semantic similarity scores for all possible pairs of 1,417 genes in the Developmental Disorder Gene2Phenotype (DDG2P) database using Human Phenotype Ontology terms. Next, for 62 probands with two molecular diagnoses in the Deciphering Developmental Disorders study, we computed semantic similarity scores between the probands' phenotypes and DDG2P phenotypes associated with the two disorders and compared the results with manual attribution of proband phenotypes to none, one, or both of the genes. We found a spectrum of phenotypic similarity for dual diagnoses, both across all DDG2P genes and within dual diagnosed probands, from phenotypically distinct through blended to indistinguishable conditions. Pairwise semantic similarity scores between two DDG2P genes were a good predictor of the extent of phenotypic blending observed in probands. Dual diagnoses involving genes linked with synergistic phenotypes can result in more extreme presentations while those involving antagonistic phenotypes have spuriously high pairwise semantic similarity scores despite a potentially milder atypical presentation. We suggest that the phenotypic contribution of two molecular diagnoses may contain discrete, synergistic, or antagonistic elements. Conceptual recognition of this phenotypic spectrum is important for making a final clinico-molecular diagnosis and providing accurate genetic counseling.

Recent positive selection can result in an excess of long identity-by-descent (IBD) haplotype segments overlapping a locus. The statistical methods that we propose here address three major objectives in studying selective sweeps: scanning for regions of interest, identifying possible sweeping alleles, and estimating a selection coefficient s. First, we implement a selection scan to locate regions with excess IBD rates. Second, we estimate the allele frequency and location of an unknown sweeping allele by aggregating over variants that are more abundant in an inferred outgroup with excess IBD rate versus the rest of the sample. Third, we propose an estimator for the selection coefficient and quantify uncertainty using the parametric bootstrap. Comparing against state-of-the-art methods in extensive simulations, we show that our methods are more precise at estimating s when s≥0.015. We also show that our 95% confidence intervals contain s in nearly 95% of our simulations. We apply these methods to study positive selection in European ancestry samples from the Trans-Omics for Precision Medicine project. We analyze eight loci where IBD rates are more than four standard deviations above the genome-wide median, including LCT where the maximum IBD rate is 35 standard deviations above the genome-wide median. Overall, we present robust and accurate approaches to study recent adaptive evolution without knowing the identity of the causal allele or using time series data.

WD repeat domain 83 opposite strand (WDR83OS) encodes the 106-aa (amino acid) protein Asterix, which heterodimerizes with CCDC47 to form the PAT (protein associated with ER translocon) complex. This complex functions as a chaperone for large proteins containing transmembrane domains to ensure proper folding. Until recently, little was known about the role of WDR83OS or CCDC47 in human disease traits. However, biallelic variants in CCDC47 were identified in four unrelated families with trichohepatoneurodevelopmental syndrome, characterized by a neurodevelopmental disorder (NDD) with liver dysfunction. Three affected siblings in an additional family share a homozygous truncating WDR83OS variant and a phenotype of NDD, dysmorphic features, and liver dysfunction. Using family-based rare variant analyses of exome sequencing (ES) data and case matching through GeneMatcher, we describe the clinical phenotypes of 11 additional individuals in eight unrelated families (nine unrelated families, 14 individuals in total) with biallelic putative truncating variants in WDR83OS. Consistent clinical features include NDD (14/14), facial dysmorphism (13/14), intractable itching (9/14), and elevated bile acids (5/6). Whereas bile acids were significantly elevated in 5/6 of individuals tested, bilirubin was normal and liver enzymes were normal to mildly elevated in all 14 individuals. In three of six individuals for whom longitudinal data were available, we observed a progressive reduction in relative head circumference. A zebrafish model lacking Wdr83os function further supports its role in the nervous system, craniofacial development, and lipid absorption. Taken together, our data support a disease-gene association between biallelic loss-of-function of WDR83OS and a neurological disease trait with hypercholanemia.

A recent publication describing the assembly of the Y chromosomes of 43 males was remarkable not only for its ambitious technical scope but also for the startling suggestion that the boundary of the pseudoautosomal region 1 (PAR1), where the human X and Y chromosomes engage in crossing-over during male meiosis, lies 500 kb distal to its previously reported location. Where is the boundary of the human PAR1? We first review the evidence that mapped the PAR boundary, or PAB, before the human genome draft sequence was produced, then examine post-genomic datasets for evidence of crossing-over between the X and Y, and lastly re-examine contiguous sequence assemblies of the PAR-NPY boundary to see whether they support a more distal PAB. We find ample evidence of X-Y crossovers throughout the 500 kb in question, some as close as 246 bp to the previously reported PAB. Our new analyses, combined with previous studies over the past 40 years, provide overwhelming evidence to support the original position and narrow the probable location of the PAB to a 201-bp window.

Large biobank samples provide an opportunity to integrate broad phenotyping, familial records, and molecular genetics data to study complex traits and diseases. We introduce Pearson-Aitken Family Genetic Risk Scores (PA-FGRS), a method for estimating disease liability from patterns of diagnoses in extended, age-censored genealogical records. We then apply the method to study a paradigmatic complex disorder, major depressive disorder (MDD), using the iPSYCH2015 case-cohort study of 30,949 MDD cases, 39,655 random population controls, and more than 2 million relatives. We show that combining PA-FGRS liabilities estimated from family records with molecular genotypes of probands improves three lines of inquiry. Incorporating PA-FGRS liabilities improves classification of MDD over and above polygenic scores, identifies robust genetic contributions to clinical heterogeneity in MDD associated with comorbidity, recurrence, and severity and can improve the power of genome-wide association studies. Our method is flexible and easy to use, and our study approaches are generalizable to other datasets and other complex traits and diseases.

Mechanisms of abdominal obesity GWAS variants have remained largely unknown. To elucidate these mechanisms, we leveraged subcutaneous adipose tissue (SAT) single nucleus RNA-sequencing and genomics data. After discovering that heritability of abdominal obesity is enriched in adipocytes, we focused on a SAT unique adipocyte marker gene, the transcription factor TBX15, and its abdominal obesity-associated deleterious missense variant, rs10494217. The allele frequency of rs10494217 revealed a north-to-south decreasing gradient, with consistent significant F_ST values observed for 25 different populations when compared to Finns, a population with a history of genetic isolation. Given the role of Tbx15 in mouse thermogenesis, the frequency may have increased as an adaptation to cold in Finns. Our selection analysis provided significant evidence of selection for the abdominal obesity risk allele T of rs10494217 in Finns, with a north-to-south decreasing trend in other populations, and demonstrated that latitude significantly predicts the allele frequency. We also discovered that the risk allele status significantly affects SAT adipocyte expression of multiple adipocyte marker genes in trans in two cohorts. Two of these trans genes have been connected to thermogenesis, supporting the thermogenic effect of the TBX15 missense variant as a possible cause of its selection. Adipose expression of one trans gene, a lncRNA, AC002066.1, was strongly associated with adipocyte size, implicating it in metabolically unhealthy adipocyte hypertrophy. In summary, the abdominal obesity variant rs10494217 was selected in Finns, and individuals with the risk allele have trans effects on adipocyte expression of genes relating to thermogenesis and adipocyte hypertrophy.

Research on brain expression quantitative trait loci (eQTLs) has illuminated the genetic underpinnings of schizophrenia (SCZ). Yet most of these studies have been centered on European populations, leading to a constrained understanding of population diversities and disease risks. To address this gap, we examined genotype and RNA-seq data from African Americans (AA, n = 158), Europeans (EUR, n = 408), and East Asians (EAS, n = 217). When comparing eQTLs between EUR and non-EUR populations, we observed concordant patterns of genetic regulatory effect, particularly in terms of the effect sizes of the eQTLs. However, 343,737 cis-eQTLs linked to 1,276 genes and 198,769 SNPs were found to be specific to non-EUR populations. Over 90% of observed population differences in eQTLs could be traced back to differences in allele frequency. Furthermore, 35% of these eQTLs were notably rare in the EUR population. Integrating brain eQTLs with SCZ signals from diverse populations, we observed a higher disease heritability enrichment of brain eQTLs in matched populations compared to mismatched ones. Prioritization analysis identified five risk genes (SFXN2, VPS37B, DENR, FTCDNL1, and NT5DC2) and three potential regulatory variants in known risk genes (CNNM2, MTRFR, and MPHOSPH9) that were missed in the EUR dataset. Our findings underscore that increasing genetic ancestral diversity is more efficient for power improvement than merely increasing the sample size within single-ancestry eQTLs datasets. Such a strategy will not only improve our understanding of the biological underpinnings of population structures but also pave the way for the identification of risk genes in SCZ.