Ben Fry

A Bivalent Chromatin Structure Marks Key Developmental Genes in Embryonic Stem Cells

The most highly conserved noncoding elements (HCNEs) in mammalian genomes cluster within regions enriched for genes encoding developmentally important transcription factors (TFs). This suggests that HCNE-rich regions may contain key regulatory controls involved in development. We explored this by examining histone methylation in mouse embryonic stem (ES) cells across 56 large HCNE-rich loci. We identified a specific modification pattern, termed bivalent domains, consisting of large regions of H3 lysine 27 methylation harboring smaller regions of H3 lysine 4 methylation. Bivalent domains tend to coincide with TF genes expressed at low levels. We propose that bivalent domains silence developmental genes in ES cells while keeping them poised for activation. We also found striking correspondences between genome sequence and histone methylation in ES cells, which become notably weaker in differentiated cells. These results highlight the importance of DNA sequence in defining the initial epigenetic landscape and suggest a novel chromatin-based mechanism for maintaining pluripotency.

Distinguishing protein-coding and noncoding genes in the human genome

Although the Human Genome Project was completed 4 years ago, the catalog of human protein-coding genes remains a matter of controversy. Current catalogs list a total of ≈24,500 putative protein-coding genes. It is broadly suspected that a large fraction of these entries are functionally meaningless ORFs present by chance in RNA transcripts, because they show no evidence of evolutionary conservation with mouse or dog. However, there is currently no scientific justification for excluding ORFs simply because they fail to show evolutionary conservation: the alternative hypothesis is that most of these ORFs are actually valid human genes that reflect gene innovation in the primate lineage or gene loss in the other lineages. Here, we reject this hypothesis by carefully analyzing the nonconserved ORFs—specifically, their properties in other primates. We show that the vast majority of these ORFs are random occurrences. The analysis yields, as a by-product, a major revision of the current human catalogs, cutting the number of protein-coding genes to ≈20,500. Specifically, it suggests that nonconserved ORFs should be added to the human gene catalog only if there is clear evidence of an encoded protein. It also provides a principled methodology for evaluating future proposed additions to the human gene catalog. Finally, the results indicate that there has been relatively little true innovation in mammalian protein-coding genes.

Population genetic study of the brain-derived neurotrophic factor (BDNF) gene.

Genetic variants in the brain-derived neurotrophic factor (BDNF) gene, predominantly the functional Val66Met polymorphism, have been associated with risk of bipolar disorder and other psychiatric disorders. However, not all studies support these findings, and overall the evidence for the association of BDNF with disease risk is weak. As differences in population genetic structure between patient samples could cause discrepant or spurious association results, we investigated this possibility by carrying out population genetic analyses of the BDNF genomic region. Substantial variation was detected in BDNF coding region single-nucleotide polymorphism (SNP) allele and haplotype frequencies between 58 global populations, with the derived Met allele of Val66Met ranging in frequency from 0 to 72% across populations. FST analyses to assess diversity in the HapMap populations determined that the Val66Met FST value was at the 99.8th percentile among all SNPs in the genome. As the BDNF population genetic differences may be due to local selection, we performed the long-range haplotype test for selection using 68 SNPs spanning the BDNF genomic region in 12 European-derived pedigrees. Evidence for positive selection was found for a high-frequency Val-carrying haplotype, with a relative extended haplotype homozygosity value above the 99th percentile compared with HapMap data (P=4.6 × 10−4). In conclusion, we observed considerable BDNF allele and haplotype diversity among global populations and evidence for positive selection at the BDNF locus. These phenomena can have a profound impact on the detection of disease susceptibility genes and must be considered in gene association studies of BDNF.

The Case for Selection at CCR5-Δ32

The C-C chemokine receptor 5, 32 base-pair deletion (CCR5-Δ32) allele confers strong resistance to infection by the AIDS virus HIV. Previous studies have suggested that CCR5-Δ32 arose within the past 1,000 y and rose to its present high frequency (5%–14%) in Europe as a result of strong positive selection, perhaps by such selective agents as the bubonic plague or smallpox during the Middle Ages. This hypothesis was based on several lines of evidence, including the absence of the allele outside of Europe and long-range linkage disequilibrium at the locus. We reevaluated this evidence with the benefit of much denser genetic maps and extensive control data. We find that the pattern of genetic variation at CCR5-Δ32 does not stand out as exceptional relative to other loci across the genome. Moreover using newer genetic maps, we estimated that the CCR5-Δ32 allele is likely to have arisen more than 5,000 y ago. While such results can not rule out the possibility that some selection may have occurred at C-C chemokine receptor 5 (CCR5), they imply that the pattern of genetic variation seen atCCR5-Δ32 is consistent with neutral evolution. More broadly, the results have general implications for the design of future studies to detect the signs of positive selection in the human genome.

Searching for signals of evolutionary selection in 168 genes related to immune function

Pathogens have played a substantial role in human evolution, with past infections shaping genetic variation at loci influencing immune function. We selected 168 genes known to be involved in the immune response, genotyped common single nucleotide polymorphisms across each gene in three population samples (CEPH Europeans from Utah, Han Chinese from Guangxi, and Yoruba Nigerians from Southwest Nigeria) and searched for evidence of selection based on four tests for non-neutral evolution: minor allele frequency (MAF), derived allele frequency (DAF), Fst versus heterozygosity and extended haplotype homozygosity (EHH). Six of the 168 genes show some evidence for non-neutral evolution in this initial screen, with two showing similar signals in independent data from the International HapMap Project. These analyses identify two loci involved in immune function that are candidates for having been subject to evolutionary selection, and highlight a number of analytical challenges in searching for selection in genome-wide polymorphism data.

Positive Selection of a Pre-Expansion CAG Repeat of the Human SCA2 Gene

A region of approximately one megabase of human Chromosome 12 shows extensive linkage disequilibrium in Utah residents with ancestry from northern and western Europe. This strikingly large linkage disequilibrium block was analyzed with statistical and experimental methods to determine whether natural selection could be implicated in shaping the current genome structure. Extended Haplotype Homozygosity and Relative Extended Haplotype Homozygosity analyses on this region mapped a core region of the strongest conserved haplotype to the exon 1 of the Spinocerebellar ataxia type 2 gene (SCA2). Direct DNA sequencing of this region of the SCA2 gene revealed a significant association between a pre-expanded allele [(CAG)8CAA(CAG)4CAA(CAG)8] of CAG repeats within exon 1 and the selected haplotype of the SCA2 gene. A significantly negative Tajimas D value (−2.20, p < 0.01) on this site consistently suggested selection on the CAG repeat. This region was also investigated in the three other populations, none of which showed signs of selection. These results suggest that a recent positive selection of the pre-expansion SCA2 CAG repeat has occurred in Utah residents with European ancestry.

Clinical and laboratory predictors of Lassa fever outcome in a dedicated treatment facility in Nigeria: a retrospective, observational cohort study

BACKGROUNDnLassa fever is a viral haemorrhagic disease endemic to west Africa. No large-scale studies exist from Nigeria, where the Lassa virus (LASV) is most diverse. LASV diversity, coupled with host genetic and environmental factors, might cause differences in disease pathophysiology. Small-scale studies in Nigeria suggest that acute kidney injury is an important clinical feature and might be a determinant of survival. We aimed to establish the demographic, clinical, and laboratory factors associated with mortality in Nigerian patients with Lassa fever, and hypothesised that LASV was the direct cause of intrinsic renal damage for a subset of the patients with Lassa fever.nnnMETHODSnWe did a retrospective, observational cohort study of consecutive patients in Nigeria with Lassa fever, who tested positive for LASV with RT-PCR, and were treated in Irrua Specialist Teaching Hospital. We did univariate and multivariate statistical analyses, including logistic regression, of all demographic, clinical, and laboratory variables available at presentation to identify the factors associated with patient mortality.nnnFINDINGSnOf 291 patients treated in Irrua Specialist Teaching Hospital between Jan 3, 2011, and Dec 11, 2015, 284 (98%) had known outcomes (died or survived) and seven (2%) were discharged against medical advice. Overall case-fatality rate was 24% (68 of 284 patients), with a 1·4 times increase in mortality risk for each 10 years of age (p=0·00017), reaching 39% (22 of 57) for patients older than 50 years. Of 284 patients, 81 (28%) had acute kidney injury and 104 (37%) had CNS manifestations and thus both were considered important complications of acute Lassa fever in Nigeria. Acute kidney injury was strongly associated with poor outcome (case-fatality rate of 60% [49 of 81 patients]; odds ratio [OR] 15, p<0·00001). Compared with patients without acute kidney injury, those with acute kidney injury had higher incidence of proteinuria (32 [82%] of 39 patients) and haematuria (29 [76%] of 38) and higher mean serum potassium (4·63 [SD 1·04] mmol/L) and lower blood urea nitrogen to creatinine ratio (8·6 for patients without clinical history of fluid loss), suggesting intrinsic renal damage. Normalisation of creatinine concentration was associated with recovery. Elevated serum creatinine (OR 1·3; p=0·046), aspartate aminotransferase (OR 1·5; p=0·075), and potassium (OR 3·6; p=0·0024) were independent predictors of death.nnnINTERPRETATIONnOur study presents detailed clinical and laboratory data for Nigerian patients with Lassa fever and provides strong evidence for intrinsic renal dysfunction in acute Lassa fever. Early recognition and treatment of acute kidney injury might significantly reduce mortality.nnnFUNDINGnGerman Research Foundation, German Center for Infection Research, Howard Hughes Medical Institute, US National Institutes of Health, and World Bank.

Introducing Processing 2.0

The Processing programming language and environment [RF07] is an open source project initiated by Casey Reas and Ben Fry at the MIT back in 2001, and it is a widely adopted tool in the current practice of creative coding and computer arts and design. Processing has the Design By Numbers project [MA01], developed at the MIT Media Lab by John Maeda, as one of its predecessors. DBNs main motivation was to introduce visual designers and artists to computational design. Processing continues with this initial goal of broadening computer literacy by making programming more accessible for creators who do not necessarily have a formal training in CS such as artists, designers, and hobbyists in general.

Investigating exhibits

For the exhibits panel I planned to have a handful of the exhibitors discuss more about the approach they used for their work, and in particular some aspects that aren’t covered in the traditional submission format employed by DIS and other conferences. Below I chose a handful of questions around some of these themes about which I’m curious, and answer them with respect to my own work in information visualization and the human genome – known as Genomic Cartography – and another project entitled Proce55ing, which is a combination of a learning environment for design and programming, and a sketchbook that supports the development of more advanced work. Investigating Exhibits