OrthoRefine: automated enhancement of prior ortholog identification via synteny.

BACKGROUND: Identifying orthologs continues to be an early and imperative step in genome analysis but remains a challenging problem. While synteny (conservation of gene order) has previously been used independently and in combination with other methods to identify orthologs, applying synteny in ortholog identification has yet to be automated in a user-friendly manner. This desire for automation and ease-of-use led us to develop OrthoRefine, a standalone program that uses synteny to refine ortholog identification. RESULTS: We developed OrthoRefine to improve the detection of orthologous genes by implementing a look-around window approach to detect synteny. We tested OrthoRefine in tandem with OrthoFinder, one of the most used software for identification of orthologs in recent years. We evaluated improvements provided by OrthoRefine in several bacterial and a eukaryotic dataset. OrthoRefine efficiently eliminates paralogs from orthologous groups detected by OrthoFinder. Using synteny increased specificity and functional ortholog identification; additionally, analysis of BLAST e-value, phylogenetics, and operon occurrence further supported using synteny for ortholog identification. A comparison of several window sizes suggested that smaller window sizes (eight genes) were generally the most suitable for identifying orthologs via synteny. However, larger windows (30 genes) performed better in datasets containing less closely related genomes. A typical run of OrthoRefine with ~ 10 bacterial genomes can be completed in a few minutes on a regular desktop PC. CONCLUSION: OrthoRefine is a simple-to-use, standalone tool that automates the application of synteny to improve ortholog detection. OrthoRefine is particularly efficient in eliminating paralogs from orthologous groups delineated by standard methods.

Biochemical functions and structure of Caenorhabditis elegans ZK177.8 protein: Aicardi-Goutières syndrome SAMHD1 dNTPase ortholog.

Mutations in sterile alpha motif domain and histidine-aspartate domain-containing protein 1 (SAMHD1) are found in a neurodevelopmental disorder, Aicardi-Goutières syndrome, and cancers, and SAMHD1, which is a deoxynucleoside triphosphate (dNTP) triphosphorylase, was identified as a myeloid-specific HIV-1 restriction factor. Here, we characterized the enzymology and structure of an SAMHD1 ortholog of Caenorhabditis elegans, ZK177.8, which also reportedly induces developmental defects upon gene knockdown. We found ZK177.8 protein is a dNTPase allosterically regulated by dGTP. The active site of ZK177.8 recognizes both 2' OH and triphosphate moieties of dNTPs but not base moiety. The dGTP activator induces the formation of the enzymatically active ZK177.8 tetramers, and ZK177.8 protein lowers cellular dNTP levels in a human monocytic cell line. Finally, ZK177.8 tetramers display very similar X-ray crystal structure with human and mouse SAMHD1s except that its lack of the canonical sterile alpha motif domain. This striking conservation in structure, function, and allosteric regulatory mechanism for the hydrolysis of the DNA building blocks supports their host developmental roles.

A Chicken Tapasin ortholog can chaperone empty HLA-B∗37:01 molecules independent of other peptide-loading components.

Human Tapasin (hTapasin) is the main chaperone of MHC-I molecules, enabling peptide loading and antigen repertoire optimization across HLA allotypes. However, it is restricted to the endoplasmic reticulum (ER) lumen as part of the protein loading complex (PLC), and therefore is highly unstable when expressed in recombinant form. Additional stabilizing co-factors such as ERp57 are required to catalyze peptide exchange in vitro, limiting uses for the generation of pMHC-I molecules of desired antigen specificities. Here, we show that the chicken Tapasin (chTapasin) ortholog can be expressed recombinantly at high yields in a stable form, independent of co-chaperones. chTapasin can bind the human HLA-B∗37:01 with low micromolar-range affinity to form a stable tertiary complex. Biophysical characterization by methyl-based NMR methods reveals that chTapasin recognizes a conserved ß2m epitope on HLA-B∗37:01, consistent with previously solved X-ray structures of hTapasin. Finally, we provide evidence that the B∗37:01/chTapasin complex is peptide-receptive and can be dissociated upon binding of high-affinity peptides. Our results highlight the use of chTapasin as a stable scaffold for protein engineering applications aiming to expand the ligand exchange function on human MHC-I and MHC-like molecules.

Inferring complex evolutionary history of the closely related East Asian wild roses in Rosa sect. Synstylae (Rosaceae) based on genomic evidence from conserved orthologs.

BACKGROUND AND AIMS: The section Synstylae in genus Rosa (Rosaceae) comprises 25-36 species and includes several major progenitors of modern rose cultivars. East Asian Synstylae species have recently diverged and are closely related, and their phylogenetic relationships remain unclear. In the present study, we employed a conserved ortholog set (COS) markers and genome-wide nuclear orthologs to elucidate their phylogenetic relationships and unravel their complex evolutionary history. METHODS: Utilising on eight Rosaceae COS (RosCOS) markers, we analysed a total of 137 accessions representing 15 East Asian Synstylae taxa to establish a robust phylogenetic framework and reconstruct ancestral areas. Furthermore, we constructed the species tree for eight representative species and estimated their divergence times based on 1,683 genome-wide orthologs. The species tree-gene tree coalescence time comparison, Patterson's D, f4-ratio, and f-branch statistics were analysed to identify incomplete lineage sorting (ILS), genetic introgression, and reticulation events using conserved ortholog data. KEY RESULTS: RosCOS markers and genome-wide orthologs effectively resolved the robust phylogeny of East Asian Rosa sect. Synstylae. Species divergence times estimated with genome-wide orthologs indicated that East Asian Synstylae species have recently diverged, with an estimated crown age of approximately 2 Mya. The rampant gene tree discordance indicated the possibility of ILS and/or genetic introgression. In the section Synstylae, deeper coalescence in the gene trees compared to the species tree suggested ILS as a source of gene tree discordance. Further, Patterson's D and f-branch statistics indicated that several lineages in the section were involved in genetic introgression. CONCLUSIONS: We have unravelled the complex evolutionary history of East Asian Rosa sect. Synstylae, including recent species divergences, ILS, and genetic introgression. Coupled with the geographical and ecological complexity of East Asia, ILS and genetic introgression may have contributed to the rapid diversification of East Asian Synstylae species by permitting adaptation to diverse environments.

Depletion of lipid storage droplet-1 delays endoreplication progression and induces cell death in Drosophila salivary gland.

Perilipins are evolutionarily conserved from insects to mammals. Drosophila lipid storage droplet-1 (LSD-1) is a lipid storage droplet membrane surface-binding protein family member and a counterpart to mammalian perilipin 1 and is known to play a role in lipolysis. However, the function of LSD-1 during specific tissue development remains under investigation. This study demonstrated the role of LSD-1 in salivary gland development. Knockdown of Lsd-1 in the salivary gland was established using the GAL4/UAS system. The third-instar larvae of knockdown flies had small salivary glands containing cells with smaller nuclei. The null mutant Drosophila also showed the same phenotype. The depletion of LSD-1 expression induced a delay of endoreplication due to decreasing CycE expression and increasing DNA damage. Lsd-1 genetically interacted with Myc in the third-instar larvae. These results demonstrate that LSD-1 is involved in cell cycle and cell death programs in the salivary gland, providing novel insight into the effects of LSD-1 in regulating salivary gland development and the interaction between LSD-1 and Myc.

Large-scale analysis of putative Euphorbiaceae R2R3-MYB transcription factors identifies a MYB involved in seed oil biosynthesis.

BACKGROUND: MYB transcription factors are widely distributed in the plant kingdom and play key roles in regulatory networks governing plant metabolism and biochemical and physiological processes. RESULTS: Here, we first determined the R2R3-MYB genes in five Euphorbiaceae genomes. The three Trp (W) residues from the first MYB domain (R2) were absolutely conserved, whereas the first W residue from the second MYB domain (R3) was preferentially mutated. The R2R3-MYBs were clustered into 48 functional subfamilies, of which 34 had both R2R3-MYBs of Euphorbiaceae species and AtMYBs, and four contained only Euphorbiaceae R2R3-MYBs. The whole-genome duplication (WGD) and/or segmental duplication (SD) played key roles in the expansion of the R2R3-MYB family. Unlike paralogous R2R3-MYB family members, orthologous R2R3-MYB members contained a higher selective pressure and were subject to a constrained evolutionary rate. VfMYB36 was specifically expressed in fruit, and its trend was consistent with the change in oil content, indicating that it might be involved in oil biosynthesis. Overexpression experiments showed that VfMYB36 could significantly provide linolenic acid (C18:3) content, which eventually led to a significant increase in oil content. CONCLUSION: Our study first provides insight into understanding the evolution and expression of R2R3-MYBs in Euphorbiaceae species, and also provides a target for the production of biomass diesel and a convenient way for breeding germplasm resources with high linolenic acid content in the future.

Peptide signaling through leucine-rich repeat receptor kinases: insight into land plant evolution.

Multicellular organisms need mechanisms for communication between cells so that they can fulfill their purpose in the organism as a whole. Over the last two decades, several small post-translationally modified peptides (PTMPs) have been identified as components of cell-to-cell signaling modules in flowering plants. Such peptides most often influence growth and development of organs not universally conserved among land plants. PTMPs have been matched to subfamily XI leucine-rich repeat receptor-like kinases with > 20 repeats. Phylogenetic analyses, facilitated by recently published genomic sequences of non-flowering plants, have identified seven clades of such receptors with a history back to the common ancestor of bryophytes and vascular plants. This raises a number of questions: When did peptide signaling arise during land plant evolution? Have orthologous peptide-receptor pairs preserved their biological functions? Has peptide signaling contributed to major innovations, such as stomata, vasculature, roots, seeds, and flowers? Using genomic, genetic, biochemical, and structural data and non-angiosperm model species, it is now possible to address these questions. The vast number of peptides that have not yet found their partners suggests furthermore that we have far more to learn about peptide signaling in the coming decades.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins.

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.

The Interactive Complex between Cytomegalovirus Kinase vCDK/pUL97 and Host Factors CDK7-Cyclin H Determines Individual Patterns of Transcription in Infected Cells.

The infection of human cytomegalovirus (HCMV) is strongly determined by the host-cell interaction in a way that the efficiency of HCMV lytic replication is dependent on the regulatory interplay between viral and cellular proteins. In particular, the activities of protein kinases, such as cyclin-dependent kinases (CDKs) and the viral CDK ortholog (vCDK/pUL97), play an important role in both viral reproduction and virus-host interaction. Very recently, we reported on the complexes formed between vCDK/pUL97, human cyclin H, and CDK7. Major hallmarks of this interplay are the interaction between cyclin H and vCDK/pUL97, which is consistently detectable across various conditions and host cell types of infection, the decrease or increase in pUL97 kinase activity resulting from cyclin H knock-down or elevated levels, respectively, and significant trans-stimulation of human CDK7 activity by pUL97 in vitro. Due to the fact that even a ternary complex of vCDK/pUL97-cyclin H-CDK7 can be detected by coimmunoprecipitation and visualized by bioinformatic structural modeling, we postulated a putative impact of the respective kinase activities on the patterns of transcription in HCMV-infected cells. Here, we undertook a first vCDK/pUL97-specific transcriptomic analysis, which combined conditions of fully lytic HCMV replication with those under specific vCDK/pUL97 or CDK7 drug-mediated inhibition or transient cyclin H knockout. The novel results were further strengthened using bioinformatic modeling of the involved multi-protein complexes. Our data underline the importance of these kinase activities for the C-terminal domain (CTD) phosphorylation-driven activation of host RNA polymerase in HCMV-infected cells. The impact of the individual experimental conditions on differentially expressed gene profiles is described in detail and discussed.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae.

BACKGROUND: Lipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species. RESULTS: To uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree. CONCLUSION: Comparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

Special Care Is Needed in Applying Phylogenetic Comparative Methods to Gene Trees with Speciation and Duplication Nodes.

How gene function evolves is a central question of evolutionary biology. It can be investigated by comparing functional genomics results between species and between genes. Most comparative studies of functional genomics have used pairwise comparisons. Yet it has been shown that this can provide biased results, as genes, like species, are phylogenetically related. Phylogenetic comparative methods should be used to correct for this, but they depend on strong assumptions, including unbiased tree estimates relative to the hypothesis being tested. Such methods have recently been used to test the "ortholog conjecture," the hypothesis that functional evolution is faster in paralogs than in orthologs. Although pairwise comparisons of tissue specificity (τ) provided support for the ortholog conjecture, phylogenetic independent contrasts did not. Our reanalysis on the same gene trees identified problems with the time calibration of duplication nodes. We find that the gene trees used suffer from important biases, due to the inclusion of trees with no duplication nodes, to the relative age of speciations and duplications, to systematic differences in branch lengths, and to non-Brownian motion of tissue specificity on many trees. We find that incorrect implementation of phylogenetic method in empirical gene trees with duplications can be problematic. Controlling for biases allows successful use of phylogenetic methods to study the evolution of gene function and provides some support for the ortholog conjecture using three different phylogenetic approaches.

Phylogenetic and Expression Analysis of Fos Transcription Factors in Zebrafish.

Members of the FOS protein family regulate gene expression responses to a multitude of extracellular signals and are dysregulated in several pathological states. Whilst mouse genetic models have provided key insights into the tissue-specific functions of these proteins in vivo, little is known about their roles during early vertebrate embryonic development. This study examined the potential of using zebrafish as a model for such studies and, more broadly, for investigating the mechanisms regulating the functions of Fos proteins in vivo. Through phylogenetic and sequence analysis, we identified six zebrafish FOS orthologues, fosaa, fosab, fosb, fosl1a, fosl1b, and fosl2, which show high conservation in key regulatory domains and post-translational modification sites compared to their equivalent human proteins. During embryogenesis, zebrafish fos genes exhibit both overlapping and distinct spatiotemporal patterns of expression in specific cell types and tissues. Most fos genes are also expressed in a variety of adult zebrafish tissues. As in humans, we also found that expression of zebrafish FOS orthologs is induced by oncogenic BRAF-ERK signalling in zebrafish melanomas. These findings suggest that zebrafish represent an alternate model to mice for investigating the regulation and functions of Fos proteins in vertebrate embryonic and adult tissues, and cancer.

Cloning and Functional Characterization of Dog OCT1 and OCT2: Another Step in Exploring Species Differences in Organic Cation Transporters.

OCT1 and OCT2 are polyspecific membrane transporters that are involved in hepatic and renal drug clearance in humans and mice. In this study, we cloned dog OCT1 and OCT2 and compared their function to the human and mouse orthologs. We used liver and kidney RNA to clone dog OCT1 and OCT2. The cloned and the publicly available RNA-Seq sequences differed from the annotated exon-intron structure of OCT1 in the dog genome CanFam3.1. An additional exon between exons 2 and 3 was identified and confirmed by sequencing in six additional dog breeds. Next, dog OCT1 and OCT2 were stably overexpressed in HEK293 cells and the transport kinetics of five drugs were analyzed. We observed strong differences in the transport kinetics between dog and human orthologs. Dog OCT1 transported fenoterol with 12.9-fold higher capacity but 14.3-fold lower affinity (higher KM) than human OCT1. Human OCT1 transported ipratropium with 5.2-fold higher capacity but 8.4-fold lower affinity than dog OCT1. Compared to human OCT2, dog OCT2 showed 10-fold lower transport of fenoterol and butylscopolamine. In conclusion, the functional characterization of dog OCT1 and OCT2 reported here may have implications when using dogs as pre-clinical models as well as for drug therapy in dogs.

Arabidopsis RING-type E3 ubiquitin ligase XBAT35.2 promotes proteasome-dependent degradation of ACD11 to attenuate abiotic stress tolerance.

Plants employ multiple mechanisms to cope with a constantly changing and challenging environment, including using the ubiquitin proteasome system (UPS) to alter their proteome to assist in initiating, modulating and terminating responses to stress. We previously reported that the ubiquitin ligase XBAT35.2 mediates the proteasome-dependent degradation of Accelerated Cell Death 11 (ACD11) to promote pathogen defense. Here, we demonstrate roles for XBAT35.2 and ACD11 in abiotic stress tolerance. As seen in response to pathogen infection, abiotic stress stabilizes XBAT35.2 and the abundance of ACD11 rose consistently with increasing concentrations of abscisic acid (ABA) and salt. Surprisingly, exposure to ABA and salt increased the stability of ACD11, and the overexpression of ACD11 improves plant survival of salt and drought stress, suggesting a role for ACD11 in promoting tolerance. Prolonged exposure to high concentrations of ABA or salt resulted in ubiquitination and the proteasome-dependent degradation of ACD11, however. The stress-induced turnover of ACD11 requires XBAT35.2, as degradation is slowed in the absence of the E3 ubiquitin ligase. Consistent with XBAT35.2 mediating the proteasome-dependent degradation of ACD11, the loss of E3 ubiquitin ligase function enhances the tolerance of salt and drought stress, whereas overexpression increases sensitivity. A model is presented where, upon the perception of abiotic stress, ACD11 abundance increases to promote tolerance. Meanwhile, XBAT35.2 accumulates and in turn promotes the degradation of ACD11 to attenuate the stress response. The results characterize XBAT35.2 as an E3 ubiquitin ligase with opposing roles in abiotic and biotic stress.

Species and population specific gene expression in blood transcriptomes of marine turtles.

BACKGROUND: Transcriptomic data has demonstrated utility to advance the study of physiological diversity and organisms' responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations. RESULTS: We generated high quality blood transcriptome assemblies for hawksbill (Eretmochelys imbricata), loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance. CONCLUSIONS: Obtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.

Evolutionary History of GLIS Genes Illuminates Their Roles in Cell Reprograming and Ciliogenesis.

The GLIS family transcription factors, GLIS1 and GLIS3, potentiate generation of induced pluripotent stem cells (iPSCs). In contrast, another GLIS family member, GLIS2, suppresses cell reprograming. To understand how these disparate roles arose, we examined evolutionary origins and genomic organization of GLIS genes. Comprehensive phylogenetic analysis shows that GLIS1 and GLIS3 originated during vertebrate whole genome duplication, whereas GLIS2 is a sister group to the GLIS1/3 and GLI families. This result is consistent with their opposing functions in cell reprograming. Glis1 evolved faster than Glis3, losing many protein-interacting motifs. This suggests that Glis1 acquired new functions under weakened evolutionary constraints. In fact, GLIS1 induces induced pluripotent stem cells more strongly. Transcriptomic data from various animal embryos demonstrate that glis1 is maternally expressed in some tetrapods, whereas vertebrate glis3 and invertebrate glis1/3 genes are rarely expressed in oocytes, suggesting that vertebrate (or tetrapod) Glis1 acquired a new expression domain and function as a maternal factor. Furthermore, comparative genomic analysis reveals that glis1/3 is part of a bilaterian-specific gene cluster, together with rfx3, ndc1, hspb11, and lrrc42. Because known functions of these genes are related to cilia formation and function, the last common ancestor of bilaterians may have acquired this cluster by shuffling gene order to establish more sophisticated epithelial tissues involving cilia. This evolutionary study highlights the significance of GLIS1/3 for cell reprograming, development, and diseases in ciliated organs such as lung, kidney, and pancreas.

Unusual orthologs shed new light on the binding mechanism of ghrelin to its receptor GHSR1a.

The gastric peptide ghrelin has important functions in energy metabolism and cellular homeostasis by activating growth hormone secretagogue receptor type 1a (GHSR1a). The N-terminal residues of ghrelin orthologs from all vertebrates are quite conserved; however, in orthologs from Cavia porcellus and Phyllostomus discolor, Ser2 and Leu5 are replaced by a smaller Ala and a positively charged Arg, respectively. In the present study, we first demonstrated that the hydrophobic Leu5 is essential for the function of human ghrelin, because Ala replacement caused an approximately 100-fold decrease in activity. However, replacement of Leu5 by an Arg residue caused much less disruption; further replacement of Ser2 by Ala almost restored full activity, although the [S2A] mutation itself showed slight detriments, implying that the positively charged Arg5 in the [S2A,L5R] mutant might form alternative interactions with certain receptor residues to compensate for the loss of the essential Leu5. To identify the responsible receptor residues, we screened GHSR1a mutants in which all conserved negatively charged residues in the extracellular regions and all aromatic residues in the ligand-binding pocket were mutated separately. According to the decrease in selectivity of the mutant receptors towards [S2A,L5R]ghrelin, we deduced that the positively charged Arg5 of the ghrelin mutant primarily interacts with the essential aromatic Phe286 at the extracellular end of the sixth transmembrane domain of GHSR1a by forming cation-π and π-π interactions. The present study provided new insights into the binding mechanism of ghrelin with its receptor, and thus would facilitate the design of novel ligands for GHSR1a.

A family with partially penetrant multicentric carpotarsal osteolysis due to gonadal mosaicism: First reported case.

Multicentric carpotarsal osteolysis (MCTO) is an autosomal dominant condition characterized by carpal-tarsal abnormalities; over half of affected individuals also develop renal disease. MCTO is caused by mutations of MAFB; however, there is no clear phenotype-genotype correlation. We describe the first reported family of variable MCTO phenotype due to mosaicism: the proband had classical skeletal features and renal involvement due to focal segmental glomerulosclerosis (FSGS), and the father had profound renal impairment due to FSGS, necessitating kidney transplantation. Mosaicism was first suspected in this family due to unequal allele ratios in the sequencing chromatograph of the initial blood sample of proband's father and confirmed by sequencing DNA extracted from the father's hair, collected from different bodily parts. This case highlights the need for a high index of clinical suspicion to detect low-level parental mosaicism, as well as a potential role for MAFB mutation screening in individuals with isolated FSGS.

Pairwise comparisons across species are problematic when analyzing functional genomic data.

There is considerable interest in comparing functional genomic data across species. One goal of such work is to provide an integrated understanding of genome and phenotype evolution. Most comparative functional genomic studies have relied on multiple pairwise comparisons between species, an approach that does not incorporate information about the evolutionary relationships among species. The statistical problems that arise from not considering these relationships can lead pairwise approaches to the wrong conclusions and are a missed opportunity to learn about biology that can only be understood in an explicit phylogenetic context. Here, we examine two recently published studies that compare gene expression across species with pairwise methods, and find reason to question the original conclusions of both. One study interpreted pairwise comparisons of gene expression as support for the ortholog conjecture, the hypothesis that orthologs tend to have more similar attributes (expression in this case) than paralogs. The other study interpreted pairwise comparisons of embryonic gene expression across distantly related animals as evidence for a distinct evolutionary process that gave rise to phyla. In each study, distinct patterns of pairwise similarity among species were originally interpreted as evidence of particular evolutionary processes, but instead, we find that they reflect species relationships. These reanalyses concretely show the inadequacy of pairwise comparisons for analyzing functional genomic data across species. It will be critical to adopt phylogenetic comparative methods in future functional genomic work. Fortunately, phylogenetic comparative biology is also a rapidly advancing field with many methods that can be directly applied to functional genomic data.

Knee Osteoarthritis Progression Is Delayed in Silent Information Regulator 2 Ortholog 1 Knock-in Mice.

Overexpression of silent information regulator 2 ortholog 1 (SIRT1) is associated with beneficial roles in aging-related diseases; however, the effects of SIRT1 overexpression on osteoarthritis (OA) progression have not yet been studied. The aim of this study was to investigate OA progression in SIRT1-KI mice using a mouse OA model. OA was induced via destabilization of the medial meniscus using 12-week-old SIRT1-KI and wild type (control) mice. OA progression was evaluated histologically based on the Osteoarthritis Research Society International (OARSI) score at 4, 8, 12, and 16 weeks after surgery. The production of SIRT1, type II collagen, MMP-13, ADAMTS-5, cleaved caspase 3, Poly (ADP-ribose) polymerase (PARP) p85, acetylated NF-κB p65, interleukin 1 beta (IL-1ß), and IL-6 was examined via immunostaining. The OARSI scores were significantly lower in SIRT1-KI mice than those in control mice at 8, 12, and 16 weeks after surgery. The proportion of SIRT1 and type II collagen-positive-chondrocytes was significantly higher in SIRT1-KI mice than that in control mice. Moreover, the proportion of MMP-13-, ADAMTS-5-, cleaved caspase 3-, PARP p85-, acetylated NF-κB p65-, IL-1ß-, and IL-6-positive chondrocytes was significantly lower in SIRT1-KI mice than that in control mice. The mechanically induced OA progression was delayed in SIRT1-KI mice compared to that in control mice. Therefore, overexpression of SIRT1 may represent a mechanism for delaying OA progression.