Selective sweep for an enhancer involucrin allele identifies skin barrier adaptation out of Africa.

The genetic modules that contribute to human evolution are poorly understood. Here we investigate positive selection in the Epidermal Differentiation Complex locus for skin barrier adaptation in diverse HapMap human populations (CEU, JPT/CHB, and YRI). Using Composite of Multiple Signals and iSAFE, we identify selective sweeps for LCE1A-SMCP and involucrin (IVL) haplotypes associated with human migration out-of-Africa, reaching near fixation in European populations. CEU-IVL is associated with increased IVL expression and a known epidermis-specific enhancer. CRISPR/Cas9 deletion of the orthologous mouse enhancer in vivo reveals a functional requirement for the enhancer to regulate Ivl expression in cis. Reporter assays confirm increased regulatory and additive enhancer effects of CEU-specific polymorphisms identified at predicted IRF1 and NFIC binding sites in the IVL enhancer (rs4845327) and its promoter (rs1854779). Together, our results identify a selective sweep for a cis regulatory module for CEU-IVL, highlighting human skin barrier evolution for increased IVL expression out-of-Africa.

Novel oxytocin receptor variants in laboring women requiring high doses of oxytocin.

BACKGROUND: Although oxytocin commonly is used to augment or induce labor, it is difficult to predict its effectiveness because oxytocin dose requirements vary significantly among women. One possibility is that women requiring high or low doses of oxytocin have variations in the oxytocin receptor gene. OBJECTIVES: To identify oxytocin receptor gene variants in laboring women with low and high oxytocin dosage requirements. STUDY DESIGN: Term, nulliparous women requiring oxytocin doses of ≤4 mU/min (low-dose-requiring, n = 83) or ≥20 mU/min (high-dose-requiring, n = 104) for labor augmentation or induction provided consent to a postpartum blood draw as a source of genomic DNA. Targeted-amplicon sequencing (coverage >30×) with MiSeq (Illumina) was performed to discover variants in the coding exons of the oxytocin receptor gene. Baseline relevant clinical history, outcomes, demographics, and oxytocin receptor gene sequence variants and their allele frequencies were compared between low-dose-requiring and high-dose-requiring women. The Scale-Invariant Feature Transform algorithm was used to predict the effect of variants on oxytocin receptor function. The Fisher exact or χ2 tests were used for categorical variables, and Student t tests or Wilcoxon rank sum tests were used for continuous variables. A P value < .05 was considered statistically significant. RESULTS: The high-dose-requiring women had greater rates of obesity and diabetes and were more likely to have undergone labor induction and required prostaglandins. High-dose-requiring women were more likely to undergo cesarean delivery for first-stage arrest and less likely to undergo cesarean delivery for nonreassuring fetal status. Targeted sequencing of the oxytocin receptor gene in the total cohort (n = 187) revealed 30 distinct coding variants: 17 nonsynonymous, 11 synonymous, and 2 small structural variants. One novel variant (A243T) was found in both the low- and high-dose-requiring groups. Three novel variants (Y106H, A240_A249del, and P197delfs*206) resulting in an amino acid substitution, loss of 9 amino acids, and a frameshift stop mutation, respectively, were identified only in low-dose-requiring women. Nine nonsynonymous variants were unique to the high-dose-requiring group. These included 3 known variants (R151C, G221S, and W228C) and 6 novel variants (M133V, R150L, H173R, A248V, G253R, and I266V). Of these, R150L, R151C, and H173R were predicted by Scale-Invariant Feature Transform algorithm to damage oxytocin receptor function. There was no statistically significant association between the numbers of synonymous and nonsynonymous substitutions in the patient groups. CONCLUSION: Obesity, diabetes, and labor induction were associated with the requirement for high doses of oxytocin. We did not identify significant differences in the prevalence of oxytocin receptor variants between low-dose-requiring and high-dose-requiring women, but novel oxytocin receptor variants were enriched in the high-dose-requiring women. We also found 3 oxytocin receptor variants (2 novel, 1 known) that were predicted to damage oxytocin receptor function and would likely increase an individual's risk for requiring a high oxytocin dose. Further investigation of oxytocin receptor variants and their effects on protein function will inform precision medicine in pregnant women.

Recent Positive Selection in Genes of the Mammalian Epidermal Differentiation Complex Locus.

The epidermal differentiation complex (EDC) is the most rapidly evolving locus in the human genome compared to that of the chimpanzee. Yet the EDC genes that are undergoing positive selection across mammals and in humans are not known. We sought to identify the positively selected genetic variants and determine the evolutionary events of the EDC using mammalian-wide and clade-specific branch- and branch-site likelihood ratio tests and a genetic algorithm (GA) branch test. Significant non-synonymous substitutions were found in filaggrin, SPRR4, LELP1, and S100A2 genes across 14 mammals. By contrast, we identified recent positive selection in SPRR4 in primates. Additionally, the GA branch test discovered lineage-specific evolution for distinct EDC genes occurring in each of the nodes in the 14-mammal phylogenetic tree. Multiple instances of positive selection for FLG, TCHHL1, SPRR4, LELP1, and S100A2 were noted among the primate branch nodes. Branch-site likelihood ratio tests further revealed positive selection in specific sites in SPRR4, LELP1, filaggrin, and repetin across 14 mammals. However, in addition to continuous evolution of SPRR4, site-specific positive selection was also found in S100A11, KPRP, SPRR1A, S100A7L2, and S100A3 in primates and filaggrin, filaggrin2, and S100A8 in great apes. Very recent human positive selection was identified in the filaggrin2 L41 site that was present in Neanderthal. Together, our results identifying recent positive selection in distinct EDC genes reveal an underappreciated evolution of epidermal skin barrier function in primates and humans.

Comparative genomics of Steinernema reveals deeply conserved gene regulatory networks.

BACKGROUND: Parasitism is a major ecological niche for a variety of nematodes. Multiple nematode lineages have specialized as pathogens, including deadly parasites of insects that are used in biological control. We have sequenced and analyzed the draft genomes and transcriptomes of the entomopathogenic nematode Steinernema carpocapsae and four congeners (S. scapterisci, S. monticolum, S. feltiae, and S. glaseri). RESULTS: We used these genomes to establish phylogenetic relationships, explore gene conservation across species, and identify genes uniquely expanded in insect parasites. Protein domain analysis in Steinernema revealed a striking expansion of numerous putative parasitism genes, including certain protease and protease inhibitor families, as well as fatty acid- and retinol-binding proteins. Stage-specific gene expression of some of these expanded families further supports the notion that they are involved in insect parasitism by Steinernema. We show that sets of novel conserved non-coding regulatory motifs are associated with orthologous genes in Steinernema and Caenorhabditis. CONCLUSIONS: We have identified a set of expanded gene families that are likely to be involved in parasitism. We have also identified a set of non-coding motifs associated with groups of orthologous genes in Steinernema and Caenorhabditis involved in neurogenesis and embryonic development that are likely part of conserved protein-DNA relationships shared between these two genera.

Origin and consequences of the relationship between protein mean and variance.

Cell-to-cell variance in protein levels (noise) is a ubiquitous phenomenon that can increase fitness by generating phenotypic differences within clonal populations of cells. An important challenge is to identify the specific molecular events that control noise. This task is complicated by the strong dependence of a protein's cell-to-cell variance on its mean expression level through a power-law like relationship (σ2∝µ1.69). Here, we dissect the nature of this relationship using a stochastic model parameterized with experimentally measured values. This framework naturally recapitulates the power-law like relationship (σ2∝µ1.6) and accurately predicts protein variance across the yeast proteome (r2 = 0.935). Using this model we identified two distinct mechanisms by which protein variance can be increased. Variables that affect promoter activation, such as nucleosome positioning, increase protein variance by changing the exponent of the power-law relationship. In contrast, variables that affect processes downstream of promoter activation, such as mRNA and protein synthesis, increase protein variance in a mean-dependent manner following the power-law. We verified our findings experimentally using an inducible gene expression system in yeast. We conclude that the power-law-like relationship between noise and protein mean is due to the kinetics of promoter activation. Our results provide a framework for understanding how molecular processes shape stochastic variation across the genome.

Regulation of the dynamic chromatin architecture of the epidermal differentiation complex is mediated by a c-Jun/AP-1-modulated enhancer.

The epidermal differentiation complex (EDC) locus comprises a syntenic and linear cluster of genes whose concomitant expression is a hallmark feature of differentiation in the developing skin epidermis. Many of the EDC proteins are cross-linked together to form the cornified envelope, an essential and discrete unit of the mammalian skin barrier. The mechanism underlying coordinate transcriptional activation of the EDC is unknown. Within the human EDC, we identified an epidermal-specific regulatory enhancer, 923, which responded to the developmental and spatiotemporal cues at the onset of epidermal differentiation in the mouse embryo. Comparative chromosomal conformation capture assays in proliferating and differentiated primary mouse keratinocytes revealed multiple physiologically sensitive chromatin interactions between the 923 enhancer and EDC gene promoters, thus depicting the dynamic chromatin topology of the EDC. We elucidate a mechanistic link between c-Jun/AP-1 and 923, whereby AP-1- and 923-mediated EDC chromatin remodeling are required for functional EDC gene activation. Thus, we identify a critical enhancer/transcription factor axis governing the dynamic regulation of the EDC chromatin architecture and gene expression and provide a framework for future studies toward understanding gene regulation in cutaneous diseases.