A Noncoding A-to-U Kozak Site Change Related to the High Transmissibility of Alpha, Delta, and Omicron VOCs.

Three prevalent SARS-CoV-2 variants of concern (VOCs) emerged and caused epidemic waves. It is essential to uncover advantageous mutations that cause the high transmissibility of VOCs. However, viral mutations are tightly linked, so traditional population genetic methods, including machine learning-based methods, cannot reliably detect mutations conferring a fitness advantage. In this study, we developed an approach based on the sequential occurrence order of mutations and the accelerated furcation rate in the pandemic-scale phylogenomic tree. We analyzed 3,777,753 high-quality SARS-CoV-2 genomic sequences and the epidemiology metadata using the Coronavirus GenBrowser. We found that two noncoding mutations at the same position (g.a28271-/u) may be crucial to the high transmissibility of Alpha, Delta, and Omicron VOCs although the noncoding mutations alone cannot increase viral transmissibility. Both mutations cause an A-to-U change at the core position -3 of the Kozak sequence of the N gene and significantly reduce the protein expression ratio of ORF9b to N. Using a convergent evolutionary analysis, we found that g.a28271-/u, S:p.P681H/R, and N:p.R203K/M occur independently on three VOC lineages, suggesting that coordinated changes of S, N, and ORF9b proteins are crucial to high viral transmissibility. Our results provide new insights into high viral transmissibility co-modulated by advantageous noncoding and nonsynonymous changes.

Coronavirus GenBrowser for monitoring the transmission and evolution of SARS-CoV-2.

Genomic epidemiology is important to study the COVID-19 pandemic, and more than two million severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic sequences were deposited into public databases. However, the exponential increase of sequences invokes unprecedented bioinformatic challenges. Here, we present the Coronavirus GenBrowser (CGB) based on a highly efficient analysis framework and a node-picking rendering strategy. In total, 1,002,739 high-quality genomic sequences with the transmission-related metadata were analyzed and visualized. The size of the core data file is only 12.20 MB, highly efficient for clean data sharing. Quick visualization modules and rich interactive operations are provided to explore the annotated SARS-CoV-2 evolutionary tree. CGB binary nomenclature is proposed to name each internal lineage. The pre-analyzed data can be filtered out according to the user-defined criteria to explore the transmission of SARS-CoV-2. Different evolutionary analyses can also be easily performed, such as the detection of accelerated evolution and ongoing positive selection. Moreover, the 75 genomic spots conserved in SARS-CoV-2 but non-conserved in other coronaviruses were identified, which may indicate the functional elements specifically important for SARS-CoV-2. The CGB was written in Java and JavaScript. It not only enables users who have no programming skills to analyze millions of genomic sequences, but also offers a panoramic vision of the transmission and evolution of SARS-CoV-2.

The protease-sensitive N-terminal polybasic region of prion protein modulates its conversion to the pathogenic prion conformer.

Conversion of normal prion protein (PrPC) to the pathogenic PrPSc conformer is central to prion diseases such as Creutzfeldt-Jakob disease and scrapie; however, the detailed mechanism of this conversion remains obscure. To investigate how the N-terminal polybasic region of PrP (NPR) influences the PrPC-to-PrPSc conversion, we analyzed two PrP mutants: ΔN6 (deletion of all six amino acids in NPR) and Met4-1 (replacement of four positively charged amino acids in NPR with methionine). We found that ΔN6 and Met4-1 differentially impacted the binding of recombinant PrP (recPrP) to the negatively charged phospholipid 1-palmitoyl-2-oleoylphosphatidylglycerol, a nonprotein cofactor that facilitates PrP conversion. Both mutant recPrPs were able to form recombinant prion (recPrPSc) in vitro, but the convertibility was greatly reduced, with ΔN6 displaying the lowest convertibility. Prion infection assays in mammalian RK13 cells expressing WT or NPR-mutant PrPs confirmed these differences in convertibility, indicating that the NPR affects the conversion of both bacterially expressed recPrP and post-translationally modified PrP in eukaryotic cells. We also found that both WT and mutant recPrPSc conformers caused prion disease in WT mice with a 100% attack rate, but the incubation times and neuropathological changes caused by two recPrPSc mutants were significantly different from each other and from that of WT recPrPSc. Together, our results support that the NPR greatly influences PrPC-to-PrPSc conversion, but it is not essential for the generation of PrPSc. Moreover, the significant differences between ΔN6 and Met4-1 suggest that not only charge but also the identity of amino acids in NPR is important to PrP conversion.

Fine human genetic map based on UK10K data set.

Recombination is a major force that shapes genetic diversity. Determination of recombination rate is important and can theoretically be improved by increasing the sample size. However, it is nearly impossible to estimate recombination rates using traditional population genetics methods when the sample size is large because these methods are highly computationally demanding. In this study, we used a refined machine learning approach to estimate the recombination rate of the human genome using the UK10K human genomic dataset with 7,562 genomic sequences and its three subsets with 200, 400 and 2,000 genomic sequences. The estimation was performed under the human Out-of-Africa demographic model. We not only obtained an accurate human genetic map, but also found that the fluctuation of estimated recombination rate is reduced along the human genome when the sample size increases. The estimated UK10K recombination rate heterogeneity is less than that estimated from its subsets. Our results demonstrate how the sample size affects the estimated recombination rate, and analyses of a larger number of genomes result in a more precise estimation of recombination rate. The accurate genetic map based on UK10K data set is also expected to benefit other human biology researches.

Adamts18 deficiency in zebrafish embryo causes defective trunk angiogenesis and caudal vein plexus formation.

ADAMTS (A Disintegrin and Metalloproteinase with Thrombospondin type I motifs) enzymes play an important role in various morphogenesis processes. To determine the functions of Adamts18 in the early stages of organogenesis, we created Adamts18 deficient zebrafish using morpholino antisense oligonucleotides (MO) to generate exon 3 skipped adamts18 mRNA transcripts. Results showed that Adamts18 deficiency in zebrafish embryos caused developmental defects, including expanded brain ventricle and hindbrain edema, eye defects, and accumulation of blood in the caudal vein. Adamts18 deficiency also led to impaired trunk angiogenesis and formation of the caudal vein plexus (CVP). Consequently, Adamts18 deficient zebrafish embryos exhibited incomplete formation of intersegment vessels (ISVs), disruption of the honeycomb structure of CVP, and reduced CVP area and loop number. Furthermore, Adamts18 deficiency resulted in impaired blood circulation in major trunk, caudal vein (CV), and common cardinal vein (CCV). These aberrant vascular phenotypes in mutant zebrafish embryos were shown to be associated with a decreased expression of multiple angiogenesis-related signaling genes, including slit/robo, dll4/Notch, cox2, and fgfr. These findings indicate the critical role of Adamts18 in the early stages of vascular network development.

Gut transcriptomic changes during hibernation in the greater horseshoe bat (Rhinolophus ferrumequinum).

BACKGROUND: The gut is the major organ for nutrient absorption and immune response in the body of animals. Although effects of fasting on the gut functions have been extensively studied in model animals (e.g. mice), little is known about the response of the gut to fasting in a natural condition (e.g. hibernation). During hibernation, animals endure the long term of fasting and hypothermia. RESULTS: Here we generated the first gut transcriptome in a wild hibernating bat (Rhinolophus ferrumequinum). We identified 1614 differentially expressed genes (DEGs) during four physiological states (Torpor, Arousal, Winter Active and Summer Active). Gene co-expression network analysis assigns 926 DEGs into six modules associated with Torpor and Arousal. Our results reveal that in response to the stress of luminal nutrient deficiency during hibernation, the gut helps to reduce food intake by overexpressing genes (e.g. CCK and GPR17) that regulate the sensitivity to insulin and leptin. At the same time, the gut contributes energy supply by overexpressing genes that increase capacity for ketogenesis (HMGCS2) and selective autophagy (TEX264). Furthermore, we identified separate sets of multiple DEGs upregulated in Torpor and Arousal whose functions are involved in innate immunity. CONCLUSION: This is the first gut transcriptome of a hibernating mammal. Our study identified candidate genes associated with regulation of food intake and enhance of innate immunity in the gut during hibernation. By comparing with previous studies, we found that two DEGs (CPE and HSPA8) were also significantly elevated during torpor in liver and brain of R. ferrumequinum and several DEGs (e.g. TXNIP and PDK1/4) were commonly upregulated during torpor in multiple tissues of different mammals. Our results support that shared expression changes may underlie the hibernation phenotype by most mammals.

A Disintegrin and Metalloproteinase with Thrombospondin Motifs 18 Deficiency Leads to Visceral Adiposity and Associated Metabolic Syndrome in Mice.

Visceral adiposity is of greater risk than obesity in s.c. adipose tissue for diabetes and cardiovascular disease. Its pathogenesis remains unclear, but it is associated with extracellular matrix (ECM) remodeling. A disintegrin and metalloproteinase with thrombospondin motifs (ADAMTSs) are a family of secreted zinc-dependent metalloproteinases that play crucial roles in development and various diseases because of their ECM remodeling activity. ADAMTS18 is an orphan ADAMTS whose function and substrate remain unclear. Herein, we showed that Adamts18 mRNA was abundantly expressed in visceral (gonadal) white adipose tissue (vWAT) during the early stage of development after birth. Adamts18 knockout (KO) mice showed increased body fat percentage and larger adipocyte size in vWAT relative to wild-type littermates. These findings may be partly attributed to ECM remodeling, especially increased expression of laminin 1 and adipokine thrombospondin 1 in vWAT. Attenuated extracellular signal-regulated kinase 1 and 2 activity, along with increased expression of adipocyte-specific transcription factors peroxisome proliferator-activated receptor-γ, CCAAT/enhancer binding protein ß, and marker gene Fabp4, was detected in vWAT of Adamts18 KO mice. Furthermore, Adamts18 KO mice showed early metabolic syndrome, including hyperlipidemia, blood glucose metabolic disorder, and hypertension. ADAMTS18 deficiency promotes atherosclerosis in apolipoprotein E-deficient mice. These results indicate a novel function of ADAMTS18 in vWAT development and associated metabolic disorders.

The Genomes of Two Bat Species with Long Constant Frequency Echolocation Calls.

Bats can perceive the world by using a wide range of sensory systems, and some of the systems have become highly specialized, such as auditory sensory perception. Among bat species, the Old World leaf-nosed bats and horseshoe bats (rhinolophoid bats) possess the most sophisticated echolocation systems. Here, we reported the whole-genome sequencing and de novo assembles of two rhinolophoid bats-the great leaf-nosed bat (Hipposideros armiger) and the Chinese rufous horseshoe bat (Rhinolophus sinicus). Comparative genomic analyses revealed the adaptation of auditory sensory perception in the rhinolophoid bat lineages, probably resulting from the extreme selectivity used in the auditory processing by these bats. Pseudogenization of some vision-related genes in rhinolophoid bats was observed, suggesting that these genes have undergone relaxed natural selection. An extensive contraction of olfactory receptor gene repertoires was observed in the lineage leading to the common ancestor of bats. Further extensive gene contractions can be observed in the branch leading to the rhinolophoid bats. Such concordance suggested that molecular changes at one sensory gene might have direct consequences for genes controlling for other sensory modalities. To characterize the population genetic structure and patterns of evolution, we re-sequenced the genome of 20 great leaf-nosed bats from four different geographical locations of China. The result showed similar sequence diversity values and little differentiation among populations. Moreover, evidence of genetic adaptations to high altitudes in the great leaf-nosed bats was observed. Taken together, our work provided a useful resource for future research on the evolution of bats.

Adaptation of the FK506 binding protein 1B to hibernation in bats.

Hibernation is an adaptive strategy used by some animals to cope with cold and food shortage. The heart rate, overall energy need, body temperature, and many other physiological functions are greatly reduced during torpor but promptly return to normal levels upon arousal. The heartbeat of torpid bats can be hundreds fold lower than that of active bats, indicating that hibernating bats have a remarkable ability to control excitation-contraction coupling in cardiac muscle. FKBP1B (calstabin 2), a peptidyl-prolyl cis-trans isomerase, is critical for the regulation of excitation-contraction coupling. Whether FKBP1B is adapted to hibernation in bats is not known. Evolutionary analyses showed that the ω values of the Fkbp1b genes of 25 mammalian species are all less than 1, and amino acid sequence alignments revealed that FKBP1B proteins are highly conserved in mammals. The expression of the Fkbp1b gene was found to be elevated at both mRNA and protein levels in two distantly related bats (Rhinolophus ferrumequinum in Yinpterochiroptera and Myotis ricketti in Yangochiroptera) during torpor. Transcription factors such as YY1 and SPs were bioinformatically determined to have a higher binding affinity to the potential regulatory regions of Fkbp1b genes in hibernating than in non-hibernating mammals. This study provides new insights into the molecular evolution of Fkbp1b in adaptation to bat hibernation.

Maintenance of neural activities in torpid Rhinolophus ferrumequinum bats revealed by 2D gel-based proteome analysis.

Bats are the only mammals capable of self-powered flying. Many bat species hibernate in winter. A reversible control of cerebral activities is critical for bats to accommodate a repeated torpor-arousal cycle during hibernation. Little is known about the molecular mechanisms that regulate neuronal activities in torpid bats. In this study, Rhinolophus ferrumequinum bat brain proteins were fractionated, and their abundance in active and torpid states was compared. Results of 2D gel-based proteomics showed that 38% of identified proteins with a significant change in abundance are involved in synaptic vesicle recycling and cytoskeletal integrity. Changes in the abundance of proteins related to RNA splicing, proteostasis, redox homeostasis, mitochondrial function, and energy metabolism were also detected. In addition, the levels of GNAO1 (guanine nucleotide-binding protein Gαo subunit), an important modulator of neuronal transmembrane signaling, were significantly increased in the insoluble protein fraction of torpid bats; this may be due to GNAO1 palmitoylation making it insoluble. Our data provide molecular evidence for the maintenance of neuronal activities in torpid bats and suggest that a reversible palmitoylation of the G protein plays a role in the regulation of neuronal activities during bat hibernation.

Adaptation of peroxisome proliferator-activated receptor alpha to hibernation in bats.

BACKGROUND: Hibernation is a survival mechanism in the winter for some animals. Fat preserved instead of glucose produced is the primary fuel during winter hibernation of mammals. Many genes involved in lipid metabolism are regulated by the peroxisome proliferator-activated receptor alpha (PPARα). The role of PPARα in hibernation of mammals remains largely unknown. Using a multidisciplinary approach, we investigated whether PPARα is adapted to hibernation in bats. RESULTS: Evolutionary analyses revealed that the ω value of Pparα of the ancestral lineage of hibernating bats in both Yinpterochiroptera and Yangochiroptera was lower than that of non-hibernating bats in Yinpterochiroptera, suggesting that a higher selective pressure acts on Pparα in hibernating bats. PPARα expression was found to be increased at both mRNA and protein levels in distantly related bats (Rhinolophus ferrumequinum and Hipposideros armiger in Yinpterochiroptera and Myotis ricketti in Yangochiroptera) during their torpid episodes. Transcription factors such as FOXL1, NFYA, NFYB, SP1, TBP, and ERG were bioinformatically determined to have a higher binding affinity to the potential regulatory regions of Pparα in hibernating than in non-hibernating mammals. Genome-wide bioinformatic analyses of 64 mammalian species showed that PPARα has more potential target genes and higher binding affinity to these genes in hibernating than in non-hibernating mammals. CONCLUSIONS: We conclude that PPARα is adapted to hibernation in bats based on the observations that Pparα has a more stringent functional constraint in the ancestral lineage of hibernating bats and a higher level of expression in hibernating than in non-hibernating bats. We also conclude that PPARα plays a very important role in hibernation as hibernators have more PPARα target genes than non-hibernators, and PPARα in hibernators has a higher binding affinity for its target genes than in non-hibernators.

A non-image-forming visual circuit mediates the innate fear of heights in male mice.

The neural basis of fear of heights remains largely unknown. In this study, we investigated the fear response to heights in male mice and observed characteristic aversive behaviors resembling human height vertigo. We identified visual input as a critical factor in mouse reactions to heights, while peripheral vestibular input was found to be nonessential for fear of heights. Unexpectedly, we found that fear of heights in naïve mice does not rely on image-forming visual processing by the primary visual cortex. Instead, a subset of neurons in the ventral lateral geniculate nucleus (vLGN), which connects to the lateral/ventrolateral periaqueductal gray (l/vlPAG), drives the expression of fear associated with heights. Additionally, we observed that a subcortical visual pathway linking the superior colliculus to the lateral posterior thalamic nucleus inhibits the defensive response to height threats. These findings highlight a rapid fear response to height threats through a subcortical visual and defensive pathway from the vLGN to the l/vlPAG.

Genomic inference of a severe human bottleneck during the Early to Middle Pleistocene transition.

Population size history is essential for studying human evolution. However, ancient population size history during the Pleistocene is notoriously difficult to unravel. In this study, we developed a fast infinitesimal time coalescent process (FitCoal) to circumvent this difficulty and calculated the composite likelihood for present-day human genomic sequences of 3154 individuals. Results showed that human ancestors went through a severe population bottleneck with about 1280 breeding individuals between around 930,000 and 813,000 years ago. The bottleneck lasted for about 117,000 years and brought human ancestors close to extinction. This bottleneck is congruent with a substantial chronological gap in the available African and Eurasian fossil record. Our results provide new insights into our ancestry and suggest a coincident speciation event.

Progressive pseudogenization: vitamin C synthesis and its loss in bats.

For the past 50 years, it was believed that all bats, like humans and guinea pigs, did not synthesize vitamin C (Vc) because they lacked activity of L-gulonolactone oxidase (GULO) in their livers. Humans and guinea pigs lack the activity due to pseudogenization of GULO in their genomes, but there is no genetic evidence to show whether such loss in bats is caused by pseudogenization. Unexpectedly, our successful molecular cloning in one frugivorous bat (Rousettus leschenaultii) and one insectivorous bat (Hipposideros armiger) ascertains that no pseudogenization occurs in these species. Furthermore, we find normal GULO protein expression using bat-specific anti-GULO polyclonal antibodies in bats, evaluated by Western blotting. Most surprisingly, GULO activity assays reveal that these two bat species have retained the ability to synthesize Vc, but at low levels compared with the mouse. It is known that bats in the genus Pteropus have lost GULO activity. We then found that functional constraints acting on the GULO of Pteropus vampyrus (which lost its function) are relaxed. These results imply that the ability to synthesize Vc in bats has not been lost completely in species as previously thought. We also suggest that the evolution of bat GULO genes can be a good model to study genetic processes associated with loss-of-function.

Association of CDH11 with Autism Spectrum Disorder Revealed by Matched-gene Co-expression Analysis and Mouse Behavioral Studies.

A large number of putative risk genes for autism spectrum disorder (ASD) have been reported. The functions of most of these susceptibility genes in developing brains remain unknown, and causal relationships between their variation and autism traits have not been established. The aim of this study was to predict putative risk genes at the whole-genome level based on the analysis of gene co-expression with a group of high-confidence ASD risk genes (hcASDs). The results showed that three gene features - gene size, mRNA abundance, and guanine-cytosine content - affect the genome-wide co-expression profiles of hcASDs. To circumvent the interference of these features in gene co-expression analysis, we developed a method to determine whether a gene is significantly co-expressed with hcASDs by statistically comparing the co-expression profile of this gene with hcASDs to that of this gene with permuted gene sets of feature-matched genes. This method is referred to as "matched-gene co-expression analysis" (MGCA). With MGCA, we demonstrated the convergence in developmental expression profiles of hcASDs and improved the efficacy of risk gene prediction. The results of analysis of two recently-reported ASD candidate genes, CDH11 and CDH9, suggested the involvement of CDH11, but not CDH9, in ASD. Consistent with this prediction, behavioral studies showed that Cdh11-null mice, but not Cdh9-null mice, have multiple autism-like behavioral alterations. This study highlights the power of MGCA in revealing ASD-associated genes and the potential role of CDH11 in ASD.

Evaluation of Three Sample Introduction Systems for Impurity Analysis of an Ultrapure Reagent Using a Scanning Mobility Particle Sizer.

The semiconductor industry continues to shrink the device sizes while applying more complex shapes and using diverse materials, which requires parallel improvements in the quality of ultrapure reagents. The need for ultrapure reagents has led to ever-higher demands for the performance of analytical instruments used to detect ultratrace impurities. In this study, nonvolatile impurities in ultrapure reagents were quantified using a scanning mobility particle sizer (SMPS). The performances of three different sample introduction systems, i.e., an electrospray (ES), an aerosol generator with a heating chamber and a Nafion desolvation membrane (NB-II), and a MicroMist nebulizer with a heated cyclonic spray chamber and a three-stage Peltier-cooled desolvation system (MM-APEX), were evaluated for the lower limit of detection of a SMPS. The MM-APEX equipped with the SMPS was able to detect NaCl additives at a concentration of 100 parts per trillion (ppt, ng/L) in ultrapure water, which was approximately 104- and 102-fold lower than those of ES and NB-II, respectively. The practical application of MM-APEX with the SMPS for commercial isopropanol samples was also studied. The results clearly demonstrate that the impurity concentrations presented by the NaCl-equivalent concentrations among different sources of isopropanol were at the ppt to parts-to-billion (ppb) scale. The SMPS system equipped with MM-APEX is capable of recognizing impurities with concentrations ranging from tens ppt to thousands of parts per million (ppm), which is beneficial for an ultratrace analysis of nonvolatile impurities in semiconductor process chemicals.

Unexpected expression of heat-activated transient receptor potential (TRP) channels in winter torpid bats and cold-activated TRP channels in summer active bats.

The ability to sense temperature changes is crucial for mammalian survival. Mammalian thermal sensing is primarily carried out by thermosensitive transient receptor potential channels (Thermo-TRPs). Some mammals hibernate to survive cold winter conditions, during which time their body temperature fluctuates dramatically. However, the underlying mechanisms by which these mammals regulate thermal responses remain unclear. Using quantitative real-time polymerase chain reaction (qRT-PCR) and the Western blotting, we found that Myotis ricketti bats had high levels of heat-activated TRPs (e.g., TRPV1 and TRPV4) during torpor in winter and cold-activated TRPs (e.g., TRPM8 and TRPC5) during active states in summer. We also found that laboratory mice had high mRNA levels of cold-activated TRPs (e.g., Trpm8 and Trpc5) under relatively hot conditions (i.e., 40 °C). These data suggest that small mammals up-regulate the expression of cold-activated TRPs even under warm or hot conditions. Binding site analysis showed that some homeobox (HOX) transcription factors (TFs) regulate the expression of hot- and cold-activated TRP genes and that some TFs of the Pit-Oct-Unc (POU) family regulate warm-sensitive and cold-activated TRP genes. The dual-luciferase reporter assay results demonstrated that TFs HOXA9, POU3F1, and POU5F1 regulate TRPC5 expression, suggesting that Thermo-TRP genes are regulated by multiple TFs of the HOX and POU families at different levels. This study provides insights into the adaptive mechanisms underlying thermal sensing used by bats to survive hibernation.

ADAMTS18 regulates vaginal opening through influencing the fusion of Mullerian duct and apoptosis of vaginal epithelial cells in mice.

The ADAMTS (A Disintegrin and Metalloproteinase with Thrombospondin Motifs) enzymes are secreted metalloproteinases with major roles in development, morphogenesis, and tissue repair via the assembly and degradation of extracellular matrix (ECM). In this study, we investigated the role of ADAMTS18 in the development of the reproductive tract in female mice by phenotyping Adamts18 knockout (Adamts18-/-) mice. The results showed that Adamst18 mRNAs were abundantly expressed in vaginal epithelial cells and muscularis cells of the developing vagina. At the time of vaginal opening (5 weeks of age), about 41 % of Adamts18-/- females showed enlarged protrusions in the upper and middle parts of the vagina, reduced vaginal length, and simultaneously exhibited vaginal atresia. 6% Adamts18-/- females exhibited vaginal septum. Histological analyses revealed that the paired Mullerian ducts in ∼33 % female Adamts18-/- embryos failed to fuse at embryonic day 15.5 (E15.5) resulting in the formation of two vaginal cavities. Results of TUNEL assay and immunohistochemistry for caspase-3 showed that the number of apoptotic cells in the terminal portion of the vagina of 5-week-old Adamts18-/- females with vaginal atresia was significantly decreased. Adamts18-/- females also showed a significant decrease in serum estradiol E2 compared to age-matched Adamts18+/+ females. Results of qRT-PCR showed that the expression level of the anti-apoptosis gene Bcl-2 was significantly increased and that of the apoptosis-related gene Epha1 was decreased in the vagina of 5-week-old Adamts18-/- females. These results suggest that ADAMTS18 regulates vaginal opening through influencing the fusion of Mullerian ducts and apoptosis of vaginal cells in mice.

New exon and accelerated evolution of placental gene Nrk occurred in the ancestral lineage of placental mammals.

INTRODUCTION: The chorioallantoic placenta is a specific organ for placental mammals. However, the adaptive events during its emergence are still poorly investigated. METHODS: We scanned the chromosome X to detect the accelerated evolution in the ancestral lineage of placental mammals, and constructed 3D protein structure models of a candidate by homology modeling. RESULTS: Eight branch-specific accelerated regions were identified. Five of these regions (P=5.61×10-11 ~ 9.03×10-8) are located in the five exons of Nik-related kinase (Nrk), which is essential in placenta development and fetoplacental induction of labor. Nrk belongs to the germinal center kinase-IV subfamily with the overall similar protein structure; however, a new exon emerged in ancestors of placental mammals and its sequence has been conserved since then. Structure modelling of NRK suggests that the accelerated exons and the placental-mammal-specific exon (as a new loop) could change the enzymatic activity and the structure of placental mammal NRK. DISCUSSION: Since the new loop is surrounded by the accelerated protein regions, it is likely that the new loop occurred and shifted the function of NRK, and then the accelerated evolution of Nrk occurred to adapt the structure change caused by the new loop in the ancestral lineage of placental mammals. Overall, this work suggests that the fundamental process of placental development and fetoplacental induction of labor has been targeted by positive Darwinian selection.

Adamts18 modulates the development of the aortic arch and common carotid artery.

Members of a disintegrin and metalloproteinases with thrombospondin motif (ADAMTS) family have been implicated in various vascular diseases. However, their functional roles in early embryonic vascular development are unknown. In this study, we showed that Adamts18 is highly expressed at E11.5-E14.5 in cells surrounding the embryonic aortic arch (AOAR) and the common carotid artery (CCA) during branchial arch artery development in mice. Adamts18 deficiency was found to cause abnormal development of AOAR, CCA, and the third and fourth branchial arch appendages, leading to hypoplastic carotid body, thymus, and variation of middle cerebral artery. Adamts18 was shown to affect the accumulation of extracellular matrix (ECM) components, in particular fibronectin (Fn), around AOAR and CCA. As a result of increased Fn accumulation, the Notch3 signaling pathway was activated to promote the differentiation of cranial neural crest cells (CNCCs) to vascular smooth muscle cells. These data indicate that Adamts18-mediated ECM homeostasis is crucial for the differentiation of CNCCs.