Genetic Variation in Human Gene Regulatory Factors Uncovers Regulatory Roles in Local Adaptation and Disease

Differences in gene regulation have been suggested to play essential roles in the evolution of phenotypic changes. Although DNA changes in cis-regulatory elements affect only the regulation of its corresponding gene, variations in gene regulatory factors (trans) can have a broader effect, because the expression of many target genes might be affected. Aiming to better understand how natural selection may have shaped the diversity of gene regulatory factors in human, we assembled a catalog of all proteins involved in controlling gene expression. We found that at least five DNA-binding transcription factor classes are enriched among genes located in candidate regions for selection, suggesting that they might be relevant for understanding regulatory mechanisms involved in human local adaptation. The class of KRAB-ZNFs, zinc-finger (ZNF) genes with a Krüppel-associated box, stands out by first, having the most genes located on candidate regions for positive selection. Second, displaying most nonsynonymous single nucleotide polymorphisms (SNPs) with high genetic differentiation between populations within these regions. Third, having 27 KRAB-ZNF gene clusters with high extended haplotype homozygosity. Our further characterization of nonsynonymous SNPs in ZNF genes located within candidate regions for selection, suggests regulatory modifications that might influence the expression of target genes at population level. Our detailed investigation of three candidate regions revealed possible explanations for how SNPs may influence the prevalence of schizophrenia, eye development, and fertility in humans, among other phenotypes. The genetic variation we characterized here may be responsible for subtle to rough regulatory changes that could be important for understanding human adaptation

Fine-scale analysis of mechanisms and controlling factors in a meiotic recombination hotspot in dogs (canis familiaris)

Meiotic recombination re-shuffles genomes from one generation to the next. In humans and most other mammals, meiotic recombination events are clustered in 1-2 kb wide recombination hotspots, whose locations are determined in trans by the protein PR-domain containing 9 (PRDM9). Mice lacking PRDM9 direct recombination to promoters and functional elements, resulting in meiotic defects. Dogs (Canis familiaris) lack a functional copy of PRDM9, yet linkage data showed that historical recombination events cluster in functional elements, suggesting that there may be a mechanism enabling controlled recombination at these locations, and in the absence of PRDM9. However nothing is known about the de-novo activity of dog recombination hotspots and the patters of recombination resolution in this PRDM9 deficient species. I investigated a dog recombination hotspot for de-novo recombination events using pooled sperm typing, and uncovered high crossover frequencies affecting up to 1 % of sperm. Frequencies can differ by one order of magnitude between dogs. Fine-scale analysis of crossover-breakpoints revealed wide distributions of breaks across up to 10 kb within the hotspot region. I further detect asymmetric breakpoint distributions between crossover orientations and crossover-associated transmission distortion, suggesting biased recombination-initiation or -repair. This work is an elaborate fine-scale dissection of a mammalian PRDM9-independent active recombination hotspot

Psoriasis drug development and GWAS interpretation through in silico analysis of transcription factor binding sites

BackgroundPsoriasis is a cytokine‐mediated skin disease that can be treated effectively with immunosuppressive biologic agents. These medications, however, are not equally effective in all patients and are poorly suited for treating mild psoriasis. To develop more targeted therapies, interfering with transcription factor (TF) activity is a promising strategy.MethodsMeta‐analysis was used to identify differentially expressed genes (DEGs) in the lesional skin from psoriasis patients (n = 237). We compiled a dictionary of 2935 binding sites representing empirically‐determined binding affinities of TFs and unconventional DNA‐binding proteins (uDBPs). This dictionary was screened to identify “psoriasis response elements” (PREs) overrepresented in sequences upstream of psoriasis DEGs.ResultsPREs are recognized by IRF1, ISGF3, NF‐kappaB and multiple TFs with helix‐turn‐helix (homeo) or other all‐alpha‐helical (high‐mobility group) DNA‐binding domains. We identified a limited set of DEGs that encode proteins interacting with PRE motifs, including TFs (GATA3, EHF, FOXM1, SOX5) and uDBPs (AVEN, RBM8A, GPAM, WISP2). PREs were prominent within enhancer regions near cytokine‐encoding DEGs (IL17A, IL19 and IL1B), suggesting that PREs might be incorporated into complex decoy oligonucleotides (cdODNs). To illustrate this idea, we designed a cdODN to concomitantly target psoriasis‐activated TFs (i.e., FOXM1, ISGF3, IRF1 and NF‐kappaB). Finally, we screened psoriasis‐associated SNPs to identify risk alleles that disrupt or engender PRE motifs. This identified possible sites of allele‐specific TF/uDBP binding and showed that PREs are disproportionately disrupted by psoriasis risk alleles.ConclusionsWe identified new TF/uDBP candidates and developed an approach that (i) connects transcriptome informatics to cdODN drug development and (ii) enhances our ability to interpret GWAS findings. Disruption of PRE motifs by psoriasis risk alleles may contribute to disease susceptibility.Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/155494/1/ctm2s4016901500545-sup-0001.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/155494/2/ctm2s4016901500545-sup-0018.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/155494/3/ctm2s4016901500545-sup-0002.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/155494/4/ctm2s4016901500545.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/155494/5/ctm2s4016901500545-sup-0009.pd

A common variant associated with dyslexia reduces expression of the KIAA0319 gene

This work was supported by the Wellcome Trust (MYD, SP, TSS, JCK, RWM, PC, SB, and APM), the Intramural Research Programs of the National Human Genome Research Institute (MYD and EDG) and National Cancer Institute (MPO), and the NIH/Ox-Cam Graduate Partnership Program (MYD).Numerous genetic association studies have implicated the KIAA0319 gene on human chromosome 6p22 in dyslexia susceptibility. The causative variant(s) remains unknown but may modulate gene expression, given that (1) a dyslexia-associated haplotype has been implicated in the reduced expression of KIAA0319, and (2) the strongest association has been found for the region spanning exon 1 of KIAA0319. Here, we test the hypothesis that variant(s) responsible for reduced KIAA0319 expression resides on the risk haplotype close to the gene's transcription start site. We identified seven single-nucleotide polymorphisms on the risk haplotype immediately upstream of KIAA0319 and determined that three of these are strongly associated with multiple reading-related traits. Using luciferase-expressing constructs containing the KIAA0319 upstream region, we characterized the minimal promoter and additional putative transcriptional regulator regions. This revealed that the minor allele of rs9461045, which shows the strongest association with dyslexia in our sample (max p-value = 0.0001), confers reduced luciferase expression in both neuronal and non-neuronal cell lines. Additionally, we found that the presence of this rs9461045 dyslexia-associated allele creates a nuclear protein-binding site, likely for the transcriptional silencer OCT-1. Knocking down OCT-1 expression in the neuronal cell line SHSY5Y using an siRNA restores KIAA0319 expression from the risk haplotype to nearly that seen from the non-risk haplotype. Our study thus pinpoints a common variant as altering the function of a dyslexia candidate gene and provides an illustrative example of the strategic approach needed to dissect the molecular basis of complex genetic traits.PostprintPeer reviewe

FE65 Binds Teashirt, Inhibiting Expression of the Primate-Specific Caspase-4

The Alzheimer disease (AD) amyloid protein precursor (APP) can bind the FE65 adaptor protein and this complex can regulate gene expression. We carried out yeast two-hybrid studies with a PTB domain of FE65, focusing on those genes that might be involved in nuclear signaling, and identified and validated Teashirt proteins as FE65 interacting proteins in neurons. Using reporter systems, we observed that FE65 could simultaneously recruit SET, a component of the inhibitor of acetyl transferase, and Teashirt, which in turn recruited histone deacetylases, to produce a powerful gene-silencing complex. We screened stable cell lines with a macroarray focusing on AD-related genes and identified CASP4, encoding caspase-4, as a target of this silencing complex. Chromatin immunoprecipitation showed a direct interaction of FE65 and Teashirt3 with the promoter region of CASP4. Expression studies in postmortem samples demonstrated decreasing expression of Teashirt and increasing expression of caspase-4 with progressive cognitive decline. Importantly, there were significant increases in caspase-4 expression associated with even the earliest neuritic plaque changes in AD. We evaluated a case-control cohort and observed evidence for a genetic association between the Teashirt genes TSHZ1 and TSHZ3 and AD, with the TSHZ3 SNP genotype correlating with expression of Teashirt3. The results were consistent with a model in which reduced expression of Teashirt3, mediated by genetic or other causes, increases caspase-4 expression, leading to progression of AD. Thus the cell biological, gene expression and genetic data support a role for Teashirt/caspase-4 in AD biology. As caspase-4 shows evidence of being a primate-specific gene, current models of AD and other neurodegenerative conditions may be incomplete because of the absence of this gene in the murine genome

A hybrid next generation transcript sequencing-based approach to identify allelic and homeolog-specific single nucleotide polymorphisms in allotetraploid white clover

peer-reviewedBackground: White clover (Trifolium repens L.) is an allotetraploid species possessing two highly collinear ancestral sub-genomes. The apparent existence of highly similar homeolog copies for the majority of genes in white clover is problematic for the development of genome-based resources in the species. This is especially true for the development of genetic markers based on single nucleotide polymorphisms (SNPs), since it is difficult to distinguish between homeolog-specific and allelic variants. Robust methods for categorising single nucleotide variants as allelic or homeolog-specific in large transcript datasets are required. We illustrate one potential approach in this study. Results: We used 454-pyrosequencing sequencing to generate ~760,000 transcript sequences from an 8th generation white clover inbred line. These were assembled and partially annotated to yield a reference transcript set comprising 71,545 sequences. We subsequently performed Illumina sequencing on three further white clover samples, generating 14 million transcript reads from a mixed sample comprising 24 divergent white clover genotypes, and 50 million reads on two further eighth generation white clover inbred lines. Mapping these reads to the reference transcript set allowed us to develop a significant SNP resource for white clover, and to partition the SNPs from the inbred lines into categories reflecting allelic or homeolog-specific variation. The potential for using haplotype reconstruction and progenitor genome comparison to assign haplotypes to specific ancestral sub-genomes of white clover is demonstrated for sequences corresponding to genes encoding dehydration responsive element binding protein and acyl-coA oxidase. Conclusions: In total, 208,854 independent SNPs in 31,715 reference sequences were discovered, approximately three quarters of which were categorised as representing allelic or homeolog-specific variation using two inbred lines. This represents a significant resource for white clover genomics and genetics studies. We discuss the potential to extend the analysis to identify a “core set” of ancestrally derived homeolog specific variants in white clover.Department of Agriculture Food & the Marine, Ireland - Research Stimulus Fund (RSF 07–566