Meta-analysis of 46,000 germline de novo mutations linked to human inherited disease.

BACKGROUND: De novo mutations (DNMs) are variants that occur anew in the offspring of noncarrier parents. They are not inherited from either parent but rather result from endogenous mutational processes involving errors of DNA repair/replication. These spontaneous errors play a significant role in the causation of genetic disorders, and their importance in the context of molecular diagnostic medicine has become steadily more apparent as more DNMs have been reported in the literature. In this study, we examined 46,489 disease-associated DNMs annotated by the Human Gene Mutation Database (HGMD) to ascertain their distribution across gene and disease categories. RESULTS: Most disease-associated DNMs reported to date are found to be associated with developmental and psychiatric disorders, a reflection of the focus of sequencing efforts over the last decade. Of the 13,277 human genes in which DNMs have so far been found, the top-10 genes with the highest proportions of DNM relative to gene size were H3-3 A, DDX3X, CSNK2B, PURA, ZC4H2, STXBP1, SCN1A, SATB2, H3-3B and TUBA1A. The distribution of CADD and REVEL scores for both disease-associated DNMs and those mutations not reported to be de novo revealed a trend towards higher deleteriousness for DNMs, consistent with the likely lower selection pressure impacting them. This contrasts with the non-DNMs, which are presumed to have been subject to continuous negative selection over multiple generations. CONCLUSION: This meta-analysis provides important information on the occurrence and distribution of disease-associated DNMs in association with heritable disease and should make a significant contribution to our understanding of this major type of mutation.

Naturally occurring genetic diseases caused by de novo variants in domestic animals.

With the advent of next-generation sequencing, an increasing number of cases of de novo variants in domestic animals have been reported in scientific literature primarily associated with clinically severe phenotypes. The emergence of new variants at each generation is a crucial aspect in understanding the pathology of early-onset diseases in animals and can provide valuable insights into similar diseases in humans. With the aim of collecting deleterious de novo variants in domestic animals, we searched the scientific literature and compiled reports on 42 de novo variants in 31 genes in domestic animals. No clear disease-associated phenotype has been established in humans for three of these genes (NUMB, ANKRD28 and KCNG1). For the remaining 28 genes, a strong similarity between animal and human phenotypes was recognized from available information in OMIM and OMIA, revealing the importance of comparative studies and supporting the use of domestic animals as natural models for human diseases, in line with the One Health approach.

A partial duplication of an X-linked gene exclusive of a primate lineage (Macaca).

Gene duplication plays a significant role in evolution. Paralogous gene copies may be lost due to the successive accumulation of deleterious mutations or remain active in the genome. In this work, a partial duplication of an X-linked region in the Macaca genus is identified and explored. Genomic comparisons reveal that the duplication encompasses the genes encoding ornithine transcarbamylase (OTC) and retinitis pigmentosa GTPase regulator (RPGR), spanning over 0.1 Mb on the chromosome 9 of Macaca. According to our analyses, the duplicated region of chromosome 9 involves partial coding sequences of both OTC and RPGR genes. Analyses of the selective pressures did not reveal significant differences in the ratio between nonsynonymous and synonymous mutations (w<1), suggesting that no selective pressures were acting in the evolutionary process. Reports for a biological role regarding some partial duplications exist in the literature, therefore, although being rare events, partial duplications of functionally important genes are worthy of study so that their impact can be explored.

Complex interactions between p.His558Arg and linked variants in the sodium voltage-gated channel alpha subunit 5 (Na V 1.5).

Common genetic polymorphisms may modify the phenotypic outcome when co-occurring with a disease-causing variant, and therefore understanding their modulating role in health and disease is of great importance. The polymorphic p.His558Arg variant of the sodium voltage-gated channel alpha subunit 5 (Na V 1.5) encoded by the SCN5A gene is a case in point, as several studies have shown it can modify the clinical phenotype in a number of cardiac diseases. To evaluate the genetic backgrounds associated with this modulating effect, we reanalysed previous electrophysiological findings regarding the p.His558Arg variant and further assessed its patterns of genetic diversity in human populations. The Na V 1.5 p.His558Arg variant was found to be in linkage disequilibrium with six other polymorphic variants that previously were also associated with cardiac traits in GWAS analyses. On account of this, incongruent reports that Arg558 allele can compensate, aggravate or have no effect on Na V 1.5, likely might have arose due to a role of p.His558Arg depending on the additional linked variants. Altogether, these results indicate a major influence of the epistatic interactions between SCN5A variants, revealing also that phenotypic severity may depend on the polymorphic background associated to each individual genome.

A Robust Assay to Monitor Ataxin-3 Amyloid Fibril Assembly.

Spinocerebellar ataxia type 3 (SCA3) is caused by the expansion of a glutamine repeat in the protein ataxin-3, which is deposited as intracellular aggregates in affected brain regions. Despite the controversial role of ataxin-3 amyloid structures in SCA3 pathology, the identification of molecules with the capacity to prevent aberrant self-assembly and stabilize functional conformation(s) of ataxin-3 is a key to the development of therapeutic solutions. Amyloid-specific kinetic assays are routinely used to measure rates of protein self-assembly in vitro and are employed during screening for fibrillation inhibitors. The high tendency of ataxin-3 to assemble into oligomeric structures implies that minor changes in experimental conditions can modify ataxin-3 amyloid assembly kinetics. Here, we determine the self-association rates of ataxin-3 and present a detailed study of the aggregation of normal and pathogenic ataxin-3, highlighting the experimental conditions that should be considered when implementing and validating ataxin-3 amyloid progress curves in different settings and in the presence of ataxin-3 interactors. This assay provides a unique and robust platform to screen for modulators of the first steps of ataxin-3 aggregation-a starting point for further studies with cell and animal models of SCA3.

The repertoire of the elongation of very long-chain fatty acids (Elovl) protein family is conserved in tambaqui (Colossoma macropomum): Gene expression profiles offer insights into the sexual differentiation process.

Elongation of very long-chain fatty acids (Elovl) proteins are critical players in the regulation of the length of a fatty acid. At present, eight members of the Elovl family (Elovl1-8), displaying a characteristic fatty acid substrate specificity, have been identified in vertebrates, including teleost fish. In general, Elovl1, Elovl3, Elovl6 and Elovl7 exhibit a substrate preference for saturated and monounsaturated fatty acids, while Elovl2, Elovl4, Elovl5 and Elovl8 use polyunsaturated fatty acids (PUFA) as substrates. PUFA elongases have received considerable attention in aquatic animals due to their involvement in the conversion of C18 PUFAs to long-chain polyunsaturated fatty acids (LC-PUFA). Here, we identified the full repertoire of elovl genes in the tambaqui Colossoma macropomum genome. A detailed phylogenetic and synteny analysis suggests a conservation of these genes among teleosts. Furthermore, based on RNAseq gene expression data, we discovered a gender bias expression of elovl genes during sex differentiation of tambaqui, toward future males. Our findings suggest a role of Elovl enzymes and fatty acid metabolism in tambaqui sexual differentiation.

Genetic Variability of the Functional Domains of Chromodomains Helicase DNA-Binding (CHD) Proteins.

In the past few years, there has been an increasing neuroscientific interest in understanding the function of mammalian chromodomains helicase DNA-binding (CHD) proteins due to their association with severe developmental syndromes. Mammalian CHDs include nine members (CHD1 to CHD9), grouped into subfamilies according to the presence of specific functional domains, generally highly conserved in evolutionary terms. Mutations affecting these domains hold great potential to disrupt protein function, leading to meaningful pathogenic scenarios, such as embryonic defects incompatible with life. Here, we analysed the evolution of CHD proteins by performing a comparative study of the functional domains of CHD proteins between orthologous and paralogous protein sequences. Our findings show that the highest degree of inter-species conservation was observed at Group II (CHD3, CHD4, and CHD5) and that most of the pathological variations documented in humans involve amino acid residues that are conserved not only between species but also between paralogs. The parallel analysis of both orthologous and paralogous proteins, in cases where gene duplications have occurred, provided extra information showing patterns of flexibility as well as interchangeability between amino acid positions. This added complexity needs to be considered when the impact of novel mutations is assessed in terms of evolutionary conservation.

Evolution and Functional Characteristics of the Novel elovl8 That Play Pivotal Roles in Fatty Acid Biosynthesis.

Elongation of very long-chain fatty acid (Elovl) proteins are key enzymes that catalyze the rate-limiting step in the fatty acid elongation pathway. The most recently discovered member of the Elovl family, Elovl8, has been proposed to be a fish-specific elongase with two gene paralogs described in teleosts. However, the biological functions of Elovl8 are still to be elucidated. In this study, we showed that in contrast to previous findings, elovl8 is not unique to teleosts, but displays a rather unique and ample phylogenetic distribution. For functional determination, we generated elovl8a (elovl8a-/-) and elovl8b (elovl8b-/-) zebrafish using CRISPR/Cas9 technology. Fatty acid composition in vivo and zebrafish liver cell experiments suggest that the substrate preference of Elovl8 overlapped with other existing Elovl enzymes. Zebrafish Elovl8a could elongate the polyunsaturated fatty acids (PUFAs) C18:2n-6 and C18:3n-3 to C20:2n-6 and C20:3n-3, respectively. Along with PUFA, zebrafish Elovl8b also showed the capacity to elongate C18:0 and C20:1. Gene expression quantification suggests that Elovl8a and Elovl8b may play a potentially important role in fatty acid biosynthesis. Overall, our results provide novel insights into the function of Elovl8a and Elovl8b, representing additional fatty acid elongases not previously described in chordates.

An ancestral nuclear receptor couple, PPAR-RXR, is exploited by organotins.

Environmental chemicals have been reported to greatly disturb the endocrine and metabolic systems of multiple animal species. A recent example involves the exploitation of the nuclear receptor (NR) heterodimeric pair composed by PPAR/RXR (peroxisome proliferator-activated receptor/retinoid X receptor), which shows lipid perturbation in mammalian species. While gene orthologues of both of these receptors have been described outside vertebrates, no functional characterization of PPAR has been carried in protostome lineages. We provide the first functional analysis of PPAR in Patella sp. (Mollusca), using model obesogens such as tributyltin (TBT), triphenyltin (TPT), and proposed natural ligands (fatty acid molecules). To gain further insights, we used site-directed mutagenesis to PPAR and replaced the tyrosine 277 by a cysteine (the human homologous amino acid and TBT anchor residue) and an alanine. Additionally, we explored the alterations in the fatty acid profiles after an exposure to the model obesogen TBT, in vivo. Our results show that TBT and TPT behave as an antagonist of Patella sp. PPAR/RXR and that the tyrosine 277 is important, but not essential in the response to TBT. Overall, these results suggest a relation between the response of the mollusc PPAR-RXR to TBT and the lipid profile alterations reported at environmentally relevant concentrations. Our findings highlight the importance of comparative analysis between protostome and deuterostome lineages to decipher the differential impact of environmental chemicals.

Compensatory epistasis explored by molecular dynamics simulations.

A non-negligible proportion of human pathogenic variants are known to be present as wild type in at least some non-human mammalian species. The standard explanation for this finding is that molecular mechanisms of compensatory epistasis can alleviate the mutations' otherwise pathogenic effects. Examples of compensated variants have been described in the literature but the interacting residue(s) postulated to play a compensatory role have rarely been ascertained. In this study, the examination of five human X-chromosomally encoded proteins (FIX, GLA, HPRT1, NDP and OTC) allowed us to identify several candidate compensated variants. Strong evidence for a compensated/compensatory pair of amino acids in the coagulation FIXa protein (involving residues 270 and 271) was found in a variety of mammalian species. Both amino acid residues are located within the 60-loop, spatially close to the 39-loop that performs a key role in coagulation serine proteases. To understand the nature of the underlying interactions, molecular dynamics simulations were performed. The predicted conformational change in the 39-loop consequent to the Glu270Lys substitution (associated with hemophilia B) appears to impair the protein's interaction with its substrate but, importantly, such steric hindrance is largely mitigated in those proteins that carry the compensatory residue (Pro271) at the neighboring amino acid position.

Common polymorphic OTC variants can act as genetic modifiers of enzymatic activity.

Understanding the role of common polymorphisms in modulating the clinical phenotype when they co-occur with a disease-causing lesion is of critical importance in medical genetics. We explored the impact of apparently neutral common polymorphisms, using the gene encoding the urea cycle enzyme, ornithine transcarbamylase (OTC), as a model system. Distinct combinations of genetic backgrounds embracing two missense polymorphisms were created in cis with the pathogenic p.Arg40His replacement. In vitro enzymatic assays revealed that the polymorphic variants were able to modulate OTC activity both in the presence or absence of the pathogenic lesion. First, we found that the combination of the minor alleles of polymorphisms p.Lys46Arg and p.Gln270Arg significantly enhanced enzymatic activity in the wild-type protein. Second, enzymatic assays revealed that the minor allele of the p.Gln270Arg polymorphism was capable of ameliorating OTC activity when combined in cis with the pathogenic p.Arg40His replacement. Structural analysis predicted that the minor allele of the p.Gln270Arg polymorphism would serve to stabilize the OTC wild-type protein, thereby corroborating the results of the experimental assays. Our findings demonstrate the potential importance of cis-interactions between common polymorphic variants and pathogenic missense mutations and illustrate how standing genetic variation can modulate protein function.

SLC35A2-CDG: Novel variant and review.

SLC35A2 encodes the X-linked transporter that carries uridine diphosphate (UDP)-galactose from the cytosol to the lumen of the Golgi apparatus and the endoplasmic reticulum. Pathogenic variants have been associated to a congenital disorder of glycosylation (CDG) with epileptic encephalopathy as a predominant feature. Among the sixty five patients described so far, a strong gender bias is observed as only seven patients are males. This work is a review and reports a SLC35A2-CDG in a male without epilepsy and with growth deficiency associated with decreased serum IGF1, minor neurological involvement, minor facial dysmorphism, and camptodactyly of fingers and toes. Sequence analysis revealed a hemizygosity for a novel de novo variant: c.233A > G (p.Lys78Arg) in SLC35A2. Further analysis of SLC35A2 sequence by comparing both orthologous and paralogous positions, revealed that not only the variant found in this study, but also most of the reported mutated positions are conserved in SLC35A2 orthologous, and many even in the paralogous SLC35A1 and SLC35A3. This is strong evidence that replacements at these positions will have a critical pathological effect and may also explain the gender bias observed among SLC35A2-CDG patients.

GBA3: a polymorphic pseudogene in humans that experienced repeated gene loss during mammalian evolution.

The gene encoding the cytosolic ß-glucosidase GBA3 shows pseudogenization due to a truncated allele (rs358231) that is polymorphic in humans. Since this enzyme is involved in the transformation of many plant ß-glycosides, this particular case of gene loss may have been influenced by dietary adaptations during evolution. In humans, apart from the inactivating allele, we found that GBA3 accumulated additional damaging mutations, implying an extensive GBA3 loss. The allelic distribution of loss-of-function alleles revealed significant differences between human populations which can be partially related with their staple diet. The analysis of mammalian orthologs disclosed that GBA3 underwent at least nine pseudogenization events. Most events of pseudogenization occurred in carnivorous lineages, suggesting a possible link to a ß-glycoside poor diet. However, GBA3 was also lost in omnivorous and herbivorous species, hinting that the physiological role of GBA3 is not fully understood and other unknown causes may underlie GBA3 pseudogenization. Such possibility relies upon a putative role in sialic acid biology, where GBA3 participates in a cellular network involving NEU2 and CMAH. Overall, our data shows that the recurrent loss of GBA3 in mammals is likely to represent an evolutionary endpoint of the relaxation of selective constraints triggered by diet-related factors.

PseudoChecker: an integrated online platform for gene inactivation inference.

The rapid expansion of high-quality genome assemblies, exemplified by ongoing initiatives such as the Genome-10K and i5k, demands novel automated methods to approach comparative genomics. Of these, the study of inactivating mutations in the coding region of genes, or pseudogenization, as a source of evolutionary novelty is mostly overlooked. Thus, to address such evolutionary/genomic events, a systematic, accurate and computationally automated approach is required. Here, we present PseudoChecker, the first integrated online platform for gene inactivation inference. Unlike the few existing methods, our comparative genomics-based approach displays full automation, a built-in graphical user interface and a novel index, PseudoIndex, for an empirical evaluation of the gene coding status. As a multi-platform online service, PseudoChecker simplifies access and usability, allowing a fast identification of disruptive mutations. An analysis of 30 genes previously reported to be eroded in mammals, and 30 viable genes from the same lineages, demonstrated that PseudoChecker was able to correctly infer 97% of loss events and 95% of functional genes, confirming its reliability. PseudoChecker is freely available, without login required, at http://pseudochecker.ciimar.up.pt.

The Echinodermata PPAR: Functional characterization and exploitation by the model lipid homeostasis regulator tributyltin.

The wide ecological relevance of lipid homeostasis modulators in the environment has been increasingly acknowledged. Tributyltin (TBT), for instance, was shown to cause lipid modulation, not only in mammals, but also in fish, molluscs, arthropods and rotifers. In vertebrates, TBT is known to interact with a nuclear receptor heterodimer module, formed by the retinoid X receptor (RXR) and the peroxisome proliferator-activated receptor (PPAR). These modulate the expression of genes involved in lipid homeostasis. In the present work, we isolated for the first time the complete coding region of the Echinodermata (Paracentrotus lividus) gene orthologues of PPAR and RXR and evaluated the ability of a model lipid homeostasis modulator, TBT, to interfere with the lipid metabolism in this species. Our results demonstrate that TBT alters the gonadal fatty acid composition and gene expression patterns: yielding sex-specific responses in fatty acid levels, including the decrease of eicosapentaenoic acid (C20:5 n-3, EPA) in males, and increase of arachidonic acid (20:4n-6, ARA) in females, and upregulation of long-chain acyl-CoA synthetase (acsl), ppar and rxr. Furthermore, an in vitro test using COS-1 cells as host and chimeric receptors with the ligand binding domain (LBD) of P. lividus PPAR and RXR shows that organotins (TBT and TPT (Triphenyltin)) suppressed activity of the heterodimer PPAR/RXR in a concentration-dependent manner. Together, these results suggest that TBT acts as a lipid homeostasis modulator at environmentally relevant concentrations in Echinodermata and highlight a possible conserved mode of action via the PPAR/RXR heterodimer.

Essential genetic findings in neurodevelopmental disorders.

Neurodevelopmental disorders (NDDs) represent a growing medical challenge in modern societies. Ever-increasing sophisticated diagnostic tools have been continuously revealing a remarkably complex architecture that embraces genetic mutations of distinct types (chromosomal rearrangements, copy number variants, small indels, and nucleotide substitutions) with distinct frequencies in the population (common, rare, de novo). Such a network of interacting players creates difficulties in establishing rigorous genotype-phenotype correlations. Furthermore, individual lifestyles may also contribute to the severity of the symptoms fueling a large spectrum of gene-environment interactions that have a key role on the relationships between genotypes and phenotypes.Herein, a review of the genetic discoveries related to NDDs is presented with the aim to provide useful general information for the medical community.

Convergent inactivation of the skin-specific C-C motif chemokine ligand 27 in mammalian evolution.

The appearance of mammalian-specific skin features was a key evolutionary event contributing for the elaboration of physiological processes such as thermoregulation, adequate hydration, locomotion, and inflammation. Skin inflammatory and autoimmune processes engage a population of skin-infiltrating T cells expressing a specific C-C chemokine receptor (CCR10) which interacts with an epidermal CC chemokine, the skin-specific C-C motif chemokine ligand 27 (CCL27). CCL27 is selectively produced in the skin by keratinocytes, particularly upon inflammation, mediating the adhesion and homing of skin-infiltrating T cells. Here, we examined the evolution and coding condition of Ccl27 in 112 placental mammalian species. Our findings reveal that a number of open reading frame inactivation events such as insertions, deletions, and start and stop codon mutations independently occurred in Cetacea, Pholidota, Sirenia, Chiroptera, and Rodentia, totalizing 18 species. The diverse habitat settings and lifestyles of Ccl27-eroded lineages probably implied distinct evolutionary triggers rendering this gene unessential. For example, in Cetacea, the rapid renewal of skin layers minimizes the need for an elaborate inflammatory mechanism, mirrored by the absence of epidermal scabs. Our findings suggest that the convergent and independent loss of Ccl27 in mammalian evolution concurred with unique adaptive roads for skin physiology.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea.

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of "gene vestiges" in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

The Singularity of Cetacea Behavior Parallels the Complete Inactivation of Melatonin Gene Modules.

: Melatonin, the hormone of darkness, is a peculiar molecule found in most living organisms. Emerging as a potent broad-spectrum antioxidant, melatonin was repurposed into extra roles such as the modulation of circadian and seasonal rhythmicity, affecting numerous aspects of physiology and behaviour, including sleep entrainment and locomotor activity. Interestingly, the pineal gland-the melatonin synthesising organ in vertebrates-was suggested to be absent or rudimentary in some mammalian lineages, including Cetacea. In Cetacea, pineal regression is paralleled by their unique bio-rhythmicity, as illustrated by the unihemispheric sleeping behaviour and long-term vigilance. Here, we examined the genes responsible for melatonin synthesis (Aanat and Asmt) and signalling (Mtnr1a and Mtnr1b) in 12 toothed and baleen whale genomes. Based on an ample genomic comparison, we deduce that melatonin-related gene modules are eroded in Cetacea.

A complete enzymatic capacity for long-chain polyunsaturated fatty acid biosynthesis is present in the Amazonian teleost tambaqui, Colossoma macropomum.

In vertebrates, the essential fatty acids (FA) that satisfy the dietary requirements for a given species depend upon its desaturation and elongation capabilities to convert the C18 polyunsaturated fatty acids (PUFA), namely linoleic acid (LA, 18:2n-6) and α-linolenic acid (ALA, 18:3n-3), into the biologically active long-chain (C20-24) polyunsaturated fatty acids (LC-PUFA), including arachidonic acid (ARA, 20:4n-6), eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3). Recent studies have established that tambaqui (Colossoma macropomum), an important aquaculture-produced species in Brazil, is a herbivorous fish that can fulfil its essential FA requirements with dietary provision C18 PUFA LA and ALA, although the molecular mechanisms underpinning such ability remained unclear. The present study aimed at cloning and functionally characterizing genes encoding key desaturase and elongase enzymes, namely fads2, elovl5 and elovl2, involved in the LC-PUFA biosynthetic pathways in tambaqui. First, a fads2-like desaturase was isolated from tambaqui. When expressed in yeast, the tambaqui Fads2 showed Δ6, Δ5 and Δ8 desaturase capacities within the same enzyme, enabling all desaturation reactions required for ARA, EPA and DHA biosynthesis. Moreover, tambaqui possesses two elongases that are bona fide orthologs of elovl5 and elovl2. Their functional characterization confirmed that they can operate towards a variety of PUFA substrates with chain lengths ranging from 18 to 22 carbons. Overall our results provide compelling evidence that demonstrates that all the desaturase and elongase activities required to convert LA and ALA into ARA, EPA and DHA are present in tambaqui within the three genes studied herein, i.e. fads2, elovl5 and elovl2.