HERC2 rs12913832 modulates human pigmentation by attenuating chromatin-loop formation between a long-range enhancer and the OCA2 promoter

Long-range gene regulation and novel therapeutic applications

An intimate relationship exists between nuclear architecture and gene activity. Unraveling the fine-scale three-dimensional structure of the genome and its impact on gene regulation is a major goal of current epigenetic research, one with direct implications for understanding the molecular mechanisms underlying human phenotypic variation and disease susceptibility. In this context, the novel revolutionary genome editing technologies and emerging new ways to manipulate genome folding offer new promises for the treatment of human disorders.

Genetic study of complex diseases in the post-GWAS era

Genome-wide association studies (GWASs) have identified thousands of genes and genetic variants (mainly SNPs) that contribute to complex diseases in humans. Functional characterization and mechanistic elucidation of these SNPs and genes action are the next major challenge. It has been well established that SNPs altering the amino acids of protein-coding genes can drastically impact protein function, and play an important role in molecular pathogenesis. Functions of regulatory SNPs can be complex and elusive, and involve gene expression regulation through the effect on RNA splicing, transcription factor binding, DNA methylation and miRNA recruitment. In the present review, we summarize the recent progress in our understanding of functional consequences of GWAS-associated non-coding regulatory SNPs, and discuss the application of systems genetics and network biology in the interpretation of GWAS findings.

Development of a novel pink-eyed dilution mouse model showing progressive darkening of the eyes and coat hair with aging

Oca2^p-cas (oculocutaneous albinism II; pink-eyed dilution castaneus) is a coat color mutant gene on mouse chromosome 7 that arose spontaneously in wild Mus musculus castaneus mice. Mice homozygous for Oca2^p-cas usually exhibit pink eyes and gray coat hair on the non-agouti genetic background, and this ordinary phenotype remains unchanged throughout life. During breeding of a mixed strain carrying this gene on the C57BL/6J background, we discovered a novel spontaneous mutation that causes darkening of the eyes and coat hair with aging. In this study, we developed a novel mouse model showing this unique phenotype. Gross observations revealed that the pink eyes and gray coat hair of the novel mutant young mice became progressively darker in color by approximately 3 months after birth. Light and transmission-electron microscopic observations revealed a marked increase in melanin pigmentation of coat hair shafts and choroid of the eye in the novel mice compared to that in the ordinary mice. Sequence analysis of Oca2^p-cas revealed a 4.1-kb deletion involving exons 15 and 16 of its wild-type gene. However, there was no sequence difference between the two types of mutant mice. Mating experiments suggested that the novel mutant phenotype was not inherited in a simple fashion, due to incomplete penetrance. The novel spontaneous mutant mouse is the first example of progressive hair darkening animals and is an essential animal model for understanding of the regulation mechanisms of melanin biosynthesis with aging.

ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters

Chromatin interactions connect distal regulatory elements to target gene promoters guiding stimulus- and lineage-specific transcription. Few factors securing chromatin interactions have so far been identified. Here, by integrating chromatin interaction maps with the large collection of transcription factor-binding profiles provided by the ENCODE project, we demonstrate that the zinc-finger protein ZNF143 preferentially occupies anchors of chromatin interactions connecting promoters with distal regulatory elements. It binds directly to promoters and associates with lineage-specific chromatin interactions and gene expression. Silencing ZNF143 or modulating its DNA-binding affinity using single-nucleotide polymorphisms (SNPs) as a surrogate of site-directed mutagenesis reveals the sequence dependency of chromatin interactions at gene promoters. We also find that chromatin interactions alone do not regulate gene expression. Together, our results identify ZNF143 as a novel chromatin-looping factor that contributes to the architectural foundation of the genome by providing sequence specificity at promoters connected with distal regulatory elements.

Chromatin interactions can connect distal regulatory elements to promoters via protein factors, but few such factors have been identified. Here, the authors show that zinc-finger protein ZNF143 is a sequence-specific chromatin-looping factor that connects promoters with distal regulatory elements.

Allele-specific transcriptional regulation of IRF4 in melanocytes is mediated by chromatin looping of the intronic rs12203592 enhancer to the IRF4 promoter

The majority of significant single-nucleotide polymorphisms (SNPs) identified with genome-wide association studies are located in non-coding regions of the genome; it is therefore possible that they are involved in transcriptional regulation of a nearby gene rather than affecting an encoded protein's function. Previously, it was demonstrated that the SNP rs12203592, located in intron 4 of the IRF4 gene, is strongly associated with human skin pigmentation and modulates an enhancer element that controls expression of IRF4. In our study, we investigated the allele-specific effect of rs12203592 on IRF4 expression in epidermal skin samples and in melanocytic cells from donors of different skin color. We focused on the characteristics and activity of the enhancer, and on long-range chromatin interactions in melanocytic cells homozygous and heterozygous for rs12203592. We found that, irrespective of the trans-activating environment, IRF4 transcription is strongly correlated with the allelic status of rs12203592, the activity of the rs12203592 enhancer and that the chromatin features depend on the rs12203592 genotype. Furthermore, we demonstrate that the rs12203592 enhancer physically interacts with the IRF4 promoter through an allele-dependent chromatin loop, and suggest that subsequent allele-specific activation of IRF4 transcription is stabilized by another allele-specific loop from the rs12203592 enhancer to an additional regulatory element in IRF4. We conclude that the non-coding SNP rs12203592 is located in a regulatory region and affects a wide range of enhancer characteristics, resulting into modulation of the enhancer's activity, its interaction with the IRF4 promoter and subsequent allele-specific transcription of IRF4. Our findings provide another example of a non-coding SNP affecting skin color by modulating enhancer-mediated transcriptional regulation.

Genome-wide transcriptome analysis of human epidermal melanocytes

Because human epidermal melanocytes (HEMs) provide critical protection against skin cancer, sunburn, and photoaging, a genome-wide perspective of gene expression in these cells is vital to understanding human skin physiology. In this study we performed high throughput sequencing of HEMs to obtain a complete data set of transcript sizes, abundances, and splicing. As expected, we found that melanocyte specific genes that function in pigmentation were among the highest expressed genes. We analyzed receptor, ion channel and transcription factor gene families to get a better understanding of the cell signaling pathways used by melanocytes. We also performed a comparative transcriptomic analysis of lightly versus darkly pigmented HEMs and found 16 genes differentially expressed in the two pigmentation phenotypes; of those, only one putative melanosomal transporter (SLC45A2) has known function in pigmentation. In addition, we found 166 transcript isoforms expressed exclusively in one pigmentation phenotype, 17 of which are genes involved in signal transduction. Our melanocyte transcriptome study provides a comprehensive view and may help identify novel pigmentation genes and potential pharmacological targets.

Genetic variation in regulatory DNA elements: the case of OCA2 transcriptional regulation

Mutations within the OCA2 gene or the complete absence of the OCA2 protein leads to oculocutaneous albinism type 2. The OCA2 protein plays a central role in melanosome biogenesis, and it is a strong determinant of the eumelanin content in melanocytes. Transcript levels of the OCA2 gene are strongly correlated with pigmentation intensities. Recent studies demonstrated that the transcriptional level of OCA2 is to a large extent determined by the noncoding SNP rs12913832 located 21.5 kb upstream of the OCA2 gene promoter. In this review, we discuss current hypotheses and the available data on the mechanism of OCA2 transcriptional regulation and how this is influenced by genetic variation. Finally, we will explore how future epigenetic studies can be used to advance our insight into the functional biology that connects genetic variation to human pigmentation.

Prostate cancer risk locus at 8q24 as a regulatory hub by physical interactions with multiple genomic loci across the genome

Chromosome 8q24 locus contains regulatory variants that modulate genetic risk to various cancers including prostate cancer (PC). However, the biological mechanism underlying this regulation is not well understood. Here, we developed a chromosome conformation capture (3C)-based multi-target sequencing technology and systematically examined three PC risk regions at the 8q24 locus and their potential regulatory targets across human genome in six cell lines. We observed frequent physical contacts of this risk locus with multiple genomic regions, in particular, inter-chromosomal interaction with CD96 at 3q13 and intra-chromosomal interaction with MYC at 8q24. We identified at least five interaction hot spots within the predicted functional regulatory elements at the 8q24 risk locus. We also found intra-chromosomal interaction genes PVT1, FAM84B and GSDMC and inter-chromosomal interaction gene CXorf36 in most of the six cell lines. Other gene regions appeared to be cell line-specific, such as RRP12 in LNCaP, USP14 in DU-145 and SMIN3 in lymphoblastoid cell line. We further found that the 8q24 functional domains more likely interacted with genomic regions containing genes enriched in critical pathways such as Wnt signaling and promoter motifs such as E2F1 and TCF3. This result suggests that the risk locus may function as a regulatory hub by physical interactions with multiple genes important for prostate carcinogenesis. Further understanding genetic effect and biological mechanism of these chromatin interactions will shed light on the newly discovered regulatory role of the risk locus in PC etiology and progression.

Cuba: exploring the history of admixture and the genetic basis of pigmentation using autosomal and uniparental markers

We carried out an admixture analysis of a sample comprising 1,019 individuals from all the provinces of Cuba. We used a panel of 128 autosomal Ancestry Informative Markers (AIMs) to estimate the admixture proportions. We also characterized a number of haplogroup diagnostic markers in the mtDNA and Y-chromosome in order to evaluate admixture using uniparental markers. Finally, we analyzed the association of 16 single nucleotide polymorphisms (SNPs) with quantitative estimates of skin pigmentation. In the total sample, the average European, African and Native American contributions as estimated from autosomal AIMs were 72%, 20% and 8%, respectively. The Eastern provinces of Cuba showed relatively higher African and Native American contributions than the Western provinces. In particular, the highest proportion of African ancestry was observed in the provinces of Guantánamo (40%) and Santiago de Cuba (39%), and the highest proportion of Native American ancestry in Granma (15%), Holguín (12%) and Las Tunas (12%). We found evidence of substantial population stratification in the current Cuban population, emphasizing the need to control for the effects of population stratification in association studies including individuals from Cuba. The results of the analyses of uniparental markers were concordant with those observed in the autosomes. These geographic patterns in admixture proportions are fully consistent with historical and archaeological information. Additionally, we identified a sex-biased pattern in the process of gene flow, with a substantially higher European contribution from the paternal side, and higher Native American and African contributions from the maternal side. This sex-biased contribution was particularly evident for Native American ancestry. Finally, we observed that SNPs located in the genes SLC24A5 and SLC45A2 are strongly associated with melanin levels in the sample.

The genome in three dimensions: a new frontier in human brain research

Less than 1.5% of the human genome encodes protein. However, vast portions of the human genome are subject to transcriptional and epigenetic regulation, and many noncoding regulatory DNA elements are thought to regulate the spatial organization of interphase chromosomes. For example, chromosomal "loopings" are pivotal for the orderly process of gene expression, by enabling distal regulatory enhancer or silencer elements to directly interact with proximal promoter and transcription start sites, potentially bypassing hundreds of kilobases of interspersed sequence on the linear genome. To date, however, epigenetic studies in the human brain are mostly limited to the exploration of DNA methylation and posttranslational modifications of the nucleosome core histones. In contrast, very little is known about the regulation of supranucleosomal structures. Here, we show that chromosome conformation capture, a widely used approach to study higher-order chromatin, is applicable to tissue collected postmortem, thereby informing about genome organization in the human brain. We introduce chromosome conformation capture protocols for brain and compare higher-order chromatin structures at the chromosome 6p22.2-22.1 schizophrenia and bipolar disorder susceptibility locus, and additional neurodevelopmental risk genes, (DPP10, MCPH1) in adult prefrontal cortex and various cell culture systems, including neurons derived from reprogrammed skin cells. We predict that the exploration of three-dimensional genome architectures and function will open up new frontiers in human brain research and psychiatric genetics and provide novel insights into the epigenetic risk architectures of regulatory noncoding DNA.

Human skin color is influenced by an intergenic DNA polymorphism regulating transcription of the nearby BNC2 pigmentation gene

Single nucleotide polymorphisms (SNPs) found to be statistically significant when associated with human diseases, and other phenotypes are most often located in non-coding regions of the genome. One example is rs10765819 located in the first intron of the BNC2 gene previously associated with (saturation of) human skin color. Here, we demonstrate that a nearby intergenic SNP (rs12350739) in high linkage disequilibrium with rs10756819 is likely the causal DNA variant for the observed BNC2 skin color association. The highly conserved region surrounding rs12350739 functions as an enhancer element regulating BNC2 transcription in human melanocytes, while the activity of this enhancer element depends on the allelic status of rs12350739. When the rs12350739-AA allele is present, the chromatin at the region surrounding rs12350739 is inaccessible and the enhancer element is only slightly active, resulting in low expression of BNC2, corresponding with light skin pigmentation. When the rs12350739-GG allele is present however, the chromatin at the region surrounding rs12350739 is more accessible and the enhancer is active, resulting in a higher expression of BNC2, corresponding with dark skin pigmentation. Overall, we demonstrate the identification of the functional DNA variant that explains the BNC2 skin color association signal, providing another important step towards further understanding human pigmentation genetics beyond statistical association. We thus deliver a clear example of how an intergenic non-coding DNA variant modulates the regulatory potential of the enhancer element it is located within, which in turn results in allele-dependent differential gene expression affecting variation in common human traits.

A molecular basis for classic blond hair color in Europeans

Hair color differences are among the most obvious examples of phenotypic variation in humans. Although genome-wide association studies (GWAS) have implicated multiple loci in human pigment variation, the causative base-pair changes are still largely unknown. Here we dissect a regulatory region of the KITLG gene (encoding KIT ligand) that is significantly associated with common blond hair color in northern Europeans. Functional tests demonstrate that the region contains a regulatory enhancer that drives expression in developing hair follicles. This enhancer contains a common SNP (rs12821256) that alters a binding site for the lymphoid enhancer-binding factor 1 (LEF1) transcription factor, reducing LEF1 responsiveness and enhancer activity in cultured human keratinocytes. Mice carrying ancestral or derived variants of the human KITLG enhancer exhibit significant differences in hair pigmentation, confirming that altered regulation of an essential growth factor contributes to the classic blond hair phenotype found in northern Europeans.

Implications of the admixture process in skin color molecular assessment

The understanding of the complex genotype-phenotype architecture of human pigmentation has clear implications for the evolutionary history of humans, as well as for medical and forensic practices. Although dozens of genes have previously been associated with human skin color, knowledge about this trait remains incomplete. In particular, studies focusing on populations outside the European-North American axis are rare, and, until now, admixed populations have seldom been considered. The present study was designed to help fill this gap. Our objective was to evaluate possible associations of 18 single nucleotide polymorphisms (SNPs), located within nine genes, and one pseudogene with the Melanin Index (MI) in two admixed Brazilian populations (Gaucho, N = 352; Baiano, N = 148) with different histories of geographic and ethnic colonization. Of the total sample, four markers were found to be significantly associated with skin color, but only two (SLC24A5 rs1426654, and SLC45A2 rs16891982) were consistently associated with MI in both samples (Gaucho and Baiano). Therefore, only these 2 SNPs should be preliminarily considered to have forensic significance because they consistently showed the association independently of the admixture level of the populations studied. We do not discard that the other two markers (HERC2 rs1129038 and TYR rs1126809) might be also relevant to admixed samples, but additional studies are necessary to confirm the real importance of these markers for skin pigmentation. Finally, our study shows associations of some SNPs with MI in a modern Brazilian admixed sample, with possible applications in forensic genetics. Some classical genetic markers in Euro-North American populations are not associated with MI in our sample. Our results point out the relevance of considering population differences in selecting an appropriate set of SNPs as phenotype predictors in forensic practice.

Genome-wide analysis of functional and evolutionary features of tele-enhancers

We investigated sequence features of enhancers separated from their target gene by at least one intermediate gene/exon (named tele-enhancers in this study) and enhancers residing inside their target gene locus. In this study, we used whole genome enhancer maps and gene expression profiles to establish a large panel of tele-enhancers. By contrasting tele-enhancers to proximal enhancers targeting heart genes, we observed that heart tele-enhancers use unique regulatory mechanisms based on the cardiac transcription factors SRF, TEAD, and NKX-2.5, whereas proximal heart enhancers rely on GATA4 instead. A functional analysis showed that tele-enhancers preferentially regulate house-keeping genes and genes with a metabolic role during heart development. In addition, tele-enhancers are significantly more conserved than their proximal counterparts. Similar trends have been observed for non-heart tissues and cell types, suggesting that our findings represent general characteristics of tele-enhancers.

Chromatin states reveal functional associations for globally defined transcription start sites in four human cell lines

BACKGROUND

Deciphering the most common modes by which chromatin regulates transcription, and how this is related to cellular status and processes is an important task for improving our understanding of human cellular biology. The FANTOM5 and ENCODE projects represent two independent large scale efforts to map regulatory and transcriptional features to the human genome. Here we investigate chromatin features around a comprehensive set of transcription start sites in four cell lines by integrating data from these two projects.

RESULTS

Transcription start sites can be distinguished by chromatin states defined by specific combinations of both chromatin mark enrichment and the profile shapes of these chromatin marks. The observed patterns can be associated with cellular functions and processes, and they also show association with expression level, location relative to nearby genes, and CpG content. In particular we find a substantial number of repressed inter- and intra-genic transcription start sites enriched for active chromatin marks and Pol II, and these sites are strongly associated with immediate-early response processes and cell signaling. Associations between start sites with similar chromatin patterns are validated by significant correlations in their global expression profiles.

CONCLUSIONS

The results confirm the link between chromatin state and cellular function for expressed transcripts, and also indicate that active chromatin states at repressed transcripts may poise transcripts for rapid activation during immune response.

Laying a solid foundation for Manhattan--'setting the functional basis for the post-GWAS era'

Genome-wide association studies (GWAS) have identified more than 8900 genetic variants, mainly single-nucleotide polymorphisms (SNPs), associated with hundreds of human traits and diseases, which define risk-associated loci. Variants that map to coding regions can affect protein sequence, translation rate, and alternative splicing, all of which influence protein function. However, the vast majority of sequence variants map to non-coding intergenic and intronic regions, and it has been much more challenging to assess the functional nature of these variants. Recent work annotating the non-coding regions of the genome has contributed to post-GWAS studies by facilitating the identification of the functional targets of risk-associated loci. Many non-coding genetic variants within risk-associated loci alter gene expression by modulating the activity of cis-regulatory elements. We review here these recent findings, discuss their implication for the post-GWAS era, and relate their importance to the interpretation of disease-associated mutations identified through whole-genome sequencing.

HBS1L-MYB intergenic variants modulate fetal hemoglobin via long-range MYB enhancers

Genetic studies have identified common variants within the intergenic region (HBS1L-MYB) between GTP-binding elongation factor HBS1L and myeloblastosis oncogene MYB on chromosome 6q that are associated with elevated fetal hemoglobin (HbF) levels and alterations of other clinically important human erythroid traits. It is unclear how these noncoding sequence variants affect multiple erythrocyte characteristics. Here, we determined that several HBS1L-MYB intergenic variants affect regulatory elements that are occupied by key erythroid transcription factors within this region. These elements interact with MYB, a critical regulator of erythroid development and HbF levels. We found that several HBS1L-MYB intergenic variants reduce transcription factor binding, affecting long-range interactions with MYB and MYB expression levels. These data provide a functional explanation for the genetic association of HBS1L-MYB intergenic polymorphisms with human erythroid traits and HbF levels. Our results further designate MYB as a target for therapeutic induction of HbF to ameliorate sickle cell and β-thalassemia disease severity.

Direct evidence for positive selection of skin, hair, and eye pigmentation in Europeans during the last 5,000 y

Pigmentation is a polygenic trait encompassing some of the most visible phenotypic variation observed in humans. Here we present direct estimates of selection acting on functional alleles in three key genes known to be involved in human pigmentation pathways--HERC2, SLC45A2, and TYR--using allele frequency estimates from Eneolithic, Bronze Age, and modern Eastern European samples and forward simulations. Neutrality was overwhelmingly rejected for all alleles studied, with point estimates of selection ranging from around 2-10% per generation. Our results provide direct evidence that strong selection favoring lighter skin, hair, and eye pigmentation has been operating in European populations over the last 5,000 y.

A polymorphism in IRF4 affects human pigmentation through a tyrosinase-dependent MITF/TFAP2A pathway

Sequence polymorphisms linked to human diseases and phenotypes in genome-wide association studies often affect noncoding regions. A SNP within an intron of the gene encoding Interferon Regulatory Factor 4 (IRF4), a transcription factor with no known role in melanocyte biology, is strongly associated with sensitivity of skin to sun exposure, freckles, blue eyes, and brown hair color. Here, we demonstrate that this SNP lies within an enhancer of IRF4 transcription in melanocytes. The allele associated with this pigmentation phenotype impairs binding of the TFAP2A transcription factor that, together with the melanocyte master regulator MITF, regulates activity of the enhancer. Assays in zebrafish and mice reveal that IRF4 cooperates with MITF to activate expression of Tyrosinase (TYR), an essential enzyme in melanin synthesis. Our findings provide a clear example of a noncoding polymorphism that affects a phenotype by modulating a developmental gene regulatory network.

Association of genetic variants with self-assessed color categories in Brazilians

The Brazilian population was formed by extensive admixture of three different ancestral roots: Amerindians, Europeans and Africans. Our previous work has shown that at an individual level, ancestry, as estimated using molecular markers, was a poor predictor of color in Brazilians. We now investigate if SNPs known to be associated with human skin pigmentation can be used to predict color in Brazilians. For that, we studied the association of fifteen SNPs, previously known to be linked with skin color, in 243 unrelated Brazilian individuals self-identified as White, Browns or Blacks from Rio de Janeiro and 212 unrelated Brazilian individuals self-identified as White or Blacks from São Paulo. The significance of association of SNP genotypes with self-assessed color was evaluated using partial regression analysis. After controlling for ancestry estimates as covariates, only four SNPs remained significantly associated with skin pigmentation: rs1426654 and rs2555364 within SLC24A5, rs16891982 at SLC45A2 and rs1042602 at TYR. These loci are known to be involved in melanin synthesis or transport of melanosomes. We found that neither genotypes of these SNPs, nor their combination with biogeographical ancestry in principal component analysis, could predict self-assessed color in Brazilians at an individual level. However, significant correlations did emerge at group level, demonstrating that even though elements other than skin, eye and hair pigmentation do influence self-assessed color in Brazilians, the sociological act of self-classification is still substantially dependent of genotype at these four SNPs.

From GWAS to function: genetic variation in sodium channel gene enhancer influences electrical patterning

The electrical activity of the heart depends on the correct interplay between key transcription factors and cis-regulatory elements, which together regulate the proper heterogeneous expression of genes encoding for ion channels and other proteins. Genome-wide association studies of ECG parameters implicated genetic variants in the genes for these factors and ion channels modulating conduction and depolarization. Here, we review recent insights into the regulation of localized expression of ion channel genes and the mechanism by which a single-nucleotide polymorphism (SNP) associated with alterations in cardiac conduction patterns in humans affects the transcriptional regulation of the sodium channel genes, SCN5A and SCN10A. The identification of regulatory elements of electrical activity genes helps to explain the impact of genetic variants in non-coding regulatory DNA sequences on regulation of cardiac conduction and the predisposition for cardiac arrhythmias.

Prediction of eye color in the Slovenian population using the IrisPlex SNPs

AIM

To evaluate the accuracy of eye color prediction based on six IrisPlex single nucleotide polymorphisms (SNP) in a Slovenian population sample.

METHODS

Six IrisPlex predictor SNPs (HERC2 - rs12913832, OCA2 - rs1800407, SLC45A2 - rs16891982 and TYR - rs1393350, SLC24A4 - rs12896399, and IRF4 - rs12203592) of 105 individuals were analyzed using single base extension approach and SNaPshot chemistry. The IrisPlex multinomial regression prediction model was used to infer eye color probabilities. The accuracy of the IrisPlex was assessed through the calculation of sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and the area under the receiver characteristic operating curves (AUC).

RESULTS

Blue eye color was observed in 44.7%, brown in 29.6%, and intermediate in 25.7% participants. Prediction accuracy expressed by the AUC was 0.966 for blue, 0.913 for brown, and 0.796 for intermediate eye color. Sensitivity was 93.6% for blue, 58.1% for brown, and 0% for intermediate eye color. Specificity was 93.1% for blue, 89.2% for brown, and 100% for intermediate eye color. PPV was 91.7% for blue and 69.2% for brown color. NPV was 94.7% for blue and 83.5% for brown eye color. These values indicate prediction accuracy comparable to that established in other studies.

CONCLUSION

Blue and brown eye color can be reliably predicted from DNA samples using only six polymorphisms, while intermediate eye color defies prediction, indicating that more research is needed to genetically predict the whole variation of eye color in humans.

Helicase-like transcription factor (Hltf) regulates G2/M transition, Wt1/Gata4/Hif-1a cardiac transcription networks, and collagen biogenesis

HLTF/Hltf regulates transcription, remodels chromatin, and coordinates DNA damage repair. Hltf is expressed in mouse brain and heart during embryonic and postnatal development. Silencing Hltf is semilethal. Seventy-four percent of congenic C57BL/6J Hltf knockout mice died, 75% within 12-24 hours of birth. Previous studies in neonatal (6-8 hour postpartum) brain revealed silencing Hltf disrupted cell cycle progression, and attenuated DNA damage repair. An RNA-Seq snapshot of neonatal heart transcriptome showed 1,536 of 20,000 total transcripts were altered (p < 0.05) - 10 up- and 1,526 downregulated. Pathway enrichment analysis with MetaCore™ showed Hltf’s regulation of the G2/M transition (p=9.726E^-15) of the cell cycle in heart is nearly identical to its role in brain. In addition, Brca1 and 12 members of the Brca1 associated genome surveillance complex are also downregulated. Activation of caspase 3 coincides with transcriptional repression of Bcl-2. Hltf loss caused downregulation of Wt1/Gata4/Hif-1a signaling cascades as well as Myh7b/miR499 transcription. Hltf-specific binding to promoters and/or regulatory regions of these genes was authenticated by ChIP-PCR. Hif-1a targets for prolyl (P4ha1, P4ha2) and lysyl (Plod2) collagen hydroxylation, PPIase enzymes (Ppid, Ppif, Ppil3) for collagen trimerization, and lysyl oxidase (Loxl2) for collagen-elastin crosslinking were downregulated. However, transcription of genes for collagens, fibronectin, Mmps and their inhibitors (Timps) was unaffected. The collective downregulation of genes whose protein products control collagen biogenesis caused disorganization of the interstitial and perivascular myocardial collagen fibrillar network as viewed with picrosirius red-staining, and authenticated with spectral imaging. Wavy collagen bundles in control hearts contrasted with collagen fibers that were thin, short and disorganized in Hltf null hearts. Collagen bundles in Hltf null hearts were tangled and fragmented. Thus, silencing Hltf during heart organogenesis compromised DNA double-strand break repair, and caused aberrant collagen biogenesis altering the structural network that transmits cardiomyocyte force into muscle contraction.

Identification and characterization of functional risk variants for colorectal cancer mapping to chromosome 11q23.1

Genome-wide association studies of colorectal cancer (CRC) have identified a number of common variants associated with modest risk, including rs3802842 at chromosome 11q23.1. Several genes map to this region but rs3802842 does not map to any known transcribed or regulatory sequences. We reasoned, therefore, that rs3802842 is not the functional single-nucleotide polymorphism (SNP), but is in linkage disequilibrium (LD) with a functional SNP(s). We performed ChIP-seq for histone modifications in SW480 and HCT-116 CRC cells, and incorporated ChIP-seq and DNase I hypersensitivity data available through ENCODE within a 137-kb genomic region containing rs3802842 on 11q23.1. We identified SNP rs10891246 in LD with rs3802842 that mapped within a bidirectional promoter region of genes C11orf92 and C11orf93. Following mutagenesis to the risk allele, the promoter demonstrated lower levels of reporter gene expression. A second SNP rs7130173 was identified in LD with rs3802842 that mapped to a candidate enhancer region, which showed strong unidirectional activity in both HCT-116 and SW480 CRC cells. The risk allele of rs7130173 demonstrated reduced enhancer activity compared with the common allele, and reduced nuclear protein binding affinity in electromobility shift assays compared with the common allele suggesting differential transcription factor (TF) binding. SNPs rs10891246 and rs7130173 are on the same haplotype, and expression quantitative trait loci (eQTL) analyses of neighboring genes implicate C11orf53, C11orf92 and C11orf93 as candidate target genes. These data imply that rs10891246 and rs7130173 are functional SNPs mapping to 11q23.1 and that C11orf53, C11orf92 and C11orf93 represent novel candidate target genes involved in CRC etiology.

Natural and orthogonal model for estimating gene-gene interactions applied to cutaneous melanoma

Epistasis, or gene-gene interaction, results from joint effects of genes on a trait; thus, the same alleles of one gene may display different genetic effects in different genetic backgrounds. In this study, we generalized the coding technique of a natural and orthogonal interaction (NOIA) model for association studies along with gene-gene interactions for dichotomous traits and human complex diseases. The NOIA model which has non-correlated estimators for genetic effects is important for estimating influence from multiple loci. We conducted simulations and data analyses to evaluate the performance of the NOIA model. Both simulation and real data analyses revealed that the NOIA statistical model had higher power for detecting main genetic effects and usually had higher power for some interaction effects than the usual model. Although associated genes have been identified for predisposing people to melanoma risk: HERC2 at 15q13.1, MC1R at 16q24.3 and CDKN2A at 9p21.3, no gene-gene interaction study has been fully explored for melanoma. By applying the NOIA statistical model to a genome-wide melanoma dataset, we confirmed the previously identified significantly associated genes and found potential regions at chromosomes 5 and 4 that may interact with the HERC2 and MC1R genes, respectively. Our study not only generalized the orthogonal NOIA model but also provided useful insights for understanding the influence of interactions on melanoma risk.

Candidate genes expressed in human islets and their role in the pathogenesis of type 1 diabetes

In type 1 diabetes (T1D), the insulin-producing β cells are destroyed by an immune-mediated process leading to complete insulin deficiency. There is a strong genetic component in T1D. Genes located in the human leukocyte antigen (HLA) region are the most important genetic determinants of disease, but more than 40 additional loci are known to significantly affect T1D risk. Since most of the currently known genetic candidates have annotated immune cell functions, it is generally considered that most of the genetic susceptibility in T1D is caused by variation in genes affecting immune cell function. Recent studies, however, indicate that most T1D candidate genes are expressed in human islets suggesting that the functions of the genes are not restricted to immune cells, but also play roles in the islets and possibly the β cells. Several candidates change expression levels within the islets following exposure to proinflammatory cytokines highlighting that these genes may be involved in the response of β cells to immune attack. In this review, the compiling evidence that many of the candidate genes are expressed in islets and β cells will be presented. Further, we perform the first systematic human islet expression analysis of all genes located in 50 T1D-associated GWAS loci using a published RNA sequencing dataset. We find that 336 out of 857 genes are expressed in human islets and that many of these interact in protein networks. Finally, the potential pathogenetic roles of some candidate genes will be discussed.

The spatial organization of the human genome

In vivo, the human genome functions as a complex, folded, three-dimensional chromatin polymer. Understanding how the human genome is spatially organized and folded inside the cell nucleus is therefore central to understanding how genes are regulated in normal development and dysregulated in disease. Established light microscopy-based approaches and more recent molecular chromosome conformation capture methods are now combining to give us unprecedented insight into this fascinating aspect of human genomics.

Exploring signatures of positive selection in pigmentation candidate genes in populations of East Asian ancestry

Background

Currently, there is very limited knowledge about the genes involved in normal pigmentation variation in East Asian populations. We carried out a genome-wide scan of signatures of positive selection using the 1000 Genomes Phase I dataset, in order to identify pigmentation genes showing putative signatures of selective sweeps in East Asia. We applied a broad range of methods to detect signatures of selection including: 1) Tests designed to identify deviations of the Site Frequency Spectrum (SFS) from neutral expectations (Tajima’s D, Fay and Wu’s H and Fu and Li’s D* and F*), 2) Tests focused on the identification of high-frequency haplotypes with extended linkage disequilibrium (iHS and Rsb) and 3) Tests based on genetic differentiation between populations (LSBL). Based on the results obtained from a genome wide analysis of 25 kb windows, we constructed an empirical distribution for each statistic across all windows, and identified pigmentation genes that are outliers in the distribution.

Results

Our tests identified twenty genes that are relevant for pigmentation biology. Of these, eight genes (ATRN, EDAR, KLHL7, MITF, OCA2, TH, TMEM33 and TRPM1,) were extreme outliers (top 0.1% of the empirical distribution) for at least one statistic, and twelve genes (ADAM17, BNC2, CTSD, DCT, EGFR, LYST, MC1R, MLPH, OPRM1, PDIA6, PMEL (SILV) and TYRP1) were in the top 1% of the empirical distribution for at least one statistic. Additionally, eight of these genes (BNC2, EGFR, LYST, MC1R, OCA2, OPRM1, PMEL (SILV) and TYRP1) have been associated with pigmentary traits in association studies.

Conclusions

We identified a number of putative pigmentation genes showing extremely unusual patterns of genetic variation in East Asia. Most of these genes are outliers for different tests and/or different populations, and have already been described in previous scans for positive selection, providing strong support to the hypothesis that recent selective sweeps left a signature in these regions. However, it will be necessary to carry out association and functional studies to demonstrate the implication of these genes in normal pigmentation variation.

Variants at the 9p21 locus and melanoma risk

Background

The influence of variants at the 9p21 locus on melanoma risk has been reported through investigation of CDKN2A variants through candidate gene approach as well as by genome wide association studies (GWAS).

Methods

In the present study we genotyped, 25 SNPs that tag 273 variants on chromosome 9p21 in 837 melanoma cases and 1154 controls from Spain. Ten SNPs were selected based on previous associations, reported in GWAS, with either melanocytic nevi or melanoma risk or both. The other 15 SNPs were selected to fine map the CDKN2A gene region.

Results

All the 10 variants selected from the GWAS showed statistically significant association with melanoma risk. Statistically significant association with melanoma risk was also observed for the carriers of the variant T-allele of rs3088440 (540 C>T) at the 3’ UTR of CDKN2A gene with an OR 1.52 (95% CI 1.14-2.04). Interaction analysis between risk associated polymorphisms and previously genotyped MC1R variants, in the present study, did not show any statistically significant association. Statistical significant association was observed for the interaction between phototypes and the rs10811629 (located in intron 5 of MTAP). The strongest association was observed between the homozygous carrier of the A–allele and phototype II with an OR of 15.93 (95% CI 5.34-47.54).

Conclusions

Our data confirmed the association of different variants at chromosome 9p21 with melanoma risk and we also found an association of a variant with skin phototypes.

Genetic analyses of the human eye colours using a novel objective method for eye colour classification

In this study, we present a new objective method for measuring the eye colour on a continuous scale that allows researchers to associate genetic markers with different shades of eye colour. With the use of the custom designed software Digital Iris Analysis Tool (DIAT), the iris was automatically identified and extracted from high resolution digital images. DIAT was made user friendly with a graphical user interface. The software counted the number of blue and brown pixels in the iris image and calculated a Pixel Index of the Eye (PIE-score) that described the eye colour quantitatively. The PIE-score ranged from -1 to 1 (brown to blue). The software eliminated the need for user based interpretation and qualitative eye colour categories. In 94% (570) of 605 analyzed eye images, the iris region was successfully extracted and a PIE-score was calculated. A very high correlation between the PIE-score and the human perception of eye colour was observed. The correlations between the PIE-scores and the six IrisPlex SNPs (HERC2 rs12913832, OCA2 rs1800407, SLC24A4 rs12896399, TYR rs1393350, SLC45A2 rs16891982 and IRF4 rs12203592) were analyzed in 570 individuals. Significant differences (p<10(-6)) in the PIE-scores of the individuals typed as HERC2 rs12913832 G (PIE=0.99) and rs12913832 GA (PIE=-0.71) or A (PIE=-0.87) were observed. We adjusted for the effect of HERC2 rs12913832 and showed that the quantitative PIE-scores were significantly associated with SNPs with minor effects (OCA2 rs1800407, SLC24A4 rs12896399 and TYR rs1393350) on the eye colour. We evaluated the two published prediction models for eye colour (IrisPlex [1] and Snipper[2]) and compared the predictions with the PIE-scores. We found good concordance with the prediction from individuals typed as HERC2 rs12913832 G. However, both methods had difficulties in categorizing individuals typed as HERC2 rs12913832 GA because of the large variation in eye colour in HERC2 rs12913832 GA individuals. With the use of the DIAT software and the PIE-score, it will be possible to automatically compare the iris colour of large numbers of iris images obtained by different studies and to perform large meta-studies that may reveal loci with small effects on the eye colour.

The role of MyoD1 and histone modifications in the activation of muscle enhancers

MyoD1 is a key regulator that orchestrates skeletal muscle differentiation through the regulation of gene expression. Although many studies have focused on its role in transcriptional control at gene promoters, less is known regarding the role of MyoD1 in the assembly of active enhancers. Here, we discuss novel data that point to the ability of MyoD1 to mediate the assembly of active enhancers that augment the transcription of genes essential for muscle development and lineage specification. Based on genome-wide studies of epigenetic marks that typify active enhancers, we recently identified the compendium of distal regulatory elements that dictate transcriptional programs during myogenesis. Superimposition of MyoD1 binding sites upon the locations of muscle enhancers revealed its unequivocal binding to a core region of nearly a third of condition-specific muscle enhancers. Further studies exploring deposition of enhancer-related epigenetic marks in myoblasts lacking MyoD1 demonstrate the dependence of muscle enhancer assembly on the presence of MyoD1. We propose a model wherein MyoD1 mediates recruitment of Set7, H3K4me1, H3K27ac, p300, and RNAP II to MyoD1-bound enhancers to establish condition-specific activation of muscle genes. Moreover, muscle enhancers are modulated through coordinated binding of transcription factors, including c-Jun, Jdp2, Meis, and Runx1, which are recruited to muscle enhancers in a MyoD1-dependent manner. Thus, MyoD1 and enhancer-associated transcription factors function coordinately to assemble and regulate enhancers, thereby augmenting expression of muscle-related genes.

Role of helicase-like transcription factor (hltf) in the G2/m transition and apoptosis in brain

HLTF participates in transcription, chromatin remodeling, DNA damage repair, and tumor suppression. Aside from being expressed in mouse brain during embryonic and postnatal development, little is known about Hltf's functional importance. Splice variant quantification of wild-type neonatal (6-8 hour postpartum) brain gave a ratio of 5:1 for Hltf isoform 1 (exons 1-25) to isoform 2 (exons 1-21 with exon 21 extended via a partial intron retention event). Western analysis showed a close correlation between mRNA and protein expression. Complete loss of Hltf caused encephalomalacia with increased apoptosis, and reduced viability. Sixty-four percent of Hltf null mice died, 48% within 12-24 hours of birth. An RNA-Seq snapshot of the neonatal brain transcriptome showed 341 of 20,000 transcripts were altered (p < 0.05) - 95 up regulated and 246 down regulated. MetaCore^TM enrichment pathway analysis revealed Hltf regulates cell cycle, cell adhesion, and TGF-beta receptor signaling. Hltf's most important role is in the G2/M transition of the cell cycle (p  =  4.672e-7) with an emphasis on transcript availability of major components in chromosome cohesion and condensation. Hltf null brains have reduced transcript levels for Rad21/Scc1, histone H3.3, Cap-E/Smc2, Cap-G/G2, and Aurora B kinase. The loss of Hltf in its yeast Rad5-like role in DNA damage repair is accompanied by down regulation of Cflar, a critical inhibitor of TNFRSF6-mediated apoptosis, and increased (p<0.0001) active caspase-3, an indicator of intrinsic triggering of apoptosis in null brains. Hltf also regulates Smad7/Bambi/Tgf-beta/Bmp5/Wnt10b signaling in brain. ChIP confirmed Hltf binding to consensus sequences in predicted (promoter Scgb3a1 gene) and previously unidentified (P-element on chromosome 7) targets. This study is the first to provide a comprehensive view of Hltf targets in brain. Moreover, it reveals how silencing Hltf disrupts cell cycle progression, and attenuates DNA damage repair.

The convergent evolution of blue iris pigmentation in primates took distinct molecular paths

How many distinct molecular paths lead to the same phenotype? One approach to this question has been to examine the genetic basis of convergent traits, which likely evolved repeatedly under a shared selective pressure. We investigated the convergent phenotype of blue iris pigmentation, which has arisen independently in four primate lineages: humans, blue-eyed black lemurs, Japanese macaques, and spider monkeys. Characterizing the phenotype across these species, we found that the variation within the blue-eyed subsets of each species occupies strongly overlapping regions of CIE L*a*b* color space. Yet whereas Japanese macaques and humans display continuous variation, the phenotypes of blue-eyed black lemurs and their sister species (whose irises are brown) occupy more clustered subspaces. Variation in an enhancer of OCA2 is primarily responsible for the phenotypic difference between humans with blue and brown irises. In the orthologous region, we found no variant that distinguishes the two lemur species or associates with quantitative phenotypic variation in Japanese macaques. Given the high similarity between the blue iris phenotypes in these species and that in humans, this finding implies that evolution has used different molecular paths to reach the same end. Am J Phys Anthropol 151:398–407, 2013.© 2013 Wiley Periodicals, Inc.

Colorful DNA polymorphisms in humans

In this review article we summarize current knowledge on how variation on the DNA level influences human pigmentation including color variation of iris, hair, and skin. We review recent progress in the field of human pigmentation genetics by focusing on the genes and DNA polymorphisms discovered to be involved in determining human pigmentation traits, their association with diseases particularly skin cancers, and their power to predict human eye, hair, and skin colors with potential utilization in forensic investigations.

TRACER: a resource to study the regulatory architecture of the mouse genome

BACKGROUND

Mammalian genes are regulated through the action of multiple regulatory elements, often distributed across large regions. The mechanisms that control the integration of these diverse inputs into specific gene expression patterns are still poorly understood. New approaches enabling the dissection of these mechanisms in vivo are needed.

RESULTS

Here, we describe TRACER (http://tracerdatabase.embl.de), a resource that centralizes information from a large on-going functional exploration of the mouse genome with different transposon-associated regulatory sensors. Hundreds of insertions have been mapped to specific genomic positions, and their corresponding regulatory potential has been documented by analysis of the expression of the reporter sensor gene in mouse embryos. The data can be easily accessed and provides information on the regulatory activities present in a large number of genomic regions, notably in gene-poor intervals that have been associated with human diseases.

CONCLUSIONS

TRACER data enables comparisons with the expression pattern of neighbouring genes, activity of surrounding regulatory elements or with other genomic features, revealing the underlying regulatory architecture of these loci. TRACER mouse lines can also be requested for in vivo transposition and chromosomal engineering, to analyse further regions of interest.

Population structure, migration, and diversifying selection in the Netherlands

Genetic variation in a population can be summarized through principal component analysis (PCA) on genome-wide data. PCs derived from such analyses are valuable for genetic association studies, where they can correct for population stratification. We investigated how to capture the genetic population structure in a well-characterized sample from the Netherlands and in a worldwide data set and examined whether (1) removing long-range linkage disequilibrium (LD) regions and LD-based SNP pruning significantly improves correlations between PCs and geography and (2) whether genetic differentiation may have been influenced by migration and/or selection. In the Netherlands, three PCs showed significant correlations with geography, distinguishing between: (1) North and South; (2) East and West; and (3) the middle-band and the rest of the country. The third PC only emerged with minimized LD, which also significantly increased correlations with geography for the other two PCs. In addition to geography, the Dutch North-South PC showed correlations with genome-wide homozygosity (r=0.245), which may reflect a serial-founder effect due to northwards migration, and also with height (♂: r=0.142, ♀: r=0.153). The divergence between subpopulations identified by PCs is partly driven by selection pressures. The first three PCs showed significant signals for diversifying selection (545 SNPs - the majority within 184 genes). The strongest signal was observed between North and South for the functional SNP in HERC2 that determines human blue/brown eye color. Thus, this study demonstrates how to increase ancestry signals in a relatively homogeneous population and how those signals can reveal evolutionary history.

Genetic architecture of skin and eye color in an African-European admixed population

Variation in human skin and eye color is substantial and especially apparent in admixed populations, yet the underlying genetic architecture is poorly understood because most genome-wide studies are based on individuals of European ancestry. We study pigmentary variation in 699 individuals from Cape Verde, where extensive West African/European admixture has given rise to a broad range in trait values and genomic ancestry proportions. We develop and apply a new approach for measuring eye color, and identify two major loci (HERC2[OCA2] P = 2.3×10⁻⁶², SLC24A5 P = 9.6×10⁻⁹) that account for both blue versus brown eye color and varying intensities of brown eye color. We identify four major loci (SLC24A5 P = 5.4×10⁻²⁷, TYR P = 1.1×10⁻⁹, APBA2[OCA2] P = 1.5×10⁻⁸, SLC45A2 P = 6×10⁻⁹) for skin color that together account for 35% of the total variance, but the genetic component with the largest effect (∼44%) is average genomic ancestry. Our results suggest that adjacent cis-acting regulatory loci for OCA2 explain the relationship between skin and eye color, and point to an underlying genetic architecture in which several genes of moderate effect act together with many genes of small effect to explain ∼70% of the estimated heritability.

Differences in skin and eye color are some of the most obvious traits that underlie human diversity, yet most of our knowledge regarding the genetic basis for these traits is based on the limited range of variation represented by individuals of European ancestry. We have studied a unique population in Cape Verde, an archipelago located off the West African coast, in which extensive mixing between individuals of Portuguese and West African ancestry has given rise to a broad range of phenotypes and ancestral genome proportions. Our results help to explain how genes work together to control the full range of pigmentary phenotypic diversity, provide new insight into the evolution of these traits, and provide a model for understanding other types of quantitative variation in admixed populations.

Identical cells with different 3D genomes; cause and consequences?

The mammalian genome is folded into topological domains, chromosomal units that probably serve to spatially accommodate enhancer-promoter interactions and control gene expression levels across cell populations. Longer-range contacts beyond topological domains are also formed, but only in subpopulations of cells. We propose a model (dog-on-a-lead model) to understand the principles behind and consequences of cell-specific remote DNA contacts and speculate that cell-specific genome topologies can cause variegated gene expression among otherwise identical cells.

Multiplexed chromosome conformation capture sequencing for rapid genome-scale high-resolution detection of long-range chromatin interactions

Chromosome conformation capture (3C) technology is a powerful and increasingly popular tool for analyzing the spatial organization of genomes. Several 3C variants have been developed (e.g., 4C, 5C, ChIA-PET, Hi-C), allowing large-scale mapping of long-range genomic interactions. Here we describe multiplexed 3C sequencing (3C-seq), a 4C variant coupled to next-generation sequencing, allowing genome-scale detection of long-range interactions with candidate regions. Compared with several other available techniques, 3C-seq offers a superior resolution (typically single restriction fragment resolution; approximately 1-8 kb on average) and can be applied in a semi-high-throughput fashion. It allows the assessment of long-range interactions of up to 192 genes or regions of interest in parallel by multiplexing library sequencing. This renders multiplexed 3C-seq an inexpensive, quick (total hands-on time of 2 weeks) and efficient method that is ideal for the in-depth analysis of complex genetic loci. The preparation of multiplexed 3C-seq libraries can be performed by any investigator with basic skills in molecular biology techniques. Data analysis requires basic expertise in bioinformatics and in Linux and Python environments. The protocol describes all materials, critical steps and bioinformatics tools required for successful application of 3C-seq technology.

Genome-wide association study implicates testis-sperm specific FKBP6 as a susceptibility locus for impaired acrosome reaction in stallions

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals. Despite compelling evidence about the genetic control over acrosome biogenesis and function, the genomics of IAR is as yet poorly understood, providing no molecular tools for diagnostics. Here we conducted Equine SNP50 Beadchip genotyping and GWAS using 7 IAR–affected and 37 control Thoroughbred stallions. A significant (P<6.75E-08) genotype–phenotype association was found in horse chromosome 13 in FK506 binding protein 6 (FKBP6). The gene belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. Direct sequencing of FKBP6 exons in cases and controls identified SNPs g.11040315G>A and g.11040379C>A (p.166H>N) in exon 4 that were significantly associated with the IAR phenotype both in the GWAS cohort (n = 44) and in a large multi-breed cohort of 265 horses. All IAR stallions were homozygous for the A-alleles, while this genotype was found only in 2% of controls. The equine FKBP6 was exclusively expressed in testis and sperm and had 5 different transcripts, of which 4 were novel. The expression of this gene in AC/AG heterozygous controls was monoallelic, and we observed a tendency for FKBP6 up-regulation in IAR stallions compared to controls. Because exon 4 SNPs had no effect on the protein structure, it is likely that FKBP6 relates to the IAR phenotype via regulatory or modifying functions. In conclusion, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in mammals. FKBP6 genotyping is recommended for the detection of IAR–susceptible individuals among potential breeding stallions. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals, and currently the molecular causes of the condition are not known. Here we report the mapping, identification, and functional analysis of a susceptibility locus for IAR in stallions. The candidate region was mapped to horse chromosome 13 by SNP genotyping and GWAS of 7 IAR affected and 44 control Thoroughbred stallions. Re-sequencing and case-control analysis of functionally relevant candidate genes in the region identified FKBP6 gene as a significantly associated locus. The association was confirmed by genotyping 265 male horses of multiple breeds. FKBP6 belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. We showed that the equine FKBP6 is exclusively and monoallelically expressed in testis and sperm and has 5 different transcripts, of which 4 were novel. Overall, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in other mammals. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Identifying and characterizing regulatory sequences in the human genome with chromatin accessibility assays

After finishing a human genome reference sequence in 2002, the genomics community has turned to the task of interpreting it. A primary focus is to identify and characterize not only protein-coding genes, but all functional elements in the genome. The effort includes both individual investigators and large-scale projects like the Encyclopedia of DNA Elements (ENCODE) project. As part of the ENCODE project, several groups have identified millions of regulatory elements in hundreds of human cell-types using DNase-seq and FAIRE-seq experiments that detect regions of nucleosome-free open chromatin. ChIP-seq experiments have also been used to discover transcription factor binding sites and map histone modifications. Nearly all identified elements are found in non-coding DNA, hypothesizing a function for previously unannotated sequence. In this review, we provide an overview of the ENCODE effort to define regulatory elements, summarize the main results, and discuss implications of the millions of regulatory elements distributed throughout the genome.

Transcription factor binding at enhancers: shaping a genomic regulatory landscape in flux

The mammalian genome is packed tightly in the nucleus of the cell. This packing is primarily facilitated by histone proteins and results in an ordered organization of the genome in chromosome territories that can be roughly divided in heterochromatic and euchromatic domains. On top of this organization several distinct gene regulatory elements on the same chromosome or other chromosomes are thought to dynamically communicate via chromatin looping. Advances in genome-wide technologies have revealed the existence of a plethora of these regulatory elements in various eukaryotic genomes. These regulatory elements are defined by particular in vitro assays as promoters, enhancers, insulators, and boundary elements. However, recent studies indicate that the in vivo distinction between these elements is often less strict. Regulatory elements are bound by a mixture of common and lineage-specific transcription factors which mediate the long-range interactions between these elements. Inappropriate modulation of the binding of these transcription factors can alter the interactions between regulatory elements, which in turn leads to aberrant gene expression with disease as an ultimate consequence. Here we discuss the bi-modal behavior of regulatory elements that act in cis (with a focus on enhancers), how their activity is modulated by transcription factor binding and the effect this has on gene regulation.

Human pigmentation genes under environmental selection

Genome-wide association studies and comparative genomics have established major loci and specific polymorphisms affecting human skin, hair and eye color. Environmental changes have had an impact on selected pigmentation genes as populations have expanded into different regions of the globe.

The HIrisPlex system for simultaneous prediction of hair and eye colour from DNA

Recently, the field of predicting phenotypes of externally visible characteristics (EVCs) from DNA genotypes with the final aim of concentrating police investigations to find persons completely unknown to investigating authorities, also referred to as Forensic DNA Phenotyping (FDP), has started to become established in forensic biology. We previously developed and forensically validated the IrisPlex system for accurate prediction of blue and brown eye colour from DNA, and recently showed that all major hair colour categories are predictable from carefully selected DNA markers. Here, we introduce the newly developed HIrisPlex system, which is capable of simultaneously predicting both hair and eye colour from DNA. HIrisPlex consists of a single multiplex assay targeting 24 eye and hair colour predictive DNA variants including all 6 IrisPlex SNPs, as well as two prediction models, a newly developed model for hair colour categories and shade, and the previously developed IrisPlex model for eye colour. The HIrisPlex assay was designed to cope with low amounts of template DNA, as well as degraded DNA, and preliminary sensitivity testing revealed full DNA profiles down to 63pg input DNA. The power of the HIrisPlex system to predict hair colour was assessed in 1551 individuals from three different parts of Europe showing different hair colour frequencies. Using a 20% subset of individuals, while 80% were used for model building, the individual-based prediction accuracies employing a prediction-guided approach were 69.5% for blond, 78.5% for brown, 80% for red and 87.5% for black hair colour on average. Results from HIrisPlex analysis on worldwide DNA samples imply that HIrisPlex hair colour prediction is reliable independent of bio-geographic ancestry (similar to previous IrisPlex findings for eye colour). We furthermore demonstrate that it is possible to infer with a prediction accuracy of >86% if a brown-eyed, black-haired individual is of non-European (excluding regions nearby Europe) versus European (including nearby regions) bio-geographic origin solely from the strength of HIrisPlex eye and hair colour probabilities, which can provide extra intelligence for future forensic applications. The HIrisPlex system introduced here, including a single multiplex test assay, an interactive tool and prediction guide, and recommendations for reporting final outcomes, represents the first tool for simultaneously establishing categorical eye and hair colour of a person from DNA. The practical forensic application of the HIrisPlex system is expected to benefit cases where other avenues of investigation, including STR profiling, provide no leads on who the unknown crime scene sample donor or the unknown missing person might be.