A systematic survey of loss-of-function variants in human protein-coding genes

A single base pair duplication in the SLC33A1 gene is associated with fetal losses and neonatal lethality in Manech Tête Rousse dairy sheep

We recently discovered that the Manech Tête Rousse (MTR) deficient homozygous haplotype 2 (MTRDHH2) probably carries a recessive lethal mutation in sheep. In this study, we fine-mapped this region through whole-genome sequencing of five MTRDHH2 heterozygous carriers and 95 non-carriers from various ovine breeds. We identified a single base pair duplication within the SLC33A1 gene, leading to a frameshift mutation and a premature stop codon (p.Arg246Alafs*3). SLC33A1 encodes a transmembrane transporter of acetyl-coenzyme A that is crucial for cellular metabolism. To investigate the lethality of this mutation in homozygous MTR sheep, we performed at-risk matings using artificial insemination (AI) between heterozygous SLC33A1 variant carriers (SLC33A1_dupG). Pregnancy was confirmed 15 days post-AI using a blood test measuring interferon Tau-stimulated MX1 gene expression. Ultrasonography between 45 and 60 days post-AI revealed a 12% reduction in AI success compared with safe matings, indicating embryonic/fetal loss. This was supported by the MX1 differential expression test suggesting fetal losses between 15 and 60 days of gestation. We also observed a 34.7% pre-weaning mortality rate in 49 lambs born from at-risk matings. Homozygous SLC33A1_dupG lambs accounted for 47% of this mortality, with deaths occurring mostly within the first 5 days without visible clinical signs. Therefore, appropriate management of SLC33A1_dupG with an allele frequency of 0.04 in the MTR selection scheme would help increase overall fertility and lamb survival.

Identification of nonsense-mediated decay inhibitors that alter the tumor immune landscape

Despite exciting developments in cancer immunotherapy, its broad application is limited by the paucity of targetable antigens on the tumor cell surface. As an intrinsic cellular pathway, nonsense-mediated decay (NMD) conceals neoantigens through the destruction of the RNA products from genes harboring truncating mutations. We developed and conducted a high throughput screen, based on the ratiometric analysis of transcripts, to identify critical mediators of NMD. This screen implicated disruption of kinase SMG1's phosphorylation of UPF1 as a potential disruptor of NMD. This led us to design a novel SMG1 inhibitor, KVS0001, that elevates the expression of transcripts and proteins resulting from truncating mutations in vivo and in vitro . Most importantly, KVS0001 concomitantly increased the presentation of immune-targetable HLA class I-associated peptides from NMD-downregulated proteins on the surface of cancer cells. KVS0001 provides new opportunities for studying NMD and the diseases in which NMD plays a role, including cancer and inherited diseases.

One Sentence Summary

Disruption of the nonsense-mediated decay pathway with a newly developed SMG1 inhibitor with in-vivo activity increases the expression of T-cell targetable cancer neoantigens resulting from truncating mutations.

Genetic/epigenetic effects in NF1 microdeletion syndrome: beyond the haploinsufficiency, looking at the contribution of not deleted genes

NF1 microdeletion syndrome, accounting for 5-11% of NF1 patients, is caused by a deletion in the NF1 region and it is generally characterized by a severe phenotype. Although 70% of NF1 microdeletion patients presents the same 1.4 Mb type-I deletion, some patients may show additional clinical features. Therefore, the contribution of several pathogenic mechanisms, besides haploinsufficiency of some genes within the deletion interval, is expected and needs to be defined. We investigated an altered expression of deletion flanking genes by qPCR in patients with type-1 NF1 deletion, compared to healthy donors, possibly contributing to the clinical traits of NF1 microdeletion syndrome. In addition, the 1.4-Mb deletion leads to changes in the 3D chromatin structure in the 17q11.2 region. Specifically, this deletion alters DNA-DNA interactions in the regions flanking the breakpoints, as demonstrated by our 4C-seq analysis. This alteration likely causes position effect on the expression of deletion flanking genes.Interestingly, 4C-seq analysis revealed that in microdeletion patients, an interaction was established between the RHOT1 promoter and the SLC6A4 gene, which showed increased expression. We performed NGS on putative modifier genes, and identified two "likely pathogenic" rare variants in RAS pathway, possibly contributing to incidental phenotypic features.This study provides new insights into understanding the pathogenesis of NF1 microdeletion syndrome and suggests a novel pathomechanism that contributes to the expression phenotype in addition to haploinsufficiency of genes located within the deletion.This is a pivotal approach that can be applied to unravel microdeletion syndromes, improving precision medicine, prognosis and patients' follow-up.

Nonlethal deleterious mutation-induced stress accelerates bacterial aging

Random mutagenesis, including when it leads to loss of gene function, is a key mechanism enabling microorganisms' long-term adaptation to new environments. However, loss-of-function mutations are often deleterious, triggering, in turn, cellular stress and complex homeostatic stress responses, called "allostasis," to promote cell survival. Here, we characterize the differential impacts of 65 nonlethal, deleterious single-gene deletions on Escherichia coli growth in three different growth environments. Further assessments of select mutants, namely, those bearing single adenosine triphosphate (ATP) synthase subunit deletions, reveal that mutants display reorganized transcriptome profiles that reflect both the environment and the specific gene deletion. We also find that ATP synthase α-subunit deleted (ΔatpA) cells exhibit elevated metabolic rates while having slower growth compared to wild-type (wt) E. coli cells. At the single-cell level, compared to wt cells, individual ΔatpA cells display near normal proliferation profiles but enter a postreplicative state earlier and exhibit a distinct senescence phenotype. These results highlight the complex interplay between genomic diversity, adaptation, and stress response and uncover an "aging cost" to individual bacterial cells for maintaining population-level resilience to environmental and genetic stress; they also suggest potential bacteriostatic antibiotic targets and -as select human genetic diseases display highly similar phenotypes, - a bacterial origin of some human diseases.

Clustering of predicted loss-of-function variants in genes linked with monogenic disease can explain incomplete penetrance

BACKGROUND

Genetic variants that severely alter protein products (e.g. nonsense, frameshift) are often associated with disease. For some genes, these predicted loss-of-function variants (pLoFs) are observed throughout the gene, whilst in others, they occur only at specific locations. We hypothesised that, for genes linked with monogenic diseases that display incomplete penetrance, pLoF variants present in apparently unaffected individuals may be limited to regions where pLoFs are tolerated. To test this, we investigated whether pLoF location could explain instances of incomplete penetrance of variants expected to be pathogenic for Mendelian conditions.

METHODS

We used exome sequence data in 454,773 individuals in the UK Biobank (UKB) to investigate the locations of pLoFs in a population cohort. We counted numbers of unique pLoF, missense, and synonymous variants in UKB in each quintile of the coding sequence (CDS) of all protein-coding genes and clustered the variants using Gaussian mixture models. We limited the analyses to genes with ≥ 5 variants of each type (16,473 genes). We compared the locations of pLoFs in UKB with all theoretically possible pLoFs in a transcript, and pathogenic pLoFs from ClinVar, and performed simulations to estimate the false-positive rate of non-uniformly distributed variants.

RESULTS

For most genes, all variant classes fell into clusters representing broadly uniform variant distributions, but genes in which haploinsufficiency causes developmental disorders were less likely to have uniform pLoF distribution than other genes (P < 2.2 × 10-6). We identified a number of genes, including ARID1B and GATA6, where pLoF variants in the first quarter of the CDS were rescued by the presence of an alternative translation start site and should not be reported as pathogenic. For other genes, such as ODC1, pLoFs were located approximately uniformly across the gene, but pathogenic pLoFs were clustered only at the end, consistent with a gain-of-function disease mechanism.

CONCLUSIONS

Our results suggest the potential benefits of localised constraint metrics and that the location of pLoF variants should be considered when interpreting variants.

Systematic analysis of nonsense variants uncovers peptide release rate as a novel modifier of nonsense-mediated mRNA decay efficiency

Nonsense variants underlie many genetic diseases. The phenotypic impact of nonsense variants is determined by Nonsense-mediated mRNA decay (NMD), which degrades transcripts with premature termination codons (PTCs). NMD activity varies across transcripts and cellular contexts via poorly understood mechanisms. Here, by leveraging human genetic datasets, we uncover that the amino acid preceding the PTC dramatically affects NMD activity in human cells. We find that glycine codons in particular support high levels of NMD and are enriched before PTCs but depleted before normal termination codons (NTCs). Gly-PTC enrichment is most pronounced in human genes that tolerate loss-of-function variants. This suggests a strong biological impact for Gly-PTC in ensuring robust elimination of potentially toxic truncated proteins from non-essential genes. Biochemical assays revealed that the peptide release rate during translation termination is highly dependent on the identity of the amino acid preceding the stop codon. This release rate is the most critical feature determining NMD activity across our massively parallel reporter assays. Together, we conclude that NMD activity is significantly modulated by the "window of opportunity" offered by translation termination kinetics. Integrating the window of opportunity model with the existing framework of NMD would enable more accurate nonsense variant interpretation in the clinic.

Putative protective genomic variation in the Lithuanian population

Genomic effect variants associated with survival and protection against complex diseases vary between populations due to microevolutionary processes. The aim of this study was to analyse diversity and distribution of effect variants in a context of potential positive selection. In total, 475 individuals of Lithuanian origin were genotyped using high-throughput scanning and/or sequencing technologies. Allele frequency analysis for the pre-selected effect variants was performed using the catalogue of single nucleotide polymorphisms. Comparison of the pre-selected effect variants with variants in primate species was carried out to ascertain which allele was derived and potentially of protective nature. Recent positive selection analysis was performed to verify this protective effect. Four variants having significantly different frequencies compared to European populations were identified while two other variants reached borderline significance. Effect variant in SLC30A8 gene may potentially protect against type 2 diabetes. The existing paradox of high rates of type 2 diabetes in the Lithuanian population and the relatively high frequencies of potentially protective genome variants against it indicate a lack of knowledge about the interactions between environmental factors, regulatory regions, and other genome variation. Identification of effect variants is a step towards better understanding of the microevolutionary processes, etiopathogenetic mechanisms, and personalised medicine.

A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites

Protein evolution is constrained by structure and function, creating patterns in residue conservation that are routinely exploited to predict structure and other features. Similar constraints should affect variation across individuals, but it is only with the growth of human population sequencing that this has been tested at scale. Now, human population constraint has established applications in pathogenicity prediction, but it has not yet been explored for structural inference. Here, we map 2.4 million population variants to 5885 protein families and quantify residue-level constraint with a new Missense Enrichment Score (MES). Analysis of 61,214 structures from the PDB spanning 3661 families shows that missense depleted sites are enriched in buried residues or those involved in small-molecule or protein binding. MES is complementary to evolutionary conservation and a combined analysis allows a new classification of residues according to a conservation plane. This approach finds functional residues that are evolutionarily diverse, which can be related to specificity, as well as family-wide conserved sites that are critical for folding or function. We also find a possible contrast between lethal and non-lethal pathogenic sites, and a surprising clinical variant hot spot at a subset of missense enriched positions.

Predicted genetic burden and frequency of phenotype-associated variants in the horse

Disease-causing variants have been identified for less than 20% of suspected equine genetic diseases. Whole genome sequencing (WGS) allows rapid identification of rare disease causal variants. However, interpreting the clinical variant consequence is confounded by the number of predicted deleterious variants that healthy individuals carry (predicted genetic burden). Estimation of the predicted genetic burden and baseline frequencies of known deleterious or phenotype associated variants within and across the major horse breeds have not been performed. We used WGS of 605 horses across 48 breeds to identify 32,818,945 variants, demonstrate a high predicted genetic burden (median 730 variants/horse, interquartile range: 613-829), show breed differences in predicted genetic burden across 12 target breeds, and estimate the high frequencies of some previously reported disease variants. This large-scale variant catalog for a major and highly athletic domestic animal species will enhance its ability to serve as a model for human phenotypes and improves our ability to discover the bases for important equine phenotypes.

Unravelling the Genetic Landscape of Hemiplegic Migraine: Exploring Innovative Strategies and Emerging Approaches

Migraine is a severe, debilitating neurovascular disorder. Hemiplegic migraine (HM) is a rare and debilitating neurological condition with a strong genetic basis. Sequencing technologies have improved the diagnosis and our understanding of the molecular pathophysiology of HM. Linkage analysis and sequencing studies in HM families have identified pathogenic variants in ion channels and related genes, including CACNA1A, ATP1A2, and SCN1A, that cause HM. However, approximately 75% of HM patients are negative for these mutations, indicating there are other genes involved in disease causation. In this review, we explored our current understanding of the genetics of HM. The evidence presented herein summarises the current knowledge of the genetics of HM, which can be expanded further to explain the remaining heritability of this debilitating condition. Innovative bioinformatics and computational strategies to cover the entire genetic spectrum of HM are also discussed in this review.

Dysregulation of mTOR signaling mediates common neurite and migration defects in both idiopathic and 16p11.2 deletion autism neural precursor cells

Autism spectrum disorder (ASD) is defined by common behavioral characteristics, raising the possibility of shared pathogenic mechanisms. Yet, vast clinical and etiological heterogeneity suggests personalized phenotypes. Surprisingly, our iPSC studies find that six individuals from two distinct ASD subtypes, idiopathic and 16p11.2 deletion, have common reductions in neural precursor cell (NPC) neurite outgrowth and migration even though whole genome sequencing demonstrates no genetic overlap between the datasets. To identify signaling differences that may contribute to these developmental defects, an unbiased phospho-(p)-proteome screen was performed. Surprisingly despite the genetic heterogeneity, hundreds of shared p-peptides were identified between autism subtypes including the mTOR pathway. mTOR signaling alterations were confirmed in all NPCs across both ASD subtypes, and mTOR modulation rescued ASD phenotypes and reproduced autism NPC-associated phenotypes in control NPCs. Thus, our studies demonstrate that genetically distinct ASD subtypes have common defects in neurite outgrowth and migration which are driven by the shared pathogenic mechanism of mTOR signaling dysregulation.

Emerging significance of butyrylcholinesterase

Butyrylcholinesterase (BChE; EC 3.1.1.8), an enzyme structurally related to acetylcholinesterase, is widely distributed in the human body. It plays a role in the detoxification of chemicals such as succinylcholine, a muscle relaxant used in anesthetic practice. BChE is well-known due to variant forms of the enzyme with little or no hydrolytic activity which exist in some endogamous communities and result in prolonged apnea following the administration of succinylcholine. Its other functions include the ability to hydrolyze acetylcholine, the cholinergic neurotransmitter in the brain, when its primary hydrolytic enzyme, acetylcholinesterase, is absent. To assess its potential roles, BChE was studied in relation to insulin resistance, type 2 diabetes mellitus, cognition, hepatic disorders, cardiovascular and cerebrovascular diseases, and inflammatory conditions. Individuals who lack the enzyme activity of BChE are otherwise healthy, until they are given drugs hydrolyzed by this enzyme. Therefore, BChE is a candidate for the study of loss-of-function mutations in humans. Studying individuals with variant forms of BChE can provide insights into whether they are protected against metabolic diseases. The potential utility of the enzyme as a biomarker for Alzheimer's disease and the response to its drug treatment can also be assessed.

Nanopore sequencing of 1000 Genomes Project samples to build a comprehensive catalog of human genetic variation

Less than half of individuals with a suspected Mendelian condition receive a precise molecular diagnosis after comprehensive clinical genetic testing. Improvements in data quality and costs have heightened interest in using long-read sequencing (LRS) to streamline clinical genomic testing, but the absence of control datasets for variant filtering and prioritization has made tertiary analysis of LRS data challenging. To address this, the 1000 Genomes Project ONT Sequencing Consortium aims to generate LRS data from at least 800 of the 1000 Genomes Project samples. Our goal is to use LRS to identify a broader spectrum of variation so we may improve our understanding of normal patterns of human variation. Here, we present data from analysis of the first 100 samples, representing all 5 superpopulations and 19 subpopulations. These samples, sequenced to an average depth of coverage of 37x and sequence read N50 of 54 kbp, have high concordance with previous studies for identifying single nucleotide and indel variants outside of homopolymer regions. Using multiple structural variant (SV) callers, we identify an average of 24,543 high-confidence SVs per genome, including shared and private SVs likely to disrupt gene function as well as pathogenic expansions within disease-associated repeats that were not detected using short reads. Evaluation of methylation signatures revealed expected patterns at known imprinted loci, samples with skewed X-inactivation patterns, and novel differentially methylated regions. All raw sequencing data, processed data, and summary statistics are publicly available, providing a valuable resource for the clinical genetics community to discover pathogenic SVs.

Inherited non-FGF23-mediated phosphaturic disorders: A kidney-centric review

Phosphate is freely filtered by the glomerulus and reabsorbed exclusively in the proximal tubule by two key transporters, NaPiIIA and NaPiIIC, encoded by SLC34A1 and SLC34A3, respectively. Regulation of these transporters occurs primarily through the hormone FGF23 and, to a lesser degree, PTH. Consequently, inherited non-FGF23 mediated phosphaturic disorders are due to generalised proximal tubular dysfunction, loss-of-function variants in SLC34A1 or SLC34A3 or excess PTH signalling. The corresponding disorders are Renal Fanconi Syndrome, Infantile Hypercalcaemia type 2, Hereditary Hypophosphataemic Rickets with Hypercalciuria and Familial Hyperparathyroidism. Several inherited forms of Fanconi renotubular syndrome (FRTS) have also been described with the underlying genes encoding for GATM, EHHADH, HNF4A and NDUFAF6. Here, we will review their pathophysiology, clinical manifestations and the implications for treatment from a kidney-centric perspective, focussing on those disorders caused by dysfunction of renal phosphate transporters. Moreover, we will highlight specific genetic aspects, as the availability of large population genetic databases has raised doubts about some of the originally proposed gene-disease associations concerning phosphate transporters or their associated proteins.

Prioritising genetic findings for drug target identification and validation

The decreasing costs of high-throughput genetic sequencing and increasing abundance of sequenced genome data have paved the way for the use of genetic data in identifying and validating potential drug targets. However, the number of identified potential drug targets is often prohibitively large to experimentally evaluate in wet lab experiments, highlighting the need for systematic approaches for target prioritisation. In this review, we discuss principles of genetically guided drug development, specifically addressing loss-of-function analysis, colocalization and Mendelian randomisation (MR), and the contexts in which each may be most suitable. We subsequently present a range of biomedical resources which can be used to annotate and prioritise disease-associated proteins identified by these studies including 1) ontologies to map genes, proteins, and disease, 2) resources for determining the druggability of a potential target, 3) tissue and cell expression of the gene encoding the potential target, and 4) key biological pathways involving the potential target. We illustrate these concepts through a worked example, identifying a prioritised set of plasma proteins associated with non-alcoholic fatty liver disease (NAFLD). We identified five proteins with strong genetic support for involvement with NAFLD: CYB5A, NT5C, NCAN, TGFBI and DAPK2. All of the identified proteins were expressed in both liver and adipose tissues, with TGFBI and DAPK2 being potentially druggable. In conclusion, the current review provides an overview of genetic evidence for drug target identification, and how biomedical databases can be used to provide actionable prioritisation, fully informing downstream experimental validation.

50,000 years of Evolutionary History of India: Insights from ~2,700 Whole Genome Sequences

India has been underrepresented in whole genome sequencing studies. We generated 2,762 high coverage genomes from India-including individuals from most geographic regions, speakers of all major languages, and tribal and caste groups-providing a comprehensive survey of genetic variation in India. With these data, we reconstruct the evolutionary history of India through space and time at fine scales. We show that most Indians derive ancestry from three ancestral groups related to ancient Iranian farmers, Eurasian Steppe pastoralists and South Asian hunter-gatherers. We uncover a common source of Iranian-related ancestry from early Neolithic cultures of Central Asia into the ancestors of Ancestral South Indians (ASI), Ancestral North Indians (ANI), Austro-asiatic-related and East Asian-related groups in India. Following these admixtures, India experienced a major demographic shift towards endogamy, resulting in extensive homozygosity and identity-by-descent sharing among individuals. At deep time scales, Indians derive around 1-2% of their ancestry from gene flow from archaic hominins, Neanderthals and Denisovans. By assembling the surviving fragments of archaic ancestry in modern Indians, we recover ~1.5 Gb (or 50%) of the introgressing Neanderthal and ~0.6 Gb (or 20%) of the introgressing Denisovan genomes, more than any other previous archaic ancestry study. Moreover, Indians have the largest variation in Neanderthal ancestry, as well as the highest amount of population-specific Neanderthal segments among worldwide groups. Finally, we demonstrate that most of the genetic variation in Indians stems from a single major migration out of Africa that occurred around 50,000 years ago, with minimal contribution from earlier migration waves. Together, these analyses provide a detailed view of the population history of India and underscore the value of expanding genomic surveys to diverse groups outside Europe.

A Comprehensive Analysis of 3 Moroccan Genomes Revealed Contributions From Both African and European Ancestries

Genetic variations in the human genome represent the differences in DNA sequence within individuals. This highlights the important role of whole human genome sequencing which has become the keystone for precision medicine and disease prediction. Morocco is an important hub for studying human population migration and mixing history. This study presents the analysis of 3 Moroccan genomes; the variant analysis revealed 6 379 606 single nucleotide variants (SNVs) and 1 050 577 small InDels. Of those identified SNVs, 219 152 were novel, with 1233 occurring in coding regions, and 5580 non-synonymous single nucleotide variants (nsSNP) variants were predicted to affect protein functions. The InDels produced 1055 coding variants and 454 non-3n length variants, and their size ranged from -49 and 49 bp. We further analysed the gene pathways of 8 novel coding variants found in the 3 genomes and revealed 5 genes involved in various diseases and biological pathways. We found that the Moroccan genomes share 92.78% of African ancestry, and 92.86% of Non-Finnish European ancestry, according to the gnomAD database. Then, population structure inference, by admixture analysis and network-based approach, revealed that the studied genomes form a mixed population structure, highlighting the increased genetic diversity in Morocco.

STIGMA: Single-cell tissue-specific gene prioritization using machine learning

Clinical exome and genome sequencing have revolutionized the understanding of human disease genetics. Yet many genes remain functionally uncharacterized, complicating the establishment of causal disease links for genetic variants. While several scoring methods have been devised to prioritize these candidate genes, these methods fall short of capturing the expression heterogeneity across cell subpopulations within tissues. Here, we introduce single-cell tissue-specific gene prioritization using machine learning (STIGMA), an approach that leverages single-cell RNA-seq (scRNA-seq) data to prioritize candidate genes associated with rare congenital diseases. STIGMA prioritizes genes by learning the temporal dynamics of gene expression across cell types during healthy organogenesis. To assess the efficacy of our framework, we applied STIGMA to mouse limb and human fetal heart scRNA-seq datasets. In a cohort of individuals with congenital limb malformation, STIGMA prioritized 469 variants in 345 genes, with UBA2 as a notable example. For congenital heart defects, we detected 34 genes harboring nonsynonymous de novo variants (nsDNVs) in two or more individuals from a set of 7,958 individuals, including the ortholog of Prdm1, which is associated with hypoplastic left ventricle and hypoplastic aortic arch. Overall, our findings demonstrate that STIGMA effectively prioritizes tissue-specific candidate genes by utilizing single-cell transcriptome data. The ability to capture the heterogeneity of gene expression across cell populations makes STIGMA a powerful tool for the discovery of disease-associated genes and facilitates the identification of causal variants underlying human genetic disorders.

Phenome- and genome-wide analyses of retinal optical coherence tomography images identify links between ocular and systemic health

The human retina is a multilayered tissue that offers a unique window into systemic health. Optical coherence tomography (OCT) is widely used in eye care and allows the noninvasive, rapid capture of retinal anatomy in exquisite detail. We conducted genotypic and phenotypic analyses of retinal layer thicknesses using macular OCT images from 44,823 UK Biobank participants. We performed OCT layer cross-phenotype association analyses (OCT-XWAS), associating retinal thicknesses with 1866 incident conditions (median 10-year follow-up) and 88 quantitative traits and blood biomarkers. We performed genome-wide association studies (GWASs), identifying inherited genetic markers that influence retinal layer thicknesses and replicated our associations among the LIFE-Adult Study (N = 6313). Last, we performed a comparative analysis of phenome- and genome-wide associations to identify putative causal links between retinal layer thicknesses and both ocular and systemic conditions. Independent associations with incident mortality were detected for thinner photoreceptor segments (PSs) and, separately, ganglion cell complex layers. Phenotypic associations were detected between thinner retinal layers and ocular, neuropsychiatric, cardiometabolic, and pulmonary conditions. A GWAS of retinal layer thicknesses yielded 259 unique loci. Consistency between epidemiologic and genetic associations suggested links between a thinner retinal nerve fiber layer with glaucoma, thinner PS with age-related macular degeneration, and poor cardiometabolic and pulmonary function with a thinner PS. In conclusion, we identified multiple inherited genetic loci and acquired systemic cardio-metabolic-pulmonary conditions associated with thinner retinal layers and identify retinal layers wherein thinning is predictive of future ocular and systemic conditions.

Genome-wide assessment of genetic diversity and transcript variations in 17 accessions of the model diatom Phaeodactylum tricornutum

Diatoms, a prominent group of phytoplankton, have a significant impact on both the oceanic food chain and carbon sequestration, thereby playing a crucial role in regulating the climate. These highly diverse organisms show a wide geographic distribution across various latitudes. In addition to their ecological significance, diatoms represent a vital source of bioactive compounds that are widely used in biotechnology applications. In the present study, we investigated the genetic and transcriptomic diversity of 17 accessions of the model diatom Phaeodactylum tricornutum including those sampled a century ago as well as more recently collected accessions. The analysis of the data reveals a higher genetic diversity and the emergence of novel clades, indicating an increasing diversity within the P. tricornutum population structure, compared to the previous study and a persistent long-term balancing selection of genes in old and newly sampled accessions. However, the study did not establish a clear link between the year of sampling and genetic diversity, thereby, rejecting the hypothesis of loss of heterozygoty in cultured strains. Transcript analysis identified novel transcript including noncoding RNA and other categories of small RNA such as PiwiRNAs. Additionally, transcripts analysis using differential expression as well as Weighted Gene Correlation Network Analysis has provided evidence that the suppression or downregulation of genes cannot be solely attributed to loss-of-function mutations. This implies that other contributing factors, such as epigenetic modifications, may play a crucial role in regulating gene expression. Our study provides novel genetic resources, which are now accessible through the platform PhaeoEpiview (https://PhaeoEpiView.univ-nantes.fr), that offer both ease of use and advanced tools to further investigate microalgae biology and ecology, consequently enriching our current understanding of these organisms.

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.

Novel compound heterozygous pathogenic variants in the SLC3A1 gene in a Chinese family with cystinuria

BACKGROUND

Cystinuria is an autosomal recessive disorder characterized by a cystine transport deficiency in the renal tubules due to mutations in two genes: SLC3A1 and SLC7A9. Cystinuria can be classified into three forms based on the genotype: type A, due to mutations in the SLC3A1 gene; type B, due to mutations in the SLC7A9 gene; and type AB, due to mutations in both genes.

METHODS

We report a 12-year-old boy from central China with cystine stones. He was from a non-consanguineous family that had no known history of genetic disease. A physical examination showed normal development and neurological behaviors. Whole-exome and Sanger sequencing were used to identify and verify the suspected pathogenic variants.

RESULTS

The compound heterozygous variants c.898_905del (p.Arg301AlafsTer6) is located in exon5 and c.1898_1899insAT (p.Asp634LeufsTer46) is located in exon10 of SLC3A1 (NM_000341.4) were deemed responsible for type A cystinuria family. The variant c.898_905del was reported in a Japanese patient in 2000, and the variant c.1898_1899insAT is novel.

CONCLUSION

A novel pathogenic heterozygous variant pair of the SLC3A1 gene was identified in a Chinese boy with type A cystinuria, enriching the mutational spectrum of the SLC3A1 gene. We attempted to find a pattern for the association between the genotype of SLC3A1 variants and the manifestations of cystinuria in patients with different onset ages. Our findings have important implications for genetic counseling and the early clinical diagnosis of cystinuria.

The Landscape of Point Mutations in Human Protein Coding Genes Leading to Pregnancy Loss

Pregnancy loss is the most frequent complication of a pregnancy which is devastating for affected families and poses a significant challenge for the health care system. Genetic factors are known to play an important role in the etiology of pregnancy loss; however, despite advances in diagnostics, the causes remain unexplained in more than 30% of cases. In this review, we aggregated the results of the decade-long studies into the genetic risk factors of pregnancy loss (including miscarriage, termination for fetal abnormality, and recurrent pregnancy loss) in euploid pregnancies, focusing on the spectrum of point mutations associated with these conditions. We reviewed the evolution of molecular genetics methods used for the genetic research into causes of pregnancy loss, and collected information about 270 individual genetic variants in 196 unique genes reported as genetic cause of pregnancy loss. Among these, variants in 18 genes have been reported by multiple studies, and two or more variants were reported as causing pregnancy loss for 57 genes. Further analysis of the properties of all known pregnancy loss genes showed that they correspond to broadly expressed, highly evolutionary conserved genes involved in crucial cell differentiation and developmental processes and related signaling pathways. Given the features of known genes, we made an effort to construct a list of candidate genes, variants in which may be expected to contribute to pregnancy loss. We believe that our results may be useful for prediction of pregnancy loss risk in couples, as well as for further investigation and revealing genetic etiology of pregnancy loss.

Whole-Exome Sequencing Among School-Aged Children With High Myopia

Importance

High myopia (HM) is one of the leading causes of visual impairment worldwide. Genetic factors are known to play an important role in the development of HM.

Objective

To identify risk variants in a large HM cohort and to examine the implications of genetic testing of schoolchildren with HM.

Design, Setting, and Participants

This cohort study retrospectively reviewed whole-exome sequencing (WES) results in 6215 schoolchildren with HM who underwent genetic testing between September 2019 and July 2020 in Wenzhou City, China. HM is defined as a spherical equivalent refraction (SER) of -6.00 diopters (D) or less. The study setting was a genetic testing laboratory and a multicenter school census. Data were analyzed from July 2021 to June 2022.

Main Outcomes and Measures

The frequency and distribution of positive germline variants, the percentage of individuals with HM in both eyes, and subsequent variant yield for common high myopia (CHM; -8.00 D ≤ SER ≤ -6.00 D), ultra myopia (UM; -10.00 D ≤ SER < -8.00 D), and extreme myopia (EM; SER < -10.00 D).

Results

Of the 6215 schoolchildren with HM, 3278 (52.74%) were male. Their mean (SD) age was 14.87 (2.02) years, including 355 students in primary school, 1970 in junior high school, and 3890 in senior high school. The mean (SD) SER was -7.51 (-1.36) D for the right eye and -7.46 (-1.34) D for the left eye. Among schoolchildren with HM, genetic testing yielded 271 potential pathogenic variants in 75 HM candidate genes in 964 diagnoses (15.52%). A total of 36 known variants were found in 490 HM participants (7.88%) and 235 protein-truncating variants (PTVs) in 506 participants (8.14%). Involved variant yield was significantly positively associated with SER (Cochran-Armitage test for trend Z = 2.5492; P = .01), which ranged from 7.66% in the CHM group, 8.70% in the UM group, to 11.90% in the EM group. We also found that primary school students with EM had the highest variant yield of PTVs (8 of 35 students [22.86%]), which was 1.77 and 4.78 times that of the UM and CHM, respectively.

Conclusions and Relevance

In this cohort study of WES for HM, several potential pathogenic variants were identified in a substantial number of schoolchildren with HM. The high variation frequency in younger students with EM can provide clues for genetic screening and clinical examinations of HM to promote long-term follow-up assessment.

Association Analyses of Predicted Loss-of-Function Variants Prioritized 15 Genes as Blood Pressure Regulators

BACKGROUND

Hypertension, clinically defined by elevated blood pressure (BP), is an important cause of mortality and morbidity worldwide. Many risk factors for hypertension are known, including a positive family history, which suggests that genetics contribute to interindividual BP variation. Genome-wide association studies (GWAS) have identified > 1000 loci associated with BP, yet the identity of the genes responsible for these associations remains largely unknown.

METHODS

To pinpoint genes that causally affect variation of BP in humans, we analyzed predicted loss-of-function (pLoF) variants in the UK Biobank whole-exome sequencing dataset (n = 454,709 participants, 6% non-European ancestry). We analyzed genetic associations between systolic or diastolic BP (SBP/DBP) and single pLoF variants (additive and recessive genetic models) as well as with the burden of very rare pLoF variants (minor allele frequency [MAF] < 0.01%).

RESULTS

Single pLoF variants in 10 genes were associated with BP (ANKDD1B, ENPEP, PNCK, BTN3A2, C1orf145 [OBSCN-AS1], CASP9, DBH, KIAA1161 [MYORG], OR4X1, and TMC3). We also found a burden of rare pLoF variants in 5 additional genes associated with BP (TTN, NOS3, FES, SMAD6, COL21A1). Except for PNCK, which is located on the X-chromosome, these genes map near variants previously associated with BP by GWAS, validating the study of pLoF variants to prioritize causal genes at GWAS loci.

CONCLUSIONS

Our study highlights 15 genes that likely modulate BP in humans, including 5 genes that harbour pLoF variants associated with lower BP.

Exploring the genetic diversity of the Japanese population: Insights from a large-scale whole genome sequencing analysis

The Japanese archipelago is a terminal location for human migration, and the contemporary Japanese people represent a unique population whose genomic diversity has been shaped by multiple migrations from Eurasia. We analyzed the genomic characteristics that define the genetic makeup of the modern Japanese population from a population genetics perspective from the genomic data of 9,287 samples obtained by high-coverage whole-genome sequencing (WGS) by the National Center Biobank Network. The dataset comprised populations from the Ryukyu Islands and other parts of the Japanese archipelago (Hondo). The Hondo population underwent two episodes of population decline during the Jomon period, corresponding to the Late Neolithic, and the Edo period, corresponding to the Early Modern era, while the Ryukyu population experienced a population decline during the shell midden period of the Late Neolithic in this region. Haplotype analysis suggested increased allele frequencies for genes related to alcohol and fatty acid metabolism, which were reported as loci that had experienced positive natural selection. Two genes related to alcohol metabolism were found to be 12,500 years out of phase with the time when they began to increase in the allele frequency; this finding indicates that the genomic diversity of Japanese people has been shaped by events closely related to agriculture and food production.

Research participants' perspectives regarding the feedback of secondary findings-A cohort from the DDD-Africa study, South Africa

Genomic researchers face an ethical dilemma regarding feedback of individual results generated from genomic studies. In the African setting, genomic research is still not widely implemented and, coupled with this, the limited African-specific guidelines on how to feedback on individual research findings. A qualitative study was performed to assess participants' expectations and preferences regarding the feedback of secondary findings from genomic research. Participants were parents of children with a developmental disorder, enrolled in the Deciphering Developmental Disorders in Africa (DDD-Africa) research project, and were purposefully selected. Three deliberative focus group discussions were conducted with 14 participants. Each deliberative focus group consisted of two separate audio-recorded interviews and presented different case scenarios for different types of secondary findings that could be theoretically detected during genomic research. We aimed to explore participants' preferences for the extent, nature, timing, and methods for receiving individual study results, specifically pertaining to secondary findings. Thematic content analysis was done, with a deductive approach to coding. Four themes emerged which included participants' perception of readiness to receive secondary findings, queries raised around who has access to research findings and feedback of findings consent, responsibilities of the researcher, and reasons for not wanting/wanting secondary findings. Overall, participants expressed that they want to receive feedback on secondary findings irrespective of disease severity and treatment availability. Lifestyle changes, early prevention or treatment, impact on future generations, and preparedness were strong motivations for wanting feedback on results. Participants felt that when the research involved minors, it was the parents' right to receive results on behalf of their children. This study provides new insights into participants' preferences around feedback on genomic research results and could serve as an important basis for creating guidelines and recommendations for feedback on genomic results in the African context.

A human mitofusin 2 mutation can cause mitophagic cardiomyopathy

Cardiac muscle has the highest mitochondrial density of any human tissue, but mitochondrial dysfunction is not a recognized cause of isolated cardiomyopathy. Here, we determined that the rare mitofusin (MFN) 2 R400Q mutation is 15-20× over-represented in clinical cardiomyopathy, whereas this specific mutation is not reported as a cause of MFN2 mutant-induced peripheral neuropathy, Charcot-Marie-Tooth disease type 2A (CMT2A). Accordingly, we interrogated the enzymatic, biophysical, and functional characteristics of MFN2 Q400 versus wild-type and CMT2A-causing MFN2 mutants. All MFN2 mutants had impaired mitochondrial fusion, the canonical MFN2 function. Compared to MFN2 T105M that lacked catalytic GTPase activity and exhibited normal activation-induced changes in conformation, MFN2 R400Q and M376A had normal GTPase activity with impaired conformational shifting. MFN2 R400Q did not suppress mitochondrial motility, provoke mitochondrial depolarization, or dominantly suppress mitochondrial respiration like MFN2 T105M. By contrast to MFN2 T105M and M376A, MFN2 R400Q was uniquely defective in recruiting Parkin to mitochondria. CRISPR editing of the R400Q mutation into the mouse Mfn2 gene induced perinatal cardiomyopathy with no other organ involvement; knock-in of Mfn2 T105M or M376V did not affect the heart. RNA sequencing and metabolomics of cardiomyopathic Mfn2 Q/Q400 hearts revealed signature abnormalities recapitulating experimental mitophagic cardiomyopathy. Indeed, cultured cardiomyoblasts and in vivo cardiomyocytes expressing MFN2 Q400 had mitophagy defects with increased sensitivity to doxorubicin. MFN2 R400Q is the first known natural mitophagy-defective MFN2 mutant. Its unique profile of dysfunction evokes mitophagic cardiomyopathy, suggesting a mechanism for enrichment in clinical cardiomyopathy.

rs67047829 genotypes of ERV3-1/ZNF117 are associated with lower body mass index in the Polish population

There is now substantial evidence that zinc-finger proteins are implicated in adiposity. Aims were to datamine for high-frequency (near-neutral selection) pretermination-codon (PTC) single-nucleotide polymorphisms (SNPs; n = 141) from a database with > 550,000 variants and analyze possible association with body mass index in a large Polish sample (n = 5757). BMI was regressed (males/females together or separately) against genetic models. Regression for rs67047829 uncovered an interaction-independent association with BMI with both sexes together: mean ± standard deviation, kg/m2: [G];[G], 25.4 ± 4.59 (n = 3650); [G](;)[A], 25.0 ± 4.28 (n = 731); [A];[A], 23.4 ± 3.60 (n = 44); additive model adjusted for age and sex: p = 4.08 × 10-5; beta: - 0.0458, 95% confidence interval (CI) - 0.0732 : - 0.0183; surviving Bonferroni correction; for males: [G];[G], 24.8 ± 4.94 (n = 1878); [G](;)[A], 24.2 ± 4.31 (n = 386); [A];[A], 22.4 ± 3.69 (n = 23); p = 4.20 × 10-4; beta: - 0.0573, CI - 0.0947 : - 0.0199. For average-height males the difference between [G];[G] and [A];[A] genotypes would correspond to ~ 6 kg, suggesting considerable protection against increased BMI. rs67047829 gives a pretermination codon in ERV3-1 which shares an exonic region and possibly promoter with ZNF117, previously associated with adiposity and type-2 diabetes. As this result occurs in a near-neutral Mendelian setting, a drug targetting ERV3-1/ZNF117 might potentially provide considerable benefits with minimal side-effects. This result needs to be replicated, followed by analyses of splice-variant mRNAs and protein expression.

Learning the kernel for rare variant genetic association test

Introduction: Compared to Genome-Wide Association Studies (GWAS) for common variants, single-marker association analysis for rare variants is underpowered. Set-based association analyses for rare variants are powerful tools that capture some of the missing heritability in trait association studies. Methods: We extend the convex-optimized SKAT (cSKAT) test set procedure which learns from data the optimal convex combination of kernels, to the full Generalised Linear Model (GLM) setting with arbitrary non-genetic covariates. We call this extended cSKAT (ecSKAT) and show that the resulting optimization problem is a quadratic programming problem that can be solved with no additional cost compared to cSKAT. Results: We show that a modified objective is related to an upper bound for the p-value through a decreasing exponential term in the objective function, indicating that optimizing this objective function is a principled way of learning the combination of kernels. We evaluate the performance of the proposed method on continuous and binary traits using simulation studies and illustrate its application using UK Biobank Whole Exome Sequencing data on hand grip strength and systemic lupus erythematosus rare variant association analysis. Discussion: Our proposed ecSKAT method enables correcting for important confounders in association studies such as age, sex or population structure for both quantitative and binary traits. Simulation studies showed that ecSKAT can recover sensible weights and achieve higher power across different sample sizes and misspecification settings. Compared to the burden test and SKAT method, ecSKAT gives a lower p-value for the genes tested in both quantitative and binary traits in the UKBiobank cohort.

The genomic footprint of whaling and isolation in fin whale populations

Twentieth century industrial whaling pushed several species to the brink of extinction, with fin whales being the most impacted. However, a small, resident population in the Gulf of California was not targeted by whaling. Here, we analyzed 50 whole-genomes from the Eastern North Pacific (ENP) and Gulf of California (GOC) fin whale populations to investigate their demographic history and the genomic effects of natural and human-induced bottlenecks. We show that the two populations diverged ~16,000 years ago, after which the ENP population expanded and then suffered a 99% reduction in effective size during the whaling period. In contrast, the GOC population remained small and isolated, receiving less than one migrant per generation. However, this low level of migration has been crucial for maintaining its viability. Our study exposes the severity of whaling, emphasizes the importance of migration, and demonstrates the use of genome-based analyses and simulations to inform conservation strategies.

Advanced variant classification framework reduces the false positive rate of predicted loss-of-function variants in population sequencing data

Predicted loss of function (pLoF) variants are often highly deleterious and play an important role in disease biology, but many pLoF variants may not result in loss of function (LoF). Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines' PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 genes associated with autosomal-recessive disease from the Genome Aggregation Database (gnomAD v.2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in a low proportion expressed across transcripts (pext) scored region, or the presence of cryptic in-frame splice rescues. Variants predicted to evade LoF or to be potential artifacts were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of pLoF variants predicted as likely not LoF/not LoF, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

Pathogenic variants of human GABRA1 gene associated with epilepsy: A computational approach

Critical for brain development, neurodevelopmental and network disorders, the GABRA1 gene encodes for the α1 subunit, an abundantly and developmentally expressed subunit of heteropentameric gamma-aminobutyric acid A receptors (GABAARs) mediating primary inhibition in the brain. Mutations of the GABAAR subunit genes including GABRA1 gene are associated with epilepsy, a group of syndromes, characterized by unprovoked seizures and diagnosed by integrative approach, that involves genetic testing. Despite the diagnostic use of genetic testing, a large fraction of the GABAAR subunit gene variants including the variants of GABRA1 gene is not known in terms of their molecular consequence, a challenge for precision and personalized medicine. Addressing this, one hundred thirty-seven GABRA1 gene variants of unknown clinical significance have been extracted from the ClinVar database and computationally analyzed for pathogenicity. Eight variants (L49H, P59L, W97R, D99G, G152S, V270G, T294R, P305L) are predicted as pathogenic and mapped to the α1 subunit's extracellular domain (ECD), transmembrane domains (TMDs) and extracellular linker. This is followed by the integration with relevant data for cellular pathology and severity of the epilepsy syndromes retrieved from the literature. Our results suggest that the pathogenic variants in the ECD of GABRA1 (L49H, P59L, W97R, D99G, G152S) will probably manifest decreased surface expression and reduced current with mild epilepsy phenotypes while V270G, T294R in the TMDs and P305L in the linker between the second and the third TMDs will likely cause reduced cell current with severe epilepsy phenotypes. The results presented in this study provides insights for clinical genetics and wet lab experimentation.

A knowledge graph approach to predict and interpret disease-causing gene interactions

BACKGROUND

Understanding the impact of gene interactions on disease phenotypes is increasingly recognised as a crucial aspect of genetic disease research. This trend is reflected by the growing amount of clinical research on oligogenic diseases, where disease manifestations are influenced by combinations of variants on a few specific genes. Although statistical machine-learning methods have been developed to identify relevant genetic variant or gene combinations associated with oligogenic diseases, they rely on abstract features and black-box models, posing challenges to interpretability for medical experts and impeding their ability to comprehend and validate predictions. In this work, we present a novel, interpretable predictive approach based on a knowledge graph that not only provides accurate predictions of disease-causing gene interactions but also offers explanations for these results.

RESULTS

We introduce BOCK, a knowledge graph constructed to explore disease-causing genetic interactions, integrating curated information on oligogenic diseases from clinical cases with relevant biomedical networks and ontologies. Using this graph, we developed a novel predictive framework based on heterogenous paths connecting gene pairs. This method trains an interpretable decision set model that not only accurately predicts pathogenic gene interactions, but also unveils the patterns associated with these diseases. A unique aspect of our approach is its ability to offer, along with each positive prediction, explanations in the form of subgraphs, revealing the specific entities and relationships that led to each pathogenic prediction.

CONCLUSION

Our method, built with interpretability in mind, leverages heterogenous path information in knowledge graphs to predict pathogenic gene interactions and generate meaningful explanations. This not only broadens our understanding of the molecular mechanisms underlying oligogenic diseases, but also presents a novel application of knowledge graphs in creating more transparent and insightful predictors for genetic research.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Clinicogenomic characteristics and synthetic lethal implications of germline homologous recombination-deficient hepatocellular carcinoma

BACKGROUNDS AND AIMS

We performed an in-depth examination of pathogenic germline variants (PGVs) and somatic variants in DNA damage response (DDR) genes in hepatocellular carcinoma (HCC) to explore their clinical and genomic impacts.

APPROACH AND RESULTS

We used a merged whole-exome or RNA sequencing data set derived from in-house ( n = 230) and The Cancer Genome Atlas ( n = 362) databases of multiethnic HCC samples. We also evaluated synthetic lethal approaches targeting mutations in homologous recombination (HR) genes using HCC cells selected from five genomic databases of cancer cell lines. A total of 110 PGVs in DDR pathways in 96 patients were selected. Of the PGV carriers, 44 were HR-altered and found to be independently associated with poorer disease-free survival after hepatectomy. The most frequently altered HR gene in both germline and somatic tissues was POLQ , and this variant was detected in 22.7% (10/44) and 23.8% (5/21) of all the corresponding carriers, respectively. PGVs in HR were significantly associated with upregulation of proliferation and replication-related genes and familial risk of HCC. Samples harboring PGVs in HR with loss of heterozygosity were most strongly correlated with the genomic footprints of deficient HR, such as mutation burden and denovoSig2 (analogous to Catalogue of Somatic Mutations in Cancer [COSMIC] 3), and poor outcome. Pharmacologic experiments with HCC cells defective in BRCA2 or POLQ suggested that tumors with this phenotype are synthetic lethal with poly(ADP-ribose) polymerase inhibitors.

CONCLUSIONS

Our findings suggest that germline HR defects in HCC tend to confer a poor prognosis and result in distinctive genomic scarring. Tests of the clinical benefits of HR-directed treatments in the affected patients are needed.

The Spectrum of the Heterozygous Effect in Biallelic Mendelian Diseases-The Symptomatic Heterozygote Issue

Heterozygous carriers of pathogenic/likely pathogenic variants in autosomal recessive disorders seem to be asymptomatic. However, in recent years, an increasing number of case reports have suggested that mild and unspecific symptoms can occur in some heterozygotes, as symptomatic heterozygotes have been identified across different disease types, including neurological, neuromuscular, hematological, and pulmonary diseases. The symptoms are usually milder in heterozygotes than in biallelic variants and occur "later in life". The status of symptomatic heterozygotes as separate entities is often disputed, and alternative diagnoses are considered. Indeed, often only a thin line exists between dual, dominant, and recessive modes of inheritance and symptomatic heterozygosity. Interestingly, recent population studies have found global disease effects in heterozygous carriers of some genetic variants. What makes the few heterozygotes symptomatic, while the majority show no symptoms? The molecular basis of this phenomenon is still unknown. Possible explanations include undiscovered deep-splicing variants, genetic and environmental modifiers, digenic/oligogenic inheritance, skewed methylation patterns, and mutational burden. Symptomatic heterozygotes are rarely reported in the literature, mainly because most did not undergo the complete diagnostic procedure, so alternative diagnoses could not be conclusively excluded. However, despite the increasing accessibility to high-throughput technologies, there still seems to be a small group of patients with mild symptoms and just one variant of autosomes in biallelic diseases. Here, we present some examples, the current state of knowledge, and possible explanations for this phenomenon, and thus argue against the existing dominant/recessive classification.

Chromatin regulators in the TBX1 network confer risk for conotruncal heart defects in 22q11.2DS

Congenital heart disease (CHD) affecting the conotruncal region of the heart, occurs in 40-50% of patients with 22q11.2 deletion syndrome (22q11.2DS). This syndrome is a rare disorder with relative genetic homogeneity that can facilitate identification of genetic modifiers. Haploinsufficiency of TBX1, encoding a T-box transcription factor, is one of the main genes responsible for the etiology of the syndrome. We suggest that genetic modifiers of conotruncal defects in patients with 22q11.2DS may be in the TBX1 gene network. To identify genetic modifiers, we analyzed rare, predicted damaging variants in whole genome sequence of 456 cases with conotruncal defects and 537 controls, with 22q11.2DS. We then performed gene set approaches and identified chromatin regulatory genes as modifiers. Chromatin genes with recurrent damaging variants include EP400, KAT6A, KMT2C, KMT2D, NSD1, CHD7 and PHF21A. In total, we identified 37 chromatin regulatory genes, that may increase risk for conotruncal heart defects in 8.5% of 22q11.2DS cases. Many of these genes were identified as risk factors for sporadic CHD in the general population. These genes are co-expressed in cardiac progenitor cells with TBX1, suggesting that they may be in the same genetic network. The genes KAT6A, KMT2C, CHD7 and EZH2, have been previously shown to genetically interact with TBX1 in mouse models. Our findings indicate that disturbance of chromatin regulatory genes impact the TBX1 gene network serving as genetic modifiers of 22q11.2DS and sporadic CHD, suggesting that there are some shared mechanisms involving the TBX1 gene network in the etiology of CHD.

eyeVarP: A computational framework for the identification of pathogenic variants specific to eye disease

PURPOSE

Disease-specific pathogenic variant prediction tools that differentiate pathogenic variants from benign have been improved through disease specificity recently. However, they have not been evaluated on disease-specific pathogenic variants compared with other diseases, which would help to prioritize disease-specific variants from several genes or novel genes. Thus, we hypothesize that features of pathogenic variants alone would provide a better model.

METHODS

We developed an eye disease-specific variant prioritization tool (eyeVarP), which applied the random forest algorithm to the data set of pathogenic variants of eye diseases and other diseases. We also developed the VarP tool and generalized pipeline to filter missense and insertion-deletion variants and predict their pathogenicity from exome or genome sequencing data, thus we provide a complete computational procedure.

RESULTS

eyeVarP outperformed pan disease-specific tools in identifying eye disease-specific pathogenic variants under the top 10. VarP outperformed 12 pathogenicity prediction tools with an accuracy of 95% in correctly identifying the pathogenicity of missense and insertion-deletion variants. The complete pipeline would help to develop disease-specific tools for other genetic disorders.

CONCLUSION

eyeVarP performs better in identifying eye disease-specific pathogenic variants using pathogenic variant features and gene features. Implementing such complete computational procedure would significantly improve the clinical variant interpretation for specific diseases.

COVID-19 severity: does the genetic landscape of rare variants matter?

Rare variants affecting host defense against pathogens may be involved in COVID-19 severity, but most rare variants are not expected to have a major impact on the course of COVID-19. We hypothesized that the accumulation of weak effects of many rare functional variants throughout the exome may contribute to the overall risk in patients with severe disease. This assumption is consistent with the omnigenic model of the relationship between genetic and phenotypic variation in complex traits, according to which association signals tend to spread across most of the genome through gene regulatory networks from genes outside the major pathways to disease-related genes. We performed whole-exome sequencing and compared the burden of rare variants in 57 patients with severe and 29 patients with mild/moderate COVID-19. At the whole-exome level, we observed an excess of rare, predominantly high-impact (HI) variants in the group with severe COVID-19. Restriction to genes intolerant to HI or damaging missense variants increased enrichment for these classes of variants. Among various sets of genes, an increased signal of rare HI variants was demonstrated predominantly for primary immunodeficiency genes and the entire set of genes associated with immune diseases, as well as for genes associated with respiratory diseases. We advocate taking the ideas of the omnigenic model into account in COVID-19 studies.

Buffering of genetic dominance by allele-specific protein complex assembly

Protein complex assembly often occurs while subunits are being translated, resulting in complexes whose subunits were translated from the same mRNA in an allele-specific manner. It has thus been hypothesized that such cotranslational assembly may counter the assembly-mediated dominant-negative effect, whereby co-assembly of mutant and wild-type subunits "poisons" complex activity. Here, we show that cotranslationally assembling subunits are much less likely to be associated with autosomal dominant relative to recessive disorders, and that subunits with dominant-negative disease mutations are significantly depleted in cotranslational assembly compared to those associated with loss-of-function mutations. We also find that complexes with known dominant-negative effects tend to expose their interfaces late during translation, lessening the likelihood of cotranslational assembly. Finally, by combining complex properties with other features, we trained a computational model for predicting proteins likely to be associated with non-loss-of-function disease mechanisms, which we believe will be of considerable utility for protein variant interpretation.

Chromosome-level, nanopore-only genome and allele-specific DNA methylation of Pallas's cat, Otocolobus manul

Pallas's cat, or the manul cat (Otocolobus manul), is a small felid native to the grasslands and steppes of central Asia. Population strongholds in Mongolia and China face growing challenges from climate change, habitat fragmentation, poaching, and other sources. These threats, combined with O. manul's zoo collection popularity and value in evolutionary biology, necessitate improvement of species genomic resources. We used standalone nanopore sequencing to assemble a 2.5 Gb, 61-contig nuclear assembly and 17097 bp mitogenome for O. manul. The primary nuclear assembly had 56× sequencing coverage, a contig N50 of 118 Mb, and a 94.7% BUSCO completeness score for Carnivora-specific genes. High genome collinearity within Felidae permitted alignment-based scaffolding onto the fishing cat (Prionailurus viverrinus) reference genome. Manul contigs spanned all 19 felid chromosomes with an inferred total gap length of less than 400 kilobases. Modified basecalling and variant phasing produced an alternate pseudohaplotype assembly and allele-specific DNA methylation calls; 61 differentially methylated regions were identified between haplotypes. Nearest features included classical imprinted genes, non-coding RNAs, and putative novel imprinted loci. The assembled mitogenome successfully resolved existing discordance between Felinae nuclear and mtDNA phylogenies. All assembly drafts were generated from 158 Gb of sequence using seven minION flow cells.

Insights into human health from phenome- and genome-wide analyses of UK Biobank retinal optical coherence tomography phenotypes

The human retina is a complex multi-layered tissue which offers a unique window into systemic health and disease. Optical coherence tomography (OCT) is widely used in eye care and allows the non-invasive, rapid capture of retinal measurements in exquisite detail. We conducted genome- and phenome-wide analyses of retinal layer thicknesses using macular OCT images from 44,823 UK Biobank participants. We performed phenome-wide association analyses, associating retinal thicknesses with 1,866 incident ICD-based conditions (median 10-year follow-up) and 88 quantitative traits and blood biomarkers. We performed genome-wide association analyses, identifying inherited genetic markers which influence the retina, and replicated our associations among 6,313 individuals from the LIFE-Adult Study. And lastly, we performed comparative association of phenome- and genome- wide associations to identify putative causal links between systemic conditions, retinal layer thicknesses, and ocular disease. Independent associations with incident mortality were detected for photoreceptor thinning and ganglion cell complex thinning. Significant phenotypic associations were detected between retinal layer thinning and ocular, neuropsychiatric, cardiometabolic and pulmonary conditions. Genome-wide association of retinal layer thicknesses yielded 259 loci. Consistency between epidemiologic and genetic associations suggested putative causal links between thinning of the retinal nerve fiber layer with glaucoma, photoreceptor segment with AMD, as well as poor cardiometabolic and pulmonary function with PS thinning, among other findings. In conclusion, retinal layer thinning predicts risk of future ocular and systemic disease. Furthermore, systemic cardio-metabolic-pulmonary conditions promote retinal thinning. Retinal imaging biomarkers, integrated into electronic health records, may inform risk prediction and potential therapeutic strategies.

Structural Dynamics Predominantly Determine the Adaptability of Proteins to Amino Acid Deletions

The insertion or deletion (indel) of amino acids has a variety of effects on protein function, ranging from disease-forming changes to gaining new functions. Despite their importance, indels have not been systematically characterized towards protein engineering or modification goals. In the present work, we focus on deletions composed of multiple contiguous amino acids (mAA-dels) and their effects on the protein (mutant) folding ability. Our analysis reveals that the mutant retains the native fold when the mAA-del obeys well-defined structural dynamics properties: localization in intrinsically flexible regions, showing low resistance to mechanical stress, and separation from allosteric signaling paths. Motivated by the possibility of distinguishing the features that underlie the adaptability of proteins to mAA-dels, and by the rapid evaluation of these features using elastic network models, we developed a positive-unlabeled learning-based classifier that can be adopted for protein design purposes. Trained on a consolidated set of features, including those reflecting the intrinsic dynamics of the regions where the mAA-dels occur, the new classifier yields a high recall of 84.3% for identifying mAA-dels that are stably tolerated by the protein. The comparative examination of the relative contribution of different features to the prediction reveals the dominant role of structural dynamics in enabling the adaptation of the mutant to mAA-del without disrupting the native fold.

A draft human pangenome reference

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals1. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.

Reconciling Mouse and Human Immunology at the Altar of Genetics

Immunity to infection has been extensively studied in humans and mice bearing naturally occurring or experimentally introduced germline mutations. Mouse studies are sometimes neglected by human immunologists, on the basis that mice are not humans and the infections studied are experimental and not natural. Conversely, human studies are sometimes neglected by mouse immunologists, on the basis of the uncontrolled conditions of study and small numbers of patients. However, both sides would agree that the infectious phenotypes of patients with inborn errors of immunity often differ from those of the corresponding mutant mice. Why is that? We argue that this important question is best addressed by revisiting and reinterpreting the findings of both mouse and human studies from a genetic perspective. Greater caution is required for reverse-genetics studies than for forward-genetics studies, but genetic analysis is sufficiently strong to define the studies likely to stand the test of time. Genetically robust mouse and human studies can provide invaluable complementary insights into the mechanisms of immunity to infection common and specific to these two species.

Advanced variant classification framework reduces the false positive rate of predicted loss of function (pLoF) variants in population sequencing data

Predicted loss of function (pLoF) variants are highly deleterious and play an important role in disease biology, but many of these variants may not actually result in loss-of-function. Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines's PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 autosomal recessive disease-genes from the Genome Aggregation Database (gnomAD, v2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in low per-base expression (pext) score regions, or the presence of cryptic splice rescues. Variants predicted to be potential artifacts or to evade LoF were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of LoF evading variants assessed, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines, and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

Effects of protein-coding variants on blood metabolite measurements and clinical biomarkers in the UK Biobank

Genome-wide association studies (GWASs) have established the contribution of common and low-frequency variants to metabolic blood measurements in the UK Biobank (UKB). To complement existing GWAS findings, we assessed the contribution of rare protein-coding variants in relation to 355 metabolic blood measurements-including 325 predominantly lipid-related nuclear magnetic resonance (NMR)-derived blood metabolite measurements (Nightingale Health Plc) and 30 clinical blood biomarkers-using 412,393 exome sequences from four genetically diverse ancestries in the UKB. Gene-level collapsing analyses were conducted to evaluate a diverse range of rare-variant architectures for the metabolic blood measurements. Altogether, we identified significant associations (p < 1 × 10^-8) for 205 distinct genes that involved 1,968 significant relationships for the Nightingale blood metabolite measurements and 331 for the clinical blood biomarkers. These include associations for rare non-synonymous variants in PLIN1 and CREB3L3 with lipid metabolite measurements and SYT7 with creatinine, among others, which may not only provide insights into novel biology but also deepen our understanding of established disease mechanisms. Of the study-wide significant clinical biomarker associations, 40% were not previously detected on analyzing coding variants in a GWAS in the same cohort, reinforcing the importance of studying rare variation to fully understand the genetic architecture of metabolic blood measurements.

Point mutations of homologs as an adaptive solution to the gene loss

Gene loss is common and influences genome evolution trajectories. Multiple adaptive strategies to compensate for gene loss have been observed, including copy number gain of paralogous genes and mutations in genes of the same pathway. By using the Ubl-specific protease 2 (ULP2) eviction model, we identify compensatory mutations in the homologous gene ULP1 by laboratory evolution and find that these mutations are capable of rescuing defects caused by the loss of ULP2. Furthermore, bioinformatics analysis of genomes of yeast gene knockout library and natural yeast isolate datasets suggests that point mutations of a homologous gene might be an additional mechanism to compensate gene loss.