Genetics in Light of Transcriptional Adaptation

Loss of cped1 does not affect bone and lean mass in zebrafish

Human genetic studies have nominated Cadherin-like and PC-esterase Domain-containing 1 ( CPED1 ) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

Neurobeachin regulates hematopoietic progenitor differentiation and survival by modulating Notch activity

ABSTRACT

Hematopoietic stem cells (HSCs) can generate all blood cells. This ability is exploited in HSC transplantation (HSCT) to treat hematologic disease. A clear understanding of the molecular mechanisms that regulate HSCT is necessary to continue improving transplant protocols. We identified the Beige and Chediak-Higashi domain-containing protein (BDCP), Neurobeachin (NBEA), as a putative regulator of HSCT. Here, we demonstrated that NBEA and related BDCPs, including LPS Responsive Beige-Like Anchor Protein (LRBA), Neurobeachin Like 1 (NBEAL1) and Lysosomal Trafficking Regulator (LYST), are required during HSCT to efficiently reconstitute the hematopoietic system of lethally irradiated mice. Nbea knockdown in mouse HSCs induced apoptosis and a differentiation block after transplantation. Nbea deficiency in hematopoietic progenitor cells perturbed the expression of genes implicated in vesicle trafficking and led to changes in NOTCH receptor localization. This resulted in perturbation of the NOTCH transcriptional program, which is required for efficient HSC engraftment. In summary, our findings reveal a novel role for NBEA in the control of NOTCH receptor turnover in hematopoietic cells and supports a model in which BDCP-regulated vesicle trafficking is required for efficient HSCT.

Mechanisms governing vesicle traffic at the Golgi apparatus

Vesicle transport at the Golgi apparatus is a well-described process, and the major protein components involved have been identified. This includes the coat proteins that function in cargo sorting and vesicle formation, and the proteins that mediate the downstream events of vesicle tethering and membrane fusion. However, despite this knowledge, there remain significant gaps in our mechanistic understanding of these processes which includes how they are coordinated in space and time. In this review we discuss recent advances that have provided new insights into the mechanisms of Golgi trafficking, focussing on vesicle formation and cargo sorting, and vesicle tethering and fusion. These studies point to a high degree of spatial organisation of trafficking components at the Golgi and indicate an inherent plasticity of trafficking. Going forward, further advancements in technology and more sophisticated functional assays are expected to yield greater understanding of the mechanisms that govern Golgi trafficking events.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or nuclear factor κB (NF-κB) inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Evaluating the association of biallelic OGDHL variants with significant phenotypic heterogeneity

BACKGROUND

Biallelic variants in OGDHL, encoding part of the α-ketoglutarate dehydrogenase complex, have been associated with highly heterogeneous neurological and neurodevelopmental disorders. However, the validity of this association remains to be confirmed. A second OGDHL patient cohort was recruited to carefully assess the gene-disease relationship.

METHODS

Using an unbiased genotype-first approach, we screened large, multiethnic aggregated sequencing datasets worldwide for biallelic OGDHL variants. We used CRISPR/Cas9 to generate zebrafish knockouts of ogdhl, ogdh paralogs, and dhtkd1 to investigate functional relationships and impact during development. Functional complementation with patient variant transcripts was conducted to systematically assess protein functionality as a readout for pathogenicity.

RESULTS

A cohort of 14 individuals from 12 unrelated families exhibited highly variable clinical phenotypes, with the majority of them presenting at least one additional variant, potentially accounting for a blended phenotype and complicating phenotypic understanding. We also uncovered extreme clinical heterogeneity and high allele frequencies, occasionally incompatible with a fully penetrant recessive disorder. Human cDNA of previously described and new variants were tested in an ogdhl zebrafish knockout model, adding functional evidence for variant reclassification. We disclosed evidence of hypomorphic alleles as well as a loss-of-function variant without deleterious effects in zebrafish variant testing also showing discordant familial segregation, challenging the relationship of OGDHL as a conventional Mendelian gene. Going further, we uncovered evidence for a complex compensatory relationship among OGDH, OGDHL, and DHTKD1 isoenzymes that are associated with neurodevelopmental disorders and exhibit complex transcriptional compensation patterns with partial functional redundancy.

CONCLUSIONS

Based on the results of genetic, clinical, and functional studies, we formed three hypotheses in which to frame observations: biallelic OGDHL variants lead to a highly variable monogenic disorder, variants in OGDHL are following a complex pattern of inheritance, or they may not be causative at all. Our study further highlights the continuing challenges of assessing the validity of reported disease-gene associations and effects of variants identified in these genes. This is particularly more complicated in making genetic diagnoses based on identification of variants in genes presenting a highly heterogenous phenotype such as "OGDHL-related disorders".

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or NF-kB inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Partial sequence identity in a 25-nucleotide long element is sufficient for transcriptional adaptation in the Caenorhabditis elegans act-5/act-3 model

Genetic robustness can be achieved via several mechanisms including transcriptional adaptation (TA), a sequence similarity-driven process whereby mutant mRNA degradation products modulate, directly or indirectly, the expression of so-called adapting genes. To identify the sequences required for this process, we utilized a transgenic approach in Caenorhabditis elegans, combining an overexpression construct for a mutant gene (act-5) and a fluorescent reporter for the corresponding adapting gene (act-3). Analyzing a series of modifications for each construct, we identified, in the 5' regulatory region of the act-3 locus, a 25-base pair (bp) element which exhibits 60% identity with a sequence in the act-5 mRNA and which, in the context of a minimal promoter, is sufficient to induce ectopic expression of the fluorescent reporter. The 25 nucleotide (nt) element in the act-5 mRNA lies between the premature termination codon (PTC) and the next exon/exon junction, suggesting the importance of this region of the mutant mRNA for TA. Additionally, we found that single-stranded RNA injections of this 25 nt element from act-5 into the intestine of wild-type larvae led to higher levels of adapting gene (act-3) mRNA. Different models have been proposed to underlie the modulation of gene expression during TA including chromatin remodeling, the inhibition of antisense RNAs, the release of transcriptional pausing, and the suppression of premature transcription termination, and our data clearly show the importance of the regulatory region of the adapting gene in this particular act-5/act-3 TA model. Our findings also suggest that RNA fragments can modulate the expression of loci exhibiting limited sequence similarity, possibly a critical observation when designing RNA based therapies.

Polyphenol oxidases regulate pollen development through modulating flavonoids homeostasis in tobacco

Pollen development is critical in plant reproduction. Polyphenol oxidases (PPOs) genes encode defense-related enzymes, but the role of PPOs in pollen development remains largely unexplored. Here, we characterized NtPPO genes, and then investigated their function in pollen via creating NtPPO9/10 double knockout mutant (cas-1), overexpression 35S::NtPPO10 (cosp) line and RNAi lines against all NtPPOs in Nicotiana tabacum. NtPPOs were abundantly expressed in the anther and pollen (especially NtPPO9/10). The pollen germination, polarity ratio and fruit weights were significantly reduced in the NtPPO-RNAi and cosp lines, while they were normal in cas-1 likely due to compensation by other NtPPO isoforms. Comparisons of metabolites and transcripts between the pollen of WT and NtPPO-RNAi, or cosp showed that decreased enzymatic activity of NtPPOs led to hyper-accumulation of flavonoids. This accumulation might reduce the content of ROS. Ca²⁺ and actin levels also decreased in pollen of the transgenic lines.Thus, the NtPPOs regulate pollen germination through the flavonoid homeostasis and ROS signal pathway. This finding provides novel insights into the native physiological functions of PPOs in pollen during reproduction.

Dopey2 and Pcdh7 orchestrate the development of embryonic neural stem cells/ progenitors in zebrafish

DOPEY2 has been shown to be associated with Down syndrome and PCDH7 might be involved in Rett syndrome and MECP2 duplication syndrome. The mechanism how both proteins play roles in these syndromes are largely unknown. Here, we show that Dopey2 and Pcdh7 balance the proliferation and differentiation of neural stem cells and progenitors during embryonic neurogenesis to generate proper size and architecture of zebrafish brains. Dopey2 and Pcdh7 mutually restricted expression of each other in zebrafish embryos. Dopey2 was responsible for the proliferation of neural stem cells/progenitors, whereas Pcdh7 was responsible for the differentiation of neural stem cells/progenitors. Both proteins were shown to orchestrate the proper development and arrangement of neural cells in zebrafish embryonic brains. The results provide an insight into mechanisms to understand how the embryonic brain is constituted and how developmental defects occur in the brains of patients with Down syndrome, Rett syndrome, or MECP2 duplication syndrome.

Multiomics and optobiotechnological approaches for the development of microalgal strain for production of aviation biofuel and biorefinery

Demand and consumption of fossil fuels is increasing daily, and oil reserves are depleting. Technological developments are required towards developing sustainable renewable energy sources and microalgae are emerging as a potential candidate for various application-driven research. Molecular understanding attained through omics and system biology approach empowering researchers to modify various metabolic pathways of microalgal system for efficient extraction of biofuel and important biomolecules. This review furnish insight into different "advanced approaches" like optogenetics, systems biology and multi-omics for enhanced production of FAS (Fatty Acid Synthesis) and lipids in microalgae and their associated challenges. These new approaches would be helpful in the path of developing microalgae inspired technological platforms for optobiorefinery, which could be explored as source material to produce biofuels and other valuable bio-compounds on a large scale.

The p190 RhoGAPs, ARHGAP35, and ARHGAP5 are implicated in GnRH neuronal development: Evidence from patients with idiopathic hypogonadotropic hypogonadism, zebrafish, and in vitro GAP activity assay

PURPOSE

The study aimed to identify novel genes for idiopathic hypogonadotropic hypogonadism (IHH).

METHODS

A cohort of 1387 probands with IHH underwent exome sequencing and de novo, familial, and cohort-wide investigations. Functional studies were performed on 2 p190 Rho GTPase-activating proteins (p190 RhoGAP), ARHGAP35 and ARHGAP5, which involved in vivo modeling in larval zebrafish and an in vitro p190A-GAP activity assay.

RESULTS

Rare protein-truncating variants (PTVs; n = 5) and missense variants in the RhoGAP domain (n = 7) in ARHGAP35 were identified in IHH cases (rare variant enrichment: PTV [unadjusted P = 3.1E-06] and missense [adjusted P = 4.9E-03] vs controls). Zebrafish modeling using gnrh3:egfp phenotype assessment showed that mutant larvae with deficient arhgap35a, the predominant ARHGAP35 paralog in the zebrafish brain, display decreased GnRH3-GFP+ neuronal area, a readout for IHH. In vitro GAP activity studies showed that 1 rare missense variant [ARHGAP35 p.(Arg1284Trp)] had decreased GAP activity. Rare PTVs (n = 2) also were discovered in ARHGAP5, a paralog of ARHGAP35; however, arhgap5 zebrafish mutants did not display significant GnRH3-GFP+ abnormalities.

CONCLUSION

This study identified ARHGAP35 as a new autosomal dominant genetic driver for IHH and ARHGAP5 as a candidate gene for IHH. These observations suggest a novel role for the p190 RhoGAP proteins in GnRH neuronal development and integrity.

Exon-dependent transcriptional adaptation by exon-junction complex proteins Y14/RNP-4 and MAGOH/MAG-1 in Caenorhabditis elegans

Transcriptional adaptation is a powerful gene regulation mechanism that can increase genetic robustness. Transcriptional adaptation occurs when a gene is mutated and is mediated by the mutant RNA, rather than by protein feedback loops. We show here that transcriptional adaptation occurs in the C. elegans clh family of Cl- channels and that it requires exon-junction complex (EJC) proteins RNP-4, MAG-1, and eiF4AIII. Depending on which exons are deleted in distinct clh-1 alleles, different clh genes are regulated in an EJC-dependent manner. Our results support the idea that different transcriptional adaptation outcomes may be directed by the differential interaction of the EJC with its target mutant RNAs.

Genetic compensation of triacylglycerol biosynthesis in the green microalga Chlamydomonas reinhardtii

Genetic compensation has been proposed to explain phenotypic differences between gene knockouts and knockdowns in several metazoan and plant model systems. With the rapid development of reverse genetic tools such as CRISPR/Cas9 and RNAi in microalgae, it is increasingly important to assess whether genetic compensation affects the phenotype of engineered algal mutants. While exploring triacylglycerol (TAG) biosynthesis pathways in the model alga Chlamydomonas reinhardtii, it was discovered that knockout of certain genes catalyzing rate-limiting steps of TAG biosynthesis, type-2 diacylglycerol acyltransferase genes (DGTTs), triggered genetic compensation under abiotic stress conditions. Genetic compensation of a DGTT1 null mutation by a related PDAT gene was observed regardless of the strain background or mutagenesis approach, for example, CRISPR/Cas 9 or insertional mutagenesis. However, no compensation was found in the PDAT knockout mutant. The effect of PDAT knockout was evaluated in a Δvtc1 mutant, in which PDAT was upregulated under stress, resulting in a 90% increase in TAG content. Knockout of PDAT in the Δvtc1 background induced a 12.8-fold upregulation of DGTT1 and a 272.3% increase in TAG content in Δvtc1/pdat1 cells, while remaining viable. These data suggest that genetic compensation contributes to the genetic robustness of microalgal TAG biosynthetic pathways, maintaining lipid and redox homeostasis in the knockout mutants under abiotic stress. This work demonstrates examples of genetic compensation in microalgae, implies the physiological relevance of genetic compensation in TAG biosynthesis under stress, and provides guidance for future genetic engineering and mutant characterization efforts.

Zebrafish models of alx-linked frontonasal dysplasia reveal a role for Alx1 and Alx3 in the anterior segment and vasculature of the developing eye

The cellular and genetic mechanisms that coordinate formation of facial sensory structures with surrounding skeletal and soft tissue elements remain poorly understood. Alx1, a homeobox transcription factor, is a key regulator of midfacial morphogenesis. ALX1 mutations in humans are linked to severe congenital anomalies of the facial skeleton (frontonasal dysplasia, FND) with malformation or absence of eyes and orbital contents (micro- and anophthalmia). Zebrafish with loss-of-function alx1 mutations develop with craniofacial and ocular defects of variable penetrance, likely due to compensatory upregulation in expression of a paralogous gene, alx3. Here we show that zebrafish alx1;alx3 mutants develop with highly penetrant cranial and ocular defects that resemble human ALX1-linked FND. alx1 and alx3 are expressed in anterior cranial neural crest (aCNC), which gives rise to the anterior neurocranium (ANC), anterior segment structures of the eye and vascular pericytes. Consistent with a functional requirement for alx genes in aCNC, alx1; alx3 mutants develop with nearly absent ANC and grossly aberrant hyaloid vasculature and ocular anterior segment, but normal retina. In vivo lineage labeling identified a requirement for alx1 and alx3 during aCNC migration, and transcriptomic analysis suggested oxidative stress response as a key target mechanism of this function. Oxidative stress is a hallmark of fetal alcohol toxicity, and we found increased penetrance of facial and ocular malformations in alx1 mutants exposed to ethanol, consistent with a protective role for alx1 against ethanol toxicity. Collectively, these data demonstrate a conserved role for zebrafish alx genes in controlling ocular and facial development, and a novel role in protecting these key midfacial structures from ethanol toxicity during embryogenesis. These data also reveal novel roles for alx genes in ocular anterior segment formation and vascular development and suggest that retinal deficits in alx mutants may be secondary to aberrant ocular vascularization and anterior segment defects. This study establishes robust zebrafish models for interrogating conserved genetic mechanisms that coordinate facial and ocular development, and for exploring gene--environment interactions relevant to fetal alcohol syndrome.

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

Genotype-Phenotype Relationships in the Context of Transcriptional Adaptation and Genetic Robustness

Genetic manipulations with a robust and predictable outcome are critical to investigate gene function, as well as for therapeutic genome engineering. For many years, knockdown approaches and reagents including RNA interference and antisense oligonucleotides dominated functional studies; however, with the advent of precise genome editing technologies, CRISPR-based knockout systems have become the state-of-the-art tools for such studies. These technologies have helped decipher the role of thousands of genes in development and disease. Their use has also revealed how limited our understanding of genotype-phenotype relationships is. The recent discovery that certain mutations can trigger the transcriptional modulation of other genes, a phenomenon called transcriptional adaptation, has provided an additional explanation for the contradicting phenotypes observed in knockdown versus knockout models and increased awareness about the use of each of these approaches. In this review, we first cover the strengths and limitations of different gene perturbation strategies. Then we highlight the diverse ways in which the genotype-phenotype relationship can be discordant between these different strategies. Finally, we review the genetic robustness mechanisms that can lead to such discrepancies, paying special attention to the recently discovered phenomenon of transcriptional adaptation. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Screening and validation of genome-edited animals

The emergence of an array of genome-editing tools in recent years has facilitated the introduction of genetic modifications directly into the embryo, increasing the ease, efficiency and catalogue of alleles accessible to researchers across a range of species. Bypassing the requirement for a selection cassette and resulting in a broad range of outcomes besides the desired allele, genome editing has altered the allele validation process both temporally and technically. Whereas traditional gene targeting relies upon selection and allows allele validation at the embryonic stem cell modification stage, screening for the presence of the intended allele now occurs in the (frequently mosaic) founder animals. Final confirmation of the edited allele can only take place at the subsequent G1 generation and the validation strategy must differentiate the desired allele from a range of unintended outcomes. Here we present some of the challenges posed by gene editing, strategies for validation and considerations for animal colony management.

Zebrafish Chromosome 14 Gene Differential Expression in the fmr1 h u2787 Model of Fragile X Syndrome

Zebrafish represent a valuable model for investigating the molecular and cellular basis of Fragile X syndrome (FXS). Reduced expression of the zebrafish FMR1 orthologous gene, fmr1, causes developmental and behavioural phenotypes related to FXS. Zebrafish homozygous for the hu2787 non-sense mutation allele of fmr1 are widely used to model FXS, although FXS-relevant phenotypes seen from morpholino antisense oligonucleotide (morpholino) suppression of fmr1 transcript translation were not observed when hu2787 was first described. The subsequent discovery of transcriptional adaptation (a form of genetic compensation), whereby mutations causing non-sense-mediated decay of transcripts can drive compensatory upregulation of homologous transcripts independent of protein feedback loops, suggested an explanation for the differences reported. We examined the whole-embryo transcriptome effects of homozygosity for fmr1 h u2787 at 2 days post fertilisation. We observed statistically significant changes in expression of a number of gene transcripts, but none from genes showing sequence homology to fmr1. Enrichment testing of differentially expressed genes implied effects on lysosome function and glycosphingolipid biosynthesis. The majority of the differentially expressed genes are located, like fmr1, on Chromosome 14. Quantitative PCR tests did not support that this was artefactual due to changes in relative chromosome abundance. Enrichment testing of the "leading edge" differentially expressed genes from Chromosome 14 revealed that their co-location on this chromosome may be associated with roles in brain development and function. The differential expression of functionally related genes due to mutation of fmr1, and located on the same chromosome as fmr1, is consistent with R.A. Fisher's assertion that the selective advantage of co-segregation of particular combinations of alleles of genes will favour, during evolution, chromosomal rearrangements that place them in linkage disequilibrium on the same chromosome. However, we cannot exclude that the apparent differential expression of genes on Chromosome 14 genes was, (if only in part), caused by differences between the expression of alleles of genes unrelated to the effects of the fmr1 h u2787 mutation and made manifest due to the limited, but non-zero, allelic diversity between the genotypes compared.

Parallel genetics of regulatory sequences using scalable genome editing in vivo

How regulatory sequences control gene expression is fundamental for explaining phenotypes in health and disease. Regulatory elements must ultimately be understood within their genomic environment and development- or tissue-specific contexts. Because this is technically challenging, few regulatory elements have been characterized in vivo. Here, we use inducible Cas9 and multiplexed guide RNAs to create hundreds of mutations in enhancers/promoters and 3' UTRs of 16 genes in C. elegans. Our software crispr-DART analyzes indel mutations in targeted DNA sequencing. We quantify the impact of mutations on expression and fitness by targeted RNA sequencing and DNA sampling. When applying our approach to the lin-41 3' UTR, generating hundreds of mutants, we find that the two adjacent binding sites for the miRNA let-7 can regulate lin-41 expression independently of each other. Finally, we map regulatory genotypes to phenotypic traits for several genes. Our approach enables parallel analysis of regulatory sequences directly in animals.

Differential intolerance to loss of function and missense mutations in genes that encode human matricellular proteins

Targeted gene disruption in mice has provided valuable insights into the functions of matricellular proteins. Apart from missense and loss of function mutations that have been associated with inherited diseases, however, their functions in humans remain unclear. The availability of deep exome sequencing data from over 140,000 individuals in the Genome Aggregation Database provided an opportunity to examine intolerance to loss of function and missense mutations in human matricellular genes. The probability of loss-of-function intolerance (pLI) differed widely within members of the thrombospondin, CYR61/CTGF/NOV (CCN), tenascin, small integrin-binding ligand N-linked glycoproteins (SIBLING), and secreted protein, acidic and rich in cysteine (SPARC) gene families. Notably, pLI values in humans had limited correlation with viability of the corresponding homozygous null mice. Among the thrombospondins, only THBS1 was highly loss-intolerant (pLI = 1). In contrast, Thbs1 is not essential for viability in mice. Several known thrombospondin-1 receptors were similarly loss-intolerant, although thrombospondin-1 is not the exclusive ligand for some of these receptors. The frequencies of missense mutations in THBS1 and the gene encoding its signaling receptor CD47 indicated conservation of some residues implicated in specific receptor binding. Deficits in missense mutations were also observed for other thrombospondin genes and for SPARC, SPOCK1, SPOCK2, TNR, and DSPP. The intolerance of THBS1 to loss of function in humans and elevated pLI values for THBS2, SPARC, SPOCK1, TNR, and CCN1 support important functions for these matricellular protein genes in humans, some of which may relate to functions in reproduction or responding to environmental stresses.