Identification of a large set of rare complete human knockouts

A Rare Cause of Primary Microcephaly: 4 New Variants in CDK5RAP2 Gene and Review of the Literature

Autosomal recessive primary microcephaly (MCPH) is a rare, genetically heterogeneous disorder characterized by congenital microcephaly, non-progressive intellectual disability, and absence of neurological abnormalities. Pathogenic variants in CDK5RAP2, linked to MCPH3, represent one of the least common causes of MCPH. Autosomal recessive primary microcephaly (MCPH) is a rare, genetically heterogeneous disorder characterized by reduced head size at birth, variable intellectual disability, and no neurological abnormalities. Among This study aimed to identify and characterize novel genetic and clinical findings in patients with CDK5RAP2 gene variants, contributing to the understanding of this rare disorder. Whole exome sequencing (WES) was conducted in 11 patients from five consanguineous families. Six CDK5RAP2 variants were identified, four of which were novel (c.4421del, c.1968G>C, c.3460C>T, c.625dup). Ten patients harbored homozygous variants, whereas one displayed compound heterozygosity. Segregation analysis confirmed carrier status in parents. Clinical evaluations aligned with typical MCPH3 features, though phenotypic variability was observed. This study expands the CDK5RAP2 variant spectrum and reinforces WES as a critical tool for diagnosing rare MCPH subtypes, guiding carrier screening, and improving genetic counseling. The novel variants highlight the genetic diversity underlying MCPH3, urging broader genomic investigations in undiagnosed cases.

Bi-allelic variants in BRF2 are associated with perinatal death and craniofacial anomalies

BACKGROUND

Variants in genes encoding multiple subunits of the RNA Polymerase III complex which synthesizes rRNAs, tRNAs, and other small RNAs were previously associated with neurological disorders, such as syndromic hypomyelination leukodystrophies, pontocerebellar hypoplasia, and cerebellofaciodental syndrome. One new such candidate is BRF2, which encodes a TFIIB-like factor that recruits the RNA polymerase III complex to type 3 promoters to initiate transcription of U6, RnaseP, and 7SK RNAs.

METHODS

We combined sequencing with functional analyses to investigate the effects of BRF2 variants.

RESULTS

We observe that a previously reported significant underrepresentation of double transmission of a splice variant results in recessive lethality in three large Icelandic families with multiple perinatal losses. Using data aggregation, we identified an additional seven individuals worldwide from three unrelated families carrying biallelic variants in BRF2. Affected individuals present a variable phenotype ranging from severe craniofacial anomalies with early death to intellectual disability with motor and speech development. In silico 3D modelling and functional analyses showed functional impairment of the identified variants, e.g., differences in target loci occupancy. Zebrafish knocked down for the orthologous brf2 presented with abnormal escape response, reduced swimming velocity and head size, and craniofacial malformations. These defects were complemented by the human wild-type but not mutated BRF2 mRNA further demonstrating their deleteriousness.

CONCLUSIONS

Overall, our results support the association of biallelic BRF2 variants with a novel neurodevelopmental disease and provide an additional link between RNA polymerase III, its targets and craniofacial anomalies.

Genomics of rare diseases in the Greater Middle East

The Greater Middle East (GME) represents a concentrated region of unparalleled genetic diversity, characterized by an abundance of distinct alleles, founder mutations and extensive autozygosity driven by high consanguinity rates. These genetic hallmarks present a unique, yet vastly untapped resource for genomic research on Mendelian diseases. Despite this immense potential, the GME continues to face substantial challenges in comprehensive data collection and analysis. This Perspective highlights the region's unique position as a natural laboratory for genetic discovery and explores the challenges that have stifled progress thus far. Importantly, we propose strategic solutions, advocating for an all-inclusive research approach. With targeted investment and focused efforts, the latent genetic wealth in the GME can be transformed into a global hub for genomic research that will redefine and advance our understanding of the human genome.

PON-P3: Accurate Prediction of Pathogenicity of Amino Acid Substitutions

Different types of information are combined during variation interpretation. Computational predictors, most often pathogenicity predictors, provide one type of information for this purpose. These tools are based on various kinds of algorithms. Although the American College of Genetics and the Association for Molecular Pathology guidelines classify variants into five categories, practically all pathogenicity predictors provide binary pathogenic/benign predictions. We developed a novel artificial intelligence-based tool, PON-P3, on the basis of a carefully selected training dataset, meticulous feature selection, and optimization. We started with 1526 features describing variations, their sequence and structural context, and parameters for the affected genes and proteins. The final random boosting method was tested and compared with a total of 23 predictors. PON-P3 performed better than recently introduced predictors, which utilize large language models or structural predictions. PON-P3 was better than methods that use evolutionary data alone or in combination with different gene and protein properties. PON-P3 classifies cases into three categories as benign, pathogenic, and variants of uncertain significance (VUSs). When binary test data were used, some metapredictors performed slightly better than PON-P3; however, in real-life situations, with patient data, those methods overpredict both pathogenic and benign cases. We predicted with PON-P3 all possible amino acid substitutions in all human proteins encoded from MANE transcripts. The method was also used to predict all unambiguous VUSs (i.e., without conflicts) in ClinVar. A total of 12.9% were predicted to be pathogenic, and 49.9% were benign.

There will always be variants of uncertain significance. Analysis of VUSs

The ACMG/AMP guidelines include five categories of which variants of uncertain significance (VUSs) have received increasing attention. Recently, Fowler and Rehm claimed that all or most VUSs could be reclassified as pathogenic or benign within few years. To test this claim, we collected validated benign, pathogenic, VUS and conflicting variants from ClinVar and LOVD and investigated differences at gene, protein, structure, and variant levels. The gene and protein features included inheritance patterns, actionability, functional categories for housekeeping, essential, complete knockout, lethality and haploinsufficient proteins, Gene Ontology annotations, and protein network properties. Structural properties included the location at secondary structural elements, intrinsically disordered regions, transmembrane regions, repeats, conservation, and accessibility. Gene features were distributions of nucleotides, their groupings, codons, and location to CpG islands. The distributions of amino acids and their groups were investigated. VUSs did not markedly differ from other variants. The only major differences were the accessibility and conservation of pathogenic variants, and reduced ratio of repeat-locating variants in VUSs. Thus, all VUSs cannot be distinguished from other types of variants. They display one form of natural biological heterogeneity. Instead of concentrating on eradicating VUSs, the community would benefit from investigating and understanding factors that contribute to phenotypic heterogeneity.

Novel mechanisms of MITF regulation identified in a mouse suppressor screen

MITF, a basic Helix-Loop-Helix Zipper (bHLHZip) transcription factor, plays vital roles in melanocyte development and functions as an oncogene. We perform a genetic screen for suppressors of the Mitf-associated pigmentation phenotype in mice and identify an intragenic Mitf mutation that terminates MITF at the K316 SUMOylation site, leading to loss of the C-end intrinsically disordered region (IDR). The resulting protein is more nuclear but less stable than wild-type MITF and retains DNA-binding ability. As a dimer, it can translocate wild-type and mutant MITF partners into the nucleus, improving its own stability thus ensuring nuclear MITF supply. smFRET analysis shows interactions between K316 SUMOylation and S409 phosphorylation sites across monomers; these interactions largely explain the observed effects. The recurrent melanoma-associated E318K mutation in MITF, which affects K316 SUMOylation, also alters protein regulation in concert with S409. This suggests that residues K316 and S409 of MITF are impacted by SUMOylation and phosphorylation, respectively, mediating effects on nuclear localization and stability through conformational changes. Our work provides a novel mechanism of genetic suppression, and an example of how apparently deleterious mutations lead to normal phenotypes.

The Spectrum of Disease-Associated Alleles in Countries with a Predominantly Slavic Population

There are more than 260 million people of Slavic descent worldwide, who reside mainly in Eastern Europe but also represent a noticeable share of the population in the USA and Canada. Slavic populations, particularly Eastern Slavs and some Western Slavs, demonstrate a surprisingly high degree of genetic homogeneity, and, consequently, remarkable contribution of recurrent alleles associated with hereditary diseases. Along with pan-European pathogenic variants with clearly elevated occurrence in Slavic people (e.g., ATP7B c.3207C>A and PAH c.1222C>T), there are at least 52 pan-Slavic germ-line mutations (e.g., NBN c.657_661del and BRCA1 c.5266dupC) as well as several disease-predisposing alleles characteristic of the particular Slavic communities (e.g., Polish SDHD c.33C>A and Russian ARSB c.1562G>A variants). From a clinical standpoint, Slavs have some features of a huge founder population, thus providing a unique opportunity for efficient genetic studies.

Exome-wide evidence of compound heterozygous effects across common phenotypes in the UK Biobank

The phenotypic impact of compound heterozygous (CH) variation has not been investigated at the population scale. We phased rare variants (MAF ∼0.001%) in the UK Biobank (UKBB) exome-sequencing data to characterize recessive effects in 175,587 individuals across 311 common diseases. A total of 6.5% of individuals carry putatively damaging CH variants, 90% of which are only identifiable upon phasing rare variants (MAF < 0.38%). We identify six recessive gene-trait associations (p < 1.68 × 10-7) after accounting for relatedness, polygenicity, nearby common variants, and rare variant burden. Of these, just one is discovered when considering homozygosity alone. Using longitudinal health records, we additionally identify and replicate a novel association between bi-allelic variation in ATP2C2 and an earlier age at onset of chronic obstructive pulmonary disease (COPD) (p < 3.58 × 10-8). Genetic phase contributes to disease risk for gene-trait pairs: ATP2C2-COPD (p = 0.000238), FLG-asthma (p = 0.00205), and USH2A-visual impairment (p = 0.0084). We demonstrate the power of phasing large-scale genetic cohorts to discover phenome-wide consequences of compound heterozygosity.

A deep catalogue of protein-coding variation in 983,578 individuals

Rare coding variants that substantially affect function provide insights into the biology of a gene1-3. However, ascertaining the frequency of such variants requires large sample sizes4-8. Here we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. In total, 23% of the Regeneron Genetics Center Million Exome (RGC-ME) data come from individuals of African, East Asian, Indigenous American, Middle Eastern and South Asian ancestry. The catalogue includes more than 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss of function (LOF), we identify 3,988 LOF-intolerant genes, including 86 that were previously assessed as tolerant and 1,153 that lack established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions that are depleted of missense variants despite being tolerant of pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this resource of coding variation from the RGC-ME dataset publicly accessible through a variant allele frequency browser.

Measuring the Efficiency of Purging by non-random Mating in Human Populations

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.

Ancestry-specific high-risk gene variant profiling unmasks diabetes-associated genes

How ancestry-associated genetic variance affects disparities in the risk of polygenic diseases and influences the identification of disease-associated genes warrants a deeper understanding. We hypothesized that the discovery of genes associated with polygenic diseases may be limited by the overreliance on single-nucleotide polymorphism (SNP)-based genomic investigation, as most significant variants identified in genome-wide SNP association studies map to introns and intergenic regions of the genome. To overcome such potential limitations, we developed a gene-constrained, function-based analytical method centered on high-risk variants (hrV) that encode frameshifts, stopgains or splice site disruption. We analyzed the total number of hrV per gene in populations of different ancestry, representing a total of 185 934 subjects. Using this analysis, we developed a quantitative index of hrV (hrVI) across 20 428 genes within each population. We then applied hrVI analysis to the discovery of genes associated with type 2 diabetes mellitus (T2DM), a polygenic disease with ancestry-related disparity. HrVI profiling and gene-to-gene comparisons of ancestry-specific hrV between the case (20 781 subjects) and control (24 440 subjects) populations in the T2DM national repository identified 57 genes associated with T2DM, 40 of which were discoverable only by ancestry-specific analysis. These results illustrate how a function-based, ancestry-specific analysis of genetic variations can accelerate the identification of genes associated with polygenic diseases. Besides T2DM, such analysis may facilitate our understanding of the genetic basis for other polygenic diseases that are also greatly influenced by environmental and behavioral factors, such as obesity, hypertension and Alzheimer's disease.

Copy-number analysis from genome sequencing data of 11,754 rare-disease parent-child trios: A model for identifying autosomal recessive human gene knockouts including a novel gene for autosomal recessive retinopathy

Purpose

In parent-child trios with genome sequencing data, we investigated inherited biallelic deletions to identify known and novel genetic disorders.

Methods

We developed a copy-number variations analysis pipeline based on autosomal genome sequencing read depth of Genomics England 100,000 Genomes Project data from 11,754 parent-child trios and additional 18,875 non-trios. A control cohort of 15,440 cancer patients provided independent deletion frequencies.

Results

Autosomal recessive (AR) modeling detected 34 distinct rare deletions that were homozygous in the proband and heterozygous in both parents. These inherited biallelic deletions were only detected in 52 trios. These "knockout" regions included 37 genes, having among them 8 with an Online Mendelian Inheritance in Man AR annotation. Deletions of NPHP1, followed by OTOA, both within segmental duplications, were the only recurrent findings explaining phenotypes in a total of 10 and 3 patients, respectively. Recurrent heterozygous NPHP1 deletions were detected in 0.3%-0.5% of controls. We reviewed "knockout" patients for the remaining 29 genes without disease associations and identified SLC66A1 as a likely novel cause for AR rod-cone dystrophy in 4 families.

Conclusion

A tailored copy-number variations analysis of genome sequencing trio data shows that biallelic inherited gene deletions are rare, with NPHP1 biallelic deletions causing nephronophthisis the leading finding. We propose SLC66A1 as a novel cause for AR retinopathy.

Precision medicine using whole genome sequencing identifies a novel dystrophin (DMD) variant for X-linked muscular dystrophy in a cat

BACKGROUND

Muscular dystrophies (MDs) are a large, heterogeneous group of degenerative muscle diseases. X-linked dystrophin-deficient MD in cats is the first genetically characterized cat model for a human disease and a few novel forms have been identified.

HYPOTHESIS/OBJECTIVES

Muscular dystrophy was suspected in a young male domestic shorthair cat. Clinical, molecular, and genetic techniques could provide a definitive diagnosis.

ANIMALS

A 1-year-old male domestic shorthair cat presented for progressive difficulty walking, macroglossia and dysphagia beginning at 6 months of age. The tongue was thickened, protruded with constant ptyalism, and thickening and rigidity of the neck and shoulders were observed.

METHODS

A complete neurological examination, baseline laboratory evaluation and biopsies of the trapezius muscle were performed with owner consent. Indirect immunofluorescence staining of muscle cryosections was performed using several monoclonal and polyclonal antibodies against dystrophy-associated proteins. DNA was isolated for genomic analyses by whole genome sequencing and comparison to DNA variants in the 99 Lives Cat Genome Sequencing dataset.

RESULTS AND CLINICAL IMPORTANCE

Aspartate aminotransferase (687 IU/L) and creatine kinase (24 830 IU/L) activities were increased and mild hypokalemia (3.7 mmol/L) was present. Biopsy samples from the trapezius muscle confirmed a degenerative and regenerative myopathy and protein alterations identified by immunohistochemistry resulted in a diagnosis of a in dystrophin-deficient form of X-linked MD. A stop gain variant (c.4849C>T; p.Gln1617Ter) dystrophin was identified by genome sequencing. Precision/genomic medicine efforts for the domestic cat and in veterinary medicine support disease variant and animal model discovery and provide opportunities for targeted treatments for companion animals.

Development of a gene edited next-generation hematopoietic cell transplant to enable acute myeloid leukemia treatment by solving off-tumor toxicity

Immunotherapy of acute myeloid leukemia (AML) has been challenging because the lack of tumor-specific antigens results in "on-target, off-tumor" toxicity. To unlock the full potential of AML therapies, we used CRISPR-Cas9 to genetically ablate the myeloid protein CD33 from healthy donor hematopoietic stem and progenitor cells (HSPCs), creating tremtelectogene empogeditemcel (trem-cel). Trem-cel is a HSPC transplant product designed to provide a reconstituted hematopoietic compartment that is resistant to anti-CD33 drug cytotoxicity. Here, we describe preclinical studies and process development of clinical-scale manufacturing of trem-cel. Preclinical data showed proof-of-concept with loss of CD33 surface protein and no impact on myeloid cell differentiation or function. At clinical scale, trem-cel could be manufactured reproducibly, routinely achieving >70% CD33 editing with no effect on cell viability, differentiation, and function. Trem-cel pharmacology studies using mouse xenograft models showed long-term engraftment, multilineage differentiation, and persistence of gene editing. Toxicology assessment revealed no adverse findings, and no significant or reproducible off-target editing events. Importantly, CD33-knockout myeloid cells were resistant to the CD33-targeted agent gemtuzumab ozogamicin in vitro and in vivo. These studies supported the initiation of the first-in-human, multicenter clinical trial evaluating the safety and efficacy of trem-cel in patients with AML (NCT04849910).

A deep catalog of protein-coding variation in 985,830 individuals

Coding variants that have significant impact on function can provide insights into the biology of a gene but are typically rare in the population. Identifying and ascertaining the frequency of such rare variants requires very large sample sizes. Here, we present the largest catalog of human protein-coding variation to date, derived from exome sequencing of 985,830 individuals of diverse ancestry to serve as a rich resource for studying rare coding variants. Individuals of African, Admixed American, East Asian, Middle Eastern, and South Asian ancestry account for 20% of this Exome dataset. Our catalog of variants includes approximately 10.5 million missense (54% novel) and 1.1 million predicted loss-of-function (pLOF) variants (65% novel, 53% observed only once). We identified individuals with rare homozygous pLOF variants in 4,874 genes, and for 1,838 of these this work is the first to document at least one pLOF homozygote. Additional insights from the RGC-ME dataset include 1) improved estimates of selection against heterozygous loss-of-function and identification of 3,459 genes intolerant to loss-of-function, 83 of which were previously assessed as tolerant to loss-of-function and 1,241 that lack disease annotations; 2) identification of regions depleted of missense variation in 457 genes that are tolerant to loss-of-function; 3) functional interpretation for 10,708 variants of unknown or conflicting significance reported in ClinVar as cryptic splice sites using splicing score thresholds based on empirical variant deleteriousness scores derived from RGC-ME; and 4) an observation that approximately 3% of sequenced individuals carry a clinically actionable genetic variant in the ACMG SF 3.1 list of genes. We make this important resource of coding variation available to the public through a variant allele frequency browser. We anticipate that this report and the RGC-ME dataset will serve as a valuable reference for understanding rare coding variation and help advance precision medicine efforts.

TMCO3, a Putative K+ :Proton Antiporter at the Golgi Apparatus, Is Important for Longitudinal Growth in Mice and Humans

Isolated short stature, defined as short stature without any other abnormalities, is a common heterogeneous condition in children. Exome sequencing identified the homozygous nonsense variant c.1832G>A/p.(Trp611*) in TMCO3 in two sisters with isolated short stature. Radiological studies, biochemical measurements, assessment of the skeletal status, and three-dimensional bone microarchitecture revealed no relevant skeletal and bone abnormalities in both sisters. The homozygous TMCO3 variant segregated with short stature in the family. TMCO3 transcript levels were reduced by ~50% in leukocyte-derived RNA of both sisters compared with controls, likely due to nonsense-mediated mRNA decay. In primary urinary cells of heterozygous family members, we detected significantly reduced TMCO3 protein levels. TMCO3 is functionally uncharacterized. We ectopically expressed wild-type TMCO3 in HeLa and ATDC5 chondrogenic cells and detected TMCO3 predominantly at the Golgi apparatus, whereas the TMCO3W611* mutant did not reach the Golgi. Coordinated co-expression of TMCO3W611* -HA and EGFP in HeLa cells confirmed intrinsic instability and/or degradation of the mutant. Tmco3 is expressed in all relevant mouse skeletal cell types. Highest abundance of Tmco3 was found in chondrocytes of the prehypertrophic zone in mouse and minipig growth plates where it co-localizes with a Golgi marker. Knockdown of Tmco3 in differentiated ATDC5 cells caused reduced and increased expression of Pthlh and Ihh, respectively. Measurement of long bones in Tmco3tm1b(KOMP)Wtsi knockout mice revealed significant shortening of forelimbs and hindlimbs. TMCO3 is a potential member of the monovalent cation:proton antiporter 2 (CPA2) family. By in silico tools and homology modeling, TMCO3 is predicted to have an N-terminal secretory signal peptide, forms a dimer localized to the membrane, and is organized in a dimerization and a core domain. The core domain contains the CPA2 motif essential for K+ binding and selectivity. Collectively, our data demonstrate that loss of TMCO3 causes growth defects in both humans and mice. © 2023 American Society for Bone and Mineral Research (ASBMR).

Novel mechanisms of MITF regulation and melanoma predisposition identified in a mouse suppressor screen

MITF, a basic-Helix-Loop-Helix Zipper (bHLHZip) transcription factor, plays vital roles in melanocyte development and functions as an oncogene. To explore MITF regulation and its role in melanoma, we conducted a genetic screen for suppressors of the Mitf-associated pigmentation phenotype. An intragenic Mitf mutation was identified, leading to termination of MITF at the K316 SUMOylation site and loss of the C-end intrinsically disordered region (IDR). The resulting protein is more nuclear but less stable than wild-type MITF and retains DNA-binding ability. Interestingly, as a dimer, it can translocate wild-type and mutant MITF partners into the nucleus, improving its own stability and ensuring an active nuclear MITF supply. Interactions between K316 SUMOylation and S409 phosphorylation sites across monomers largely explain the observed effects. Notably, the recurrent melanoma-associated E318K mutation in MITF, which affects K316 SUMOylation, also alters protein regulation in concert with S409, unraveling a novel regulatory mechanism with unexpected disease insights.

Exome-wide evidence of compound heterozygous effects across common phenotypes in the UK Biobank

Exome-sequencing association studies have successfully linked rare protein-coding variation to risk of thousands of diseases. However, the relationship between rare deleterious compound heterozygous (CH) variation and their phenotypic impact has not been fully investigated. Here, we leverage advances in statistical phasing to accurately phase rare variants (MAF ~ 0.001%) in exome sequencing data from 175,587 UK Biobank (UKBB) participants, which we then systematically annotate to identify putatively deleterious CH coding variation. We show that 6.5% of individuals carry such damaging variants in the CH state, with 90% of variants occurring at MAF < 0.34%. Using a logistic mixed model framework, systematically accounting for relatedness, polygenic risk, nearby common variants, and rare variant burden, we investigate recessive effects in common complex diseases. We find six exome-wide significant (P < 1.68 × 10 - 7 ) and 17 nominally significant (P < 5.25 × 10 - 5 ) gene-trait associations. Among these, only four would have been identified without accounting for CH variation in the gene. We further incorporate age-at-diagnosis information from primary care electronic health records, to show that genetic phase influences lifetime risk of disease across 20 gene-trait combinations (FDR < 5%). Using a permutation approach, we find evidence for genetic phase contributing to disease susceptibility for a collection of gene-trait pairs, including FLG-asthma (P = 0.00205 ) and USH2A-visual impairment (P = 0.0084 ). Taken together, we demonstrate the utility of phasing large-scale genetic sequencing cohorts for robust identification of the phenome-wide consequences of compound heterozygosity.

Accurate rare variant phasing of whole-genome and whole-exome sequencing data in the UK Biobank

Phasing involves distinguishing the two parentally inherited copies of each chromosome into haplotypes. Here, we introduce SHAPEIT5, a new phasing method that quickly and accurately processes large sequencing datasets and applied it to UK Biobank (UKB) whole-genome and whole-exome sequencing data. We demonstrate that SHAPEIT5 phases rare variants with low switch error rates of below 5% for variants present in just 1 sample out of 100,000. Furthermore, we outline a method for phasing singletons, which, although less precise, constitutes an important step towards future developments. We then demonstrate that the use of UKB as a reference panel improves the accuracy of genotype imputation, which is even more pronounced when phased with SHAPEIT5 compared with other methods. Finally, we screen the UKB data for loss-of-function compound heterozygous events and identify 549 genes where both gene copies are knocked out. These genes complement current knowledge of gene essentiality in the human genome.

Tozorakimab (MEDI3506): an anti-IL-33 antibody that inhibits IL-33 signalling via ST2 and RAGE/EGFR to reduce inflammation and epithelial dysfunction

Interleukin (IL)-33 is a broad-acting alarmin cytokine that can drive inflammatory responses following tissue damage or infection and is a promising target for treatment of inflammatory disease. Here, we describe the identification of tozorakimab (MEDI3506), a potent, human anti-IL-33 monoclonal antibody, which can inhibit reduced IL-33 (IL-33red) and oxidized IL-33 (IL-33ox) activities through distinct serum-stimulated 2 (ST2) and receptor for advanced glycation end products/epidermal growth factor receptor (RAGE/EGFR complex) signalling pathways. We hypothesized that a therapeutic antibody would require an affinity higher than that of ST2 for IL-33, with an association rate greater than 107 M-1 s-1, to effectively neutralize IL-33 following rapid release from damaged tissue. An innovative antibody generation campaign identified tozorakimab, an antibody with a femtomolar affinity for IL-33red and a fast association rate (8.5 × 107 M-1 s-1), which was comparable to soluble ST2. Tozorakimab potently inhibited ST2-dependent inflammatory responses driven by IL-33 in primary human cells and in a murine model of lung epithelial injury. Additionally, tozorakimab prevented the oxidation of IL-33 and its activity via the RAGE/EGFR signalling pathway, thus increasing in vitro epithelial cell migration and repair. Tozorakimab is a novel therapeutic agent with a dual mechanism of action that blocks IL-33red and IL-33ox signalling, offering potential to reduce inflammation and epithelial dysfunction in human disease.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

Direct observation of the evolution of cell-type-specific microRNA expression signatures supports the hematopoietic origin model of endothelial cells

The evolution of specialized cell-types is a long-standing interest of biologists, but given the deep time-scales very difficult to reconstruct or observe. microRNAs have been linked to the evolution of cellular complexity and may inform on specialization. The endothelium is a vertebrate-specific specialization of the circulatory system that enabled a critical new level of vasoregulation. The evolutionary origin of these endothelial cells is unclear. We hypothesized that Mir-126, an endothelial cell-specific microRNA may be informative. We here reconstruct the evolutionary history of Mir-126. Mir-126 likely appeared in the last common ancestor of vertebrates and tunicates, which was a species without an endothelium, within an intron of the evolutionary much older EGF Like Domain Multiple (Egfl) locus. Mir-126 has a complex evolutionary history due to duplications and losses of both the host gene and the microRNA. Taking advantage of the strong evolutionary conservation of the microRNA among Olfactores, and using RNA in situ hybridization, we localized Mir-126 in the tunicate Ciona robusta. We found exclusive expression of the mature Mir-126 in granular amebocytes, supporting a long-proposed scenario that endothelial cells arose from hemoblasts, a type of proto-endothelial amoebocyte found throughout invertebrates. This observed change of expression of Mir-126 from proto-endothelial amoebocytes in the tunicate to endothelial cells in vertebrates is the first direct observation of the evolution of a cell-type in relation to microRNA expression indicating that microRNAs can be a prerequisite of cell-type evolution.

Ribonuclease inhibitor 1 (RNH1) deficiency cause congenital cataracts and global developmental delay with infection-induced psychomotor regression and anemia

Ribonuclease inhibitor 1, also known as angiogenin inhibitor 1, encoded by RNH1, is a ubiquitously expressed leucine-rich repeat protein, which is highly conserved in mammalian species. Inactivation of rnh1 in mice causes an embryonically lethal anemia, but the exact biological function of RNH1 in humans remains unknown and no human genetic disease has so far been associated with RNH1. Here, we describe a family with two out of seven siblings affected by a disease characterized by congenital cataract, global developmental delay, myopathy and psychomotor deterioration, seizures and periodic anemia associated with upper respiratory tract infections. A homozygous splice-site variant (c.615-2A > C) in RNH1 segregated with the disease. Sequencing of RNA derived from patient fibroblasts and cDNA analysis of skeletal muscle mRNA showed aberrant splicing with skipping of exon 7. Western blot analysis revealed a total lack of the RNH1 protein. Functional analysis revealed that patient fibroblasts were more sensitive to RNase A exposure, and this phenotype was reversed by transduction with a lentivirus expressing RNH1 to complement patient cells. Our results demonstrate that loss-of-function of RNH1 in humans is associated with a multiorgan developmental disease with recessive inheritance. It may be speculated that the infection-induced deterioration resulted from an increased susceptibility toward extracellular RNases and/or other inflammatory responses normally kept in place by the RNase inhibitor RNH1.

Molecular Genetic Characteristics of FANCI, a Proposed New Ovarian Cancer Predisposing Gene

FANCI was recently identified as a new candidate ovarian cancer (OC)-predisposing gene from the genetic analysis of carriers of FANCI c.1813C>T; p.L605F in OC families. Here, we aimed to investigate the molecular genetic characteristics of FANCI, as they have not been described in the context of cancer. We first investigated the germline genetic landscape of two sisters with OC from the discovery FANCI c.1813C>T; p.L605F family (F1528) to re-affirm the plausibility of this candidate. As we did not find other conclusive candidates, we then performed a candidate gene approach to identify other candidate variants in genes involved in the FANCI protein interactome in OC families negative for pathogenic variants in BRCA1, BRCA2, BRIP1, RAD51C, RAD51D, and FANCI, which identified four candidate variants. We then investigated FANCI in high-grade serous ovarian carcinoma (HGSC) from FANCI c.1813C>T carriers and found evidence of loss of the wild-type allele in tumour DNA from some of these cases. The somatic genetic landscape of OC tumours from FANCI c.1813C>T carriers was investigated for mutations in selected genes, copy number alterations, and mutational signatures, which determined that the profiles of tumours from carriers were characteristic of features exhibited by HGSC cases. As other OC-predisposing genes such as BRCA1 and BRCA2 are known to increase the risk of other cancers including breast cancer, we investigated the carrier frequency of germline FANCI c.1813C>T in various cancer types and found overall more carriers among cancer cases compared to cancer-free controls (p = 0.007). In these different tumour types, we also identified a spectrum of somatic variants in FANCI that were not restricted to any specific region within the gene. Collectively, these findings expand on the characteristics described for OC cases carrying FANCI c.1813C>T; p.L605F and suggest the possible involvement of FANCI in other cancer types at the germline and/or somatic level.

Systematic errors in annotations of truncations, loss-of-function and synonymous variants

Description of genetic phenomena and variations requires exact language and concepts. Vast amounts of variation data are produced with next-generation sequencing pipelines. The obtained variations are automatically annotated, e.g., for their functional consequences. These tools and pipelines, along with systematic nomenclature, mainly work well, but there are still some problems in nomenclature, organization of some databases, misuse of concepts and certain practices. Therefore, systematic errors prevent correct annotation and often preclude further analysis of certain variation types. Problems and solutions are described for presumed protein truncations, variants that are claimed to be of loss-of-function based on the type of variation, and synonymous variants that are not synonymous and lead to sequence changes or to missing protein.

Founder vs. non-founder BRCA1/2 pathogenic alleles: the analysis of Belarusian breast and ovarian cancer patients and review of other studies on ethnically homogenous populations

The spectrum of BRCA1/2 mutations demonstrates significant interethnic variations. We analyzed for the first time the entire BRCA1/2 coding region in 340 Belarusian cancer patients with clinical signs of BRCA1/2-related disease, including 168 women with bilateral and/or early-onset breast cancer (BC), 104 patients with ovarian cancer and 68 subjects with multiple primary malignancies involving BC and/or OC. BRCA1/2 pathogenic alleles were detected in 98 (29%) women, with 67 (68%) of these being represented by founder alleles. Systematic comparison with other relevant studies revealed that the founder effect observed in Belarus is among the highest estimates observed worldwide. These findings are surprising, given that the population of Belarus did not experience geographic or cultural isolation throughout history.

Diverse monogenic subforms of human spermatogenic failure

Non-obstructive azoospermia (NOA) is the most severe form of male infertility and typically incurable. Defining the genetic basis of NOA has proven challenging, and the most advanced classification of NOA subforms is not based on genetics, but simple description of testis histology. In this study, we exome-sequenced over 1000 clinically diagnosed NOA cases and identified a plausible recessive Mendelian cause in 20%. We find further support for 21 genes in a 2-stage burden test with 2072 cases and 11,587 fertile controls. The disrupted genes are primarily on the autosomes, enriched for undescribed human "knockouts", and, for the most part, have yet to be linked to a Mendelian trait. Integration with single-cell RNA sequencing data shows that azoospermia genes can be grouped into molecular subforms with synchronized expression patterns, and analogs of these subforms exist in mice. This analysis framework identifies groups of genes with known roles in spermatogenesis but also reveals unrecognized subforms, such as a set of genes expressed across mitotic divisions of differentiating spermatogonia. Our findings highlight NOA as an understudied Mendelian disorder and provide a conceptual structure for organizing the complex genetics of male infertility, which may provide a rational basis for disease classification.

Identification of PCSK9-like human gene knockouts using metabolomics, proteomics, and whole-genome sequencing in a consanguineous population

Natural human knockouts of genes associated with desirable outcomes, such as PCSK9 with low levels of LDL-cholesterol, can lead to the discovery of new drug targets and treatments. Rare loss-of-function variants are more likely to be found in the homozygous state in consanguineous populations, and deep molecular phenotyping of blood samples from homozygous carriers can help to discriminate between silent and functional variants. Here, we combined whole-genome sequencing with proteomics and metabolomics for 2,935 individuals from the Qatar Biobank (QBB) to evaluate the power of this approach for finding genes of clinical and pharmaceutical interest. As proof-of-concept, we identified a homozygous carrier of a very rare PCSK9 variant with extremely low circulating PCSK9 levels and low LDL. Our study demonstrates that the chances of finding such variants are about 168 times higher in QBB compared with GnomAD and emphasizes the potential of consanguineous populations for drug discovery.

Developmental disruption and restoration of brain synaptome architecture in the murine Pax6 neurodevelopmental disease model

Neurodevelopmental disorders of genetic origin delay the acquisition of normal abilities and cause disabling phenotypes. Nevertheless, spontaneous attenuation and even complete amelioration of symptoms in early childhood and adolescence can occur in many disorders, suggesting that brain circuits possess an intrinsic capacity to overcome the deficits arising from some germline mutations. We examined the molecular composition of almost a trillion excitatory synapses on a brain-wide scale between birth and adulthood in mice carrying a mutation in the homeobox transcription factor Pax6, a neurodevelopmental disorder model. Pax6 haploinsufficiency had no impact on total synapse number at any age. By contrast, the molecular composition of excitatory synapses, the postnatal expansion of synapse diversity and the acquisition of normal synaptome architecture were delayed in all brain regions, interfering with networks and electrophysiological simulations of cognitive functions. Specific excitatory synapse types and subtypes were affected in two key developmental age-windows. These phenotypes were reversed within 2-3 weeks of onset, restoring synapse diversity and synaptome architecture to the normal developmental trajectory. Synapse subtypes with rapid protein turnover mediated the synaptome remodeling. This brain-wide capacity for remodeling of synapse molecular composition to recover and maintain the developmental trajectory of synaptome architecture may help confer resilience to neurodevelopmental genetic disorders.

Neutrophil extracellular traps in systemic autoimmune and autoinflammatory diseases

Systemic autoimmune diseases are characterized by the failure of the immune system to differentiate self from non-self. These conditions are associated with significant morbidity and mortality, and they can affect many organs and systems, having significant clinical heterogeneity. Recent discoveries have highlighted that neutrophils, and in particular the neutrophil extracellular traps that they can release upon activation, can have central roles in the initiation and perpetuation of systemic autoimmune disorders and orchestrate complex inflammatory responses that lead to organ damage. Dysregulation of neutrophil cell death can lead to the modification of autoantigens and their presentation to the adaptive immune system. Furthermore, subsets of neutrophils that seem to be more prevalent in patients with systemic autoimmune disorders can promote vascular damage and increased oxidative stress. With the emergence of new technologies allowing for improved assessments of neutrophils, the complexity of neutrophil biology and its dysregulation is now starting to be understood. In this Review, we provide an overview of the roles of neutrophils in systemic autoimmune and autoinflammatory diseases and address putative therapeutic targets that may be explored based on this new knowledge.

Incomplete Penetrance and Variable Expressivity: From Clinical Studies to Population Cohorts

The same genetic variant found in different individuals can cause a range of diverse phenotypes, from no discernible clinical phenotype to severe disease, even among related individuals. Such variants can be said to display incomplete penetrance, a binary phenomenon where the genotype either causes the expected clinical phenotype or it does not, or they can be said to display variable expressivity, in which the same genotype can cause a wide range of clinical symptoms across a spectrum. Both incomplete penetrance and variable expressivity are thought to be caused by a range of factors, including common variants, variants in regulatory regions, epigenetics, environmental factors, and lifestyle. Many thousands of genetic variants have been identified as the cause of monogenic disorders, mostly determined through small clinical studies, and thus, the penetrance and expressivity of these variants may be overestimated when compared to their effect on the general population. With the wealth of population cohort data currently available, the penetrance and expressivity of such genetic variants can be investigated across a much wider contingent, potentially helping to reclassify variants that were previously thought to be completely penetrant. Research into the penetrance and expressivity of such genetic variants is important for clinical classification, both for determining causative mechanisms of disease in the affected population and for providing accurate risk information through genetic counseling. A genotype-based definition of the causes of rare diseases incorporating information from population cohorts and clinical studies is critical for our understanding of incomplete penetrance and variable expressivity. This review examines our current knowledge of the penetrance and expressivity of genetic variants in rare disease and across populations, as well as looking into the potential causes of the variation seen, including genetic modifiers, mosaicism, and polygenic factors, among others. We also considered the challenges that come with investigating penetrance and expressivity.

Rare and population-specific functional variation across pig lines

BACKGROUND

It is expected that functional, mainly missense and loss-of-function (LOF), and regulatory variants are responsible for most phenotypic differences between breeds and genetic lines of livestock species that have undergone diverse selection histories. However, there is still limited knowledge about the existing missense and LOF variation in commercial livestock populations, in particular regarding population-specific variation and how it can affect applications such as across-breed genomic prediction.

METHODS

We re-sequenced the whole genome of 7848 individuals from nine commercial pig lines (average sequencing coverage: 4.1×) and imputed whole-genome genotypes for 440,610 pedigree-related individuals. The called variants were categorized according to predicted functional annotation (from LOF to intergenic) and prevalence level (number of lines in which the variant segregated; from private to widespread). Variants in each category were examined in terms of their distribution along the genome, alternative allele frequency, per-site Wright's fixation index (F_ST), individual load, and association to production traits.

RESULTS

Of the 46 million called variants, 28% were private (called in only one line) and 21% were widespread (called in all nine lines). Genomic regions with a low recombination rate were enriched with private variants. Low-prevalence variants (called in one or a few lines only) were enriched for lower allele frequencies, lower F_ST, and putatively functional and regulatory roles (including LOF and deleterious missense variants). On average, individuals carried fewer private deleterious missense alleles than expected compared to alleles with other predicted consequences. Only a small subset of the low-prevalence variants had intermediate allele frequencies and explained small fractions of phenotypic variance (up to 3.2%) of production traits. The significant low-prevalence variants had higher per-site F_ST than the non-significant ones. These associated low-prevalence variants were tagged by other more widespread variants in high linkage disequilibrium, including intergenic variants.

CONCLUSIONS

Most low-prevalence variants have low minor allele frequencies and only a small subset of low-prevalence variants contributed detectable fractions of phenotypic variance of production traits. Accounting for low-prevalence variants is therefore unlikely to noticeably benefit across-breed analyses, such as the prediction of genomic breeding values in a population using reference populations of a different genetic background.

Beyond Mendelian Inheritance: Genetic Buffering and Phenotype Variability

Understanding the way genes work amongst individuals and across generations to shape form and function is a common theme for many genetic studies. The recent advances in genetics, genome engineering and DNA sequencing reinforced the notion that genes are not the only players that determine a phenotype. Due to physiological or pathological fluctuations in gene expression, even genetically identical cells can behave and manifest different phenotypes under the same conditions. Here, we discuss mechanisms that can influence or even disrupt the axis between genotype and phenotype; the role of modifier genes, the general concept of genetic redundancy, genetic compensation, the recently described transcriptional adaptation, environmental stressors, and phenotypic plasticity. We furthermore highlight the usage of induced pluripotent stem cells (iPSCs), the generation of isogenic lines through genome engineering, and sequencing technologies can help extract new genetic and epigenetic mechanisms from what is hitherto considered 'noise'.

Population-level deficit of homozygosity unveils CPSF3 as an intellectual disability syndrome gene

Predicting the pathogenicity of biallelic missense variants can be challenging. Here, we use a deficit of observed homozygous carriers of missense variants, versus an expected number in a set of 153,054 chip-genotyped Icelanders, to identify potentially pathogenic genotypes. We follow three missense variants with a complete deficit of homozygosity and find that their pathogenic effect in homozygous state ranges from severe childhood disease to early embryonic lethality. One of these variants is in CPSF3, a gene not previously linked to disease. From a set of clinically sequenced Icelanders, and by sequencing archival samples targeted through the Icelandic genealogy, we find four homozygous carriers. Additionally, we find two homozygous carriers of Mexican descent of another missense variant in CPSF3. All six homozygous carriers of missense variants in CPSF3 show severe intellectual disability, seizures, microcephaly, and abnormal muscle tone. Here, we show how the absence of certain homozygous genotypes from a large population set can elucidate causes of previously unexplained recessive diseases and early miscarriage.

Characterising a homozygous two-exon deletion in UQCRH: comparing human and mouse phenotypes

Mitochondrial disorders are clinically and genetically diverse, with isolated complex III (CIII) deficiency being relatively rare. Here, we describe two affected cousins, presenting with recurrent episodes of severe lactic acidosis, hyperammonaemia, hypoglycaemia and encephalopathy. Genetic investigations in both cases identified a homozygous deletion of exons 2 and 3 of UQCRH, which encodes a structural complex III (CIII) subunit. We generated a mouse model with the equivalent homozygous Uqcrh deletion (Uqcrh-/- ), which also presented with lactic acidosis and hyperammonaemia, but had a more severe, non-episodic phenotype, resulting in failure to thrive and early death. The biochemical phenotypes observed in patient and Uqcrh-/- mouse tissues were remarkably similar, displaying impaired CIII activity, decreased molecular weight of fully assembled holoenzyme and an increase of an unexpected large supercomplex (SXL ), comprising mostly of one complex I (CI) dimer and one CIII dimer. This phenotypic similarity along with lentiviral rescue experiments in patient fibroblasts verifies the pathogenicity of the shared genetic defect, demonstrating that the Uqcrh-/- mouse is a valuable model for future studies of human CIII deficiency.

The genetic structure of the Turkish population reveals high levels of variation and admixture

The construction of population-based variomes has contributed substantially to our understanding of the genetic basis of human inherited disease. Here, we investigated the genetic structure of Turkey from 3,362 unrelated subjects whose whole exomes (n = 2,589) or whole genomes (n = 773) were sequenced to generate a Turkish (TR) Variome that should serve to facilitate disease gene discovery in Turkey. Consistent with the history of present-day Turkey as a crossroads between Europe and Asia, we found extensive admixture between Balkan, Caucasus, Middle Eastern, and European populations with a closer genetic relationship of the TR population to Europeans than hitherto appreciated. We determined that 50% of TR individuals had high inbreeding coefficients (≥0.0156) with runs of homozygosity longer than 4 Mb being found exclusively in the TR population when compared to 1000 Genomes Project populations. We also found that 28% of exome and 49% of genome variants in the very rare range (allele frequency < 0.005) are unique to the modern TR population. We annotated these variants based on their functional consequences to establish a TR Variome containing alleles of potential medical relevance, a repository of homozygous loss-of-function variants and a TR reference panel for genotype imputation using high-quality haplotypes, to facilitate genome-wide association studies. In addition to providing information on the genetic structure of the modern TR population, these data provide an invaluable resource for future studies to identify variants that are associated with specific phenotypes as well as establishing the phenotypic consequences of mutations in specific genes.

Natural human knockouts and Mendelian disorders: deep phenotyping in Italian isolates

Whole genome sequencing (WGS) allows the identification of human knockouts (HKOs), individuals in whom loss of function (LoF) variants disrupt both alleles of a given gene. HKOs are a valuable model for understanding the consequences of genes function loss. Naturally occurring biallelic LoF variants tend to be significantly enriched in "genetic isolates," making these populations specifically suited for HKO studies. In this work, a meticulous WGS data analysis combined with an in-depth phenotypic assessment of 947 individuals from three Italian genetic isolates led to the identification of ten biallelic LoF variants in ten OMIM genes associated with known autosomal recessive diseases. Notably, only a minority of the identified HKOs (C7, F12, and GPR68 genes) displayed the expected phenotype. For most of the genes, instead, (ACADSB, FANCL, GRK1, LGI4, MPO, PGAM2, and RP1L1), the carriers showed none or few of the signs and symptoms typically associated with the related diseases. Of particular interest is a case presenting with a FANCL biallelic LoF variant and a positive diepoxybutane test but lacking a full Fanconi anemia phenotypic spectrum. Identifying KO subjects displaying expected phenotypes suggests that the lack of correct genetic diagnoses may lead to inappropriate and delayed treatment. In contrast, the presence of HKOs with phenotypes deviating from the expected patterns underlines how LoF variants may be responsible for broader phenotypic spectra. Overall, these results highlight the importance of in-depth phenotypical characterization to understand the role of LoF variants and the advantage of studying these variants in genetic isolates.

The genetic architecture of age-related hearing impairment revealed by genome-wide association analysis

Age-related hearing impairment (ARHI) is the most common sensory disorder in older adults. We conducted a genome-wide association meta-analysis of 121,934 ARHI cases and 591,699 controls from Iceland and the UK. We identified 21 novel sequence variants, of which 13 are rare, under either additive or recessive models. Of special interest are a missense variant in LOXHD1 (MAF = 1.96%) and a tandem duplication in FBF1 covering 4 exons (MAF = 0.22%) associating with ARHI (OR = 3.7 for homozygotes, P = 1.7 × 10-22 and OR = 4.2 for heterozygotes, P = 5.7 × 10-27, respectively). We constructed an ARHI genetic risk score (GRS) using common variants and showed that a common variant GRS can identify individuals at risk comparable to carriers of rare high penetrance variants. Furthermore, we found that ARHI and tinnitus share genetic causes. This study sheds a new light on the genetic architecture of ARHI, through several rare variants in both Mendelian deafness genes and genes not previously linked to hearing.

A data set of variants derived from 1455 clinical and research exomes is efficient in variant prioritization for early-onset monogenic disorders in Indians

Given the genomic uniqueness, a local data set is most desired for Indians, who are underrepresented in existing public databases. We hypothesize patients with rare monogenic disorders and their family members can provide a reliable source of common variants in the population. Exome sequencing (ES) data from families with rare Mendelian disorders was aggregated from five centers in India. The dataset was refined by excluding related individuals and removing the disease-causing variants (refined cohort). The efficiency of these data sets was assessed in a new set of 50 exomes against gnomAD and GenomeAsia. Our original cohort comprised 1455 individuals from 1203 families. The refined cohort had 836 unrelated individuals that retained 1,251,064 variants with 181,125 population-specific and 489,618 common variants. The allele frequencies from our cohort helped to define 97,609 rare variants in gnomAD and 44,520 rare variants in GenomeAsia as common variants in our population. Our variant dataset provided an additional 1.7% and 0.1% efficiency for prioritizing heterozygous and homozygous variants respectively for rare monogenic disorders. We observed additional 19 genes/human knockouts. We list carrier frequency for 142 recessive disorders. This is a large and useful resource of exonic variants for Indians. Despite limitations, datasets from patients are efficient tools for variant prioritization in a resource-limited setting.

Bite of the wolf: innate immune responses propagate autoimmunity in lupus

The etiopathogenesis of systemic lupus erythematosus (SLE), a clinically heterogeneous multisystemic syndrome that derives its name from the initial characterization of facial lesions that resemble the bite of a wolf, is considered a complex, multifactorial interplay between underlying genetic susceptibility factors and the environment. Prominent pathogenic factors include the induction of aberrant cell death pathways coupled with defective cell death clearance mechanisms that promote excessive externalization of modified cellular and nuclear debris with subsequent loss of tolerance to a wide variety of autoantigens and innate and adaptive immune dysregulation. While abnormalities in adaptive immunity are well recognized and are key to the pathogenesis of SLE, recent findings have emphasized fundamental roles of the innate immune system in the initiation and propagation of autoimmunity and the development of organ damage in this disease. This Review focuses on recent discoveries regarding the role of components of the innate immune system, specifically neutrophils and interferons, in promoting various aspects of lupus pathogenesis, with potential implications for novel therapeutic strategies.

Plasticity of Carbohydrate Transport at the Blood-Brain Barrier

Neuronal function is highly energy demanding, requiring efficient transport of nutrients into the central nervous system (CNS). Simultaneously the brain must be protected from the influx of unwanted solutes. Most of the energy is supplied from dietary sugars, delivered from circulation via the blood-brain barrier (BBB). Therefore, selective transporters are required to shuttle metabolites into the nervous system where they can be utilized. The Drosophila BBB is formed by perineural and subperineurial glial cells, which effectively separate the brain from the surrounding hemolymph, maintaining a constant microenvironment. We identified two previously unknown BBB transporters, MFS3 (Major Facilitator Superfamily Transporter 3), located in the perineurial glial cells, and Pippin, found in both the perineurial and subperineurial glial cells. Both transporters facilitate uptake of circulating trehalose and glucose into the BBB-forming glial cells. RNA interference-mediated knockdown of these transporters leads to pupal lethality. However, null mutants reach adulthood, although they do show reduced lifespan and activity. Here, we report that both carbohydrate transport efficiency and resulting lethality found upon loss of MFS3 or Pippin are rescued via compensatory upregulation of Tret1-1, another BBB carbohydrate transporter, in Mfs3 and pippin null mutants, while RNAi-mediated knockdown is not compensated for. This means that the compensatory mechanisms in place upon mRNA degradation following RNA interference can be vastly different from those resulting from a null mutation.

How wide is the application of genetic big data in biomedicine

In the era of big data, massive genetic data, as a new industry, has quickly swept almost all industries, especially the pharmaceutical industry. As countries around the world start to build their own gene banks, scientists study the data to explore the origins and migration of humans. Moreover, big data encourage the development of cancer therapy and bring good news to cancer patients. Big datum has been involved in the study of many diseases, and it has been found that analyzing diseases at the gene level can lead to more beneficial treatment options than ordinary treatments. This review will introduce the development of extensive data in medical research from the perspective of big data and tumor, neurological and psychiatric diseases, cardiovascular diseases, other applications and the development direction of big data in medicine.

Duchenne Muscular Dystrophy and Early Onset Hypertrophic Cardiomyopathy associated with Mutations in Dystrophin and Hypertrophic Cardiomyopathy-Associated Genes

Duchenne muscular dystrophy (DMD) is a progressive muscular damage disorder caused by mutations in dystrophin gene. Cardiomyopathy may first be evident after 10 years of age and increases in incidence with age. We present a boy diagnosed at 18 months with a rare phenotype of DMD in association with early-onset hypertrophic cardiomyopathy (HCM). The cause of DMD is a deletion of exons 51–54 of dystrophin gene. The cause of HCM was verified by whole exome sequencing. Novel missense variations in two genes: MAP2K5 inherited from the mother and ACTN2 inherited from the father, or de novo. The combination of MAP2K5, ACTN2, and dystrophin mutations, could be causing the HCM in our patient. This is the second patient diagnosed, at relatively young age, with DMD and HCM, with novel variations in genes known to cause HCM. This study demonstrates the need for genetic diagnosis to elucidate the underlying pathology of HCM.

Less Is More, Natural Loss-of-Function Mutation Is a Strategy for Adaptation

Gene gain and loss are crucial factors that shape the evolutionary success of diverse organisms. In the past two decades, more attention has been paid to the significance of gene gain through gene duplication or de novo genes. However, gene loss through natural loss-of-function (LoF) mutations, which is prevalent in the genomes of diverse organisms, has been largely ignored. With the development of sequencing techniques, many genomes have been sequenced across diverse species and can be used to study the evolutionary patterns of gene loss. In this review, we summarize recent advances in research on various aspects of LoF mutations, including their identification, evolutionary dynamics in natural populations, and functional effects. In particular, we discuss how LoF mutations can provide insights into the minimum gene set (or the essential gene set) of an organism. Furthermore, we emphasize their potential impact on adaptation. At the genome level, although most LoF mutations are neutral or deleterious, at least some of them are under positive selection and may contribute to biodiversity and adaptation. Overall, we highlight the importance of natural LoF mutations as a robust framework for understanding biological questions in general.

Genetic compensation prevents myopathy and heart failure in an in vivo model of Bag3 deficiency

Mutations in the molecular co-chaperone Bcl2-associated athanogene 3 (BAG3) are found to cause dilated cardiomyopathy (DCM), resulting in systolic dysfunction and heart failure, as well as myofibrillar myopathy (MFM), which is characterized by protein aggregation and myofibrillar disintegration in skeletal muscle cells. Here, we generated a CRISPR/Cas9-induced Bag3 knockout zebrafish line and found the complete preservation of heart and skeletal muscle structure and function during embryonic development, in contrast to morpholino-mediated knockdown of Bag3. Intriguingly, genetic compensation, a process of transcriptional adaptation which acts independent of protein feedback loops, was found to prevent heart and skeletal muscle damage in our Bag3 knockout model. Proteomic profiling and quantitative real-time PCR analyses identified Bag2, another member of the Bag protein family, significantly upregulated on a transcript and protein level in bag3^-/- mutants. This implied that the decay of bag3 mutant mRNA in homozygous bag3^-/- embryos caused the transcriptional upregulation of bag2 expression. We further demonstrated that morpholino-mediated knockdown of Bag2 in bag3^-/- embryos evoked severe functional and structural heart and skeletal muscle defects, which are similar to Bag3 morphants. However, Bag2 knockdown in bag3^+/+ or bag3^+/- embryos did not result in (cardio-)myopathy. Finally, we found that inhibition of the nonsense-mediated mRNA decay (NMD) machinery by knockdown of upf1, an essential NMD factor, caused severe heart and skeletal muscle defects in bag3^-/- mutants due to the blockade of transcriptional adaptation of bag2 expression. Our findings provide evidence that genetic compensation might vitally influence the penetrance of disease-causing bag3 mutations in vivo.

One form of genetic compensation is described as transcriptional adaptation of gene expression triggered by deleterious gene mutations. Although the precise molecular mechanism that induces genetic compensation needs to be defined, it represents a powerful biological phenomenon that warrants genetic robustness. We find that antisense-mediated knockdown of Bag3 in zebrafish embryos causes heart failure and myopathy. By contrast, CRISPR/Cas9-induced depletion of Bag3 does not result in the abrogation of heart and skeletal muscle function in zebrafish embryos. We find here that transcriptional activation of the Bag family member bag2 is capable of restoring heart and skeletal muscle function in bag3 mutant embryos, whereas this compensatory mechanism is not present in the bag3 morphants. Furthermore, we show that nonsense-mediated decay of bag3 mRNA is the molecular trigger for the compensatory upregulation of bag2. Our study provides evidence that genetic compensation via transcriptional adaptation is a vital modulator of disease peculiarity and penetrance in bag3 mutant zebrafish and that this biological phenomenon might also be active in certain human BAG3 mutation carriers.

Genetic background and transplantation outcomes: insights from genome-wide association studies

PURPOSE OF REVIEW

The current review summarizes recent advances in the genetic studies of transplantation outcomes, including new genome-wide association studies for acute rejection, allograft survival, pharmacogenomics, and common transplant comorbidities.

RECENT FINDINGS

Genetic studies of kidney transplantation outcomes have begun to address the question of genetic compatibility beyond human leukocyte antigens, including the role of genome-wide mismatches in missense variants, and the 'genomic collision' hypothesis under which the risk of rejection may be increased in recipients homozygous for loss-of-function variants with grafts from nonhomozygous donors. In recent pilot studies, missense mismatch scores for transmembrane and secreted proteins were associated with antibodies against the mismatched peptides and reduced allograft survival. A 'genomic collision' at the LIMS1 locus involving a common deletion near LIMS1 gene was associated with anti-LIMS1 antibody response and increased risk of rejection. Additional genetic factors under active investigation include genome-wide polygenic risk scores for renal function and apolipoprotein L1 risk genotypes in African-American kidney donors. Due to the heterogeneity and complexity of clinical outcomes, new genome-wide association studies for rejection, allograft survival, and specific transplant comorbidities will require larger multicenter meta-analyses.

SUMMARY

Genetic compatibilities between donor and recipient represent an important determinant of rejection and long-term allograft survival. Genetic background of transplant donors may be additionally predictive of allograft function, while recipient's genomes are likely determinant of a wide range of transplantation outcomes, from rejection susceptibility to pharmacogenetics and various comorbidities related to prolonged immunosuppression.

Functional effects of protein variants

Genetic and other variations frequently affect protein functions. Scientific articles can contain confusing descriptions about which function or property is affected, and in many cases the statements are pure speculation without any experimental evidence. To clarify functional effects of protein variations of genetic or non-genetic origin, a systematic conceptualisation and framework are introduced. This framework describes protein functional effects on abundance, activity, specificity and affinity, along with countermeasures, which allow cells, tissues and organisms to tolerate, avoid, repair, attenuate or resist (TARAR) the effects. Effects on abundance discussed include gene dosage, restricted expression, mis-localisation and degradation. Enzymopathies, effects on kinetics, allostery and regulation of protein activity are subtopics for the effects of variants on activity. Variation outcomes on specificity and affinity comprise promiscuity, specificity, affinity and moonlighting. TARAR mechanisms redress variations with active and passive processes including chaperones, redundancy, robustness, canalisation and metabolic and signalling rewiring. A framework for pragmatic protein function analysis and presentation is introduced. All of the mechanisms and effects are described along with representative examples, most often in relation to diseases. In addition, protein function is discussed from evolutionary point of view. Application of the presented framework facilitates unambiguous, detailed and specific description of functional effects and their systematic study.

Sequence Variants in TAAR5 and Other Loci Affect Human Odor Perception and Naming

Olfactory receptor (OR) genes in humans form a special class characterized by unusually high DNA sequence diversity, which should give rise to differences in perception and behavior. In the largest genome-wide association study to date based on olfactory testing, we investigated odor perception and naming with smell tasks performed by 9,122 Icelanders, with replication in a separate sample of 2,204 individuals. We discovered an association between a low-frequency missense variant in TAAR5 and reduced intensity rating of fish odor containing trimethylamine (p.Ser95Pro, p_combined = 5.6 × 10^-15). We demonstrate that TAAR5 genotype affects aversion to fish odor, reflected by linguistic descriptions of the odor and pleasantness ratings. We also discovered common sequence variants in two canonical olfactory receptor loci that associate with increased intensity and naming of licorice odor (trans-anethole: lead variant p.Lys233Asn in OR6C70, p_combined = 8.8 × 10^-16 and p_combined = 1.4 × 10^-9) and enhanced naming of cinnamon (trans-cinnamaldehyde; intergenic variant rs317787-T, p_combined = 5.0 × 10^-17). Together, our results show that TAAR5 genotype variation influences human odor responses and highlight that sequence diversity in canonical OR genes can lead to enhanced olfactory ability, in contrast to the view that greater tolerance for mutations in the human OR repertoire leads to diminished function.

Transcriptional adaptation: a mechanism underlying genetic robustness

Mutations play a crucial role in evolution as they provide the genetic variation that allows evolutionary change. Although some mutations in regulatory elements or coding regions can be beneficial, a large number of them disrupt gene function and reduce fitness. Organisms utilize several mechanisms to compensate for the damaging consequences of genetic perturbations. One such mechanism is the recently identified process of transcriptional adaptation (TA): during this event, mutations that cause mutant mRNA degradation trigger the transcriptional modulation of so-called adapting genes. In some cases, for example when one (or more) of the upregulated genes is functionally redundant with the mutated gene, this process compensates for the loss of the mutated gene's product. Notably, unlike other mechanisms underlying genetic robustness, TA is not triggered by the loss of protein function, an observation that has prompted studies into the machinery of TA and the contexts in which it functions. Here, we review the discovery and current understanding of TA, and discuss how its main features appear to be conserved across species. In light of these findings, we also speculate on the importance of TA in the context of human disease, and provide some recommendations for genome-editing strategies that should be more effective.