Search

PURPOSE: To describe the clinical, electrophysiological and genetic spectrum of inherited retinal diseases associated with variants in the PRPH2 gene. METHODS: A total of 241 patients from 168 families across 15 sites in 9 countries with pathogenic or likely pathogenic variants in PRPH2 were included. Records were reviewed for age at symptom onset, visual acuity, full-field ERG, fundus colour photography, fundus autofluorescence (FAF), and SD-OCT. Images were graded into six phenotypes. Statistical analyses were performed to determine genotype-phenotype correlations. RESULTS: The median age at symptom onset was 40 years (range, 4-78 years). FAF phenotypes included normal (5%), butterfly pattern dystrophy, or vitelliform macular dystrophy (11%), central areolar choroidal dystrophy (28%), pseudo-Stargardt pattern dystrophy (41%), and retinitis pigmentosa (25%). Symptom onset was earlier in retinitis pigmentosa as compared with pseudo-Stargardt pattern dystrophy (34 vs 44 years; P = 0.004). The median visual acuity was 0.18 logMAR (interquartile range, 0-0.54 logMAR) and 0.18 logMAR (interquartile range 0-0.42 logMAR) in the right and left eyes, respectively. ERG showed a significantly reduced amplitude across all components (P < 0.001) and a peak time delay in the light-adapted 30-Hz flicker and single-flash b-wave (P < 0.001). Twenty-two variants were novel. The central areolar choroidal dystrophy phenotype was associated with 13 missense variants. The remaining variants showed marked phenotypic variability. CONCLUSIONS: We described six distinct FAF phenotypes associated with variants in the PRPH2 gene. One FAF phenotype may have multiple ERG phenotypes, demonstrating a discordance between structure and function. Given the vast spectrum of PRPH2 disease our findings are useful for future clinical trials.

The Phenome-Wide Association Study (PheWAS) is increasingly used to broadly screen for potential treatment effects, e.g., IL6R variant as a proxy for IL6R antagonists. This approach offers an opportunity to address the limited power in clinical trials to study differential treatment effects across patient subgroups. However, limited methods exist to efficiently test for differences across subgroups in the thousands of multiple comparisons generated as part of a PheWAS. In this study, we developed an approach that maximizes the power to test for heterogeneous genotype-phenotype associations and applied this approach to an IL6R PheWAS among individuals of African (AFR) and European (EUR) ancestries. We identified 29 traits with differences in IL6R variant-phenotype associations, including a lower risk of type 2 diabetes in AFR (OR 0.96) vs EUR (OR 1.0, p-value for heterogeneity = 8.5 × 10-3), and higher white blood cell count (p-value for heterogeneity = 8.5 × 10-131). These data suggest a more salutary effect of IL6R blockade for T2D among individuals of AFR vs EUR ancestry and provide data to inform ongoing clinical trials targeting IL6 for an expanding number of conditions. Moreover, the method to test for heterogeneity of associations can be applied broadly to other large-scale genotype-phenotype screens in diverse populations.

Genetic variants in the SLC6A1 gene can cause a broad phenotypic disease spectrum by altering the protein function. Thus, systematically curated clinically relevant genotype-phenotype associations are needed to understand the disease mechanism and improve therapeutic decision-making. We aggregated genetic and clinical data from 172 individuals with likely pathogenic/pathogenic (lp/p) SLC6A1 variants and functional data for 184 variants (14.1% lp/p). Clinical and functional data were available for a subset of 126 individuals. We explored the potential associations of variant positions on the GAT1 3D structure with variant pathogenicity, altered molecular function and phenotype severity using bioinformatic approaches. The GAT1 transmembrane domains 1, 6 and extracellular loop 4 (EL4) were enriched for patient over population variants. Across functionally tested missense variants (n = 156), the spatial proximity from the ligand was associated with loss-of-function in the GAT1 transporter activity. For variants with complete loss of in vitro GABA uptake, we found a 4.6-fold enrichment in patients having severe disease versus non-severe disease (P = 2.9 × 10-3, 95% confidence interval: 1.5-15.3). In summary, we delineated associations between the 3D structure and variant pathogenicity, variant function and phenotype in SLC6A1-related disorders. This knowledge supports biology-informed variant interpretation and research on GAT1 function. All our data can be interactively explored in the SLC6A1 portal (https://slc6a1-portal.broadinstitute.org/).

BRAT1 biallelic variants are associated with rigidity and multifocal seizure syndrome, lethal neonatal (RMFSL), and neurodevelopmental disorder associating cerebellar atrophy with or without seizures syndrome (NEDCAS). To date, forty individuals have been reported in the literature. We collected clinical and molecular data from 57 additional cases allowing us to study a large cohort of 97 individuals and draw phenotype-genotype correlations. Fifty-nine individuals presented with BRAT1-related RMFSL phenotype. Most of them had no psychomotor acquisition (100%), epilepsy (100%), microcephaly (91%), limb rigidity (93%), and died prematurely (93%). Thirty-eight individuals presented a non-lethal phenotype of BRAT1-related NEDCAS phenotype. Seventy-six percent of the patients in this group were able to walk and 68% were able to say at least a few words. Most of them had cerebellar ataxia (82%), axial hypotonia (79%) and cerebellar atrophy (100%). Genotype-phenotype correlations in our cohort revealed that biallelic nonsense, frameshift or inframe deletion/insertion variants result in the severe BRAT1-related RMFSL phenotype (46/46; 100%). In contrast, genotypes with at least one missense were more likely associated with NEDCAS (28/34; 82%). The phenotype of patients carrying splice variants was variable: 41% presented with RMFSL (7/17) and 59% with NEDCAS (10/17).

Cornelia de Lange Syndrome (CdLS) is a rare, dominantly inherited multisystem developmental disorder characterized by highly variable manifestations of growth and developmental delays, upper limb involvement, hypertrichosis, cardiac, gastrointestinal, craniofacial, and other systemic features. Pathogenic variants in genes encoding cohesin complex structural subunits and regulatory proteins (NIPBL, SMC1A, SMC3, HDAC8, and RAD21) are the major pathogenic contributors to CdLS. Heterozygous or hemizygous variants in the genes encoding these five proteins have been found to be contributory to CdLS, with variants in NIPBL accounting for the majority (>60%) of cases, and the only gene identified to date that results in the severe or classic form of CdLS when mutated. Pathogenic variants in cohesin genes other than NIPBL tend to result in a less severe phenotype. Causative variants in additional genes, such as ANKRD11, EP300, AFF4, TAF1, and BRD4, can cause a CdLS-like phenotype. The common role that these genes, and others, play as critical regulators of developmental transcriptional control has led to the conditions they cause being referred to as disorders of transcriptional regulation (or DTRs). Here, we report the results of a comprehensive molecular analysis in a cohort of 716 probands with typical and atypical CdLS in order to delineate the genetic contribution of causative variants in cohesin complex genes as well as novel candidate genes, genotype-phenotype correlations, and the utility of genome sequencing in understanding the mutational landscape in this population.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a rare, irreversible, and progressive disease of the lungs. Common genetic variants, in addition to nongenetic factors, have been consistently associated with IPF. Rare variants identified by candidate gene, family-based, and exome studies have also been reported to associate with IPF. However, the extent to which rare variants, genome-wide, may contribute to the risk of IPF remains unknown. Objectives: We used whole-genome sequencing to investigate the role of rare variants, genome-wide, on IPF risk. Methods: As part of the Trans-Omics for Precision Medicine Program, we sequenced 2,180 cases of IPF. Association testing focused on the aggregated effect of rare variants (minor allele frequency ⩽0.01) within genes or regions. We also identified individual rare variants that are influential within genes and estimated the heritability of IPF on the basis of rare and common variants. Measurements and Main Results: Rare variants in both TERT and RTEL1 were significantly associated with IPF. A single rare variant in each of the TERT and RTEL1 genes was found to consistently influence the aggregated test statistics. There was no significant evidence of association with other previously reported rare variants. The SNP heritability of IPF was estimated to be 32% (SE = 3%). Conclusions: Rare variants within the TERT and RTEL1 genes and well-established common variants have the largest contribution to IPF risk overall. Efforts in risk profiling or the development of therapies for IPF that focus on TERT, RTEL1, common variants, and environmental risk factors are likely to have the largest impact on this complex disease.

Genetic association studies of many heritable traits resulting from physiological testing often have modest sample sizes due to the cost and burden of the required phenotyping. This reduces statistical power and limits discovery of multiple genetic associations. We present a strategy to leverage pleiotropy between traits to both discover new loci and to provide mechanistic hypotheses of the underlying pathophysiology. Specifically, we combine a colocalization test with a locus-level test of pleiotropy. In simulations, we show that this approach is highly selective for identifying true pleiotropy driven by the same causative variant, thereby improves the chance to replicate the associations in underpowered validation cohorts and leads to higher interpretability. Here, as an exemplar, we use Obstructive Sleep Apnea (OSA), a common disorder diagnosed using overnight multi-channel physiological testing. We leverage pleiotropy with relevant cellular and cardio-metabolic phenotypes and gene expression traits to map new risk loci in an underpowered OSA GWAS. We identify several pleiotropic loci harboring suggestive associations to OSA and genome-wide significant associations to other traits, and show that their OSA association replicates in independent cohorts of diverse ancestries. By investigating pleiotropic loci, our strategy allows proposing new hypotheses about OSA pathobiology across many physiological layers. For example, we identify and replicate the pleiotropy across the plateletcrit, OSA and an eQTL of DNA primase subunit 1 (PRIM1) in immune cells. We find suggestive links between OSA, a measure of lung function (FEV1/FVC), and an eQTL of matrix metallopeptidase 15 (MMP15) in lung tissue. We also link a previously known genome-wide significant peak for OSA in the hexokinase 1 (HK1) locus to hematocrit and other red blood cell related traits. Thus, the analysis of pleiotropic associations has the potential to assemble diverse phenotypes into a chain of mechanistic hypotheses that provide insight into the pathogenesis of complex human diseases.

BackgroundA key goal of precision medicine is to disaggregate common, complex diseases into discrete molecular subtypes. Rare coding variants in the low-density lipoprotein receptor gene (LDLR) are identified in 1% to 2% of coronary artery disease (CAD) patients, defining a molecular subtype with risk driven by hypercholesterolemia.MethodsTo search for additional subtypes, we compared the frequency of rare, predicted loss-of-function and damaging missense variants aggregated within a given gene in 41 081 CAD cases versus 217 115 controls.ResultsRare variants in LDLR were most strongly associated with CAD, present in 1% of cases and associated with 4.4-fold increased CAD risk. A second subtype was characterized by variants in endothelial nitric oxide synthase gene (NOS3), a key enzyme regulating vascular tone, endothelial function, and platelet aggregation. A rare predicted loss-of-function or damaging missense variants in NOS3 was present in 0.6% of cases and associated with 2.42-fold increased risk of CAD (95% CI, 1.80-3.26; P=5.50×10-9). These variants were associated with higher systolic blood pressure (+3.25 mm Hg; [95% CI, 1.86-4.65]; P=5.00×10-6) and increased risk of hypertension (adjusted odds ratio 1.31; [95% CI, 1.14-1.51]; P=2.00×10-4) but not circulating cholesterol concentrations, suggesting that, beyond lipid pathways, nitric oxide synthesis is a key nonlipid driver of CAD risk.ConclusionsBeyond LDLR, we identified an additional nonlipid molecular subtype of CAD characterized by rare variants in the NOS3 gene.

BackgroundInherited retinal dystrophies describe a heterogeneous group of retinal diseases that lead to the irreversible degeneration of rod and cone photoreceptors and eventual blindness. Recessive loss-of-function mutations in Tubulin Tyrosine Ligase Like 5 (TTLL5) represent a recently described cause of inherited cone-rod and cone dystrophy. This study describes the unusual phenotypes of three patients with autosomal recessive mutations in TTLL5. Examination of these patients included funduscopic evaluation, spectral-domain optical coherence tomography, short-wavelength autofluorescence, and full-field electroretinography (ffERG). Genetic diagnoses were confirmed using whole exome capture. Protein modeling of the identified variants was performed to explore potential genotype-phenotype correlations.ResultsGenetic testing revealed five novel variants in TTLL5 in three unrelated patients with retinal dystrophy. Clinical imaging demonstrated features of sectoral cone-rod dystrophy and cone dystrophy, with phenotypic variability seen across all three patients. One patient also developed high-frequency hearing loss during a similar time period as the onset of retinal disease, potentially suggestive of a syndromic disorder. Retinal structure findings were corroborated with functional measures including ffERG findings that supported these diagnoses. Modeling of the five variants suggest that they cause different effects on protein function, providing a potential reason for genotype-phenotype correlation in these patients.ConclusionsThe authors report retinal phenotypic findings in three unrelated patients with novel mutations causing autosomal recessive TTLL5-mediated retinal dystrophy. These findings broaden the understanding of the phenotypes associated with TTLL5-mediated retinal disease and suggest that mutations in TTLL5 should be considered as a potential cause of sectoral retinal dystrophy in addition to cone-rod and cone dystrophies.

Biobanks facilitate genome-wide association studies (GWASs), which have mapped genomic loci across a range of human diseases and traits. However, most biobanks are primarily composed of individuals of European ancestry. We introduce the Global Biobank Meta-analysis Initiative (GBMI)-a collaborative network of 23 biobanks from 4 continents representing more than 2.2 million consented individuals with genetic data linked to electronic health records. GBMI meta-analyzes summary statistics from GWASs generated using harmonized genotypes and phenotypes from member biobanks for 14 exemplar diseases and endpoints. This strategy validates that GWASs conducted in diverse biobanks can be integrated despite heterogeneity in case definitions, recruitment strategies, and baseline characteristics. This collaborative effort improves GWAS power for diseases, benefits understudied diseases, and improves risk prediction while also enabling the nomination of disease genes and drug candidates by incorporating gene and protein expression data and providing insight into the underlying biology of human diseases and traits.

We assessed genotype-phenotype correlations among the visual, auditory, and olfactory phenotypes of 127 participants with Usher syndrome (USH2) (n =80) or nonsyndromic autosomal recessive retinitis pigmentosa (ARRP) (n = 47) due to USH2A variants, using clinical data and molecular diagnostics from the Rate of Progression in USH2A Related Retinal Degeneration (RUSH2A) study. USH2A truncating alleles were associated with USH2 and had a dose-dependent effect on hearing loss severity with no effect on visual loss severity within the USH2 subgroup. A group of missense alleles in an interfibronectin domain appeared to be hypomorphic in ARRP. These alleles were associated with later age of onset, larger visual field area, better sensitivity thresholds, and better electroretinographic responses. No effect of genotype on the severity of olfactory deficits was observed. This study unveils a unique, tissue-specific USH2A allelic hierarchy with important prognostic implications for patient counseling and treatment trial endpoints. These findings may inform clinical care or research approaches in others with allelic disorders or pleiotropic phenotypes.

Background and objectivesTo study human leukocyte antigen (HLA) allele associations in anti-leucine-rich glioma-inactivated 1 (LGI1) encephalitis.MethodsA multiethnic cohort of 269 patients with anti-LGI1 encephalitis and 1,359 controls was included. Four-digit HLA sequencing and genome wide association single-nucleotide polymorphism typing imputation (0.99 concordance) were used for HLA typing. Significance of primary and secondary associations was tested using χ2, Fisher exact tests, or logistic regression with the control of population stratification covariates when applicable.ResultsDRB1*07:01 and DQA1*02:01, 2 alleles in strong linkage disequilibrium, were associated with the disease (90% vs 24%, OR = 27.8, p < 10e-50) across ethnicity independent of variation at DRB3 and DQB1, 2 flanking HLA loci. DRB1*07:01 homozygosity was associated with a doubling of risk (OR = 2.1, p = 0.010), suggesting causality. DRB1*07:01 negative subjects were younger (p = 0.003) and more frequently female (p = 0.015). Three patients with malignant thymomas did not carry DRB1*07:01, whereas patients with other tumors had high DRB1*07:01 frequency, suggesting that the presence of tumors other than thymomas may be coincidental and not causal. In both DRB1*07:01 heterozygous individuals and DRB1*07:01 negative subjects, DRB1*04:02 was associated with anti-LGI1 encephalitis, indicating an independent effect of this allele (OR = 6.85, p = 4.57 × 10-6 and OR = 8.93, p = 2.50 × 10-3, respectively). DRB1*04:02 was also independently associated with younger age at onset (β = -6.68, p = 9.78 × 10-3). Major histocompatibility complex peptide-binding predictions using LGI1-derived peptides revealed divergent binding propensities for DRB1*04:02 and DRB1*07:01 alleles, suggesting independent pathogenic mechanisms.DiscussionIn addition to the established primary DRB1*07:01 association in anti-LGI1 encephalitis, we observe a secondary effect of DRB1*04:02 with lower age at onset. Our study provides evidence for secondary effects within HLA locus that correlate with clinical phenotypes in anti-LGI1 encephalitis.

To determine the impact of genetic variants on the brain, we used genetically informed brain atlases in genome-wide association studies of regional cortical surface area and thickness in 39,898 adults and 9136 children. We uncovered 440 genome-wide significant loci in the discovery cohort and 800 from a post hoc combined meta-analysis. Loci in adulthood were largely captured in childhood, showing signatures of negative selection, and were linked to early neurodevelopment and pathways associated with neuropsychiatric risk. Opposing gradations of decreased surface area and increased thickness were associated with common inversion polymorphisms. Inferior frontal regions, encompassing Broca's area, which is important for speech, were enriched for human-specific genomic elements. Thus, a mixed genetic landscape of conserved and human-specific features is concordant with brain hierarchy and morphogenetic gradients.

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia and is characterized by chromosomal aberrations including 13q, 11q, and 17p deletions and a trisomy of chromosome 12 (T12). 13q deletions are often associated with 11q and 17p deletions in aggressive cases. Conversely, T12 CLLs show a variable prognosis, and association with 13q deletions is uncommon. The miR-15a/16-1 cluster is the functional target of 13q deletions, leading to BCL2 overexpression. Chromosomal aberrations in CLL are associated with prognosis, and their identification is carried out by fluorescence in situ hybridization (FISH). Since standard FISH only detects large deletions, we investigated the presence of undetected microdeletions targeting miR-15a/16-1 in CLL cases. We found that ∼34% of CLL samples show an unreported loss of the miR-15a/16-1 locus regardless of their cytogenetic profile. Interestingly, 15 out of 39 (∼39%) of all CLLs with T12, carry microdeletions of miR-15a/16-1, indicating that, in patients with T12, miR-15a/16-1 are mostly inactivated by microdeletions. In addition, ∼40% of CLL cases bearing T12, 17p-, and 11q- showed unidentified microdeletions of miR-15a/16-1, suggesting that miR-15a/16-1 loss cooperates with such chromosomal alterations in CLL. These data may have clinical relevance for the successful stratification of patients for treatment.

As the master regulator in utero, the placenta is core to the Developmental Origins of Health and Disease (DOHaD) hypothesis but is historically understudied. To identify placental gene-trait associations (GTAs) across the life course, we perform distal mediator-enriched transcriptome-wide association studies (TWAS) for 40 traits, integrating placental multi-omics from the Extremely Low Gestational Age Newborn Study. At [Formula: see text], we detect 248 GTAs, mostly for neonatal and metabolic traits, across 176 genes, enriched for cell growth and immunological pathways. In aggregate, genetic effects mediated by placental expression significantly explain 4 early-life traits but no later-in-life traits. 89 GTAs show significant mediation through distal genetic variants, identifying hypotheses for distal regulation of GTAs. Investigation of one hypothesis in human placenta-derived choriocarcinoma cells reveal that knockdown of mediator gene EPS15 upregulates predicted targets SPATA13 and FAM214A, both associated with waist-hip ratio in TWAS, and multiple genes involved in metabolic pathways. These results suggest profound health impacts of placental genomic regulation in developmental programming across the life course.

Although the cohort-level accuracy of polygenic risk scores (PRSs)-estimates of genetic value at the individual level-has been widely assessed, uncertainty in PRSs remains underexplored. In the present study, we show that Bayesian PRS methods can estimate the variance of an individual's PRS and can yield well-calibrated credible intervals via posterior sampling. For 13 real traits in the UK Biobank (n = 291,273 unrelated 'white British'), we observe large variances in individual PRS estimates which impact interpretation of PRS-based stratification; averaging across traits, only 0.8% (s.d. = 1.6%) of individuals with PRS point estimates in the top decile have corresponding 95% credible intervals fully contained in the top decile. We provide an analytical estimator for the expectation of individual PRS variance as a function of SNP heritability, number of causal SNPs and sample size. Our results showcase the importance of incorporating uncertainty in individual PRS estimates into subsequent analyses.

Dilated cardiomyopathy (DCM) is a life-threatening form of heart disease that is typically characterized by progressive thinning of the ventricular walls, chamber dilation, and systolic dysfunction. Multiple mutations in the gene encoding filamin C (FLNC), an actin-binding cytoskeletal protein in cardiomyocytes, have been found in patients with DCM. However, the mechanisms that lead to contractile impairment and DCM in patients with FLNC variants are poorly understood. To determine how FLNC regulates systolic force transmission and DCM remodeling, we used an inducible, cardiac-specific FLNC-knockout (icKO) model to produce a rapid onset of DCM in adult mice. Loss of FLNC reduced systolic force development in single cardiomyocytes and isolated papillary muscles but did not affect twitch kinetics or calcium transients. Electron and immunofluorescence microscopy showed significant defects in Z-disk alignment in icKO mice and altered myofilament lattice geometry. Moreover, a loss of FLNC induces a softening myocyte cortex and structural adaptations at the subcellular level that contribute to disrupted longitudinal force production during contraction. Spatially explicit computational models showed that these structural defects could be explained by a loss of inter-myofibril elastic coupling at the Z-disk. Our work identifies FLNC as a key regulator of the multiscale ultrastructure of cardiomyocytes and therefore plays an important role in maintaining systolic mechanotransmission pathways, the dysfunction of which may be key in driving progressive DCM.

Genetically regulated gene expression has helped elucidate the biological mechanisms underlying complex traits. Improved high-throughput technology allows similar interrogation of the genetically regulated proteome for understanding complex trait mechanisms. Here, we used the Trans-omics for Precision Medicine (TOPMed) Multi-omics pilot study, which comprises data from Multi-Ethnic Study of Atherosclerosis (MESA), to optimize genetic predictors of the plasma proteome for genetically regulated proteome-wide association studies (PWAS) in diverse populations. We built predictive models for protein abundances using data collected in TOPMed MESA, for which we have measured 1,305 proteins by a SOMAscan assay. We compared predictive models built via elastic net regression to models integrating posterior inclusion probabilities estimated by fine-mapping SNPs prior to elastic net. In order to investigate the transferability of predictive models across ancestries, we built protein prediction models in all four of the TOPMed MESA populations, African American (n = 183), Chinese (n = 71), European (n = 416), and Hispanic/Latino (n = 301), as well as in all populations combined. As expected, fine-mapping produced more significant protein prediction models, especially in African ancestries populations, potentially increasing opportunity for discovery. When we tested our TOPMed MESA models in the independent European INTERVAL study, fine-mapping improved cross-ancestries prediction for some proteins. Using GWAS summary statistics from the Population Architecture using Genomics and Epidemiology (PAGE) study, which comprises ∼50,000 Hispanic/Latinos, African Americans, Asians, Native Hawaiians, and Native Americans, we applied S-PrediXcan to perform PWAS for 28 complex traits. The most protein-trait associations were discovered, colocalized, and replicated in large independent GWAS using proteome prediction model training populations with similar ancestries to PAGE. At current training population sample sizes, performance between baseline and fine-mapped protein prediction models in PWAS was similar, highlighting the utility of elastic net. Our predictive models in diverse populations are publicly available for use in proteome mapping methods at https://doi.org/10.5281/zenodo.4837327.

BackgroundPTEN is a well-established risk gene for autism spectrum disorder (ASD). Yet, little is known about how PTEN mutations and associated molecular processes influence neurobehavioral function in mutation carriers with (PTEN-ASD) and without ASD (PTEN no-ASD). The primary aim of the present study was to examine group differences in peripheral blood-derived PTEN pathway protein levels between PTEN-ASD, PTEN no-ASD, and idiopathic macrocephalic ASD patients (macro-ASD). Secondarily, associations between protein levels and neurobehavioral functions were examined in the full cohort.MethodsPatients were recruited at four tertiary medical centers. Peripheral blood-derived protein levels from canonical PTEN pathways (PI3K/AKT and MAPK/ERK) were analyzed using Western blot analyses blinded to genotype and ASD status. Neurobehavioral measures included standardized assessments of global cognitive ability and multiple neurobehavioral domains. Analysis of variance models examined group differences in demographic, neurobehavioral, and protein measures. Bivariate correlations, structural models, and statistical learning procedures estimated associations between molecular and neurobehavioral variables. To complement patient data, Western blots for downstream proteins were generated to evaluate canonical PTEN pathways in the PTEN-m3m4 mouse model.ResultsParticipants included 61 patients (25 PTEN-ASD, 16 PTEN no-ASD, and 20 macro-ASD). Decreased PTEN and S6 were observed in both PTEN mutation groups. Reductions in MnSOD and increases in P-S6 were observed in ASD groups. Elevated neural P-AKT/AKT and P-S6/S6 from PTEN murine models parallel our patient observations. Patient PTEN and AKT levels were independently associated with global cognitive ability, and p27 expression was associated with frontal sub-cortical functions. As a group, molecular measures added significant predictive value to several neurobehavioral domains over and above PTEN mutation status.LimitationsSample sizes were small, precluding within-group analyses. Protein and neurobehavioral data were limited to a single evaluation. A small number of patients were excluded with invalid protein data, and cognitively impaired patients had missing data on some assessments.ConclusionsSeveral canonical PTEN pathway molecules appear to influence the presence of ASD and modify neurobehavioral function in PTEN mutation patients. Protein assays of the PTEN pathway may be useful for predicting neurobehavioral outcomes in PTEN patients. Future longitudinal analyses are needed to replicate these findings and evaluate within-group relationships between protein and neurobehavioral measures.Trial registrationClinicalTrials.gov Identifier NCT02461446.

BackgroundPost-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells.ResultsHere, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees.ConclusionsThis study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

BackgroundWe aimed to define the clinical and variant spectrum and to provide novel molecular insights into the DHX30-associated neurodevelopmental disorder.MethodsClinical and genetic data from affected individuals were collected through Facebook-based family support group, GeneMatcher, and our network of collaborators. We investigated the impact of novel missense variants with respect to ATPase and helicase activity, stress granule (SG) formation, global translation, and their effect on embryonic development in zebrafish. SG formation was additionally analyzed in CRISPR/Cas9-mediated DHX30-deficient HEK293T and zebrafish models, along with in vivo behavioral assays.ResultsWe identified 25 previously unreported individuals, ten of whom carry novel variants, two of which are recurrent, and provide evidence of gonadal mosaicism in one family. All 19 individuals harboring heterozygous missense variants within helicase core motifs (HCMs) have global developmental delay, intellectual disability, severe speech impairment, and gait abnormalities. These variants impair the ATPase and helicase activity of DHX30, trigger SG formation, interfere with global translation, and cause developmental defects in a zebrafish model. Notably, 4 individuals harboring heterozygous variants resulting either in haploinsufficiency or truncated proteins presented with a milder clinical course, similar to an individual harboring a de novo mosaic HCM missense variant. Functionally, we established DHX30 as an ATP-dependent RNA helicase and as an evolutionary conserved factor in SG assembly. Based on the clinical course, the variant location, and type we establish two distinct clinical subtypes. DHX30 loss-of-function variants cause a milder phenotype whereas a severe phenotype is caused by HCM missense variants that, in addition to the loss of ATPase and helicase activity, lead to a detrimental gain-of-function with respect to SG formation. Behavioral characterization of dhx30-deficient zebrafish revealed altered sleep-wake activity and social interaction, partially resembling the human phenotype.ConclusionsOur study highlights the usefulness of social media to define novel Mendelian disorders and exemplifies how functional analyses accompanied by clinical and genetic findings can define clinically distinct subtypes for ultra-rare disorders. Such approaches require close interdisciplinary collaboration between families/legal representatives of the affected individuals, clinicians, molecular genetics diagnostic laboratories, and research laboratories.

BackgroundAngelman Syndrome (AS) is a rare genetic disorder characterized by impaired communication, motor and balance deficits, intellectual disabilities, recurring seizures and abnormal sleep patterns. The genetic cause of AS is neuronal-specific loss of expression of UBE3A (ubiquitin-protein ligase E6-AP), an imprinted gene. Seizure and sleep disorders are highly prevalent (> 80%) in the AS population. The present experiments were designed to identify translational, neurophysiological outcome measures in a model of AS.MethodsWe used the exon-2 deletion mouse (Ube3a-del) on a C57BL/6J background to assess seizure, sleep and electrophysiological phenotypes. Seizure susceptibility has been reported in Ube3a-del mice with a variety of seizure induction methods. Here, we provoked seizures by a single high-dose injection of 80 mg/kg pentylenetetrazole. Novel experiments included the utilization of wireless telemetry devices to acquire global electroencephalogram (EEG) and neurophysiological data on electrographic seizures, power spectra, light-dark cycles, sleep stages and sleep spindles in Ube3a-del and WT mice.ResultsUbe3a-del mice exhibited reduced seizure threshold compared to WT. EEG illustrated that Ube3a-del mice had increased epileptiform spiking activity and delta power, which corroborates findings from other laboratories and recapitulates clinical reports in AS. This is the first report to use a cortical surface-based recording by a wireless telemetry device over tethered/fixed head-mount depth recordings. Less time in both paradoxical and slow-wave sleep, longer latencies to paradoxical sleep stages and total less sleep time in Ube3a-del mice were observed compared to WT. For the first time, we detected fewer sleep spindles in the AS mouse model.LimitationsThis study was limited to the exon 2 deletion mouse model, and future work will investigate the rat model of AS, containing a complete Ube3a deletion and pair EEG with behavior.ConclusionsOur data enhance rigor and translatability as our study provides important corroboration of previous reports on epileptiform and elevated delta power. For the first time we report neurophysiological phenotypes collected via translational methodology. Furthermore, this is the first report of reduced sleep spindles, a critical marker of memory consolidation during sleep, in an AS model. Our results are useful outcomes for therapeutic testing.

BackgroundSleep-disordered breathing is a common disorder associated with significant morbidity. The genetic architecture of sleep-disordered breathing remains poorly understood. Through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, we performed the first whole-genome sequence analysis of sleep-disordered breathing.MethodsThe study sample was comprised of 7988 individuals of diverse ancestry. Common-variant and pathway analyses included an additional 13,257 individuals. We examined five complementary traits describing different aspects of sleep-disordered breathing: the apnea-hypopnea index, average oxyhemoglobin desaturation per event, average and minimum oxyhemoglobin saturation across the sleep episode, and the percentage of sleep with oxyhemoglobin saturation

Many computational pipelines exist for the detection of differentially expressed genes. However, computational pipelines for functional gene detection rarely exist. We developed a new computational pipeline for functional gene identification from transcriptome profiling data. Key features of the pipeline include batch effect correction, clustering optimization by gap statistics, gene ontology analysis of clustered genes, and literature analysis for functional gene discovery. By leveraging this pipeline on RNA-seq datasets from two mouse retinal development studies, we identified 7 candidate genes involved in the formation of the photoreceptor outer segment. The expression of top three candidate genes (Pde8b, Laptm4b, and Nr1h4) in the outer segment of the developing mouse retina were experimentally validated by immunohistochemical analysis. This computational pipeline can accurately predict novel functional gene for a specific biological process, e.g., development of the outer segment and synapses of the photoreceptor cells in the mouse retina. This pipeline can also be useful to discover functional genes for other biological processes and in other organs and tissues.

A frameshift deletion variant in the Wnt pathway gene dishevelled 2 (DVL2) is associated with a truncated, kinked tail ("screw tail") in English Bulldogs, French Bulldogs and Boston Terriers. These breeds are also characterized by distinctive morphological traits, including a wide head, flat face and short-limbed dwarfism, which are characteristic of Robinow syndrome in humans, caused by defects in genes such as DVL1 and DVL3. Based on these phenotypic and genetic similarities, it has previously been hypothesized that the canine DVL2 variant results in a syndromic phenotype called the Robinow-like syndrome. In our study, we investigated the distribution of the DVL2 variant in 1954 dogs from 15 breeds, identifying breeds with allele variation and enabling the dissection of the genotype-phenotype correlation for the first time. With CT examinations in American Staffordshire Terriers, we confirmed that the DVL2 allele is associated with caudal vertebral malformations and a brachycephalic phenotype. We also hypothesize that the variant may be linked to additional health conditions, including brachycephalic obstructive airway syndrome and congenital heart defects. Altogether, our study strengthens the role of DVL2 as one of the contributors to the "bulldog type" morphology and features on the spectrum of human Robinow syndrome.

We have previously suggested a central role for mitochondria in the observed sex differences in metabolic traits. However, the mechanisms by which sex differences affect adipose mitochondrial function and metabolic syndrome are unclear. Here we show that in both mice and humans, adipose mitochondrial functions are elevated in females and are strongly associated with adiposity, insulin resistance and plasma lipids. Using a panel of diverse inbred strains of mice, we identify a genetic locus on mouse chromosome 17 that controls mitochondrial mass and function in adipose tissue in a sex- and tissue-specific manner. This locus contains Ndufv2 and regulates the expression of at least 89 mitochondrial genes in females, including oxidative phosphorylation genes and those related to mitochondrial DNA content. Overexpression studies indicate that Ndufv2 mediates these effects by regulating supercomplex assembly and elevating mitochondrial reactive oxygen species production, which generates a signal that increases mitochondrial biogenesis.

Genotype imputation using a reference panel that combines high-density array data and publicly available whole genome sequence consortium variant data is potentially a cost-effective method to increase the density of extant lower-density array datasets. In this study, three datasets (two Border Collie; one Italian Spinone) generated using a legacy array (Illumina CanineHD, 173 662 SNPs) were utilised to assess the feasibility and accuracy of this approach and to gather additional evidence for the efficacy of canine genotype imputation. The cosmopolitan reference panels used to impute genotypes comprised dogs of 158 breeds, mixed breed dogs, wolves and Chinese indigenous dogs, as well as breed-specific individuals genotyped using the Axiom Canine HD array. The two Border Collie reference panels comprised 808 individuals including 79 Border Collies and 426 326 or 426 332 SNPs; and the Italian Spinone reference panel comprised 807 individuals including 38 Italian Spinoni and 476 313 SNPs. A high accuracy for imputation was observed, with the lowest accuracy observed for one of the Border Collie datasets (mean R2  = 0.94) and the highest for the Italian Spinone dataset (mean R2  = 0.97). This study's findings demonstrate that imputation of a legacy array study set using a reference panel comprising both breed-specific array data and multi-breed variant data derived from whole genomes is effective and accurate. The process of canine genotype imputation, using the valuable growing resource of publicly available canine genome variant datasets alongside breed-specific data, is described in detail to facilitate and encourage use of this technique in canine genetics.

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase statistical power. Harmonization of the data usually requires substantial effort due to heterogeneity in phenotype definitions, study design, data collection procedures, and data-set organization. Here we describe a centralized system for phenotype harmonization that includes input from phenotype domain and study experts, quality control, documentation, reproducible results, and data-sharing mechanisms. This system was developed for the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program, which is generating genomic and other -omics data for more than 80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants (recruited in 1948-2012) from up to 17 studies per phenotype. Here we discuss challenges in this undertaking and how they were addressed. The harmonized phenotype data and associated documentation have been submitted to National Institutes of Health data repositories for controlled access by the scientific community. We also provide materials to facilitate future harmonization efforts by the community, which include 1) the software code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify, or extend these harmonizations to additional studies, and 2) the results of labeling thousands of phenotype variables with controlled vocabulary terms.

The phenotypic variability associated with pathogenic variants in Lysine Acetyltransferase 6B (KAT6B, a.k.a. MORF, MYST4) results in several interrelated syndromes including Say-Barber-Biesecker-Young-Simpson Syndrome and Genitopatellar Syndrome. Here we present 20 new cases representing 10 novel KAT6B variants. These patients exhibit a range of clinical phenotypes including intellectual disability, mobility and language difficulties, craniofacial dysmorphology, and skeletal anomalies. Given the range of features previously described for KAT6B-related syndromes, we have identified additional phenotypes including concern for keratoconus, sensitivity to light or noise, recurring infections, and fractures in greater numbers than previously reported. We surveyed clinicians to qualitatively assess the ways families engage with genetic counselors upon diagnosis. We found that 56% (10/18) of individuals receive diagnoses before the age of 2 years (median age = 1.96 years), making it challenging to address future complications with limited accessible information and vast phenotypic severity. We used CRISPR to introduce truncating variants into the KAT6B gene in model cell lines and performed chromatin accessibility and transcriptome sequencing to identify key dysregulated pathways. This study expands the clinical spectrum and addresses the challenges to management and genetic counseling for patients with KAT6B-related disorders.

Behaviors and disorders related to self-regulation, such as substance use, antisocial behavior and attention-deficit/hyperactivity disorder, are collectively referred to as externalizing and have shared genetic liability. We applied a multivariate approach that leverages genetic correlations among externalizing traits for genome-wide association analyses. By pooling data from ~1.5 million people, our approach is statistically more powerful than single-trait analyses and identifies more than 500 genetic loci. The loci were enriched for genes expressed in the brain and related to nervous system development. A polygenic score constructed from our results predicts a range of behavioral and medical outcomes that were not part of genome-wide analyses, including traits that until now lacked well-performing polygenic scores, such as opioid use disorder, suicide, HIV infections, criminal convictions and unemployment. Our findings are consistent with the idea that persistent difficulties in self-regulation can be conceptualized as a neurodevelopmental trait with complex and far-reaching social and health correlates.

Behaviors and disorders related to self-regulation, such as substance use, antisocial behavior and attention-deficit/hyperactivity disorder, are collectively referred to as externalizing and have shared genetic liability. We applied a multivariate approach that leverages genetic correlations among externalizing traits for genome-wide association analyses. By pooling data from ~1.5 million people, our approach is statistically more powerful than single-trait analyses and identifies more than 500 genetic loci. The loci were enriched for genes expressed in the brain and related to nervous system development. A polygenic score constructed from our results predicts a range of behavioral and medical outcomes that were not part of genome-wide analyses, including traits that until now lacked well-performing polygenic scores, such as opioid use disorder, suicide, HIV infections, criminal convictions and unemployment. Our findings are consistent with the idea that persistent difficulties in self-regulation can be conceptualized as a neurodevelopmental trait with complex and far-reaching social and health correlates.

Childhood aggressive behavior (AGG) has a substantial heritability of around 50%. Here we present a genome-wide association meta-analysis (GWAMA) of childhood AGG, in which all phenotype measures across childhood ages from multiple assessors were included. We analyzed phenotype assessments for a total of 328 935 observations from 87 485 children aged between 1.5 and 18 years, while accounting for sample overlap. We also meta-analyzed within subsets of the data, i.e., within rater, instrument and age. SNP-heritability for the overall meta-analysis (AGGoverall) was 3.31% (SE = 0.0038). We found no genome-wide significant SNPs for AGGoverall. The gene-based analysis returned three significant genes: ST3GAL3 (P = 1.6E-06), PCDH7 (P = 2.0E-06), and IPO13 (P = 2.5E-06). All three genes have previously been associated with educational traits. Polygenic scores based on our GWAMA significantly predicted aggression in a holdout sample of children (variance explained = 0.44%) and in retrospectively assessed childhood aggression (variance explained = 0.20%). Genetic correlations (rg) among rater-specific assessment of AGG ranged from rg = 0.46 between self- and teacher-assessment to rg = 0.81 between mother- and teacher-assessment. We obtained moderate-to-strong rgs with selected phenotypes from multiple domains, but hardly with any of the classical biomarkers thought to be associated with AGG. Significant genetic correlations were observed with most psychiatric and psychological traits (range [Formula: see text]: 0.19-1.00), except for obsessive-compulsive disorder. Aggression had a negative genetic correlation (rg = ~-0.5) with cognitive traits and age at first birth. Aggression was strongly genetically correlated with smoking phenotypes (range [Formula: see text]: 0.46-0.60). The genetic correlations between aggression and psychiatric disorders were weaker for teacher-reported AGG than for mother- and self-reported AGG. The current GWAMA of childhood aggression provides a powerful tool to interrogate the rater-specific genetic etiology of AGG.

Key messageGenetic modification of spike architecture is essential for improving wheat yield. Newly identified loci for the 'Miracle wheat' phenotype on chromosomes 1AS and 2BS have significant effects on spike traits. The wheat (Triticum ssp.) inflorescence, also known as a spike, forms an unbranched inflorescence in which the inflorescence meristem generates axillary spikelet meristems (SMs) destined to become sessile spikelets. Previously, we identified the putatively causative mutation in the branched headt (bht) gene (TtBH-A1) of tetraploid wheat (T. turgidum convar. compositum (L.f.) Filat.) responsible for the loss of SM identity, converting the non-branching spike to a branched wheat spike. In the current study, we performed whole-genome quantitative trait loci (QTL) analysis using 146 recombinant inbred lines (RILs) derived from a cross between spike-branching wheat ('Miracle wheat') and an elite durum wheat cultivar showing broad phenotypic variation for spike architecture. Besides the previously found gene at the bht-A1 locus on the short arm of chromosome 2A, we also mapped two new modifier QTL for spike-branching on the short arm of chromosome 1A, termed bht-A2, and 2BS. Using biparental mapping population and GWAS in 302 diverse accessions, the 2BS locus was highly associated with coding sequence variation found at the homoeo-allele of TtBH-B1 (bht-B1). Thus, RILs that combined both bht-A1 and bht-B1 alleles showed additive genetic effects leading to increased penetrance and expressivity of the supernumerary spikelet and/or mini-spike formation.

Hereditary Diffuse Gastric Cancer (HDGC) syndrome is associated with CDH1 germline likely pathogenic/pathogenic variants. Carriers of CDH1 germline likely pathogenic/pathogenic variants are predisposed to diffuse gastric cancer and lobular breast cancer. This study aims to classify the CDH1 c.[715G>A] missense variant identified in a diffuse gastric cancer prone family by performing splicing studies. RT-PCR and subsequent cloning experiments were performed to investigate whether this variant completely disrupts normal splicing. This variant preferentially abolishes normal splicing through activation of a cryptic 3' acceptor splice site within exon 6 of CDH1, presumably leading to a premature protein truncation within first extracellular domain repeat of E-cadherin protein. Our results contributed to evidence necessary to resolve pathogenicity classification of this variant, indicating that this variant is to be classified as pathogenic.

The primary goal of the unrelated population HLA diversity (UPHD) component of the 17th International HLA and Immunogenetics Workshop was to characterize HLA alleles at maximum allelic-resolution in worldwide populations and re-evaluate patterns of HLA diversity across populations. The UPHD project included HLA genotype and sequence data, generated by various next-generation sequencing methods, from 4,240 individuals collated from 12 different countries. Population data included well-defined large datasets from the USA and smaller samples from Europe, Australia, and Western Asia. Allele and haplotype frequencies varied across populations from distant geographical regions. HLA genetic diversity estimated at 2- and 4-field allelic resolution revealed that diversity at the majority of loci, particularly for European-descent populations, was lower at the 2-field resolution. Several common alleles with identical protein sequences differing only by intronic substitutions were found in distinct haplotypes, revealing a more detailed characterization of linkage between variants within the HLA region. The examination of coding and non-coding nucleotide variation revealed many examples in which almost complete biunivocal relations between common alleles at different loci were observed resulting in higher linkage disequilibrium. Our reference data of HLA profiles characterized at maximum resolution from many populations is useful for anthropological studies, unrelated donor searches, transplantation, and disease association studies.

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10-31).

Plasma lipids are known heritable risk factors for cardiovascular disease, but increasing evidence also supports shared genetics with diseases of other organ systems. We devised a comprehensive three-phase framework to identify new lipid-associated genes and study the relationships among lipids, genotypes, gene expression and hundreds of complex human diseases from the Electronic Medical Records and Genomics (347 traits) and the UK Biobank (549 traits). Aside from 67 new lipid-associated genes with strong replication, we found evidence for pleiotropic SNPs/genes between lipids and diseases across the phenome. These include discordant pleiotropy in the HLA region between lipids and multiple sclerosis and putative causal paths between triglycerides and gout, among several others. Our findings give insights into the genetic basis of the relationship between plasma lipids and diseases on a phenome-wide scale and can provide context for future prevention and treatment strategies.

Wiedemann-Steiner syndrome (WSS) is an autosomal dominant disorder caused by monoallelic variants in KMT2A and characterized by intellectual disability and hypertrichosis. We performed a retrospective, multicenter, observational study of 104 individuals with WSS from five continents to characterize the clinical and molecular spectrum of WSS in diverse populations, to identify physical features that may be more prevalent in White versus Black Indigenous People of Color individuals, to delineate genotype-phenotype correlations, to define developmental milestones, to describe the syndrome through adulthood, and to examine clinicians' differential diagnoses. Sixty-nine of the 82 variants (84%) observed in the study were not previously reported in the literature. Common clinical features identified in the cohort included: developmental delay or intellectual disability (97%), constipation (63.8%), failure to thrive (67.7%), feeding difficulties (66.3%), hypertrichosis cubiti (57%), short stature (57.8%), and vertebral anomalies (46.9%). The median ages at walking and first words were 20 months and 18 months, respectively. Hypotonia was associated with loss of function (LoF) variants, and seizures were associated with non-LoF variants. This study identifies genotype-phenotype correlations as well as race-facial feature associations in an ethnically diverse cohort, and accurately defines developmental trajectories, medical comorbidities, and long-term outcomes in individuals with WSS.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.

Prioritizing genes for translation to therapeutics for common diseases has been challenging. Here, we propose an approach to identify drug targets with high probability of success by focusing on genes with both gain of function (GoF) and loss of function (LoF) mutations associated with opposing effects on phenotype (Bidirectional Effect Selected Targets, BEST). We find 98 BEST genes for a variety of indications. Drugs targeting those genes are 3.8-fold more likely to be approved than non-BEST genes. We focus on five genes (IGF1R, NPPC, NPR2, FGFR3, and SHOX) with evidence for bidirectional effects on stature. Rare protein-altering variants in those genes result in significantly increased risk for idiopathic short stature (ISS) (OR = 2.75, p = 3.99 × 10-8). Finally, using functional experiments, we demonstrate that adding an exogenous CNP analog (encoded by NPPC) rescues the phenotype, thus validating its potential as a therapeutic treatment for ISS. Our results show the value of looking for bidirectional effects to identify and validate drug targets.

Autosomal genetic analyses of blood lipids have yielded key insights for coronary heart disease (CHD). However, X chromosome genetic variation is understudied for blood lipids in large sample sizes. We now analyze genetic and blood lipid data in a high-coverage whole X chromosome sequencing study of 65,322 multi-ancestry participants and perform replication among 456,893 European participants. Common alleles on chromosome Xq23 are strongly associated with reduced total cholesterol, LDL cholesterol, and triglycerides (min P = 8.5 × 10-72), with similar effects for males and females. Chromosome Xq23 lipid-lowering alleles are associated with reduced odds for CHD among 42,545 cases and 591,247 controls (P = 1.7 × 10-4), and reduced odds for diabetes mellitus type 2 among 54,095 cases and 573,885 controls (P = 1.4 × 10-5). Although we observe an association with increased BMI, waist-to-hip ratio adjusted for BMI is reduced, bioimpedance analyses indicate increased gluteofemoral fat, and abdominal MRI analyses indicate reduced visceral adiposity. Co-localization analyses strongly correlate increased CHRDL1 gene expression, particularly in adipose tissue, with reduced concentrations of blood lipids.

The International Mouse Phenotyping Consortium reports the generation of new mouse mutant strains for over 5,000 genes, including 2,850 novel null, 2,987 novel conditional- ready, and 4,433 novel reporter alleles.

The most recent updated classification of inborn errors of immunity/primary immunodeficiencies, compiled by the International Union of Immunological Societies Expert Committee, was published in January 2020. Within days of completing this report, it was already out of date, evidenced by the frequent publication of genetic variants proposed to cause novel inborn errors of immunity. As the next formal report from the IUIS Expert Committee will not be published until 2022, we felt it important to provide the community with a brief update of recent contributions to the field of inborn errors of immunity. Herein, we highlight studies that have identified 26 additional monogenic gene defects that reach the threshold to represent novel causes of immune defects.

Studies of acute myeloid leukemia rely on DNA sequencing and immunophenotyping by flow cytometry as primary tools for disease characterization. However, leukemia tumor heterogeneity complicates integration of DNA variants and immunophenotypes from separate measurements. Here we introduce DAb-seq, a technology for simultaneous capture of DNA genotype and cell surface phenotype from single cells at high throughput, enabling direct profiling of proteogenomic states in tens of thousands of cells. To demonstrate the approach, we analyze the disease of three patients with leukemia over multiple treatment timepoints and disease recurrences. We observe complex genotype-phenotype dynamics that illustrate the subtlety of the disease process and the degree of incongruity between blast cell genotype and phenotype in different clinical scenarios. Our results highlight the importance of combined single-cell DNA and protein measurements to fully characterize the heterogeneity of leukemia.

BackgroundThe Chr17q12-21.2 region is the strongest and most consistently associated region with asthma susceptibility. The functional genes or single nucleotide polymorphisms (SNPs) are not obvious due to linkage disequilibrium.ObjectivesWe sought to comprehensively investigate whole-genome sequence and RNA sequence from human bronchial epithelial cells to dissect functional genes/SNPs for asthma severity in the Severe Asthma Research Program.MethodsExpression quantitative trait loci analysis (n = 114), correlation analysis (n = 156) of gene expression and asthma phenotypes, and pathway analysis were performed in bronchial epithelial cells and replicated. Genetic association for asthma severity (426 severe vs 531 nonsevere asthma) and longitudinal asthma exacerbations (n = 273) was performed.ResultsMultiple SNPs in gasdermin B (GSDMB) associated with asthma severity (odds ratio, >1.25) and longitudinal asthma exacerbations (P < .05). Expression quantitative trait loci analyses identified multiple SNPs associated with expression levels of post-GPI attachment to proteins 3, GSDMB, or gasdermin A (3.1 × 10-9 -4). Higher expression levels of GSDMB correlated with asthma and greater number of exacerbations (P < .05). Expression levels of GSDMB correlated with genes involved in IFN signaling, MHC class I antigen presentation, and immune system pathways (false-discovery rate-adjusted P < .05). rs1031458 and rs3902920 in GSDMB colocalized with IFN regulatory factor binding sites and associated with GSDMB expression, asthma severity, and asthma exacerbations (P < .05).ConclusionsBy using a unique set of gene expression data from lung cells obtained using bronchoscopy from comprehensively characterized subjects with asthma, we show that SNPs in GSDMB associated with asthma severity, exacerbations, and GSDMB expression levels. Furthermore, its expression levels correlated with asthma exacerbations and antiviral pathways. Thus, GSDMB is a functional gene for both asthma susceptibility and severity.

Rationale: The 17q12-21.1 locus is one of the most highly replicated genetic associations with asthma. Individuals of African descent have lower linkage disequilibrium in this region, which could facilitate identifying causal variants.Objectives: To identify functional variants at 17q12-21.1 associated with early-onset asthma among African American individuals.Methods: We evaluated African American participants from SAPPHIRE (Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity) (n = 1,940), SAGE II (Study of African Americans, Asthma, Genes and Environment) (n = 885), and GCPD-A (Study of the Genetic Causes of Complex Pediatric Disorders-Asthma) (n = 2,805). Associations with asthma onset at ages under 5 years were meta-analyzed across cohorts. The lead signal was reevaluated considering haplotypes informed by genetic ancestry (i.e., African vs. European). Both an expression-quantitative trait locus analysis and a phenome-wide association study were performed on the lead variant.Measurements and Main Results: The meta-analyzed results from SAPPHIRE, SAGE II, and the GCPD-A identified rs11078928 as the top association for early-onset asthma. A haplotype analysis suggested that the asthma association partitioned most closely with the rs11078928 genotype. Genetic ancestry did not appear to influence the effect of this variant. In the expression-quantitative trait locus analysis, rs11078928 was related to alternative splicing of GSDMB (gasdermin-B) transcripts. The phenome-wide association study of rs11078928 suggested that this variant was predominantly associated with asthma and asthma-associated symptoms.Conclusions: A splice-acceptor polymorphism appears to be a causal variant for asthma at the 17q12-21.1 locus. This variant appears to have the same magnitude of effect in individuals of African and European descent.

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.

The integration of proteomic, transcriptomic, and genetic variant annotation data will improve our understanding of genotype-phenotype associations. Due, in part, to challenges associated with accurate inter-database mapping, such multi-omic studies have not extended to chemoproteomics, a method that measures the intrinsic reactivity and potential "druggability" of nucleophilic amino acid side chains. Here, we evaluated mapping approaches to match chemoproteomic-detected cysteine and lysine residues with their genetic coordinates. Our analysis revealed that database update cycles and reliance on stable identifiers can lead to pervasive misidentification of labeled residues. Enabled by this examination of mapping strategies, we then integrated our chemoproteomics data with computational methods for predicting genetic variant pathogenicity, which revealed that codons of highly reactive cysteines are enriched for genetic variants that are predicted to be more deleterious and allowed us to identify and functionally characterize a new damaging residue in the cysteine protease caspase-8. Our study provides a roadmap for more precise inter-database mapping and points to untapped opportunities to improve the predictive power of pathogenicity scores and to advance prioritization of putative druggable sites.

Along with early-onset cancers, multiple primary cancers (MPCs) are likely resulting from increased genetic susceptibility; however, the associated predisposition genes or prevalence of the pathogenic variants genes in MPC patients are often unknown. We screened 71 patients with MPC of the stomach, colorectal, and endometrium, sequencing 65 cancer predisposition genes. A subset of 19 patients with early-onset MPC of stomach and colorectum were further evaluated for variants in cancer related genes using both normal and tumor whole exome sequencing. Among 71 patients with MPCs, variants classified to be pathogenic were observed in 15 (21.1%) patients and affected Lynch Syndrome (LS) genes: MLH1 (n = 10), MSH6 (n = 2), PMS2 (n = 2), and MSH2 (n = 1). All carriers had tumors with high microsatellite instability and 13 of them (86.7%) were early-onset, consistent with LS. In 19 patients with early-onset MPCs, loss of function (LoF) variants in RECQL5 were more prevalent in non-LS MPC than in matched sporadic cancer patients (OR = 31.6, 2.73-1700.6, p = 0.001). Additionally, there were high-confidence LoF variants at FANCG and CASP8 in two patients accompanied by somatic loss of heterozygosity in tumor, respectively. The results suggest that genetic screening should be considered for synchronous cancers and metachronous MPCs of the LS tumor spectrum, particularly in early-onset. Susceptibility variants in non-LS genes for MPC patients may exist, but evidence for their role is more elusive than for LS patients.

ObjectiveCertain copy number variants (CNVs) greatly increase the risk of autism. The authors conducted a genetics-first study to investigate whether heterogeneity in the clinical presentation of autism is underpinned by specific genotype-phenotype relationships.MethodsThis international study included 547 individuals (mean age, 12.3 years [SD=4.2], 54% male) who were ascertained on the basis of having a genetic diagnosis of a rare CNV associated with high risk of autism (82 16p11.2 deletion carriers, 50 16p11.2 duplication carriers, 370 22q11.2 deletion carriers, and 45 22q11.2 duplication carriers), as well as 2,027 individuals (mean age, 9.1 years [SD=4.9], 86% male) with autism of heterogeneous etiology. Assessments included the Autism Diagnostic Interview-Revised and IQ testing.ResultsThe four genetic variant groups differed in autism symptom severity, autism subdomain profile, and IQ profile. However, substantial variability was observed in phenotypic outcome in individual genetic variant groups (74%-97% of the variance, depending on the trait), whereas variability between groups was low (1%-21%, depending on the trait). CNV carriers who met autism criteria were compared with individuals with heterogeneous autism, and a range of profile differences were identified. When clinical cutoff scores were applied, 54% of individuals with one of the four CNVs who did not meet full autism diagnostic criteria had elevated levels of autistic traits.ConclusionsMany CNV carriers do not meet full diagnostic criteria for autism but nevertheless meet clinical cutoffs for autistic traits. Although profile differences between variants were observed, there is considerable variability in clinical symptoms in the same variant.