Regulation of DNA methylation dictates Cd4 expression during the development of helper and cytotoxic T cell lineages.

During development, progenitor cells with binary potential give rise to daughter cells that have distinct functions. Heritable epigenetic mechanisms then lock in gene-expression programs that define lineage identity. Regulation of the gene encoding the T cell-specific coreceptor CD4 in helper and cytotoxic T cells exemplifies this process, with enhancer- and silencer-regulated establishment of epigenetic memory for stable gene expression and repression, respectively. Using a genetic screen, we identified the DNA-methylation machinery as essential for maintaining silencing of Cd4 in the cytotoxic lineage. Furthermore, we found a requirement for the proximal enhancer in mediating the removal of DNA-methylation marks from Cd4, which allowed stable expression of Cd4 in helper T cells. Our findings suggest that stage-specific methylation and demethylation events in Cd4 regulate its heritable expression in response to the distinct signals that dictate lineage 'choice' during T cell development.

Returning integrated genomic risk and clinical recommendations: The eMERGE study.

PURPOSE: Assessing the risk of common, complex diseases requires consideration of clinical risk factors as well as monogenic and polygenic risks, which in turn may be reflected in family history. Returning risks to individuals and providers may influence preventive care or use of prophylactic therapies for those individuals at high genetic risk. METHODS: To enable integrated genetic risk assessment, the eMERGE (electronic MEdical Records and GEnomics) network is enrolling 25,000 diverse individuals in a prospective cohort study across 10 sites. The network developed methods to return cross-ancestry polygenic risk scores, monogenic risks, family history, and clinical risk assessments via a genome-informed risk assessment (GIRA) report and will assess uptake of care recommendations after return of results. RESULTS: GIRAs include summary care recommendations for 11 conditions, education pages, and clinical laboratory reports. The return of high-risk GIRA to individuals and providers includes guidelines for care and lifestyle recommendations. Assembling the GIRA required infrastructure and workflows for ingesting and presenting content from multiple sources. Recruitment began in February 2022. CONCLUSION: Return of a novel report for communicating monogenic, polygenic, and family history-based risk factors will inform the benefits of integrated genetic risk assessment for routine health care.

Genome-wide DNA methylation analysis implicates enrichment of interferon pathway in African American patients with Systemic Lupus Erythematosus and European Americans with lupus nephritis.

Systemic Lupus Erythematosus (SLE) is a chronic, multisystem, inflammatory autoimmune disease that disproportionately affects women. Trends in SLE prevalence and clinical course differ by ancestry, with those of African American ancestry presenting with more active, severe and rapidly progressive disease than European Americans. Previous research established altered epigenetic signatures in SLE patients compared to controls. However, the contribution of aberrant DNA methylation (DNAm) to the risk of SLE by ancestry and differences among patients with SLE-associated Lupus Nephritis (LN) has not been well described. We evaluated the DNA methylomes of 87 individuals including 41 SLE patients, with and without LN, and 46 controls enrolled in an ancestry diverse, well-characterized cohort study of established SLE (41 SLE patients [20 SLE-LN+, 21 SLE-LN-] and 46 sex-, race- and age-matched controls; 55% African American, 45% European American). Participants were genotyped using the Infinium Global Diversity Array (GDA), and genetic ancestry was estimated using principal components. Genome-wide DNA methylation was initially measured using the Illumina MethylationEPIC 850K Beadchip array followed by methylation-specific qPCR to validate the methylation status at putative loci. Differentially Methylated Positions (DMP) were identified using a case-control approach adjusted for ancestry. We identified a total of 51 DMPs in CpGs among SLE patients compared to controls. Genes proximal to these CpGs were highly enriched for involvement in type I interferon signaling. DMPs among European American SLE patients with LN were similar to African American SLE patients with and without LN. Our findings were validated using an orthogonal, methyl-specific PCR for three SLE-associated DMPs near or proximal to MX1, USP18, and IFITM1. Our study confirms previous reports that DMPs in CpGs associated with SLE are enriched in type I interferon genes. However, we show that European American SLE patients with LN have similar DNAm patterns to African American SLE patients irrespective of LN, suggesting that aberrant DNAm alters activity of type I interferon pathway leading to more severe disease independent of ancestry.

Genetic associations with age at dementia onset in the PSEN1 E280A Colombian kindred.

INTRODUCTION: Genetic associations with Alzheimer's disease (AD) age at onset (AAO) could reveal genetic variants with therapeutic applications. We present a large Colombian kindred with autosomal dominant AD (ADAD) as a unique opportunity to discover AAO genetic associations. METHODS: A genetic association study was conducted to examine ADAD AAO in 340 individuals with the PSEN1 E280A mutation via TOPMed array imputation. Replication was assessed in two ADAD cohorts, one sporadic early-onset AD study and four late-onset AD studies. RESULTS: 13 variants had p<1×10-7 or p<1×10-5 with replication including three independent loci with candidate associations with clusterin including near CLU. Other suggestive associations were identified in or near HS3ST1, HSPG2, ACE, LRP1B, TSPAN10, and TSPAN14. DISCUSSION: Variants with suggestive associations with AAO were associated with biological processes including clusterin, heparin sulfate, and amyloid processing. The detection of these effects in the presence of a strong mutation for ADAD reinforces their potentially impactful role.

Skewed allelic expression on X chromosome associated with aberrant expression of XIST on systemic lupus erythematosus lymphocytes.

A common feature of autoimmune diseases, including systemic lupus erythematosus (SLE), is an increased prevalence in women. However, the molecular basis for sex disparity in SLE remains poorly understood. To examine the role of X-linked transcription in SLE adaptive immune cells, we performed RNA-seq in T cell and B cell subsets from either healthy donors or patients with SLE. Analyses of allelic expression (AE) profiles identified a pattern of increased allelic imbalance across the entire X chromosome in SLE lymphocytes. X-linked genes exhibiting AE in SLE had an extensive overlap with genes known to escape X chromosome inactivation (XCI). XIST RNA was overexpressed in SLE patients. Differential XIST expression correlated with AE profiles more positively at X-linked genes than the genome-wide background. Analysis of three independent RNA-seq data verified the XIST-associated skewed AE on X chromosome in SLE. Integrative analyses of DNA methylation profiles showed an increased variability of DNA methylation levels at these AE-related X-linked genes. In cultured lymphoblastic cells, knockdown of XIST specifically altered allelic imbalance patterns between X chromosomes. Our study provides genetic evidence that upregulation of XIST accompanied with more skewed allelic expression on X chromosome is associated with the pathogenesis of SLE and may provide mechanistic insights into the increased incidence of SLE in females.

Genetic influences on susceptibility to rheumatoid arthritis in African-Americans.

Large meta-analyses of rheumatoid arthritis (RA) susceptibility in European (EUR) and East Asian (EAS) populations have identified >100 RA risk loci, but genome-wide studies of RA in African-Americans (AAs) are absent. To address this disparity, we performed an analysis of 916 AA RA patients and 1392 controls and aggregated our data with genotyping data from >100 000 EUR and Asian RA patients and controls. We identified two novel risk loci that appear to be specific to AAs: GPC5 and RBFOX1 (PAA < 5 × 10-9). Most RA risk loci are shared across different ethnicities, but among discordant loci, we observed strong enrichment of variants having large effect sizes. We found strong evidence of effect concordance for only 3 of the 21 largest effect index variants in EURs. We used the trans-ethnic fine-mapping algorithm PAINTOR3 to prioritize risk variants in >90 RA risk loci. Addition of AA data to those of EUR and EAS descent enabled identification of seven novel high-confidence candidate pathogenic variants (defined by posterior probability > 0.8). In summary, our trans-ethnic analyses are the first to include AAs, identified several new RA risk loci and point to candidate pathogenic variants that may underlie this common autoimmune disease. These findings may lead to better ways to diagnose or stratify treatment approaches in RA.

Identifying rare, medically relevant variation via population-based genomic screening in Alabama: opportunities and pitfalls.

PURPOSE: To evaluate the effectiveness and specificity of population-based genomic screening in Alabama. METHODS: The Alabama Genomic Health Initiative (AGHI) has enrolled and evaluated 5369 participants for the presence of pathogenic/likely pathogenic (P/LP) variants using the Illumina Global Screening Array (GSA), with validation of all P/LP variants via Sanger sequencing in a CLIA-certified laboratory before return of results. RESULTS: Among 131 variants identified by the GSA that were evaluated by Sanger sequencing, 67 (51%) were false positives (FP). For 39 of the 67 FP variants, a benign/likely benign variant was present at or near the targeted P/LP variant. Variants detected within African American individuals were significantly enriched for FPs, likely due to a higher rate of nontargeted alternative alleles close to array-targeted P/LP variants. CONCLUSION: In AGHI, we have implemented an array-based process to screen for highly penetrant genetic variants in actionable disease genes. We demonstrate the need for clinical validation of array-identified variants in direct-to-consumer or population testing, especially for diverse populations.

CYP450 polymorphisms and clinical pharmacogenetics of ibuprofen after lower third molar extraction.

PURPOSE: This study hypothesized that drugs accumulate in the bloodstream of poor-metabolizing patients and may have more adverse effects and different pain perceptions and aimed to investigate the influence of CYP450 polymorphisms on acute postoperative pain, swelling, and trismus controlled by ibuprofen (600 mg) in 200 volunteers after dental extraction. In addition, surgical outcomes can determine pain, edema, and trismus and indicate inflammatory reactions after oral surgeries. METHODS: Genetic sequencing was performed to identify CYP450 polymorphisms and the surgical parameters evaluated: pre and postoperative swelling, trismus, and temperature; self-reported postoperative pain with visual analog scale (VAS); rescue medication consumed; and severity of adverse reactions. RESULTS: A multiple linear regression model with independent variables [single nucleotide polymorphisms (SNPs), BMI (body mass index), duration, and difficulty of surgery] and dependent variables [postoperative pain by sum of pain intensity difference (SPID), trismus, and swelling] was used for analysis. The duration of surgery was a predictor for pain at 8 h and 96 h after surgery, and BMI was a predictor for both swelling and trismus on the 2nd postoperative day. When evaluating CYP2C8 and C9 genotyped SNPs, it was observed that normal metabolizers showed higher pain levels than the intermediate/poor metabolizers on the postoperative periods as compared with time 0 h. In another analysis, the poor metabolizers for CYP2C8 and C9 presented lower levels of postoperative pain after 8 h and used rescue medication earlier than normal metabolizers. CONCLUSION: Ibuprofen 600 mg was very effective in controlling inflammatory pain after lower third molar surgeries, without relevant adverse reactions; although in a very subtle way, patients with poor metabolism had higher levels of pain in the first hours, and no longer after 8 h, and used pain relief medication earlier. TRIAL REGISTRATION: The study was registered with ClinicalTrials.gov ID (NCT03169127), on March 16th, 2017.

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Genetic studies of body mass index yield new insights for obesity biology.

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in up to 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 × 10(-8)), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for ∼2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous system in obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

Blood Leukocyte DNA Methylation Predicts Risk of Future Myocardial Infarction and Coronary Heart Disease.

BACKGROUND: DNA methylation is implicated in coronary heart disease (CHD), but current evidence is based on small, cross-sectional studies. We examined blood DNA methylation in relation to incident CHD across multiple prospective cohorts. METHODS: Nine population-based cohorts from the United States and Europe profiled epigenome-wide blood leukocyte DNA methylation using the Illumina Infinium 450k microarray, and prospectively ascertained CHD events including coronary insufficiency/unstable angina, recognized myocardial infarction, coronary revascularization, and coronary death. Cohorts conducted race-specific analyses adjusted for age, sex, smoking, education, body mass index, blood cell type proportions, and technical variables. We conducted fixed-effect meta-analyses across cohorts. RESULTS: Among 11 461 individuals (mean age 64 years, 67% women, 35% African American) free of CHD at baseline, 1895 developed CHD during a mean follow-up of 11.2 years. Methylation levels at 52 CpG (cytosine-phosphate-guanine) sites were associated with incident CHD or myocardial infarction (false discovery rate<0.05). These CpGs map to genes with key roles in calcium regulation (ATP2B2, CASR, GUCA1B, HPCAL1), and genes identified in genome- and epigenome-wide studies of serum calcium (CASR), serum calcium-related risk of CHD (CASR), coronary artery calcified plaque (PTPRN2), and kidney function (CDH23, HPCAL1), among others. Mendelian randomization analyses supported a causal effect of DNA methylation on incident CHD; these CpGs map to active regulatory regions proximal to long non-coding RNA transcripts. CONCLUSION: Methylation of blood-derived DNA is associated with risk of future CHD across diverse populations and may serve as an informative tool for gaining further insight on the development of CHD.

DNA Methylation Analysis Identifies Loci for Blood Pressure Regulation.

Genome-wide association studies have identified hundreds of genetic variants associated with blood pressure (BP), but sequence variation accounts for a small fraction of the phenotypic variance. Epigenetic changes may alter the expression of genes involved in BP regulation and explain part of the missing heritability. We therefore conducted a two-stage meta-analysis of the cross-sectional associations of systolic and diastolic BP with blood-derived genome-wide DNA methylation measured on the Infinium HumanMethylation450 BeadChip in 17,010 individuals of European, African American, and Hispanic ancestry. Of 31 discovery-stage cytosine-phosphate-guanine (CpG) dinucleotides, 13 replicated after Bonferroni correction (discovery: N = 9,828, p < 1.0 × 10-7; replication: N = 7,182, p < 1.6 × 10-3). The replicated methylation sites are heritable (h2 > 30%) and independent of known BP genetic variants, explaining an additional 1.4% and 2.0% of the interindividual variation in systolic and diastolic BP, respectively. Bidirectional Mendelian randomization among up to 4,513 individuals of European ancestry from 4 cohorts suggested that methylation at cg08035323 (TAF1B-YWHAQ) influences BP, while BP influences methylation at cg00533891 (ZMIZ1), cg00574958 (CPT1A), and cg02711608 (SLC1A5). Gene expression analyses further identified six genes (TSPAN2, SLC7A11, UNC93B1, CPT1A, PTMS, and LPCAT3) with evidence of triangular associations between methylation, gene expression, and BP. Additional integrative Mendelian randomization analyses of gene expression and DNA methylation suggested that the expression of TSPAN2 is a putative mediator of association between DNA methylation at cg23999170 and BP. These findings suggest that heritable DNA methylation plays a role in regulating BP independently of previously known genetic variants.

Epigenome-wide association study of rheumatoid arthritis identifies differentially methylated loci in B cells.

Epigenetic regulation of immune cell types could be critical for the development and maintenance of autoimmune diseases like rheumatoid arthritis (RA). B cells are highly relevant in RA, since patients express autoantibodies and depleting this cell type is a successful therapeutic approach. Epigenetic variation, such as DNA methylation, may mediate the pathogenic activity of B cells. In this study, we performed an epigenome-wide association study (EWAS) for RA with three different replication cohorts, to identify disease-specific alterations in DNA methylation in B cells. CpG methylation in isolated B lymphocytes was assayed on the Illumina HumanMethylation450 BeadChip in a discovery cohort of RA patients (N = 50) and controls (N = 75). Differential methylation was observed in 64 CpG sites (q < 0.05). Six biological pathways were also differentially methylated in RA B cells. Analysis in an independent cohort of patients (N = 15) and controls (N = 15) validated the association of 10 CpG sites located on 8 genes CD1C, TNFSF10, PARVG, NID1, DHRS12, ITPK1, ACSF3 and TNFRSF13C, and 2 intergenic regions. Differential methylation at the CBL signaling pathway was replicated. Using an additional case-control cohort (N = 24), the association between RA risk and CpGs cg18972751 at CD1C (P = 2.26 × 10-9) and cg03055671 at TNFSF10 (P = 1.67 × 10-8) genes was further validated. Differential methylation at genes CD1C, TNFSF10, PARVG, NID1, DHRS12, ITPK1, ACSF3, TNFRSF13C and intergenic region chr10p12.31 was replicated in a cohort of systemic lupus erythematosus (SLE) patients (N = 47) and controls (N = 56). Our results highlight genes that may drive the pathogenic activity of B cells in RA and suggest shared methylation patterns with SLE.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels.

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Menopause accelerates biological aging.

Although epigenetic processes have been linked to aging and disease in other systems, it is not yet known whether they relate to reproductive aging. Recently, we developed a highly accurate epigenetic biomarker of age (known as the "epigenetic clock"), which is based on DNA methylation levels. Here we carry out an epigenetic clock analysis of blood, saliva, and buccal epithelium using data from four large studies: the Women's Health Initiative (n = 1,864); Invecchiare nel Chianti (n = 200); Parkinson's disease, Environment, and Genes (n = 256); and the United Kingdom Medical Research Council National Survey of Health and Development (n = 790). We find that increased epigenetic age acceleration in blood is significantly associated with earlier menopause (P = 0.00091), bilateral oophorectomy (P = 0.0018), and a longer time since menopause (P = 0.017). Conversely, epigenetic age acceleration in buccal epithelium and saliva do not relate to age at menopause; however, a higher epigenetic age in saliva is exhibited in women who undergo bilateral oophorectomy (P = 0.0079), while a lower epigenetic age in buccal epithelium was found for women who underwent menopausal hormone therapy (P = 0.00078). Using genetic data, we find evidence of coheritability between age at menopause and epigenetic age acceleration in blood. Using Mendelian randomization analysis, we find that two SNPs that are highly associated with age at menopause exhibit a significant association with epigenetic age acceleration. Overall, our Mendelian randomization approach and other lines of evidence suggest that menopause accelerates epigenetic aging of blood, but mechanistic studies will be needed to dissect cause-and-effect relationships further.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study.

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Fine-mapping of lipid regions in global populations discovers ethnic-specific signals and refines previously identified lipid loci.

Genome-wide association studies have identified over 150 loci associated with lipid traits, however, no large-scale studies exist for Hispanics and other minority populations. Additionally, the genetic architecture of lipid-influencing loci remains largely unknown. We performed one of the most racially/ethnically diverse fine-mapping genetic studies of HDL-C, LDL-C, and triglycerides to-date using SNPs on the MetaboChip array on 54,119 individuals: 21,304 African Americans, 19,829 Hispanic Americans, 12,456 Asians, and 530 American Indians. The majority of signals found in these groups generalize to European Americans. While we uncovered signals unique to racial/ethnic populations, we also observed systematically consistent lipid associations across these groups. In African Americans, we identified three novel signals associated with HDL-C (LPL, APOA5, LCAT) and two associated with LDL-C (ABCG8, DHODH). In addition, using this population, we refined the location for 16 out of the 58 known MetaboChip lipid loci. These results can guide tailored screening efforts, reveal population-specific responses to lipid-lowering medications, and aid in the development of new targeted drug therapies.

Signatures of post-zygotic structural genetic aberrations in the cells of histologically normal breast tissue that can predispose to sporadic breast cancer.

Sporadic breast cancer (SBC) is a common disease without robust means of early risk prediction in the population. We studied 282 females with SBC, focusing on copy number aberrations in cancer-free breast tissue (uninvolved margin, UM) outside the primary tumor (PT). In total, 1162 UMs (1-14 per breast) were studied. Comparative analysis between UM(s), PT(s), and blood/skin from the same patient as a control is the core of the study design. We identified 108 patients with at least one aberrant UM, representing 38.3% of cases. Gains in gene copy number were the principal type of mutations in microscopically normal breast cells, suggesting that oncogenic activation of genes via increased gene copy number is a predominant mechanism for initiation of SBC pathogenesis. The gain of ERBB2, with overexpression of HER2 protein, was the most common aberration in normal cells. Five additional growth factor receptor genes (EGFR, FGFR1, IGF1R, LIFR, and NGFR) also showed recurrent gains, and these were occasionally present in combination with the gain of ERBB2. All the aberrations found in the normal breast cells were previously described in cancer literature, suggesting their causative, driving role in pathogenesis of SBC. We demonstrate that analysis of normal cells from cancer patients leads to identification of signatures that may increase risk of SBC and our results could influence the choice of surgical intervention to remove all predisposing cells. Early detection of copy number gains suggesting a predisposition toward cancer development, long before detectable tumors are formed, is a key to the anticipated shift into a preventive paradigm of personalized medicine for breast cancer.

Hypermethylation of MIR21 in CD4+ T cells from patients with relapsing-remitting multiple sclerosis associates with lower miRNA-21 levels and concomitant up-regulation of its target genes.

BACKGROUND: Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system caused by genetic and environmental factors. DNA methylation, an epigenetic mechanism that controls genome activity, may provide a link between genetic and environmental risk factors. OBJECTIVE: We sought to identify DNA methylation changes in CD4+ T cells in patients with relapsing-remitting (RR-MS) and secondary-progressive (SP-MS) disease and healthy controls (HC). METHODS: We performed DNA methylation analysis in CD4+ T cells from RR-MS, SP-MS, and HC and associated identified changes with the nearby risk allele, smoking, age, and gene expression. RESULTS: We observed significant methylation differences in the VMP1/MIR21 locus, with RR-MS displaying higher methylation compared to SP-MS and HC. VMP1/MIR21 methylation did not correlate with a known MS risk variant in VMP1 or smoking but displayed a significant negative correlation with age and the levels of mature miR-21 in CD4+ T cells. Accordingly, RR-MS displayed lower levels of miR-21 compared to SP-MS, which might reflect differences in age between the groups, and healthy individuals and a significant enrichment of up-regulated miR-21 target genes. CONCLUSION: Disease-related changes in epigenetic marking of MIR21 in RR-MS lead to differences in miR-21 expression with a consequence on miR-21 target genes.

Recurrent evolution of melanism in South American felids.

Morphological variation in natural populations is a genomic test bed for studying the interface between molecular evolution and population genetics, but some of the most interesting questions involve non-model organisms that lack well annotated reference genomes. Many felid species exhibit polymorphism for melanism but the relative roles played by genetic drift, natural selection, and interspecies hybridization remain uncertain. We identify mutations of Agouti signaling protein (ASIP) or the Melanocortin 1 receptor (MC1R) as independent causes of melanism in three closely related South American species: the pampas cat (Leopardus colocolo), the kodkod (Leopardus guigna), and Geoffroy's cat (Leopardus geoffroyi). To assess population level variation in the regions surrounding the causative mutations we apply genomic resources from the domestic cat to carry out clone-based capture and targeted resequencing of 299 kb and 251 kb segments that contain ASIP and MC1R, respectively, from 54 individuals (13-21 per species), achieving enrichment of ~500-2500-fold and ~150x coverage. Our analysis points to unique evolutionary histories for each of the three species, with a strong selective sweep in the pampas cat, a distinctive but short melanism-specific haplotype in the Geoffroy's cat, and reduced nucleotide diversity for both ancestral and melanism-bearing chromosomes in the kodkod. These results reveal an important role for natural selection in a trait of longstanding interest to ecologists, geneticists, and the lay community, and provide a platform for comparative studies of morphological variation in other natural populations.