Guidelines for investigating causality of sequence variants in human disease

Recessive variants in MYO1C as a potential novel cause of proteinuric kidney disease

Background

Steroid-resistant nephrotic syndrome is the second leading cause of chronic kidney disease among patients <25 years of age. Through whole exome sequencing, identification of >65 monogenic causes has rendered insights into disease mechanisms of nephrotic syndrome.

Methods

To elucidate novel monogenic causes of NS, we combined homozygosity mapping with ES in a worldwide cohort of 1649 pediatric patients with NS.

Results

We identified homozygous missense variants in MYO1C in two unrelated children with nephrotic syndrome (c.292C>T, p.R98W; c.2273 A>T, p.K758M). We evaluated publicly available kidney single-cell RNA sequencing datasets and found MYO1Cto be predominantly expressed in podocytes. We then performed structural modeling in molecular viewer PyMol using the super function aligning shared regions within both partial structures of MYO1C (4byf and 4r8g). In both structures, calmodulin, a common regulator of myosin activity, is shown to bind to the IQ motif. At both residue sites (K758; R98), there are ion-ion interactions stabilizing intradomain and ligand interactions: R98 binds to nearby D220 within the Myosin Motor Domain and K758 binds to E14 on a calmodulin molecule. Variants of these charged residues to non-charged amino acids could ablate these ionic interactions, weakening protein structure and function establishing the impact of these variants.

Conclusion

We here identified recessive variants in MYO1C as a potential novel cause of nephrotic syndrome in children.

DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

Podocyte-specific Nup160 knockout mice develop nephrotic syndrome and glomerulosclerosis

More than 60 monogenic genes mutated in steroid-resistant nephrotic syndrome (SRNS) have been identified. Our previous study found that mutations in nucleoporin 160 kD (NUP160) are implicated in SRNS. The NUP160 gene encodes a component of the nuclear pore complex. Recently, two siblings with homozygous NUP160 mutations presented with SRNS and a nervous system disorder. However, replication of nephrotic syndrome (NS)-associated phenotypes in a mammalian model following loss of Nup160 is needed to prove that NUP160 mutations cause SRNS. Here, we generated a podocyte-specific Nup160 knockout (Nup160podKO) mouse model using CRISPR/Cas9 and Cre/loxP technologies. We investigated NS-associated phenotypes in these Nup160podKO mice. We verified efficient abrogation of Nup160 in Nup160podKO mice at both the DNA and protein levels. We showed that Nup160podKO mice develop typical signs of NS. Nup160podKO mice exhibited progression of proteinuria to average albumin/creatinine ratio (ACR) levels of 15.06 ± 2.71 mg/mg at 26 weeks, and had lower serum albumin levels of 13.13 ± 1.34 g/l at 30 weeks. Littermate control mice had urinary ACR mean values of 0.03 mg/mg and serum albumin values of 22.89 ± 0.34 g/l at the corresponding ages. Further, Nup160podKO mice exhibited glomerulosclerosis compared with littermate control mice. Podocyte-specific Nup160 knockout in mice led to NS and glomerulosclerosis. Thus, our findings strongly support that mutations in NUP160 cause SRNS. The newly generated Nup160podKO mice are a reliable mammalian model for future study of the pathogenesis of NUP160-associated SRNS.

Implication of transcription factor FOXD2 dysfunction in syndromic congenital anomalies of the kidney and urinary tract (CAKUT)

Congenital anomalies of the kidney and urinary tract (CAKUT) are the predominant cause for chronic kidney disease below age 30 years. Many monogenic forms have been discovered due to comprehensive genetic testing like exome sequencing. However, disease-causing variants in known disease-associated genes only explain a proportion of cases. Here, we aim to unravel underlying molecular mechanisms of syndromic CAKUT in three unrelated multiplex families with presumed autosomal recessive inheritance. Exome sequencing in the index individuals revealed three different rare homozygous variants in FOXD2, encoding a transcription factor not previously implicated in CAKUT in humans: a frameshift in the Arabic and a missense variant each in the Turkish and the Israeli family with segregation patterns consistent with autosomal recessive inheritance. CRISPR/Cas9-derived Foxd2 knockout mice presented with a bilateral dilated kidney pelvis accompanied by atrophy of the kidney papilla and mandibular, ophthalmologic, and behavioral anomalies, recapitulating the human phenotype. In a complementary approach to study pathomechanisms of FOXD2-dysfunction-mediated developmental kidney defects, we generated CRISPR/Cas9-mediated knockout of Foxd2 in ureteric bud-induced mouse metanephric mesenchyme cells. Transcriptomic analyses revealed enrichment of numerous differentially expressed genes important for kidney/urogenital development, including Pax2 and Wnt4 as well as gene expression changes indicating a shift toward a stromal cell identity. Histology of Foxd2 knockout mouse kidneys confirmed increased fibrosis. Further, genome-wide association studies suggest that FOXD2 could play a role for maintenance of podocyte integrity during adulthood. Thus, our studies help in genetic diagnostics of monogenic CAKUT and in understanding of monogenic and multifactorial kidney diseases.

Clinical-Genomic Analysis of 1261 Patients with Ehlers-Danlos Syndrome Outlines an Articulo-Autonomic Gene Network (Entome)

Systematic evaluation of 80 history and 40 history findings diagnosed 1261 patients with Ehlers-Danlos syndrome (EDS) by direct or online interaction, and 60 key findings were selected for their relation to clinical mechanisms and/or management. Genomic testing results in 566 of these patients supported EDS relevance by their differences from those in 82 developmental disability patients and by their association with general rather than type-specific EDS findings. The 437 nuclear and 79 mitochondrial DNA changes included 71 impacting joint matrix (49 COL5), 39 bone (30 COL1/2/9/11), 22 vessel (12 COL3/8VWF), 43 vessel-heart (17FBN1/11TGFB/BR), 59 muscle (28 COL6/12), 56 neural (16 SCN9A/10A/11A), and 74 autonomic (13 POLG/25porphyria related). These genes were distributed over all chromosomes but the Y, a network analogized to an 'entome' where DNA change disrupts truncal mechanisms (skin constraint, neuromuscular support, joint vessel flexibility) and produces a mirroring cascade of articular and autonomic symptoms. The implied sequences of genes from nodal proteins to hypermobility to branching tissue laxity or dysautonomia symptoms would be ideal for large language/artificial intelligence analyses.

Discovery of novel genetic syndromes in Latin America: Opportunities and challenges

Latin America (LatAm) has a rich and historically significant role in delineating both novel and well-documented genetic disorders. However, the ongoing advancements in the field of human genetics pose challenges to the relatively slow adaption of LatAm in the field. Here, we describe past and present contributions of LatAm to the discovery of novel genetic disorders, often referred as novel gene-disease associations (NGDA). We also describe the current methodologies for discovery of NGDA, taking into account the latest developments in genomics. We provide an overview of opportunities and challenges for NGDA research in LatAm considering the steps currently performed to identify and validate such associations. Given the multiple and diverse needs of populations and countries in LatAm, it is imperative to foster collaborations amongst patients, indigenous people, clinicians and scientists. Such collaborative effort is essential for sustaining and enhancing the LatAm´s contributions to the field of NGDA.

Genetic variant interpretation for the neurologist - A pragmatic approach in the next-generation sequencing era in childhood epilepsy

Genetic advances over the past decade have enhanced our understanding of the genetic landscape of childhood epilepsy. However a major challenge for clinicians ha been understanding the rationale and systematic approach towards interpretation of the clinical significance of variant(s) detected in their patients. As the clinical paradigm evolves from gene panels to whole exome or whole genome testing including rapid genome sequencing, the number of patients tested and variants identified per patient will only increase. Each step in the process of variant interpretation has limitations and there is no single criterion which enables the clinician to draw reliable conclusions on a causal relationship between the variant and disease without robust clinical phenotyping. Although many automated online analysis software tools are available, these carry a risk of misinterpretation. This guideline provides a pragmatic, real-world approach to variant interpretation for the child neurologist. The focus will be on ascertaining aspects such as variant frequency, subtype, inheritance pattern, structural and functional consequence with regard to genotype-phenotype correlations, while refraining from mere interpretation of the classification provided in a genetic test report. It will not replace the expert advice of colleagues in clinical genetics, however as genomic investigations become a first-line test for epilepsy, it is vital that neurologists and epileptologists are equipped to navigate this landscape.

A Multilayered Post-Genome-Wide Association Study Analysis Pipeline Defines Functional Variants and Target Genes for Systemic Lupus Erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWASs) have revealed multiple SLE susceptibility loci and associated single-nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown.

METHODS

Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci (cis- and trans-expression quantitative trait loci) with promoter capture high-throughput capture chromatin conformation (PCHi-C), allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, quantitative polymerase chain reaction, and Western blot.

RESULTS

Anchoring on 452 index SNPs, we selected 9,931 high linkage disequilibrium (r2 > 0.8) SNPs and defined 182 independent non-human leukocyte antigen (HLA) SLE loci. The 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated reference SNP identifier 57668933 (rs57668933) as a functional variant regulating multiple targets, including SLE-risk gene ELF1 in B cells.

CONCLUSION

We demonstrate and validate post-GWAS strategies for using multidimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms to guide experimental characterization.

The first case of intellectual disability caused by novel compound heterozygosity for NUDT2 variants

NUDT2 is an enzyme important for maintaining the intracellular level of the diadenosine tetraphosphate (Ap4A). Bi-allelic loss of function variants in NUDT2 has recently been reported as a rare cause of intellectual disability (ID). Herein, we describe a Chinese girl with ID, attention deficit hyperactivity disorder (ADHD), and motor delays with abnormal walking posture and difficulty climbing stairs, who bears compound heterozygous variants c.34 C > T (p.R12*) and c.194T > G (p.I65R) in NUDT2. Homozygous variants c.34 C > T (p.R12*) or c.186del (p.A63Qfs*3) in NUDT2 were previously reported to cause ID. This is the first patient with ID due to compound heterozygous variants in NUDT2 and p.I65R is a novel missense variant. This study enriched the genotype and phenotype of NUDT2-related ID and supported the critical developmental involvement of NUDT2.

Effects of LDLR variants rs5928, rs750518671 and rs879254797 on protein structure and functional activity in HepG2 cells transfected with CRISPR/Cas9 constructs

Familial Hypercholesterolemia (FH) is a genetic disorder associated with premature atherosclerosis and increased risk of cardiovascular diseases. LDLR deleterious mutations are associated with FH, however the role of some missense variants in FH pathogenicity remains to be elucidated. This study explored the predictive impact of LDLR missense variants on protein structure and investigated their functional effects on LDLR expression in HepG2 cells transfected with CRISPR/Cas9 constructs. FH (n = 287) and non-FH patients (n = 45) were selected, and lipid profile was obtained from medical records. LDLR variants were identified using an exon-targeted gene sequencing strategy, considering its cost-effective to increase accuracy in the identification step of the most likely FH-related variants in a less laborious process. LDLR variants were selected based on conflicting pathogenicity results found in Clinvar, in silico prediction tools, affected LDLR domains, and less common variants considering minor allele frequency < 0.05. Molecular modeling studies were used to predict the effects of LDLR missense variants on protein structure. Recombinant LDLR variants were constructed using CRISPR/Cas9 system and were used to transfect HepG2 cells. Functional assays in transfected cells were performed to assess LDLR expression using flow cytometry and western blotting, and LDLR activity using flow cytometry and confocal microscopy. The variants rs121908039 (c.551G>A, p.C184Y), rs879254797 (c.1118G>A, p.G373D), rs28941776 (c.1646G>A, p.G549D), rs750518671 (c.2389G>C, p.V797L), rs5928 (c.2441G>A, p.R814Q) and rs137853964 (c.2479G>A, p.V827I) were selected for molecular docking analysis. The p.C184Y exhibited a favorable energy change for protein stability due to its interaction with EGF-A/EGF-B regions; p.G373D and p.G549D displayed intermediate energy changes; and p.R814Q and p.V827I showed smaller energy changes. The results of functional assays showed that p.G373D, p.V797L and p.R814Q reduced LDLR expression and activity (p < 0.05). Microscopic analysis of the p.V797L and p.G373D variants revealed altered lipid localization and accumulation in transfected HepG2 cells. Carriers of p.G549D, p.V797L and p.R814Q had higher LDL cholesterol levels than non-FH group, and (p < 0.05). p.G373D and p.G549D were associated with clinical manifestations of FH. In conclusion, the p.C184Y, p.G373D, p.G549D and p.R814Q variants alter protein stability and intramolecular interactions, while p.V797L has a minimal impact on protein stability, and p.V827I has no significant intramolecular interactions. p.G373D, p.V767L and p.R814Q are associated with impaired LDLR expression and activity.

A new variant in the ZCCHC8 gene: diverse clinical phenotypes and expression in the lung

Introduction

Pulmonary fibrosis is a severe disease which can be familial. A genetic cause can only be found in ∼40% of families. Searching for shared novel genetic variants may aid the discovery of new genetic causes of disease.

Methods

Whole-exome sequencing was performed in 152 unrelated patients with a suspected genetic cause of pulmonary fibrosis from the St Antonius interstitial lung disease biobank. Variants of interest were selected by filtering for novel, potentially deleterious variants that were present in at least three unrelated pulmonary fibrosis patients.

Results

The novel c.586G>A p.(E196K) variant in the ZCCHC8 gene was observed in three unrelated patients: two familial patients and one sporadic patient, who was later genealogically linked to one of the families. The variant was identified in nine additional relatives with pulmonary fibrosis and other telomere-related phenotypes, such as pulmonary arterial venous malformations, emphysema, myelodysplastic syndrome, acute myeloid leukaemia and dyskeratosis congenita. One family showed incomplete segregation, with absence of the variant in one pulmonary fibrosis patient who carried a PARN variant. The majority of ZCCHC8 variant carriers showed short telomeres in blood. ZCCHC8 protein was located in different lung cell types, including alveolar type 2 (AT2) pneumocytes, the culprit cells in pulmonary fibrosis. AT2 cells showed telomere shortening and increased DNA damage, which was comparable to patients with sporadic pulmonary fibrosis and those with pulmonary fibrosis carrying a telomere-related gene variant, respectively.

Discussion

The ZCCHC8 c.586G>A variant confirms the involvement of ZCCHC8 in pulmonary fibrosis and short-telomere syndromes and underlines the importance of including the ZCCHC8 gene in diagnostic gene panels for these diseases.

Dissecting the Immune System through Gene Regulation

The immune system plays a dual role in human health, functioning both as a protector against pathogens and, at times, as a contributor to disease. This feature emphasizes the importance to uncover the underlying causes of its malfunctions, necessitating an in-depth analysis in both pathological and physiological conditions to better understand the immune system and immune disorders. Recent advances in scientific technology have enabled extensive investigations into gene regulation, a crucial mechanism governing cellular functionality. Studying gene regulatory mechanisms within the immune system is a promising avenue for enhancing our understanding of immune cells and the immune system as a whole. The gene regulatory mechanisms, revealed through various methodologies, and their implications in the field of immunology are discussed in this chapter.

A probabilistic graphical model for estimating selection coefficient of missense variants from human population sequence data

Accurately predicting the effect of missense variants is a central problem in interpretation of genomic variation. Commonly used computational methods does not capture the quantitative impact on fitness in populations. We developed MisFit to estimate missense fitness effect using biobank-scale human population genome data. MisFit jointly models the effect at molecular level ( d ) and population level (selection coefficient, s ), assuming that in the same gene, missense variants with similar d have similar s . MisFit is a probabilistic graphical model that integrates deep neural network components and population genetics models efficiently with inductive bias based on biological causality of variant effect. We trained it by maximizing probability of observed allele counts in 236,017 European individuals. We show that s is informative in predicting frequency across ancestries and consistent with the fraction of de novo mutations given s . Finally, MisFit outperforms previous methods in prioritizing missense variants in individuals with neurodevelopmental disorders.

The effects of pathogenic and likely pathogenic variants for inherited hemostasis disorders in 140 214 UK Biobank participants

Rare genetic diseases affect millions, and identifying causal DNA variants is essential for patient care. Therefore, it is imperative to estimate the effect of each independent variant and improve their pathogenicity classification. Our study of 140 214 unrelated UK Biobank (UKB) participants found that each of them carries a median of 7 variants previously reported as pathogenic or likely pathogenic. We focused on 967 diagnostic-grade gene (DGG) variants for rare bleeding, thrombotic, and platelet disorders (BTPDs) observed in 12 367 UKB participants. By association analysis, for a subset of these variants, we estimated effect sizes for platelet count and volume, and odds ratios for bleeding and thrombosis. Variants causal of some autosomal recessive platelet disorders revealed phenotypic consequences in carriers. Loss-of-function variants in MPL, which cause chronic amegakaryocytic thrombocytopenia if biallelic, were unexpectedly associated with increased platelet counts in carriers. We also demonstrated that common variants identified by genome-wide association studies (GWAS) for platelet count or thrombosis risk may influence the penetrance of rare variants in BTPD DGGs on their associated hemostasis disorders. Network-propagation analysis applied to an interactome of 18 410 nodes and 571 917 edges showed that GWAS variants with large effect sizes are enriched in DGGs and their first-order interactors. Finally, we illustrate the modifying effect of polygenic scores for platelet count and thrombosis risk on disease severity in participants carrying rare variants in TUBB1 or PROC and PROS1, respectively. Our findings demonstrate the power of association analyses using large population datasets in improving pathogenicity classifications of rare variants.

Changes in Treg and Breg cells in a healthy pediatric population

The interpretation of clinical diagnostic results in suspected inborn errors of immunity, including Tregopathies, is hampered by the lack of age-stratified reference values for regulatory T cells (Treg) in the pediatric population and a consensus on which Treg immunophenotype to use. Regulatory B cells (Breg) are an important component of the regulatory system that have been poorly studied in the pediatric population. We analyzed (1) the correlation between the three immunophenotypic definitions of Treg (CD4+CD25hiCD127low, CD4+CD25hiCD127lowFoxP3+, CD4+CD25hiFoxP3+), and with CD4+CD25hi and (2) the changes in Treg and Breg frequencies and their maturation status with age. We performed peripheral blood immunophenotyping of Treg and Breg (CD19+CD24hiCD38hi) by flow cytometry in 55 healthy pediatric controls. We observed that Treg numbers varied depending on the definition used, and the frequency ranged between 3.3-9.7% for CD4+CD25hiCD127low, 0.07-1.6% for CD4+CD25hiCD127lowFoxP3+, and 0.24-2.83% for CD4+CD25hiFoxP3+. The correlation between the three definitions of Treg was positive for most age ranges, especially between the two intracellular panels and with CD4+CD25hi vs CD4+CD25hiCD127low. Treg and Breg frequencies tended to decline after 7 and 3 years onwards, respectively. Treg's maturation status increased with age, with a decline of naïve Treg and an increase in memory/effector Treg from age 7 onwards. Memory Breg increased progressively from age 3 onwards. In conclusion, the number of Treg frequencies spans a wide range depending on the immunophenotypic definition used despite a good level of correlation exists between them. The decline in numbers and maturation process with age occurs earlier in Breg than in Treg.

Novel genetic resources associated with sucrose and stachyose content through genome-wide association study in soybean (Glycine max (L.) Merr.)

The nutritional value of soybean [Glycine max (L.) Merr.] for animals is influenced by soluble carbohydrates, such as sucrose and stachyose. Although sucrose is nutritionally desirable, stachyose is an antinutrient causing diarrhea and flatulence in non-ruminant animals. We conducted a genome-wide association study of 220 soybean accessions using 21,317 single nucleotide polymorphisms (SNPs) from the SoySNP50K iSelect Beadchip data to identify significant SNPs associated with sucrose and stachyose content. Seven significant SNPs were identified for sucrose content across chromosomes (Chrs.) 2, 8, 12, 17, and 20, while thirteen significant SNPs were identified for stachyose content across Chrs. 2, 5, 8, 9, 10, 13, 14, and 15. Among those significant SNPs, three sucrose-related SNPs on Chrs. 8 and 17 were novel, while twelve stachyose-related SNPs on Chrs. 2, 5, 8, 9, 10, 13, 14, and 15 were novel. Based on Phytozome, STRING, and GO annotation, 17 and 24 candidate genes for sucrose and stachyose content, respectively, were highly associated with the carbohydrate metabolic pathway. Among these, the publicly available RNA-seq Atlas database highlighted four candidate genes associated with sucrose (Glyma.08g361200 and Glyma.17g258100) and stachyose (Glyma.05g025300 and Glyma.13g077900) content, which had higher gene expression levels in developing seed and multiple parts of the soybean plant. The results of this study will extend knowledge of the molecular mechanism and genetic basis underlying sucrose and stachyose content in soybean seed. Furthermore, the novel candidate genes and SNPs can be valuable genetic resources that soybean breeders may utilize to modify carbohydrate profiles for animal and human usage.

Variants of IL6, IL10, FCN2, RNASE3, IL12B and IL17B loci are associated with Schistosoma mansoni worm burden in the Albert Nile region of Uganda

BACKGROUND

Individuals genetically susceptible to high schistosomiasis worm burden may contribute disproportionately to transmission and could be prioritized for control. Identifying genes involved may guide development of therapy.

METHODOLOGY/PRINCIPAL FINDINGS

A cohort of 606 children aged 10-15 years were recruited in the Albert Nile region of Uganda and assessed for Schistosoma mansoni worm burden using the Up-Converting Particle Lateral Flow (UCP-LF) test detecting circulating anodic antigen (CAA), point-of-care Circulating Cathodic Antigen (POC-CCA) and Kato-Katz tests. Whole genome genotyping was conducted on 326 children comprising the top and bottom 25% of worm burden. Linear models were fitted to identify variants associated with worm burden in preselected candidate genes. Expression quantitative trait locus (eQTL) analysis was conducted for candidate genes with UCP-LF worm burden included as a covariate. Single Nucleotide Polymorphism loci associated with UCP-LF CAA included IL6 rs2066992 (OR = 0.43, p = 0.0006) and rs7793163 (OR = 2.0, p = 0.0007); IL21 SNP kgp513476 (OR 1.79, p = 0.0025) and IL17B SNP kgp708159 (OR = 0.35, p = 0.0028). A haplotype in the IL10 locus was associated with lower worm burden (OR = 0.53, p = 0.015) and overlapped SNPs rs1800896, rs1800871 and rs1800872. Significant haplotypes (p<0.05, overlapping significant SNP) associated with worm burden were observed in IL6 and the Th17 pathway IL12B and IL17B genes. There were significant eQTL in the IL6, IL5, IL21, IL25 and IFNG regions.

CONCLUSIONS

Variants associated with S. mansoni worm burden were in IL6, FCN2, RNASE3, IL10, IL12B and IL17B gene loci. However only eQTL associations remained significant after Bonferroni correction. In summary, immune balance, pathogen recognition and Th17 pathways may play a role in modulating Schistosoma worm burden. Individuals carrying risk variants may be targeted first in allocation of control efforts to reduce the burden of schistosomiasis in the community.

Genetic Defects in Early-Onset Inflammatory Bowel Disease

Inflammatory bowel disease (IBD) represents a spectrum of disease, which is characterized by chronic gastrointestinal inflammation. Monogenic mutations driving IBD pathogenesis are more highly represented in early-onset compared to adult-onset disease. The pathogenic genes which dysregulate host immune responses in monogenic IBD affect both the innate (ie, intestinal barrier, phagocytes) and adaptive immune systems (ie, T cells, B cells). Advanced genomic and targeted functional testing can improve clinical decision making and present increased opportunities for precision medicine approaches in this important patient population.

Severe kidney dysfunction in sialidosis mice reveals an essential role for neuraminidase 1 in reabsorption

Sialidosis is an ultra-rare multisystemic lysosomal disease caused by mutations in the neuraminidase 1 (NEU1) gene. The severe type II form of the disease manifests with a prenatal/infantile or juvenile onset, bone abnormalities, severe neuropathology, and visceromegaly. A subset of these patients present with nephrosialidosis, characterized by abrupt onset of fulminant glomerular nephropathy. We studied the pathophysiological mechanism of the disease in 2 NEU1-deficient mouse models, a constitutive Neu1-knockout, Neu1ΔEx3, and a conditional phagocyte-specific knockout, Neu1Cx3cr1ΔEx3. Mice of both strains exhibited terminal urinary retention and severe kidney damage with elevated urinary albumin levels, loss of nephrons, renal fibrosis, presence of storage vacuoles, and dysmorphic mitochondria in the intraglomerular and tubular cells. Glycoprotein sialylation in glomeruli, proximal distal tubules, and distal tubules was drastically increased, including that of an endocytic reabsorption receptor megalin. The pool of megalin bearing O-linked glycans with terminal galactose residues, essential for protein targeting and activity, was reduced to below detection levels. Megalin levels were severely reduced, and the protein was directed to lysosomes instead of the apical membrane. Together, our results demonstrated that desialylation by NEU1 plays a crucial role in processing and cellular trafficking of megalin and that NEU1 deficiency in sialidosis impairs megalin-mediated protein reabsorption.

Newborn genetic screening is highly effective for high-risk infants: A single-centre study in China

Background

Newborn genetic screening (NBGS) is promising for early detection of genetic diseases in newborns. However, little is known about its clinical effectiveness in special groups like high-risk infants. To address this gap, we aimed to investigate the impact of NBGS on high-risk infants.

Methods

We screened 10 334 healthy newborns from the general maternity unit and 886 high-risk infants from the neonatal ward using both traditional newborn screening (tNBS) and NBGS, and collected clinical data from electronic medical records.

Results

We found that high-risk infants had a higher proportion of eutocia (P < 0.01) and prematurity (P < 0.01). For high-risk infants vs healthy newborns screened by tNBS, the primary screening positive rate was 3.84% vs 1.31%, the false positive rate (FPR) was 3.62% vs 1.18% (P < 0.001), and the positive predictive value (PPV) was 5.88% vs 8.27%. For NBGS vs tNBS in high-risk infants, the primary screening positive rate was 0.54% vs 3.68%, the FPR was 0.22% vs 3.47%, and the PPV was 60.00% vs 5.88%.

Conclusions

We found that combined newborn screening can effectively reduce the FPR caused by the high-risk symptoms and improve the PPV in high-risk infants, sufficient for more accurately showing the true status of the disease.

Genome Sequencing of Consanguineous Family Implicates Ubiquitin-Specific Protease 53 (USP53) Variant in Psychosis/Schizophrenia: Wild-Type Expression in Murine Hippocampal CA 1-3 and Granular Dentate with AMPA Synapse Interactions

Psychosis is a severe mental disorder characterized by abnormal thoughts and perceptions (e.g., hallucinations) occurring quintessentially in schizophrenia and in several other neuropsychiatric disorders. Schizophrenia is widely considered as a neurodevelopmental disorder that onsets during teenage/early adulthood. A multiplex consanguineous Pakistani family was afflicted with severe psychosis and apparent autosomal recessive transmission. The first-cousin parents and five children were healthy, whereas two teenage daughters were severely affected. Structured interviews confirmed the diagnosis of DSM-V schizophrenia. Probands and father underwent next-generation sequencing. All available relatives were subjected to confirmatory Sanger sequencing. Homozygosity mapping and directed a priori filtering identified only one rare variant [MAF < 5(10)-5] at a residue conserved across vertebrates. The variant was a non-catalytic deubiquitinase, USP53 (p.Cys228Arg), predicted in silico as damaging. Genome sequencing did not identify any other potentially pathogenic single nucleotide variant or structural variant. Since the literature on USP53 lacked relevance to mental illness or CNS expression, studies were conducted which revealed USP53 localization in regions of the hippocampus (CA 1-3) and granular dentate. The staining pattern was like that seen with GRIA2/GluA2 and GRIP2 antibodies. All three proteins coimmunoprecipitated. These findings support the glutamate hypothesis of schizophrenia as part of the AMPA-R interactome. If confirmed, USP53 appears to be one of the few Mendelian variants potentially causal to a common-appearing mental disorder that is a rare genetic form of schizophrenia.

A multilayered post-GWAS analysis pipeline defines functional variants and target genes for systemic lupus erythematosus (SLE)

Objectives

Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWAS) have revealed multiple SLE susceptibility loci and associated single nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown.

Methods

Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci ( cis - and trans -eQTLs) with promoter-capture Hi-C, allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, qPCR, and Western blot.

Results

Anchoring on 452 index SNPs, we selected 9,931 high-linkage disequilibrium (r 2 >0.8) SNPs and defined 182 independent non-HLA SLE loci. 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated rs57668933 as a functional variant regulating multiple targets, including SLE risk gene ELF1 , in B-cells.

Conclusion

We demonstrate and validate post-GWAS strategies for utilizing multi-dimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms, to guide experimental characterization.

The effect of ascertainment on penetrance estimates for rare variants: Implications for establishing pathogenicity and for genetic counselling

Next-generation sequencing has led to an explosion of genetic findings for many rare diseases. However, most of the variants identified are very rare and were also identified in small pedigrees, which creates challenges in terms of penetrance estimation and translation into genetic counselling in the setting of cascade testing. We use simulations to show that for a rare (dominant) disorder where a variant is identified in a small number of small pedigrees, the penetrance estimate can both have large uncertainty and be drastically inflated, due to underlying ascertainment bias. We have developed PenEst, an app that allows users to investigate the phenomenon across ranges of parameter settings. We also illustrate robust ascertainment corrections via the LOD (logarithm of the odds) score, and recommend a LOD-based approach to assessing pathogenicity of rare variants in the presence of reduced penetrance.

Deficiency in the neural cell adhesion molecule 2 (NCAM2) reduces axonal levels of beta-site amyloid precursor protein cleaving enzyme 1 (BACE1), affects axonal organization in the hippocampus, and leads to behavioral deficits

The neural cell adhesion molecule 2 (NCAM2) regulates axonal organization in the central nervous system via mechanisms that have remained poorly understood. We now show that NCAM2 increases axonal levels of beta-site amyloid precursor protein cleaving enzyme 1 (BACE1), a protease that regulates axonal guidance. In brains of NCAM2-deficient mice, BACE1 levels are reduced in hippocampal mossy fiber projections, and the infrapyramidal bundle of these projections is shortened. This abnormal axonal organization correlates with impaired short-term spatial memory and cognitive flexibility in NCAM2-deficient male and female mice. Self-grooming, rearing, digging and olfactory acuity are increased in NCAM2-deficient male mice, when compared with littermate wild-type mice of the same sex. NCAM2-deficient female mice also show increased self-grooming, but are reduced in rearing, and do not differ from female wild-type mice in olfactory acuity and digging behavior. Our results indicate that errors in axonal guidance and organization caused by impaired BACE1 function can underlie the manifestation of neurodevelopmental disorders, including autism as found in humans with deletions of the NCAM2 gene.

Leveraging genetic diversity to understand monogenic Parkinson's disease's landscape in AfrAbia

Parkinson's disease may be caused by a single highly deleterious and penetrant pathogenic variant in 5-10% of cases (monogenic). Research into these mutational disorders yields important pathophysiological insights. This article examines the phenotype, genotype, pathophysiology, and geographic and ethnic distribution of genetic forms of disease. Well established Parkinson's disease (PD) causal variants can follow an autosomal dominant (SNCA, LRRK2, and VPS35) and autosomal recessive pattern of inheritance (PRKN, PINK1, and DJ). Parkinson's disease is a worldwide condition, yet the AfrAbia population is understudied in this regard. No prevalence or incidence investigations have been conducted yet. Few studies on genetic risk factors for PD in AfrAbia communities have been reported which supported the notion that the prevalence and incidence rates of PD in AfrAbia are generally lower than those reported for European and North American populations. There have been only a handful of documented genetic studies of PD in AfrAbia and very limited cohort and case-control research studies on PD have been documented. In this article, we provide a summary of prior conducted research on monogenic PD in Africa and highlight data gaps and promising new research directions. We emphasize that monogenic Parkinson's disease is influenced by distinctions in ethnicity and geography, thereby reinforcing the need for global initiatives to aggregate large numbers of patients and identify novel candidate genes. The current article increases our knowledge of the genetics of Parkinson's disease (PD) and helps to further our knowledge on the genetic factors that contribute to PD, such as the lower penetrance and varying clinical expressivity of known genetic variants, particularly in AfrAbian PD patients.

Toward Identification of Functional Sequences and Variants in Noncoding DNA

Understanding the noncoding part of the genome, which encodes gene regulation, is necessary to identify genetic mechanisms of disease and translate findings from genome-wide association studies into actionable results for treatments and personalized care. Here we provide an overview of the computational analysis of noncoding regions, starting from gene-regulatory mechanisms and their representation in data. Deep learning methods, when applied to these data, highlight important regulatory sequence elements and predict the functional effects of genetic variants. These and other algorithms are used to predict damaging sequence variants. Finally, we introduce rare-variant association tests that incorporate functional annotations and predictions in order to increase interpretability and statistical power.

Recessive CHRM5 variant as a potential cause of neurogenic bladder

Neurogenic bladder is caused by disruption of neuronal pathways regulating bladder relaxation and contraction. In severe cases, neurogenic bladder can lead to vesicoureteral reflux, hydroureter, and chronic kidney disease. These complications overlap with manifestations of congenital anomalies of the kidney and urinary tract (CAKUT). To identify novel monogenic causes of neurogenic bladder, we applied exome sequencing (ES) to our cohort of families with CAKUT. By ES, we have identified a homozygous missense variant (p.Gln184Arg) in CHRM5 (cholinergic receptor, muscarinic, 5) in a patient with neurogenic bladder and secondary complications of CAKUT. CHRM5 codes for a seven transmembrane-spanning G-protein-coupled muscarinic acetylcholine receptor. CHRM5 is shown to be expressed in murine and human bladder walls and is reported to cause bladder overactivity in Chrm5 knockout mice. We investigated CHRM5 as a potential novel candidate gene for neurogenic bladder with secondary complications of CAKUT. CHRM5 is similar to the cholinergic bladder neuron receptor CHRNA3, which Mann et al. published as the first monogenic cause of neurogenic bladder. However, functional in vitro studies did not reveal evidence to strengthen the status as a candidate gene. Discovering additional families with CHRM5 variants could help to further assess the genes' candidate status.

Leveraging AI Advances and Online Tools for Structure-Based Variant Analysis

Understanding how a gene variant affects protein function is important in life science, as it helps explain traits or dysfunctions in organisms. In a clinical setting, this understanding makes it possible to improve and personalize patient care. Bioinformatic tools often only assign a pathogenicity score, rather than providing information about the molecular basis for phenotypes. Experimental testing can furnish this information, but this is slow and costly and requires expertise and equipment not available in a clinical setting. Conversely, mapping a gene variant onto the three-dimensional (3D) protein structure provides a fast molecular assessment free of charge. Before 2021, this type of analysis was severely limited by the availability of experimentally determined 3D protein structures. Advances in artificial intelligence algorithms now allow confident prediction of protein structural features from sequence alone. The aim of the protocols presented here is to enable non-experts to use databases and online tools to investigate the molecular effect of a genetic variant. The Basic Protocol relies only on the online resources AlphaFold, Protein Structure Database, and UniProt. Alternate Protocols document the usage of the Protein Data Bank, SWISS-MODEL, ColabFold, and PyMOL for structure-based variant analysis. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: 3D Mapping based on UniProt and AlphaFold Alternate Protocol 1: Using experimental models from the PDB Alternate Protocol 2: Using information from homology modeling with SWISS-MODEL Alternate Protocol 3: Predicting 3D structures with ColabFold Alternate Protocol 4: Structure visualization and analysis with PyMOL.

A Clinical Qualification Protocol Highlights Overlapping Genomic Influences and Neuro-Autonomic Mechanisms in Ehlers-Danlos and Long COVID-19 Syndromes

A substantial fraction of the 15% with double-jointedness or hypermobility have the traditionally ascertained joint-skeletal, cutaneous, and cardiovascular symptoms of connective tissue dysplasia and its particular manifestation as Ehlers-Danlos syndrome (EDS). The holistic ascertainment of 120 findings in 1261 EDS patients added neuro-autonomic symptoms like headaches, muscle weakness, brain fog, chronic fatigue, dyspnea, and bowel irregularity to those of arthralgia and skin laxity, 15 of these symptoms shared with those of post-infectious SARS-CoV-2 (long COVID-19). Underlying articulo-autonomic mechanisms guided a clinical qualification protocol that qualified DNA variants in 317 genes as having diagnostic utility for EDS, six of them identical (F2-LIFR-NLRP3-STAT1-T1CAM1-TNFRSF13B) and eighteen similar to those modifying COVID-19 severity/EDS, including ADAMTS13/ADAMTS2-C3/C1R-IKBKG/IKBKAP-PIK3C3/PIK3R1-POLD4/POLG-TMPRSS2/TMPRSS6-WNT3/WNT10A. Also, contributing to EDS and COVID-19 severity were forty and three genes, respectively, impacting mitochondrial functions as well as parts of an overlapping gene network, or entome, that are hypothesized to mediate the cognitive-behavioral, neuro-autonomic, and immune-inflammatory alterations of connective tissue in these conditions. The further characterization of long COVID-19 natural history and genetic predisposition will be necessary before these parallels to EDS can be carefully delineated and translated into therapies.

Deep structured learning for variant prioritization in Mendelian diseases

Effective computer-aided or automated variant evaluations for monogenic diseases will expedite clinical diagnostic and research efforts of known and novel disease-causing genes. Here we introduce MAVERICK: a Mendelian Approach to Variant Effect pRedICtion built in Keras. MAVERICK is an ensemble of transformer-based neural networks that can classify a wide range of protein-altering single nucleotide variants (SNVs) and indels and assesses whether a variant would be pathogenic in the context of dominant or recessive inheritance. We demonstrate that MAVERICK outperforms all other major programs that assess pathogenicity in a Mendelian context. In a cohort of 644 previously solved patients with Mendelian diseases, MAVERICK ranks the causative pathogenic variant within the top five variants in over 95% of cases. Seventy-six percent of cases were solved by the top-ranked variant. MAVERICK ranks the causative pathogenic variant in hitherto novel disease genes within the first five candidate variants in 70% of cases. MAVERICK has already facilitated the identification of a novel disease gene causing a degenerative motor neuron disease. These results represent a significant step towards automated identification of causal variants in patients with Mendelian diseases.

Sub-region analysis of DMD gene in cases with idiopathic generalized epilepsy

Gene sub-region encoded protein domain is the basic unit for protein structure and function. The DMD gene is the largest coding gene in humans, with its phenotype relevant to idiopathic generalized epilepsy. We hypothesized variants clustered in sub-regions of idiopathic generalized epilepsy genes and investigated the relationship between the DMD gene and idiopathic generalized epilepsy. Whole exome sequencing was performed in 106 idiopathic generalized epilepsy individuals. DMD variants were filtered with variant type, allele frequency, in silico prediction, hemizygous or homozygous status in the population, inheritance mode, and domain location. Variants located at the sub-regions were selected by the subRVIS software. The pathogenicity of variants was evaluated by the American College of Medical Genetics and Genomics criteria. Articles on functional studies related to epilepsy for variants clustered protein domains were reviewed. In sub-regions of the DMD gene, two variants were identified in two unrelated cases with juvenile absence epilepsy or juvenile myoclonic epilepsy. The pathogenicity of both variants was uncertain significance. Allele frequency of both variants in probands with idiopathic generalized epilepsy reached statistical significance compared with the population (Fisher's test, p = 2.02 × 10-6, adjusted α = 4.52 × 10-6). The variants clustered in the spectrin domain of dystrophin, which binds to glycoprotein complexes and indirectly affects ion channels contributing to epileptogenesis. Gene sub-region analysis suggests a weak association between the DMD gene and idiopathic generalized epilepsy. Functional analysis of gene sub-region helps infer the pathogenesis of idiopathic generalized epilepsy.

Innovative Family-Based Genetically Informed Series of Analyses of Whole-Exome Data Supports Likely Inheritance for Grammar in Children with Specific Language Impairment

Individuals with specific language impairment (SLI) struggle with language acquisition despite average non-verbal intelligence and otherwise typical development. One SLI account focuses on grammar acquisition delay. The current study aimed to detect novel rare genetic variants associated with performance on a grammar assessment, the Test of Early Grammatical Impairment (TEGI), in English-speaking children. The TEGI was selected due to its sensitivity and specificity, consistently high heritability estimates, and its absence from all but one molecular genetic study. We performed whole exome sequencing (WES) in eight families with SLI (n = 74 total) and follow-up Sanger sequencing in additional unrelated probands (n = 146). We prioritized rare exonic variants shared by individuals with low TEGI performance (n = 34) from at least two families under two filtering workflows: (1) novel and (2) previously reported candidate genes. Candidate variants were observed on six new genes (PDHA2, PCDHB3, FURIN, NOL6, IQGAP3, and BAHCC1), and two genes previously reported for overall language ability (GLI3 and FLNB). We specifically suggest PCDHB3, a protocadherin gene, and NOL6 are critical for ribosome synthesis, as they are important targets of SLI investigation. The proposed SLI candidate genes associated with TEGI performance emphasize the utility of precise phenotyping and family-based genetic study.

Invited review: Good practices in genome-wide association studies to identify candidate sequence variants in dairy cattle

Genotype data from dairy cattle selection programs have greatly facilitated GWAS to identify variants related to economic traits. Results can enhance the accuracy of genomic prediction, analyze more complex models that go beyond additive effects, elucidate the genetic architecture of a trait, and finally, decipher the underlying biology of traits. The entire process, comprising data generation, quality control, statistical analyses, interpretation of association results, and linking results to biology should be designed and executed to minimize the generation of false-positive and false-negative associations and misleading links to biological processes. This review aims to provide general guidelines for data analysis that address data quality control, association tests, adjustment for population stratification, and significance evaluation to improve the reliability of conclusions. We also provide guidance on post-GWAS strategy and the interpretation of results. These guidelines are tailored to dairy cattle, which are characterized by long-range linkage disequilibrium, large half-sib families, and routinely collected phenotypes, requiring different approaches than those applied in human GWAS. We discuss common limitations and challenges that have been overlooked in the analysis and interpretation of GWAS to identify candidate sequence variants in dairy cattle.

The landscape of tolerated genetic variation in humans and primates

Personalized genome sequencing has revealed millions of genetic differences between individuals, but our understanding of their clinical relevance remains largely incomplete. To systematically decipher the effects of human genetic variants, we obtained whole-genome sequencing data for 809 individuals from 233 primate species and identified 4.3 million common protein-altering variants with orthologs in humans. We show that these variants can be inferred to have nondeleterious effects in humans based on their presence at high allele frequencies in other primate populations. We use this resource to classify 6% of all possible human protein-altering variants as likely benign and impute the pathogenicity of the remaining 94% of variants with deep learning, achieving state-of-the-art accuracy for diagnosing pathogenic variants in patients with genetic diseases.

Identification and validation of differentially expressed genes for targeted therapy in NSCLC using integrated bioinformatics analysis

Background

Despite the high prevalence of lung cancer, with a five-year survival rate of only 23%, the underlying molecular mechanisms of non-small cell lung cancer (NSCLC) remain unknown. There is a great need to identify reliable candidate biomarker genes for early diagnosis and targeted therapeutic strategies to prevent cancer progression.

Methods

In this study, four datasets obtained from the Gene Expression Omnibus were evaluated for NSCLC- associated differentially expressed genes (DEGs) using bioinformatics analysis. About 10 common significant DEGs were shortlisted based on their p-value and FDR (DOCK4, ID2, SASH1, NPR1, GJA4, TBX2, CD24, HBEGF, GATA3, and DDR1). The expression of significant genes was validated using experimental data obtained from TCGA and the Human Protein Atlas database. The human proteomic data for post- translational modifications was used to interpret the mutations in these genes.

Results

Validation of DEGs revealed a significant difference in the expression of hub genes in normal and tumor tissues. Mutation analysis revealed 22.69%, 48.95%, and 47.21% sequence predicted disordered regions of DOCK4, GJA4, and HBEGF, respectively. The gene-gene and drug-gene network analysis revealed important interactions between genes and chemicals suggesting they could act as probable drug targets. The system-level network showed important interactions between these genes, and the drug interaction network showed that these genes are affected by several types of chemicals that could serve as potential drug targets.

Conclusions

The study demonstrates the importance of systemic genetics in identifying potential drug- targeted therapies for NSCLC. The integrative system- level approach should contribute to a better understanding of disease etiology and may accelerate drug discovery for many cancer types.

Whole Exome-Trio Analysis Reveals Rare Variants Associated with Congenital Pouch Colon

Anorectal malformations (ARM) are individually common, but Congenital Pouch Colon (CPC) is a rare anorectal anomaly that causes a dilated pouch and communication with the genitourinary tract. In this work, we attempted to identify de novo heterozygous missense variants, and further discovered variants of unknown significance (VUS) which could provide insights into CPC manifestation. From whole exome sequencing (WES) performed earlier, the trio exomes were analyzed from those who were admitted to J.K. Lon Hospital, SMS Medical College, Jaipur, India, between 2011 and 2017. The proband exomes were compared with the unaffected sibling/family members, and we sought to ask whether any variants of significant interest were associated with the CPC manifestation. The WES data from a total of 64 samples including 16 affected neonates (11 male and 5 female) with their parents and unaffected siblings were used for the study. We examined the role of rare allelic variation associated with CPC in a 16 proband/parent trio family, comparing the mutations to those of their unaffected parents/siblings. We also performed RNA-Seq as a pilot to find whether or not the genes harboring these mutations were differentially expressed. Our study revealed extremely rare variants, viz., TAF1B, MUC5B and FRG1, which were further validated for disease-causing mutations associated with CPC, further closing the gaps of surgery by bringing intervention in therapies.

Genetic Influences on the Developing Young Brain and Risk for Neuropsychiatric Disorders

Imaging genetics provides an opportunity to discern associations between genetic variants and brain imaging phenotypes. Historically, the field has focused on adults and adolescents; very few imaging genetics studies have focused on brain development in infancy and early childhood (from birth to age 6 years). This is an important knowledge gap because developmental changes in the brain during the prenatal and early postnatal period are regulated by dynamic gene expression patterns that likely play an important role in establishing an individual's risk for later psychiatric illness and neurodevelopmental disabilities. In this review, we summarize findings from imaging genetics studies spanning from early infancy to early childhood, with a focus on studies examining genetic risk for neuropsychiatric disorders. We also introduce the Organization for Imaging Genomics in Infancy (ORIGINs), a working group of the ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) consortium, which was established to facilitate large-scale imaging genetics studies in infancy and early childhood.

The landscape of tolerated genetic variation in humans and primates

Personalized genome sequencing has revealed millions of genetic differences between individuals, but our understanding of their clinical relevance remains largely incomplete. To systematically decipher the effects of human genetic variants, we obtained whole genome sequencing data for 809 individuals from 233 primate species, and identified 4.3 million common protein-altering variants with orthologs in human. We show that these variants can be inferred to have non-deleterious effects in human based on their presence at high allele frequencies in other primate populations. We use this resource to classify 6% of all possible human protein-altering variants as likely benign and impute the pathogenicity of the remaining 94% of variants with deep learning, achieving state-of-the-art accuracy for diagnosing pathogenic variants in patients with genetic diseases.

One Sentence Summary

Deep learning classifier trained on 4.3 million common primate missense variants predicts variant pathogenicity in humans.

Genetic and Molecular Contributors in Eosinophilic Esophagitis

PURPOSE OF REVIEW

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic and molecular etiology. The interest of the scientific community in EoE has grown considerably over the past three decades, and the understanding of the genetic and molecular mechanisms involved in this disease has greatly increased.

RECENT FINDINGS

In this article, we aim to provide both historic aspects and updates on the recent genetic and molecular advances in the understanding of EoE. Although EoE is a relatively newly described disorder, much progress has been made toward identifying the genetic and molecular factors contributing to the disease pathogenesis by a variety of approaches with next-generation sequencing technologies, including genome-wide association study, whole exome sequencing, and bulk and single-cell RNA sequencing. This review highlights the multifaceted impacts of various findings that have shaped the current molecular and genetic landscape of EoE, providing insights that facilitate further understanding of the disease process.

Identifying Novel Risk Conferring Genes Involved in Glycosylation Processes with Familial Schizophrenia in an Indian Cohort: Prediction of ADAMTS9 gene Variant for Structural Stability

Schizophrenia is a complex neuropsychiatric disorder and heritability is as high as 80% making it the most heritable mental disorder. Although GWAS has identified numerous variants, the pathophysiology is still elusive. Here, an attempt was made to identify genetic risk factors in familial cases of schizophrenia that are associated with a common causative pathway. To achieve this objective, exome sequencing was done in 4 familial cases and identified six unique coding variants in five genes. Among these genes, PIGQ gene has two pathogenic variants, one nonsense and in-frame deletion. One missense variant in GALNT16 and one in GALNT5 have variable damaging score, however, the other variants, in ADAMTS9 and in LTBP4 have the highest damaging score. Further analysis showed that the variant of LTBP4 was not present in the functional domain. The other missense variant in the ADAMTS9 gene was found to be significant and was present in the thrombospondin repeat motif, one of the important motifs. Detailed molecular dynamics simulation study on this variant showed a damaging effect on structural stability. Since, all these genes culminated into the glycosylation process, it was evident that an aberrant glycosylation process may be one of the risk factors. Although, extracellular matrix formation through glycosylation have been shown to be associated, the involvement of ADAMTS9 and PIGQ gene mediated glycosylation has not been reported. In this paper, a novel link between ADAMTS9 and PIGQ gene with schizophrenia have been reported. Therefore, this novel observation has contributed immensely to the existing knowledge on risk factor of Schizophrenia.

Retinitis Pigmentosa: Current Clinical Management and Emerging Therapies

Retinitis pigmentosa (RP) comprises a group of inherited retinal dystrophies characterized by the degeneration of rod photoreceptors, followed by the degeneration of cone photoreceptors. As a result of photoreceptor degeneration, affected individuals experience gradual loss of visual function, with primary symptoms of progressive nyctalopia, constricted visual fields and, ultimately, central vision loss. The onset, severity and clinical course of RP shows great variability and unpredictability, with most patients already experiencing some degree of visual disability in childhood. While RP is currently untreatable for the majority of patients, significant efforts have been made in the development of genetic therapies, which offer new hope for treatment for patients affected by inherited retinal dystrophies. In this exciting era of emerging gene therapies, it remains imperative to continue supporting patients with RP using all available options to manage their condition. Patients with RP experience a wide variety of physical, mental and social-emotional difficulties during their lifetime, of which some require timely intervention. This review aims to familiarize readers with clinical management options that are currently available for patients with RP.

Implication of FOXD2 dysfunction in syndromic congenital anomalies of the kidney and urinary tract (CAKUT)

Background

Congenital anomalies of the kidney and urinary tract (CAKUT) are the predominant cause for chronic kidney disease below 30 years of age. Many monogenic forms have been discovered mainly due to comprehensive genetic testing like exome sequencing (ES). However, disease-causing variants in known disease-associated genes still only explain a proportion of cases. Aim of this study was to unravel the underlying molecular mechanism of syndromic CAKUT in two multiplex families with presumed autosomal recessive inheritance.

Methods and Results

ES in the index individuals revealed two different rare homozygous variants in FOXD2, a transcription factor not previously implicated in CAKUT in humans: a frameshift in family 1 and a missense variant in family 2 with family segregation patterns consistent with autosomal-recessive inheritance. CRISPR/Cas9-derived Foxd2 knock-out (KO) mice presented with bilateral dilated renal pelvis accompanied by renal papilla atrophy while extrarenal features included mandibular, ophthalmologic, and behavioral anomalies, recapitulating the phenotype of humans with FOXD2 dysfunction. To study the pathomechanism of FOXD2-dysfunction-mediated developmental renal defects, in a complementary approach, we generated CRISPR/Cas9-mediated KO of Foxd2 in ureteric-bud-induced mouse metanephric mesenchyme cells. Transcriptomic analyses revealed enrichment of numerous differentially expressed genes important in renal/urogenital development, including Pax2 and Wnt4 as well as gene expression changes indicating a cell identity shift towards a stromal cell identity. Histology of Foxd2 KO mouse kidneys confirmed increased fibrosis. Further, GWAS data (genome-wide association studies) suggests that FOXD2 could play a role for maintenance of podocyte integrity during adulthood.

Conclusions

In summary, our data implicate that FOXD2 dysfunction is a very rare cause of autosomal recessive syndromic CAKUT and suggest disturbances of the PAX2-WNT4 cell signaling axis contribute to this phenotype.

Hacking hematopoiesis - emerging tools for examining variant effects

Hematopoiesis is a continuous process of blood and immune cell production. It is orchestrated by thousands of gene products that respond to extracellular signals by guiding cell fate decisions to meet the needs of the organism. Although much of our knowledge of this process comes from work in model systems, we have learned a great deal from studies on human genetic variation. Considerable insight has emerged from studies on presumed monogenic blood disorders, which continue to provide key insights into the mechanisms critical for hematopoiesis. Furthermore, the emergence of large-scale biobanks and cohorts has uncovered thousands of genomic loci associated with blood cell traits and diseases. Some of these blood cell trait-associated loci act as modifiers of what were once thought to be monogenic blood diseases. However, most of these loci await functional validation. Here, we discuss the validation bottleneck and emerging methods to more effectively connect variant to function. In particular, we highlight recent innovations in genome editing, which have paved the path forward for high-throughput functional assessment of loci. Finally, we discuss existing barriers to progress, including challenges in manipulating the genomes of primary hematopoietic cells.

Hypothesis-free phenotype prediction within a genetics-first framework

Cohort-wide sequencing studies have revealed that the largest category of variants is those deemed 'rare', even for the subset located in coding regions (99% of known coding variants are seen in less than 1% of the population. Associative methods give some understanding how rare genetic variants influence disease and organism-level phenotypes. But here we show that additional discoveries can be made through a knowledge-based approach using protein domains and ontologies (function and phenotype) that considers all coding variants regardless of allele frequency. We describe an ab initio, genetics-first method making molecular knowledge-based interpretations for exome-wide non-synonymous variants for phenotypes at the organism and cellular level. By using this reverse approach, we identify plausible genetic causes for developmental disorders that have eluded other established methods and present molecular hypotheses for the causal genetics of 40 phenotypes generated from a direct-to-consumer genotype cohort. This system offers a chance to extract further discovery from genetic data after standard tools have been applied.

Genetic Variants in ARHGEF6 Cause Congenital Anomalies of the Kidneys and Urinary Tract in Humans, Mice, and Frogs

BACKGROUND

About 40 disease genes have been described to date for isolated CAKUT, the most common cause of childhood CKD. However, these genes account for only 20% of cases. ARHGEF6, a guanine nucleotide exchange factor that is implicated in biologic processes such as cell migration and focal adhesion, acts downstream of integrin-linked kinase (ILK) and parvin proteins. A genetic variant of ILK that causes murine renal agenesis abrogates the interaction of ILK with a murine focal adhesion protein encoded by Parva , leading to CAKUT in mice with this variant.

METHODS

To identify novel genes that, when mutated, result in CAKUT, we performed exome sequencing in an international cohort of 1265 families with CAKUT. We also assessed the effects in vitro of wild-type and mutant ARHGEF6 proteins, and the effects of Arhgef6 deficiency in mouse and frog models.

RESULTS

We detected six different hemizygous variants in the gene ARHGEF6 (which is located on the X chromosome in humans) in eight individuals from six families with CAKUT. In kidney cells, overexpression of wild-type ARHGEF6 -but not proband-derived mutant ARHGEF6 -increased active levels of CDC42/RAC1, induced lamellipodia formation, and stimulated PARVA-dependent cell spreading. ARHGEF6-mutant proteins showed loss of interaction with PARVA. Three-dimensional Madin-Darby canine kidney cell cultures expressing ARHGEF6-mutant proteins exhibited reduced lumen formation and polarity defects. Arhgef6 deficiency in mouse and frog models recapitulated features of human CAKUT.

CONCLUSIONS

Deleterious variants in ARHGEF6 may cause dysregulation of integrin-parvin-RAC1/CDC42 signaling, thereby leading to X-linked CAKUT.

Functional characterization of GATA6 genetic variants associated with mild congenital heart defects

Damaging GATA6 variants can cause moderate congenital heart defects. With the application of next-generation sequencing approaches, various novel GATA6 variants with unknown significance have been identified from a broad spectrum of congenital heart defects. However, functional assessment for distinct GATA6 variants from different severity of congenital heart defects, especially from mild defects, is lacking, which hinders our understanding of the genotype-phenotype correlations and underlying mechanisms. Here, we assessed the functional consequences of nine rare GATA6 variants, which had been implicated as the most significant variants associated with mild congenital heart defects using the largest case and control cohort. We examined the effects of these variants on subcellular localization, transcriptional activity, and protein interactions in 293T or AC16 cells and their ability to rescue heart malformation in gata6 zebrafish mutant. We found that two of these nine variants, Q120X and S424I, significantly decreased transcriptional activity. Additionally, Q120X altered subcellular localization. Consistent with the in vitro results, the in vivo results showed that Q120X and S424I lost their potency to rescue ventricular malformation in gata6 -/- embryos. The results indicated that Q120X and S424I are pathogenic in mild congenital heart defects. Further, the inconsistence of severely impaired Q120X function and mild CHDs phenotype suggested the complexity of the genotype-phenotype correlation between the GATA6 variant and heart phenotype, which may help to inform prenatal genetic counseling and pre-implantation genotyping for congenital heart defects.

Informing variant assessment using structured evidence from prior classifications (PS1, PM5, and PVS1 sequence variant interpretation criteria)

PURPOSE

This study aimed to explore whether evidence of pathogenicity from prior variant classifications in ClinVar could be used to inform variant interpretation using the American College of Medical Genetics and Genomics/Association for Molecular Pathology clinical guidelines.

METHODS

We identified distinct single-nucleotide variants (SNVs) that are either similar in location or in functional consequence to pathogenic variants in ClinVar and analyzed evidence in support of pathogenicity using 3 interpretation criteria.

RESULTS

Thousands of variants, including many in clinically actionable disease genes (American College of Medical Genetics and Genomics secondary findings v3.0), have evidence of pathogenicity from existing variant classifications, accounting for 2.5% of nonsynonymous SNVs within ClinVar. Notably, there are many variants with uncertain or conflicting classifications that cause the same amino acid substitution as other pathogenic variants (PS1, N = 323), variants that are predicted to cause different amino acid substitutions in the same codon as pathogenic variants (PM5, N = 7692), and loss-of-function variants that are present in genes in which many loss-of-function variants are classified as pathogenic (PVS1, N = 3635). Most of these variants have similar computational predictions of pathogenicity and splicing effect as their associated pathogenic variants.

CONCLUSION

Broadly, for >1.4 million SNVs exome wide, information from previously classified variants could be used to provide evidence of pathogenicity. We have developed a pipeline to identify variants meeting these criteria that may inform interpretation efforts.

Genetics of circadian rhythms and sleep in human health and disease

Circadian rhythms and sleep are fundamental biological processes integral to human health. Their disruption is associated with detrimental physiological consequences, including cognitive, metabolic, cardiovascular and immunological dysfunctions. Yet many of the molecular underpinnings of sleep regulation in health and disease have remained elusive. Given the moderate heritability of circadian and sleep traits, genetics offers an opportunity that complements insights from model organism studies to advance our fundamental molecular understanding of human circadian and sleep physiology and linked chronic disease biology. Here, we review recent discoveries of the genetics of circadian and sleep physiology and disorders with a focus on those that reveal causal contributions to complex diseases.