UHRF1P contributes to IL-17A-mediated systemic lupus erythematosus via UHRF1-MAP4K3 axis.

Inflammatory T cells contribute to the pathogenesis of autoimmune diseases such as systemic lupus erythematosus (SLE). Analysis of the T-cell transcriptomics data of two independent SLE patient cohorts by three machine learning models revealed the pseudogene UHRF1P as a novel SLE biomarker. The pseudogene-encoded UHRF1P protein was overexpressed in peripheral blood T cells of SLE patients. The UHRF1P protein lacks the amino-terminus of its parental UHRF1 protein, resulting in missing the proteasome-binding ubiquitin-like (Ubl) domain of UHRF1. T-cell-specific UHRF1P transgenic mice manifested the induction of IL-17A and autoimmune inflammation. Mechanistically, UHFR1P prevented UHRF1-induced Lys48-linked ubiquitination and degradation of MAP4K3 (GLK), which is a kinase known to induce IL-17A. Consistently, IL-17A induction and autoimmune phenotypes of UHRF1P transgenic mice were obliterated by MAP4K3 knockout. Collectively, UHRF1P overexpression in T cells inhibits the E3 ligase function of its parental UHRF1 and induces autoimmune diseases.

Convergent evolution of kidney sizes and supraorbital salt glands for birds living in saline habitats.

Only a small number of avian species inhabit salty environments. To understand how they adapted, we examined the evolution of kidney sizes, supraorbital salt glands (SSGs), and the utilization of salty habitats across 230 species spanning 25 avian orders. Phylogenetic analysis indicates that SSGs, large kidneys, and thriving in salty habitats emerged convergently in birds. Transition rate analysis reveals that species possessing SSGs and large kidneys tended to move from low-to high-salinity environments, while others moved in the opposite direction. However, habitat salinity also influenced kidney evolution; lineages residing in high-salinity environments tended to develop larger kidneys than those in low-salinity environments. Our findings suggest that SSGs and large kidneys may have evolved through adaptation to high salinity. Overall, habitat conditions and physiological traits influenced avian adaptation to salty environments in a reciprocal manner. These results shed the new light on the evolutionary mechanisms underlying functional diversity in birds.

Degeneration of the Olfactory System in a Murid Rodent that Evolved Diurnalism.

In mammalian research, it has been debated what can initiate an evolutionary tradeoff between different senses, and the phenomenon of sensory tradeoff in rodents, the most abundant mammalian clade, is not evident. The Nile rat (Arvicanthis niloticus), a murid rodent, recently adapted to a diurnal niche through an evolutionary acquisition of daylight vision with enhanced visual acuity. As such, this model provides an opportunity for a cross-species investigation where comparative morphological and multi-omic analyses of the Nile rat are made with its closely related nocturnal species, e.g. the mouse (Mus musculus) and the rat (Rattus norvegicus). Thus, morphological examinations were performed, and evolutionary reductions in relative sizes of turbinal bone surfaces, the cribriform plate, and the olfactory bulb were discovered in Nile rats. Subsequently, we compared multiple murid genomes, and profiled olfactory epithelium transcriptomes of mice and Nile rats at various ages with RNA sequencing. The results further demonstrate that, in comparison with mouse olfactory receptor (OR) genes, Nile rat OR genes have experienced less frequent gain, more frequent loss, and more frequent expression reduction during their evolution. Furthermore, functional degeneration of coding sequences in the Nile rat lineage was found in OR genes, yet not in other genes. Taken together, these results suggest that acquisition of improved vision in the Nile rat has been accompanied by degeneration of both olfaction-related anatomical structures and OR gene repertoires, consistent with the hypothesis of an olfaction-vision tradeoff initiated by the switch from a nocturnal to a diurnal lifestyle in mammals.

Phenotypic heterogeneity in human genetic diseases: ultrasensitivity-mediated threshold effects as a unifying molecular mechanism.

Phenotypic heterogeneity is very common in genetic systems and in human diseases and has important consequences for disease diagnosis and treatment. In addition to the many genetic and non-genetic (e.g., epigenetic, environmental) factors reported to account for part of the heterogeneity, we stress the importance of stochastic fluctuation and regulatory network topology in contributing to phenotypic heterogeneity. We argue that a threshold effect is a unifying principle to explain the phenomenon; that ultrasensitivity is the molecular mechanism for this threshold effect; and discuss the three conditions for phenotypic heterogeneity to occur. We suggest that threshold effects occur not only at the cellular level, but also at the organ level. We stress the importance of context-dependence and its relationship to pleiotropy and edgetic mutations. Based on this model, we provide practical strategies to study human genetic diseases. By understanding the network mechanism for ultrasensitivity and identifying the critical factor, we may manipulate the weak spot to gently nudge the system from an ultrasensitive state to a stable non-disease state. Our analysis provides a new insight into the prevention and treatment of genetic diseases.

Long-read-based Genome Assembly of Drosophila gunungcola Reveals Fewer Chemosensory Genes in Flower-breeding Species.

Drosophila gunungcola exhibits reproductive activities on the fresh flowers of several plant species and is an emerging model to study the co-option of morphological and behavioral traits in male courtship display. Here, we report a near-chromosome-level genome assembly that was constructed based on long-read PacBio sequencing data (with ∼66× coverage) and annotated with the assistant from RNA-seq transcriptome data of whole organisms at various developmental stages. A nuclear genome of 189 Mb with 13,950 protein-coding genes and a mitogenome of 17.5 kb were acquired. Few interchromosomal rearrangements were found in the comparisons of synteny with Drosophila elegans, its sister species, and Drosophila melanogaster, suggesting that the gene compositions on each Muller element are evolutionarily conserved. Loss events of several OR and IR genes in D. gunungcola and D. elegans were revealed when orthologous genomic regions were compared across species in the D. melanogaster species group. This high-quality reference genome will facilitate further comparative studies on traits related to the evolution of sexual behavior and diet specialization.

Weak gene-gene interaction facilitates the evolution of gene expression plasticity.

BACKGROUND: Individual organisms may exhibit phenotypic plasticity when they acclimate to different conditions. Such plastic responses may facilitate or constrain the adaptation of their descendant populations to new environments, complicating their evolutionary trajectories beyond the genetic blueprint. Intriguingly, phenotypic plasticity itself can evolve in terms of its direction and magnitude during adaptation. However, we know little about what determines the evolution of phenotypic plasticity, including gene expression plasticity. Recent laboratory-based studies suggest dominance of reversing gene expression plasticity-plastic responses that move the levels of gene expression away from the new optima. Nevertheless, evidence from natural populations is still limited. RESULTS: Here, we studied gene expression plasticity and its evolution in the montane and lowland populations of an elevationally widespread songbird-the Rufous-capped Babbler (Cyanoderma ruficeps)-with reciprocal transplant experiments and transcriptomic analyses; we set common gardens at altitudes close to these populations' native ranges. We confirmed the prevalence of reversing plasticity in genes associated with altitudinal adaptation. Interestingly, we found a positive relationship between magnitude and degree of evolution in gene expression plasticity, which was pertinent to not only adaptation-associated genes but also the whole transcriptomes from multiple tissues. Furthermore, we revealed that genes with weaker expressional interactions with other genes tended to exhibit stronger plasticity and higher degree of plasticity evolution, which explains the positive magnitude-evolution relationship. CONCLUSIONS: Our experimental evidence demonstrates that species may initiate their adaptation to new habitats with genes exhibiting strong expression plasticity. We also highlight the role of expression interdependence among genes in regulating the magnitude and evolution of expression plasticity. This study illuminates how the evolution of phenotypic plasticity in gene expression facilitates the adaptation of species to challenging environments in nature.

Epigenomic signatures on paralogous genes reveal underappreciated universality of active histone codes adopted across animals.

The results of conventional gene-based analyses which combine epigenome and transcriptome data, including those conducted by the ENCODE/modENCODE projects, suggest various histone modifications performing regulatory functions in controlling mRNA expression (referred to as a histone code) in several model animals. While some histone codes were found to be universally adopted across organisms, "species-specific" histone codes have also been defined. We found that the characterization of these histone codes was confounded by factors (e.g. gene essentiality, expression breadth) that are independent of, but correlated with, gene expression levels. Hence, we attempted to decode histone marks in mouse (Mus musculus), fly (Drosophila melanogaster), and worm (Caenorhabditis elegans) genomes by examining ratios of RNA sequencing (and chromatin immunoprecipitation sequencing) intensities between paralog genes to remove confounding effects that would otherwise be present in a gene-based approach. With this paralog-based approach, associations between four histone modifications (H3K4me3, H3K27ac, H3K9ac, and H3K36me3) and gene expression are substantially revised. For example, we demonstrate that H3K27ac and H3K9ac represent universal active marks in promoters, rather than worm-specific marks as previously reported. Second, acting regions of the studied active marks that are common across species (and across a wide range of tissues at different developmental stages) were found to extend beyond the previously defined regions. Thus, it appears that the active histone codes analyzed have a universality that has previously been underappreciated. Our results suggested that these universal codes, including those previously considered species-specific, could have an ancient origin, and are important in regulating animal gene expression abundance.

Characterization of TMAO productivity from carnitine challenge facilitates personalized nutrition and microbiome signatures discovery.

The capability of gut microbiota in degrading foods and drugs administered orally can result in diversified efficacies and toxicity interpersonally and cause significant impact on human health. Production of atherogenic trimethylamine N-oxide (TMAO) from carnitine is a gut microbiota-directed pathway and varies widely among individuals. Here, we demonstrated a personalized TMAO formation and carnitine bioavailability from carnitine supplements by differentiating individual TMAO productivities with a recently developed oral carnitine challenge test (OCCT). By exploring gut microbiome in subjects characterized by TMAO producer phenotypes, we identified 39 operational taxonomy units that were highly correlated to TMAO productivity, including Emergencia timonensis, which has been recently discovered to convert γ-butyrobetaine to TMA in vitro. A microbiome-based random forest classifier was therefore constructed to predict the TMAO producer phenotype (AUROC = 0.81) which was then validated with an external cohort (AUROC = 0.80). A novel bacterium called Ihubacter massiliensis was also discovered to be a key microbe for TMA/TMAO production by using an OCCT-based humanized gnotobiotic mice model. Simply combining the presence of E. timonensis and I. massiliensis could account for 43% of high TMAO producers with 97% specificity. Collectively, this human gut microbiota phenotype-directed approach offers potential for developing precision medicine and provides insights into translational research. Video Abstract.

Comparative Transcriptomics across Nematode Life Cycles Reveal Gene Expression Conservation and Correlated Evolution in Adjacent Developmental Stages.

Nematodes are highly abundant animals with diverse habitats and lifestyles. Some are free living whereas others parasitize animals or plants, and among the latter, infection abilities change across developmental stages to infect hosts and complete life cycles. To determine the relationship between transcriptome evolution and morphological divergences among nematodes, we compared 48 transcriptomes of different developmental stages across eight nematode species. The transcriptomes were clustered broadly into embryo, larva, and adult stages, with the developmental plastic stages were separated from common larval stages within the larval branch. This suggests that development was the major determining factor after lifestyle changes, such as parasitism, during transcriptome evolution. Such patterns were partly accounted for by tissue-specific genes-such as those in oocytes and the hypodermis-being expressed at different proportions. Although nematodes typically have 3-5 larval stages, the transcriptomes for these stages were found to be highly correlated within each species, suggesting high similarity among larval stages across species. For the Caenorhabditis elegans-Caenorhabditis briggsae and Strongyloides stercoralis-Strongyloides venezuelensis comparisons, we found that ∼50% of genes were expressed at multiple stages, whereas half of their orthologs were also expressed in multiple but different stages. Such frequent changes in expression have resulted in concerted transcriptome evolution across adjacent stages, thus generating species-specific transcriptomes over the course of nematode evolution. Our study provides a first insight into the evolution of nematode transcriptomes beyond embryonic development.

Novel eye genes systematically discovered through an integrated analysis of mouse transcriptomes and phenome.

In the last few decades, reverse genetic and high throughput approaches have been frequently applied to the mouse (Mus musculus) to understand how genes function in tissues/organs and during development in a mammalian system. Despite these efforts, the associated phenotypes for the majority of mouse genes remained to be fully characterized. Here, we performed an integrated transcriptome-phenome analysis by identifying coexpressed gene modules based on tissue transcriptomes profiled with each of various platforms and functionally interpreting these modules using the mouse phenotypic data. Consequently, >15,000 mouse genes were linked with at least one of the 47 tissue functions that were examined. Specifically, our approach predicted >50 genes previously unknown to be involved in mice (Mus musculus) visual functions. Fifteen genes were selected for further analysis based on their potential biomedical relevance and compatibility with further experimental validation. Gene-specific morpholinos were introduced into zebrafish (Danio rerio) to target their corresponding orthologs. Quantitative assessments of phenotypes of developing eyes confirmed predicted eye-related functions of 13 out of the 15 genes examined. These novel eye genes include: Adal, Ankrd33, Car14, Ccdc126, Dhx32, Dkk3, Fam169a, Grifin, Kcnj14, Lrit2, Ppef2, Ppm1n, and Wdr17. The results highlighted the potential for this phenome-based approach to assist the experimental design of mutating and phenotyping mouse genes that aims to fully reveal the functional landscape of mammalian genomes.

Reduced Translational Efficiency of Eukaryotic Genes after Duplication Events.

Control of gene expression has been found to be predominantly determined at the level of protein translation. However, to date, reduced expression from duplicated genes in eukaryotes for dosage maintenance has only been linked to transcriptional control involving epigenetic mechanisms. Here, we hypothesize that dosage maintenance following gene duplication also involves regulation at the protein level. To test this hypothesis, we compared transcriptome and proteome data of yeast models, Saccharomyces cerevisiae and Schizosaccharomyces pombe, and worm models, Caenorhabditis elegans and Caenorhabditis briggsae, to investigate lineage-specifically duplicated genes. Duplicated genes in both eukaryotic models exhibited a reduced protein-to-mRNA abundance ratio. Moreover, dosage sensitive genes, represented by genes encoding protein complex subunits, reduced their protein-to-mRNA abundance ratios more significantly than the other genes after duplication events. An analysis of ribosome profiling (Ribo-Seq) data further showed that reduced translational efficiency was more prominent for dosage sensitive genes than for the other genes. Meanwhile, no difference in protein degradation rate was associated with duplication events. Translationally repressed duplicated genes were also more likely to be inhibited at the level of transcription. Taken together, these results suggest that translation-mediated dosage control is partially contributed by natural selection and it enhances transcriptional control in maintaining gene dosage after gene duplication events during eukaryotic genome evolution.

Genome-Wide Microarray Analysis Suggests Transcriptomic Response May Not Play a Major Role in High- to Low-Altitude Acclimation in Harvest Mouse (Micromys minutus).

The harvest mouse (Micromys minutus) is a small rodent species with a wide range of vertical distribution in Taiwan, extending from the sea level to 3100 m altitude. This species has recently suffered from habitat loss in high-altitude areas due to orchard cultivation, which may have resulted in mouse migration from high to low altitude. To investigate whether there is any physiological mechanism involved in altitude acclimation, rat cDNA microarray was used to compare transcriptomic patterns of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site. Of the 23,188 genes being analyzed, 47 (33 up-regulated and 14 down-regulated) were found to have differential expression (fold change > 4 or < -4, ANOVA p < 0.05). However, after multiple testing correction with a false discovery rate (FDR), only the result for Tnfrsf12a was found to be statistically significant (fold change = 13, FDR p < 0.05). The result was confirmed by quantitative polymerase chain reaction (q-PCR). The expression of Tnfrsf12a possibly relates to the skeletal muscle biology and thus can be correlated with altitude acclimation. However, finding only one gene transcript with significant alteration suggests that transcriptomic response may not play a major role in high- to low-altitude acclimation in harvest mouse.

Standing genetic variation as the predominant source for adaptation of a songbird.

What kind of genetic variation contributes the most to adaptation is a fundamental question in evolutionary biology. By resequencing genomes of 80 individuals, we inferred the origin of genomic variants associated with a complex adaptive syndrome involving multiple quantitative traits, namely, adaptation between high and low altitudes, in the vinous-throated parrotbill (Sinosuthora webbiana) in Taiwan. By comparing these variants with those in the Asian mainland population, we revealed standing variation in 24 noncoding genomic regions to be the predominant genetic source of adaptation. Parrotbills at both high and low altitudes exhibited signatures of recent selection, suggesting that not only the front but also the trailing edges of postglacial expanding populations could be subjected to environmental stresses. This study verifies and quantifies the importance of standing variation in adaptation in a cohort of genes, illustrating that the evolutionary potential of a population depends significantly on its preexisting genetic diversity. These findings provide important context for understanding adaptation and conservation of species in the Anthropocene.

Identification of TMAO-producer phenotype and host-diet-gut dysbiosis by carnitine challenge test in human and germ-free mice.

OBJECTIVE: The gut microbiota-derived metabolite, trimethylamine N-oxide (TMAO) plays an important role in cardiovascular disease (CVD). The fasting plasma TMAO was shown as a prognostic indicator of CVD incident in patients and raised the interest of intervention targeting gut microbiota. Here we develop a clinically applicable method called oral carnitine challenge test (OCCT) for TMAO-related therapeutic drug efforts assessment and personalising dietary guidance. DESIGN: A pharmacokinetic study was performed to verify the design of OCCT protocol. The OCCT was conducted in 23 vegetarians and 34 omnivores to validate gut microbiota TMAO production capacity. The OCCT survey was integrated with gut microbiome, host genotypes, dietary records and serum biochemistry. A humanised gnotobiotic mice study was performed for translational validation. RESULTS: The OCCT showed better efficacy than fasting plasma TMAO to identify TMAO producer phenotype. The omnivores exhibited a 10-fold higher OR to be high TMAO producer than vegetarians. The TMAO-associated taxa found by OCCT in this study were consistent with previous animal studies. The TMAO producer phenotypes were also reproduced in humanised gnotobiotic mice model. Besides, we found the faecal CntA gene was not associated with TMAO production; therefore, other key relevant microbial genes might be involved. Finally, we demonstrated the urine TMAO exhibited a strong positive correlation with plasma TMAO (r=0.92, p<0.0001) and improved the feasibility of OCCT. CONCLUSION: The OCCT can be used to identify TMAO-producer phenotype of gut microbiota and may serve as a personal guidance in CVD prevention and treatment. TRIAL REGISTRATION NUMBER: NCT02838732; Results.

modPhEA: model organism Phenotype Enrichment Analysis of eukaryotic gene sets.

MOTIVATION: Genome-scale phenotypic data are available for many model organisms, yet existing tools to functionally interpret gene sets from these phenotypic data are largely based on mutagenesis-derived phenotypes observed in mouse or human. RESULTS: Data from both mutagenesis and knockdown experiments are incorporated into modPhEA to allow users to perform enrichment analyses based on phenotypes observed in budding yeast (Saccharomyces cerevisiae), roundworm (Caenorhabditis elegans), fruit fly (Drosophila melanogaster), zebrafish (Danio rerio), mouse (Mus musculus) and humans (Homo sapiens). The phenotypes analysed can be customized to investigate complex traits and gene sets from any fully sequenced animal or fungal genome are also supported by modPhEA. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://evol.nhri.org.tw/modPhEA/. CONTACT: liaoby@nhri.org.tw. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Recruitment of histone modifications to assist mRNA dosage maintenance after degeneration of cytosine DNA methylation during animal evolution.

Following gene duplication, mRNA expression of the duplicated gene is reduced to maintain mRNA dosage. In mammals, this process is achieved with increased cytosine DNA methylation of the promoters of duplicated genes to suppress transcriptional initiation. However, not all animal species possess a full apparatus for cytosine DNA methylation. For such species, such as the roundworm (Caenorhabditis elegans, "worm" hereafter) or fruit fly (Drosophila melanogaster, "fly" hereafter), it is unclear how reduced expression of duplicated genes has been achieved evolutionarily. Here, we hypothesize that in the absence of a classical cytosine DNA methylation pathway, histone modifications play an increasing role in maintaining mRNA dosage following gene duplication. We initially verified that reduced gene expression of duplicated genes had occurred in the worm, fly, and mouse (Mus musculus). Next, several histone marks, with the capacity to control mRNA abundance in the models studied, were examined. In the worm and fly, but not in the mouse, multiple histone modifications were found to assist mRNA dosage maintenance following gene duplication events and the possible involvement of adenine DNA methylation in this process was excluded. Furthermore, the histone marks and acting regions that mediated the reduction in duplicated gene expression were found to be largely organism specific. Thus, it appears that many of the histone marks that maintain mRNA dosage were independently recruited during the evolution of worms and flies to compensate for the loss of cytosine DNA methylation machinery from their genomes.

Functionalities of expressed messenger RNAs revealed from mutant phenotypes.

Total messenger RNAs mRNAs that are produced from a given gene under a certain set of conditions include both functional and nonfunctional transcripts. The high prevalence of nonfunctional mRNAs that have been detected in cells has raised questions regarding the functional implications of mRNA expression patterns and divergences. Phenotypes that result from the mutagenesis of protein-coding genes have provided the most straightforward descriptions of gene functions, and such data obtained from model organisms have facilitated investigations of the functionalities of expressed mRNAs. Mutant phenotype data from mouse tissues have revealed various attributes of functional mRNAs, including tissue-specificity, strength of expression, and evolutionary conservation. In addition, the role that mRNA expression evolution plays in driving morphological evolution has been revealed from studies designed to exploit morphological and physiological phenotypes of mouse mutants. Investigations into yeast essential genes (defined by an absence of colony growth after gene deletion) have further described gene regulatory strategies that reduce protein expression noise by mediating the rates of transcription and translation. In addition to the functional significance of expressed mRNAs as described in the abovementioned findings, the functionalities of other type of RNAs (i.e., noncoding RNAs) remain to be characterized with systematic mutations and phenotyping of the DNA regions that encode these RNA molecules. WIREs RNA 2016, 7:416-427. doi: 10.1002/wrna.1329 For further resources related to this article, please visit the WIREs website.

Proteins with Highly Evolvable Domain Architectures Are Nonessential but Highly Retained.

The functions of proteins are usually determined by domains, and the sequential order in which domains are connected to make up a protein chain is known as the domain architecture. Here, we constructed evolutionary networks of protein domain architectures in species from three major life lineages (bacteria, fungi, and metazoans) by connecting any two architectures between which an evolutionary event could be inferred by a model that assumes maximum parsimony. We found that proteins with domain architectures with a higher level of evolvability, indicated by a greater number of connections in the evolutionary network, are present in a wider range of species. However, these proteins tend to be less essential to the organism, are duplicated more often during evolution, have more isoforms, and, intriguingly, tend to be associated with functional categories important for organismal adaptation. These results reveal the presence, in many genomes, of genes coding for a core set of nonessential proteins that have a highly evolvable domain architecture and thus a repertoire of genetic materials accessible for organismal adaptation.

Metabolic characteristics of dominant microbes and key rare species from an acidic hot spring in Taiwan revealed by metagenomics.

BACKGROUND: Microbial diversity and community structures in acidic hot springs have been characterized by 16S rRNA gene-based diversity surveys. However, our understanding regarding the interactions among microbes, or between microbes and environmental factors, remains limited. RESULTS: In the present study, a metagenomic approach, followed by bioinformatics analyses, were used to predict interactions within the microbial ecosystem in Shi-Huang-Ping (SHP), an acidic hot spring in northern Taiwan. Characterizing environmental parameters and potential metabolic pathways highlighted the importance of carbon assimilatory pathways. Four distinct carbon assimilatory pathways were identified in five dominant genera of bacteria. Of those dominant carbon fixers, Hydrogenobaculum bacteria outcompeted other carbon assimilators and dominated the SHP, presumably due to their ability to metabolize hydrogen and to withstand an anaerobic environment with fluctuating temperatures. Furthermore, most dominant microbes were capable of metabolizing inorganic sulfur-related compounds (abundant in SHP). However, Acidithiobacillus ferrooxidans was the only species among key rare microbes with the capability to fix nitrogen, suggesting a key role in nitrogen cycling. In addition to potential metabolic interactions, based on the 16S rRNAs gene sequence of Nanoarchaeum-related and its potential host Ignicoccus-related archaea, as well as sequences of viruses and CRISPR arrays, we inferred that there were complex microbe-microbe interactions. CONCLUSIONS: Our study provided evidence that there were numerous microbe-microbe and microbe-environment interactions within the microbial community in an acidic hot spring. We proposed that Hydrogenobaculum bacteria were the dominant microbial genus, as they were able to metabolize hydrogen, assimilate carbon and live in an anaerobic environment with fluctuating temperatures.