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.

Background Worldwide, respiratory syncytial virus (RSV) infections are among the most common causes of infant hospitalization. Host genetic factors influencing the risk and severity of RSV infection are not well known. Methods We conducted a genome-wide association study (GWAS) to investigate single-nucleotide polymorphisms (SNPs) associated with severe RSV infections using a nested case-control design based on 2 Danish cohorts. We compared SNPs from 1786 children hospitalized with RSV to 45 060 controls without an RSV-coded hospitalization. We performed gene-based testing, tissue enrichment, gene-set enrichment, and a meta-analysis of the 2 cohorts. Finally, an analysis of potential associations between the severity of RSV infection and genetic markers was performed. Results We did not detect any significant genome-wide associations between SNPs and RSV infection or the severity of RSV. We did find potential loci associated with RSV infections on chromosome 5 in 1 cohort but failed to replicate any signals in both cohorts. Conclusions Despite being the largest GWAS of severe RSV infection, we did not detect any genome-wide significant loci. This may be an indication of a lack of power or an absence of signal. Future studies might include mild illness and need to be larger to detect any significant associations. [ABSTRACT FROM AUTHOR]

Copyright of Arquivos Brasileiros de Oftalmologia is the property of Arquivos Brasileiros de Oftalmologia and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

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.

Introduction Neuronal ceroid lipofuscinoses (CLNs) are a group of lysosomal storage disorders of genetic origin, characterized by progressive neurodegeneration and intracellular accumulation of autofluorescent lipopigment. Thirteen genes related to CLNs are currently described, showing genetic and allelic heterogeneity, most of them with an autosomal recessive pattern. Due to the few descriptions of cases related to CLNs in Brazil, it is necessary to describe the phenotypic and genotypic characteristics of these patients. This study aims to evaluate the genotypic profile and correlate it with the phenotypic characteristics of patients with CLN in a children's hospital. Methods This study was performed as a descriptive cross-sectional study with analysis of medical records, imaging, and laboratory tests of patients who had a confirmed molecular diagnosis of CLN. Results The sample consisted of 11 patients from nine families with different subtypes of CLNs (CLN2, 5, 6, 7, and 8), with CLN2 being the most prevalent in the study. A total of 16 mutation variants were identified in genes associated with the five CLNs described in this study, with typical and atypical clinical phenotypes depending on the subtype and its variants. Conclusion Novel mutations identified in the patients in this study showed phenotypes of rapid and severe progression in the CLN2 patient and similar characteristics in CLN6 and CLN7 patients, as previously described in the literature. [ABSTRACT FROM AUTHOR]

Background/Purpose: Congenital nephrotic syndrome (CNS) is one of the important causes of end-stage kidney disease in children. Studies on the genotype, phenotype, and clinical outcome in infants with CNS caused by genetic mutations are scarce. Methods: We analyzed the genetic background, clinical manifestations, treatment response, and prognosis of pediatric patients with CNS in Taiwan. Results: Fifteen infants with CNS were enrolled, and 11 patients of median age 21 (interquartile range 3∼44) days caused by genetic mutations from 10 unrelated families were included in the study. Of the eleven patients, 9 had extra-renal manifestations including microcephaly, facial dysmorphism, and skeletal anomalies. More than two-thirds of the patients had disease onset before 1 month of age. Diffuse meningeal sclerosis was the most common histological characteristic. Whole exome sequencing followed by direct Sanger sequence revealed mutations in OSGEP (R247Q), WT1 (R366H and R467Q), LAMB2 (Q1209∗ and c. 5432–5451 19 bp deletion), NUP93 (D302V), and LAGE3 (c.188+1G > A). Three of the variants were novel. Corticosteroids and/or immunosuppressants were administered in 2 patients, but both were refractory to treatment. During the mean 3.5 years of follow-up, all but two died of uremia and sepsis. The two survivors reached end-stage kidney disease and required peritoneal dialysis, and one of them underwent uneventful renal transplantation. Conclusion: The majority of patients with CNS in Taiwan were caused by OSGEP followed by WT1 mutation. R247Q is the hotspot mutation of OSGEP in Taiwan. CNS patients in Taiwan suffer from significant morbidity and mortality.

Congenital nephrotic syndrome (CNS) is one of the important causes of end-stage kidney disease in children. Studies on the genotype, phenotype, and clinical outcome in infants with CNS caused by genetic mutations are scarce. We analyzed the genetic background, clinical manifestations, treatment response, and prognosis of pediatric patients with CNS in Taiwan. Fifteen infants with CNS were enrolled, and 11 patients of median age 21 (interquartile range 3∼44) days caused by genetic mutations from 10 unrelated families were included in the study. Of the eleven patients, 9 had extra-renal manifestations including microcephaly, facial dysmorphism, and skeletal anomalies. More than two-thirds of the patients had disease onset before 1 month of age. Diffuse meningeal sclerosis was the most common histological characteristic. Whole exome sequencing followed by direct Sanger sequence revealed mutations in OSGEP (R247Q), WT1 (R366H and R467Q), LAMB2 (Q1209∗ and c. 5432–5451 19 bp deletion), NUP93 (D302V), and LAGE3 (c.188+1G > A). Three of the variants were novel. Corticosteroids and/or immunosuppressants were administered in 2 patients, but both were refractory to treatment. During the mean 3.5 years of follow-up, all but two died of uremia and sepsis. The two survivors reached end-stage kidney disease and required peritoneal dialysis, and one of them underwent uneventful renal transplantation. The majority of patients with CNS in Taiwan were caused by OSGEP followed by WT1 mutation. R247Q is the hotspot mutation of OSGEP in Taiwan. CNS patients in Taiwan suffer from significant morbidity and mortality. [ABSTRACT FROM AUTHOR]

Background: Periodontitis and benign prostatic hyperplasia (BPH) are all common chronic diseases with higher incidence in middle-aged and old men. Several studies have indicated a potential association between periodontitis and BPH, although the findings remain inconclusive. However, there is no mendelian randomization (MR) studies to assess this association. Methods: The 40 men who had received health check-ups were included in an epidemiological study. Genetic data of BPH (13118 cases and 72799 controls) and periodontitis (3046 cases and 195395 controls) from FinnGen project was used to perform two-sample MR analysis. The inverse-variance weighted (IVW) model was identified as the primary analytical method, with MR Egger, weighted median, simple mode, and weighted mode serving as additional approaches. Results: The epidemiological analysis demonstrated a lack of statistically significant differences in the prevalence of clinical BPH between severe periodontitis group and non-severe periodontitis group. Similarly, no statistically significant differences were found in the prevalence of severe periodontitis among individuals with clinical BPH compared to those without. Additionally, Among the five models utilized in MR analysis, including the IVW model, no evidence of a causal link between periodontitis and BPH was observed. Conclusion: The findings from our epidemiological investigation and MR analysis do not provide support for a causal relationship between periodontitis and BPH. [ABSTRACT FROM AUTHOR]

Identifying genotype-by-environment interaction (GEI) is challenging because the GEI analysis generally has low power. Large-scale consortium-based studies are ultimately needed to achieve adequate power for identifying GEI. We introduce Multi-Trait Analysis of Gene–Environment Interactions (MTAGEI), a powerful, robust, and computationally efficient framework to test gene–environment interactions on multiple traits in large data sets, such as the UK Biobank (UKB). To facilitate the meta-analysis of GEI studies in a consortium, MTAGEI efficiently generates summary statistics of genetic associations for multiple traits under different environmental conditions and integrates the summary statistics for GEI analysis. MTAGEI enhances the power of GEI analysis by aggregating GEI signals across multiple traits and variants that would otherwise be difficult to detect individually. MTAGEI achieves robustness by combining complementary tests under a wide spectrum of genetic architectures. We demonstrate the advantages of MTAGEI over existing single-trait-based GEI tests through extensive simulation studies and the analysis of the whole exome sequencing data from the UKB. [ABSTRACT FROM AUTHOR]

Innate and adaptive immunity play a crucial role in allergic bronchopulmonary aspergillosis (ABPA) pathogenesis. We performed next-generation sequencing using the Illumina TruSight One panel (4,811 human disease-associated genes, at least 20 × coverage) and selected 22 known immune genes (toll-like receptors (TLRs), C-type lectin, interleukin-4 receptor, and others). We included ABPA (n = 18), asthma without ABPA (n = 12), and healthy controls (n = 8). We analyzed 3011 SNPs from 22 genes and identified 145 SNPs (13 genes) that were present only in the disease groups and absent in controls. The SNP frequency overall was significantly higher in ABPA than in asthmatics (89/145 [61.4%] vs. 56/145 [38.6%], p = 0.0001). The SNP frequency in the TLR10 gene was also significantly higher in ABPA than in asthma (p = 0.017). Association analysis further revealed three genes having significant associations. Of these, NOS3 and HLA-DQB1 are associated with antimicrobial activity and adaptive immunity. More extensive studies are required to confirm our findings. [ABSTRACT FROM AUTHOR]

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.

Aims: We evaluated the involvement of lncRNAs in the development of pathologies associated with chronic hyperglycaemia in rat models in a model of type 1, type 2 and gestational diabetes. Methods: Reports were searched in Dialnet, Scielo, HINARI, Springer, ClinicalKey, OTseeker, PubMed and different grey literature databases with any restrictions. Bibliography databases will be searched from their inception to December 2022. Results: Thirty‐seven studies met our criteria, and they had the following characteristics: original experimental studies on diabetes, the lncRNAs were extracted or measured from tissues of specific areas and the results were expressed in terms of standard measures by RT‐PCR. In most studies, both primary and secondary outcomes were mentioned. On the other hand, we found a total of nine diabetic complications, being retinopathy, nephropathy and neuropathy the most representatives. Additionally, it was found that MALAT1, H19, NEAT1 and TUG1 are the most studied lncRNAs about these complications in rats. On the other hand, the lncRNAs with the highest rate of change were MSTRG.1662 (17.85; 13.78, 21.93), ENSRNOT00000093120_Aox3 (7.13; 5.95, 8.31) and NONRATG013497.2 (−5.55; −7.18, −3.93). Conclusions: This review found a significant involvement of lncRNAs in the progression of pathologies associated with chronic hyperglycaemia in rat models, and further studies are needed to establish their potential as biomarkers and therapeutic targets for diabetes. [ABSTRACT FROM AUTHOR]

Introduction: Joubert syndrome (JS) arises from defects of primary cilia resulting in potential malformations of the brain, kidneys, eyes, liver, and limbs. Several of the 35+ genes associated with JS have recognized genotype/phenotype correlations, but most genes have not had enough reported individuals to draw meaningful conclusions. Methods: A PubMed literature review identified 688 individuals with JS across 32 genes and 112 publications to bolster known genotype/phenotype relationships and identify new correlations. All included patients had the "molar tooth sign" and a confirmed genetic diagnosis. Individuals were categorized by age, ethnicity, sex and the presence of developmental disability/intellectual disability, hypotonia, abnormal eye movements, ataxia, visual impairment, renal impairment, polydactyly, and liver abnormalities. Results: Most genes demonstrated unique phenotypic profiles. Grouping proteins based on physiologic interactions established stronger phenotypic relationships that reflect known ciliary pathophysiology. Age‐stratified data demonstrated that end‐organ disease is progressive in JS. Most genes demonstrated a significant skew towards having variants with either residual protein function or no residual protein function. Conclusion: This cohort demonstrates that clinically meaningful genotype/phenotype relationships exist within most JS‐related genes and can be referenced to allow for more personalized clinical care. [ABSTRACT FROM AUTHOR]

Currently, only the general control of the risk factors is known to prevent lacunar cerebral infarction, but it is unknown which type of medication for controlling the risk factors has a causal relationship with reducing the risk of lacunar infarction. To unlock this medical mystery, drug-target Mendelian randomization analysis was applied to estimate the effect of common antihypertensive agents, hypolipidemic agents, and hypoglycemic agents on lacunar stroke. Lacunar stroke data for the transethnic analysis were derived from meta-analyses comprising 7338 cases and 254,798 controls. We have confirmed that genetic variants mimicking calcium channel blockers were found to most stably prevent lacunar stroke. The genetic variants at or near HMGCR, NPC1L1, and APOC3 were predicted to decrease lacunar stroke incidence in drug-target MR analysis. These variants mimic the effects of statins, ezetimibe, and antisense anti-apoC3 agents, respectively. Genetically proxied GLP1R agonism had a marginal effect on lacunar stroke, while a genetically proxied improvement in overall glycemic control was associated with reduced lacunar stroke risk. Here, we show that certain categories of drugs currently used in clinical practice can more effectively reduce the risk of stroke. Repurposing several drugs with well-established safety and low costs for lacunar stroke prevention should be given high priority when doctors are making decisions in clinical practice. This may contribute to healthier brain aging. [ABSTRACT FROM AUTHOR]

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 The clinical heterogeneity in major depressive disorder (MDD), variable treatment response, and conflicting findings limit the ability of genomics toward the discovery of evidence-based diagnosis and treatment regimen. This study attempts to curate all genetic association findings to evaluate potential variants for clinical translation. Methods We systematically reviewed all candidates and genome-wide association studies for both MDD susceptibility and antidepressant response, independently, using MEDLINE, particularly to identify replicated findings. These variants were evaluated for functional consequences using different in silico tools and further estimated their diagnostic predictability by calculating positive predictive values. Results A total of 217 significantly associated studies comprising 1200 variants across 545 genes and 128 studies including 921 variants across 412 genes were included with MDD susceptibility and antidepressant response, respectively. Although the majority of associations were confirmed by a single study, we identified 31 and 18 replicated variants (in at least 2 studies) for MDD and antidepressant response. Functional annotation of these 31 variants predicted 20% coding variants as deleterious/damaging and 80.6% variants with regulatory effect. Similarly, the response-related 18 variants revealed 25% coding variant as damaging and 88.2% with substantial regulatory potential. Finally, we could calculate the diagnostic predictability of 19 and 5 variants whose positive predictive values ranges from 0.49 to 0.66 for MDD and 0.36 to 0.66 for response. Conclusions The replicated variants presented in our data are promising for disease diagnosis and improved response outcomes. Although these quantitative assessment measures are solely directive of available observational evidence, robust homogenous validation studies are required to strengthen these variants for molecular diagnostic application. [ABSTRACT FROM AUTHOR]

Background Two new SNPs have been recently associated to Alzheimer’s disease in African American populations: FCGRIIB rs1050501 C/T, and PILRA rs1859788 A/G. The risk of Alzheimer’s disease in FCGRIIB C and PILRA A allele carriers is three times higher than in non-carriers. However, the association between these and other single nucleotide polymorphisms (SNPs) has not been assessed. Methods Linkage disequilibrium analysis, with r= 0.8 as a threshold value, was used to impute new candidate SNPs, on genomic data from both genes in 26 populations worldwide (n= 2504) from the 1000Genomes database. Results Four SNPs (rs13376485, rs3767640, rs3767639 and rs3767641) were linked to rs1050501 and one (rs2405442) to rs1859788 in the whole sample. Conclusions Five novel SNPs could be associated with Alzheimer’s disease susceptibility and play a causal role, even if none of them are exon variants since their potential roles in the regulation of gene expression.

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.

Abstract Background Observational studies have shown an association between higher bilirubin levels and improved respiratory health outcomes. Targeting higher bilirubin levels has been proposed as a novel therapeutic strategy in COPD. However, bilirubin levels are influenced by multiple intrinsic and extrinsic factors, and these observational studies are prone to confounding. Genetic analyses are one approach to overcoming residual confounding in observational studies. Objectives To test associations between a genetic determinant of bilirubin levels and respiratory health outcomes. Methods COPDGene participants underwent genotyping at the baseline visit. We confirmed established associations between homozygosity for rs6742078 and higher bilirubin, and between higher bilirubin and decreased risk of acute respiratory events within this cohort. For our primary analysis, we used negative binomial regression to test associations between homozygosity for rs6742078 and rate of acute respiratory events. Results 8,727 participants (n = 6,228 non-Hispanic white and 2,499 African American) were included. Higher bilirubin was associated with decreased rate of acute respiratory events [incidence rate ratio (IRR) 0.85, 95% CI 0.75 to 0.96 per SD increase in bilirubin intensity]. We did not find significant associations between homozygosity for rs6742078 and acute respiratory events (IRR 0.94, 95% CI 0.70 to 1.25 for non-Hispanic white and 1.09, 95% CI 0.91 to 1.31 for African American participants). Conclusions A genetic determinant of higher bilirubin levels was not associated with better respiratory health outcomes. These results do not support targeting higher bilirubin levels as a therapeutic strategy in COPD.

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.

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).