CheekAge: a next-generation buccal epigenetic aging clock associated with lifestyle and health.

Epigenetic aging clocks are computational models that predict age using DNA methylation information. Initially, first-generation clocks were developed to make predictions using CpGs that change with age. Over time, next-generation clocks were created using CpGs that relate to both age and health. Since existing next-generation clocks were constructed in blood, we sought to develop a next-generation clock optimized for prediction in cheek swabs, which are non-invasive and easy to collect. To do this, we collected MethylationEPIC data as well as lifestyle and health information from 8045 diverse adults. Using a novel simulated annealing approach that allowed us to incorporate lifestyle and health factors into training as well as a combination of CpG filtering, CpG clustering, and clock ensembling, we constructed CheekAge, an epigenetic aging clock that has a strong correlation with age, displays high test-retest reproducibility across replicates, and significantly associates with a plethora of lifestyle and health factors, such as BMI, smoking status, and alcohol intake. We validated CheekAge in an internal dataset and multiple publicly available datasets, including samples from patients with progeria or meningioma. In addition to exploring the underlying biology of the data and clock, we provide a free online tool that allows users to mine our methylomic data and predict epigenetic age.

Mode and Tempo of 3D Genome Evolution in Drosophila.

Topologically associating domains (TADs) are thought to play an important role in preventing gene misexpression by spatially constraining enhancer-promoter contacts. The deleterious nature of gene misexpression implies that TADs should, therefore, be conserved among related species. Several early studies comparing chromosome conformation between species reported high levels of TAD conservation; however, more recent studies have questioned these results. Furthermore, recent work suggests that TAD reorganization is not associated with extensive changes in gene expression. Here, we investigate the evolutionary conservation of TADs among 11 species of Drosophila. We use Hi-C data to identify TADs in each species and employ a comparative phylogenetic approach to derive empirical estimates of the rate of TAD evolution. Surprisingly, we find that TADs evolve rapidly. However, we also find that the rate of evolution depends on the chromatin state of the TAD, with TADs enriched for developmentally regulated chromatin evolving significantly slower than TADs enriched for broadly expressed, active chromatin. We also find that, after controlling for differences in chromatin state, highly conserved TADs do not exhibit higher levels of gene expression constraint. These results suggest that, in general, most TADs evolve rapidly and their divergence is not associated with widespread changes in gene expression. However, higher levels of evolutionary conservation and gene expression constraints in TADs enriched for developmentally regulated chromatin suggest that these TAD subtypes may be more important for regulating gene expression, likely due to the larger number of long-distance enhancer-promoter contacts associated with developmental genes.

A set of common buccal CpGs that predict epigenetic age and associate with lifespan-regulating genes.

Epigenetic aging clocks are computational models that use DNA methylation sites to predict age. Since cheek swabs are non-invasive and painless, collecting DNA from buccal tissue is highly desirable. Here, we review 11 existing clocks that have been applied to buccal tissue. Two of these were exclusively trained on adults and, while moderately accurate, have not been used to capture health-relevant differences in epigenetic age. Using 130 common CpGs utilized by two or more existing buccal clocks, we generate a proof-of-concept predictor in an adult methylomic dataset. In addition to accurately estimating age (r = 0.95 and mean absolute error = 3.88 years), this clock predicted that Down syndrome subjects were significantly older relative to controls. A literature and database review of CpG-associated genes identified numerous genes (e.g., CLOCK, ELOVL2, and VGF) and molecules (e.g., alpha-linolenic acid, glycine, and spermidine) reported to influence lifespan and/or age-related disease in model organisms.

Two hundred and five newly assembled mitogenomes provide mixed evidence for rivers as drivers of speciation for Amazonian primates.

Mitochondrial DNA remains a cornerstone for molecular ecology, especially for study species from which high-quality tissue samples cannot be easily obtained. Methods using mitochondrial markers are usually reliant on reference databases, but these are often incomplete. Furthermore, available mitochondrial genomes often lack crucial metadata, such as sampling location, limiting their utility for many analyses. Here, we assembled 205 new mitochondrial genomes for platyrrhine primates, most from the Amazon and with known sampling locations. We present a dated mitogenomic phylogeny based on these samples along with additional published platyrrhine mitogenomes, and use this to assess support for the long-standing riverine barrier hypothesis (RBH), which proposes that river formation was a major driver of speciation in Amazonian primates. Along the Amazon, Negro, and Madeira rivers, we found mixed support for the RBH. While we identified divergences that coincide with a river barrier, only some occur synchronously and also overlap with the proposed dates of river formation. The most compelling evidence is for the Amazon river potentially driving speciation within bearded saki monkeys (Chiropotes spp.) and within the smallest extant platyrrhines, the marmosets and tamarins. However, we also found that even large rivers do not appear to be barriers for some primates, including howler monkeys (Alouatta spp.), uakaris (Cacajao spp.), sakis (Pithecia spp.), and robust capuchins (Sapajus spp.). Our results support a more nuanced, clade-specific effect of riverine barriers and suggest that other evolutionary mechanisms, besides the RBH and allopatric speciation, may have played an important role in the diversification of platyrrhines.

3D genome evolution and reorganization in the Drosophila melanogaster species group.

Topologically associating domains, or TADs, are functional units that organize chromosomes into 3D structures of interacting chromatin. TADs play an important role in regulating gene expression by constraining enhancer-promoter contacts and there is evidence that deletion of TAD boundaries leads to aberrant expression of neighboring genes. While the mechanisms of TAD formation have been well-studied, current knowledge on the patterns of TAD evolution across species is limited. Due to the integral role TADs play in gene regulation, their structure and organization is expected to be conserved during evolution. However, more recent research suggests that TAD structures diverge relatively rapidly. We use Hi-C chromosome conformation capture to measure evolutionary conservation of whole TADs and TAD boundary elements between D. melanogaster and D. triauraria, two early-branching species from the melanogaster species group which diverged ∼15 million years ago. We find that the majority of TADs have been reorganized since the common ancestor of D. melanogaster and D. triauraria, via a combination of chromosomal rearrangements and gain/loss of TAD boundaries. TAD reorganization between these two species is associated with a localized effect on gene expression, near the site of disruption. By separating TADs into subtypes based on their chromatin state, we find that different subtypes are evolving under different evolutionary forces. TADs enriched for broadly expressed, transcriptionally active genes are evolving rapidly, potentially due to positive selection, whereas TADs enriched for developmentally-regulated genes remain conserved, presumably due to their importance in restricting gene-regulatory element interactions. These results provide novel insight into the evolutionary dynamics of TADs and help to reconcile contradictory reports related to the evolutionary conservation of TADs and whether changes in TAD structure affect gene expression.

Comparing the selective landscape of TLR7 and TLR8 across primates reveals unique sites under positive selection in Alouatta.

Among primates, susceptibility to yellow fever (YFV), a single-stranded (ss) RNA virus, ranges from complete resistance to high susceptibility. Howler monkeys (genus Alouatta) are the most susceptible to YFV. In order to identify Alouatta-specific genetic factors that may be responsible for their susceptibility, we collected skin samples from howler monkey museum specimens of the species A. caraya and A. guariba clamitans. We compared the rate of nonsynonymous to synonymous (dN/dS) changes of Toll-like receptor (TLR) 7 and TLR8, the two genes responsible for detecting all ssRNA viruses, across the Primate order. Overall, we found that the TLR7 gene is under stronger purifying selection in howler monkeys compared to other New World and Old World primates, but TLR8 is under the same selective pressure. When we evaluated dN/dS at each codon, we found six codons under positive selection in Alouatta TLR8 and two codons under positive selection in TLR7. The changes in TLR7 are unique to A. guariba clamitans and are found in functionally important regions likely to affect detection of ssRNA viruses by TLR7/TLR8, as well as downstream signaling. These amino acid differences in A. guariba clamitans may play a role in YFV susceptibility. These results have implications for identifying genetic factors affecting YFV susceptibility in primates.

Positively selected variants in functionally important regions of TLR7 in Alouatta guariba clamitans with yellow fever virus exposure in Northern Argentina.

BACKGROUND: In 2007-2009, a major yellow fever virus (YFV) outbreak in Northern Argentina decimated the local howler monkey (Alouatta) population. AIMS: To evaluate whether the surviving howler monkeys possess advantageous genetic variants inherited from monkeys alive prior to the YFV outbreak, we explored the relationship between Toll-like receptor (TLR) 7 and TLR8 gene variation and YFV susceptibility. METHODS: We used samples from Alouatta individuals in Misiones, Argentina alive before the YFV outbreak, individuals that died during the outbreak, and individuals that survived the outbreak and are alive today. We measured genetic divergence between Alouatta YFV exposure groups and evaluated Alouatta-specific substitutions for functional consequences. RESULTS: We did not find different allele frequencies in the post-YFV exposure Alouatta group compared to the pre-exposure group. We identified three nonsynonymous variants in TLR7 in Alouatta guariba clamitans. Two of these substitutions are under positive selection in functionally important regions of the gene. DISCUSSION AND CONCLUSIONS: Our results did not indicate that surviving howler monkey spossess advantageous genetic variants at greater frequency than those alive before the YFV outbreak. However, the positively selected unique coding differences in A. guariba clamitans are in the region important in pathogen detection which may affect YFV resistance. Morework is necessary to fully explore this hypothesis.