Controlled Noise: Evidence of epigenetic regulation of Single-Cell expression variability.

MOTIVATION: Understanding single-cell expression variability (scEV) or gene expression noise among cells of the same type and state is crucial for delineating population-level cellular function. While epigenetic mechanisms are widely implicated in gene expression regulation, a definitive link between chromatin accessibility and scEV remains elusive. Recent advances in single-cell techniques enable the study of single-cell multiomics data that include the simultaneous measurement of scATAC-seq and scRNA-seq within individual cells, presenting an unprecedented opportunity to address this gap. RESULTS: This paper introduces an innovative testing pipeline to investigate the association between chromatin accessibility and scEV. With single-cell multiomics data of scATAC-seq and scRNA-seq, the pipeline hinges on comparing the prediction performance of scATAC-seq data on gene expression levels between highly variable genes (HVGs) and non-highly variable genes (non-HVGs). Applying this pipeline to paired scATAC-seq and scRNA-seq data from human hematopoietic stem and progenitor cells, we observed a significantly superior prediction performance of scATAC-seq data for HVGs compared to non-HVGs. Notably, there was substantial overlap between well-predicted genes and HVGs. The gene pathways enriched from well-predicted genes are highly pertinent to cell type-specific functions. Our findings support the notion that scEV largely stems from cell-to-cell variability in chromatin accessibility, providing compelling evidence for the epigenetic regulation of scEV and offering promising avenues for investigating gene regulation mechanisms at the single-cell level. AVAILABILITY: The source code and data used in this paper can be found at https://github.com/SiweiCui/EpigeneticControlOfSingle-CellExpressionVariability. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of a copper-responsive small molecule inhibitor of uropathogenic Escherichia coli.

Urinary tract infections (UTIs) are a major global health problem and are caused predominantly by uropathogenic Escherichia coli (UPEC). UTIs are a leading cause of prescription antimicrobial use. Incessant increase in antimicrobial resistance in UPEC and other uropathogens poses a serious threat to the current treatment practices. Copper is an effector of nutritional immunity that impedes the growth of pathogens during infection. We hypothesized that copper would augment the toxicity of select small molecules against bacterial pathogens. We conducted a small molecule screening campaign with a library of 51,098 molecules to detect hits that inhibit a UPEC ΔtolC mutant in a copper-dependent manner. A molecule, denoted as E. coli inhibitor or ECIN, was identified as a copper-responsive inhibitor of wild-type UPEC strains. Our gene expression and metal content analysis results demonstrate that ECIN works in concert with copper to exacerbate Cu toxicity in UPEC. ECIN has a broad spectrum of activity against pathogens of medical and veterinary significance including Acinetobacter baumannii, Pseudomonas aeruginosa, and methicillin-resistant Staphylococcus aureus. Subinhibitory levels of ECIN eliminate UPEC biofilm formation. Transcriptome analysis of UPEC treated with ECIN reveals induction of multiple stress response systems. Furthermore, we demonstrate that L-cysteine rescues the growth of UPEC exposed to ECIN. In summary, we report the identification and characterization of a novel copper-responsive small molecule inhibitor of UPEC.IMPORTANCEUrinary tract infection (UTI) is a ubiquitous infectious condition affecting millions of people annually. Uropathogenic Escherichia coli (UPEC) is the predominant etiological agent of UTI. However, UTIs are becoming increasingly difficult to resolve with antimicrobials due to increased antimicrobial resistance in UPEC and other uropathogens. Here, we report the identification and characterization of a novel copper-responsive small molecule inhibitor of UPEC. In addition to E. coli, this small molecule also inhibits pathogens of medical and veterinary significance including Acinetobacter baumannii, Pseudomonas aeruginosa, and methicillin-resistant Staphylococcus aureus.

Suppression of FOXO1 attenuates inflamm-aging and improves liver function during aging.

The liver is a key metabolic organ that maintains whole-body nutrient homeostasis. Aging-induced liver function alterations contribute to systemic susceptibility to aging-related diseases. However, the molecular mechanisms of liver aging remain insufficiently understood. In this study, we performed bulk RNA-Seq and single-cell RNA-Seq analyses to investigate the underlying mechanisms of the aging-induced liver function changes. We found that liver inflammation, glucose intolerance, and liver fat deposition were aggravated in old mice. Aging significantly increased pro-inflammation in hepatic macrophages. Furthermore, we found that Kupffer cells (KCs) were the major driver to induce pro-inflammation in hepatic macrophages during aging. In KCs, aging significantly increased pro-inflammatory levels; in monocyte-derived macrophages (MDMs), aging had a limited effect on pro-inflammation but led to a functional quiescence in antigen presentation and phagosome process. In addition, we identified an aging-responsive KC-specific (ARKC) gene set that potentially mediates aging-induced pro-inflammation in KCs. Interestingly, FOXO1 activity was significantly increased in the liver of old mice. FOXO1 inhibition by AS1842856 significantly alleviated glucose intolerance, hepatic steatosis, and systemic inflammation in old mice. FOXO1 inhibition significantly attenuated aging-induced pro-inflammation in KCs partially through downregulation of ARKC genes. However, FOXO1 inhibition had a limited effect on aging-induced functional quiescence in MDMs. These results indicate that aging induces pro-inflammation in liver mainly through targeting KCs and FOXO1 is a key player in aging-induced pro-inflammation in KCs. Thus, FOXO1 could be a potential therapeutic target for the treatment of age-associated chronic diseases.

Hepatocyte Adenosine Kinase Promotes Excessive Fat Deposition and Liver Inflammation.

BACKGROUND & AIMS: Nonalcoholic fatty liver disease is highly associated with obesity and progresses to nonalcoholic steatohepatitis when the liver develops overt inflammatory damage. While removing adenosine in the purine salvage pathway, adenosine kinase (ADK) regulates methylation reactions. We aimed to study whether hepatocyte ADK functions as an obesogenic gene/enzyme to promote excessive fat deposition and liver inflammation. METHODS: Liver sections of human subjects were examined for ADK expression using immunohistochemistry. Mice with hepatocyte-specific ADK disruption or overexpression were examined for hepatic fat deposition and inflammation. Liver lipidomics, hepatocyte RNA sequencing (RNA-seq), and single-cell RNA-seq for liver nonparenchymal cells were performed to analyze ADK regulation of hepatocyte metabolic responses and hepatocyte-nonparenchymal cells crosstalk. RESULTS: Whereas patients with nonalcoholic fatty liver disease had increased hepatic ADK levels, mice with hepatocyte-specific ADK disruption displayed decreased hepatic fat deposition on a chow diet and were protected from diet-induced excessive hepatic fat deposition and inflammation. In contrast, mice with hepatocyte-specific ADK overexpression displayed increased body weight and adiposity and elevated degrees of hepatic steatosis and inflammation compared with control mice. RNA-seq and epigenetic analyses indicated that ADK increased hepatic DNA methylation and decreased hepatic Ppara expression and fatty acid oxidation. Lipidomic and single-cell RNA-seq analyses indicated that ADK-driven hepatocyte factors, due to mitochondrial dysfunction, enhanced macrophage proinflammatory activation in manners involving increased expression of stimulator of interferon genes. CONCLUSIONS: Hepatocyte ADK functions to promote excessive fat deposition and liver inflammation through suppressing hepatocyte fatty acid oxidation and producing hepatocyte-derived proinflammatory mediators. Therefore, hepatocyte ADK is a therapeutic target for managing obesity and nonalcoholic fatty liver disease.

Aryl Hydrocarbon Receptor (AhR) Signaling in Colonic Cells and Tumors.

The aryl hydrocarbon receptor (AhR) is overexpressed in many tumor types and exhibits tumor-specific tumor promoter and tumor suppressor-like activity. In colon cancer, most but not all studies suggest that the AhR exhibits tumor suppressor activity which is enhanced by AhR ligands acting as agonists. Our studies investigated the role of the AhR in colon tumorigenesis using wild-type and AhR-knockout mice, the inflammation model of colon tumorigenesis using mice treated with azoxymethane (AOM)/dextran sodium sulfate (DSS) and APCS580/+; KrasG12D/+ mice all of which form intestinal tumors. The effects of tissue-specific AhR loss in the intestine of the tumor-forming mice on colonic stem cells, organoid-initiating capacity, colon tumor formation and mechanisms of AhR-mediated effects were investigated. Loss of AhR enhanced stem cell and tumor growth and in the AOM/DSS model AhR-dependent suppression of FOXM1 and downstream genes was important for AhR-dependent anticancer activity. Furthermore, the effectiveness of interleukin-22 (IL22) in colonic epithelial cells was also dependent on AhR expression. IL22 induced phosphorylation of STAT3, inhibited colonic organoid growth, promoted colonic cell proliferation in vivo and enhanced DNA repair in AOM/DSS-induced tumors. In this mouse model, the AhR suppressed SOCS3 expression and enhanced IL22-mediated activation of STAT3, whereas the loss of the AhR increased levels of SOCS3 which in turn inhibited IL22-induced STAT3 activation. In the APCS580/+; KrasG12D/+ mouse model, the loss of the AhR enhanced Wnt signaling and colon carcinogenesis. Results in both mouse models of colon carcinogenesis were complemented by single cell transcriptomics on colonic intestinal crypts which also showed that AhR deletion promoted expression of FOXM1-regulated genes in multiple colonic cell subtypes. These results support the role of the AhR as a tumor suppressor-like gene in the colon.

Ablation of long noncoding RNA MALAT1 activates antioxidant pathway and alleviates sepsis in mice.

The metastasis-associated lung adenocarcinoma transcript1 (MALAT1) is a long noncoding RNA (lncRNA) and is known for its role in cancer development and prognosis. In this study, we report that MALAT1 plays an important role in regulating acute inflammatory responses in sepsis. In patient samples, MALAT1 expression was positively correlated with severity of sepsis. In cultured macrophages, LPS treatment significantly induced MALAT1 expression, while genetic ablation of MALAT1 greatly reduced proinflammatory cytokine levels. Furthermore, MALAT1-ablated mice had significantly increased survival rates in cecal ligation and puncture (CLP)-induced sepsis and LPS-induced endotoxemia. One novel and salient feature of MALAT1-ablated mice is greatly reduced ROS level in macrophages and other cell types and increased glutathione/oxidized glutathione (GSH/GSSG) ratio in macrophages, suggesting an increased antioxidant capacity. We showed a mechanism for MALAT1 ablation leading to enhanced antioxidant capacity is through activation of methionine cycle by epitranscriptomical regulation of methionine adenosyltransferase 2A (MAT2A). MAT2A 3'UTR can be methylated by METTL16 which was known to directly bind to MALAT1. MALAT1 ablation was found to reduce methylation in MAT2A hairpin1 and increase MAT2A protein levels. Our results suggest a MALAT1-METTL16-MAT2A interactive axis which may be targeted for treatments of sepsis.

Transcriptomic Profiling of Gene Expression Associated with Granulosa Cell Tumor Development in a Mouse Model.

Ovarian granulosa cell tumors (GCTs) are rare sex cord-stromal tumors, accounting for ~5% ovarian tumors. The etiology of GCTs remains poorly defined. Genetically engineered mouse models are potentially valuable for understanding the pathogenesis of GCTs. Mice harboring constitutively active TGFß signaling (TGFBR1-CA) develop ovarian GCTs that phenocopy several hormonal and molecular characteristics of human GCTs. To determine molecular alterations in the ovary upon TGFß signaling activation, we performed transcriptomic profiling of gene expression associated with GCT development using ovaries from 1-month-old TGFBR1-CA mice and age-matched controls. RNA-sequencing and bioinformatics analysis coupled with the validation of select target genes revealed dysregulations of multiple cellular events and signaling molecules/pathways. The differentially expressed genes are enriched not only for known GCT-related pathways and tumorigenic events but also for signaling events potentially mediated by neuroactive ligand-receptor interaction, relaxin signaling, insulin signaling, and complements in TGFBR1-CA ovaries. Additionally, a comparative analysis of our data in mice with genes dysregulated in human GCTs or granulosa cells overexpressing a mutant FOXL2, the genetic hallmark of adult GCTs, identified some common genes altered in both conditions. In summary, this study has revealed the molecular signature of ovarian GCTs in a mouse model that harbors the constitutive activation of TGFBR1. The findings may be further exploited to understand the pathogenesis of a class of poorly defined ovarian tumors.

Single-cell RNA Sequencing Reveals How the Aryl Hydrocarbon Receptor Shapes Cellular Differentiation Potency in the Mouse Colon.

Despite recent progress recognizing the importance of aryl hydrocarbon receptor (Ahr)-dependent signaling in suppressing colon tumorigenesis, its role in regulating colonic crypt homeostasis remains unclear. To assess the effects of Ahr on intestinal epithelial cell heterogeneity and functional phenotypes, we utilized single-cell transcriptomics and advanced analytic strategies to generate a high-quality atlas for colonic intestinal crypts from wild-type and intestinal-specific Ahr knockout mice. Here we observed the promotive effects of Ahr deletion on Foxm1-regulated genes in crypt-associated canonical epithelial cell types and subtypes of goblet cells and deep crypt-secretory cells. We also show that intestinal Ahr deletion elevated single-cell entropy (a measure of differentiation potency or cell stemness) and RNA velocity length (a measure of the rate of cell differentiation) in noncycling and cycling Lgr5+ stem cells. In general, intercellular signaling cross-talk via soluble and membrane-bound factors was perturbed in Ahr-null colonocytes. Taken together, our single-cell RNA sequencing analyses provide new evidence of the molecular function of Ahr in modulating putative stem cell driver genes, cell potency lineage decisions, and cell-cell communication in vivo. PREVENTION RELEVANCE: Our mouse single-cell RNA sequencing analyses provide new evidence of the molecular function of Ahr in modulating colonic stemness and cell-cell communication in vivo. From a cancer prevention perspective, Ahr should be considered a therapeutic target to recalibrate remodeling of the intestinal stem cell niche.

Adipose tissue inflammation and systemic insulin resistance in mice with diet-induced obesity is possibly associated with disruption of PFKFB3 in hematopoietic cells.

Obesity-associated inflammation in white adipose tissue (WAT) is a causal factor of systemic insulin resistance; however, precisely how immune cells regulate WAT inflammation in relation to systemic insulin resistance remains to be elucidated. The present study examined a role for 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3 (PFKFB3) in hematopoietic cells in regulating WAT inflammation and systemic insulin sensitivity. Male C57BL/6J mice were fed a high-fat diet (HFD) or low-fat diet (LFD) for 12 weeks and examined for WAT inducible 6-phosphofructo-2-kinase (iPFK2) content, while additional HFD-fed mice were treated with rosiglitazone and examined for PFKFB3 mRNAs in WAT stromal vascular cells (SVC). Also, chimeric mice in which PFKFB3 was disrupted only in hematopoietic cells and control chimeric mice were also fed an HFD and examined for HFD-induced WAT inflammation and systemic insulin resistance. In vitro, adipocytes were co-cultured with bone marrow-derived macrophages and examined for adipocyte proinflammatory responses and insulin signaling. Compared with their respective levels in controls, WAT iPFK2 amount in HFD-fed mice and WAT SVC PFKFB3 mRNAs in rosiglitazone-treated mice were significantly increased. When the inflammatory responses were analyzed, peritoneal macrophages from PFKFB3-disrputed mice revealed increased proinflammatory activation and decreased anti-inflammatory activation compared with control macrophages. At the whole animal level, hematopoietic cell-specific PFKFB3 disruption enhanced the effects of HFD feeding on promoting WAT inflammation, impairing WAT insulin signaling, and increasing systemic insulin resistance. In vitro, adipocytes co-cultured with PFKFB3-disrupted macrophages revealed increased proinflammatory responses and decreased insulin signaling compared with adipocytes co-cultured with control macrophages. These results suggest that PFKFB3 disruption in hematopoietic cells only exacerbates HFD-induced WAT inflammation and systemic insulin resistance.

Single-cell RNA sequencing of a European and an African lymphoblastoid cell line.

In biomedical research, lymphoblastoid cell lines (LCLs), often established by in vitro infection of resting B cells with Epstein-Barr virus, are commonly used as surrogates for peripheral blood lymphocytes. Genomic and transcriptomic information on LCLs has been used to study the impact of genetic variation on gene expression in humans. Here we present single-cell RNA sequencing (scRNA-seq) data on GM12878 and GM18502-two LCLs derived from the blood of female donors of European and African ancestry, respectively. Cells from three samples (the two LCLs and a 1:1 mixture of the two) were prepared separately using a 10x Genomics Chromium Controller and deeply sequenced. The final dataset contained 7,045 cells from GM12878, 5,189 from GM18502, and 5,820 from the mixture, offering valuable information on single-cell gene expression in highly homogenous cell populations. This dataset is a suitable reference for population differentiation in gene expression at the single-cell level. Data from the mixture provide additional valuable information facilitating the development of statistical methods for data normalization and batch effect correction.

Long noncoding RNA MALAT1 regulates generation of reactive oxygen species and the insulin responses in male mice.

The metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) is a long noncoding RNA and its overexpression is associated with the development of many types of malignancy. MALAT1 null mice show no overt phenotype. However, in transcriptome analysis of MALAT1 null mice we found significant upregulation of nuclear factor-erythroid 2 p45-related factor 2 (Nrf2) regulated antioxidant genes including Nqo1 and Cat with significant reduction in reactive oxygen species (ROS) and greatly reduced ROS-generated protein carbonylation in hepatocyte and islets. We performed lncRNA pulldown assay using biotinylated antisense oligonucleotides against MALAT1 and found MALAT1 interacted with Nrf2, suggesting Nrf2 is transcriptionally regulated by MALAT1. Exposure to excessive ROS has been shown to cause insulin resistance through activation of c-Jun N-terminal kinase (JNK) which leads to inhibition of insulin receptor substrate 1 (IRS-1) and insulin-induced phosphorylation of serine/threonine kinase Akt. We found MALAT1 ablation suppressed JNK activity with concomitant insulin-induced activation of IRS-1 and phosphorylation of Akt suggesting MALAT1 regulated insulin responses. MALAT1 null mice exhibited sensitized insulin-signaling response to fast-refeeding and glucose/insulin challenges and significantly increased insulin secretion in response to glucose challenge in isolated MALAT1 null islets, suggesting an increased insulin sensitivity. In summary, we demonstrate that MALAT1 plays an important role in regulating insulin sensitivity and has the potential as a therapeutic target for the treatment of diabetes as well as other diseases caused by excessive exposure to ROS.

Integrating multiple fitting regression and Bayes decision for cancer diagnosis with transcriptomic data from tumor-educated blood platelets.

The application of machine learning in cancer diagnostics has shown great promise and is of importance in clinic settings. Here we consider applying machine learning methods to transcriptomic data derived from tumor-educated platelets (TEPs) from individuals with different types of cancer. We aim to define a reliability measure for diagnostic purposes to increase the potential for facilitating personalized treatments. To this end, we present a novel classification method called MFRB (for Multiple Fitting Regression and Bayes decision), which integrates the process of multiple fitting regression (MFR) with Bayes decision theory. MFR is first used to map multidimensional features of the transcriptomic data into a one-dimensional feature. The probability density function of each class in the mapped space is then adjusted using the Gaussian probability density function. Finally, the Bayes decision theory is used to build a probabilistic classifier with the estimated probability density functions. The output of MFRB can be used to determine which class a sample belongs to, as well as to assign a reliability measure for a given class. The classical support vector machine (SVM) and probabilistic SVM (PSVM) are used to evaluate the performance of the proposed method with simulated and real TEP datasets. Our results indicate that the proposed MFRB method achieves the best performance compared to SVM and PSVM, mainly due to its strong generalization ability for limited, imbalanced, and noisy data.

Adaptive human CDKAL1 variants underlie hormonal response variations at the enteroinsular axis.

Recent analyses have identified positively selected loci that explain differences in immune responses, body forms, and adaptations to extreme climates, but variants that describe adaptations in energy-balance regulation remain underexplored. To identify variants that confer adaptations in energy-balance regulation, we explored the evolutionary history and functional associations of candidate variants in 207 genes. We screened single nucleotide polymorphisms in genes that had been associated with energy-balance regulation for unusual genetic patterns in human populations, followed by studying associations among selected variants and serum levels of GIP, insulin, and C-peptide in pregnant women after an oral glucose tolerance test. Our analysis indicated that 5' variants in CDKAL1, CYB5R4, GAD2, and PPARG are marked with statistically significant signals of gene-environment interactions. Importantly, studies of serum hormone levels showed that variants in CDKAL1 are associated with glucose-induced GIP and insulin responses (p<0.05). On the other hand, a GAD2 variant exhibited a significant association with glucose-induced C-peptide response. In addition, simulation analysis indicated that a type 2 diabetes risk variant in CDKAL1 (rs7754840) was selected in East Asians ∼6,900 years ago. Taken together, these data indicated that variants in CDKAL1 and GAD2 were targets of prior environmental selection. Because the selection of the CDKAL1 variant overlapped with the selection of a cluster of GIP variants in the same population ∼11,800 to 2,000 years ago, we speculate that these regulatory genes at the human enteroinsular axis could be highly responsive to environmental selection in recent human history.

Population-level expression variability of mitochondrial DNA-encoded genes in humans.

Human mitochondria contain multiple copies of a circular genome made up of double-stranded DNA (mtDNA) that encodes proteins involved in cellular respiration. Transcript abundance of mtDNA-encoded genes varies between human individuals, yet the level of variation in the general population has not been systematically assessed. In the present study, we revisited large-scale RNA sequencing data generated from lymphoblastoid cell lines of HapMap samples of European and African ancestry to estimate transcript abundance and quantify expression variation for mtDNA-encoded genes. In both populations, we detected up to over 100-fold difference in mtDNA gene expression between individuals. The marked variation was not due to differences in mtDNA copy number between individuals, but was shaped by the transcription of hundreds of nuclear genes. Many of these nuclear genes were co-expressed with one another, resulting in a module-enriched co-expression network. Significant correlations in expression between genes of the mtDNA and nuclear genomes were used to identify factors involved with the regulation of mitochondrial functions. In conclusion, we determined the baseline amount of variability in mtDNA gene expression in general human populations and cataloged a complete set of nuclear genes whose expression levels are correlated with those of mtDNA-encoded genes. Our findings will enable the integration of information from both mtDNA and nuclear genetic systems, and facilitate the discovery of novel regulatory pathways involving mitochondrial functions.

Widespread divergence of the CEACAM/PSG genes in vertebrates and humans suggests sensitivity to selection.

In mammals, carcinoembryonic antigen cell adhesion molecules (CEACAMs) and pregnancy-specific glycoproteins (PSGs) play important roles in the regulation of pathogen transmission, tumorigenesis, insulin signaling turnover, and fetal-maternal interactions. However, how these genes evolved and to what extent they diverged in humans remain to be investigated specifically. Based on syntenic mapping of chordate genomes, we reveal that diverging homologs with a prototypic CEACAM architecture-including an extracellular domain with immunoglobulin variable and constant domain-like regions, and an intracellular domain containing ITAM motif-are present from cartilaginous fish to humans, but are absent in sea lamprey, cephalochordate or urochordate. Interestingly, the CEACAM/PSG gene inventory underwent radical divergence in various vertebrate lineages: from zero in avian species to dozens in therian mammals. In addition, analyses of genetic variations in human populations showed the presence of various types of copy number variations (CNVs) at the CEACAM/PSG locus. These copy number polymorphisms have 3-80% frequency in select populations, and encompass single to more than six PSG genes. Furthermore, we found that CEACAM/PSG genes contain a significantly higher density of nonsynonymous single nucleotide polymorphism (SNP) compared to the chromosome average, and many CEACAM/PSG SNPs exhibit high population differentiation. Taken together, our study suggested that CEACAM/PSG genes have had a more dynamic evolutionary history in vertebrates than previously thought. Given that CEACAM/PSGs play important roles in maternal-fetal interaction and pathogen recognition, these data have laid the groundwork for future analysis of adaptive CEACAM/PSG genotype-phenotypic relationships in normal and complicated pregnancies as well as other etiologies.

Identification of metabolic modifiers that underlie phenotypic variations in energy-balance regulation.

OBJECTIVE: Although recent studies have shown that human genomes contain hundreds of loci that exhibit signatures of positive selection, variants that are associated with adaptation in energy-balance regulation remain elusive. We reasoned that the difficulty in identifying such variants could be due to heterogeneity in selection pressure and that an integrative approach that incorporated experiment-based evidence and population genetics-based statistical judgments would be needed to reveal important metabolic modifiers in humans. RESEARCH DESIGN AND METHODS: To identify common metabolic modifiers that underlie phenotypic variation in diabetes-associated or obesity-associated traits in humans, or both, we screened 207 candidate loci for regulatory single nucleotide polymorphisms (SNPs) that exhibited evidence of gene-environmental interactions. RESULTS: Three SNPs (rs3895874, rs3848460, and rs937301) at the 5' gene region of human GIP were identified as prime metabolic-modifier candidates at the enteroinsular axis. Functional studies have shown that GIP promoter reporters carrying derived alleles of these three SNPs (haplotype GIP(-1920A)) have significantly lower transcriptional activities than those with ancestral alleles at corresponding positions (haplotype GIP(-1920G)). Consistently, studies of pregnant women who have undergone a screening test for gestational diabetes have shown that patients with a homozygous GIP(-1920A/A) genotype have significantly lower serum concentrations of glucose-dependent insulinotropic polypeptide (GIP) than those carrying an ancestral GIP(-1920G) haplotype. After controlling for a GIPR variation, we showed that serum glucose concentrations of patients carrying GIP(-1920A/A) homozygotes are significantly higher than that of those carrying an ancestral GIP(-1920G) haplotype (odds ratio 3.53). CONCLUSIONS: Our proof-of-concept study indicates that common regulatory GIP variants impart a difference in GIP and glucose metabolism. The study also provides a rare example that identified the common variant-common phenotypic variation pattern based on evidence of moderate gene-environmental interactions.

Context-dependent robustness to 5' splice site polymorphisms in human populations.

There has been growing evidence for extensive diversity of alternative splicing in human populations. Genetic variants within the 5' splice site can cause splicing differences among human individuals and constitute an important class of human disease mutations. In this study, we explored whether natural variations of splicing could reveal important signals of 5' splice site recognition. In seven lymphoblastoid cell lines of Asian, European and African ancestry, we identified 1174 single nucleotide polymorphisms (SNPs) within the consensus 5' splice site. We selected 129 SNPs predicted to significantly alter the splice site activity, and quantitatively examined their splicing impact in the seven individuals. Surprisingly, outside of the essential GT dinucleotide position, only ∼14% of the tested SNPs altered splicing. Bioinformatic and minigene analyses identified signals that could modify the impact of 5' splice site polymorphisms, most notably a strong 3' splice site and the presence of intronic motifs downstream of the 5' splice site. Strikingly, we found that the poly-G run, a known intronic splicing enhancer, was the most significantly enriched motif downstream of exons unaffected by 5' splice site SNPs. In TRIM62, the upstream 3' splice site and downstream intronic poly-G runs functioned redundantly to protect an exon from its 5' splice site polymorphism. Collectively, our study reveals widespread context-dependent robustness to 5' splice site polymorphisms in human transcriptomes. Consequently, certain exons are more susceptible to 5' splice site mutations. Additionally, our work demonstrates that genetic diversity of alternative splicing can provide significant insights into the splicing code of mammalian cells.

Adaptive selection of an incretin gene in Eurasian populations.

Diversities in human physiology have been partially shaped by adaptation to natural environments and changing cultures. Recent genomic analyses have revealed single nucleotide polymorphisms (SNPs) that are associated with adaptations in immune responses, obvious changes in human body forms, or adaptations to extreme climates in select human populations. Here, we report that the human GIP locus was differentially selected among human populations based on the analysis of a nonsynonymous SNP (rs2291725). Comparative and functional analyses showed that the human GIP gene encodes a cryptic glucose-dependent insulinotropic polypeptide (GIP) isoform (GIP55S or GIP55G) that encompasses the SNP and is resistant to serum degradation relative to the known mature GIP peptide. Importantly, we found that GIP55G, which is encoded by the derived allele, exhibits a higher bioactivity compared with GIP55S, which is derived from the ancestral allele. Haplotype structure analysis suggests that the derived allele at rs2291725 arose to dominance in East Asians ∼8100 yr ago due to positive selection. The combined results suggested that rs2291725 represents a functional mutation and may contribute to the population genetics observation. Given that GIP signaling plays a critical role in homeostasis regulation at both the enteroinsular and enteroadipocyte axes, our study highlights the importance of understanding adaptations in energy-balance regulation in the face of the emerging diabetes and obesity epidemics.

Natural selection retains overrepresented out-of-frame stop codons against frameshift peptides in prokaryotes.

BACKGROUND: Out-of-frame stop codons (OSCs) occur naturally in coding sequences of all organisms, providing a mechanism of early termination of translation in incorrect reading frame so that the metabolic cost associated with frameshift events can be reduced. Given such a functional significance, we expect statistically overrepresented OSCs in coding sequences as a result of a widespread selection. Accordingly, we examined available prokaryotic genomes to look for evidence of this selection. RESULTS: The complete genome sequences of 990 prokaryotes were obtained from NCBI GenBank. We found that low G+C content coding sequences contain significantly more OSCs and G+C content at specific codon positions were the principal determinants of OSC usage bias in the different reading frames. To investigate if there is overrepresentation of OSCs, we modeled the trinucleotide and hexanucleotide biases of the coding sequences using Markov models, and calculated the expected OSC frequencies for each organism using a Monte Carlo approach. More than 93% of 342 phylogenetically representative prokaryotic genomes contain excess OSCs. Interestingly the degree of OSC overrepresentation correlates positively with G+C content, which may represent a compensatory mechanism for the negative correlation of OSC frequency with G+C content. We extended the analysis using additional compositional bias models and showed that lower-order bias like codon usage and dipeptide bias could not explain the OSC overrepresentation. The degree of OSC overrepresentation was found to correlate negatively with the optimal growth temperature of the organism after correcting for the G+C% and AT skew of the coding sequence. CONCLUSIONS: The present study uses approaches with statistical rigor to show that OSC overrepresentation is a widespread phenomenon among prokaryotes. Our results support the hypothesis that OSCs carry functional significance and have been selected in the course of genome evolution to act against unintended frameshift occurrences. Some results also hint that OSC overrepresentation being a compensatory mechanism to make up for the decrease in OSCs in high G+C organisms, thus revealing the interplay between two different determinants of OSC frequency.

Genomic and experimental evidence for a potential sexual cycle in the pathogenic thermal dimorphic fungus Penicillium marneffei.

All meiotic genes (except HOP1) and genes encoding putative pheromone processing enzymes, pheromone receptors and pheromone response pathways proteins in Aspergillus fumigatus and Aspergillus nidulans and a putative MAT-1 alpha box mating-type gene were present in the Penicillium marneffei genome. A putative MAT-2 high-mobility group mating-type gene was amplified from a MAT-1 alpha box mating-type gene-negative P. marneffei strain. Among 37 P. marneffei patient strains, MAT-1 alpha box and MAT-2 high-mobility group mating-type genes were present in 23 and 14 isolates, respectively. We speculate that P. marneffei can potentially be a heterothallic fungus that does not switch mating type.