Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

A gene-rich, transcriptionally active environment and the pre-deposition of repressive marks are predictive of susceptibility to KRAB/KAP1-mediated silencing.

AbstractBACKGROUND: KRAB-ZFPs (Krüppel-associated box domain-zinc finger proteins) are vertebrate-restricted transcriptional repressors encoded in the hundreds by the mouse and human genomes. They act via an essential cofactor, KAP1, which recruits effectors responsible for the formation of facultative heterochromatin. We have recently shown that KRAB/KAP1 can mediate long-range transcriptional repression through heterochromatin spreading, but also demonstrated that this process is at times countered by endogenous influences.METHOD: To investigate this issue further we used an ectopic KRAB-based repressor. This system allowed us to tether KRAB/KAP1 to hundreds of euchromatic sites within genes, and to record its impact on gene expression. We then correlated this KRAB/KAP1-mediated transcriptional effect to pre-existing genomic and chromatin structures to identify specific characteristics making a gene susceptible to repression.RESULTS: We found that genes that were susceptible to KRAB/KAP1-mediated silencing carried higher levels of repressive histone marks both at the promoter and over the transcribed region than genes that were insensitive. In parallel, we found a high enrichment in euchromatic marks within both the close and more distant environment of these genes.CONCLUSION: Together, these data indicate that high levels of gene activity in the genomic environment and the pre-deposition of repressive histone marks within a gene increase its susceptibility to KRAB/KAP1-mediated repression.

Characteristic bimodal profiles of RNA polymerase II at thousands of active mammalian promoters.

BACKGROUND: In mammals, ChIP-seq studies of RNA polymerase II (PolII) occupancy have been performed to reveal how recruitment, initiation and pausing of PolII may control transcription rates, but the focus is rarely on obtaining finely resolved profiles that can portray the progression of PolII through sequential promoter states. RESULTS: Here, we analyze PolII binding profiles from high-coverage ChIP-seq on promoters of actively transcribed genes in mouse and humans. We show that the enrichment of PolII near transcription start sites exhibits a stereotypical bimodal structure, with one peak near active transcription start sites and a second peak 110 base pairs downstream from the first. Using an empirical model that reliably quantifies the spatial PolII signal, gene by gene, we show that the first PolII peak allows for refined positioning of transcription start sites, which is corroborated by mRNA sequencing. This bimodal signature is found both in mouse and humans. Analysis of the pausing-related factors NELF and DSIF suggests that the downstream peak reflects widespread pausing at the +1 nucleosome barrier. Several features of the bimodal pattern are correlated with sequence features such as CpG content and TATA boxes, as well as the histone mark H3K4me3. CONCLUSIONS: We thus show how high coverage DNA sequencing experiments can reveal as-yet unnoticed bimodal spatial features of PolII accumulation that are frequent at individual mammalian genes and reminiscent of transcription initiation and pausing. The initiation-pausing hypothesis is corroborated by evidence from run-on sequencing and immunoprecipitation in other cell types and species.

Mechanisms of reproductive isolation between an ant species of hybrid origin and one of its parents.

The establishment of new species by hybridization is difficult because it requires the development of reproductive isolation (RI) in sympatry to escape the homogenizing effects of gene flow from the parental species. Here we investigated the role of two pre- and two postzygotic mechanisms of RI in a system comprising two interdependent Pogonomyrmex harvester ant lineages (the H1 and H2 lineages) of hybrid origin and one of their parental species (P. rugosus). Similar to most other ants, P. rugosus is characterized by an environmental system of caste determination with female brood developing either into queens or workers depending on nongenetic factors. By contrast, there is a strong genetic component to caste determination in the H1 and H2 lineages because the developmental fate of female brood depends on the genetic origin of the parents, with interlineage eggs developing into workers and intralineage eggs developing into queens. The study of a mixed mating aggregation revealed strong differences in mating flight timing between P. rugosus and the two lineages as a first mechanism of RI. A second important prezygotic mechanism was assortative mating. Laboratory experiments also provided support for one of the two investigated mechanisms of postzygotic isolation. The majority of offspring produced from the few matings between P. rugosus and the lineages aborted at the egg stage. This hybrid inviability was under maternal influence, with hybrids produced by P. rugosus queens being always inviable whereas a small proportion of H2 lineage queens produced large numbers of adult hybrid offspring. Finally, we found no evidence that genetic caste determination acted as a second postzygotic mechanism reducing gene flow between P. rugosus and the H lineages. The few viable P. rugosus-H hybrids were not preferentially shunted into functionally sterile workers but developed into both workers and queens. Overall, these results reveal that the nearly complete (99.5%) RI between P. rugosus and the two hybrid lineages stems from the combination of two typical prezygotic mechanisms (mating time divergence and assortative mating) and one postzygotic mechanism (hybrid inviability).

Ethnic minority-majority asymmetry in national attitudes around the world: A multilevel analysis

Using data from the International Social Survey Programme, this research investigated asymmetric attitudes of ethnic minorities and majorities towards their country and explored the impact of human development, ethnic diversity, and social inequality as country-level moderators of national attitudes. In line with the general hypothesis of ethnic asymmetry, we found that ethnic, linguistic, and religious majorities were more identified with the nation and more strongly endorsed nationalist ideology than minorities (H1, 33 countries). Multilevel analyses revealed that this pattern of asymmetry was moderated by country-level characteristics: the difference between minorities and majorities was greatest in ethnically diverse countries and in egalitarian, low inequality contexts. We also observed a larger positive correlation between ethnic subgroup identification and both national identification and nationalism for majorities than for minorities (H2, 20 countries). A stronger overall relationship between ethnic and national identification was observed in countries with a low level of human development. The greatest minority-majority differences in the relationship between ethnic identification and national attitudes were found in egalitarian countries with a strong welfare state tradition.

Novel H3K4me3 marks are enriched at human- and chimpanzee-specific cytogenetic structures.

Human and chimpanzee genomes are 98.8% identical within comparable sequences. However, they differ structurally in nine pericentric inversions, one fusion that originated human chromosome 2, and content and localization of heterochromatin and lineage-specific segmental duplications. The possible functional consequences of these cytogenetic and structural differences are not fully understood and their possible involvement in speciation remains unclear. We show that subtelomeric regions-regions that have a species-specific organization, are more divergent in sequence, and are enriched in genes and recombination hotspots-are significantly enriched for species-specific histone modifications that decorate transcription start sites in different tissues in both human and chimpanzee. The human lineage-specific chromosome 2 fusion point and ancestral centromere locus as well as chromosome 1 and 18 pericentric inversion breakpoints showed enrichment of human-specific H3K4me3 peaks in the prefrontal cortex. Our results reveal an association between plastic regions and potential novel regulatory elements.

SNCA, MAPT, and GSK3B in Parkinson disease: a gene-gene interaction study.

Background and purpose:  Recent evidence suggests that variation in the SNCA, MAPT, and GSK3B genes interacts in affecting risk for Parkinson disease (PD). In the current study, we attempt to validate previously published findings, evaluating gene-gene interactions between SNCA, MAPT, and GSK3B in association with PD. Methods:  Three Caucasian PD patient-control series from the United States, Ireland, and Norway (combined n = 1020 patients and 1095 controls) were genotyped for SNCA rs356219, MAPT H1/H2-discriminating SNP rs1052553, and GSK3B rs334558 and rs6438552. Results:  Our findings indicate that as previously reported, the SNCA rs356219-G allele and MAPT rs1052553 (H1 haplotype) were both associated with an increased risk of PD, whilst contrary to previous reports, GSK3B variants were not. No pair-wise interaction was observed between SNCA, MAPT, and GSK3B; the risk effects of SNCA rs356219-G and MAPT rs1052553-H1 were seen in a similar manner across genotypes of other variants, with no evidence suggesting synergistic, antagonistic, or deferential effects. Conclusions:  In the Caucasian patient-control series examined, risk for PD was influenced by variation in SNCA and MAPT but not GSK3B. Additionally, those three genes did not interact in determining disease risk.

Implication of DNA methylation and PAX5 factor in the transcriptional regulation of hTERT, the human telomerase reverse transcriptase gene

RESUME La télomérase est une enzyme dite "d'immortalité" qui permet aux cellules de maintenir la longueur de leurs télomères, ce qui confère une capacité de réplication illimitée aux cellules reproductrices et cancéreuses. A l'inverse, les cellules somatiques normales, qui n'expriment pas la télomérase, ont une capacité de réplication limitée. La sous-unité catalytique de la télomérase, hTERT, est définie comme le facteur limitant l'activité télomérasique. Entre activateurs et répresseurs, le rôle de la méthylation de l'ADN et de l'acétylation des histones, de nombreux modèles ont été suggérés. La découverte de l'implication de CTCF dans la régulation transcriptionnelle de hTERT explique en partie le mécanisme de répression de la télomérase dans la plupart des cellules somatiques et sa réactivation dans les cellules tumorales. Dans les cellules télomérase-positives, l'activité inhibitrice de CTCF est bloquée par un mécanisme dépendent ou non de la méthylation. Dans la plupart des carcinomes, une hyperméthylation de la région 5' de hTERT bloque l'effet inhibiteur de CTCF, alors qu'une petite région hypométhylée permet un faible niveau de transcription du gène. Nous avons démontré que la protéine MBD2 se lie spécifiquement sur la région 5' méthylée de hTERT dans différentes lignées cellulaires et qu'elle est impliquée dans la répression partielle de la transcription de hTERT dans les cellules tumorales méthylées. Par contre, nous avons montré que dans les lymphocytes B normaux et néoplasiques, la régulation de hTERT est indépendante de la méthylation. Dans ces cellules, le facteur PAX5 se lie sur la région 5' de hTERT en aval du site d'initiation de la traduction (ATG). L'expression exogène de PAX5 dans les cellules télomérase-négatives active la transcription de hTERT, alors que la répression de PAX5 dans les cellules lymphomateuses inhibe la transcription du gène. PAX5 est donc directement impliqué dans l'activation de l'expression de hTERT dans les lymphocytes B exprimant la télomérase. Ces résultats révèlent des différences entre les niveaux de méthylation de hTERT dans les cellules de carcinomes et les lymphocytes B exprimant la télomérase. La méthylation de hTERT en tant que biomarqueur de cancer a été évaluée, puis appliquée à la détection de métastases. Nous avons ainsi montré que la méthylation de hTERT est positivement corrélée au diagnostic cytologique dans les liquides céphalorachidiens. Nos résultats conduisent à un modèle de régulation de hTERT, qui aide à comprendre comment la transcription de ce gène est régulée par CTCF, avec un mécanisme lié ou non à la méthylation du gène hTERT. La méthylation de hTERT s'est aussi révélée être un nouveau et prometteur biomarqueur de cancer. SUMMARY Human telomerase is an "immortalizing" enzyme that enables cells to maintain telomere length, allowing unlimited replicative capacity to reproductive and cancer cells. Conversely, normal somatic cells that do not express telomerase have a finite replicative capacity. The catalytic subunit of telomerase, hTERT, is defined as the limiting factor for telomerase activity. Between activators and repressors, and the role of DNA methylation and histone acetylation, an abundance of hTERT regulatory models have been suggested. The discovery of the implication of CTCF in the transcriptional regulation of hTERT in part explained the mechanism of silencing of telomerase in most somatic cells and its reactivation in neoplastic cells. In telomerase-positive cells, the inhibitory activity of CTCF is blocked by methylation-dependent and -independent mechanisms. In most carcinoma cells, hypermethylation of the hTERT 5' region has been shown to block the inhibitory effect of CTCF, while a short hypomethylated region allows a low transcription level of the gene. We have demonstrated that MBD2 protein specifically binds the methylated 5' region of hTERT in different cell lines and is therefore involved in the partial repression of hTERT transcription in methylated tumor cells. In contrast, we have shown that in normal and neoplastic B cells, hTERT regulation is methylation-independent. The PAX5 factor has been shown to bind to the hTERT 5'region downstream of the ATG translational start site. Ectopic expression of PAX5 in telomerase-negative cells or repression of PAX5 expression in B lymphoma cells respectively activated and repressed hTERT transcription. Thus, PAX5 is strongly implicated in hTERT expression activation in telomerase-positive B cells. These results reveal differences between the hTERT methylation patterns in telomerase-positive carcinoma cells and telomerase-positive normal B cells. The potential of hTERT methylation as a cancer biomarker was evaluated and applied to the detection of metastasis. We have shown that hTERT methylation correlates with the cytological diagnosis in cerebrospinal fluids. Our results suggest a model of hTERT gene regulation, which helps us to better understand how hTERT transcription is regulated by CTCF in methylation-dependant and independent mechanisms. Our data also indicate that hTERT methylation is a promising new cancer biomarker.

Extensive remodeling of DC function by rapid maturation-induced transcriptional silencing.

The activation, or maturation, of dendritic cells (DCs) is crucial for the initiation of adaptive T-cell mediated immune responses. Research on the molecular mechanisms implicated in DC maturation has focused primarily on inducible gene-expression events promoting the acquisition of new functions, such as cytokine production and enhanced T-cell-stimulatory capacity. In contrast, mechanisms that modulate DC function by inducing widespread gene-silencing remain poorly understood. Yet the termination of key functions is known to be critical for the function of activated DCs. Genome-wide analysis of activation-induced histone deacetylation, combined with genome-wide quantification of activation-induced silencing of nascent transcription, led us to identify a novel inducible transcriptional-repression pathway that makes major contributions to the DC-maturation process. This silencing response is a rapid primary event distinct from repression mechanisms known to operate at later stages of DC maturation. The repressed genes function in pivotal processes--including antigen-presentation, extracellular signal detection, intracellular signal transduction and lipid-mediator biosynthesis--underscoring the central contribution of the silencing mechanism to rapid reshaping of DC function. Interestingly, promoters of the repressed genes exhibit a surprisingly high frequency of PU.1-occupied sites, suggesting a novel role for this lineage-specific transcription factor in marking genes poised for inducible repression.

Impact of Sirtuin 2 knockout on innate immune responses

Purpose/Objective: Histone deacetylases (HDACs) deacetylate histones and transcriptional regulators thereby affecting numerous biological functions. Seven mammalian sirtuins (SIRT1-7) constitute the NAD-dependent class III subfamily of HDACs. Sirtuins are the center of great interest due to their regulatory role in the control of metabolism, ageing and age-related diseases. Up to now, little is known about the influence of sirtuins on immune responses, and nothing about the role of SIRT2. The aim of the study was to analyze the influence of SIRT2 knockout on immune cell development and innate immune responses in vitro and in vivo. Materials and methods: SIRT2 germline knockout were produced on a C57BL/6J background. The cellularity of thymus and spleen was assessed by flow cytometry (n = 3). Bone marrow derived macrophages (BMDMs) and dendritic cells (BMDCs) and splenocytes were stimulated with LPS, Pam3CSK4 lipopeptide, CpG ODN, E. coli, S. aureus, TSST-1, SEB, anti-CD3+ CD28 and concanavalin A (n = 3_8). TNF, IL-2, IL-6, IL-12p40 and IFNc production, SIRT1_7 and CD40 expression, and proliferation were quantified by real time-PCR, ELISA, flow cytometry and H3-thymidine incorporation. Mice (n = 6_16) were challenged with LPS, TNF/D-galactosamine, E. coli and K. pneumonia titrated to cause either mild or severe infections or shock. Blood was collected to quantify cytokines and bacteria. Mortality was checked regularly. Results: SIRT2 is the most expressed sirtuin in macrophages and myeloid DCs. To test whether SIRT2 impacts on innate immune responses, we generated SIRT2 germline knockout mice. SIRT2-/- mice born at the expected Mendelian ratio and develop normally. The proportions and absolute numbers of DN1-4, DP and SP thymocytes, and of T-cells (DN and SP, naı¨ve and memory), B-cells (immature and mature), DCs (cDCs and pDCs) and granulocytes in the spleen are similar in SIRT2+/+ and SIRT2-/- mice. SIRT2+/+ and SIRT2-/- BMDMs, BMDCs and splenocytes produce cytokines (RNA and protein), upregulate CD40, and proliferate to the same extent. SIRT2+/+ and SIRT2-/- mice respond similarly (cytokine blood levels, bacterial counts and mortality) to non-severe and lethal endotoxemia, E. coli peritonitis, K. pneumonia pneumonia and TNF-induced shock. Conclusions: SIRT2 knockout has no dramatic impact on the development of immune cells and on innate immune responses in vitro and in vivo. Considering that SIRT2 may participate to control metabolic homeostasis, we are currently assessing the impact of SIRT2 deficiency on innate immune responses under metabolic stress.

Expression pattern of sirtuins in innate immune cells and influence of Sirtuin 1/2 inhibition on innate immune responses to Toll-like receptor ligands and to infection

Purpose/Objective: The family of histone deacetylases comprises 18 members in mammals, among which seven sirtuins (SIRT1-7). Sirtuins are NADP-dependent enzymes that have been involved in the control of cell metabolism, proliferation and survival. The expression pattern of sirtuins and their influence on host response to microbial infection remain largely unknown. The aim of the study was to analyze the expression of SIRT1-7 and to address the effects of SIRT1/2 inhibition on innate immune responses in vitro and in vivo.. Materials and methods: in vitro: Bone marrow (BM), BM-derived macrophages (BMDMs) and dendritic cells (BMDCs) and RAW 264.7 and J774.1 macrophage cell lines were stimulated for 0, 2, 6 and 18 h with LPS, Pam3CSK4 and CpG ODN. SIRT1-7 mRNA was quantified by real time-PCR. TNF was measured by ELISA. In vivo: BALB/c mice were challenged with LPS (350 lg i.p.) with or without a SIRT1/2 inhibitor. Blood and organs were collected after 0, 1, 4, 8 and 24 h to quantify SIRT1-7 and TNF. Mortality was assessed daily. Results: Bone marrow, macrophages and DCs express, in order of abundance, SIRT2 > > SIRT1, SIRT3 and SIRT6 > SIRT4, SIRT5 and SIRT7. Microbial products decrease the expression of all sirtuins except SIRT6 in a time dependent manner in BMDMs (0_24 h). SIRT2 is the most expressed sirtuin also in the liver, kidney (together with SIRT3) and spleen. Upon LPS challenge, SIRT1, SIRT3, SIRT4 and SIRT7 mRNA levels decrease in the liver (from 4 h to 24 h), whereas SIRT1-7 mRNA levels decrease within 1 h in both kidney and spleen. Pharmacological inhibition of SIRT1/2 decreases TNF production by macrophages stimulated with LPS, Pam3CSK4 and CpG ODN (n = 6; P < 0.001). In agreement, prophylactic treatment with a SIRT1/2 inhibitor decreases TNF production (n = 8; P = 0.04) and increases survival (n = 13, P = 0.03) of mice challenged with LPS. Conclusions: Sirtuins are expressed in innate immune cells. Inhibition of SIRT1/2 activity decreases cytokine production by macrophages and protects from endotoxemia, suggesting that sirtuin inhibitors may represent novel adjunctive therapy for treating inflammatory disorders such as sepsis.

Copy number variation and chromatin structure

The functional consequences of structural variation in the human genome range from adaptation, to phenotypic variation, to predisposition to diseases. Copy number variation (CNV) was shown to influence the phenotype by modifying, in a somewhat dose-dependent manner, the expression of genes that map within them, as well as that of genes located on their flanks. To assess the possible mechanism(s) behind this neighboring effect, we compared histone modification status of cell lines from patients affected by Williams-Beuren, Williams-Beuren region duplication, Smith-Magenis or DiGeorge Syndrome and control individuals using a high-throughput version of chromatin immuno-precipitation method (ChIP), called ChlP-seq. We monitored monomethylation of lysine K20 on histone H4 and trimethylation of lysine K27 on histone H3, as proxies for open and condensed chromatin, respectively. Consistent with the changes in expression levels observed for multiple genes mapping on the entire length of chromosomes affected by structural variants, we also detected regions with modified histone status between samples, up- and downstream from the critical regions, up to the end of the rearranged chromosome. We also gauged the intrachromosomal interactions of these cell lines utilizing chromosome conformation capture (4C-seq) technique. We observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together, possibly forming an interacting cluster with each other and the WBSCR. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We conclude, that large genomic rearrangements can lead to changes in the state of the chromatin spreading far away from the critical region, thus possibly affecting expression globally and as a result modifying the phenotype of the patients. - Les conséquences fonctionnelles des variations structurelles dans le génome humain sont vastes, allant de l'adaptation, en passant par les variations phénotypiques, aux prédispositions à certaines maladies. Il a été démontré que les variations du nombre de copies (CNV) influencent le phénotype en modifiant, d'une manière plus ou moins dose-dépendante, l'expression des gènes se situant à l'intérieur de ces régions, mais également celle des gènes se trouvant dans les régions flanquantes. Afin d'étudier les mécanismes possibles sous-jacents à cet effet de voisinage, nous avons comparé les états de modification des histones dans des lignées cellulaires dérivées de patients atteints du syndrome de Williams-Beuren, de la duplication de la région Williams-Beuren, du syndrome de Smith-Magenis ou du syndrome de Di- George et d'individus contrôles en utilisant une version haut-débit de la méthode d'immunoprécipitation de la chromatine (ChIP), appelée ChIP-seq. Nous avons suivi la mono-méthylation de la lysine K20 sur l'histone H4 et la tri-méthylation de la lysine K27 sur l'histone H3, marqueurs respectifs de la chromatine ouverte et fermée. En accord avec les changements de niveaux d'expression observés pour de multiples gènes tout le long des chromosomes affectés par les CNVs, nous avons aussi détecté des régions présentant des modifications d'histones entre les échantillons, situées de part et d'autre des régions critiques, jusqu'aux extrémités du chromosome réarrangé. Nous avons aussi évalué les interactions intra-chromosomiques ayant lieu dans ces cellules par l'utilisation de la technique de capture de conformation des chromosomes (4C-seq). Nous avons observé qu'un groupe de gènes flanquants la région critique du syndrome de Williams-Beuren (WBSCR) forment souvent une boucle, constituant un groupe d'interactions privilégiées entre ces gènes et la WBSCR. La délétion de la WBSCR perturbe l'expression de ce groupe de gènes flanquants, mais également les interactions à grande échelle entre eux et la région réarrangée. Nous en concluons que les larges réarrangements génomiques peuvent aboutir à des changements de l'état de la chromatine pouvant s'étendre bien plus loin que la région critique, affectant donc potentiellement l'expression de manière globale et ainsi modifiant le phénotype des patients.

Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons.

Phosphate homeostasis was studied in a monocotyledonous model plant through the characterization of the PHO1 gene family in rice (Oryza sativa). Bioinformatics and phylogenetic analysis showed that the rice genome has three PHO1 homologs, which cluster with the Arabidopsis (Arabidopsis thaliana) AtPHO1 and AtPHO1;H1, the only two genes known to be involved in root-to-shoot transfer of phosphate. In contrast to the Arabidopsis PHO1 gene family, all three rice PHO1 genes have a cis-natural antisense transcript located at the 5 ' end of the genes. Strand-specific quantitative reverse transcription-PCR analyses revealed distinct patterns of expression for sense and antisense transcripts for all three genes, both at the level of tissue expression and in response to nutrient stress. The most abundantly expressed gene was OsPHO1;2 in the roots, for both sense and antisense transcripts. However, while the OsPHO1;2 sense transcript was relatively stable under various nutrient deficiencies, the antisense transcript was highly induced by inorganic phosphate (Pi) deficiency. Characterization of Ospho1;1 and Ospho1;2 insertion mutants revealed that only Ospho1;2 mutants had defects in Pi homeostasis, namely strong reduction in Pi transfer from root to shoot, which was accompanied by low-shoot and high-root Pi. Our data identify OsPHO1;2 as playing a key role in the transfer of Pi from roots to shoots in rice, and indicate that this gene could be regulated by its cis-natural antisense transcripts. Furthermore, phylogenetic analysis of PHO1 homologs in monocotyledons and dicotyledons revealed the emergence of a distinct clade of PHO1 genes in dicotyledons, which include members having roles other than long-distance Pi transport.

Coordinated effects of sequence variation on DNA binding, chromatin structure, and transcription.

DNA sequence variation has been associated with quantitative changes in molecular phenotypes such as gene expression, but its impact on chromatin states is poorly characterized. To understand the interplay between chromatin and genetic control of gene regulation, we quantified allelic variability in transcription factor binding, histone modifications, and gene expression within humans. We found abundant allelic specificity in chromatin and extensive local, short-range, and long-range allelic coordination among the studied molecular phenotypes. We observed genetic influence on most of these phenotypes, with histone modifications exhibiting strong context-dependent behavior. Our results implicate transcription factors as primary mediators of sequence-specific regulation of gene expression programs, with histone modifications frequently reflecting the primary regulatory event.