Polypyrimidine tract‐binding protein is required for the repression of gene expression by all‐trans retinoic acid

All‐trans retinoic acid is a key regulator of early development. High concentrations of retinoic acid interfere with differentiation and migration of neural crest cells. Here we report that a dinucleotide repeat in the cis‐element of Snail2 (previously known as Slug) gene plays a role in repression by all‐trans retinoic acid. We analyzed the cis‐acting regulatory regions of the Xenopus Snail2 gene, whose expression is repressed by all‐trans retinoic acid. The analysis identified a TG/CA repeat as a necessary element for the repression. By performing a yeast one‐hybrid screen, we found that a polypyrimidine tract‐binding protein (PTB), which is known to be a regulator of the alternative splicing of pre‐messenger RNA, binds to the TG/CA repeat. Overexpression and knockdown experiments for PTB in HEK293 cells and Xenopus embryos indicated that PTB is required for repression by retinoic acid. The green fluorescent protein‐PTB fusion protein was localized in the nucleus of 293T cells. In situ hybridization for PTB in Xenopus embryos showed that PTB is expressed at the regions including neural crest at the early stages. Our results indicate that PTB plays a role in the repression of gene expression by retinoic acid through binding to the TG/CA repeats.     

Natural sequence variants of yeast environmental sensors confer cell‐to‐cell expression variability

Living systems may have evolved probabilistic bet hedging strategies that generate cell‐to‐cell phenotypic diversity in anticipation of environmental catastrophes, as opposed to adaptation via a deterministic response to environmental changes. Evolution of bet hedging assumes that genotypes segregating in natural populations modulate the level of intraclonal diversity, which so far has largely remained hypothetical. Using a fluorescent P_met17‐GFP reporter, we mapped four genetic loci conferring to a wild yeast strain an elevated cell‐to‐cell variability in the expression of MET17, a gene regulated by the methionine pathway. A frameshift mutation in the Erc1p transmembrane transporter, probably resulting from a release of laboratory strains from negative selection, reduced P_met17‐GFP expression variability. At a second locus, cis‐regulatory polymorphisms increased mean expression of the Mup1p methionine permease, causing increased expression variability in trans. These results demonstrate that an expression quantitative trait locus (eQTL) can simultaneously have a deterministic effect in cis and a probabilistic effect in trans. Our observations indicate that the evolution of transmembrane transporter genes can tune intraclonal variation and may therefore be implicated in both reactive and anticipatory strategies of adaptation.     

Influence of gene position on the expression divergence of oxidative response genes in intraspecific yeast

Phenotypic variation can arise from differences in the protein coding sequence and in the regulatory elements. However, little is known about the contribution of regulatory difference to the expression divergence, especially the cis and trans regulatory variation to the expression divergence in intraspecific populations. In this study, we used two different yeast strains, BY4743 and RM11‐1a/α, to study the regulatory variation to the expression divergence between BY and RM under oxidative stress condition. Our results indicated that the expression divergence of BY and RM is mainly due to trans regulatory variations under both normal and oxidative stress conditions. However, cis regulatory variation seems to play a very important role in oxidative stress response in yeast because 36% of genes showed an increase in cis regulatory variation effect compared with 13% of genes that showed an increase in trans regulatory variation effect after oxidative stress. Our data also indicated that genes located on the longer arm of the chromosomes are more susceptible to cis variation effect under oxidative stress than genes on the shorter arm of the chromosomes.     

Complex patterns of cis‐regulatory polymorphisms in ebony underlie standing pigmentation variation in Drosophila melanogaster

Pigmentation traits in adult Drosophila melanogaster were used in this study to investigate how phenotypic variations in continuous ecological traits can be maintained in a natural population. First, pigmentation variation in the adult female was measured at seven different body positions in 20 strains from the Drosophila melanogaster Genetic Reference Panel (DGRP) originating from a natural population in North Carolina. Next, to assess the contributions of cis‐regulatory polymorphisms of the genes involved in the melanin biosynthesis pathway, allele‐specific expression levels of four genes were quantified by amplicon sequencing using a 454 GS Junior. Among those genes, ebony was significantly associated with pigmentation intensity of the thoracic segment. Detailed sequence analysis of the gene regulatory regions of this gene indicated that many different functional cis‐regulatory alleles are segregating in the population and that variations outside the core enhancer element could potentially play important roles in the regulation of gene expression. In addition, a slight enrichment of distantly associated SNP pairs was observed in the ~10 kb cis‐regulatory region of ebony, which suggested the presence of interacting elements scattered across the region. In contrast, sequence analysis in the core cis‐regulatory region of tan indicated that SNPs within the region are significantly associated with allele‐specific expression level of this gene. Collectively, the data suggest that the underlying genetic differences in the cis‐regulatory regions that control intraspecific pigmentation variation can be more complex than those of interspecific pigmentation trait differences, where causal genetic changes are typically confined to modular enhancer elements.     

Species‐specific coordinated gene expression and trans‐regulation of larval color pattern in three swallowtail butterflies

SUMMARY The diversity of butterfly larval color pattern has been attracted to people since Darwin's time; however, its molecular mechanisms still remain largely unknown. Larval body markings often differ completely between closely related species under natural selection. The final instar larvae of the swallowtail butterflies Papilio xuthus and Papilio polytes show a green camouflage pattern, whereas those of Papilio machaon show a warning color pattern, although P. xuthus and P. machaon are closely related species. To identify the genes that contribute to species divergence, we compared the expression pattern of eight pigment‐associated genes between three Papilio species. The spatial expression pattern of melanin‐related genes coincided with the species‐specific cuticular markings. We newly found that the combination of bilin‐binding protein and yellow‐related gene (YRG) correlated perfectly with larval blue, yellow, and green coloration. To distinguish whether the interspecific differences in pigment‐associated genes are caused by cis‐regulatory changes or distribution differences in trans‐regulatory proteins, we compared species‐specific mRNA expression in an F1 hybrid specimen. Px‐YRG and Pp‐YRG showed a similar expression pattern, suggesting that the change in expression of YRG is caused mainly by changes in the distribution of trans‐regulatory proteins. Our findings shed light on the gene regulatory networks for butterfly larval color pattern.     

Predominant contribution of cis‐regulatory divergence in the evolution of mouse alternative splicing

Divergence of alternative splicing represents one of the major driving forces to shape phenotypic diversity during evolution. However, the extent to which these divergences could be explained by the evolving cis‐regulatory versus trans‐acting factors remains unresolved. To globally investigate the relative contributions of the two factors for the first time in mammals, we measured splicing difference between C57BL/6J and SPRET/EiJ mouse strains and allele‐specific splicing pattern in their F1 hybrid. Out of 11,818 alternative splicing events expressed in the cultured fibroblast cells, we identified 796 with significant difference between the parental strains. After integrating allele‐specific data from F1 hybrid, we demonstrated that these events could be predominately attributed to cis‐regulatory variants, including those residing at and beyond canonical splicing sites. Contrary to previous observations in Drosophila, such predominant contribution was consistently observed across different types of alternative splicing. Further analysis of liver tissues from the same mouse strains and reanalysis of published datasets on other strains showed similar trends, implying in general the predominant contribution of cis‐regulatory changes in the evolution of mouse alternative splicing.     

Transcribed microsatellite allele lengths are often correlated with gene expression in natural sunflower populations

Microsatellites are common in genomes of most eukaryotic species. Due to their high mutability, an adaptive role for microsatellites has been considered. However, little is known concerning the contribution of microsatellites towards phenotypic variation. We used populations of the common sunflower (Helianthus annuus) at two latitudes to quantify the effect of microsatellite allele length on phenotype at the level of gene expression. We conducted a common garden experiment with seed collected from sunflower populations in Kansas and Oklahoma followed by an RNA‐Seq experiment on 95 individuals. The effect of microsatellite allele length on gene expression was assessed across 3,325 microsatellites that could be consistently scored. Our study revealed 479 microsatellites at which allele length significantly correlates with gene expression (eSTRs). When irregular allele sizes not conforming to the motif length were removed, the number of eSTRs rose to 2,379. The percentage of variation in gene expression explained by eSTRs ranged from 1%–86% when controlling for population and allele‐by‐population interaction effects at the 479 eSTRs. Of these eSTRs, 70.4% are in untranslated regions (UTRs). A gene ontology (GO) analysis revealed that eSTRs are significantly enriched for GO terms associated with cis‐ and trans‐regulatory processes. Our findings suggest that a substantial number of transcribed microsatellites can influence gene expression.     

Simple inheritance,complex regulation: Supergene‐mediated fire ant queen polymorphism

The fire ant Solenopsis invicta exists in two alternate social forms: monogyne nests contain a single reproductive queen and polygyne nests contain multiple reproductive queens. This colony‐level social polymorphism corresponds with individual differences in queen physiology, queen dispersal patterns and worker discrimination behaviours, all evidently regulated by an inversion‐based supergene that spans more than 13 Mb of a “social chromosome,” contains over 400 protein‐coding genes and rarely undergoes recombination. The specific mechanisms by which this supergene influences expression of the many distinctive features that characterize the alternate forms remain almost wholly unknown. To advance our understanding of these mechanisms, we explore the effects of social chromosome genotype and natal colony social form on gene expression in queens sampled as they embarked on nuptial flights, using RNA‐sequencing of brains and ovaries. We observe a large effect of natal social form, that is, of the social/developmental environment, on gene expression profiles, with similarly substantial effects of genotype, including: (a) supergene‐associated gene upregulation, (b) allele‐specific expression and (c) pronounced extra‐supergene trans‐regulatory effects. These findings, along with observed spatial variation in differential and allele‐specific expression within the supergene region, highlight the complex gene regulatory landscape that emerged following divergence of the inversion‐mediated Sb haplotype from its homologue, which presumably largely retained the ancestral gene order. The distinctive supergene‐associated gene expression trajectories we document at the onset of a queen’s reproductive life expand the known record of relevant molecular correlates of a complex social polymorphism and point to putative genetic factors underpinning the alternate social syndromes.     

Role of 5′‐ and 3′‐untranslated regions of mRNAs in human diseases

Protein synthesis is often regulated at the level of initiation of translation, making it a critical step. This regulation occurs by both the cis‐regulatory elements, which are located in the 5′‐ and 3′‐UTRs (untranslated regions), and trans‐acting factors. A breakdown in this regulation machinery can perturb cellular metabolism, leading to various physiological abnormalities. The highly structured UTRs, along with features such as GC‐richness, upstream open reading frames and internal ribosome entry sites, significantly influence the rate of translation of mRNAs. In this review, we discuss how changes in the cis‐regulatory sequences of the UTRs, for example, point mutations and truncations, influence expression of specific genes at the level of translation. Such modifications may tilt the physiological balance from healthy to diseased states, resulting in conditions such as hereditary thrombocythaemia, breast cancer, fragile X syndrome, bipolar affective disorder and Alzheimer's disease. This information tends to establish the crucial role of UTRs, perhaps as much as that of coding sequences, in health and disease.     

The genetic basis of evolutionary change in gene expression levels

The regulation of gene expression is an important determinant of organismal phenotype and evolution. However, the widespread recognition of this fact occurred long after the synthesis of evolution and genetics. Here, we give a brief sketch of thoughts regarding gene regulation in the history of evolution and genetics. We then review the development of genome-wide studies of gene regulatory variation in the context of the location and mode of action of the causative genetic changes. In particular, we review mapping of the genetic basis of expression variation through expression quantitative trait locus studies and measuring the cis/trans component of expression variation in allele-specific expression studies. We conclude by proposing a systematic integration of ideas that combines global mapping studies, cis/trans tests and modern population genetics methodologies, in order to directly estimate the forces acting on regulatory variation within and between species.     

Effects of cis and trans regulatory variations on the expression divergence of heat shock response genes between yeast strains

Phenotypic variation among individuals in a population can be due to DNA sequence variation in protein coding regions or in regulatory elements. Recently, many studies have indicated that mutations in regulatory elements may be the major cause of phenotypic evolution. However, the mechanisms for evolutionary changes in gene expression are still not well understood. Here, we studied the relative roles of cis and trans regulatory changes in Saccharomyces cerevisiae cells to cope with heat stress. It has been found that the expression level of ~ 300 genes was induced at least two fold and that of ~ 500 genes was repressed at least two fold in response to heat shock. From the former set of genes, we randomly selected 65 genes that showed polymorphism(s) between the BY and RM strains for pyrosequencing analysis to explore the relative contributions of cis and trans regulatory variations to the expression divergence between BY and RM. Our data indicated that the expression divergence between BY and RM was mainly due to trans regulatory variations under either the normal condition or the heat stress condition. However, the relative contribution of trans regulatory variation was decreased from 76.9% to 61.5% after the heat shock stress. These results indicated that the cis regulatory variation may play an important role in the adaption to heat stress. In our data, 43.1% (28 genes) of the 65 genes showed the same trend of cis or trans variation effect after the heat shock stress, 35.4% (23 genes) showed an increased cis variation effect and 21.5% (14 genes) showed an increased trans variation effect after the heat shock stress. Thus, our data give insights into the relative roles of cis and trans variations in response to heat shock in yeast.