Genome-wide linkage and positional association analyses identify associations of novel AFF3 and NTM genes with triglycerides: the GenSalt study

Whole genome sequencing of mouse lines divergently selected for fatness (FLI) and leanness (FHI) revealed several genetic variants as candidates for novel obesity genes

BACKGROUND

Analysing genomes of animal model organisms is widely used for understanding the genetic basis of complex traits and diseases, such as obesity, for which only a few mouse models exist, however, without their lean counterparts.

OBJECTIVE

To analyse genetic differences in the unique mouse models of polygenic obesity (Fat line) and leanness (Lean line) originating from the same base population and established by divergent selection over more than 60 generations.

METHODS

Genetic variability was analysed using WGS. Variants were identified with GATK and annotated with Ensembl VEP. g.Profiler, WebGestalt, and KEGG were used for GO and pathway enrichment analysis. miRNA seed regions were obtained with miRPathDB 2.0, LncRRIsearch was used to predict targets of identified lncRNAs, and genes influencing adipose tissue amount were searched using the IMPC database.

RESULTS

WGS analysis revealed 6.3 million SNPs, 1.3 million were new. Thousands of potentially impactful SNPs were identified, including within 24 genes related to adipose tissue amount. SNP density was highest in pseudogenes and regulatory RNAs. The Lean line carries SNP rs248726381 in the seed region of mmu-miR-3086-3p, which may affect fatty acid metabolism. KEGG analysis showed deleterious missense variants in immune response and diabetes genes, with food perception pathways being most enriched. Gene prioritisation considering SNP GERP scores, variant consequences, and allele comparison with other mouse lines identified seven novel obesity candidate genes: 4930441H08Rik, Aff3, Fam237b, Gm36633, Pced1a, Tecrl, and Zfp536.

CONCLUSION

WGS revealed many genetic differences between the lines that accumulated over the selection period, including variants with potential negative impacts on gene function. Given the increasing availability of mouse strains and genetic polymorphism catalogues, the study is a valuable resource for researchers to study obesity.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Plasma PFOA and PFOS Levels, DNA Methylation, and Blood Lipid Levels: A Pilot Study

Exposure to perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) is associated with blood lipids in adults, but the underlying mechanisms remain unclear. This pilot study aimed to investigate the associations between PFOA or PFOS and epigenome-wide DNA methylation and assess the mediating effect of DNA methylation on the PFOA/PFOS-blood lipid association. We measured plasma PFOA/PFOS and leukocyte DNA methylation in 98 patients enrolled from the hospital between October 2018 and August 2019. The median plasma PFOA/PFOS levels were 0.85 and 2.29 ng/mL. Plasma PFOA and PFOS levels were significantly associated with elevated total cholesterol (TC) and low-density lipoprotein cholesterol (LDL) levels. There were 63/87 CpG positions and 8/11 differentially methylated regions (DMRs) associated with plasma PFOA/PFOS levels, respectively. In addition, 5 CpG positions (annotated to AFF3, CREB5, NRG2, USF2, and intergenic region) and one DMR annotated to IRF6 may mediate the association between plasma PFOA/PFOS and LDL levels (mediated proportion from 7.29 to 46.77%); two CpG positions may mediate the association between plasma PFOA/PFOS and TC levels (annotated to CREB5 and USF2, mediated proportion is around 30%). The data suggest that PFOA/PFOS exposure alters DNA methylation. More importantly, the association of PFOA/PFOS with lipid indicators was partly mediated by DNA methylation changes in lipid metabolism-related genes.

Genetic Contributors of Incident Stroke in 10,700 African Americans With Hypertension: A Meta-Analysis From the Genetics of Hypertension Associated Treatments and Reasons for Geographic and Racial Differences in Stroke Studies

Background: African Americans (AAs) suffer a higher stroke burden due to hypertension. Identifying genetic contributors to stroke among AAs with hypertension is critical to understanding the genetic basis of the disease, as well as detecting at-risk individuals.

Methods: In a population comprising over 10,700 AAs treated for hypertension from the Genetics of Hypertension Associated Treatments (GenHAT) and Reasons for Geographic and Racial Differences in Stroke (REGARDS) studies, we performed an inverse variance-weighted meta-analysis of incident stroke. Additionally, we tested the predictive accuracy of a polygenic risk score (PRS) derived from a European ancestral population in both GenHAT and REGARDS AAs aiming to evaluate cross-ethnic performance.

Results: We identified 10 statistically significant (p < 5.00E-08) and 90 additional suggestive (p < 1.00E-06) variants associated with incident stroke in the meta-analysis. Six of the top 10 variants were located in an intergenic region on chromosome 18 (LINC01443-LOC644669). Additional variants of interest were located in or near the COL12A1, SNTG1, PCDH7, TMTC1, and NTM genes. Replication was conducted in the Warfarin Pharmacogenomics Cohort (WPC), and while none of the variants were directly validated, seven intronic variants of NTM proximal to our target variants, had a p-value <5.00E-04 in the WPC. The inclusion of the PRS did not improve the prediction accuracy compared to a reference model adjusting for age, sex, and genetic ancestry in either study and had lower predictive accuracy compared to models accounting for established stroke risk factors. These results demonstrate the necessity for PRS derivation in AAs, particularly for diseases that affect AAs disproportionately.

Conclusion: This study highlights biologically plausible genetic determinants for incident stroke in hypertensive AAs. Ultimately, a better understanding of genetic risk factors for stroke in AAs may give new insight into stroke burden and potential clinical tools for those among the highest at risk.

Whole-Genome Methylation Analysis Revealed ART-Specific DNA Methylation Pattern of Neuro- and Immune-System Pathways in Chinese Human Neonates

The DNA methylation of human offspring can change due to the use of assisted reproductive technology (ART). In order to find the differentially methylated regions (DMRs) in ART newborns, cord blood maternal cell contamination and parent DNA methylation background, which will add noise to the real difference, must be removed. We analyzed newborns’ heel blood from six families to identify the DMRs between ART and natural pregnancy newborns, and the genetic model of methylation was explored, meanwhile we analyzed 32 samples of umbilical cord blood of infants born with ART and those of normal pregnancy to confirm which differences are consistent with cord blood data. The DNA methylation level was lower in ART-assisted offspring at the whole genome-wide level. Differentially methylated sites, DMRs, and cord blood differentially expressed genes were enriched in the important pathways of the immune system and nervous system, the genetic patterns of DNA methylation could be changed in the ART group. A total of three imprinted genes and 28 housekeeping genes which were involved in the nervous and immune systems were significant different between the two groups, six of them were detected both in heel blood and cord blood. We concluded that there is an ART-specific DNA methylation pattern involved in neuro- and immune-system pathways of human ART neonates, providing an epigenetic basis for the potential long-term health risks in ART-conceived neonates.

Adverse Cardiovascular Outcomes and Antihypertensive Treatment: A Genome-Wide Interaction Meta-Analysis in the International Consortium for Antihypertensive Pharmacogenomics Studies

We sought to identify genome-wide variants influencing antihypertensive drug response and adverse cardiovascular outcomes, utilizing data from four randomized controlled trials in the International Consortium for Antihypertensive Pharmacogenomics Studies (ICAPS). Genome-wide antihypertensive drug-single nucleotide polymorphism (SNP) interaction tests for four drug classes (β-blockers, n = 9,195; calcium channel blockers (CCBs), n = 10,511; thiazide/thiazide-like diuretics, n = 3,516; ACE-inhibitors/ARBs, n = 2,559) and cardiovascular outcomes (incident myocardial infarction, stroke, or death) were analyzed among patients with hypertension of European ancestry. Top SNPs from the meta-analyses were tested for replication of cardiovascular outcomes in an independent Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) study (n = 21,267), blood pressure (BP) response in independent ICAPS studies (n = 1,552), and ethnic validation in African Americans from the Genetics of Hypertension Associated Treatment study (GenHAT; n = 5,115). One signal reached genome-wide significance in the β-blocker-SNP interaction analysis (rs139945292, Interaction P = 1.56 × 10-8 ). rs139945292 was validated through BP response to β-blockers, with the T-allele associated with less BP reduction (systolic BP response P = 6 × 10-4 , Beta = 3.09, diastolic BP response P = 5 × 10-3 , Beta = 1.53). The T-allele was also associated with increased adverse cardiovascular risk within the β-blocker treated patients' subgroup (P = 2.35 × 10-4 , odds ratio = 1.57, 95% confidence interval = 1.23-1.99). The locus showed nominal replication in CHARGE, and consistent directional trends in β-blocker treated African Americans. rs139945292 is an expression quantitative trait locus for the 50 kb upstream gene NTM (neurotrimin). No SNPs attained genome-wide significance for any other drugs classes. Top SNPs were located near CALB1 (CCB), FLJ367777 (ACE-inhibitor), and CES5AP1 (thiazide). The NTM region is associated with increased risk for adverse cardiovascular outcomes and less BP reduction in β-blocker treated patients. Further investigation into this region is warranted.

A genomic prediction model for racecourse starts in the Thoroughbred horse

Durability traits in Thoroughbred horses are heritable, economically valuable and may affect horse welfare. The aims of this study were to test the hypotheses that (i) durability traits are heritable and (ii) genetic data may be used to predict a horse's potential to have a racecourse start. Heritability for the phenotype 'number of 2- and 3-year-old starts' was estimated to be h m 2  = 0.11 ± 0.02 (n = 4499). A genome-wide association study identified SNP contributions to the trait. The neurotrimin (NTM), opioid-binding protein/cell adhesion molecule like (OPCML) and prolylcarboxypeptidase (PRCP) genes were identified as candidate genes associated with the trait. NTM functions in brain development and has been shown to have been selected during the domestication of the horse. PRCP is an established expression quantitative trait locus involved in the interaction between voluntary exercise and body composition in mice. We hypothesise that variation at these loci contributes to the motivation of the horse to exercise, which may influence its response to the demands of the training and racing environment. A random forest with mixed effects (RFME) model identified a set of SNPs that contributed to 24.7% of the heritable variation in the trait. In an independent validation set (n = 528 horses), the cohort with high genetic potential for a racecourse start had significantly fewer unraced horses (16% unraced) than did low (27% unraced) potential horses and had more favourable race outcomes among those that raced. Therefore, the information from SNPs included in the model may be used to predict horses with a greater chance of a racecourse start.

IGF1 Gene Is Associated With Triglyceride Levels In Subjects With Family History Of Hypertension From The SAPPHIRe And TWB Projects

Chromosome 12q23-q24 has been linked to triglyceride (TG) levels by previous linkage studies, and it contains the Insulin-like growth factor 1 (IGF1) gene. We investigated the association between IGF1 and TG levels using two independent samples collected in Taiwan. First, based on 954 siblings in 397 families from the Stanford Asian Pacific Program in Hypertension and Insulin Resistance (SAPPHIRe), we found that rs978458 was associated with TG levels (β = -0.049, p = 0.0043) under a recessive genetic model. Specifically, subjects carrying the homozygous genotype of the minor allele had lower TG levels, compared with other subjects. Then, a series of stratification analyses in a large sample of 13,193 unrelated subjects from the Taiwan biobank (TWB) project showed that this association appeared in subjects with a family history (FH) of hypertension (β = -0.045, p = 0.0000034), but not in subjects without such an FH. A re-examination of the SAPPHIRe sample confirmed that this association appeared in subjects with an FH of hypertension (β = -0.068, p = 0.0025), but not in subjects without an FH. The successful replication in two independent samples indicated that IGF1 is associated with TG levels in subjects with an FH of hypertension in Taiwan.

Deficit in emotional learning in neurotrimin knockout mice

Neurotrimin (Ntm) belongs to the IgLON family of cell adhesion molecules with Lsamp, Obcam and kilon that regulate the outgrowth of neurites mostly by forming heterodimers. IgLONs have been associated with psychiatric disorders, intelligence, body weight, heart disease and tumours. This study provides an initial behavioural and pharmacological characterization of the phenotype of Ntm-deficient mice. We expected to see at least some overlap with the phenotype of Lsamp-deficient mice as Ntm and Lsamp are the main interaction partners in the IgLON family and are colocalized in some brain regions. However, Ntm-deficient mice displayed none of the deviations in behaviour that we have previously shown in Lsamp-deficient mice, but differently from Lsamp-deficient mice, had a deficit in emotional learning in the active avoidance task. The only overlap was decreased sensitivity to the locomotor stimulating effect of amphetamine in both knockout models. Thus, despite being interaction partners, on the behavioural level Lsamp seems to play a much more central role than Ntm and the roles of these two proteins seem to be complementary rather than overlapping.

A Combination of CD28 (rs1980422) and IRF5 (rs10488631) Polymorphisms Is Associated with Seropositivity in Rheumatoid Arthritis: A Case Control Study

Introduction

The aim of the study was to analyse genetic architecture of RA by utilizing multiparametric statistical methods such as linear discriminant analysis (LDA) and redundancy analysis (RDA).

Methods

A total of 1393 volunteers, 499 patients with RA and 894 healthy controls were included in the study. The presence of shared epitope (SE) in HLA-DRB1 and 11 SNPs (PTPN22 C/T (rs2476601), STAT4 G/T (rs7574865), CTLA4 A/G (rs3087243), TRAF1/C5 A/G (rs3761847), IRF5 T/C (rs10488631), TNFAIP3 C/T (rs5029937), AFF3 A/T (rs11676922), PADI4 C/T (rs2240340), CD28 T/C (rs1980422), CSK G/A (rs34933034) and FCGR3A A/C (rs396991), rheumatoid factor (RF), anti–citrullinated protein antibodies (ACPA) and clinical status was analysed using the LDA and RDA.

Results

HLA-DRB1, PTPN22, STAT4, IRF5 and PADI4 significantly discriminated between RA patients and healthy controls in LDA. The correlation between RA diagnosis and the explanatory variables in the model was 0.328 (Trace = 0.107; F = 13.715; P = 0.0002). The risk variants of IRF5 and CD28 genes were found to be common determinants for seropositivity in RDA, while positivity of RF alone was associated with the CTLA4 risk variant in heterozygous form. The correlation between serologic status and genetic determinants on the 1st ordinal axis was 0.468, and 0.145 on the 2nd one (Trace = 0.179; F = 6.135; P = 0.001). The risk alleles in AFF3 gene together with the presence of ACPA were associated with higher clinical severity of RA.

Conclusions

The association among multiple risk variants related to T cell receptor signalling with seropositivity may play an important role in distinct clinical phenotypes of RA. Our study demonstrates that multiparametric analyses represent a powerful tool for investigation of mutual relationships of potential risk factors in complex diseases such as RA.