Caseins are among cardinal proteins that evolved in the lineage leading to mammals. In milk, caseins and calcium phosphate (CaP) form a huge complex called casein micelle. By forming the micelle, milk maintains high CaP concentrations, which help altricial mammalian neonates to grow bone and teeth. Two types of caseins are known. Ca-sensitive caseins (α(s)- and β-caseins) bind Ca but precipitate at high Ca concentrations, whereas Ca-insensitive casein (κ-casein) does not usually interact with Ca but instead stabilizes the micelle. Thus, it is thought that these two types of caseins are both necessary for stable micelle formation. Both types of caseins show high substitution rates, which make it difficult to elucidate the evolution of caseins. Yet, recent studies have revealed that all casein genes belong to the secretory calcium-binding phosphoprotein (SCPP) gene family that arose by gene duplication. In the present study, we investigated exon-intron structures and phylogenetic distributions of casein and other SCPP genes, particularly the odontogenic ameloblast-associated (ODAM) gene, the SCPP-Pro-Gln-rich 1 (SCPPPQ1) gene, and the follicular dendritic cell secreted peptide (FDCSP) gene. The results suggest that contemporary Ca-sensitive casein genes arose from a putative common ancestor, which we refer to as CSN1/2. The six putative exons comprising CSN1/2 are all found in SCPPPQ1, although ODAM also shares four of these exons. By contrast, the five exons of the Ca-insensitive casein gene are all reminiscent of FDCSP. The phylogenetic distribution of these genes suggests that both SCPPPQ1 and FDCSP arose from ODAM. We thus argue that all casein genes evolved from ODAM via two different pathways; Ca-sensitive casein genes likely originated directly from SCPPPQ1, whereas the Ca-insensitive casein genes directly differentiated from FDCSP. Further, expression of ODAM, SCPPPQ1, and FDCSP was detected in dental tissues, supporting the idea that both types of caseins evolved as Ca-binding proteins. Based on these findings, we propose two alternative hypotheses for micelle formation in primitive milk. The conserved biochemical characteristics in caseins and their immediate ancestors also suggest that many slight genetic modifications have created modern caseins, proteins vital to the sustained success of mammals.

Enamelin (ENAM) has been shown to be a crucial protein for enamel formation and mineralization. Previous molecular analyses have indicated a probable origin early in vertebrate evolution, which is supported by the presence of enamel/enameloid tissues in early vertebrates. In contrast to these hypotheses, ENAM was only characterized in mammals. Our aims were to 1) look for ENAM in representatives of nonmammalian tetrapods, 2) search for a pseudogene in the chicken genome, and 3) see whether the new sequences could bring new information on ENAM evolution. Using in silico approach and polymerase chain reaction, we obtained and characterized the messenger RNA sequences of ENAM in a frog, a lizard, and a crocodile; the genomic DNA sequences of ENAM in a frog and a lizard; and the putative sequence of chicken ENAM pseudogene. The comparison with mammalian ENAM sequences has revealed 1) the presence of an additional coding exon, named exon 8b, in sauropsids and marsupials, 2) a simpler 5'-untranslated region in nonmammalian ENAMs, 3) many sequence variations in the large exons while there are a few conserved regions in small exons, and 4) 25 amino acids that have been conserved during 350 million years of tetrapod evolution and hence of crucial biological importance. The chicken pseudogene was identified in a region that was not expected when considering the gene synteny in mammals. Together with the location of lizard ENAM in a homologous region, this result indicates that enamel genes were probably translocated in an ancestor of the sauropsid lineage. This study supports the origin of ENAM earlier in vertebrate evolution, confirms that tooth loss in modern birds led to the invalidation of enamel genes, and adds information on the important role played by, for example, the phosphorylated serines and the glycosylated asparagines for correct ENAM functions.

No abstract

Calcitonin (CT) is one of the hormones involved in vertebrate calcium regulation. It has been proposed to act as a hypocalcemic factor, but the regulatory pathways remain to be clarified. We investigated the CT/calcitonin gene-related peptide (CGRP) family in zebrafish and its potential involvement in calcium homeostasis. We identified the presence of four receptors: CTR, CRLR1, CRLR2, and CRLR3. From the phylogenetic analysis, together with the effect observed after CT and CGRP overexpression, we concluded that CTR appears to be a CT receptor and CRLR1 a CGRP receptor. The distribution of these two receptors shows a major presence in the central nervous system and in tissues involved in ionoregulation. Zebrafish embryos kept in high-Ca 2þ -concentration medium showed upregulation of CT and CTR expression and downregulation of the epithelial calcium channel (ECaC). Embryos injected with CT morpholino (CALC MO) incubated in high-Ca 2þ medium, showed downregulation of CTR together with upregulation on ECaC mRNA expression. In contrast, overexpression of CT cRNA induced the downregulation of ECaC mRNA synthesis, concomitant with the downregulation in the calcium content after 30 hours postfertilization. At 4 days postfertilization, CT cRNA injection induced upregulation of hypercalcemic factors, with subsequent increase in the calcium content. These results suggest that CT acts as a hypocalcemic factor in calcium regulation, probably through inhibition of ECaC synthesis. ß

Genetic approaches have shown that several genes could modify caries susceptibility; AmelogeninX (AMELX) has been repeatedly designated. Here, we hypothesized that AMELX mutations resulting in discrete changes of enamel microstructure may be found in children with a severe caries phenotype. In parallel, possible AMELX mutations that could explain resistance to caries may be found in caries-free patients. In this study, coding exons of AMELX and exon-intron boundaries were sequenced in 399 individuals with extensive caries (250) or caries-free (149) individuals from nine French hospital groups. No mutation responsible for a direct change of amelogenin function was identified. Seven single-nucleotide polymorphisms (SNPs) were found, 3 presenting a high allele frequency, and 1 being detected for the first time. Three SNPs were located in coding regions, 2 of them being non-synonymous. Both evolutionary and statistical analyses showed that none of these SNPs was associated with caries susceptibility, suggesting that AMELX is not a gene candidate in our studied population.

The neuropeptide gonadotropin-releasing hormone (GnRH) plays an important role in the control of reproductive functions. Vertebrates possess multiple GnRH forms that are classified into three main groups, namely GnRH1, GnRH2, and GnRH3. In order to gain more insights into the GnRH gene family in vertebrates, we sought to identify which paralogs of this family are present in cartilaginous fish. For this purpose, we searched the genomes and/or transcriptomes of three representative species of this group, the small-spotted catshark, Scyliorhinus canicula, the whale shark, Rhincodon typus and the elephant shark Callorhinchus milii. In each species, we report the identification of three GnRH genes. In catshark and whale shark, phylogenetic and synteny analysis showed that these three genes correspond to GnRH1, GnRH2, and GnRH3. In both species, GnRH1 was found to encode a novel form of GnRH whose primary structure was determined as follows: QHWSFDLRPG. In elephant shark, the three genes correspond to GnRH1a and GnRH1b, two copies of the GnRH1 gene, plus GnRH2. 3D structure prediction of the chondrichthyan GnRH-associated peptides (GAPs) revealed that catshark GAP1, GAP2, and elephant shark GAP2 peptides exhibit a helix-loop-helix (HLH) structure. This structure observed for many osteichthyan GAP1 and GAP2, may convey GAP biological activity. This HLH structure could not be observed for elephant shark GAP1a and GAP1b. As for all other GAP3 described so far, no typical 3D HLH structure was observed for catshark nor whale shark GAP3. RT-PCR analysis revealed that GnRH1, GnRH2, and GnRH3 genes are differentially expressed in the catshark brain. GnRH1 mRNA appeared predominant in the diencephalon while GnRH2 and GnRH3 mRNAs seemed to be most abundant in the mesencephalon and telencephalon, respectively. Taken together, our results show that the GnRH gene repertoire of the vertebrate ancestor was entirely conserved in the chondrichthyan lineage but that the GnRH3 gene was probably lost in holocephali. They also suggest that the three GnRH neuronal systems previously described in the brain of bony vertebrates are also present in cartilaginous fish.

Background We previously reported that impaired type I IFN activity, due to inborn errors of TLR3- and TLR7-dependent type I interferon (IFN) immunity or to autoantibodies against type I IFN, account for 15–20% of cases of life-threatening COVID-19 in unvaccinated patients. Therefore, the determinants of life-threatening COVID-19 remain to be identified in ~ 80% of cases. Methods We report here a genome-wide rare variant burden association analysis in 3269 unvaccinated patients with life-threatening COVID-19, and 1373 unvaccinated SARS-CoV-2-infected individuals without pneumonia. Among the 928 patients tested for autoantibodies against type I IFN, a quarter (234) were positive and were excluded. Results No gene reached genome-wide significance. Under a recessive model, the most significant gene with at-risk variants was TLR7, with an OR of 27.68 (95%CI 1.5–528.7, P = 1.1 × 10−4) for biochemically loss-of-function (bLOF) variants. We replicated the enrichment in rare predicted LOF (pLOF) variants at 13 influenza susceptibility loci involved in TLR3-dependent type I IFN immunity (OR = 3.70[95%CI 1.3–8.2], P = 2.1 × 10−4). This enrichment was further strengthened by (1) adding the recently reported TYK2 and TLR7 COVID-19 loci, particularly under a recessive model (OR = 19.65[95%CI 2.1–2635.4], P = 3.4 × 10−3), and (2) considering as pLOF branchpoint variants with potentially strong impacts on splicing among the 15 loci (OR = 4.40[9%CI 2.3–8.4], P = 7.7 × 10−8). Finally, the patients with pLOF/bLOF variants at these 15 loci were significantly younger (mean age [SD] = 43.3 [20.3] years) than the other patients (56.0 [17.3] years; P = 1.68 × 10−5). Conclusions Rare variants of TLR3- and TLR7-dependent type I IFN immunity genes can underlie life-threatening COVID-19, particularly with recessive inheritance, in patients under 60 years old.