Human knockout research: new horizons and opportunities

NanoRanger enables rapid single-base-pair resolution of genomic disorders

BACKGROUND

Delineating base-resolution breakpoints of complex rearrangements is crucial for an accurate clinical understanding of pathogenic variants and for carrier screening within family networks or the broader population. However, despite advances in genetic testing using short-read sequencing (SRS), this task remains costly and challenging.

METHODS

This study addresses the challenges of resolving missing disease-causing breakpoints in complex genomic disorders with suspected homozygous rearrangements by employing multiple long-read sequencing (LRS) strategies, including a novel and efficient strategy named nanopore-based rapid acquisition of neighboring genomic regions (NanoRanger). NanoRanger does not require large amounts of ultrahigh-molecular-weight DNA and stands out for its ease of use and rapid acquisition of large genomic regions of interest with deep coverage.

FINDINGS

We describe a cohort of 16 familial cases, each harboring homozygous rearrangements that defied breakpoint determination by SRS and optical genome mapping (OGM). NanoRanger identified the breakpoints with single-base-pair resolution, enabling accurate determination of the carrier status of unaffected family members as well as the founder nature of these genomic lesions and their frequency in the local population. The resolved breakpoints revealed that repetitive DNA, gene regulatory elements, and transcription activity contribute to genome instability in these novel recessive rearrangements.

CONCLUSIONS

Our data suggest that NanoRanger greatly improves the success rate of resolving base-resolution breakpoints of complex genomic disorders and expands access to LRS for the benefit of patients with Mendelian disorders.

FUNDING

M.L. is supported by KAUST Baseline Award no. BAS/1/1080-01-01 and KAUST Research Translation Fund Award no. REI/1/4742-01.

Identification of novel genetic variations in ABCB6 and GRN genes associated with HIV-associated lipodystrophy

Protease inhibitors (PIs) are associated with an incidence of lipodystrophy among people living with HIV(PLHIV). Lipodystrophiesare characterised by the loss of adipose tissue. Evidence suggests that a patient's lipodystrophy phenotype is influenced by genetic mutation, age, gender, and environmental and genetic factors, such as single-nucleotide variants (SNVs). Pathogenic variants are considered to cause a more significant loss of adipose tissue compared to non-pathogenic. Lipid metabolising enzymes and transporter genes have a role in regulating lipoprotein metabolism and have been associated with lipodystrophy in HIV-infected patients (LDHIV). The long-term effect of the lipodystrophy syndrome is related to cardiovascular diseases (CVDs). Hence, we determined the SNVs of lipid metabolising enzymes and transporter genes in a total of 48 patient samples, of which 24 were with and 24 were without HIV-associated lipodystrophy (HIVLD) using next-generation sequencing. A panel of lipid metabolism, transport and elimination genes were sequenced. Three novel heterozygous non-synonymous variants at exon 8 (c.C1400A:p.S467Y, c.G1385A:p.G462E, and c.T1339C:p.S447P) in the ABCB6 gene were identified in patients with lipodystrophy. One homozygous non-synonymous SNV (exon5:c.T358C:p.S120P) in the GRN gene was identified in patients with lipodystrophy. One novelstop-gain SNV (exon5:c.C373T:p.Q125X) was found in the GRN gene among patients without lipodystrophy. Patients without lipodystrophy had one homozygous non-synonymous SNV (exon9:c.G1462T:p.G488C) in the ABCB6 gene. Our findings suggest that novel heterozygous non-synonymous variants in the ABCB6 gene may contribute to defective protein production, potentially intensifying the severity of lipodystrophy. Additionally, identifying a stop-gain SNV in the GRN gene among patients without lipodystrophy implies a potential role in the development of HIVLD.

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.

Diagnostic implications of pitfalls in causal variant identification based on 4577 molecularly characterized families

Despite large sequencing and data sharing efforts, previously characterized pathogenic variants only account for a fraction of Mendelian disease patients, which highlights the need for accurate identification and interpretation of novel variants. In a large Mendelian cohort of 4577 molecularly characterized families, numerous scenarios in which variant identification and interpretation can be challenging are encountered. We describe categories of challenges that cover the phenotype (e.g. novel allelic disorders), pedigree structure (e.g. imprinting disorders masquerading as autosomal recessive phenotypes), positional mapping (e.g. double recombination events abrogating candidate autozygous intervals), gene (e.g. novel gene-disease assertion) and variant (e.g. complex compound inheritance). Overall, we estimate a probability of 34.3% for encountering at least one of these challenges. Importantly, our data show that by only addressing non-sequencing-based challenges, around 71% increase in the diagnostic yield can be expected. Indeed, by applying these lessons to a cohort of 314 cases with negative clinical exome or genome reports, we could identify the likely causal variant in 54.5%. Our work highlights the need to have a thorough approach to undiagnosed diseases by considering a wide range of challenges rather than a narrow focus on sequencing technologies. It is hoped that by sharing this experience, the yield of undiagnosed disease programs globally can be improved.

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.

Family history is key to the interpretation of exome sequencing in the prenatal context: Unexpected diagnosis of Basal Cell Nevus Syndrome

OBJECTIVE

To reach a molecular diagnosis for a family with two consecutive fetuses presenting with multiple congenital anomalies.

METHOD

The two fetuses underwent prenatal ultrasound, autopsy, radiologic and genetic investigation. Genetic analysis included karyotype and array-CGH for both fetuses and trio-based whole exome sequencing (WES) only for the second fetus.

RESULTS

WES results, initially focusing on recessive or dominant de novo variants, were negative. However, as a result of new relevant information regarding family history, the variant c.648_651dup in the PTCH1 gene was identified as causative of the fetal phenotype.

CONCLUSION

This case further highlights how WES data analysis and interpretation strongly rely on family history and robust genotype-phenotype correlation. This is even more relevant in the prenatal setting, where access to fetal phenotype is limited and prenatal recognition of many morbid genes is not fully explored. We also provide a detailed description of the prenatal manifestations of Basal Cell Nevus Syndrome. This article is protected by copyright. All rights reserved.

ASHG 2020 Curt Stern Award introduction: Fowzan Sami Alkuraya

This article is based on the address given by the author at the 2020 virtual meeting of the American Society of Human Genetics (ASHG) on October 26, 2020. The video of the original address can be found at the ASHG website.

2020 Curt Stern Award address: a more perfect clinical genome-how consanguineous populations contribute to the medical annotation of the human genome

This article is based on the address given by the author at the 2020 virtual meeting of the American Society of Human Genetics (ASHG) on October 26, 2020. The video of the original address can be found at the ASHG website.

Common homozygosity for predicted loss-of-function variants reveals both redundant and advantageous effects of dispensable human genes

Humans homozygous or hemizygous for variants predicted to cause a loss of function (LoF) of the corresponding protein do not necessarily present with overt clinical phenotypes. We report here 190 autosomal genes with 207 predicted LoF variants, for which the frequency of homozygous individuals exceeds 1% in at least one human population from five major ancestry groups. No such genes were identified on the X and Y chromosomes. Manual curation revealed that 28 variants (15%) had been misannotated as LoF. Of the 179 remaining variants in 166 genes, only 11 alleles in 11 genes had previously been confirmed experimentally to be LoF. The set of 166 dispensable genes was enriched in olfactory receptor genes (41 genes). The 41 dispensable olfactory receptor genes displayed a relaxation of selective constraints similar to that observed for other olfactory receptor genes. The 125 dispensable nonolfactory receptor genes also displayed a relaxation of selective constraints consistent with greater redundancy. Sixty-two of these 125 genes were found to be dispensable in at least three human populations, suggesting possible evolution toward pseudogenes. Of the 179 LoF variants, 68 could be tested for two neutrality statistics, and 8 displayed robust signals of positive selection. These latter variants included a known FUT2 variant that confers resistance to intestinal viruses, and an APOL3 variant involved in resistance to parasitic infections. Overall, the identification of 166 genes for which a sizeable proportion of humans are homozygous for predicted LoF alleles reveals both redundancies and advantages of such deficiencies for human survival.

Lessons Learned from Large-Scale, First-Tier Clinical Exome Sequencing in a Highly Consanguineous Population

We report the results of clinical exome sequencing (CES) on >2,200 previously unpublished Saudi families as a first-tier test. The predominance of autosomal-recessive causes allowed us to make several key observations. We highlight 155 genes that we propose to be recessive, disease-related candidates. We report additional mutational events in 64 previously reported candidates (40 recessive), and these events support their candidacy. We report recessive forms of genes that were previously associated only with dominant disorders and that have phenotypes ranging from consistent with to conspicuously distinct from the known dominant phenotypes. We also report homozygous loss-of-function events that can inform the genetics of complex diseases. We were also able to deduce the likely causal variant in most couples who presented after the loss of one or more children, but we lack samples from those children. Although a similar pattern of mostly recessive causes was observed in the prenatal setting, the higher proportion of loss-of-function events in these cases was notable. The allelic series presented by the wealth of recessive variants greatly expanded the phenotypic expression of the respective genes. We also make important observations about dominant disorders; these observations include the pattern of de novo variants, the identification of 74 candidate dominant, disease-related genes, and the potential confirmation of 21 previously reported candidates. Finally, we describe the influence of a predominantly autosomal-recessive landscape on the clinical utility of rapid sequencing (Flash Exome). Our cohort's genotypic and phenotypic data represent a unique resource that can contribute to improved variant interpretation through data sharing.

Autozygome and high throughput confirmation of disease genes candidacy

Purpose

Establishing links between Mendelian phenotypes and genes enables the proper interpretation of variants therein. Autozygome, a rich source of homozygous variants, has been successfully utilized for the high throughput identification of novel autosomal recessive disease genes. Here, we highlight the utility of the autozygome for the high throughput confirmation of previously published tentative links to diseases.

Methods

Autozygome and exome analysis of patients with suspected Mendelian phenotypes. All variants were classified according to the American College of Medical Genetics and Genomics guidelines.

Results

We highlight 30 published candidate genes (ACTL6B, ADAM22, AGTPBP1, APC, C12orf4, C3orf17 (NEPRO), CENPF, CNPY3, COL27A1, DMBX1, FUT8, GOLGA2, KIAA0556, LENG8, MCIDAS, MTMR9, MYH11, QRSL1, RUBCN, SLC25A42, SLC9A1, TBXT, TFG, THUMPD1, TRAF3IP2, UFC1, UFM1, WDR81, XRCC2, ZAK) in which we identified homozygous likely deleterious variants in patients with compatible phenotypes. We also identified homozygous likely deleterious variants in 18 published candidate genes (ABCA2, ARL6IP1, ATP8A2, CDK9, CNKSR1, DGAT1, DMXL2, GEMIN4, HCN2, HCRT, MYO9A, PARS2, PLOD3, PREPL, SCLT1, STX3, TXNRD2, WIPI2) although the associated phenotypes are sufficiently different from the original reports that they represent phenotypic expansion or potentially distinct allelic disorders.

Conclusions

Our results should facilitate the timely relabeling of these candidate disease genes in relevant databases to improve the yield of clinical genomic sequencing.

Human genetics of infectious diseases: Unique insights into immunological redundancy

For almost any given human-tropic virus, bacterium, fungus, or parasite, the clinical outcome of primary infection is enormously variable, ranging from asymptomatic to lethal infection. This variability has long been thought to be largely determined by the germline genetics of the human host, and this is increasingly being demonstrated to be the case. The number and diversity of known inborn errors of immunity is continually increasing, and we focus here on autosomal and X-linked recessive traits underlying complete deficiencies of the encoded protein. Schematically, four types of infectious phenotype have been observed in individuals with such deficiencies, each providing information about the redundancy of the corresponding human gene, in terms of host defense in natural conditions. The lack of a protein can confer vulnerability to a broad range of microbes in most, if not all patients, through the disruption of a key immunological component. In such cases, the gene concerned is of low redundancy. However, the lack of a protein may also confer vulnerability to a narrow range of microbes, sometimes a single pathogen, and not necessarily in all patients. In such cases, the gene concerned is highly redundant. Conversely, the deficiency may be apparently neutral, conferring no detectable predisposition to infection in any individual. In such cases, the gene concerned is completely redundant. Finally, the lack of a protein may, paradoxically, be advantageous to the host, conferring resistance to one or more infections. In such cases, the gene is considered to display beneficial redundancy. These findings reflect the current state of evolution of humans and microbes, and should not be considered predictive of redundancy, or of a lack of redundancy, in the distant future. Nevertheless, these observations are of potential interest to present-day biologists testing immunological hypotheses experimentally and physicians managing patients with immunological or infectious conditions.

Using ALoFT to determine the impact of putative loss-of-function variants in protein-coding genes

Variants predicted to result in the loss of function of human genes have attracted interest because of their clinical impact and surprising prevalence in healthy individuals. Here, we present ALoFT (annotation of loss-of-function transcripts), a method to annotate and predict the disease-causing potential of loss-of-function variants. Using data from Mendelian disease-gene discovery projects, we show that ALoFT can distinguish between loss-of-function variants that are deleterious as heterozygotes and those causing disease only in the homozygous state. Investigation of variants discovered in healthy populations suggests that each individual carries at least two heterozygous premature stop alleles that could potentially lead to disease if present as homozygotes. When applied to de novo putative loss-of-function variants in autism-affected families, ALoFT distinguishes between deleterious variants in patients and benign variants in unaffected siblings. Finally, analysis of somatic variants in >6500 cancer exomes shows that putative loss-of-function variants predicted to be deleterious by ALoFT are enriched in known driver genes.Variants causing loss of function (LoF) of human genes have clinical implications. Here, the authors present a method to predict disease-causing potential of LoF variants, ALoFT (annotation of Loss-of-Function Transcripts) and show its application to interpreting LoF variants in different contexts.

Molecular autopsy in maternal-fetal medicine

PurposeThe application of genomic sequencing to investigate unexplained death during early human development, a form of lethality likely enriched for severe Mendelian disorders, has been limited.MethodsIn this study, we employed exome sequencing as a molecular autopsy tool in a cohort of 44 families with at least one death or lethal fetal malformation at any stage of in utero development. Where no DNA was available from the fetus, we performed molecular autopsy by proxy, i.e., through parental testing.ResultsPathogenic or likely pathogenic variants were identified in 22 families (50%), and variants of unknown significance were identified in further 15 families (34%). These variants were in genes known to cause embryonic or perinatal lethality (ALPL, GUSB, SLC17A5, MRPS16, THSD1, PIEZO1, and CTSA), genes known to cause Mendelian phenotypes that do not typically include embryonic lethality (INVS, FKTN, MYBPC3, COL11A2, KRIT1, ASCC1, NEB, LZTR1, TTC21B, AGT, KLHL41, GFPT1, and WDR81) and genes with no established links to human disease that we propose as novel candidates supported by embryonic lethality of their orthologs or other lines of evidence (MS4A7, SERPINA11, FCRL4, MYBPHL, PRPF19, VPS13D, KIAA1109, MOCS3, SVOPL, FEN1, HSPB11, KIF19, and EXOC3L2).ConclusionOur results suggest that molecular autopsy in pregnancy losses is a practical and high-yield alternative to traditional autopsy, and an opportunity for bringing precision medicine to the clinical practice of perinatology.

GWAS signals revisited using human knockouts

PurposeGenome-wide association studies (GWAS) have been instrumental to our understanding of the genetic risk determinants of complex traits. A common challenge in GWAS is the interpretation of signals, which are usually attributed to the genes closest to the polymorphic markers that display the strongest statistical association. Naturally occurring complete loss of function (knockout) of these genes in humans can inform GWAS interpretation by unmasking their deficiency state in a clinical context.MethodsWe exploited the unique population structure of Saudi Arabia to identify novel knockout events in genes previously highlighted in GWAS using combined autozygome/exome analysis.ResultsWe report five families with homozygous truncating mutations in genes that had only been linked to human disease through GWAS. The phenotypes observed in the natural knockouts for these genes (TRAF3IP2, FRMD3, RSRC1, BTBD9, and PXDNL) range from consistent with, to unrelated to, the previously reported GWAS phenotype.ConclusionWe expand the role of human knockouts in the medical annotation of the human genome, and show their potential value in informing the interpretation of GWAS of complex traits.

Autozygosity reveals recessive mutations and novel mechanisms in dominant genes: implications in variant interpretation

PURPOSE

The purpose of this study is to describe recessive alleles in strictly dominant genes. Identifying recessive mutations in genes for which only dominant disease or risk alleles have been reported can expand our understanding of the medical relevance of these genes both phenotypically and mechanistically. The Saudi population is enriched for autozygosity, which enhances the homozygous occurrence of alleles, including pathogenic alleles in genes that have been associated only with a dominant inheritance pattern.

METHODS

Exome sequencing of patients from consanguineous families with likely recessive phenotypes was performed. In one family, the genotype of the deceased children was inferred from their parents due to lack of available samples.

RESULTS

We describe the identification of 11 recessive variants (5 of which are reported here for the first time) in 11 genes for which only dominant disease or risk alleles have been reported. The observed phenotypes for these recessive variants were novel (e.g., FBN2-related myopathy and CSF1R-related brain malformation and osteopetrosis), typical (e.g., ACTG2-related visceral myopathy), or an apparently healthy state (e.g., PDE11A), consistent with the corresponding mouse knockout phenotypes.

CONCLUSION

Our results show that, in the era of genomic sequencing and "reverse phenotyping," recessive variants in dominant genes should not be dismissed based on perceived "incompatibility" with the patient's phenotype before careful consideration.Genet Med advance online publication 06 April 2017.

Clinical genomics expands the morbid genome of intellectual disability and offers a high diagnostic yield

Intellectual disability (ID) is a measurable phenotypic consequence of genetic and environmental factors. In this study, we prospectively assessed the diagnostic yield of genomic tools (molecular karyotyping, multi-gene panel and exome sequencing) in a cohort of 337 ID subjects as a first-tier test and compared it with a standard clinical evaluation performed in parallel. Standard clinical evaluation suggested a diagnosis in 16% of cases (54/337) but only 70% of these (38/54) were subsequently confirmed. On the other hand, the genomic approach revealed a likely diagnosis in 58% (n=196). These included copy number variants in 14% (n=54, 15% are novel), and point mutations revealed by multi-gene panel and exome sequencing in the remaining 43% (1% were found to have Fragile-X). The identified point mutations were mostly recessive (n=117, 81%), consistent with the high consanguinity of the study cohort, but also X-linked (n=8, 6%) and de novo dominant (n=19, 13%). When applied directly on all cases with negative molecular karyotyping, the diagnostic yield of exome sequencing was 60% (77/129). Exome sequencing also identified likely pathogenic variants in three novel candidate genes (DENND5A, NEMF and DNHD1) each of which harbored independent homozygous mutations in patients with overlapping phenotypes. In addition, exome sequencing revealed de novo and recessive variants in 32 genes (MAMDC2, TUBAL3, CPNE6, KLHL24, USP2, PIP5K1A, UBE4A, TP53TG5, ATOH1, C16ORF90, SLC39A14, TRERF1, RGL1, CDH11, SYDE2, HIRA, FEZF2, PROCA1, PIANP, PLK2, QRFPR, AP3B2, NUDT2, UFC1, BTN3A2, TADA1, ARFGEF3, FAM160B1, ZMYM5, SLC45A1, ARHGAP33 and CAPS2), which we highlight as potential candidates on the basis of several lines of evidence, and one of these genes (SLC39A14) was biallelically inactivated in a potentially treatable form of hypermanganesemia and neurodegeneration. Finally, likely causal variants in previously published candidate genes were identified (ASTN1, HELZ, THOC6, WDR45B, ADRA2B and CLIP1), thus supporting their involvement in ID pathogenesis. Our results expand the morbid genome of ID and support the adoption of genomics as a first-tier test for individuals with ID.

Opportunities and challenges of whole-genome and -exome sequencing

Recent advances in the development of sequencing technologies provide researchers with unprecedented possibilities for genetic analyses. In this review, we will discuss the history of genetic studies and the progress driven by next-generation sequencing (NGS), using complex inflammatory bowel diseases as an example. We focus on the opportunities, but also challenges that researchers are facing when working with NGS data to unravel the genetic causes underlying diseases.

Revisiting the morbid genome of Mendelian disorders

BACKGROUND

The pathogenicity of many Mendelian variants has been challenged by large-scale sequencing efforts. However, many rare and benign "disease mutations" are difficult to analyze due to their rarity. The Saudi Arabian variome is enriched for homozygosity due to inbreeding, a key advantage that can be exploited for the critical examination of previously published variants.

RESULTS

We collated all "disease-related mutations" listed in the Human Gene Mutation Database (HGMD) and ClinVar, including "variants of uncertain significance" (VOUS). We find that the use of public databases including 1000 Genomes, ExAC, and Kaviar can reclassify many of these variants as likely benign. Our Saudi Human Genome Program (SHGP) can reclassify many variants that are rare in public databases. Furthermore, SGPD allows us to observe many previously reported variants in the homozygous state and our extensive phenotyping of participants makes it possible to demonstrate the lack of phenotype for these variants, thus challenging their pathogenicity despite their rarity. We also find that 18 VOUS BRCA1 and BRCA2 variants that are listed in BRCA Exchange are present at least once in the homozygous state in patients who lack features of Fanconi anemia. Reassuringly, we could reciprocally demonstrate that none of those labeled as "pathogenic" were observed in the homozygous statue in individuals who lack Fanconi phenotype in our database.

CONCLUSION

Our study shows the importance of revisiting disease-related databases using public resources as well as of population-specific resources to improve the specificity of the morbid genome of Mendelian diseases in humans.

High-throughput discovery of novel developmental phenotypes

Approximately one-third of all mammalian genes are essential for life. Phenotypes resulting from knockouts of these genes in mice have provided tremendous insight into gene function and congenital disorders. As part of the International Mouse Phenotyping Consortium effort to generate and phenotypically characterize 5,000 knockout mouse lines, here we identify 410 lethal genes during the production of the first 1,751 unique gene knockouts. Using a standardized phenotyping platform that incorporates high-resolution 3D imaging, we identify phenotypes at multiple time points for previously uncharacterized genes and additional phenotypes for genes with previously reported mutant phenotypes. Unexpectedly, our analysis reveals that incomplete penetrance and variable expressivity are common even on a defined genetic background. In addition, we show that human disease genes are enriched for essential genes, thus providing a dataset that facilitates the prioritization and validation of mutations identified in clinical sequencing efforts.

NGS-based reverse genetic screen for common embryonic lethal mutations compromising fertility in livestock

We herein report the result of a large-scale, next generation sequencing (NGS)-based screen for embryonic lethal (EL) mutations in Belgian beef and New Zealand dairy cattle. We estimated by simulation that cattle might carry, on average, ∼0.5 recessive EL mutations. We mined exome sequence data from >600 animals, and identified 1377 stop-gain, 3139 frame-shift, 1341 splice-site, 22,939 disruptive missense, 62,399 benign missense, and 92,163 synonymous variants. We show that cattle have a comparable load of loss-of-function (LoF) variants (defined as stop-gain, frame-shift, or splice-site variants) as humans despite having a more variable exome. We genotyped >40,000 animals for up to 296 LoF and 3483 disruptive missense, breed-specific variants. We identified candidate EL mutations based on the observation of a significant depletion in homozygotes. We estimated the proportion of EL mutations at 15% of tested LoF and 6% of tested disruptive missense variants. We confirmed the EL nature of nine candidate variants by genotyping 200 carrier × carrier trios, and demonstrating the absence of homozygous offspring. The nine identified EL mutations segregate at frequencies ranging from 1.2% to 6.6% in the studied populations and collectively account for the mortality of ∼0.6% of conceptuses. We show that EL mutations preferentially affect gene products fulfilling basic cellular functions. The resulting information will be useful to avoid at-risk matings, thereby improving fertility.

A null mutation in TNIK defines a novel locus for intellectual disability

Intellectual disability (ID) is one of the most common disabilities and, although many genes have been implicated in its etiology, the genetic heterogeneity of ID continues to expand. The purpose of the study was to describe a novel autosomal recessive non-syndromic ID locus. Autozygome and linkage analysis, and exome sequencing followed by RNA and protein analysis of the candidate disease gene were performed. We describe two multiplex consanguineous families with non-syndromic ID phenotype, which maps to a critical linkage locus on 3q26. Exome sequencing of the index in each family revealed the same homozygous truncating mutation in TNIK that results in complete loss of the protein. TNIK is a kinase with a well-established role in dendrite development and synaptic transmission. The phenotype we observe in human patients who lack TNIK is consistent with the previously published Tnik (-/-) phenotype in the murine model. Our data strongly implicate TNIK deficiency in the causation of ID in humans.

Whole-exome sequencing in an isolated population from the Dalmatian island of Vis

We have whole-exome sequenced 176 individuals from the isolated population of the island of Vis in Croatia in order to describe exonic variation architecture. We found 290 577 single nucleotide variants (SNVs), 65% of which are singletons, low frequency or rare variants. A total of 25 430 (9%) SNVs are novel, previously not catalogued in NHLBI GO Exome Sequencing Project, UK10K-Generation Scotland, 1000Genomes Project, ExAC or NCBI Reference Assembly dbSNP. The majority of these variants (76%) are singletons. Comparable to data obtained from UK10K-Generation Scotland that were sequenced and analysed using the same protocols, we detected an enrichment of potentially damaging variants (non-synonymous and loss-of-function) in the low frequency and common variant categories. On average 115 (range 93–140) genotypes with loss-of-function variants, 23 (15–34) of which were homozygous, were identified per person. The landscape of loss-of-function variants across an exome revealed that variants mainly accumulated in genes on the xenobiotic-related pathways, of which majority coded for enzymes. The frequency of loss-of-function variants was additionally increased in Vis runs of homozygosity regions where variants mainly affected signalling pathways. This work confirms the isolate status of Vis population by means of whole-exome sequence and reveals the pattern of loss-of-function mutations, which resembles the trails of adaptive evolution that were found in other species. By cataloguing the exomic variants and describing the allelic structure of the Vis population, this study will serve as a valuable resource for future genetic studies of human diseases, population genetics and evolution in this population.

Testing the role of predicted gene knockouts in human anthropometric trait variation

Although the role of complete gene inactivation by two loss-of-function mutations inherited in trans is well-established in recessive Mendelian diseases, we have not yet explored how such gene knockouts (KOs) could influence complex human phenotypes. Here, we developed a statistical framework to test the association between gene KOs and quantitative human traits. Our method is flexible, publicly available, and compatible with common genotype format files (e.g. PLINK and vcf). We characterized gene KOs in 4498 participants from the NHLBI Exome Sequence Project (ESP) sequenced at high coverage (>100×), 1976 French Canadians from the Montreal Heart Institute Biobank sequenced at low coverage (5.7×), and >100 000 participants from the Genetic Investigation of ANthropometric Traits (GIANT) Consortium genotyped on an exome array. We tested associations between gene KOs and three anthropometric traits: body mass index (BMI), height and BMI-adjusted waist-to-hip ratio (WHR). Despite our large sample size and multiple datasets available, we could not detect robust associations between specific gene KOs and quantitative anthropometric traits. Our results highlight several limitations and challenges for future gene KO studies in humans, in particular when there is no prior knowledge on the phenotypes that might be affected by the tested gene KOs. They also suggest that gene KOs identified with current DNA sequencing methodologies probably do not strongly influence normal variation in BMI, height, and WHR in the general human population.

GOLGA2, encoding a master regulator of golgi apparatus, is mutated in a patient with a neuromuscular disorder

Golgi apparatus (GA) is a membrane-bound organelle that serves a multitude of critical cellular functions including protein secretion and sorting, and cellular polarity. Many Mendelian diseases are caused by mutations in genes encoding various components of GA. GOLGA2 encodes GM130, a necessary component for the assembly of GA as a single complex, and its deficiency has been found to result in severe cellular phenotypes. We describe the first human patient with a homozygous apparently loss of function mutation in GOLGA2. The phenotype is a neuromuscular disorder characterized by developmental delay, seizures, progressive microcephaly, and muscular dystrophy. Knockdown of golga2 in zebrafish resulted in severe skeletal muscle disorganization and microcephaly recapitulating loss of function human phenotype. Our data suggest an important developmental role of GM130 in humans and zebrafish.

TLE6 mutation causes the earliest known human embryonic lethality

BACKGROUND

Embryonic lethality is a recognized phenotypic expression of individual gene mutations in model organisms. However, identifying embryonic lethal genes in humans is challenging, especially when the phenotype is manifested at the preimplantation stage.

RESULTS

In an ongoing effort to exploit the highly consanguineous nature of the Saudi population to catalog recessively acting embryonic lethal genes in humans, we have identified two families with a female-limited infertility phenotype. Using autozygosity mapping and whole exome sequencing, we map this phenotype to a single mutation in TLE6, a maternal effect gene that encodes a member of the subcortical maternal complex in mammalian oocytes. Consistent with the published phenotype of mouse Tle6 mutants, embryos from female patients who are homozygous for the TLE6 mutation fail to undergo early cleavage, with resulting sterility. The human mutation abrogates TLE6 phosphorylation, a step that is reported to be critical for the PKA-mediated progression of oocyte meiosis II. Furthermore, the TLE6 mutation impairs its binding to components of the subcortical maternal complex.

CONCLUSION

In this first report of a human defect in a member of the subcortical maternal subcritical maternal complex, we show that the TLE6 mutation is gender-specific and leads to the earliest known human embryonic lethality phenotype.

Revisiting disease genes based on whole-exome sequencing in consanguineous populations

Assigning a causal role for genes in disease states is one of the most significant medical applications of human genetics research. The requirement for at least two different pathogenic alleles in the same gene in individuals with a similar phenotype to assign a causal link has not always been fully adhered to, and we now know that even two alleles may not necessarily constitute sufficient evidence. Autozygosity is a rich source of natural "knockout" events by virtue of rendering ancestral loss-of-function (LOF) variants homozygous. In this study, we exploit this phenomenon by examining 523 exomes enriched for autozygosity to call into question previously published disease links for several genes based on the identification of confirmed homozygous LOF variants in the absence of the purported diseases. This study highlights an additional advantage of consanguineous populations in the quest to improve the medical annotation of the human genome.

Loss of gene function and evolution of human phenotypes

Humans have acquired many distinct evolutionary traits after the human-chimpanzee divergence. These phenotypes have resulted from genetic changes that occurred in the human genome and were retained by natural selection. Comparative primate genome analyses reveal that loss-of-function mutations are common in the human genome. Some of these gene inactivation events were revealed to be associated with the emergence of advantageous phenotypes and were therefore positively selected and fixed in modern humans (the "less-ismore" hypothesis). Representative cases of human gene inactivation and their functional implications are presented in this review. Functional studies of additional inactive genes will provide insight into the molecular mechanisms underlying acquisition of various human-specific traits.

Identification of embryonic lethal genes in humans by autozygosity mapping and exome sequencing in consanguineous families

Background

Identifying genetic variants that lead to discernible phenotypes is the core of Mendelian genetics. An approach that considers embryonic lethality as a bona fide Mendelian phenotype has the potential to reveal novel genetic causes, which will further our understanding of early human development at a molecular level. Consanguineous families in which embryonic lethality segregates as a recessive Mendelian phenotype offer a unique opportunity for high throughput novel gene discovery as has been established for other recessive postnatal phenotypes.

Results

We have studied 24 eligible families using autozygosity mapping and whole-exome sequencing. In addition to revealing mutations in genes previously linked to embryonic lethality in severe cases, our approach revealed seven novel candidate genes (THSD1, PIGC, UBN1, MYOM1, DNAH14, GALNT14, and FZD6). A founder mutation in one of these genes, THSD1, which has been linked to vascular permeability, accounted for embryonic lethality in three of the study families. Unlike the other six candidate genes, we were able to identify a second mutation in THSD1 in a family with a less severe phenotype consisting of hydrops fetalis and persistent postnatal edema, which provides further support for the proposed link between this gene and embryonic lethality.

Conclusions

Our study represents an important step towards the systematic analysis of “embryonic lethal genes” in humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0681-6) contains supplementary material, which is available to authorized users.

Natural human knockouts and the era of genotype to phenotype

Complete loss of gene function in humans by naturally occurring biallelic loss-of-function mutations (human knockout) is not a new concept. However, the recent identification of human knockouts along the entire spectrum of health and disease by next-generation sequencing promises to unlock their full potential to accelerate the medical and functional annotation of the human genome.